
Deeply bound KNN, KNNN and KNNNN states are discussed. The effective force
exerted by the K meson on the nucleons is calculated with static nucleons. Next
the binding energies are obtained by solving the Schrodinger equation or by
variational calculations.
The dominant attraction comes from the Swave Lambda(1405) and an additional
contribution is due to Sigma(1385). The latter state is formed at the nuclear
peripheries and absorbs a sizable piece of the binding energy. It also
generates new branches of quasibound states. The lowest binding energies based
on a phenomenological KN input fall into the 4080 MeV range for KNN, 90150
MeV for KNNN and 120220 MeV for Kalpha systems. The uncertainties are due to
unknown KN interactions in the distant subthreshold energy region.

The coupled etaN, piN, gammaN, 2piN system is described by a Kmatrix
method. The parameters in this model are adjusted to get an optimal fit to
piN>piN, piN>etaN, gammaN>piN and gammaN>etaN data in an energy
range of about 100 MeV or so each side of the etathreshold. Compared with our
earlier analysis, we now utilize recent Crystal Ball data. However, the outcome
confirms our previous result that the etanucleon scattering length (a) is
large with a value of 0.91(6)+i0.27(2) fm.

This is a followup to our earlier work for the charge (vector) and matter
(scalar) distributions for Swave states in a heavylight meson, where the
heavy quark is static and the light quark has a mass about that of the strange
quark. The calculation is again carried out with dynamical fermions on a
16^3x24 lattice with a lattice spacing of about 0.14 fm. It is shown that
several features of the S and Pwave distributions are in qualitative
agreement with what one expects from a simple onebody Dirac equation
interpretation.

In this Chapter QCD interactions between a quark and an antiquark are
discussed. In the heavy quark limit these potentials can be related to
quarkonia and $1/m$ corrections can be systematically determined. Excitations
of the ground state potential provide an entry point into the phenomenology of
quarkgluon hybrids. The shortdistance behaviour of nonperturbative
potentials can serve as a test of resummation and convergence of perturbative
expansions. Torelons and potentials between nonfundamental colour charges
offer a window into the origin of the confinement mechanism and relate to
effective string descriptions of low energy aspects of QCD.

In our earlier work, the charge and matter radial distributions of
heavylight mesons were measured on a 16^3 times 24 lattice with a lattice
spacing of a = 0.17 fm and a light quark mass about that of the strange quark.
Several major improvements have now been made: 1) Dynamical fermions are used
with a approximately 0.14 fm; 2) More gauge configurations are included (78 vs
20); 3) Offaxis, in addition to onaxis, insertions are made; 4) The data
analysis is much more complete. In particular, distributions involving excited
states are extracted.

The etadeuteron scattering amplitude is calculated with a recent Kmatrix
model of etaN interactions. The existence of a narrow virtual state in the
etadeuteron system is inferred. This state couples strongly to the etad
channel but its effect on final state interactions in the pn > etad reaction
is shown to be rather weak. Fairly large etaN scattering lengths are allowed
by the experimental data.

The coupled etaN, piN, gammaN system is described by a Kmatrix method.
The parameters in this model are adjusted to get an optimal fit to piN >pi
N$, piN >etaN$, gammaN>piN and gammaN >etaN data in an energy
range of about 100MeV or so each side of the eta threshold. The outcome is the
appearance of two solutions one which has an etanucleon scattering length (a)
of about 1.0 fm and a second with a~0.2fm. However, the second solution has an
unconventional nonLorentian form for the Tmatrix in the region of the
1535(20)MeV and 1650(30)MeV Swave piN resonances.

Recent lattice calculations on the interaction energy of two heavylight
mesons (Q^2\bar{q}^2) are interpreted in terms of the potential of the
corresponding single heavylight meson (Q\bar{q}). This model leads to a large
overestimate of the binding compared with the lattice data  unless the basic
Q\bar{q} potential is modified to become a fourquark potential.

The potential between staticlight mesons forming a mesonmeson or a
mesonantimeson system is calculated in quenched and unquenched SU(3) gauge
theory. We use the SheikholeslamiWohlert action and statistical estimators of
light quark propagators with maximal variance reduction. The dependence of the
potentials on the light quark spin and isospin and the effect of meson exchange
is investigated. Our main motivation is exploration of bound states of two
mesons and string breaking. The latter also involves the twoquark potential
and the correlation between twoquark and twomeson states.

By comparing lattice data for the two heavylight meson system (Q^2 qbar^2)
with a standard manybody approach employing only interquark potentials, it is
shown that the use of unmodified twoquark potentials leads to a gross
overestimate of the binding energy.

The coupled pi+N, eta+N, gamma+N systems are described by a Kmatrix method.
The parameters in this model are adjusted to get an optimal fit to pi+N >
pi+N, pi+N > eta+N, gamma+N > pi+N and gamma+N > eta+N data in an energy
range of about 100 MeV each side of the etathreshold. The coupling of photons
to the N(1535) state is extracted and also an alternative to the current
S11(gamma+N > pi+N) amplitudes suggested. Expansions are given for the etaeta
and gammaeta amplitudes in terms of the etamomentum. Effects of interference
of this state with background potential interactions are discussed and
experimental consequences are indicated.

