
We discuss effective field theory treatments of the problem of three
particles interacting via shortrange forces. One case of such a system is
neutrondeuteron scattering at low energies. We demonstrate that in attractive
channels the renormalizationgroup evolution of the 1+2 scattering amplitude
may be complicated by the presence of eigenvalues greater than unity in the
kernel. We also show that these eigenvalues can be removed from the kernel by
one subtraction, resulting in an equation which is renormalizationgroup
invariant. A unique solution for 1+2 scattering phase shifts is then obtained.
We give an explicit demonstration of our procedure for both the case of three
spinless bosons and the case of the doublet channel in nd scattering. After the
contribution of the twobody effective range is included in the effective field
theory, it gives a good description of the nd doublet phase shifts below
deuteron breakup threshold.

We present a comprehensive analysis of deep inelastic scattering from He3
and H3, focusing in particular on the extraction of the free neutron structure
function, F_2^n. Nuclear corrections are shown to cancel to within 12% for the
isospinweighted ratio of He3 to H3 structure functions, which leads to more
than an order of magnitude improvement in the current uncertainty on the
neutron to proton ratio F_2^n/F_2^p at large x. Theoretical uncertainties
originating from the nuclear wave function, including possible nonnucleonic
components, are evaluated. Measurement of the He3 and H3 structure functions
will, in addition, determine the magnitude of the EMC effect in all A < 4
nuclei.

The BetheSalpeter (BS) amplitude for $\pi N$ scattering is evaluated at the
off mass shell points corresponding to the Low Energy Theorems (LET) based on
PCAC and current algebra. The results suggest a way of maintaining consistency
between BS equation and LET.

A covariant model of elastic pionnucleon scattering based on the
BetheSalpeter equation is presented. We obtain a good description of the S
and Pwave phase shifts up to 360 MeV laboratory energy. We also compare
results from the Kmatrix approach and several 3dimensional quasipotential
equations to the BetheSalpeter equation.

We discuss a possible measurement of the ratio of the nucleon structure
functions, F2n/F2p, and the ratio of the up to down quark distributions, u/d,
at large x, by performing deep inelastic electron scattering from the 3H and
3He mirror nuclei with the 11 GeV upgraded beam of Jefferson Lab. The
measurement is expected to be almost free of nuclear effects, which introduce a
significant uncertainty in the extraction of these two ratios from deep
inelastic scattering off the proton and deuteron. The results are expected to
test perturbative and nonperturbative mechanisms of spinflavor symmetry
breaking in the nucleon, and constrain the structure function parameterizations
needed for the interpretation of high energy ep, pp and ppbar collider data.
The precision of the expected data can also allow for testing competing
parameterizations of the nuclear EMC effect and provide valuable constraints on
models of its dynamical origin.

We study deep inelastic scattering on isospin asymmetric nuclei. In
particular, the difference of the nuclear structure functions and the Gottfried
sum rule for the lightest mirror nuclei, 3He and 3H, are investigated. It is
found that such systems can provide significant information on charge symmetry
breaking and flavor asymmetry in the nuclear medium. Furthermore, we propose a
new method to extract the neutron structure function from radioactive isotopes
far from the line of stability. We also discuss the flavor asymmetry in the
DrellYan process with isospin asymmetric nuclei.

In a recent paper Yang and Bodek extracted the free neutron structure
function at large x by extrapolating the density dependence of the nuclear EMC
ratios to the deuteron. We clarify why this approach is illdefined for light
nuclei, and introduces a large theoretical bias into the extraction of F_2n at
large x.

Dimensional regularization is applied to the LippmannSchwinger equation for
a separable potential which gives rise to logarithmic singularities in the Born
series. For this potential a subtraction at a fixed energy can be used to
renormalize the amplitude and produce a finite solution to the integral
equation for all energies. This can be done either algebraically or
numerically. In the latter case dimensional regularization can be implemented
by solving the integral equation in a lower number of dimensions, fixing the
potential strength, and computing the phase shifts, while taking the limit as
the number of dimensions approaches three. We demonstrate that these steps can
be carried out in a numerically stable way, and show that the results thereby
obtained agree with those found when the renormalization is performed
algebraically to four significant figures.

A relativistic description of pionnucleon scattering based on the
fourdimensional BetheSalpeter equation is presented. The kernel of the
equation consists of s and uchannel nucleon and delta pole diagrams, as well
as rho and sigma exchange in the tchannel. The BetheSalpeter equation is
solved by means of a Wick rotation, and good fits are obtained to the s and
pwave phase shifts up to 360 MeV pion laboratory energy. The coupling
constants determined by the fits are consistent with the commonly accepted
values in the literature.

To understand the difference between the value of the nn and LambdaLambda
twobody matrix elements, we examine the role of the coupling between the
LambdaLambda and XiN channels, first within the framework of SU(3), and
second in the numerical results for the binding energy of LambdaLambda 6He. We
find that it is essential to include the coupled channel nature of the BB
interaction as we proceed from strangeness zero to the strangeness 1 and 2
channels if we are to understand the relative magnitude of the twobody matrix
elements in the different strangeness channels.

A covariant model of elastic pionnucleon scattering based on the
BetheSalpeter equation is presented. The kernel consists of s and uchannel
nucleon and delta poles, along with rho and sigma exchange in the tchannel. A
good fit is obtained to the s and pwave phase shifts up to the twopion
production threshold.

The reaction $\Xi^d\to n\Lambda\Lambda$ is studied within the framework of
the Faddeev equations as a possible tool to gain insight into the final state
$\Lambda$$\Lambda$ interaction. The neutron differential energy spectrum
gives a final state interaction that is sensitive to both the
$\Lambda$$\Lambda$ amplitude at threshold, and the coupling between the
$\Lambda$$\Lambda$ and $\Xi$$N$ channels. The latter is a result of
interference between two mechanisms for the production of the final state,
which suggests that this reaction could give a measure of flavor SU(3)
violation in the twobaryon system.

