• We study multi-hadron systems with a single heavy quark (charm or bottom) in the limit of heavy quark mass. The spin degeneracy of the states with quantum numbers $(j+1/2)^{P}$ and $(j-1/2)^{P}$ for $j \neq 0$, known in a normal hadron, can be generalized to multi-hadron systems. The spin degeneracy is the universal phenomena for any multi-hadron systems with a single heavy quark, irrespective of their internal structures, including compact multi-quarks, hadronic molecules and exotic nuclei. We demonstrate the spin degeneracy in the hadronic systems formed by a heavy hadron effective theory; $P^{(\ast)}N$ states with a $P^{(\ast)}=\bar{D}^{(\ast)}$, $B^{(\ast)}$ meson and a nucleon $N$, and a $P^{(\ast)}$ meson in nuclear matter.
  • We consider heavy mesons with charm and bottom in nuclear medium. We construct the effective Lagrangian with axial-vector coupling including 1/M-corrections for the heavy meson mass $M$ by following the velocity-rearrangement invariance. As an application, we consider heavy mesons, $\bar{\mathrm{D}}$ and $\bar{\mathrm{D}}^{\ast}$ mesons for charm and $\mathrm{B}$ and $\mathrm{B}^{\ast}$ mesons for bottom, bound in nuclear matter, and we discuss their in-medium masses modified by the interaction with nucleons via pion exchanges including the 1/M-corrections. The mass modifications are affected by the gluon dynamics in nuclear medium. By comparison with the heavy quark effective theory, we find that the effects of scale anomaly become suppressed in nuclear medium. We also find that the contributions from the chromoelectric gluon are enhanced in nuclear medium, while those from the chromomagnetic gluon are reduced. We propose to use heavy mesons as probes to research the gluon fields in nuclear medium in experimental studies.
  • We discuss the dynamics in finite density medium including a heavy impurity particle (hadron or quark) with a heavy flavor, charm and bottom, at zero temperature. As a system, we consider a $\bar{D}$ ($B$) meson embedded in nuclear matter as a heavy impurity boson with SU(2) isospin symmetry. As another system, we also consider a charm (bottom) quark embedded in quark matter as a heavy impurity fermion with SU(3) color symmetry. We suppose a vector current interaction with SU($n$) symmetry ($n \ge 2$) for the fermion composing the Fermi surface and the embedded heavy impurity particle, and calculate the scattering amplitude perturbatively for the small coupling constant up to one-loop level. We obtain that the scattering amplitude has a logarithmic enhancement in the large mass limit of the heavy impurity particle, and show that the perturbative calculation breaks down for any small coupling constant in this limit.
  • We study heavy hadron spectroscopy near heavy meson thresholds. We employ heavy pseudoscalar meson P and heavy vector meson P* as effective degrees of freedom and consider meson exchange potentials between them. All possible composite states which can be constructed from the P and P* mesons are studied up to the total angular momentum J <= 2. We consider, as exotic states, isosinglet states with exotic J^{PC} quantum numbers and isotriplet states. We solve numerically the Schr\"odinger equation with channel-couplings for each state. We found B(*)barB(*) molecule states for I^G(J^{PC}) = 1^+(1^{+-}) correspond to the masses of twin resonances Zb(10610) and Zb(10650). We predict several possible B(*)barB(*) bound and/or resonant states in other channels. On the other hand, there are no B(*)barB(*) bound and/or resonant states whose quantum numbers are exotic.
  • We discuss the mass modifications of $\bar{\mathrm{D}}$ and $\bar{\mathrm{D}}^{\ast}$ ($\mathrm{B}$ and $\mathrm{B}^{\ast}$) mesons in nuclear medium. The heavy quark symmetry for $\bar{\mathrm{D}}$ and $\bar{\mathrm{D}}^{\ast}$ ($\mathrm{B}$ and $\mathrm{B}^{\ast}$) mesons is adopted, and the interaction between a $\bar{\mathrm{D}}$ or $\bar{\mathrm{D}}^{\ast}$ ($\mathrm{B}$ or $\mathrm{B}^{\ast}$) meson and a nucleon is supplied from the pion exchange. We find the negative mass shifts for $\bar{\mathrm{D}}$ meson and $\mathrm{B}$ meson, and hence that the $\bar{\mathrm{D}}$ and $\mathrm{B}$ mesons are bound in the nuclear medium. As applications, we consider the atomic nuclei with a $\bar{\mathrm{D}}$ meson, $^{40}_{\bar{\mathrm{D}}}\mathrm{Ca}$ and $^{208}_{\bar{\mathrm{D}}}\mathrm{Pb}$, and investigate the energy levels of the $\bar{\mathrm{D}}$ meson in each nucleus. We also discuss the mass shifts in the isospin asymmetric nuclear medium, and present the possible phenomenon about distribution of isospin density around a $\bar{\mathrm{D}}$ or $\mathrm{B}$ meson in nuclear medium. We find that the mass shifts of $\bar{\mathrm{D}}^{\ast}$ and $\mathrm{B}^{\ast}$ mesons have large imaginary parts, which would prevent precise study of the energy levels of $\bar{\mathrm{D}}^{\ast}$ and $\mathrm{B}^{\ast}$ mesons in nuclei.
  • We consider a novel class of hard branching hadronic processes a+b->c+d+e, where hadrons c and d have large and nearly opposite transverse momenta and large invariant energy which is a finite fraction of the total invariant energy. We use color transparency logic to argue that these processes can be used to study quark generalized parton distributions (GPDs) for baryons and mesons in hadron collisions hence complementing and adding to the studies of GPDs in the exclusive DIS (deep inelastic scattering) processes. We propose that a number of GPDs can be investigated in hadron facilities such as J-PARC (Japan Proton Accelerator Research Complex) facility and GSI-FAIR (Gesellschaft fur Schwerionenforschung -Facility for Antiproton and Ion Research) project. In this work, the GPDs for the nucleon and for the N->Delta transition are studied in the reaction N+N->N+pi+B, where N, pi, and B are a nucleon, a pion, and a baryon (nucleon or Delta), respectively, with a large momentum transfer between B (or pi) and the incident nucleon. In particular, the Efremov-Radyushkin-Brodsky-Lepage (ERBL) region of the GPDs can be measured in such exclusive reactions. We estimate the cross section of the processes N+N->N+pi+B by using current models for relevant GPDs and information about large angle pi N reactions. We find that it will be feasible to measure these cross sections at the high-energy hadron facilities, and get novel information about the nucleon structure, for example, contributions of quark orbital angular momenta to the nucleon spin. The studies of N->Delta transition GPDs could be valuable also for investigating electromagnetic properties of the transition.
  • The semi-relativistic quark potential model is surprisingly powerful for heavy-light systems if the bound state equation is treated correctly using 1/m_Q expansion with heavy quark mass m_Q. We elucidate the reasons why our semi-relativistic model succeeds in predicting and reproducing all the mass spectra of heavy-light systems so far reported, D/D_s/B/B_s, by reviewing and comparing recent experimental data with the results of our model and others. Especially the mass spectra of the so-called D_{sJ}, i.e., D_{s0}^* and D_{s1}', are successfully reproduced only by our model but not by other models.
  • We discuss internal structure of an exotic hadron by using fragmentation functions. The fragmentation functions for the f_0(980) meson are obtained by a global analysis of e^++e^- \to f_0+X data. Quark configuration of the f_0(980) could be determined by peak positions and second moments of the obtained fragmentation functions.
  • It is proposed that fragmentation functions should be used to identify exotic hadrons. As an example, fragmentation functions of the scalar meson f_0(980) are investigated. It is pointed out that the second moments and functional forms of the u- and s-quark fragmentation functions can distinguish the tetraquark structure from $q\bar q$. By the global analysis of f_0 (980) production data in electron-positron annihilation, its fragmentation functions and their uncertainties are determined. It is found that the current available data are not sufficient to determine its internal structure, while precise data in future should be able to identify exotic quark configurations.
  • Fragmentation functions of pion, kaon, and nucleon are determined by global analyses of hadron-production data in $e^+e^-$ annihilation. It is particularly important that uncertainties of the fragmentation functions are estimated for the first time. We found that light-quark and gluon fragmentation functions have large uncertainties, so that one should be careful in using these functions for hadron-production processes in heavy-ion collisions and lepton scattering. The analysis is extended to possible exotic hadron search by fragmentation functions. We found that internal structure of $f_0 (980)$, such as $s\bar s$ or tetraquark configuration, can be determined by noting differences between favored and disfavored fragmentation functions.
  • Fragmentation functions are determined for the pion, kaon, and proton by analyzing charged-hadron production data in electron-positron annihilation. It is important that uncertainties of the determined fragmentation functions are estimated in this analysis. Analysis results indicate that gluon and light-quark functions have large uncertainties especially at small Q^2. We find that next-to-leading-order (NLO) uncertainties are significantly reduced in comparison with leading-order (LO) ones in the pion and kaon. The fragmentation functions are very different in various analysis groups. However, all the recent functions are roughly within the estimated uncertainties, which indicates that they are consistent with each other. We provide a code for calculating the fragmentation functions and their uncertainties at a given kinematical point of z and Q^2 by a user.
  • Fragmentation functions and their uncertainties are determined for pion, kaon, and proton by a global $\chi^2$ analysis of charged-hadron production data in electron-positron annihilation and by the Hessian method for error estimation. It is especially important that the uncertainties of the fragmentation functions are estimated in this analysis. The results indicate that the fragmentation functions, especially gluon and light-quark fragmentation functions, have large uncertainties at small $Q^2$. There are large differences between widely-used functions by KKP (Kniehl, Kramer, and P\"otter) and Kretzer; however, they are compatible with each other and also with our functions if the uncertainties are taken into account. We find that determination of the fragmentation functions is improved in next-to-leading-order (NLO) analyses for the pion and kaon in comparison with leading-order ones. Such a NLO improvement is not obvious in the proton. Since the uncertainties are large at small $Q^2$, the uncertainty estimation is very important for analyzing hadron-production data at small $Q^2$ or $p_T$ ($Q^2, p_T^2 << M_Z^2$) in lepton scattering and hadron-hadron collisions. A code is available for general users for calculating obtained fragmentation functions.
  • Fragmentation functions are determined for pions, kaons, and nucleons by a global analysis of charged-hadron production data in electron-positron annihilation. The optimum functions are obtained in both leading order (LO) and next-to-leading order (NLO) of alpha_s. It is important that uncertainties of the fragmentation functions are estimated in this work by the Hessian method. We found that the uncertainties are large at small Q^2 and that they are generally reduced in the NLO in comparison with the LO ones. We supply a code for calculating the fragmentation functions and their uncertainties for the pions, kaons, and nucleons at given z and Q^2.
  • Recent discovery of D_sJ particles, which are considered to be a great trouble by experimentalists as well as theorists, had already been resolved by our potential model proposed some time ago by two of us (T.M. and T.M.), in which the Hamiltonian and wave functions are expanded in 1/(heavy quark mass) respecting heavy quark symmetry. Using our model, we explain how narrow states like D_sJ can be realized, predict properties of 0^+ and 1^+ states of B and B_s heavy mesons, and interpret how global SU(3) symmetry seems to be recovered for these 0^+ and 1^+ heavy mesons.
  • The analysis\cite{ohkuma:1} of helicity transfer to $\Lambda_{c}^{+}$ in polarized proton collisions is extended to the proton helicity correlations with the $\Lambda_{c}^{+}$ transverse polarization in the production plane (parameter $D_{LS}$). The available spin transfer observables for the collisions of {\em two} longitudinally polarized protons are evaluated. It is shown that, in the central region at $\Lambda_{c}^{+}$ transverse momenta of a few GeV/c, $D_{LS}$ parameters are of about the same size as the helicity-to-helicity correlations. The methodical issue of using spin transfers for cross-checks of systematic errors in cross-section $A_{LL}$ measurements at polarized proton colliders is also briefly discussed.
  • Predicted masses of J^P=0^+ and 1^+ states of D_{sJ} and D by a potential model proposed some time ago by two of us (T.M. and T.M.), in which the Hamiltonian and wave functions are expanded in 1/m_Q with m_Q heavy quark mass respecting heavy quark symmetry, have recently been confirmed by BaBar and Belle experiments within one percent accuracy. In this Letter, decay modes of 0^+ and 1^+ states of B and B_s mesons are discussed using the predicted masses of this model.
  • In heavy quark effective theory, heavy mesons which contain a heavy quark (or antiquark) are classified by $s_l^{pi_l}$, i.e., the total angular momentum $s_l$ and the parity $pi_l$ of the light quark degrees of freedom around a static heavy quark. In this case, however, one needs to separately estimate the parity other than the angular momentum of a light quark to describe heavy mesons. A new operator $K$ was proposed some time ago by two of us (T.M. and T.M.). In this Letter, we show that the quantum number $k$ of this operator is enough to describe both the total angular momentum of the light quark degrees of freedom and the parity of a heavy meson, and derive a simple relation between $k$ and $s_l^{pi_l}$.
  • We have calculated the cross section of $\psi'$ photoproduction at small-$p_{T}$ regions in polarized $\gamma p$ collisions for testing the color-octet model, which is based on NRQCD factorization formalism. We show that the measurement of the two-spin asymmetry of this process might be effective to constrain the value of NRQCD matrix elements.
  • To test a possible color-octet contribution of heavy quarkonium productions, we propose $\psi$'-productions for small-$p_T$ regions, where the color-singlet $\psi$' state cannot be produced by the gluon-gluon fusion owing to charge conjugation. We also point out that these processes are very useful to deduce informations of spin-dependent gluon distributions $\Delta g$.