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 S-wave 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 quasi-bound states. The lowest binding energies based on a phenomenological KN input fall into the 40-80 MeV range for KNN, 90-150 MeV for KNNN and 120-220 MeV for K-alpha systems. The uncertainties are due to unknown KN interactions in the distant subthreshold energy region.

The coupled eta-N, pi-N, gamma-N, 2pi-N system is described by a K-matrix 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 or so each side of the eta-threshold. Compared with our earlier analysis, we now utilize recent Crystal Ball data. However, the outcome confirms our previous result that the eta-nucleon scattering length (a) is large with a value of 0.91(6)+i0.27(2) fm.

This is a follow-up to our earlier work for the charge (vector) and matter (scalar) distributions for S-wave states in a heavy-light 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 P-wave distributions are in qualitative agreement with what one expects from a simple one-body Dirac equation interpretation.

In this Chapter QCD interactions between a quark and an anti-quark 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 quark-gluon hybrids. The short-distance behaviour of non-perturbative potentials can serve as a test of resummation and convergence of perturbative expansions. Torelons and potentials between non-fundamental 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 heavy-light 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) Off-axis, in addition to on-axis, insertions are made; 4) The data analysis is much more complete. In particular, distributions involving excited states are extracted.

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

The coupled eta-N, pi-N, gamma-N system is described by a K-matrix 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 100MeV or so each side of the eta threshold. The outcome is the appearance of two solutions one which has an eta-nucleon scattering length (a) of about 1.0 fm and a second with a~0.2fm. However, the second solution has an unconventional non-Lorentian form for the T-matrix in the region of the 1535(20)MeV and 1650(30)MeV S-wave pi-N resonances.

Recent lattice calculations on the interaction energy of two heavy-light mesons (Q^2\bar{q}^2) are interpreted in terms of the potential of the corresponding single heavy-light 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 four-quark potential.

The potential between static-light mesons forming a meson-meson or a meson-antimeson system is calculated in quenched and unquenched SU(3) gauge theory. We use the Sheikholeslami-Wohlert 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 two-quark potential and the correlation between two-quark and two-meson states.

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

The coupled pi+N, eta+N, gamma+N systems are described by a K-matrix 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 eta-threshold. 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 eta-eta and gamma-eta amplitudes in terms of the eta-momentum. 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 heavy-light mesons is measured in quenched SU(3) using stochastic estimates of light quark propagators with the Sheikholeslami-Wohlert 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 four-quark case is equivalent to a qbar q qbar q system in SU(2) and is relevant to meson-meson interactions. By subtracting two-body flux tubes from the four-quark distribution we isolate the flux contribution connected with the four-body 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 flux-tube 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 beta-functions 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 'dip-like' structure characteristic of the Isgur-Paton model.

The coupled eta N, pi N system is described by a K-matrix 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 center-of-mass, 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 flux-tubes 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 flux-tube with E_u symmetry. We compare our results with hadronic string and flux-tube models. Sum rules are used to extract generalised beta-functions, 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 flux-tubes -- are applied in three ways: 1) To extract generalised beta-functions: 2) As a consistency check for the use of excited gluon flux-tubes and as an estimate of the quark self-energies: 3) To extract approximate sum rules using a simplified form of V(R). Also the flux-tube profiles are compared with hadronic string and flux-tube models.

We use lattice sum rules for the static quark potential to determine the beta-function for symmetric and asymmetric lattices non-perturbatively. We also study the colour field distributions in excited gluonic states.

Eta-deuteron 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 two-body 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.