
We study multihadron 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 $(j1/2)^{P}$ for $j \neq 0$, known in a normal
hadron, can be generalized to multihadron systems. The spin degeneracy is the
universal phenomena for any multihadron systems with a single heavy quark,
irrespective of their internal structures, including compact multiquarks,
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 axialvector coupling including
1/Mcorrections for the heavy meson mass $M$ by following the
velocityrearrangement 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 inmedium masses modified by the interaction with
nucleons via pion exchanges including the 1/Mcorrections. 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 oneloop 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 channelcouplings 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 JPARC
(Japan Proton Accelerator Research Complex) facility and GSIFAIR (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
EfremovRadyushkinBrodskyLepage (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 highenergy 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 semirelativistic quark potential model is surprisingly powerful for
heavylight 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
semirelativistic model succeeds in predicting and reproducing all the mass
spectra of heavylight 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 socalled 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 squark fragmentation functions can distinguish the
tetraquark structure from $q\bar q$. By the global analysis of f_0 (980)
production data in electronpositron 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 hadronproduction data in $e^+e^$ annihilation. It is particularly
important that uncertainties of the fragmentation functions are estimated for
the first time. We found that lightquark and gluon fragmentation functions
have large uncertainties, so that one should be careful in using these
functions for hadronproduction processes in heavyion 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 chargedhadron production data in electronpositron annihilation. It
is important that uncertainties of the determined fragmentation functions are
estimated in this analysis. Analysis results indicate that gluon and
lightquark functions have large uncertainties especially at small Q^2. We find
that nexttoleadingorder (NLO) uncertainties are significantly reduced in
comparison with leadingorder (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 chargedhadron production
data in electronpositron 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 lightquark
fragmentation functions, have large uncertainties at small $Q^2$. There are
large differences between widelyused 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
nexttoleadingorder (NLO) analyses for the pion and kaon in comparison with
leadingorder 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 hadronproduction data at small $Q^2$ or $p_T$ ($Q^2,
p_T^2 << M_Z^2$) in lepton scattering and hadronhadron 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 chargedhadron production data in electronpositron
annihilation. The optimum functions are obtained in both leading order (LO) and
nexttoleading 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
helicitytohelicity correlations. The methodical issue of using spin transfers
for crosschecks of systematic errors in crosssection $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
coloroctet model, which is based on NRQCD factorization formalism. We show
that the measurement of the twospin asymmetry of this process might be
effective to constrain the value of NRQCD matrix elements.

To test a possible coloroctet contribution of heavy quarkonium productions,
we propose $\psi$'productions for small$p_T$ regions, where the colorsinglet
$\psi$' state cannot be produced by the gluongluon fusion owing to charge
conjugation. We also point out that these processes are very useful to deduce
informations of spindependent gluon distributions $\Delta g$.