with the commissioning of the large hadron collider ( lhc ) , it will soon be possible to test models of multiparticle production in hadron - hadron collisions up to an energy of 14 tev . charged particle densities , @xmath6 , especially the values at mid - rapidity and their dependence on center - of - mass energy @xmath7 ( c.m.s . )
are important for understanding the mechanism of hadron production and the interplay of soft and hard scattering contributions in the lhc energy range .
the rate of parton - parton and multi parton - parton ( mpi ) scattering are strongly correlated to the observed particle multiplicity ( related also to _ initial entropy _ and _ initial energy density _ generated in the collision process ) .
new data on inclusive charged particle distributions from the lhc in @xmath8 collisions have become available @xcite .
these results complement previous data on @xmath8 and @xmath9 collisions taken at lower energies @xmath7 = 0.02 - 1.96 tev @xcite .
many of these measurements have been used to constrain phenomenological models of soft - hadronic interactions and to predict properties at higher energies @xcite .
the recent lhc data may lead to a better theoretical understanding based on a quantum chromodynamics ( qcd ) approach @xcite .
the heavy ion jet interacting generator ( hijing ) @xcite and hijing / bb v1.10 models @xcite have been used extensively to study particle production in @xmath8 collisions and to determine the physical properties of the ultra - dense matter produced in relativistic heavy - ion collisions . in the lund @xcite and dual - parton ( dpm ) @xcite models
multi qcd strings or flux tubes were proposed to describe soft multiparticle production in longitudinal color fields .
color exchange between high @xmath10 partons in the projectile and target create confined color flux tubes of tension ( @xmath11 gev / fm ) , that must neutralize through pair production or color singlet hadronization approximately uniformly in rapidity . in nucleus - nucleus collisions , the @xmath12 enhancement of the local parton density of high @xmath10 partons allows for higher color casimir representations to be excited . those flux tubes with stronger longitudinal color fields than in average @xmath8 reactions
have been called color ropes @xcite and naturally have higher string tensions @xcite . recently
, an extension of color glass condensate ( cgc ) theory has proposed a more detailed dynamical `` glasma '' model @xcite of color ropes .    in the hijing model
@xcite the soft beam jet fragmentation is modeled by simple diquark - quark strings as in the lund model with multiple gluon kinks induced by soft gluon radiations .
hard collisions are included with standard perturbative qcd ( pqcd ) as programmed in the pythia generator @xcite .
however , hijing differs from pythia by the incusion of a geometric scaling multiple jet production model .
thus this model contains both longitudinal field induced _ soft _ beam jet multiparticle production and collinear factorized pqcd based _ hard _ multiple jet production for @xmath13 gev / c .    a systematic comparison with data on @xmath8 and @xmath9 collisions in a wide energy range @xcite revealed that minijet production and fragmentation as implemented in the hijing model provide a simultaneous and consistent explanation of several effects : the inclusive spectra at moderate transverse momentum ( @xmath14 ) ; the energy dependence of the central rapidity density ; the two particle correlation function ; and the degree of violation of koba - nielsen - olesen ( kno ) scaling @xcite , @xcite up to tevatron energy ( @xmath15 tev ) .
however , the model failed to describe the dependence of the mean value of transverse momentum ( @xmath16 ) on charged particle multiplicity ( @xmath17 ) .
wang argued @xcite that by requiring high @xmath17 within a limited pseudorapidity ( @xmath18 ) range one necessarily biases the data towards higher @xmath14 minijets , hence the observed increase of @xmath16 versus ( vs. ) @xmath6 @xcite .
this effect has also been associated with the presence of transverse flow of the hadronic matter @xcite and was proposed as possibly due to _ quark - gluon plasma _
( qgp ) formation already in @xmath8 collisions .
initial states of color gauge fields produced in high - energy heavy - ion collisions have also recently been discussed in ref .
decay of a strong color electric field ( scf ) ( @xmath19=@xmath20 v / m ) due to the schwinger mechanisms @xcite plays an important role at the initial stage of heavy - ion collisions at ultra - relativistic energies .
a thermalization scenario based on the analogy between schwinger mechanisms and the hawking - unruh effect has been proposed @xcite .
it was also suggested that the back - reaction and screening effects of quark and anti - quark pairs on external electric field could even lead to the phenomenon of plasma oscillations @xcite .
recently , the schwinger mechanism has been revisited @xcite and pair production in time - dependent electric fields has been studied @xcite .
it was concluded that particles with large momentum were likely to have been created earlier , and for very short temporal widths ( @xmath21 , where the compton time @xmath22 ) and as a consequence the schwinger formula could underestimate the reachable particle number density . in previous papers
, we have shown that the dynamics of strangeness production in @xmath8 and au + au collisions at relativistic heavy ion collider ( rhic ) energies deviates considerably from calculations based on schwinger - like estimates for homogeneous and constant color fields and point to the contribution of fluctuations of transient strong color fields ( scf ) @xcite .    in this paper
we extend our study of the dynamic consequences of scf in the framework of the hijing / bb v2.0 model @xcite to particle production in hadron - hadron collisions at lhc energies .
we explore dynamical effects associated with long range coherent fields ( _ i.e _ , strong longitudinal color fields , scf ) , including baryon junctions and loops @xcite , with emphasis on the new observables measured in @xmath8 collisions at the lhc by the alice , atlas and cms collaborations .
our study aims to investigate a broad set of observables sensitive to the dynamics of the collisions , covering both longitudinal and transverse degree of freedom .
in addition , this study is intended to provide the @xmath8 baseline to future extrapolations to lhc studies of proton - nucleus ( @xmath23 ) and nucleus - nucleus ( @xmath24 ) collisions .
in hijing / bb v2.0 , in addition to conventional quark - diquark longitudinal electric fields , novel color flux topology junction anti - junction ( jj ) loops are also implemented . in a dual superconductor model of color confinement for the three - quark positioning in a y - geometry ,
the flux tubes converge first toward the center of the triangle and there is also another component that runs in the opposite direction ( like a _ mercedes star _ ) @xcite .
unlike the conventional diquark - quark implemented in lund and the hijing model @xcite , the hijing / bb v1.10 @xcite model allows the diquark - quark to split with the three independent flux lines tied together with an @xmath25 a junction and terminating on three delocalized fundamental casimir quarks . we introduced @xcite a new version ( v2.0 ) of hijing / bb that differs from hijing / bb v1.10 @xcite in its implementation of additional more complex flux topologies via junction anti - junction ( jj ) loops .
we parametrize the probability that a junction loop occurs in the string .
moreover , we enhance the intrinsic ( anti)diquark - quark @xmath14 kick ( by a factor @xmath26 = 3 ) of all ( @xmath27-qq ) strings that contain one or multiple jj loops .
the reason for this is the mechanism behind the dynamic of diquark - quark breaking ( see ref .
@xcite for details ) .    in string fragmentation phenomenology
, it has been proposed that the observed strong enhancement of strange particle production in nuclear collisions could be naturally explained via strong color field ( scf ) effects @xcite . for a uniform chromoelectric flux tube with field ( _ e _ ) , the pair production rate @xcite per unit volume for
a ( light)heavy quark ( @xmath28 ) is given by @xmath29 where @xmath30 ( diquark ) , @xmath31 ( strange ) , @xmath32 ( charm ) or @xmath33 ( bottom ) .
the _ current quark masses _ are @xmath34 = 0.45 gev @xcite , @xmath35 gev , @xmath36 gev , and @xmath37 gev @xcite .
the _ constituent quark masses _ of light non - strange quarks are @xmath38 = 0.23 gev , of the strange quark is @xmath39=0.35 gev @xcite , and of the diquark is @xmath40 gev @xcite .
note that @xmath41 is the effective string tension in low energy @xmath8 reactions with @xmath42 gev / fm and @xmath43 , @xmath44 are the second order casimir operators ( see ref .
an enhanced rate for spontaneous pair production is naturally associated with `` _ strong chromo - electric fields _ '' , such that @xmath45 1 _ at least some of the time_. in a strong longitudinal color electric field , the heavier flavor suppression factor @xmath46 varies with string tension via the well known schwinger formula @xcite , since @xmath47 for @xmath48 , @xmath31 , @xmath32 or @xmath33 and @xmath49 , @xmath50 . in our calculations
, we assume that @xmath51 = 0.5 gev , @xmath52 = 0.28 gev , @xmath53 = 1.30 gev .
therefore , the above formula implies a suppression of heavier quark production according to @xmath54 : @xmath50 : @xmath55 : @xmath31 : @xmath32 @xmath56 1 : 1 : 0.02 : 0.3 : 10@xmath57 for the vacuum string tension @xmath1 = 1 gev / fm . for a color rope , on the other hand , if the _
average string tension _
value @xmath58 increases , the suppression factors @xmath46 increase .
we show below that this dynamical mechanism improves considerably the description of the strange meson / hyperon data at the tevatron and at lhc energies .
saturation physics is based on the observation that small - x hadronic and nuclear wave functions , and , thus the scattering cross sections as well , are described by the same internal momentum scale known as the _ saturation scale _ ( @xmath59 ) .
a recent analysis of @xmath8 data up to lhc 7 tev has shown that , with the @xmath60 factorized ( glr ) gluon fusion approximation @xcite , the growth of the @xmath5 can be accounted for if the saturation scale grows with c.m.s . energy as @xmath61 with @xmath62 .
the saturation scale is also increasing with atomic number as @xmath63 @xcite .
it was argued that the effective string tension ( @xmath58 ) of color ropes should scale with @xmath64 @xcite .
however , in hijing the string / rope fragmentation is the only soft source of multiparticle production and multiple minijets provide a semi - hard additional source that is computable within collinear factorized standard pqcd with initial and final radiation ( dglap evolution @xcite ) . in order to achieve a quantitative description , within our hijing / bb framework
we will show that combined effects of hard and soft sources of multiparticle production can reproduce the available data in the range @xmath65 tev only with a reduced dependence of the effective string tension on @xmath7 .
we find that the data can be well reproduced taking @xmath66 where @xmath67 = 1 gev / fm is the vacuum string tension value , @xmath68 = 1 gev@xmath69 is a scale factor , and @xmath70 is adjusted to give @xmath71 gev / fm at the rhic energy @xmath72 tev .
our phenomenological @xmath73 is compared to @xmath74 in fig .
[ fig : kappa_vs_cgc ] , where @xmath58 = 1.40
gev / fm at @xmath7 = 0.017 tev increases to @xmath58 = 3.14 gev / fm at @xmath7 = 14 tev .    [ h ! ]   in @xmath8 collisions ( eq .
[ eq : kappa_sup ] ) .
the values @xmath75 ( eq .
[ eq : larry10 ] ) are from cgc model fits to lhc data from ref .
[ fig : kappa_vs_cgc ] , title="fig:",height=226 ] 0.5 cm    the energy dependence of the string tension leads to a variation of the diquark / quark suppression factors , as well as the enhanced intrinsic transverse momentum @xmath60 .
these include i ) the ratio of production rates of diquark - quark to quark pairs ( diquark - quark suppression factor ) , @xmath76 ; ii ) the ratio of production rates of strange to non - strange quark pairs ( strangeness suppression factor ) , @xmath77 ; iii ) the extra suppression associated with a diquark containing a strange quark compared to the normal suppression of strange quark ( @xmath78 ) , @xmath79 ; iv ) the suppression of spin 1 diquarks relative to spin 0 ones ( apart from the factor of 3 enhancement of the former based on counting the number of spin states ) , @xmath80 ; and v ) the ( anti)quark ( @xmath81 ) and ( anti)diquark ( @xmath82 ) gaussian width .
as an example we plot in fig .
[ fig : mg - fig1s_u ] the energy dependence of the suppression factor @xmath83 , when the string tension values @xmath73 are taken from eq .
[ eq : kappa_sup ] in comparison with the values predicted using @xmath75 ( eq .  [ eq : larry10 ] from the cgc model fit @xcite ) .
[ h ! ] , using @xmath73 from eq .
[ eq : kappa_sup ] and using @xmath74 from ref .
@xcite are compared .
( note that this schwinger suppression factor is not used in the cgc model ) .
[ fig : mg - fig1s_u ] , title="fig:",height=226 ] 0.5 cm    the contributions of multiple jets to the multiplicity distributions in @xmath8 and @xmath9 collisions have been studied in detail in ref .  @xcite . within the hijing model
, one assumes that nucleon - nucleon collisions at high energy can be divided into soft and hard processes with at least one pair of jet production with @xmath84 .
a cut - off scale @xmath85 in the transverse momentum of the final jet production has to be introduced below which the interaction is considered non - perturbative and can be characterized by a finite soft parton cross section @xmath86 .
the inclusive jet cross section @xmath87 at leading order ( lo ) @xcite is @xmath88 where , @xmath89 depends on the parton - parton cross section @xmath90 and parton distribution functions ( pdf ) @xmath91 .
the summation runs over all parton species ; @xmath92 and @xmath93 are the rapidities of the scattered partons ; @xmath94 and @xmath95 are the light - cone momentum fractions carried by the initial partons .
the factor @xmath96 accounts for the next - to - leading order ( nlo ) corrections to the leading order jet cross section . in the default hijing model ( v.1.383 ) , the duke - owens parameterization @xcite of pdfs in nucleons is used . with the duke - owens parameterization of pdfs ,
an energy independent cut - off scale @xmath972 gev/@xmath32 and a constant soft parton cross section @xmath98 mb are sufficient to reproduce the experimental data on total and inelastic cross sections and the hadron central rapidity density in @xmath99 collisions @xcite .
our results have been obtained using the same set of parameters for hard scatterings as in the latest version of hijing ( v1.383 ) .
charged hadron multiplicity measurements are the first results of the lhc physics program .
data for @xmath8 collisions were reported by the alice , atlas , and cms collaborations @xcite .
the new data at mid - rapidity for non single diffractive interactions ( nsd ) and inelastic scattering ( inel ) are shown in fig .
[ fig : pp_dn_ch_0 ] , which includes also similar results at lower energies .
the main result is an observed sizeable increase of the central pseudorapidity density with c.m.s .
energy .
the main contribution to the multiplicity comes from soft interactions with only a small component originating from hard scattering of the partonic constituents of the proton .
in contrast to the higher @xmath14 regime , well described by pqcd , particle production in soft collisions is generally modeled phenomenologically to describe the different @xmath8 scattering processes : elastic scattering ( el ) , single diffractive ( sd ) , double diffractive dissociation ( dd ) , and inelastic non - diffractive scattering ( nd ) . experimentally , minimum bias events are a close approximation of nsd interactions , _
i.e. _ , @xmath100 , where @xmath101 is the total cross section .
the selection of nsd events is energy dependent and differs somewhat for different experimental triggers .
the event selection of inelastic processes ( inel ) includes sd interactions : @xmath102 .
the data therefore must be corrected for the sd component , involving model dependent calculations .
the results of the model calculations for both nsd ( left panel ) and inel ( right panel ) are also shown in fig .
[ fig : pp_dn_ch_0 ] .
solid lines depict results including the strong color field ( scf ) effects whereas the results without scf effects are shown as dashed lines . without scf effects
the model strongly overestimates the central charged particle density .
the absolute value and energy increase of the central rapidity density are well reproduced assuming the energy dependent string tension given in eq .
[ eq : kappa_sup ] . at higher lhc energies ( 2.36 and 7 tev ) , a discrepancy of 10 - 15 % is observed .
[ h ! ]     and
@xmath9 interactions for non - single - diffractive ( nsd ) ( left panel ) and inelastic ( inel ) ( right panel ) interactions as a function of c.m.s . energy .
the solid and dashed lines are the results with and without scf , respectively .
the data are from refs .
@xcite ( left panel ) and from refs .
@xcite ( right panel ) . only statistical error bars are shown .
the open stars at 7 and 14 tev ( alice fit ) , are obtained by a power law fit to lower energy data from refs .
@xcite , @xcite .
[ fig : pp_dn_ch_0 ] , title="fig:",height=151 ] 0.5 cm    as the colliding energy increases , the rate of multiple parton interactions ( mpi ) also increases , producing a rise in the central multiplicity .
the increase with energy in our phenomenology is due to the interplay of the increased mini - jet production in high colliding energy with scf effects . for an increase of strangeness suppression factors due to an increase of string tension with energy ( @xmath103 ) ,
the model predicts a decrease of produced pions due to energy conservation .
lower values of @xmath73 imply smaller values for strangeness suppression factors , therefore a higher multiplicity of mesons ( mostly pions ) .    [ t ! ]     + @xmath104 ( @xmath105 ) collisions at various c.m.s .
the solid histograms are the results with scf and jj loops .
the data are from refs .
@xcite ( ua5 ) , @xcite ( star ) , @xcite ( cdf ) , @xcite ( alice ) , @xcite ( cms ) .
only statistical error bars are shown .
[ fig : pp_dn_ch_eta ] , title="fig:",width=325,height=325 ] 0.5 cm    changing the effective value @xmath73 = 2.89 gev / fm to @xmath73 = 2.0 gev / fm results in an increased multiplicity of 11% at 7 tev where the effect is greatest .
in addition , the multiplicity depends also on the value of the cut - off parameter @xmath85 .
low values of @xmath85 imply high rates of parton - parton scattering and hence high levels of particle multiplicity . evidently for increasing values of @xmath85
the opposite is expected . at 7 tev , where the effect is also the largest , changing @xmath85 by @xmath106 gev results in a change of only @xmath107 for the central pseudorapidity density .
data on the charged particle pseudorapidity distribution are also available over a limited @xmath18 range @xcite @xcite .
these are presented in fig .
[ fig : pp_dn_ch_eta ] where we also include results obtained at 0.2 tev by the ua5 @xcite and star @xcite collaborations , and with cdf results @xcite obtained at the tevatron for @xmath9 collisions at 1.8 tev .
consistent with the discussion above , a scenario with scf effects ( solid histograms ) reproduces the measured multiplicity distributions well .
at all energies considered , theoretical calculations predict a central dip at mid - rapidity that is consistent with the observations . at 0.9 and 2.36
tev , the shape of the distribution measured by the alice collaboration is very well reproduced , while the cms results show a much flatter distributions than calculated ( also the case at 7 tev ) .
data over a larger rapidity range are needed to determine the shape of the falling density in the fragmentation region . for completeness , predictions at 10 and 14 tev , the higher lhc energies , are also shown .
the measured transverse momentum distributions for nsd events over an energy range @xmath108 tev are shown in fig .
[ fig : pp_ch_pt_all ] .
these recent measurements are performed in the central rapidity region and cover a wide @xmath14 range ( @xmath109 gev / c ) , where both hard and soft processes are expected to contribute .
the data of atlas and cms are measured in larger pseudorapidity intervals ( @xmath110 ) .
in contrast , alice and cdf measurements are in a very central region ( @xmath111 and @xmath112 , respectively ) .
the calculation takes into account the difference in acceptance , but as can be concluded from fig .
[ fig : pp_dn_ch_eta ] , this difference in pseudorapidity range has a negligible effect on the measured cross section .
[ h ! ]     , title="fig:",width=325,height=325 ] 0.5 cm    the model calculations including scf effects describe the data well at @xmath113 tev , but lead to a harder spectrum that observed at higher energy .
the model gives a fair description of the spectral shape at low @xmath14 but overestimate the data at high @xmath14 .
the discrepancy is highest ( up to a factor of three ) at @xmath114 tev .
we do not understand the source of this discrepancy and await higher @xmath14 data to draw a firm conclusion .
however , in our phenomenology this could indicate that jet quenching , _
i.e. _ suppression of high @xmath14 particles like that observed at rhic energies in nucleus - nucleus collisions , could also appear in @xmath8 collisions in events with large multiplicity .
this overestimation of high @xmath14 yield leads also to a similar overestimation of the mean transverse momentum ( @xmath16 ) and of the correlation of mean @xmath16 as a function of n@xmath115 , which we do not discuss here .    within our model
we generate events with different numbers of mini - jets up to some maximum values .
high numbers of mini - jets lead to higher multiplicity events . for events with ten mini - jets ( the maximum assumed in our calculation ) ,
the central charged particle pseudorapidity density could increase up to @xmath56 20 and the total multiplicity could be greater than 150 at tevatron energy ( 1.8 tev ) .
the measurements of two particle correlations over the entire azimuth could reveal the jet structure related to high @xmath14 particles @xcite .
the study of these correlations in events with high multiplicity could help us to draw a firm conclusion with regard to a possible jet quenching phenomenon in @xmath8 collisions .
the @xmath8 single particle inclusive @xmath14 spectra measurements are important for understanding collision dynamics , since the various particles show different systematic behavior , as observed at rhic energy @xcite . detailed theoretical predictions for single inclusive hadron production ( including strange hyperons ) are discussed in this section .
baryon - to - meson ratios ( b / m ) are experimental observables that can be used at the lhc for investigating multi - parton interactions and helping to understand the underlying physics @xcite .
[ h ! ]   = 1.8 tev of baryon and mesons transverse momentum at mid - rapidity ( @xmath116 ) ( a ) and their ratios for single strange particles ( @xmath117)/@xmath118 ( b ) in minimum bias @xmath9 collisions .
the solid and dashed lines have the same meaning as in fig .
[ fig : pp_dn_ch_0 ] . experimental results at @xmath7 = 1.8 tev are from ref .
@xcite ( cdf collaboration ) .
error bars include only the statistical uncertainties .
the ratios ( b ) have been calculated by us , dividing the spectra in the panel ( a ) .
[ fig : pp_ks0_l0_18 ] , title="fig:",height=151 ] 0.5 cm    unexpectedly high b / m ratios observed in @xmath119 collisions have been discussed in terms of recombination and coalescence mechanisms @xcite .
such high ratios at intermediate @xmath14 were also reported in @xmath8 collisions at rhic @xcite and at the tevatron @xcite . in @xmath8 collisions , however , a coalescence / hadronization scenario is not favored due to low phase space density in the final state .
our hijing / bb model , with scf effects included , provides an alternative dynamical explanation of the heavy - ion data at rhic energies .
we have shown that the model also predicts an increasing yield of ( multi)strange particles , thereby better describing the experimental data @xcite .
[ h ! ]   ( a ) and their ratio ( b ) for minimum bias @xmath9 collisions at @xmath7 = 7.0 tev .
the solid and dashed lines have the same meaning as in fig .
[ fig : pp_dn_ch_0 ] .
the experimental results ( a ) are from ref .
@xcite(cms collaboration ) .
error bars include only the statistical uncertainties .
the ratios ( b ) have been calculated by us , dividing the spectra in the panel ( a ) .
[ fig : pp_ks0_l0_7 ] , title="fig:",height=151 ] 0.5 cm    figure  [ fig : pp_ks0_l0_18 ] ( b ) shows a comparison of model predictions with cdf experimental data @xcite of the strange baryon - to - meson ratio @xmath120 at 1.8 tev .
the particle @xmath14 spectra are shown in the panel ( a ) .
the measured ratio is fairly well described within our phenomenology .
the larger string tension parameterization results in a predicted increase of the ratio @xmath120 by a factor of @xmath56 10 at the tevatron energy .
the @xmath14 spectra show that the increased ratio is due almost entirely to an increase of the @xmath121 cross section that is well described with @xmath122 .
the model underestimates the kaon production by 15 - 25 % .
the model predictions at 7 tev , currently the maximum energy where there are data @xcite , are shown in fig .
[ fig : pp_ks0_l0_7 ] ( b ) .
the model gives a good description of hyperon production ( @xmath123 ) , for which an increase by a factor of ten is still predicted if scf effects are considered .
however , our calculations underestimate by approximately a factor of two the yield of @xmath118 and , as a consequence , result in a higher ratio than that observed .
this result needs further investigation ( theoretical and experimental ) on heavy flavour production at this energy .
we note that the models pythia @xcite and energy - conserving partons off - shell remnants and spliting of partons ladders ( epos ) @xcite can not reproduce the observed high b / m ratio ( see fig . 6 and fig . 7 from ref .
the preliminary data reported by the cms collaboration indicate high hyperon yield .
comparisons with pythia results show that this model significantly underestimates the hyperon yields in @xmath8 collisions at 0.9 and 7 tev @xcite .
the strange particle ratios could also be the manifestation of new collective phenomena . in the epos model such an increase
is obtained if the production of a _ mini - plasma _ is considered in @xmath8 collisions @xcite , @xcite .
if confirmed by future measurements , the study of these observables could open a perspective on new physics in @xmath8 interactions .
similar conclusions are obtained from the study of the proton / pion ( @xmath124 , @xmath125 ) ratios where data exist at lower energies @xcite .
these data , shown in figure  [ fig : pp_pi_pbar_18 ] , are limited to low @xmath126 gev / c .
adding scf effects results in a very sizable increase of the ratio and our calculations provide a good description of the data in the measured range .
however , as the calculations indicate , to draw a final conclusion measurements at intermediate and high @xmath14 are needed .
the solid and dashed lines have the same meaning as in fig .
[ fig : pp_dn_ch_0 ] . experimental results for @xmath127 at @xmath7 = 0.54 tev ( left panel ) and at @xmath7 = 1.8 tev ( right panel ) are from ref .
@xcite ( e735 collaboration ) .
the results at mid - rapidity at @xmath7 = 0.54 tev ( left panel ) are from ref .
@xcite ( ua2 collaboration ) .
error bars include only the statistical uncertainties .
[ fig : pp_pi_pbar_18 ] , title="fig:",height=151 ] 0.5 cm    to the extent that the lhc experiments are able to identify hadron species , such data will provide vital input to validate this interpretation .
the model predictions at lhc energies for the @xmath14 dependence of the @xmath128 ratio are shown in fig .
[ fig : pp_pi_pbar_14 ] . an enhancement up to the highest lhc energy and a weak energy dependence , with a saturation that sets in for a c.m.s .
energy @xmath129 tev , is predicted .
note that preliminary data at 0.9 tev reported by the alice collaboration for @xmath14 spectra of pions ( @xmath130 ) and protons ( @xmath104 )
@xcite cover to @xmath131 gev / c .
the model results , with scf effects included ( dot - dashed histogram in fig .  [ fig : pp_pi_pbar_14 ] ) are consistent with a @xmath132 ratio derived from the spectra reported by alice at 0.9 tev in ref .
@xcite .
[ h ! ] ) for minimum bias events at mid - rapidity at lhc energies .
the upper curves correspond to calculations which include the effects of scf and jj loops .
the lower curves corresponds to calculations without scf effects .
the results are at 0.9 tev ( dot dashed - histograms ) , at 7 tev ( dashed histograms ) , at 10 tev ( dotted histograms ) and at 14 tev ( solid histograms ) .
[ fig : pp_pi_pbar_14 ] , title="fig:",height=226 ] 0.5 cm    in our approach , the dynamical mechanism that leads to such high values of b / m ratios is scf appearing at the initial stage of the interaction .
the scf mechanism strongly modifies the fragmentation processes ( strangeness suppression factors ) and thus results in a huge increase of ( strange)baryons .
this interpretation is also supported by more sophisticated theoretical calculations , in a scenario in which a time - dependent pulse for the initial strength of the color field is considered @xcite .
the large enhancement of the baryon - to - meson ratios demonstrates that scf could play an important role in multiparticle production in @xmath8 collisions at lhc energies and that high energy density fluctuations can reach very high densities , potentially comparable to those reached in central au + au collisions at rhic energies @xcite .      from the study of the baryon - anti - baryon asymmetry one can learn about the mechanism of baryon number transport .
baryon production via the conventional default quark - diquark mechanisms in the lund string model are known to be inadequate even in @xmath133 phenomenology .
this is one of the main reasons for our continued exploration of alternative baryon junction mechanisms .
the details of the new implementation of jj loops in hijing / bb v2.0 are described in ref .
@xcite sec ii a. in hijing / bb v2.0 the main two mechanisms for baryon production are quark - diquark ( @xmath134 ) string fragmentation and junction anti - junction ( jj ) loops @xcite in which baryons are produced approximately in pairs . in a junction loop
a color flux line splits at some intermediate point into two flux lines at one junction and then the flux lines fuse back at an anti - junction somewhere further along the original flux line .
the distance in rapidity between these points is chosen via a regge distribution  @xcite .
we assume that , out of the non single diffractive nucleon - nucleon ( nn ) interaction cross section ( @xmath135 ) , a fraction @xmath136 of the events excite a junction loop .
the probability that the incident baryon has a jj loop in p(a)+a collisions after @xmath137 simulated binary collisions is given by    @xmath138    where @xmath139 = 17 mb , @xmath140 is parametrized in the model , and the total inelastic nucleon - nucleon cross sections are calculated @xcite .
these cross sections imply that a junction loop occurs in @xmath8 collisions with a rather high probability ( at rhic energy @xmath141 for @xmath142 = 42 mb ) .
taking a constant value for @xmath139 results in a decrease with energy of the probability @xmath143 , due to a faster increase with energy of @xmath142 relative to @xmath140 .
the actual probability is modified also by string fragmentation processes for which we consider a threshold cutoff mass @xmath144 gev / c@xmath69 in order to have enough kinematic phase space to produce bb pairs .
we have investigated the sensitivity of the results to the value of parameters @xmath139 and @xmath145 and found no significant variation on pseudorapidity distributions of charged particles for 15 mb @xmath146 25 mb and for 4 gev / c@xmath69 @xmath147 6 gev / c@xmath69 .
baryon number transport is quantified in terms of the rapidity loss ( @xmath148 , where @xmath149 and @xmath150 are the rapidity of incoming beam and outgoing baryon , respectively ) and has been discussed within our model phenomenology for @xmath119 collisions in refs .
it was shown that hijing / bb v1.0 overestimate the stopping power and give a mild energy dependence of net - baryons at mid - rapidity .
the energy dependence of net - baryons at mid - rapidity per participant pair within hijing / bb v2.0 is proportional to @xmath151 , similar to the dependence predicted in ref .
@xcite , with the assumption that jj is the dominant mechanism .
this dependence is obtained assuming the following parameters @xcite : @xmath152 gev@xmath69 ( the usual parameter of regge theory ) , @xmath153 ( the reggeon intercept of the trajectory ) , and @xmath154 ( where @xmath155 ) for the pomeron intercept .    [ h ! ]     distributions ( a ) at mid - rapidity for the @xmath156 ratio at @xmath7 = 0.9 tev ; rapidity distributions are shown in part ( b ) .
the results are from three possible scenarios : without contributions from scf and jj loops ( dotted histograms ) , including only the effect of jj loops ( dashed histograms ) and including both effects ( scf and jj loops ) ( solid histograms ) .
the data are from ref .
@xcite ( alice collaboration ) . only statistical errors are shown .
[ fig : pp_rap_pbar_pt ] , title="fig:",height=151 ] 0.5 cm    recently , the alice collaboration reported results @xcite on mid - rapidity anti - proton - to - proton ratio in @xmath8 collisions at @xmath157 and 7 tev and equivalently the proton - anti - proton asymmetry , @xmath158 .
these data could be used to constrain regge inspired model descriptions of baryon asymmetry .
the authors state that , within statistical errors , the observed @xmath156 ratio shows no dependence on transverse momentum or rapidity in the limited measured acceptance ( @xmath159 ; @xmath160 gev / c ) . in fig .
[ fig : pp_rap_pbar_pt ] ( a ) are compared the hijing / bb v2.0 model predictions with the published data at 0.9 tev .
our model predicts negligible dependence on @xmath14 for @xmath161 gev / c , and a slight @xmath14 dependence at higher @xmath14 , where both effects ( jj loops and scf ) could contribute .
the full calculation is shown by the solid histogram .
the dotted line is the prediction without jj loops and scf effects while the dashed line includes the effect of jj loops only . over
the measured ranges the rapidity distribution of the ratio @xmath156 distribution is not sensitive to the various scenarios presented .
the scenario with combined effects results in a wider rapidity distribution at @xmath162 ( b ) .
the narrow structure observed near @xmath163 has no physical significance and we believe that it is likely due to a numerical artifact of our current implementation of fragmentation scheme .
separate proton and anti - proton rapidity distributions are , however , much more sensitive to scf effects , as seen in fig .
[ fig : pp_p_pbar_rap ] .
the model predicts a substantial increase ( by a factor of @xmath164 ) for p(@xmath105 ) , when scf are taken into account . due to the high cutoff mass @xmath165 gev
/ c@xmath69 the effect of jj loops is very small over the entire rapidity region .
analysis of the 7 tev data leads to similar conclusions , but shows less effect at high @xmath14 .
[ h ! ]     = 0.9 tev .
the histograms have the same meaning as in fig .
[ fig : pp_rap_pbar_pt ] .
[ fig : pp_p_pbar_rap ] , title="fig:",height=151 ] 0.5 cm    over the measured range the feed - down - corrected data for the @xmath156 ratio rises from @xmath166 at 0.9 tev to @xmath167 at 7 tev @xcite .
although the measured mid - rapidity ratio is close to unity there is a small but significant excess of protons over anti - protons corresponding to an asymmetry of @xmath168 and @xmath169 at @xmath170 and 7 tev , respectively . within our model using a regge intercept of @xmath171 = 1/2 results in a proton - anti - proton asymmetry @xmath172 = 0.032 ( 0.002 ) at 0.9 tev ( 7 tev ) .
the result overestimates the value of asymmetry at 0.9 tev .
these values show that a fraction of the baryon number associated with the beam particles is transported over rapidity intervals of less than seven units .
however , a better understanding of the baryon transport would require measuring the rapidity dependence of the asymmetry over a larger range .    in ref .
@xcite we discuss the predicted rapidity correlation length for production of baryon - anti - baryon pairs within our model . the predicted rapidity correlation length ( @xmath173 depends upon the value of the regge intercept @xmath174 @xcite .
a value of @xmath175 leads to rapidity correlations in the range @xmath176 , while a value @xmath177 @xcite is associated with infinite range rapidity correlations .
this kind of analysis and future measurements at lhc energies will help us to determine better the specific parameters characterizing jj loops used in the calculations of baryon number transport and/or baryon - anti - baryon asymmetry while helping to establish the validity of regge inspired models .
we have studied the influence of strong longitudinal color fields and of possible multi - gluon dynamics ( _ gluon junctions _ ) in particle production in @xmath8 collisions , with a focus on rhic , to tevatron , and lhc energies .
we have investigated a set of observables sensitive to the dynamics of the collisions , covering both longitudinal and transverse degree of freedom . a detailed comparison with newly available experimental data from the lhc has been performed .
we found that the inclusion of the multiple minijet source limits the growth of the string tension @xmath73 to be approximately only linear as a function of saturation scale @xmath59 ( eq .  [ eq : kappa_sup ] ) , in contrast to recent approaches @xcite where @xmath73 scales as @xmath75 ( eq .  [ eq : larry10 ] ) .
the reason for this is that in the cgc model the collinear factorized minijet mechanism is suppressed by geometric scaling to much higher @xmath14 .
future measurements at lhc energies ( @xmath7 = 7 and 14 tev ) , extended to high @xmath14 , will help to clarify the validity of this mechanism .
we have shown that scf could play an important role in particle production at mid - rapidity in @xmath8 collisions .
our calculations show that high - energy density fluctuations in @xmath8 collisions at lhc can reach densities comparable to those reached in central nuclear ( @xmath119 ) collisions at rhic .
a large enhancement of the ( strange)baryon - to - meson ratios that persists up to the highest lhc energy can be explained as an effect of scf that appears at the initial stage of the interaction .
the mechanisms of hadron production are very sensitive to the early phase of the collisions , when fluctuations of the color field strength are highest .
strong color field effects are modeled by varying the effective string tension that controls the @xmath178 and @xmath179 pair creation rates and strangeness suppression factors .
scf , therefore , may modify the fragmentation processes with a resultant huge increase of ( strange)baryons .
we show that both jj loops and scf effects could play an important role in baryon ( anti - baryon ) production at mid - rapidity in @xmath8 collisions at lhc energies . introducing a new jj loop algorithm in the framework of hijing
/ bb v2.0 leads to a consistent and significant improvement in the description of the recent experimental results for proton - anti - proton and for baryon - anti - baryon asymmetry in comparison to the older versions hijing / b or hijing / bb v1.0 @xcite .
we have shown that baryon number transport is suppressed for @xmath180 , a result that is confirmed by recent alice measurements @xcite .
the present study is limited to the effect of initial state baryon production via possible junction dynamics in strong fields .
it would be very useful to consider a generalization of back reaction effects @xcite to the case not only of pair production relevant for mesons but to the more difficult three string junction configurations needed to describe baryon production .    a greater sensitivity to scf effects
is expected also for open charm and bottom production @xcite .
in particular , measurements of rapidity and @xmath14 distributions for particles involving charm and bottom quark , would provide an important test of the relevance of scf fluctuations , helping us to determine values of the suppression factors @xmath46 ( where q = qq , @xmath31 , @xmath32 , @xmath33 ) , which have strong dependence on the main parameters of qcd ( the constituent and current quark masses ) and on the system size . even though the success of this procedure has been clearly illustrated here , a fuller understanding of particle production and especially of ( multi)strange particles in ultra - relativistic @xmath8 collisions at the lhc
remains an exciting open question , and will continue to challenge many theoretical ideas .
0.2 cm    * acknowledgments : * we thank s. das gupta and s. jeon for useful discussions and continued support .
we thank p. levai for helpful discussions and suggestions throughout this project .
vtp acknowledges computer facilities at columbia university , new york , where part of these calculations were performed .
this work was supported by the natural sciences and engineering research council of canada .
this work was also supported by the division of nuclear science , of the u. s. department of energy under contract no .
de - ac03 - 76sf00098 and de - fg02 - 93er-40764 .
t.  pierog , s.  porteboeuf , i.  karpenko and k.  werner , to appear in the proceedings of 45th rencontre de moriond , qcd and high energy interactions , 13 - 20 march 2010 , la thuile , aosta valley , italy , arxiv:1005.4526 [ hep - ph ] .
e.  eichten , i.  hinchliffe , k.  d.  lane and c.  quigg , rev .  mod .
* 56 * , 579 ( 1984 ) [ addendum - ibid .
* 58 * , 1065 ( 1986 ) ] .
d.  w.  duke and j.  f.  owens , phys .  rev .
d * 30 * , 49 ( 1984 ) .
d.  acosta  _ et al .
_ ,  cdf  collaboration ,  phys .
d * 65 * , 072005 ( 2002 ) .
h.  ricaud , a.  kalweit and a.  maire , j.  phys .
g * 37 * , 094049 ( 2010 ) .
r.  j.  fries , b.  muller , c.  nonaka and s.  a.  bass , phys .
c * 68 * , 044902 ( 2003 ) .
v.  greco , c.  m.  ko and p.  levai , phys .
c * 68 * , 034904 ( 2003 ) .
r.  c.  hwa and c.  b.  yang , phys .
c * 67 * , 034902 ( 2003 ) .
b.  i.  abelev _ et al . _ [ star collaboration ] , phys .
c * 75 * , 064901 ( 2007 ) .
d.  e.  acosta _ et al .
_ [ cdf collaboration ] , phys .
d * 72 * , 052001 ( 2005 ) .
cms collaboration , cern cms physics analysis summary qcd-10 - 007 ( 2010 ) , http://cdsweb.cern.ch/record/1279344 ; keith a. ulmer for the cms collaboration , arxiv:1012.1605 [ hep - ex ] .
k.  werner , t.  hirano , i.  karpenko , t.  pierog , s.  porteboeuf , m.  bleicher and s.  haussler , nucl .
suppl .   * 196 * , 36 ( 2009 ) .
k.  werner , i.  karpenko and t.  pierog , j.  phys .
ser .   * 230 * , 012026 ( 2010 ) .