relativistic heavy ion collision experiments at the bnl - ags ( au(10.7  agev)+au ) and at the cern - sps ( pb(160  agev)+pb ) have yielded a large variety of fascinating data .
various observables like the strong @xmath2 suppression , enhanced yield of intermediate mass dilepton pairs , enhanced ( anti-)hyperon yields , the creation of antimatter clusters and strong transverse flow seem to indicate the formation of very dense and highly excited matter @xcite . since these observables are connected in a non - trivial way it is a tempting task for theoreticians to model high energy heavy ion collisions in a consistent way and simultaneously predict this wide range of observables from a few hundred mev up to several thousand gev per nucleon at lhc .
bear in mind , however , that up to now there is no unique theoretical description of the underlying hadron@xmath3hadron interactions , with their vastly different characteristics at different incident energies and in different kinematic intervals .
perturbative quantum chromodynamics ( pqcd ) can be applied to describe hard processes , i.e. processes with large four - momentum , @xmath4 , transfer .
but pqcd is formally inappropriate for the description of the soft interactions because of the absence of the large @xmath5scale .
therefore , low@xmath6 collisions are described in terms of phenomenological models .
early on , multiple production of secondaries in relativistic hadronic collisions has been described within the hydrodynamic approach @xcite . then regge theory @xcite and
multiperipheral models have been developed to understand the phenomenology of the soft interactions .
they avoid the difficulties attributed to the statistical models .
an inconvenient point of this approach is the large number of free parameters , which have to be fixed by comparison to experiment .
subsequently , various qcd - motivated quark@xmath3parton models have been introduced .
consequently , a vast variety of models for hadronic- and nuclear collisions have been developed .
they may be subdivided into macroscopic ( statistical and hydrodynamical ) models @xcite and microscopic ( string- , transport- , cascade- , etc . )
models like , e.g. urqmd @xcite , which is applied in the present paper , fritiof @xcite , venus@xcite , qgsm @xcite , rqmd @xcite and others @xcite including the parton cascade approach @xcite . in the hydrodynamical ( thermal ) model
one assumes local ( global ) equilibrium - the dynamics is characterized by the equation of state employed .
the microscopic models describe subsequent individual hadron@xmath3hadron collisions .    for low and intermediate energies
hadron@xmath3hadron and nucleus@xmath3nucleus collisions are described in terms of the interactions between hadrons and their excited states , resonances , i.e. on the ( quasi-)particle level . at high energies the quark and gluon degrees of freedom can not be neglected .
then the concept of color string excitations is introduced with their subsequent fragmentation into hadrons . in lead@xmath3on@xmath3lead collisions at the full
sps energy one finds in the urqmd model that the ten most frequent hadron@xmath3hadron collision types namely , @xmath7 ( in decreasing order of frequency ) describe only 50% of the total h - h collisions .
the inclusion of an additional 120 h - h collision types allows modeling up to 90% of the collisions predicted in the urqmd model , while several thousand different h - h combinations are needed to cover more than 99% of the total number of h - h collisions .
since only a few of these cross sections are measured , one relies heavily on extrapolations ( and transformation , e.g. via the detailed balance principle ) of known processes .
therefore , here we want to present and analyze the detailed elementary h - h input used in the transport model urqmd .
the h - h predictions are a necessary basis for our understanding of the dynamics of the complex heavy ion reactions .
this paper is structured as follows : a brief description of the basic principles of the urqmd model is given in sec .
section  [ sec3 ] presents urqmd - results of the different h - h cross - sections for different reactions and a comparison with the available experimental data .
the additive quark model ( aqm ) is used to calculate unknown cross - sections .
urqmd is the first microscopic model which attempts to include the color coherent phenomena .
the implications of the effects of color opacity and color transparency in the model are discussed .
the treatment of the formation and decay of resonances and strings is described in detail in sec .
the importance of finite size effects in the fragmentation of strings is demonstrated .
section  [ sec5 ] discusses the generation of the transverse momentum of particles in the model . in sec .
[ sec6 ] several predictions of observables for elementary channels are presented , which are especially interesting for the upcoming proton - proton run of the na49 collaboration .
finally , a summary and conclusions are given .
the urqmd - model @xcite is a microscopic transport theory based on the covariant propagation of all hadrons on classical trajectories in combination with stochastic binary scatterings , color string formation and resonance decay .
it represents a monte carlo solution of a large set of coupled partial integro - differential equations for the time evolution of the various phase space densities @xmath8 of particle species @xmath9 etc . , which non - relativistically assumes the boltzmann form : + + = st f_i(x , p ) , where @xmath10 and @xmath11 are the position and momentum of the particle , respectively , and @xmath12 denotes the collision ( or rather source- ) term of these particle species , which are connected to any other particle species @xmath13 .    the exchange of electric and baryonic charge , strangeness and four momentum in the @xmath14-channel is considered for baryon - baryon ( bb ) collisions at low energies , while meson - baryon ( mb ) and meson - meson ( mm ) interactions are treated via the formation and decay of resonances , i.e. the @xmath15-channel reactions .
@xmath14-channel reactions for mb and mm collisions are taken into account from @xmath16  gev on increasing to the only mb , mm interaction type above @xmath17  gev . for nucleus - nucleus collisions
the soft binary and ternary interactions between nucleons can be described by the real part of the in - medium g - matrix , which is approximated by a non - relativistic density - dependent skyrme potential of the form v^sk = t_1 _ i j ( x_i - x_j ) + t_2 _ i j k ( x_i - x_j ) ( x_j - x_k ) , where @xmath18 denotes the coordinate variable in the quantum phase space . the first term simulates the attractive potential of the nn - interaction , and the second one yields the saturation .
in addition , yukawa and coulomb potentials are implemented in the model .
the potentials allow to calculate the equation of state of the interacting many body system , as long as it is dominated by nucleons .
note that these potential interactions are only used in the model for baryons / nucleons with relative momenta @xmath19 of less than 2@xmath20gev/@xmath21 .
for the hadronic collisions discussed here , the potential interactions are omitted .
further details of the application of the urqmd model to heavy - ion reactions may be found in @xcite .
this framework allows to bridge with one concise model the entire available range of energies from the sis energy region ( @xmath22gev ) to the rhic energy ( @xmath23gev ) . at the highest energies
, a huge number of different particle species can be produced .
the model should allow for subsequent rescatterings .
the collision term in the urqmd model includes more than fifty baryon species and five meson nonetts ( 45 mesons ) .
the baryons and baryon resonances included in the urqmd are listed in table [ tab1 ] .
in addition , their antiparticles have been implemented using charge - conjugation to assure full baryon - antibaryon symmetry .
figures [ fig1 ] and [ fig2 ] depict the implemented meson multiplets : pseudo - scalar , vector , scalar , pseudo - vector and ( not shown in the figs . ) the tensor mesons as well as the heavy vector meson resonances @xmath24 , @xmath25 , @xmath26 , and @xmath27 .
extremely heavy meson resonances ( @xmath28  gev ) are not explicitly implemented , however they may be important when investigating , e.g. the dynamics of @xmath29 correlations in future experiments .    all particles can be produced in hadron - hadron collisions and can interact further with each other .
the different decay channels all nucleon- , @xmath30- and hyperon - resonances up to 2.25@xmath20gev/@xmath31 mass as well as the meson ( e.g. k@xmath32 ) decays etc . are implemented . at higher energies
we take advantage of the hadron universality and use a string model for the decay of intermediate states .
the cross - sections of various hadronic processes as well as the formation and fragmentation of the strings are discussed in the subsequent chapters .
a basic input into the microscopic transport models are the particle species and -energy dependent cross - sections of hadron - hadron interactions .
the total cross - sections are interpreted geometrically .
a collision between two hadrons will occur if @xmath33 , where @xmath34 and @xmath35 are the impact parameter of the hadrons and the total cross - section of the two hadrons , respectively . in the urqmd model the total cross - section @xmath35 depends on the isospins of colliding particles , their flavor and the c.m .
however , partial cross - sections are then used to calculate the relative weights for the different channels .
only a small fraction of all possible hadronic cross - sections has been measured . in the following sections , we compare the urqmd cross - sections with experimental data .
if no data are available , the additive quark model and detailed balance arguments are used to extrapolate such unknown observables .      the total bb cross - section of the reaction @xmath36 has the general form _ tot^bb ( ) ( 2s_d + 1 ) ( 2s_e + 1 ) |m|^2 , with the spins of the particles , @xmath37 , momenta of the pairs of particles , @xmath38 , in the two - particle rest frame , and the matrix element @xmath39 . the matrix element @xmath39
, however , can take on a very complicated form and may be in general a function of of all the particle s quantum numbers as well as it s momenta and the c.m .
energy .
if high quality experimental data on the respective cross section exists , a phenomenological fit to the respective data is by far the most accurate approach for implementing the cross section .
otherwise , we have to rely on simplified assumption for the matrix element and employ general symmetries , like the principle of detailed balance ( see section  [ res_sec ] ) .
let us start by investigating the total cross - section of proton - proton collisions from a beam momentum of 0.1@xmath20gev/@xmath21 up to @xmath40gev/@xmath21 as shown in fig .
[ fig3 ] . the total and inelastic cross - sections of the @xmath41 reaction
are well measured in this energy region @xcite .
one finds a complex structure in this cross - section : local minima at 700@xmath20mev/@xmath21 ( @xmath42gev ) and 100@xmath20gev/@xmath21 ( @xmath43gev ) , the maximum at 2@xmath20gev/@xmath21 ( @xmath44gev ) and a rise above 100@xmath20gev/@xmath21 .
note that the steep rise in the data below 300@xmath20mev/@xmath21 is due to soft coulomb interaction of the protons and is taken care of via the potential interaction .
the structure in the @xmath41 cross - section is mainly due to the inelastic channels which are shown in fig .
one clearly sees the @xmath30 excitation with its increasing cross - section at low energies .
the different partial cross sections depicted in fig .
[ fig4 ] are discussed in section  [ res_sec ] .
a detailed comparison of prominent outgoing channels is depicted in fig .
[ ppc ] . here
we show calculations of exclusive ( @xmath45 ) and inclusive ( @xmath46 ) cross sections for the production of neutral mesons @xmath47 as a function of the excess energies @xmath48 . here
@xmath49 is the energy of the production threshold calculated as @xmath50 with the proton mass @xmath51 and the pole mass of the meson @xmath52 . in the case of @xmath53 mesons
we count only those with masses within @xmath54  mev around the pole mass to compare with data .
the exclusive @xmath55 production just above threshold @xcite is overestimated by a factor of about two .
note that above 3.5  gev the exclusive cross sections become less important because the string and multiple decay channels open and allow for multiple resonance production . for upcoming gsi - sis experiments
relevant @xmath56 values are below 4  gev , where so far no data on the inclusive channels are available and therefore rely heavily on extrapolations .    in fig .
[ pbar ] the cross section of @xmath57 reactions is shown .
good agreement with the data @xcite is found over a large energy range .    at higher energies ,
the contributions of the different nucleonic resonances decrease and give way to the excitation of color strings , which is the dominant process at high energies in our model .
the total cross section above the resonance region is given by the cern - hera parameterization @xcite as shown in table [ tab2 ] .
partonic pqcd scattering is not included into the urqmd model in the present version .
the difference between the total and the elastic cross - section is taken as the inelastic cross - section .
the mb cross - sections are dominated by the formation of @xmath15-channel resonances , i.e. the formation of a transient state of mass @xmath58 , containing the total c.m .
energy of the two incoming hadrons . on the quark level
such a process implies that a quark from the baryon annihilates an antiquark from the incoming meson . below 2.2 gev c.m .
energy intermediate resonance states get excited .
the total cross - section of these reactions are given by the expression : @xmath59 which depends on the total decay width @xmath60 , on the partial decay width @xmath61 and on the c.m .
energy @xmath56 . at higher energies
the quark - antiquark annihilation processes become less important .
there , @xmath14-channel excitations of the hadrons dominate , where the exchange of mesons and pomeron exchange determines the total cross - section of the mb interaction @xcite .
figures [ fig5 ] and [ fig6 ] show the cross - section of pion - proton reactions at different energies . in fig .
[ fig5 ] ( @xmath62 ) one probes predominantly the creation of the @xmath63 ( @xmath64 ) resonance .
note that the low energy s - wave @xmath65 scattering is not included into the urqmd fit .
the resonance peak at @xmath66  gev / c is from the @xmath67 resonances . in comparison fig .
[ fig6 ] ( @xmath68 ) depicts many strong uncharged non - strange baryon resonances , e.g. the @xmath69 etc .
the total cross - section in the intermediate energy regime is therefore the sum of the individual excitation modes of baryon resonances , the s - wave at lower energies is left out .
let us now investigate collisions of strange mesons with baryons . for @xmath70 mesons
strange @xmath15-channel resonances can be formed on non - strange baryons due to the annihilation of the @xmath71-quark .
a comparison of these processes from urqmd with the experimental data is shown in fig .
the formation of hyperon resonances is clearly visible at lower energies , while the universal @xmath14-channel reaction dominates the high energy tail .
[ fig8 ] shows the cross - section of @xmath72-mesons ( @xmath73 ) on protons . in this case , the formation of resonances is forbidden , since the @xmath74-quark can not be annihilated by non - strange baryons ( for strange baryons the formation of resonances is still possible ) .
here we use only the elastic channel and the @xmath14-channel excitation of both particles . the cross section at very high energies is given by the cern - hera parametrization as shown in table [ tab2 ] .      due to the fact that the experimental preparation of meson beams and targets is restricted to @xmath75 s and k s ,
only very little is known about mm collisions in general . for the description of heavy ion collisions the importance of this
channel increases with energy : at 1 a@xmath20gev beam energy we find that the production of new hadrons ( mostly pions ) is only a ten percent effect . at ags energies ( 10 a@xmath20gev )
the amount of mesons roughly equals the number of incoming nucleons . going on to the sps ( 160 a@xmath20gev )
the picture changes drastically : the produced particles dominate the reactions , while the incoming nucleons have dropped to a 15% admixture in particle density and multiplicity @xcite .    to describe the total meson - meson reaction cross - sections , we make use of the additive quark model ( see below ) and the principle of detailed balance , which assumes the reversibility of the particle interactions .    fig .
[ fig9 ] compares the calculated cross - section of @xmath76 scattering to experimental data @xcite .
the spectrum is dominated by the formation of the @xmath77 with a mass of 770  mev , the other two small peaks belong to the @xmath78 and @xmath79 resonances .
the @xmath78 resonance is not visible in the data since the experimental analysis of meson - meson scattering is model dependent .
figure  [ fig10 ] shows the implemented elastic @xmath80 cross - section , which remains constant in the whole energy range where data are available @xcite .
strangeness production in the meson - meson channel is possible , e.g. via the reaction @xmath81 as shown in fig .
[ fig11 ] .
rescattering of strange mesons is implemented via the resonance formation ( cf .
[ fig12 ] dominated by the @xmath82 resonance @xcite ) , or elastically as depicted in fig .
[ fig13 ] @xcite .
finally , we predict the cross - sections of @xmath83 , @xmath84 and @xmath85-reactions ( fig .
[ fig14 ] ) which are of utmost importance for the production of thermal photons and dileptons . at higher energies
other meson resonances can be formed .
to model mm interactions above the resonance region @xmath86  gev we use the rescaled total @xmath65 cross section : _
tot^mm(s > 1.7)=_tot^p(s ) .
this is justified , since at high energies the total cross section is given by quark counting . in the energy region from @xmath87  gev to @xmath88  gev @xmath15-channel interactions are taken into account , while from @xmath89  gev on @xmath14-channel excitation of both mesons becomes the mm interaction process of increasing importance in the model .
the cross section for high energetic reactions are taken from the aqm - rescaled @xmath90 paramatrization by the cern - hera group ( see table [ tab2 ] ) .
the physics of baryon@xmath3antibaryon interactions has been an area of much theoretical and experimental activity for a rather long period .
it is well - known that at energies @xmath91gev/@xmath21 an important contribution to the total interaction cross - section comes from the process of annihilation , where only mesons are left in the final state . though the earlier experiments on @xmath92annihilation revealed a number of differences from the non - annihilation channels , it is not clearly understood whether these differences arose simply due to the kinematic restrictions on the available phase space , or whether they are related to dynamical differences between the non - annihilation and annihilation mechanisms .
the experimental results obtained in @xcite by comparison of @xmath41 with non - annihilation @xmath93 interactions at 32 gev/@xmath21 support the conclusion of equivalence of @xmath41 and non - annihilation @xmath93 interaction processes .
still , the nature of the baryon annihilation is subject to theoretical discussions . in the framework of the quark models based on topological @xmath94 expansion ,
the annihilation process is associated with the annihilation of string junctions , i.e. the point where strings are connected , such that the baryons have a @xmath95-form . in this case three
@xmath96 strings are formed .
the theory also allows for diagrams where the string junctions and one or two of the valence quarks can annihilate , corresponding to the formation of two strings or one string .
other theories consider an annihilation mechanism without invoking the string junction hypothesis .
this intriguing question has yet to be clarified ( for review see , e.g. @xcite , and references herein ) .    to avoid the difficulties attributed to these theoretical approaches
, urqmd is adjusted to known experimental data .
the total @xmath97 cross - section is shown in fig .
[ fig15 ] , as well as the annihilation and the elastic cross - sections .
the urqmd parameterizations depicted by lines are taken from koch and dover @xcite : _ ann^|p
p = _ 0^n , with @xmath98 mb , @xmath99 , a = 50 mev and b = 0.6 .
the @xmath100 cross - section does not differ significantly from the @xmath97 cross - section @xcite , hence they are set equal in the urqmd model .    at higher energies ,
cern - hera parameterizations @xcite are used for the total and elastic channel : ^|p p_tot , el(p)=a + b p^n + cln^2(p ) + dln(p ) , with the laboratory momentum @xmath11 in gev/@xmath21 and the cross - section @xmath101 in mb .
the parameters of the fit are listed in table  [ tab2 ] .    below @xmath102  gev
/ c the following parameterization is used : _
tot(p ) = \ {    ll 75.0 + 43.1 p^-1 + 2.6 p^-2 -3.9 p & ; 0.3<p<5 + 271.6 & ; p < 0.3    .    _
el(p ) = \ {    ll 31.6 + 18.3 p^-1 -1.1 p^-2 - 3.8 p & ; 0.3<p < 5 + 78.6 & ; p < 0.3    .
the sum of annihilation and elastic cross - sections do not yield the total cross - section : = _ tot-_el-_ann in urqmd this difference @xmath103 is interpreted as the diffractive cross - section which describes the excitation at least one of the collision particles to a resonance or to a string via pomeron exchange .
the annihilation of baryon@xmath3antibaryon pairs proceeds in the urqmd model according to _ rearrangement _ diagrams . here
the formation of two @xmath104-strings of equal energies in the c.m .
system is assumed while the remaining constituent quarks are rearranged into newly produced hadrons .
the generalization of the @xmath93 cross - section towards all possible antibaryon@xmath3baryon collisions can be done in different ways :    1 .   the anti - baryon baryon cross - section at a given c.m .
energy @xmath56 is equal to the @xmath93 annihilation cross - section at the same @xmath56 : _
|bb|_=_|pp| _ 2 .
the anti - baryon baryon cross - section at a given relative momentum @xmath105 is equal to the @xmath93 annihilation cross - section at the same @xmath105 ( fig .
[ fig16 ] ) : _
|bb|_p_rel=_|pp|_p_rel    in the urqmd we have chosen the first parameterization . since @xmath106 , the annihilation cross - section drops rapidly with rising particle mass .
the different treatment this cross - section can lead to systematic shifts in the antibaryon distributions for massive systems .
therefore , anti - baryon production off nuclei may be used to solve this ambiguity .
unknown cross - sections are calculated on the basis of the additive quark model ( aqm ) @xcite , which assumes the existence of dressed valence quarks , interacting very weakly inside of the hadron . at the phenomenological level
, the aqm gives a correct quantitative and qualitative description of , e.g. , the asymmetry of c.m .
spectra of secondaries in meson@xmath3nucleon and photon@xmath3nucleon reactions .
the aqm has predicted also the important role of resonances for the multiple production of particles in hadronic interactions .
to apply the aqm to calculations of heavy ion reactions , one needs to know the cross - sections of the quark interactions , which can be evaluated from the quark masses .
then , the unknown total cross - section of the high energy reaction can be calculated assuming a 40% reduced @xmath15-quark cross - section ( compared to that of @xmath107- and @xmath34-quark ) .
the elastic cross - section is derived from regge theory @xcite : _ elastic = 0.039 _ total^ [ mb ] , where _ total = 40()^m_1+m_2 ( 1 - 0.4 ) ( 1 - 0.4 ) [ mb ] , where @xmath108 for particle @xmath109 being a meson ( a baryon ) and @xmath110 is the number of strange quarks in particle @xmath109 .
this formula results from the high energy reactions , therefore there is no difference between antiparticles and particles . for bb - reactions
no additional energy dependence is employed in collisions involving strange baryons .
non - strange baryon cross - sections are not treated via the additive quark model , they have an explicit energy dependence in line with experimental data .
the mb and mm cross - sections are rescaled via : _
the hyperon - nucleon cross - section , which is taken in the urqmd model from the additive quark model , is in good agreement with the data above @xmath111  mev ( cf .
[ fig17 ] ) .
the total cross - sections calculated for baryon - baryon , meson - baryon and meson - meson interactions are listed in tables  [ tab3]@xmath3[tab5 ] .
quantum chromodynamics has important applications of the dynamical role of color degrees of freedom to the strong interactions at ultrarelativistic energies ( for a review , see @xcite and references therein ) .
the theory is presented in ref .
@xcite in detail , here we just sketch the main ideas of color optics and coherent phenomena in high energy physics .
hadrons are dynamical objects which come in fock space of configurations of very different spatial sizes . at high energies , incident
hadronic quark@xmath3gluon configurations can be considered frozen as a result of lorentz time dilation . due to the long coherence length at such high energies one can apply geometrical color optics .
small objects produced in hard processes with high momentum transfer @xmath4 then have reduced interaction cross sections . in processes with moderate @xmath4
such compact objects , which are a coherent superposition of eigenstates of the qcd hamiltonian , should gain size .
when the quark@xmath3gluon configuration is large , it will lead to an increased interaction cross - section of the hadron with the nuclear medium .
therefore , the fluctuations of the hadron s spatial extent give rise to the color transparency and color opacity phenomena : when a small object is produced , it interacts only very weakly with other hadrons due to color screening .
moreover , since - at sufficiently high energies - the small - sized configuration of this object is frozen , the nuclear medium appears to be transparent for such hadrons ( color transparency ) .
in contrast , hadronic configurations which are larger than average interact with larger cross - section , giving rise to color opacity .
nucleus - nucleus collisions provide a tool to investigate the effect of color transparency , for instance , in the production of leading nucleons . on the other hand , the complementary color opacity effect , i.e. large - sized configurations , can cause stronger stopping and significant fluctuations in the transverse energy of secondaries in central reactions .    a first step to investigate these qcd effects within a microscopic transport model
is made by incorporating the color fluctuations in the elementary hadron@xmath3hadron reactions in the urqmd model .
thus , one needs to know the probability @xmath112 that a given configuration interacts with a nucleon with a total cross - section @xmath101 .
it is convenient to consider moments of the distribution : ^0 & = & |^0 p ( ) = 1 , + ^1 & = & |^1 p ( ) =| , etc . , where @xmath113 denotes the average cross - section .
the second moment @xmath114 can be determined from the diffractive dissociation experiments .
in addition , further information can be obtained from qcd , which implies @xcite : p()^n_q-2 , for @xmath115 , where @xmath116 is the number of valence quarks .
thus , for the nucleon and the pion distributions it follows for @xmath115 : p_n ( ) & & , + p
_ ( ) & & constant
. from these arguments , @xmath112 can be construct .
[ fig18 ] shows the resulting broad @xmath112 distribution for proton projectiles and the even broader one for the pions @xcite .
the effect of the color fluctuations on proton - proton interactions at different impact parameters , @xmath117 , is shown in fig .
[ fig19 ] for the urqmd model .
the charged pion multiplicity distribution decreases monotonically with rising @xmath117 .
@xmath118 has a non - vanishing tail for @xmath119fm , in contrast to the abrupt geometrical edge of the distribution as calculated in the static geometric models without color fluctuations .
the production and decay of resonances is the most important h - h reaction channel below @xmath120  gev for bb and 3  gev for mm and mb reactions .
baryon resonances are produced in two different ways , namely    * _ hard production _ : n+n@xmath121n ,  @xmath122 ,  n@xmath32n , etc . *
_ soft production _ : @xmath123p@xmath124 ,  k@xmath125+p@xmath126 ...    the formation of @xmath15-channel resonances is fitted to measured data , e.g. in the reaction @xmath36 we use the general form _
tot^bb ( ) ( 2s_d + 1 ) ( 2s_e + 1 ) |m(m_d , m_e)|^2 , with the spins of the particles , @xmath37 , momenta of the pairs of particles , @xmath38 , in the two - particle rest frame , and the matrix element @xmath127 , which here depends only on the masses of the outgoing hadrons , @xmath128 .
there are six channels of the excitation of non - strange resonances in the urqmd model , namely @xmath129 .
the @xmath130 is explicitly listed , whereas higher excitations of the @xmath30 resonance have been denoted as @xmath131 . for each of these 6 channels specific assumptions
are made with respect to the form of the matrix element @xmath132 , and the free parameters are adjusted to the available experimental data , when available :    1 .
@xmath133 excitation : @xmath134 with @xmath135 gev , @xmath136 gev and @xmath137 .
note that this form of the matrix element has been adjusted to fit the data shown in figure [ pp - nd1232 ] .
2 .   @xmath138 , @xmath139 , @xmath140 and @xmath141 excitation : @xmath142 with @xmath143 for @xmath138 , @xmath144 for @xmath139 and @xmath145 for @xmath140 . in the case of @xmath141
there are insufficient data available , therefore we use the same matrix element and parameters as in the case of @xmath140 . since @xmath146 is valid for all above cases , the matrix element can not diverge .
@xmath147 excitation : @xmath148 with @xmath149 .
figure  [ pp - nd1232 ] shows the fit of the urqmd @xmath150 cross section to experimental measurements @xcite .
the measurements refer to the @xmath151 exit channel and have been rescaled to match the full isospin - summed cross section . in the case of the exclusive @xmath130 cross section the quality of the data and thus also
the quality of the resulting fit is very good . for higher resonance excitations
this is unfortunately no longer the case and additional measurements are needed to clarify the situation .
one has to keep in mind , however , that the experimental extraction of exclusive resonance production cross sections is only possible via two- or three - particle correlations ( e.g. a pion - nucleon correlation in the case of the @xmath30 ) which introduces large systematic errors , especially for broad resonances .    in figure  [ pp - ppstar ] the urqmd cross sections for the processes @xmath152 , @xmath153 @xmath154 and @xmath155
are compared to data @xcite .
one single parameter has been used to fit all four cross sections .
the data situation is not as good as in the case of the @xmath130 resonance , some ambiguities are visible which results in the quality of the fit being not as good as in the previous case .
the parameters for the other classes are fitted in the same fashion .
the cross section for exclusive @xmath156 can be seen in figure  [ pp-2d1232 ] .
the data points @xcite hint at a resonance like structure which can not be reproduced with the urqmd ansatz for resonance - excitation cross sections .
however , the data deviates considerably from other cross sections for resonance excitation ( e.g. @xmath157 ) . considering hadron - universality and the similarities in all other resonance excitation cross sections
this casts a certain doubt on the accuracy of the measurement of the resonance - like peak .
figure  [ pp - nd1920 ] shows the urqmd fit for the exclusive @xmath158 production .
the same matrix element is used for the entire class of @xmath159 reactions . in the case of exclusive @xmath160 production
the matrix element has been fitted to the @xmath161 exit channel ( see figure  [ pp - d1232n1680 ] ) .
the extrapolation to the case of @xmath162 production can be seen in figure  [ pp - d1232n1520 ] . for the exclusive @xmath163 reaction class
the data situation is unsatisfactory , therefore we used the same matrix - element as in the @xmath164 case .    the decay of the resonances proceeds according to the branching ratios compiled by the particle data group @xcite .
the resonance decay products have isotropical distributions in the rest frame of the resonance .
if the resonance decays into @xmath165 particles , then the corresponding @xmath166body phase space is used to calculate their @xmath167 momenta stochastically .
it is necessary to note that a consistent description of angular momentum distributions points to a rather intricate problem of transport theory itself : if one considers the whole scattering interaction to be described by one single quantum mechanical process there are correlations between the final and the initial stage .
for instance the angular distribution of the final particles with respect to the axes of the incoming momenta in the cms system . however , a fitting of the angular distributions to experimental data may conflict with the basic assumption of transport theories that the multiple scattering processes can be considered to be of markovian type , i.e. after each scattering process or resonance formation the outgoing particles completely forget about their entrance channels . in the case of a spin 0 resonance there is no preferred direction for the emission of the final particles , while for spin 1 ( and other ) the different magnetic quantum numbers are statistically occupied , so that also in these cases there is no preferred angle of emission .
all produced particles are able to rescatter within the nuclear medium , therefore the excitation of resonances by the annihilation of mesons on baryons included as depicted in fig .
[ fig5 ] for the reaction @xmath168 .    also the @xmath169 channel ( fig .
[ fig6 ] ) shows a rich structure of baryon resonances .
the total meson - baryon cross section is given by formula  [ mbbreitwig ] .
there , the total and partial decay widths also define the inverse reaction , i.e. the different decay - channels of the respective resonance .
thus , the principle of detailed balance is applied .
based on this principle we calculate all resonance formation cross sections from the measured decay properties of the possible resonances up to c.m .
energies of 2.25  gev/@xmath31 for baryon resonance and 1.7  gev/@xmath31 in the case of mm and mb reactions . above this
energy collisions are modeled by the formation of @xmath15-channel string or , at higher energies ( beginning at @xmath170  gev ) , by one / two @xmath14-channel strings . in the strangeness channel
elastic collisions are possible for those meson - baryon combinations which are not able to form a resonance , while the creation of @xmath14-channel strings is always possible at sufficiently large energies ( c.f . fig .
[ fig7 ] for the formation of hyperon resonances and fig .
[ fig8 ] for the non - resonant channel ) . at high collision energies
both cross section become equal due to quark counting rules .    in more general terms , the principle of detailed balance can be derived by assuming time - reversal invariance of the interaction hamiltonian and can be formulated in the following way : [ detbal ] ( yx ) p_y^2 g_y = ( xy ) p_x^2 g_x , with @xmath171 denoting the c.m .- momenta of the particles and @xmath172 being the spin - isospin degeneracy factors .
thus , if the cross section of the reaction @xmath173 is known , the back reaction @xmath174 can be easily obtained .
this principle is in urqmd widely applied for the calculation of baryon - resonance absorption cross sections , such as @xmath175 . for a detailed discussion of the application of the principle of detailed balance to resonance absorption and @xmath56-dependent decay widths we refer to  @xcite .    as was mentioned above ,
not only baryon - baryon and meson - baryon collisions have to be included in the proposed scheme . at high energies and in massive aa systems meson - meson collisions
may dominate the multiple production of secondaries . unfortunately , there are only few channels for which the experimental information exists , like the process of @xmath176 scattering ( fig .
[ fig9 ] ) , which is fairly described by the urqmd model .      both ,
for the high energy regime and for baryon@xmath3antibaryon annihilation we apply a string model ( similar to the lund model @xcite ) to describe the inelastic reactions .
the constituents , quark and diquark ( or anti - quarks ) , of the incoming hadron also define the predominant emission patterns of the events . the amount of stopping in nucleus - nucleus is strongly correlated to the detailed dynamics of the diquark in the hadronic medium @xcite .
recently different additional mechanisms of baryon number transport in nucleus - nucleus reactions have been investigated : baryonic junctions as suggested by @xcite which yield an enhanced hyperon and proton production cross section at central rapidities @xcite ; di - quark breaking due to interaction with the hadronic medium as predicted by @xcite .
the di - quark breaking component is also taken into account in the urqmd model ( di - quark breaking probability@xmath177 ) .
however , the above mechanisms are of minor importance in the urqmd approach since rescattering of the leading di - quark with hadrons is explicitly taken into account .
since gluons are massless particles with spin @xmath178 , the static strong interaction between quarks at small distances @xmath179 may be described by a potential @xmath180 . at sufficiently large distances the color field between two quarks or anti
- quarks transforms into the color string . denoting the string tension as @xmath181 one defines the linear string potential @xmath182 between the quarks / diquarks located at @xmath183 and @xmath184 , respectively
this form of the potential is chosen from heuristic considerations , based on the quark confinement hypothesis , and is supported by lattice qcd calculations @xcite .
the transverse directions have not to be taken into account , since they are negligibly small compared to the longitudinal excitation of the hadron string . hence
, we get the dynamics of the quark system ( with quark momenta @xmath185 and @xmath186 ) from the hamiltonian @xmath187 @xmath188 which leads to the following equations of motion for the massless endpoints of the string : @xmath189 a change in momentum is directly related to the sign of @xmath190 , while the direction of propagation is defined by the sign of the momentum @xmath191 of the quark .
this results in a typical  yo - yo  type evolution of the quark system .
if the momentum transfer is large enough , the excitation of the string may exceed some critical limit .
after that it will be energetically favorable to break the string into pieces by producing @xmath104-pairs from the vacuum .
each of the produced @xmath104-pairs will have small relative momenta in their rest frame .
owing to the fact that the color string is uniformly stretched , the hadrons produced as a result of the string fragmentation will be uniformly distributed within the kinematically allowed interval between @xmath192 and @xmath193 .
the probability of the pair production process has been calculated by schwinger for the case of an infinite homogeneous electrical field .
his result can be used to motivate the decay of qcd color field between the quarks .
note that qcd is a non - abelian theory , therefore the color field need not be homogeneous and it is definitely not infinite .
this leads to the modifications of string decays which will be discussed later .
the probability @xmath194 for the creation of a quark - antiquark pair with mass @xmath195 in a color field with a string tension @xmath181 is : @xmath196 where a typical value for @xmath181 is 1  gev / fm .
this relation is motivated by schwinger s qed - based result for particle - antiparticle creation in a strong electric field ( see the discussion of formula  [ schw ] in section  [ finite ] ) .
the relative production probabilities of the different quark flavors are adjusted to @xmath197-data : @xmath198 the production of charmed ( and heavier ) quarks is strongly suppressed in the string picture , hence they are exclusively produced in hard qcd processes .
the elementary diquark is introduced to allow for baryon - antibaryon production in the string .
a schematic view of the decaying string is shown in fig .  [ fig20 ] @xmath3 a non - strange baryon decays into a hyperon , a kaon and a pion .    to decide which type of hadron is produced from the quark configuration that is created in the color field we choose in the case of a produced    ( * i * ) _ baryon _ @xmath3 the octet and decuplet with equal probabilities ;    ( * ii * ) _ meson _ @xmath3 the meson nonet with a probability proportional to the spin degeneracy and inverse mean mass @xmath195 , @xmath199 the singlet states are mixtures of @xmath200 , @xmath201 and @xmath202 .
they are projected onto su@xmath203 hadrons with the flavor mixing angles from the quadratic gell - mann - okubo mass formula @xcite . for the scalar mesons
this formula is not applicable , here an ideal mixing angle ( @xmath204 ) is assumed ( the mixing angles are depicted in table [ mixang ] ) .    the field@xmath3feynman fragmentation mechanism @xcite , which allows the independent string decay from both ends of the string is used in the urqmd model .
the string break - up is treated iteratively : string @xmath205 hadron + smaller string .
the conservation laws are fulfilled .
the diquark is permitted to convert into mesons via the breaking of the diquark link , thus transporting the baryon number into central rapidities .    on both sides of the fragmenting string
the new particles are formed randomly .
if a resonance is produced , its mass is determined by a breit - wigner mass distribution .
the transverse momentum is assigned to this particle according to a thermal momentum distribution , resulting in a temperature of 170  mev .
after that the fragmentation function determines the fraction of the longitudinal momentum of the string transferred to the hadron .
this procedure can be described in a covariant manner by the light cone variables defined as : @xmath206 the light cone momentum @xmath207 given to the newly produced hadron is : @xmath208 the fragmentation of a baryonic string reads : @xmath209    the main input is the fragmentation function which yields the probability distribution @xmath210 .
this function regulates the fraction of energy and momentum given to the produced hadron in the stochastic fragmentation of the color string . for newly produced particles the field - feynman function @xcite :
@xmath211 is used .
@xmath212 drops rapidly with increasing @xmath213 ( fig .
[ fig21 ] ) .
therefore , the longitudinal momenta of e.g. produced antibaryons ( fig .  [ fig22 ] ) and pions ( fig .
[ fig23 ] ) are small ( they stick to central rapidities ) , in line with the experimental data .
the rapidity spectra of these particles have a characteristic gaussian - like shape , in contrast to the baryon spectra in pp , as it is clearly seen in figure [ fig22 ] .
the proton is on average less stopped , since it is build up from the leading diquark in the string ( leading particle effect ) .
[ fig24 ] compares the @xmath214 distribution of protons and @xmath215 s for the feynman scaling variable @xmath216 measured in pp reactions at 205  gev/@xmath21 .
the data on leading baryons can only be reproduced when a modified fragmentation function is used for the leading baryons ( cf .
[ fig21 ] , dashed curve ) .
this leading baryon fragmentation function is of gaussian form : @xmath217 } \quad,\ ] ] with parameters @xmath218 and @xmath219 .
it is obvious that modeling the momentum loss in elementary collisions has a strong influence on the rapidity spectra of particles produced in heavy - ion collisions .
this can be seen if one compares the spectra of p , @xmath215 s ( created from the leading baryon ) and @xmath220 s or mesons ( created from newly produced quarks ) in lead - lead collisions at the sps energies @xcite .
let us now discuss finite size effects in the process of string fragmentation .
the string is essentially a color field which connects two color charges , the @xmath221 $ ] quark and the @xmath222 $ ] di - quark or anti - quark , at the ends of the string .
quantum electrodynamics ( qed ) predicts spontaneous particle@xmath3antiparticle creation in strong electric fields @xcite .
this effect should hold for particle creation in a strong color fields .
schwinger s qed result w_=^_n=1 [ schw ] is often adopted to the case of color fields by equating @xmath223 with the string tension @xmath181 .
however , for strings , several important assumptions , which lead to schwinger s result are not fulfilled : firstly , the color field is not infinitely extended .
it is bound radially by the interaction length of the gluons and longitudinally roughly by the @xmath221 $ ] and @xmath222 $ ] endpoints .
secondly , the two endpoints of the string move ( with close to the speed of light ) in opposite directions . finally , the assumption of a constant field strength @xmath181 seems to be fulfilled along the string , but may become invalid in hot and dense matter @xcite .    the influence of the finite radial size of the string can be studied by solving the dirac equation for the ( newly produced ) pairs in a finite volume . the field is restricted to a cylindrical volume of length @xmath224 and and area @xmath225 , where @xmath226 is the cylinder radius . along the cylinder
axis a homogeneous color field is assumed .
the boundary condition on the surface of the cylinder leads to enclosure of the color charges .
taking @xmath227 neglects the longitudinal direction .
the constraint on the cylinder surface is given by linear mit - boundary conditions @xcite , thus we have to solve the equations : @xmath228 the mit model allows for an analytic solution@xcite . the occupation numbers in the limits @xmath229
yield the pair creation probability as @xcite w_r=^_n=1 \ { _ n , } .
thus the string radius is directly related to the pair creation rate @xmath230 , with @xmath231 being the n@xmath232 momentum eigenvector in the solution of eqs .
( [ dirac1 ] and [ dirac2 ] ) for a given projection @xmath233 of the corresponding bessel functions .
the expression in curly brackets gives the deviation from the schwinger formula .
for @xmath234 this second factor vanishes and one is back to the infinite case .
this result allows to calculate the strangeness suppression @xmath235 for different string radii as shown in fig .
[ fig26 ] .
the overall pair creation rate has been calculated as a function of the longitudinal size of the color field for the corresponding dirac equation @xcite in terms of confluent hypergeometric functions .
the resulting particle production rate per volume d@xmath236 and time interval d@xmath14 , @xmath237 at @xmath238 is shown in fig .
[ fig27 ] as a function of @xmath213 , the distance from the string center , for different string lengths @xmath224 . a clear depletion of particle production near the endpoints of the string endings is visible ( eq.[schw ] ) . for short strings
the particle production gets enhanced toward the string center .
the above discussed finite size effects have a strong influence on the results , especially on the production of heavy quarks and particles like anti - baryons , -hyperons , etc . up to now
these corrections are included only in the non - strange antibaryon sector .
this effect is very important especially in the case of the @xmath215-particle , since it is mostly created at the string ends in the fragmentation of a leading proton .
indeed , it has been reported that most of the transport models which use a string fragmentation scheme based on the schwinger formalism tend to overestimate the @xmath215 s @xcite .
therefore , this finite size effects shall be included for strange baryons , in order to correct the overestimated number of @xmath215 s and @xmath239 s .
in hot and dense nuclear matter , most hadrons suffer interactions from the many surrounding particles . as a result ,
the effective mass may change with density .
many dynamical properties of hadrons are modified in the medium .
in - medium two - body scattering cross - sections may differ from the free space values .
those effects can be studied in the framework of relativistic transport theory , i.e. the relativistic boltzmann - uehling - uhlenbeck ( rbuu ) equation .
this type of transport equation has been used extensively to the study of relativistic heavy - ion collisions and turned out to be very successful .
the urqmd approach uses an analytical expression for the differential cross - section of in - medium @xmath240 elastic scattering derived from the collision term of the rbuu equation @xcite to determine the scattering angles between the outgoing particles in elementary hadron - hadron collisions .
it is assumed that the angular distributions for all relevant two - body processes are similar modified in an analogous manner .
they are approximated by the differential in - medium @xmath240 elastic scattering cross - section : @xmath241 with the direct term @xmath242 and the exchange term @xmath243 the ( pseudo-)scalar and vector coupling constants are @xmath244 , @xmath245 , and @xmath246 and @xmath247 is the in - medium mass , @xmath248 are the mandelstam variables .
the in - medium single - particle energy is given by @xmath249 the formula for the differential cross section of in - medium @xmath240 elastic scattering can be used for elementary hadron - hadron collision if it is scaled by @xmath250 where @xmath251 and @xmath252 denote the effective masses of the incoming hadrons .
furthermore , we take into account finite size effects of the hadrons and part of the short range correlation by introducing a phenomenological form factor at each vertex . for the baryon - baryon - meson vertex
the common form @xmath253 is used , where @xmath254 is the cut - off mass for meson @xmath255 .
+ the total energy and the masses of the incoming hadrons serve as input for calculating the angular distribution .
it is worth to stress again that eq.([angeq ] ) is used to calculate only the angular distributions for all elementary elastic two - body processes but not the corresponding total cross sections . the inverse slopes ( temperatures ) as calculated in the urqmd from the transverse momentum spectra of pions - by fitting the @xmath256 distribution with an exponential - compare well with values extracted from thermal fits to data , @xmath257 and @xmath197 . ]
@xcite ( fig .
[ fig28 ] ) .
the urqmd model reproduces nicely the total , elastic and inelastic cross - sections of numerous hadronic reactions .
the model also predicts the particle multiplicities ( i.e. the inclusive cross - sections ) as well as the ( lorentz - invariant ) cross - sections , which may come in the form of @xmath214- , rapidity- , or transverse momentum distributions .
the abundances of the most important particle species produced in @xmath41 collisions at 12 gev/@xmath21 @xcite are listed in table  [ tab12gev ] - the model predictions are in line with the data within 15% .
the yields of various particle species produced in @xmath41 collisions at 205 gev/@xmath21 ( @xmath258 gev ) @xcite are listed in table  [ tab6 ] .
it is easy to see that the model predictions lie generally well within the 10% of the data except for the strange baryons .
table  [ tab7 ] presents the results on particle production in @xmath41 interactions at @xmath259 gev @xcite .
again , the agreement between the experimental and theoretical results is good .
as discussed above , the production of @xmath215(@xmath220 ) s is overestimated by a factor of 2 - 3 due to the neglect of finite size effects in the strings .
the rapidity spectra of pions and @xmath214-distributions of baryons , as obtained from the urqmd model for @xmath41 interactions at 205 gev/@xmath21 ( fig .
[ fig22 ] and [ fig23 ] , respectively ) have been discussed already in the previous section .
correlations between the transverse and longitudinal momenta of charged pions , produced in the same reaction , can be studied @xcite .
the transverse momenta of both positively and negatively charged pions have evident minima at @xmath260 as shown in fig .
[ fig29 ] ( the  sea - gull
 effect ) . then the transverse momentum increase nearly linearly with the rising longitudinal momenta .