The Nijmegen OBE potential D is SU(3) rotated to model the strangeness 2
sector of the baryonbaryon force. Soft core repulsion is introduced to
regularize the singular nature of the OBE functions. The strength of the
$^1$S$_0$ $t=0$ interaction is adjusted to model three scenarios: (1) a bound
state, (2) a narrow virtual, or antibound, state, and (3) an unbound state in
which the force is weakly attractive. Using a separable approximation to these
potentials, the binding energy of $^{6}_{\Lambda\Lambda}$He is calculated in an
$\Lambda\Lambda\alpha$ model. The resultant binding energies suggest that the
strength of the $^1$S$_0$ $t=0$ $\Lambda\Lambda$ and the $nn$ interactions
should be similar, if the coupling between the $\Lambda\Lambda$ and $\Xi N$
channels is taken into consideration.

The $\Lambda\Lambda$ interaction resulting from the SU(3) rotation of the
Swave component of the Nijmegen OBE potential $D$ is used to calculate the
binding energy of $^{ 6}_{\Lambda\Lambda}$He as a $\Lambda\Lambda\alpha$
threebody system, and the neutron differential energy spectrum for the
reaction $\Xi^ + d \to n + \Lambda + \Lambda$.

Equaltime Green's function is used to derive a threedimensional integral
equation from the BetheSalpeter equation. The resultant equation, in the
absence of antiparticles, is identical to the use of timeordered diagrams,
and has been used within the framework of $\phi^2\sigma$ coupling to study the
role of energy dependence and nonlocality when the twobody potential is the
sum of $\sigma$exchange and crossed $\sigma$ exchange. The results show that
nonlocality and energy dependence make a substantial contribution to both the
onshell and offshell amplitudes.

We derive a set of coupled fourdimensional integral equations for the
$NN\pi NN$ system using our modified version of the Taylor method of
classificationofdiagrams. These equations are covariant, obey two and
threebody unitarity and contain subtraction terms which eliminate the
doublecounting present in some previous fourdimensional $NN\pi NN$
equations. The equations are then recast into a from convenient for computation
by grouping the subtraction terms together and obtaining a set of twofragment
scattering equations for the amplitudes of interest.

By calculating the contribution of the $\pi\pi$ threebody force to the
threenucleon binding energy in terms of the $\pi N$ amplitude using
perturbation theory, we are able to determine the importance of the energy
dependence and the contribution of the different partial waves of the $\pi N$
amplitude to the threenucleon force. A separable representation of the
nonpole $\pi N$ amplitude allows us to write the threenucleon force in terms
of the amplitude for $NN\rightarrow NN^*$, propagation of the $NNN^*$ system,
and the amplitude for $NN^*\rightarrow NN$, with $N^*$ being the $\pi N$
quasiparticle amplitude in a given state. The division of the $\pi N$
amplitude into a pole and nonpole gives a procedure for the determination of
the $\pi NN$ form factor within the model. The total contribution of the
threebody force to the binding energy of the triton for the separable
approximation to the Paris nucleonnucleon potential (PEST) is found to be very
small mainly as a result of the energy dependence of the $\pi N$ amplitude, the
cancellation between the $S$ and $P$wave $\pi N$ amplitudes, and the soft
$\pi NN$ form factor.

By calculating the contribution of the $\pi\pi$ threebody force to the
threenucleon binding energy in terms of the $\pi N$ amplitude using
perturbation theory, we are able to determine the contribution of the different
$\pi N$ partial waves to the threenucleon force. The division of the $\pi N$
amplitude into a pole and nonpole gives a unique procedure for the
determination of the $\pi NN$ form factor in the model. The total contribution
of the threebody force to the binding energy of the triton is found to be very
small.

The derivation of scattering equations connecting the amplitudes obtained
from diagrammatic expansions is of interest in many branches of physics. One
method for deriving such equations is the classificationofdiagrams technique
of Taylor. However, as we shall explain in this paper, there are certain points
of Taylor's method which require clarification. Firstly, it is not clear
whether Taylor's original method is equivalent to the simpler
classificationofdiagrams scheme used by Thomas, Rinat, Afnan and Blankleider
(TRAB). Secondly, when the Taylor method is applied to certain problems in a
timedependent perturbation theory it leads to the overcounting of some
diagrams. This paper first restates Taylor's method, in the process uncovering
reasons why certain diagrams might be doublecounted in the Taylor method. It
then explores how far Taylor's method is equivalent to the simpler TRAB method.
Finally, it examines precisely why the doublecounting occurs in Taylor's
method, and derives corrections which compensate for this doublecounting.

This paper is the first in a series of three which attempt to resolve the
difficulties that have plagued the $NN\pi NN$ problem for the past ten years.
Various theoretical inconsistencies in the current formulation have been
pointed out and this work aims to eliminate these inconsistencies and so, we
hope, produce agreement with experiment. This is to be done by using covariant
perturbation theory, in which these inconsistencies are not present. The
covariant perturbation theory is developed starting from a model Lagrangian, in
order to fix notation and phases. It is shown that both oldfashioned
"timeordered" perturbation theory and the convolution integral of Kvinikhidze
and Blankleider may be recovered from the covariant perturbation theory when
certain approximations are made. The connection of these results with the work
of Klein, L\'evy, Macke and Kadyshevsky is discussed. Two forthcoming papers
will pursue this covariant calculation in the $NN\pi NN$ system using the
model and perturbation scheme developed in this paper and derive fully
covariant $NN\pi NN$ equations without the double counting problems present in
previous covariant equations.