In order to understand the binding of four static quarks, flux distributions
corresponding to these binding energies are studied in quenched SU(2) and also
related to a model for the energies. The potential relevant to string breaking
between two heavylight mesons is measured in quenched SU(3) using stochastic
estimates of light quark propagators with the SheikholeslamiWohlert action.

The optical potential for the $\eta$meson is calculated. This is based on a
new determination of a large $\eta$N scattering length and recent experiments
of $\eta$ production in two nucleon collisions. These determine the $\eta$ 
nuclear potential well to be about 45 MeV deep and able to support S, P, D and
F quasibound states. Some of the higher angular momentum states are fairly
narrow and offer a chance of experimental detection.

The full spatial distribution of the color fields of two and four static
quarks is measured in lattice SU(2) field theory at separations up to 1 fm at
beta=2.4. The fourquark case is equivalent to a qbar q qbar q system in SU(2)
and is relevant to mesonmeson interactions. By subtracting twobody flux tubes
from the fourquark distribution we isolate the flux contribution connected
with the fourbody binding energy. This contribution is further studied using a
model for the binding energies. Lattice sum rules for two and four quarks are
used to verify the results.

Action and energy fluxtube profiles are computed, in SU(2) with
beta=2.4,2.5, for two quarks up to 1 fm apart and for which the colour fields
are in their ground state (A_1g) and the first (E_u) and higher (A'_1g) excited
gluonic states. When these profiles are integrated over all space, a scaling
comparison is made between the beta=2.4 and 2.5 data. Using sum rules, these
integrated forms also permit an estimate to be made of generalised
betafunctions giving b(2.4)=0.312(15), b(2.5)=0.323(9), f(2.4)=0.65(1) and
f(2.5)=0.68(1). When the profiles are integrated only over planes transverse to
the interquark line and assuming underlying string features, scaling
comparisons are again made near the centres of the interquark line for the
largest interquark distances. For the A'_{1g} case, some of the profiles
exhibit a 'diplike' structure characteristic of the IsgurPaton model.

The coupled eta N, pi N system is described by a Kmatrix method. The
parameters in this model are adjusted to get an optimal fit to pi N > pi N, pi
N > eta N and gamma N > eta N data in an energy range of about 100MeV each
side of the eta threshold. In the notation T^(1)+iq_(eta) = 1/a + (r_0)/2
q^2_(eta) + s q^4_(eta), q_(eta) being the momentum in the eta N
centerofmass, the resulting effective range parameters for eta N scattering
are found to be a(fm) = 0.75(4)i0.27(3), r_0(fm) = 1.50(13)i0.24(4) and
s(fm^3) = 0.10(2)i0.01(1)

The spatial distribution of the action and energy in the colour fields of
fluxtubes is studied in lattice SU(2) field theory for static quarks at
separations up to 1 fm. Special attention is paid to the structure of the
colour fields associated with an excited fluxtube with E_u symmetry. We
compare our results with hadronic string and fluxtube models. Sum rules are
used to extract generalised betafunctions, to describe the expected colour
field behaviour and to cross check the methods used.

Sum rules  relating the static quark potential V(R) to the spatial
distribution of the action and energy in the colour fields of fluxtubes  are
applied in three ways: 1) To extract generalised betafunctions: 2) As a
consistency check for the use of excited gluon fluxtubes and as an estimate of
the quark selfenergies: 3) To extract approximate sum rules using a simplified
form of V(R). Also the fluxtube profiles are compared with hadronic string and
fluxtube models.

We use lattice sum rules for the static quark potential to determine the
betafunction for symmetric and asymmetric lattices nonperturbatively. We also
study the colour field distributions in excited gluonic states.

Etadeuteron scattering lengths are calculated. A summation of the multiple
scattering series is carried out and the result is checked against more
involved calculations. The necessity to go beyond the fixed nucleon
approximation is emphasized. It is shown that a quasibound or virtual state in
the $\eta$deuteron system may occur within the range of $\eta$nucleon
scattering lengths suggested by other experiments.

A first attempt to understand hadron dynamics at low energies in terms of the
fundamental quark and gluon degrees of freedom incorporates the effects of the
gluonic field into a potential depending only on the spatial positions of the
quarks, which are considered in the infinite mass limit. A suitable framework
for calculating such potentials between static quarks, i.e.\ a generalization
of the Wilson loop will be discussed.
Making a connection with recent Monte Carlo lattice simulations for the
lowest two energies of a system of two quarks and two antiquarks, the static
$qq\bar{q}\bar{q}$potential will be calculated in perturbation theory to
fourth order. The result will be shown to be exactly equal to the prediction of
a straightforward twobody approach, which in Monte Carlo lattice simulations
has been found to be a reasonable approximation for very small interquark
distances.

Using multiple scattering theory the scattering lengths of $\eta$ mesons on
helium nuclei are calculated and checked against final state $\eta$
interactions from the $pd \rightarrow \eta ^3$He and $dd \rightarrow \eta ^4$He
reactions. The existence of an $\eta^4$He quasibound state is indicated.