the agreement with the experimental data is fair .
the single event correlation between the transverse momenta of the produced pions and the multiplicity of the negatively charged hadrons in the event is presented in fig .
[ fig30 ] .
the measured distribution is reproduced nicely for @xmath261 s , while for the @xmath262 s , 20% deviations of the calculated spectrum from the data are evident at high ( @xmath263 ) event multiplicities .    to probe the ability of the urqmd model to describe hadron - hadron collision even at the today highest available bombarding energies for nuclei , we compare the calculated he+he collision at isr with data @xcite as shown in fig .
[ hehe ] .
it is not surprising , that the light helium system is transparent .
a baryon free area of 3 units in rapidity is produced .
the urqmd model prediction describes the data fairly well .
the calculated particle yields can be increased by 15% if one includes multiple jet creation into the model description @xcite .
the urqmd model reproduces the cross - sections and spectra of particles in hadronic collisions fairly well .
since hadronic interactions build up the basic input to simulate the hadron - nucleus and nucleus - nucleus interactions in the model , it is interesting to examine the applicability of the urqmd model to these reactions .
the full comparison with the experimental @xmath264 and @xmath265 data is an ongoing program , which is not completed yet @xcite .
the ability of the model to reproduce e.g. dilepton yields in @xmath41 collisions , which is of interest in high energy physics is shown here as an example .
[ fig31 ] shows the urqmd calculations of the dilepton spectrum for @xmath11+be ( which serves as a substitute for pp reactions ) at 450 gev/@xmath21 .
dilepton sources considered here are dalitz decays ( @xmath266 , @xmath55 and @xmath267 ) and direct vector meson decays ( @xmath77 , @xmath267 and @xmath268 ) .
dalitz decays of heavier meson and baryon resonances are included explicitly via their emission of @xmath77 mesons ( assuming vector meson dominance ) . to avoid double counting ,
the @xmath77 mesons from @xmath55 s , and @xmath267 s are excluded from the @xmath77 contribution .
pion annihilation is included dynamically into the contribution of decaying @xmath77 mesons via the channel @xmath269 .
the calculation of dilepton yields without modifications of the @xmath77 mass pole is compatible with the ceres data @xcite .
the basic hadronic interaction processes incorporated into the urqmd approach are given .
the implemented cross - sections of various hadron - hadron collisions , as well as their extrapolations into the high energy region are presented .
the model treatment of excitation and decay of intermediate objects like resonances and strings is described .
the importance of the finite size effects for the process of string fragmentation is discussed .
the urqmd model is a microscopic transport model which allows to study optionally color fluctuations , i.e. color coherence phenomena , as color transparency and color opacity , as well as the expansion of small wave packets from the point of the production .
the model predictions are compared with the available experimental data on particle yields in hadronic interactions for a wide range of c.m .
the urqmd model treats the elementary processes reasonably well .
more accurate data on proton - proton , proton - neutron , as well as meson - baryon collisions are needed to improve the extrapolation to nucleus - nucleus interactions at high energies , where already in the present model - about @xmath270 elementary hadron - hadron reactions are possible .
this work was supported by the graduiertenkolleg fr theoretische und experimentelle schwerionenphysik , frankfurt@xmath3giessen , the bundesministerium fr bildung und forschung , the gesellschaft fr schwerionenforschung , darmstadt , deutsche forschungsgemeinschaft , the alexander von humboldt@xmath3stiftung , the josef buchmann - stiftung , daad and in part by doe grant de - fg02 - 96er40945 .
quark matter96 , heidelberg , germany , nucl .
phys . a * 610 * ( 1996 ) ; quark matter97 , tsukuba , japan , nucl .
a * 638 * ( 1998 ) ; s.  a.  bass , m. gyulassy , h.  stcker , w.  greiner , j. phys .
g ( 1998 ) in print ; strange quark matter 97 , j. phys .
g * 23 * ( 1997 ) ; strange quark matter 98 , j. phys .
g ( 1998 ) .
landau , e.m .
lifshitz , fluid mechanics ( pergamon , oxford , 1959 ) ; h.  stcker , w.  greiner , phys . rep . * 137 * , 277 ( 1986 ) ; r.b .
clare , d.  strottman , phys .
rep . * 141 * , 177 ( 1986 ) ; e.  schnedermann , j.  sollfrank , u.  heinz , phys .
c * 48 * , 2462 ( 1993 ) ; j.  brachmann , a.  dumitru , j.  maruhn , h.  stcker , w.  greiner , d.  rischke , nucl .
a * 619 * , 391 ( 1997 ) .
bass , m.  belkacem , m.  bleicher , m.  brandstetter , l.  bravina , c.  ernst , l.  gerland , m.  hofmann , s.  hofmann , j.  konopka , g.  mao , l.  neise , s.  soff , c.  spieles , h.  weber , l.a .
winckelmann , h.  stcker , w.  greiner , ch .
hartnack , j.  aichelin , n.  amelin , prog .
. phys . * 41 * , 225 ( 1998 ) .
zabrodin , l.v .
bravina , o.l .
kodolova , n.a .
kruglov , a.s .
proskuryakov , l.i .
sarycheva , m.y .
bogolyubsky , m.s .
levitsky , v.v .
maximov , a.a .
minaenko , a.m.  moiseev , s.v .
chekulaev , phys .
d * 52 * , 1316 ( 1995 ) .
t.  kafka , r.  engelmann , r.  godbole , j.  hanlon , m.  pratap , h.  wahl , r.  singer , y.  cho , t.  fields , l.g .
hyman , l.  voyvodic , r.  walker , j.  whitmore , phys .
d * 16 * , 1261 ( 1977 ) ; g.  charlton et al .
. lett . * 30 * , 574 ( 1973 ) .
( 75,40)(0,0 ) ( 3,16)@xmath262 ( 17,16)@xmath266 ( 26,16)@xmath261 ( 17,13)@xmath271 ( 11,13)@xmath55 ( 8,27)@xmath272 ( 22,27)@xmath72 ( 8,2)@xmath273 ( 22,2)@xmath274 ( 0,15)(1,0)38 ( 15,0)(0,1)38 ( 5,15)(1,2)5 ( 10,25)(1,0)10 ( 20,25)(1,-2)5 ( 25,15)(-1,-2)5 ( 20,5)(-1,0)10 ( 10,5)(-1,2)5 ( 10,16)(0,-1)2 ( 20,16)(0,-1)2 ( 15,15 ) ( 15,15 ) ( 15,15 ) ( 5,15 ) ( 10,25 ) ( 20,25 ) ( 25,15 ) ( 20,5 ) ( 10,5 ) ( 34,13)@xmath275 ( 18,32)@xmath95 ( 43,16)@xmath276 ( 57,16)@xmath277 ( 66,16)@xmath278 ( 57,13)@xmath279 ( 51,13)@xmath267 ( 48,27)@xmath280 ( 62,27)@xmath281 ( 48,2)@xmath282 ( 62,2)@xmath283 ( 40,15)(1,0)38 ( 55,0)(0,1)38 ( 45,15)(1,2)5 ( 50,25)(1,0)10 ( 60,25)(1,-2)5 ( 65,15)(-1,-2)5 ( 60,5)(-1,0)10 ( 50,5)(-1,2)5 ( 50,16)(0,-1)2 ( 60,16)(0,-1)2 ( 55,15 ) ( 55,15 ) ( 55,15 ) ( 45,15 ) ( 50,25 ) ( 60,25 ) ( 65,15 ) ( 60,5 ) ( 50,5 ) ( 74,13)@xmath275 ( 58,32)@xmath95    ( 75,40)(0,0 ) ( 3,16)@xmath284 ( 17,16)@xmath285 ( 26,16)@xmath286 ( 17,13)@xmath287 ( 11,13)@xmath288 ( 8,27)@xmath289 ( 22,27)@xmath290 ( 8,2)@xmath291 ( 22,2)@xmath292 ( 0,15)(1,0)38 ( 15,0)(0,1)38 ( 5,15)(1,2)5 ( 10,25)(1,0)10 ( 20,25)(1,-2)5 ( 25,15)(-1,-2)5 ( 20,5)(-1,0)10 ( 10,5)(-1,2)5 ( 10,16)(0,-1)2 ( 20,16)(0,-1)2 ( 15,15 ) ( 15,15 ) ( 15,15 ) ( 5,15 ) ( 10,25 ) ( 20,25 ) ( 25,15 ) ( 20,5 ) ( 10,5 ) ( 34,13)@xmath275 ( 18,32)@xmath95 ( 43,16)@xmath293 ( 57,16)@xmath294 ( 66,16)@xmath295 ( 57,13)@xmath296 ( 51,13)@xmath297 ( 48,27)@xmath298 ( 62,27)@xmath299 ( 48,2)@xmath300 ( 62,2)@xmath301 ( 40,15)(1,0)38 ( 55,0)(0,1)38 ( 45,15)(1,2)5 ( 50,25)(1,0)10 ( 60,25)(1,-2)5 ( 65,15)(-1,-2)5 ( 60,5)(-1,0)10 ( 50,5)(-1,2)5 ( 50,16)(0,-1)2 ( 60,16)(0,-1)2 ( 55,15 ) ( 55,15 ) ( 55,15 ) ( 45,15 ) ( 50,25 ) ( 60,25 ) ( 65,15 ) ( 60,5 ) ( 50,5 ) ( 74,13)@xmath275 ( 58,32)@xmath95