the intuitive feeling that hadron matter should be maximally perturbed in the high energy extremely inelastic collisions was a main reason of our effort to consider such processes .
it was a hope to observe new dynamical phenomena , or new degrees of freedom , unattainable in other ordinary hadron reactions .
this paper presents a first attempt to describe the particularity of considered processes , to give a review of existing models prediction and , at the end , we will offer the field - theoretical formalism for hadron inelastic processes .
thus , considering hadrons mean multiplicity @xmath3 as a natural scale of the produced hadrons multiplicity @xmath4 at given cm energy @xmath5 , we would assume that n>>. 0.1at the same time we wish to have n<<n_max=/m , 0.2where @xmath6 gev is the characteristic hadron mass .
the last restriction is introduced to weaken the unphysical constraints from the finite , for given @xmath7 , phase space volume .
we should assume therefore that @xmath7 is high enough .
the multiple production is the process of colliding particles where kinetic energy is dissipated into the mass of produced particles .
then one may validate that the entropy @xmath8 accedes its maximum in the domain ( ) since the multiplicity @xmath4 characterizes the rate of stochastization , i.e. the level of incident energy dissipation over existing ( free ) degrees of freedom .
there is also another quantitative definition of our reactions .
let @xmath9 be the energy of the fastest particle in the given frame and let @xmath10 be the total incident energy in the same frame .
then the difference ( @xmath11 ) is the energy spent on production of the less energetic particles .
it is useful to consider the inelasticity coefficient = = 1 - 1 .
0.2ait defines the portion of spent energy .
therefore , we wish to consider processes with 1 - 1 . 0.2cfollowing ( ) we will assume that 1->>. 0.2dtherefore , the kinetic energy of produced particles in our processes can not be arbitrarily small .
using thermodynamical terminology , we wish investigate the production and properties of comparatively ` cold ' multi - hadron ( mostly of @xmath12-mesons ) state .
we would like to note from the very beginning that we have only a qualitative scenario of such states which may be produced and the review of the corresponding why s may be considered as the main purpose of this paper .
although at the end of this paper ( appendix k ) we will describe principal features of possible field - theoretical solution of our problem .
it should be noted the absence of any authentic experimental information concerning discussed processes .
moreover , actually the hadron inelastic interactions with a set peculiar to them of unsolved theoretical problems will be considered .
nevertheless we suggest to work in this field in spite of these difficulties because the system with extremal properties may be more transparent since the asymptotics always simplify a picture .
we would demonstrate this idea and will try to put it in the basis of developed theoretical methods .
absence of experimental information about so high inelastic hadron processes is the consequence of the smallness of corresponding cross sections . besides this it was unclear for what purpose the experimental efforts should be done .
we would like to convince the reader that the discussed problem is interesting and important .
for instance , we will discuss a possibility that asymptotics over @xmath4 may replace in definite sense the asymptotics over @xmath5 .
a short address to experimentalists will be given in concluding sec.4.1.2 .
we hope that the paper would be useful both for theorists and experimentalists . by this reason
the main text of this paper will contain only the qualitative discussion of the problems and results . the quantitative proof , formulation of pure theoretical methods etc .
are added in the appendices
. considered extremal problem is a good theoretical laboratory and is described in the appendices and theoretical methods may be applied for other physical problems .
we would like to point out that the special technique was built for the problem discussed , see sec.2 .
_ having the very high multiplicity ( vhm ) state it is natural to use the thermodynamical methods .
_ we will offer for this purpose the real - time @xmath0-matrix interpretation of thermodynamics .
it can be shown in what quantitative conditions it will coincide with simpler canonical imaginary - time matsubara formalism
. we will give also the generalization of the real - time finite - temperature perturbation theory in the case of local temperature @xmath13 distribution , when @xmath14 .
this will allow to use the thermodynamical description if the system is far from equilibrium .
 _ the particle spectrum in the vhm region is soft .
it is just a situation when the collective phenomena should be important . _ to describe this phenomena , the decomposition on correlators will be adopted .
the origin of this decomposition lie in the mayer s ` group decomposition ' . in multiple production physics
this decomposition is known also as the ` multi - component description ' .
it based on the idea that the multiple production process may include various mechanisms .
in sec.3 we will investigate model predictions for vhm region .
we would like to note two main conclusions :     existing multiperipheral - type models are unable to describe the vhm region .
the infrared region of the pqcd becomes important even if the constraint ( ) is taken into account .
in sec.2.1 we will try to formulate the phenomenology of our problem , i.e. the way the vhm processes may be described and what type of phenomena one may expect .
the importance of thermodynamical methods will became evident and we will offer in sec.2.2 a general description of the corresponding formalism .    it is important to note that we may classify the possible asymptotics over @xmath4 .
we will find that there exist only three classes of asymptotics .
this will simplify consideration definitely restricting the possibilities .
in sec.2.3 we will use the thermodynamical language to give a physical interpretation of these classes of asymptotics .
we will see in result that in our choice of the vhm final state this should lead to reorganization of multiple production dynamics : we will get out of the habitual multiperipheral picture in the vhm domain .
moreover , one may assume that the semiclassical approximation becomes exact in the vhm domain .
this naturally leads to the idea to search such a scheme of calculation which depends on the choice of final state
. quantitative description of this idea may be realized as is described in appendix k.      the vhm production phenomena includes two sub - problems .
first of all it is the dynamical problem of incident energy degradation into the secondary particle energies and the second one is the description of the final state
.    we will start discussion in sec.2.1.1 from the second part of the problem to explain that the statistical methods are essential for us .    in sec.2.1.2
we will try to outline at least qualitatively the main mechanisms of hadron production .
the peculiarity of hadron production phenomena consists in the presence of hidden constraints , the consequence of local non - abelian gauge symmetry .
the constraints may prevent thermalization and the incident energy dissipation is confused in this case .
just the ` confusing ' effect is dominant in the hadron multiple production processes if @xmath15 .
the expected change of dominant mechanism of hadron production is discussed sec.2.1.3 .
it is important that , in spite of hidden constraints , the system may freely evolve to define the vhm state .
such a vhm state should be in equilibrium .
formal definition of the ` equilibrium ' notion will be given in sec.2.2.2 .
the problem described contains small parameters @xmath16 and @xmath17 . to have the possibility of estimation of contributions in accordance with these parameters
one should include them into the formalism .
this becomes possible if and only if the integral quantities are calculated .
so , the multiplicity @xmath4 is an index only if the multiple production amplitudes @xmath18 are considered .
but the cross section @xmath19 is a nontrivial @xmath20 of @xmath4 .
we will calculated by this logic mostly integrals of @xmath21 , living out the consideration the amplitudes , see also sec.2.2 , where the first realization of this idea is offered ( the @xmath22 attempt to realize this idea one may find in ) .
this is a general methodological feature of our consideration .
the multiple production cross section @xmath19 falls down rapidly in the discussed very high multiplicity ( vhm ) domain ( ) and for this reason the multiplicities @xmath23 are not accessible experimentally . at
the lhc energy @xmath24 is valid and we will assume that @xmath25 is just the discussed vhm region ( @xmath26 at the lhc energy ) .
we will explain later why n~^2 0.1is chosen for the definition of vhm region .    generally speaking , having the state of a large number of particles , it is reasonable to depart from an exact definition of the number @xmath4 of created particles , their individual energies @xmath27 , momenta @xmath28 , etc . since they can not be defined exactly by experiment .
indeed , for instance , full reconstruction of kinematics is a practically impossible task because of neutral particles , neutrinos , the more so as @xmath29 is considered .
we suppose that nothing will happen if @xmath4 is measured with @xmath30 accuracy since @xmath31 is easily attainable in the vhm region . besides it is practically impossible to deal with theory which operates by the @xmath32 ( @xmath33 ! ) variables .
artificial reduction of the set of the necessary variables may lead to a temporary success only .
indeed , the last thirty years of multiple production physics development was based on the inclusive approach , when the measured quantities ( cross sections ) depend on a few dynamical parameters only .
but later on the experiment and its fractal analyses shows that the situation is not so simple , also as , for instance , for the classical turbulence .
so , the event - by - event experimental data shows that the particle density fluctuation is unexpectedly large and the fractal dimension @xmath34 is not equal to zero .
we know that if the fractal dimension is non - trivial , then the system is extremely ` non - regular ' .
so , @xmath350.3 for the perimeter of great britain and @xmath350.5 for norway .
the discrepancy marks the fact that the shore of norway is much more broken then of great britain .
it is noticeable that the fractal dimension @xmath34 crucially depends on the type of reaction , incident energy and so on .
it is evident that one may choose from @xmath32 an arbitrary finite set of variables to characterize the multiple production process .
but the fractal analyses shows that such approach would lead to the same effect as if one may hear , for example , only the first violin of mahler s music .
so , it is important to understand when the restricted set of dynamical variables will admit to describe the process ( state ) @xmath36 .
the same problem was solved in statistical physics , where the ` rough ' description by a restricted number of ( thermodynamical ) parameters is a basis of its success , see the discussion of rough variables description e.g. in the review of uhlenbeck .
we will search the same solution desiring to build a complete theory of the vhm hadron reactions .
we want to note that just vhm process may be in this sense ` simple ' : at all evidence , the system becomes ` quiet ' in the vhm region and by this reason its ` rough ' thermodynamical description is available .
it seems natural , therefore , to start investigation of multiple production phenomena from the ( extremely rare ) vhm processes .
the dominant inelastic hadron processes at @xmath15 are saturated by production of low transverse momentum hadrons .
one of the approaches explains this phenomena by the nonperturbative effect of quarks created from vacuum .
corresponding dynamics looks as follows . at the expanse of transverse kinetic motion color charges may separated at large distances . nevertheless the transverse motion is suppressed since separation leads to increasing polarization of vacuum , because of confinement phenomenon .
then , as in qed , the vacuum becomes unstable in regard to the tunneling creation of real fermions .
just on their creation , the transverse kinetic energy is spent and , as a result , particles can not have , with exponential accuracy , high transverse energies .
this picture is attractive being simple and transparent , but despite numerous efforts , there is no quantitative description of this phenomenon till now .
briefly , the problem is connected with the unknown mechanism of strong coloured electric fields formation among distant colored charges , see also .
one may use other terms .
the soft channel of multiple production means the long - range correlation among hadron coloured constituents . under this special correlation
the non - abelian gauge field theory conservation laws constraint were implied .
they are important in dynamics since each conservation law decreases the number of independent degrees of freedom at least on one unity ( this may explain why hadrons @xmath37 ) , i.e. it has the nonperturbative effect .
moreover , in so - called integrable systems each independent integral of motion ( in involution ) reduces the number of degrees of freedom on two units . as a result
there is no thermalization in such systems and the corresponding mean multiplicity @xmath38 should be equal to zero , see appendix k.    existence of multiple production , @xmath39 , testifies the statement that the thermalization phenomena exist in hadron processes , i.e. the system of yang - mills fields is not completely integrable .
but the most probable process with @xmath15 did not lead to the final state with maximal entropy since @xmath40 , i.e. the definite restrictions on the dissipation dynamics should be taken into account .
such problems , being intermediate , are mostly complicated ones
.    the quantitative theory of this phenomena may lead to deep revision of the main notions of existing quantum field theory , see appendix k.
so , the dynamical display of hidden conservation laws of the hadron system are probably unstable since we expect that the system is not completely integrable , solitary field configurations @xmath41
. then the quantum theory should be able to describe quantum excitations of these fields , i.e. to count the fluctuations of ` curved ' manifolds .
the canonical perturbation theory methods , formulated in terms of creation and annihilation of particles in the external field @xmath41 , are too complicated , see and references cited therein .
for this reason , existing calculations usually do not exceed the semiclassical approximation .
we hope that , as described in appendix k , the quantization scheme would be able to solve this problem ( see also the example described there ) .
another approach assumes that the special ` @xmath42-channel ' ladder type feynman diagrams are able to describe the @xmath15 region .
this approach did not take into account the confinement nonperturbative effects introducing the hardly controllable supposition that the free quarks and gluons may form a complete set of states .
formally this is right , but at all evidence , the decomposition on this fock basis is realized in the non - abelian gauge field theories on zero measure ( see appendix k ) .
nevertheless one may reject this argument assuming that the process is happening on a sufficiently small distances .
corresponding contribution came from the so called ` hard pomeron ' .
but the intrinsic problems of the accuracy of chosen logarithmic approximations , the understanding of the so called non - logarithmic corrections , of the fate of the infrared divergences stay unsolved till now .
so , the multiple production of soft hadron phenomena seems unsolvable on the day - to - day level of understanding and we may with pure conscience move it away .
this is why the hadron inelastic reactions lost some popularity , migrating the last two decades to the class of ` non - interesting ' problems . yet , in a number of modern fundamental experiments , multiple production plays , at least , the role of background to the investigated phenomena and for this reason we should be ready for the quantitative estimation of it .
our hope to describe such a complicated problem as the multiple production phenomenon in the vhm conditions is based on the following idea . at the very
beginning of this century , a couple p. and t.euhrenfest , had offered a model to visualize boltzmann s interpretation of the irreversibility phenomena in statistics .
the model is extremely simple and fruitful .
it considers the two boxes with @xmath43 numerated balls . choosing the label of the balls @xmath44 one must take the ball with the corresponding label from one box and put it into another one .
one may repeat this action an arbitrary number of times @xmath42 .
starting from the highly ` nonequilibrium ' state with all balls in one box , @xmath45 , it is seen stationary with @xmath42 tendency to equalization of the number of balls in the boxes fig.1 .
the stationarity means that the number of balls in the other box rises @xmath46 at least on an early stage of the process .
this signifies presence of an @xmath47 flow ( of balls ) toward the preferable ( equilibrium ) state .
one can hope that this model reflects a physical reality of nonequilibrium processes with initial state very far from equilibrium .
a theory of such processes with ( irreversible ) flow toward a state with maximal entropy should be sufficiently simple being close to the stationary markovian
.    the vhm production process may be , at least on an early stage , stationary markovian .
if this is so then one may neglect long - range effects , non - perturbative as well , since they are not markovian as follows from the experience described in sec.2.1.2 .
this is possible if the vhm process is happening so fast ( being the short - range phenomenon ) that the confinement forces became ` frozen ' .
it can be shown that the quantitative reason for this phenomena is a fast ( exponential ) reduction with @xmath4 of the soft channel contribution into the hadron production process .
so , we expect a change of the multiple production dynamics in the vhm region .    thus , the main input idea consists from two general propositions .
the first of them is following :    ( * i * ) _ the hadron vhm production processes should be close to the stationary markovian .
_    ` freezing ' the confinement constraints , the entropy @xmath8 may exceed for given energy @xmath5 its available maximum in the vhm domain .
then one can assume that the vhm final state is in ` equilibrium ' , or is close to it .
so ,    ( * ii * ) _ the vhm final state should be close to equilibrium _
+ is our second basic proposition
.    we would select and appreciate particle physics models in accordance with these propositions .
the field - theoretical description of statistical systems at a finite temperature is based usually on the formal analogy between imaginary time and inverse temperature @xmath48 .
this analogy is formulated by schwinger as the ` euclidean postulate ' and it assumes that ( i ) the system is in equilibrium , i.e. it should allow the arbitrary rearrangement of states of temporal sequence in the described process , and ( ii ) there is not special space - time long - range correlations among states of the process , i.e , for instance , the symmetry constraints should not play a crucial role . we do not know _ ad hoc _ whether or not to apply the ` euclidean postulate ' for given @xmath4 and @xmath7 , even if ( ) is satisfied .
for this reason we are forced to formulate the theory in the natural real - time terms .    the first important quantitative attempt to build the real - time finite - temperature field theory discovered the formal problem of the so called ` pinch - singularities ' .
further investigation of the theory has allowed to demonstrate the cancellation mechanism of these unphysical singularities .
this is attained by doubling the degrees of freedom : the green functions of the theory represent @xmath49 @xmath50 @xmath49 matrix .
it surely makes the theory more complicated , but the operator formalism of the thermo - field dynamics shows the unavoidable character of this complication .
notice that the canonical real - time finite temperature field - theoretical description of the statistical systems based on the kubo - martin - schwinger ( kms ) boundary condition for a field is : ( t)=(t - i ) .
22.1 it , without fail , leads to the @xmath51 fluctuation - dissipation conditions ( see also ) . by this reason
it can not be applied in our case , where the dissipation problem is solved .
the origin of this boundary condition is shown in appendix a.    we will use a more natural microcanonical formalism for particle physics .
the thermodynamical ` rough ' variables are introduced in these approach as the lagrange multipliers of corresponding conservation laws .
the physical meaning of these ` rough ' variables are defined by the corresponding equations of state .
we shall use the @xmath0-matrix approach which is natural for the description of the time evolution .
( the @xmath0-matrix description was used also in . ) for this purpose the amplitudes = a_nm(p , q ) 22.2 of the @xmath52- into @xmath4-particles transition will be introduced .
the in- and out - states must be composed from mass - shell particles . using these amplitudes we will calculate & r_nm(p , q)=|a_nm(p , q)|^2= + & = < p_1,p_2,
... ,p_n|q_1,q_2, ... ,q_m >
< q_1,q_2, ...
,p_n>. 22.3this will lead to the doubling of the degrees of freedom .
the temperature description will be introduced ( see also ) noting that , for instance , & d_n=|a_nm(p , q)|^2d_n(p ) , + & d_n(p)=_1^n ,  ( p)= ( p^2 + m^2)^1/2 , 22.4 is the differential measure of the final state .
it is a first example where the usefulness of the probability - like quantity @xmath53 is seen .
the measure ( ) is defined on the energy - momentum shell : _
i=1^np_i = p .
22.4ait should be underlined that @xmath54 are the translationally invariant amplitudes and four equalities : _
i=1^np_i=_i=1^mq_i 22.4bare obeyed identically .
so , ( ) are the constraints and to take them into account one may multiply @xmath55 on @xmath56 where @xmath57 is the time - like 4-vector . it is evident that integration over @xmath57 with factor @xmath58 gives the constraints ( ) .
one may simplify the calculation assuming that all calculations are performed , for example , in the cm frame @xmath59 .
then one may ignore the space components considering @xmath60 .
this is the equivalent of the assumption that only the energy conservation law is important .
the last step is the substitution @xmath61 , where @xmath62 is our lagrange multiplier . to define its physical meaning
one should solve the equation of state : e =- d_n_k=1^ne^-(p_k)-d_n ( ) .
3such a definition of temperature as the lagrange multiplier of the energy conservation law is obvious for microcanonical description
.    the initial - state temperature will be introduced in the same way , taking into account ( ) .
so , we will construct the two - temperature theory .
it is impossible to use the kms boundary condition in such a two temperatures description ( the equation of state can be applied at the very end of the calculations ) .
it should be notice that the ` density matrix ' @xmath63 , defined in ( [ 22.3 ] ) , describes the ` closed - path motion ' in the functional space .
so , if = < n , out|e^is(_+)|m , in > o1and ^*= < q_1,q_2, ...
,p_n>= + = < m , in|e^-is(_-)|n , out > o2then , by definition , _+(_)=_-(_)= ( _ ) , o3with some ` turning - point ' fields @xmath64 , where @xmath65 is the remote hypersurface . the value of @xmath66 specifies the environment of the system . we will show that ( ) coincides with the kms boundary condition in some special cases . here
consequences of the vacuum boundary condition : (_)=0 03aare analyzed .
one should admit also that below boundary conditions are not unique : one can consider arbitrary organization of the environment of the considered system .
the @xmath0-matrix interpretation is able to show the way as an arbitrary boundary condition may be adopted .
this should extend the potentialities of the real - time finite - temperature field - theoretical methods .
it seems useful to illustrate the above microcanonical approach by the simplest example , see also . by definition ,
the @xmath4 particles production cross section _
n(s)=d_n(p)(q_1+q_2-^n_i=1p_i)|a_n(p , q)|^2 , 18where @xmath67 is the ordinary @xmath4 particle production amplitude in accelerator experiments .    considering the fourier transform of energy - momentum conservation @xmath68-function
one can introduce the generating function @xmath69 , see and references cited therein .
we may find in the result that @xmath70 is defined by the equality : _
n(e)=^+i_-ie^e_n ( ) ,   e=(q_1)+(q_2 ) , 3awhere _
n()= \{^n_i=1 } |a_n|^2=d_n ( ) .
19most probable value of @xmath62 in ( ) is defined by the equation of state ( ) . inserting ( ) into ( )
we find the expression ( ) if the momentum conservation shell is neglected .
the last one is possible since the cross sections are always measured in the definite frame .
let us consider the simplest example of noninteracting particles : @xmath71 where @xmath72 is the bessel function . inserting this expression into ( )
we can find that in the nonrelativistic case ( @xmath73 ) @xmath74 i.e. , we find the well known equality : e_kin = t , o4where @xmath75 is the mean kinetic energy and @xmath76 is the temperature ( boltzmann constant was taken equal to one ) .
it is important to note that the equation ( ) has a unique real solution @xmath77 rising with @xmath4 and decreasing with @xmath7 .
the expansion of the integral ( ) near @xmath77 unavoidably gives an asymptotic series with zero convergence radii since @xmath78 is the essentially nonlinear function of @xmath62 , see also sec.2.2.2 .
this means that , generally speaking , fluctuations in the vicinity of @xmath77 may be arbitrarily high and in this case @xmath77 has no physical sense .
but if the fluctuations are gaussian , then @xmath78 coincides with the partition function of the @xmath4 particle state and @xmath77 may be interpreted as the inverse temperature .
we will put the observation of this important fact in the basis of our thermodynamical description of the vhm region .
the notion of ` equilibrium ' over some parameter @xmath79 in our understanding is a requirement that the fluctuations in the vicinity of its mean value , @xmath80 , have a gaussian character .
notice , in this case , one can use this variable for a ` rough ' description of the system .
we would like to show now that the corresponding equilibrium condition would have the meaning of the correlations relaxation condition of bogolyubov , see also .
let us define the conditions when the fluctuations in the vicinity of @xmath81 are gaussian .
firstly , to estimate the integral ( ) in the vicinity of the extremum , @xmath81 , we should expand @xmath82 over @xmath62 :    _ n(+_c)=_n(_c)-+^2 _ n(_c)-^3_n(_c)+ ...
21cand , secondly , expand the exponent in the integral ( ) over , for instance , @xmath83 etc . in the result , if higher terms in ( ) are neglected , @xmath84-th term of the perturbation series _
n , k~\{}^k ( ) .
22ctherefore , because of euler s @xmath85 function , the perturbation theory near @xmath81 leads to the asymptotic series .
the supposition to define this series formally , for instance , in the borel sense is not interesting from physical point of view .
indeed , such a formal solution assumes that the fluctuations near @xmath81 may be arbitrarily high .
then , by this reason , the value of @xmath81 loses its significance : arbitrary values of @xmath86 are important in this case .
nevertheless it is important to know that our asymptotic series exists in some definite sense , i.e. we can calculate the integral over @xmath62 by expanding it over @xmath86 .
therefore , if the considered series is asymptotic , we may estimate it by first term if /^3_c
< < ( ^2_n(_c)/^2_c)^3/2 . 23cone of the possible solutions of this condition is /^3_c0 .
24cif this condition is satisfied , then the fluctuations are gaussian with dispersion @xmath87 see ( ) .    let us consider now ( ) carefully .
we will find by computing derivatives that this condition means the following approximate equality : -3 + 2 0 , 25cwhere @xmath88 means the @xmath84-th derivative .
for identical particles , ^(k)_n(_c)=n^k(-1)^kd_n(_c)^k_i=1(q_i ) , 26ctherefore , the left hand side of ( ) is the 3-point correlator @xmath89 since @xmath90 is a density of states for given @xmath62 : & k_3d_3(q)(<^3_i=1(q_i)>__c- 3<^2_i=1(q_i)>__c < ( q_3)>__c+ . + & + .2 ^ 3_i=1<(q_i)>__c ) , 27cwhere the index @xmath81 means that averaging is performed with the boltzmann factor @xmath91 .
notice , in distinction with bogolyubov , @xmath89 is the @xmath92 correlation function .
so , in our interpretation , one can introduce the notion of temperature @xmath93 if and only if the @xmath94 energy flows , measured by the corresponding correlation functions , are to die out .    as a result , to have all the fluctuations in the vicinity of @xmath81 gaussian , we should have @xmath95 , @xmath96 . notice , as follows from ( ) , the set of minimal conditions actually looks as follows : k_m<<k_2 ,  m3 . 27dif
the experiment confirms this conditions then , independently from the number of produced particles , the final state may be described with high enough accuracy by one parameter @xmath81 and the energy spectrum of particles is gaussian . in this conditions
one may return to the statistical and the hydrodynamical models .
considering @xmath81 as a physical ( measurable ) quantity , we are forced to assume that both the total energy of the system , @xmath97 , and the conjugate to it , variable @xmath81 , may be measured simultaneously with high accuracy .
we would like to show now that the ordinary big partition function of the statistical system coincide with _
n , m_(_1,z_1;_2,z_2 ) r_nm(p , q)=(,z ) , densmwhere @xmath98 is defined by ( ) . the summation and integration are performed with constraints that the mean energy of particles in the initial(final ) state is @xmath99 . on
may interpret @xmath100 in the first approximation as the temperature and @xmath101 as the activity for initial(final ) state .    direct calculation , see appendix b , gives the following expression for generating functional : ( , z)=e^-i*n*(_i^*_j)r_0 ( ) , 22.22awhere the particle number operator ( @xmath102 ) & * n*(_i^*_j)=- dxdx(_+(x)d_+-(x - x,_2,z_2)_-(x)- + & -_-(x)d_-+(x - x,_1,z_1)_+(x ) ) p22and r_0()=z(_+)z^ * (-_-) , p23where @xmath103 is defined in ( ): @xmath104 and , for the vacuum boundary condition @xmath105 , d_+-(x - x,,z)=-id_1(q)e^iq(x - x)e^-(q)z(q ) 22.23a + d_-+(x - x,,z)=id_1(q)e^-iq(x - x)e^-(q)z(q ) 22.24aare respectively the positive and negative frequency correlation functions at @xmath106 .
it is evident , r_nm(p , q)=. _
k=1^n\{e^_1(p_k ) } _ k=1^m\{e^_2(q_k ) } ( , z)|_z_i=0 .
p2notice , defining @xmath98 through the generating functional we extract the boltzmann factors @xmath107 since the energy - momentum conservation @xmath68-functions were extracted from amplitudes @xmath54 .
we suppose that @xmath103 may be computed perturbatively . as a result , ( @xmath108 is the variational derivative ) r(,z)= e^-iv(-i_+)+iv(-i_- ) e^ dx dxj_i ( x)d_ik(x - x;,z)j_k ( x ) , 22.29awhere @xmath109 is the matrix green function .
these green functions are defined on the mills time contours @xmath110 in the complex time plane ( @xmath111 ) , see fig.2 .
this definition of the time contours coincides with the keldysh time contour .
the generating functional ( ) has the same structure as the generating functional of niemi and semenoff .
the difference is only in the definition of green functions @xmath112 .
this choice is a consequence of the boundary condition ( ) .
so , if ( ) is used , then the green function is defined by eq .
notice also that if @xmath113 then a new green function obey kms boundary condition , see ( ) .
following niemi and semenoff one can write ( ) in the form : ( ) = d_nse^is_ns( ) , nisewhere the functional measure @xmath114 and the action @xmath115 are defined on the closed complex time contour @xmath116 , see fig.3 .
the choice of initial time @xmath117 and @xmath118 is arbitrary
. then one can perform shifts : @xmath119 and @xmath120 . in result ,
( * i * ) if @xmath113 , ( * ii * ) if contributions from imaginary parts @xmath121 and @xmath122 of the contour @xmath116 have disappeared in this limit , ( * iii * ) if the integral ( ) may be calculated perturbatively then this integral is a compact form of the representation ( ) .
notice that the requirements ( * i * ) - ( * iii * ) are the equivalent of the euclidean postulate of schwinger . in this frame one can consider another limit @xmath123 .
then the @xmath116 contour reduces to the matsubara imaginary time contour , fig.4 .
later on we will use this @xmath0-matrix interpretation of thermodynamics .
but one should take in mind that corresponding results will hide assumptions ( * i * ) - ( * iii * ) .
we would like to mention the ambivalent role of external particles in our @xmath0-matrix interpretation of thermodynamics . in the ordinary matsubara formalism
the temperature is measured assuming that the system under consideration is in equilibrium with the thermostat , i.e. temperature is the energy characteristics of @xmath124 particles . in our definition
the temperature is the mean energy of produced , i.e. non - interacting , particles .
it can be shown that both definitions lead to the same result .
explanation of this coincidence is the following .
let us consider the point of particle production as the coordinate of fictitious ` particle ' .
this ` particle ' interacts since the connected contributions into the amplitudes @xmath125 only are considered , and has the equal to produced particles momentum and so on .
the set of these ` particles ' form a system .
interaction among these ` particles ' may be described by corresponding correlation functions , see sec.2.3.3 .
let us consider now the limit @xmath123
. in this limit ( ) reduces to ( ) = d_me^-s_m( ) , matswhere the imaginary time measure @xmath126 and action @xmath127 are defined on the matsubara time contour . the periodic boundary condition ( )
should be used calculating integral ( ) .
the rules and corresponding problems in the integral ( ) can be calculated as described in many textbooks , see also .    in the limit
considered the time was eliminated in the formalism and the integral in ( ) performed over all states of the ` particles ' system with the weight @xmath128 .
notice , the doubling of degrees of freedom has disappeared and our fictitious ` particles ' became real ones .
on the other hand , the produced particles may be considered as the probes through which we measure the interacting fields . as was mentioned above
, their mean energy defines the temperature , if the energy correlations are relaxed .
if even one of the conditions ( * i * ) - ( * iii * ) is not satisfied then one can not reduce our @xmath0-matrix formalism to the imaginary time matsubara theory .
then one can ask : is there any possibility , staying in the frame of @xmath0-matrix formalism , to conserve the statistics formalism .
this question is discussed in the appendix c , where the wigner functions approach is applied .
it may be shown that the formalism may be generalized to describe the kinetic phase of the nonequilibrium process , where the temperature should have the local meaning .
the comparison with the ` local equilibrium hypothesis ' is discussed at the end of appendix c.      our further consideration will based on the model independent ( formal ) classification of asymptotics .      we will consider the generating function : t(s , z)=_n=1^n_max z^n_n(s ) ,  s=(p_1+p_2)^2>>m^2 ,   n_max=/m .
21.3this step is natural since the number of particles is not conserved in our problem .
so , the total cross section and the averaged multiplicity will be : _
tot(s)=t(s,1)=_n_n(s ) ,   = _ n n ( _ n(s)/_tot)= .t(s , z)|_z=1 . 21.4    at the same time , the inverse mellin transform gives _
n=.t(s , z)|_z=0= t(s , z)= e^(-nz+t(s , z ) ) .
21.5the essential values of @xmath129 in this integral are defined by the equation ( of state ) : n = zt(z , s ) .
21.6taking into account the definition of the mean multiplicity @xmath3 , given in ( ) , we can conclude that the solution of ( ) @xmath130 is equal to one at @xmath131 .
therefore , @xmath132 are essential in the vhm domain .
the asymptotics over n ( @xmath133 is assumed ) are governed by the smallest solution @xmath130 of ( ) because of the asymptotic estimation of the integral ( ): _ n(s)e^-nz_c(n , s ) .
21.7    let us assume that in the vhm region and at high energies , @xmath134 , there exist such a value of @xmath135 that we can neglect in ( ) the dependence on the upper boundary @xmath136 .
this formal trick with the thermodynamical limit allows to consider @xmath137 as the nontrivial function of @xmath129 for finite @xmath7 .
then , it follows from ( ) that z_c(n , s)z_s  at  n , 21.8where @xmath138 is the leftmost singularity of @xmath137 in the right half plane of complex @xmath129 .
one can say that the singularity of @xmath137 attracts @xmath135 if @xmath139vhm .
we will put this observation in the basis of vhm processes phenomenology .
we would like to underline once more that actually @xmath137 is regular for arbitrary finite @xmath129 if @xmath7 is finite . but
@xmath135 behaves in the vhm domain as if it is attracted by the ( imaginary ) singularity @xmath138 .
and just this @xmath135 defines @xmath70 in the vhm domain .
we want to note that actually the energy @xmath5 should be high enough to use such an estimation .
one can notice from the estimation ( ) that @xmath70 weakly depends on the character of the singularity .
therefore it is enough to classify only the possible positions of @xmath138 .
we may distinguish following possibilities : & ( * a * )  z_s= :  _ n < o(e^-n ) + & ( * b * )  z_s=1 :  _ n >
o(e^-n ) + & ( * c * )  1<z_s < :  _ n = o(e^-n ) , 21.9i.e . , following this classification , the cross section may decrease faster ( * a * ) , slower ( * b * ) , or as ( * c * ) an arbitrary power of @xmath140
. it is evident , if all these possibilities may be realized in nature , then we should expect the asymptotics ( * b * ) .    as was explained in sec.2.2.1
, @xmath70 has the meaning of the @xmath4 particle partition function in the energy representation
. then @xmath137 should be the ` big partition function ' . taking this interpretation into account ,
as follows from lee - yang theorem , @xmath137 can not be singular at @xmath141 .    at the same time , the direct calculations based on the physically acceptable interaction potentials give the following restriction from above : ( * d * )  _ n < o(1/n ) 21.2this means that @xmath70 should decrease faster than any power of @xmath142 .
it should be noted that our classification predicts rough ( asymptotic ) behavior only and did not exclude local increase of the cross section @xmath70 .
one may notice that -=z_c(n , s)+o(1/n ) .
21.10using thermodynamical terminology , the asymptotics of @xmath70 is governed by the physical value of the activity @xmath135 .
one can introduce also the chemical potential @xmath143 .
it defines the work needed for one particle creation , @xmath144 , where @xmath145 is the produced particles mean energy .
so , one may introduce the chemical potential if and only if @xmath146 and @xmath135 may be used as the ` rough ' variables .
then the above formulated classification has a natural explanation .
so , ( * a * ) means that the system is stable with reference to particle production and the activity @xmath135 is the increasing function of @xmath4 , the asymptotics ( * b * ) may realized if and only if the system is unstable . in this case
@xmath135 is the decreasing function .
the asymptotics ( * c * ) is not realized in equilibrium thermodynamics
.    we will show that the asymptotics ( * a * ) reflects the multiperipheral processes kinematics : created particles form jets moving in the cm frame with different velocities along the incoming particles directions , i.e. with restricted transverse momentum , see sec.3.1.1 .
the asymptotics ( * b * ) assumes the condensation - like phenomena , see sec.3.3
. the third type asymptotics ( * c * ) is predicted by stationary markovian processes with the pqcd jets kinematics , see sec.3.2.2 .
the dis kinematics may be considered as the intermediate , see sec.3.2.1 .
this interpretation of classes ( ) allows to conclude that we should expect reorganization of production dynamics in the vhm region : the soft channel ( * a * ) of particle production should yield a place to the hard dynamics ( * c * ) , if the ground state of the investigated system is stable with reference to the particle production .
otherwise we will have asymptotics ( * b * ) .
let us consider the system with several correlation scales .
for example , in statistics one should distinguish correlation length among particles ( molecules ) and correlation length among droplets if the two - phase region is considered . in particle physics
, one should distinguish in this sense correlation among particles produced in result of resonance decay and correlations among resonances .
in pqcd one may distinguish correlations of particles in jet and correlation among jets .
there exist many model description of this physical picture . in statistics
mayer s group decomposition is well known . in particle physics one should note also the many - component formalism . we will consider the generating functions ( functionals ) formalism considering mostly jet correlations . in many respects
it overlaps the above mentioned approaches .
the generating function @xmath137 may be written in the form : t(z , s)=_k=1^c_k(s)= _ l=1^z^lb_l , g1where the coefficients @xmath147 are the moments of the multiplicity distribution p_n(s)=_n(s)/_tot(s ) .
g2so , c_1(s)=_n np_n(s)= ,  c_2(s)=_n n(n-1)p_n(s)-^2 g3and so on . using the connection with the inclusive distribution functions
@xmath148 : t(z , s)=_k=1^d_k(q ) f_k(q_1,q_2, ... ,q_k;s ) , g4it is easy to find that & c_1(s)=d_1(q ) f_1(q;s)=|f_1(s ) , + & c_2(s)=d_2(q ) \{f_2(q_1,q_2;s)-f_1(q_1;s)f_1(q_2;s)}=    _ \{k}_l=0^(_i=1^lk_i - k ) _ i=1^l\ { } , g6where @xmath149 and @xmath150 one may invert formulae ( ): |f_k(s)= _ \{n_k}=0^(_k=1^kn_k - l ) _
k=1^()^n_k g7    the mayer s group coefficients @xmath151 in ( ) have the following connection with @xmath147 : b_l(s)=_k=0^c_k+l(s ) .
g23    it seems useful to illustrate the effectiveness of the generating function method by the following example .
we will consider the transformation ( multiplicity @xmath152 activity @xmath129 ) to show the origin of the koba - nielsen - olesen scaling ( kno - scaling ) .
if @xmath153 , @xmath154 , then @xmath70 is described by the poisson formulae : _
k8it corresponds to the case of absence of correlations .
let us consider more weak assumption : c_m(s)=_m(c_1(s))^m , k9where @xmath155 is the energy independent constant , see also , where a generalization of kno scaling on the semi - inclusive processes was offered .
then t(z , s)=_m=1\{(z-1)|n(s)}^m .
k10to find the consequences of this assumption , let us find the most probable values of @xmath129 .
the equation of state @xmath156 has solution @xmath157 increasing with @xmath4 since @xmath137 is an increasing function of @xmath129 , if and only if , @xmath137 is nonsingular at finite @xmath129 . as was mentioned above , the last condition has deep physical meaning and practically assumes the absence of the first order phase transition .
let us introduce a new variable : = ( z-1)|n(s ) .
k12the corresponding equation of state looks as follows : = ( 1+)t ( ) .
k13so , with @xmath158 accuracy , one can assume that _
c(n/|n(s ) ) . k14are essential .
it follows from this estimation that such scaling dependence is rightful at least in the neighborhood of @xmath106 , i.e. in vicinity of main contributions into @xmath159 .
this gives : |n(s)_n(s)=_tot(s)(n/|n(s ) ) , k15where ( n/|n(s))t(_c(n/|n(s)))\{n/|n(s)_c(n/|n(s))}o(e^-n ) k16is the unknown function .
the asymptotic estimation follows from the fact that @xmath160 should be a nondecreasing function of @xmath4 , as follows from nonsingularity of @xmath137 .    the estimation ( ) is right at least at @xmath161 .
the range validity of @xmath4 , where solution of ( ) is acceptable , depends from exact form of @xmath137 . indeed , if @xmath162 , @xmath163 , then ( ) is right at all values of @xmath4 and it is enough to have the condition @xmath161 .
but if @xmath164 , @xmath163 , then ( ) is acceptable if and only if @xmath165 .
representation ( ) shows that just @xmath38 is the natural scale of multiplicity @xmath4 .
this representation was offered first as a reaction on the so called feynman scaling for inclusive cross section : f_k(q_1,q_2, ... ,q_k)~^k_i=1 .
k17    as follows from estimation ( ) , the limiting kno prediction assumes that @xmath166 . in this regime @xmath137
should be singular at @xmath167 .
the normalization condition @xmath168 gives : @xmath169 , where @xmath170 is the constant .
notice , such behavior of the big partition function @xmath137 is natural for stationary markovian processes described by logistic equations . in the field theory such equation describes the qcd jets .
we wish to generalize expansion ( ) to take into account the possibility of many - component structure of the multiple production processes .
let us consider particle production through the generation , for instance , of jets . in this case decay of a particle of high virtuality @xmath171 forms a jet of lower virtuality particles .
it is evident that one should distinguish correlation among particles in the jet , and correlation among jets .
let @xmath172 be the @xmath173 that @xmath174-th jet of mass @xmath175 includes @xmath176 particles , @xmath177 , where _
i=1^n_jn_i = n g8the jets are the result of particles decay .
then let us assume that @xmath178 defines the mean number of jets of mass @xmath175 and momentum @xmath179 : & ^(1)_n()=_n_j_\{n}_n_j ( _
i=1^n_jn_i - n ) + & _ i=1^n_j \ { e^-(p_i)|n_1(m_i , p_i)_n_i(m_i ) } , g9where @xmath180 . notice the boltzmann factor @xmath107 , where @xmath181 is the jets energy , play the same role as the corresponding factor in ( ) and introduced to take into account the energy conservation law .
we consider the vhm domain and for this reason @xmath182 is assumed .
it is useful to avoid the particles number conservation law ( ) . for this purpose
we will introduce ^(1)(,z)=_n^(1)_n()=\ { e^-(p)|n_1(m , p)(t(z , m)-1 ) } , g10where t(z , m)=_n z^n_n(m ) . g11comparing ( ) with ( )
we may conclude that @xmath183 plays the role of activity of jets .
then the generalization is evident : & ( , z)=_k^(k)(,z)= + & = \{_k_i=1^k\{dm_id_1(p_i ) e^-(p_i)(t(z , m_i)-1)}. + & .|n_k(m_1,p_1, ... ,m_k , p_k ) } , g12where @xmath184 has the same meaning as @xmath147 , i.e. @xmath184 is the correlation function of @xmath84 jets .
let us consider the following ` bootstrap ' regime when @xmath185 is defined by the equation : ( , z ) e^-(p)|n_1(m , p)t(z , m ) .
g13inserting here the strict expression ( ) we find a nonlinear equation for @xmath183 .
the solution of ( ) assumes that e^-(p)|n_1 t > > \{e^-(p ) |n_1t}^2 g14and e^-|n_1 t > > ( e^-t ) ( e^-t ) |n_2
g15    to solve eq .
( ) in the vhm region , where the leftmost singularity over z is important , let us consider the @xmath186 : t(z , m)= ,  _ 0>0 , g16where @xmath187 is the polynomial function of @xmath129 , @xmath188 . using the normalization condition : |n_j=.t(z , m)|_z=1 g17we can find : a(m)_0=|n_j-(1,m ) ,  (1,m).(z , m)|_z=1 . g18
the partition function of the jet @xmath183 defined by @xmath186 ( ) is singular at z_s(m)=1 +
. g19this singularity would be significant in the vhm region if @xmath189 is decreasing function of @xmath52 .
this means an assumption that @xmath190 so , in first approximation we will choose a(m)=|n_j(m ) .
g20this choice may be confirmed by concrete model calculations .
taking into account the energy conservation law , conditions ( ) and ( ) are satisfied if \{-n } <
g21at @xmath139vhm .
therefore , ( ) obey eq .
( ) with exponential accuracy in the vhm region , i.e. if n > > = .
g22we assume here that one can find so large @xmath4 and @xmath7 that with exponential accuracy the factors @xmath191 did not play an important role .
but at low energies the condition ( ) is important and the factors @xmath191 should be taken into account .
notice now important consequence of our ` bootstrap ' solution : it means that we can leave production of the heavy jets only , if @xmath192 . on the other hand ,
let us choose @xmath193 , where @xmath194 is the function of @xmath7 , and consider @xmath195 .
then the condition ( ) defines @xmath196 : if z_0 > > g22then we are able to obey the inequalities ( ) and ( ) .
the jet mean multiplicity , see sec.3.2.2 , |n_j(s)~. g26then = \{e^/-1}^-1~<<|n_j(s ) .
g50at @xmath195 .
therefore , ( ) may be satisfied outside the vhm domain .
let us compare now the solutions of the equation of state . inserting ( ) into ( )
we can find for a jet of mass @xmath5 that z^1_c = z^1_s-=1 + - .
g24the two - jet contribution of the masses @xmath197 gives : z^2_c = z^2_s-=1 + - . g25at arbitrary finite energies @xmath198 and , as follows from ( ) , they decrease @xmath199 with energy .    noting the normalization condition , @xmath200 , and assuming that the vacuum is stable , i.e. @xmath201
, we can conclude that     if @xmath202 then @xmath138 attracts @xmath130 , i.e. @xmath203 , and if @xmath204 then this contributions should be significant in @xmath159 ;     if @xmath195 , then @xmath204 , and if @xmath4 satisfy the inequality ( ) , or if @xmath196 satisfy the inequality ( ) , then considered contributions are significant in @xmath159 .
+ it is the ( energy - multiplicity ) asymptotics equivalence principle .
one of the simplest consequences of this principle is the prediction that the mean transverse momentum of created particles should increase with multiplicity at sufficiently high energies .
this principle is the consequence of independence of contributions in the vhm domain on the type of singularity in the complex @xmath129 plane and of the energy conservation law .
just the last one shifts the two - jet singularity to the right and @xmath205 .
we would like to notice also that this effect , when the mostly ` energetic population ' survives was described mathematically by v.volterra .
intuitively evident is that one may find the ` energetic population ' searching the vhm one , or , it is the same , giving it a rich supply , i.e. to give the population enough energy .
this is our ( energy - multiplicity ) equivalence ( ( @xmath206)-equivalence ) principle .
notice , if the amount of supply is too high then few populations may grow .
this is the case when the difference @xmath198 tends to zero at high energies .
we would like to note that singular at finite @xmath129 partition functions was predicted in the @xmath207-theory , in qcd jets , in the generalized bose - einstein distribution model . in all of this models decay of the essentially nonequilibrium initial state ( highly virtual
parton , heavy resonance , etc ) was described .
one may distinguish the phases of the media by a characteristic correlation length .
then the phase transition may be considered as the process of changing correlation length .
our ` bootstrap ' solution predicts just such phenomena : at low multiplicities the long - range correlations among light jets is dominant .
the ` bootstrap ' solution predicts that for vhm processes just the short - range correlations among particles of the heavy jet become dominant .
the ( @xmath206)-equivalence means that this transition is a pure dynamical effect .
multiple production phenomena was first observed more than seventy years ago . during this time
a vast experimental information was accumulated concerning hadron inelastic interactions , see the review papers .
now we know that at high energy @xmath208 : + ( i ) the total cross section @xmath209 of hadron interactions is enhanced almost completely by the inelastic channels ; + ( ii ) the mean multiplicity of produced hadrons @xmath3 slowly ( logarithmically ) grows with @xmath208 ; + ( iii ) the interaction radii of hadrons @xmath210 slowly ( logarithmically ) increase with energy . ; + ( iv ) the multiplicity distribution @xmath19 is wider than the poisson distribution . ; + ( v ) the mean value of transverse momentum @xmath84 of produced hadrons is restricted and is independent of the incident energy @xmath5 and produced particle multiplicity @xmath4 ; + ( vi ) the one - particle energy spectrum @xmath211 .
first of all , ( i ) means that the high energy hadron interaction may be considered as the ordinary dissipation process . in this process the kinetic energy of incident particles is spent in produced particle mass formation .
the vhm process takes place in the vacuum and then it was assumed on the early stages that the multiple production phenomena reflects a natural tendency of the excited hadron system to get to equilibrium with the environment . in this way one can introduce as a first approximation the model that the excited hadron system evolves without any restrictions . in this model
we should have @xmath212 .
the dissipation is maximal in this case and the entropy @xmath8 exceeds its maximum .
this simple model has definite popularity up to 70-th .
but the experimental data ( ii ) and ( iii ) prohibit this model and it was forgotten .    choosing the model we would like to hope that the considered model    takes into account experimental conditions ( i ) - ( vi ) in the @xmath15 domain ;     has natural asymptotics over multiplicity to the vhm region .
+ it is necessary to remember also that     new channels of hadron production may arise in the vhm region .
+ it is impossible to understand all possibilities without those offered in sec.2.3.2 in the classification of asymptotics .
thus , we will observe predictions of     _ multiperipheral models _
, + distinguishing the soft pomeron models , see .     _ the dual - resonance model _
+ predictions for the vhm region are described also .
it can be shown that this models predict asymptotics ( a ) if @xmath213 .
just this result explains why vhm domain is defined by the condition ( ) .
we are forced since it allows to include pqcd , forbidden by the multiperipheral models , considering     _ hard pomeron model _
+ production of mini - jets .
but we will find using monte carlo simulations that the pqcd pomeron is unable to adopt the hard channels of hadron production .
then     _ the deep inelastic processes _
+ for vhm region will be considered to generalize the dglap kinematics in the case of heavy qcd jets production .
the analysis shows that transition to the vhm leads to the necessity to include low-@xmath214 sub - processes . as a result
we get out of the range of pqcd validity .
 _ multiple reduction of jets . _
+ we will see , that at very high energies in the vhm region the heavy jets creation should be a dominant process if the vacuum is stable with reference to the particle production .
we will consider also decay of the ` false vacuum ' to describe the consequence of     _ phase transition _ + in the vhm domain .
this channel is hardly seen for @xmath15 since the confinement constraints may prevent cooling of the system up to phase transitions condition .
later on multiple production physics was developed on the basis of experimental observation ( iv ) .
the regge pole model naturally explains this experimental data and , at the same time , absorbs all experimental information , ( i)-(vi ) . at the very beginning ,
adopting the regge poles notion without its microscopical explanation , this description was self - consistent .
the efforts to extend the regge pole model to the relativistic hadron reactions was ended by the reggeon diagram technique , see and references cited therein , and it was used later to construct the perturbation theory for @xmath70 .
it was shown that the multiplicity distribution is wider then the poissonian because of the multi - pomeron exchanges .
the leading energy asymptotics pomeron contribution reflects the created particle kinematics described in appendix d , where the available kinematical scenario in the frame of pqcd are described .
so , the longitudinal momentum of produced particles is large and is strictly ordered .
at the same time , particles transverse momentum is restricted .
let us consider the inelasticity coefficient introduced in ( ) @xmath215 , where @xmath9 is the energy of the fastest particle in the laboratory frame .
then the strict ordering of particles in the pomeron kinematics means that @xmath216 is independent of the index of the particle .
so , if the fastest particle has the energy @xmath217 , then the following particle should have the energy @xmath218 , and so on .
following to this law , @xmath219-th particle would have the energy @xmath220 . in the laboratory frame
the energy should degrade to @xmath221 . inserting here the above formulated estimation of @xmath222
we can find that if the number of produced particles is n_0 ,  n_0=-(1-)^2>0 ,  = ( s / m^2 ) , 1.1athen we may expect the total degradation of energy .
this degradation is the necessary condition noting that the total cross section of slowly moving particles may depend only slightly on energy and is seems necessary for natural explanation of the weak dependence of the hadron cross sections on the energy
. this consideration would be lorentz - covariant if one can find the slowly moving particle in an arbitrary frame .
resulting estimation of mean multiplicity have good qualitative experimental confirmation .
notice that it was assumed deriving ( ) that the energy degrades step by step .
in other words , if we introduce a time of degradation , then the time @xmath223 is needed for complete degradation of energy .
assuming the random walk in the normal to incident particle plane , we can conclude that the points of particle production are located on a disk ( in the moving frame ) of radii @xmath224 .
this means that the interaction radii should grow with the energy of the colliding particles .
if @xmath225 is the cross section to observe particle @xmath226 inclusively in the @xmath227 and @xmath228 particles collision , then it was found experimentally that the ratio = = 1.2is universal .
this may be interpreted as the direct evidence of fact that the hadron interactions have a large - distance character , i.e. that the interaction radii should be large .
this picture assumes that the probability to have total degradation of energy is ~e^-b^2/4 , 1.3where @xmath229 is some dimensional constant ( the slope of regge trajectory ) and @xmath230 is the 2-dimensional impact parameter .
this formulae has also the explanation connected to the vacuum instability with reference to the real particle production in the strong colour electric field .
the above picture has natural restrictions .
we can assume that each of the produced particles may be the source of above described @xmath42-channel cascade of the energy degradation .
this means that in the frame of the pomeron phenomenology , we are able to describe the production of n <
^2 1.4particles only . if @xmath231 then the density of particles in the diffraction disk becomes large and ( a ) one should introduce short - distance interactions , or ( b ) rise interaction radii .
it will be shown that just ( a ) is preferable .
we will build the perturbation theory in the phenomenological frames ( i ) - ( vi ) .
considering the system with variable number of particles the generating function t(z , s)=_n z^n _ n(s ) i1would be useful .
one can use also the decomposition : t(z , s)=_tot(s)\{(z-1)c_1(s)+(z-1)^2c_2(s)+ ... } , i2where , by definition , c_1(s)==.t(z , s)|_z=1 i3is the mean multiplicity , c_2(s)=.t(z , s)|_z=1 i4is the second binomial momentum , and so on .    our idea is to assume that all @xmath232 , @xmath154 may be calculated perturbatively choosing p(0,s)=e^(z-1 ) i5as the born approximation ` superpropagator ' .
it is evident that ( ) leads to poisson distribution .
then , having in mind ( ii ) , ( iii ) and ( v ) we will use following @xmath186 : p(q , s)=e^(o)-q^2se^(z-1 ) , i6where the transverse momentum @xmath233 is conjugate to the impact parameter @xmath230 .
so , the born term ( ) is a fourier transform of the simple product of ( ) and ( ) .
it contains only one free parameter , the pomeron intercept @xmath234 . on the phenomenological level
it is not important to know the dynamical ( microscopical ) origin of ( ) .    for our purpose
the laplace transform of @xmath235 would be useful . if @xmath236 , then (  , q^2)=_0^de^-p(q^2,s)= .
i8it is the propagator of two - dimensional field theory with mass squared @xmath237 knowing the gribov s reggeon calculus completed by the abramovski - gribov - kancheli ( agk ) cutting rules one can investigate the consequences of this approach .
the lla approximation of the pqcd gives = ( 0)-1=_s0.55 ,  _ s=0.2 .
i10but radiative corrections give @xmath238 .
we will call this solution as the bfkl model .
the quantitative origin of the restriction ( ) is following .
the contribution of the diagram with @xmath239 pomeron exchange gives , since the diffraction radii increase with @xmath7 , see ( ) , mean value of the impact parameter decreasing with @xmath239 : @xmath240 where @xmath241 . on the other hand , the number of necessary pomeron exchanges @xmath242 since one pomeron gives maximal contribution ( with factorial accuracy ) at @xmath243 . in result , |^2~a. i11therefore , if the transverse momentum of created particles is a restricted quantity , i.e. @xmath244 , where @xmath245 is a constant , then the mechanism of particle production is valid up to n~^2 .
o8    following our general idea , it will be enough for us find the position of singularity over @xmath129 .
analysis shows that ( ) predict the singularity at infinity .    in appendix e the gribov s reggeon diagram technique with cut pomerons is described and ( ) is derived .
it can be shown using this technique that the model with the critical pomeron , @xmath246 , is inconsistent from the physical point of view .    as was mentioned above
, the model with @xmath247 is natural for the pqcd .
the concrete value of @xmath248 will not be important for us .
we will assume only that 0<0@xmath249>0@xmath250|^2~b^24(+(n/))@xmath251z@xmath252n~n_max~s@xmath253e^-n@xmath254~s^-1/2@xmath255(m)@xmath256m@xmath257m@xmath258^r(m)@xmath259^r=1/2@xmath260^r_n(m)@xmath261m@xmath262n@xmath263|n^r(m)@xmath264^r_n(m)@xmath265|n^r(m)@xmath266m
/ m_0@xmath267n@xmath268n|n^2(s)= |n_0^r ( /m_0)@xmath269|n(s)<<n<|n^2(s)@xmath270q|q@xmath271n@xmath272d_ab(x , q^2;n)@xmath273d_ab(x , q^2)@xmath274b@xmath275q^2<0@xmath276a@xmath277~@xmath278>>@xmath279_s()<<1@xmath280q^2@xmath279x@xmath281(1/x)>>1@xmath282d_ab(x , q^2)@xmath283(1/x)@xmath284|q^2|@xmath285~(1/x)/|q^2|@xmath286@xmath287@xmath288 ^ 2<<k_i^2<<-q^2@xmath289k_i^2>0@xmath290~(1/)@xmath291k_i^2>>^2@xmath292|n_j@xmath293|k|@xmath294|n_j @xmath295@xmath296||k_i|@xmath297|k_i|@xmath298k_i@xmath299x@xmath300x0@xmath301t_ab(x , q^2;z)@xmath302n@xmath303z_c@xmath304n@xmath305z_c@xmath306z_c(x , q^2;n)@xmath307n@xmath308d_ab(x , q^2;n)@xmath309t_ab(x , q^2;z)@xmath310z=1@xmath311d_ab(x , q^2;n)@xmath312c_ab(x , q^2;n)@xmath313n@xmath314c_ab(x , q^2;n)~n@xmath315c_ab(x , q^2;n)=const.@xmath316n@xmath317nn_p@xmath318n@xmath319c_ab(x , q^2;n)@xmath320|_ab(x , q^2;n)@xmath321|n_j()@xmath322@xmath323@xmath324n>>|n_g@xmath325(,z)@xmath326x0@xmath327_s(-q^2/^2)~1@xmath279_s(1/x)~1@xmath328|q^2|la^2@xmath329@xmath330_s(^2)<<1@xmath331(^3)_6@xmath332|q^2|la^2@xmath333z_s-1~(1/|n)@xmath334a^(1j)(n , m)@xmath335n@xmath336|n_j(m)@xmath337m@xmath338c_j@xmath339n/|n@xmath340m@xmath341m_1@xmath279m_2@xmath342m_1=m_2=m/2@xmath343m_2<<m_1m@xmath344a^(1j)@xmath345n_1+n_2=n@xmath346@xmath347@xmath348|n_j(m)@xmath349m@xmath350a^1j@xmath336a^2j@xmath351z@xmath352s@xmath353t(z , s)@xmath354z=1@xmath355z_c@xmath356s@xmath279n@xmath357t@xmath358s@xmath359s@xmath360n@xmath361n@xmath362_n(s)@xmath363_n(s)@xmath364n@xmath365w_n(x_1,x_2, ... ,x_n)@xmath366x_k=(u , q)_k@xmath367s@xmath368s@xmath369(|n(s)/n)<<1@xmath370|n(s)@xmath371n@xmath372s@xmath373_max@xmath374(_max / s)0@xmath375s@xmath376n@xmath377n_max=/m@xmath378@xmath336t(z , s)@xmath379z@xmath380z>1@xmath381t(z , s)@xmath382n_max@xmath383t(z , s)@xmath384n@xmath385z@xmath386t(z)@xmath387z@xmath388t(z)@xmath389z=1@xmath390z@xmath391n@xmath392z@xmath393n@xmath394n=@xmath395z=1@xmath396t(z)@xmath397z=1@xmath398z=@xmath399z@xmath400_n@xmath401_n\{-cn^}@xmath402t(z)@xmath403z=1@xmath4041 , 1.3awhere @xmath130 is smallest solution of the equation of state n = zt(z ) .
1.2a    therefore , to have the singularity at @xmath106 , we should consider @xmath405 as a decreasing function of @xmath4 . on the other hand , at constant temperature
, @xmath406 is the chemical potential , i.e. is the work necessary for creation of one particle .
so , the singularity at @xmath106 means that the system is unstable : the less work is necessary for creation of one more particle if @xmath407 is a decreasing function of @xmath4 .
the physical explanation of this phenomena is the following , see also .
the generating function @xmath408 has following expansion : t(z)=\{_lz^lb_l } , 1.4awhere @xmath151 are known as the mayer s group coefficients .
they can be expressed through the inclusive correlation functions and may be used to describe formation of droplets of correlated particles , see sec.2.3.3 .
so , if droplet consist of @xmath409 particles , then b_l~e^-l^(d-1)/d 1.5is the mean number of such droplets .
here @xmath410 is the surface energy of @xmath411-dimensional droplet .
inserting this estimation into ( ) , t(z)~_l e^(l - l^(d-1)/d ) ,  = z. 1.6the first term in the exponent @xmath412 is the volume energy of the droplet and being positive it tries to enlarge the droplet .
the second surface term @xmath413 tries to shrink it .
therefore , the singularity at @xmath106 is the consequence of instability : at @xmath132 the volume energy abundance leads to unlimited growth of the droplet .    in conclusion
we wish to formulate once more the main assumptions .
it was assumed first of all that the system under consideration is in equilibrium .
this condition may be naturally reached in the statistics , where one can wait the arbitrary time till the system becomes in equilibrium .
note , in the critical domain , the time of relaxation @xmath414 , @xmath415 , @xmath416 , @xmath417 is the critical temperature
.    we can not give the guarantee that in the high energy hadron collisions the final state system is in equilibrium .
the reason of this uncertainty is the finite time the inelastic processes and presence of hidden ( confinement ) constraints on the dynamics .
but if the confinement forces are frozen in the vhm domain , i.e. the production process is ` fast ' , then the equilibrium may be reached .
we may formulate the quantitative conditions , when the equilibrium is satisfied .
one should have the gaussian energy spectra of created particles .
if this condition is hardly investigated in the experiment , then one should consider the relaxation of ` long - range ' correlations .
this excludes the usage of relaxation condition for the ` short - range ' ( i.e. resonance ) correlation    ( ii ) .
the second condition consists in the requirement that the system should be in the critical domain , where the ( equilibrium ) fluctuations of the system become high . having no theory of hadron interaction at high energies we can not define where lies the ` critical domain ' and even whether it exists or not .
but , having the vhm ` cold ' final state , we can hope that the critical domain is achieved .
the quantitative realization of this picture is given in appendix j. it is important to note that used there the semiclassical approximation is rightfully in the vhm domain .
it seems useful to start the discussion of models outlined the main problems , from authors point of view .
* a. * _ soft colour parton problems _
the infrared region of soft colour parton interactions is a very important problem of high energy hadron dynamics .
such fundamental questions as the infrared divergences of pqcd , collective phenomena in the coloured particles system and symmetry breaking are the phenomena of the infrared domain .
the standard ( most popular ) hadron theory considers pqcd at small distances ( in the scale of @xmath418gev ) as the exact theory .
this statement is confirmed by a number of experiments , namely deep - inelastic scattering data , hard jets observation .
but the pqcd predictions have a finite range of validity since the non - perturbative effects should be taken into account at distances larger then @xmath419 .
it is natural to assume , building the complete theory , that at large distances the non - perturbative effects @xmath420 the perturbative ones . as a result pqcd loses its predictability @xmath421 by the non - perturbative effects .
notice , the pqcd running coupling constant @xmath422 becomes infinite at @xmath423 and we do not know what happens with pqcd if @xmath424 .
there are few possibilities .
for instance , there is a suspicion that at @xmath425 the properties of theory changed so drastically ( being defined on a new vacuum ) that even the notions of pqcd @xmath426 .
this means that pqcd should be truncated from below on the ` fundamental ' scale @xmath427 .
it seems natural that this infrared cut - off would influence the soft hadrons emission .
the new possibility is described in appendix k. this strict formalism allows to conclude that pure pqcd contributions are realized on zero measure , i.e. it is the phenomenological theory only .
the successive approach shows that the yang - mills theory should be described in terms of ( action , angle)-like variables .
the last one means that the self - consistent description excludes such notions as the ` gluon ' . as a result of
this substitution new perturbation theory would be free from infrared divergences , i.e. there is not necessity to introduce the infrared cutoff parameter @xmath427 .
( moreover , in the sector of vector fields ( without quarks ) the theory is ultraviolet stable . )
it seems important for this reason to investigate experimentally just vhm events , where the soft colour partons production is dominant .
it is important to try to raise the role of pqcd in the ` forbidden ' area of large distances .
the vhm processes are at highly unusual condition , where the non - perturbative effects must be negligible .
* b*. _ dissipation problems _    the highly nonequilibrium states decay ( thermalization ) which means in the pqcd terms that the process of vhm formation should be enhanced , at least in asymptotics over multiplicity and energy , by jets .
it is the general conclusion of nonequilibrium thermodynamics and it means that the very nonequilibrium initial state tends to equilibrium ( thermalized ) as fast as possible .
the entropy @xmath428 of a system is proportional to the number of created particles and , therefore , @xmath428 should tend to its maximum in the vhm region .
but the maximum of entropy testifies also the equilibrium of the system .
* c*. _ collective phenomena _
we should underline that the collective phenomena may take place if and only if the particles interaction energy is comparable to the kinetic one .
the vhm system considered may be ` cold ' and ` equilibrium ' . for this reason the vhm state
is mostly adopted for investigation of collective phenomena .
one of possible states in which the collective effects , see a.g . , may be important is the ` could coloured plasma ' .
the fundamental interest presents the problem of vacuum structure of yang - mills theory .
for instance , if the process of cooling is ` fast ' , since the dissipation process of vhm final state formation should be as fast as possible , then one may consider formation of vacuum domains with various properties
. then decay of these domains may lead to large fluctuations , for instance , of the isotopic spin .
another important question is the collective phenomena in the vhm final state .
the last one may be created ` perturbatively ' , for instance , by the formation of heavy jets . then the colour charges should be confined .
there are various predictions about this process .
one of them predicts that there should be a first order phase transition .
now we can ask : what can models say concerning the above problems ?    * a*. _ soft coloured partons production _
the multiperipheral models predict fast decreasing of topological cross section in the vhm domain @xmath429 , @xmath430 . at the same time
the mean transverse momentum should @xmath431 in this domain since the interaction radii should ` increase ' with @xmath4 .
the bfkl pomeron predicts the same asymptotics , @xmath430 , but the pqcd jet predicts @xmath166 .
the @xmath22 attempt to insert into the bfkl pomeron the production of particles via ( mini)jets seems impossible .
this ` insertion ' can be done into the dis ladder but investigation of the lla kinematics in the vhm domain allows the conclusion that the ` low - x ' contributions should be taken into account .
all this experience allows the assumption that in the vhm domain no ` @xmath42-channel ladder ' diagrams play sufficient role .
this means the existence of a transition to the processes with jet dominance .
the pqcd is unable to predict the transition mechanism .
* b*. _ transition into equilibrium _
if the @xmath42-channel ladders are ` destroyed ' in the vhm region , then jets , despite the small factor @xmath432 in the cross sections , are the only mechanism of particle production in the vhm domain .
dominance of heavy jets in the vhm domain may naturally explain the tendency towards equilibrium .
but the description of thermalization in terms of jets of massless gluons production destroy this hope : the jet contribution @xmath166 assumes ` bremsstrahlung ' of @xmath433 gluons .
this prevents the equilibrium since ordering without fail introduce the non - relaxing correlations .
* c*. _ collective phenomena _
considering the collective phenomena , we proposer to distinguish ( a ) the collective phenomena connected with the vacuum and ( b ) the collective phenomena produced in the vhm system . following the experience of sec.3.3 we can conclude that the signal of vacuum instability is inequality : @xmath434 .
the case ( b ) will not effect the cross section @xmath70 .
but if the system became equilibrium in the vhm domain then the collective phenomena may be investigated using ordinary thermodynamical methods . for instance , noting that @xmath435 is the free energy one can measure the thermal capacity ( _ c , z)=c(_c , z ) . o23then
, comparing capacities of hadron and @xmath436-quanta systems we can say whether or not the phase transition happened .
the connection of the equilibrium and relaxation of correlations is well known .
continuing this idea , if the vhm system is in equilibrium one may assume that the colour charges in the pre - confinement phase of vhm event form the plasma .
one should note here that the expected plasma is ` cold ' and ` dense ' .
for this reason no long - range confinement forces would act among colour charges . then , being ` cold ' , in such a system various , collective phenomena may be important .
the experimental possibilities in the vhm domain are not clear up to now .
nevertheless first steps toward formulation of trigger system was done , see . by this reason we would like to restrict ourselves by following two general questions .
it seems that these questions are mostly important being in the very beginning of vhm theory .
0.3 cm * i*. _ for what values of multiplicity at a given energy the vhm processes become hard _ ?
+ the answer to this question depends on the value of the incident energy .
if we know the answer then it will appear possible to estimate     the role of multiperipheral contributions ,     the jet production rate ,     the role of vacuum instabilities .
it seems that the experimental answer to this question is absent since produced particles are soft , theirs mean transverse and longitudinal momenta have the same value . in our understanding
this question means : the total transverse energy may be extremely high .
it is interesting also to search the heavy jets , i.e. to observe the fluctuations of particle density in the event - by - event experiment , but this program seems vague since , for all evidence , fractal dimensions tend to zero with increasing multiplicities .    0.3 cm * ii*. _ for what values of multiplicity does the vhm final state become equilibrium _ ?
+ we hope that having answer on this question we would be able     to investigate the status of pqcd ,     to observe the phase transition phenomena directly ,     to estimate the role of confinement constraints .
the equilibrium means that the energy correlation functions mean values are small .
it is interesting also , for example , the charge equilibrium , when the mean value of charge correlation functions are small .
notice , the effect of the phase space boundary may lead to ` equilibrium ' . indeed ,
if @xmath437 and @xmath438 then one may neglect the momentum dependence of the amplitudes @xmath439 . in this case
the momentum dependence is defined by the boltzmann exponent @xmath440 , @xmath441 , only and we get naturally to the gaussian law for momentum distribution .
correlators should be small in this case since there is no interactions among particles ( @xmath439 are constants ) .
but our question assumes that we investigate the possibility of equilibrium when @xmath133 and @xmath442 .
there are various approaches to build the real - time finite - temperature field theories of schwinger - keldysh type ( e.g. ) .
all of them use various tricks for analytical continuation of imaginary - time matsubara formalism to real time .
the basis of the approaches is the introduction of the matsubara field operator _m ( * x * , ) = e^h_s ( * x * ) e^-h , 32where @xmath443 is the interaction - picture operator introduced instead of the habitual heisenberg operator @xmath444 eq .
( ) introduces the averaging over the gibbs ensemble instead of averaging over zero - temperature vacuum states .
if the interaction switched on at the moment @xmath117 adiabatically and switched off at @xmath118 then there is the unitary transformation : (x , t)=u(t_i , t_f ) u(t_i , t)_s ( x)u(t , t_i ) . a32introducing the complex mills time contours to connect @xmath117 to @xmath42 , @xmath42 to @xmath118 and @xmath118 to @xmath117 we form a ` closed - time ' contour @xmath445 ( the end - points of the contours @xmath446 together )
this allows to write the last equality ( ) in the compact form : @xmath447 where @xmath448 is the time - ordering on the contour @xmath445 operator .
the generating functional @xmath449 of correlation ( green ) functions has the form : @xmath450 where @xmath451 means averaging over the initial state .
if the initial correlations have a little effect , we can perform averaging over the gibbs ensemble .
this is the main assumption of the formalism : the generating functional of the green functions @xmath449 has the form in this case : @xmath452 with @xmath453 . in accordance with ( ) we have : @xmath454 and , as a result , z(j)=d e^i_c _ d^4 x \ { l(x ) + j(x)(x ) } 33where the path integration is performed with kms periodic boundary condition : @xmath455 in ( ) the contour @xmath456 connects @xmath117 to @xmath118 , @xmath118 to @xmath117 and @xmath117 to @xmath457
. therefore it contains an imaginary - time matsubara part @xmath117 to @xmath457 . a more symmetrical formulation uses the following realization : @xmath117 to @xmath118 , @xmath118 to @xmath458 , @xmath458 to @xmath459 and @xmath460 to @xmath457 ( e.g. ) .
this case also contains the imaginary - time parts of the time contour .
therefore , eq . ( ) presents the analytical continuation of the matsubara generating functional to real times .
one can note that if this analytical continuation is possible for @xmath449 then representation ( ) gives good recipe of regularization of frequency integrals in the matsubara perturbation theory , e.g. .
but it gives nothing new for our problem since the matsubara formalism is a formalism for equilibrium states only .
the starting point of our calculations is the @xmath4- into @xmath52-particles transition amplitude @xmath125 , the derivation of which is well known procedure in the lehmann - symanzik - zimmermann ( lsz ) reduction formalism framework , see also .
the @xmath461-point green function @xmath462 are introduced for this purpose through the generating functional @xmath463 :    g_nm(x , y)=(-i)^n+m_k=1^n(x_k)_k=1^m ( y_k)z_j , 22.5where ( x)= ,  * and z_j = de^is_j ( ) . 22.6 the action s_j ( ) = s_0()-v()+dxj(x)(x ) , 22.7 where @xmath464 is the free part and @xmath465 describes the interactions . at the end of the calculations one can put @xmath466 .    to provide the convergence of the integral ( [ 22.6 ] ) over the scalar field @xmath467
the action @xmath468 must contain a positive imaginary part .
usually for this purpose feynman s @xmath469-prescription is used .
but it is better for us to use the integral on the mills complex time contour @xmath470 .
for example , c _ :  tt+i , + 0 ,   -t+o5and after all the calculations return the time contour on the real axis putting @xmath471 .    in eq .
( [ 22.6 ] ) the integration is performed over all field configurations with standard vacuum boundary condition : d^4 x _ ( ^)= _ _ d_^=0 , 22.8 which leads to zero contribution from the surface term .
let us introduce now field @xmath472 through the equation : - = j(x ) 22.9 and perform the shift @xmath473 in integral ( [ 22.6 ] ) , conserving boundary condition ( [ 22.8 ] ) .
considering @xmath472 as the probe field created by the source : ( x)=dy g_0 ( x - y)j(y ) , + ( ^2 + m^2)_x g_0 ( x - y)=(x - y ) , 22.10 only the connected green function @xmath474 will be interesting for us .
therefore , g_nm^c(x , y)=(-i)^n+m_k=1^n(x_k ) _ k=1^m(y_k)z ( ) , 22.11 where z()=de^is()-iv(+ ) 22.12 is the new generating functional .    to calculate the nontrivial elements of the @xmath0-matrix we must put the external particles on the mass shell .
formally this procedure means amputation of the external legs of @xmath474 and further multiplication on the free particle wave functions . as a result
the amplitude of @xmath4- into @xmath52-particles transition @xmath125 in the momentum representation has the form : a_nm(q , p)=(-i)^n+m_k=1^n(q_k ) _ k=1^m^ * ( p_k ) z ( ) .
22.13 here we introduce the particle distraction operator ( q)=dxe^-iqx(x ) ,  ( x)=. 22.14    supposing that the momentum of particles are insufficient for us the probability of @xmath4- into @xmath52-particles transition is defined by the integral : r_nm = d_n ( q ) d_m ( p ) ^(4)(_k=1^nq_k - _ k=1^mp_k ) |a_nm|^2 , 22.15 where d_n(q)=_k=1^nd(q_k)= _ k=1^n , 22.16is the lorentz - invariant phase space element .
we assume that the energy - momentum conservation @xmath68-function was extracted from the amplitude .
note that @xmath475 is the divergent quantity . to avoid this problem with trivial divergence ,
connected integration over reference frame , let us divide the energy - momentum fixing @xmath68-function into two parts : ^(4)(q_k -p_k)=d^4 p ^(4)(p - q_k)^(4)(p - p_k ) 22.17 and consider a new quantity : r(p)=_n , md_n ( q ) d_m ( p ) ^(4)(p-_k=1^n q_k ) ^(4)(p-_k=1^n p_k ) |a_nm|^2 22.18defined on the energy momentum shell ( ) . here we suppose that the numbers of particles are not fixed .
it is not too hard to see that , up to a phase space volume , r = d^4p  r(p ) 22.19 is the imaginary part of the amplitude @xmath476 .
therefore , computing @xmath477 the standard renormalization procedure can be applied and the new divergences will not arise in our formalism .
the fourier transformation of @xmath68-functions in ( [ 22.18 ] ) allows one to write @xmath478 in the form : r(p)= e^ip(_1+_2 ) ( ) , 22.20where ( ) = _ n , m _ k=1^m\{d(q_k)e^-i_1 q_k }
_ k=1^n\{d(p_k)e^-i_2 p_k } |a_nm|^2 . 22.21
introducing the fourier - transformed probability @xmath479 we assume that the phase - space volume is not fixed exactly , i.e. it is proposed that the 4-vector @xmath235 is fixed with some accuracy if @xmath480 are fixed .
the energy and momentum in our approach are still locally conserved quantities since an amplitude @xmath125 is translationally invariant .
so , we can perform the transformation : _
1 q_k = ( _ 1 -_1 ) q_k + _ 1 q_k ( _ 1 -_1 ) q_k + _ 1 p  since 4-momenta are conserved . the choice of @xmath481 fixes the reference frame .
this degree of freedom of the theory was considered in .
inserting ( [ 22.13 ] ) into ( [ 22.21 ] ) we find that & ( ) = \{-idx dx ( _ + ( x)d_+-(x - x,_2)_-(x)- + & -_-(x)d_-+(x - x,_1)_+(x ) ) } z(_+)z^ * (-_-) , 22.22where @xmath482 and @xmath483 are the positive and negative frequency correlation functions correspondingly :
d_+-(x - x,)=-id_1(q)e^iq(x - x- ) 22.23 describes the process of particles creation at the time @xmath484 and its absorption at @xmath485 , @xmath486 , and @xmath57 is the cm 4-coordinate .
function d_-+(x - x,)=id_1(q)e^-iq(x - x+ ) 22.24 describes the opposite process , @xmath487 .
these functions obey the homogeneous equations : ( ^2 + m^2)_x g_+-=(^2 + m^2)_x g_-+=0 22.25since the propagation of on mass - shell particles is described .
we suppose that @xmath103 may be computed perturbatively . for this purpose
the following transformations will be used ( @xmath488 at @xmath489 ) : & e^-iv()= e^-idx ( x)(x ) e^idx j(x)(x ) e^-iv()= + & = e^dx ( x)(x ) e^-iv()= + & = e^-iv(-i ) e^idx j(x)(x ) , 22.26 where @xmath490 was defined in ( [ * ] ) and @xmath491 in ( [ 22.14 ] ) . at the end of the calculations , the auxiliary variables
@xmath492 , @xmath493 can be taken equal to zero . using the first equality in ( [ 22.26 ] ) we find that z()= e^-idx ( x)(x ) e^-iv(+ ) e^-dx dx j(x)d_++(x - x)j(x ) , 22.27where @xmath494 is the causal green function : ( ^2 + m^2)_x g_++ ( x - y)=(x - y ) .
22.28 inserting ( [ 22.27 ] ) into ( [ 22.22 ] ) after simple manipulations with differential operators , see ( [ 22.26 ] ) , we may find the expression : & ( ) = e^-iv(-i_+)+iv(-i_-)\ { dx dx + & ( j_+ ( x)d_+-(x - x,_1)j_- ( x)- j_- ( x)d_-+(x - x,_2)j_+ ( x)- + & - j_+ ( x)d_++(x - x)j_+ ( x)+ j_- ( x)d_(x - x)j_- ( x ) ) } , 22.29where d_=(d_++)^ * 22.30 is the anticausal green function .    considering the system with a large number of particles
, we can simplify the calculations choosing the cm frame @xmath495 .
it is useful also to rotate the contours of integration over @xmath496 for the result , omitting the unnecessary constant , we will consider @xmath497 .
external particles play a double role in the @xmath0-matrix approach : their interactions create and annihilate the system under consideration and , on the other hand , they are probes through which the measurement of a system is performed . since @xmath498 are the conjugate to the particles energies quantities we will interpret them as the inverse temperatures in the initial ( @xmath499 ) and final ( @xmath500 ) states of interacting fields .
they are the ` good ' parameters if and only if the energy correlations are relaxed .
0.2 cm    * kubo - martin - schwinger boundary condition *    0.2 cm    the simplest ( minimal ) choice of @xmath501 assumes that the system under consideration is surrounded by black - body radiation .
this interpretation restores niemi - semenoff s formulation of the real - time finite temperature field theory .    indeed , as follows from ( ) , the generating functional @xmath479 is defined by corresponding generating functional & _ 0 ( _ ) = z(_+)z^ * (-_-)
= d_+ d_- e^is_0(_+)-is_0(_- ) + &
e^-iv(_+ + _ + ) + iv(__- ) , 22.35 see ( [ 22.22 ] ) .
the fields @xmath502 and @xmath503 were defined on the time contours @xmath470 and @xmath504 .
as was mentioned above , see ( ) , the path integral ( [ 22.35 ] ) describes the closed path motion in the space of fields @xmath467 .
we want to use this fact and introduce a more general boundary condition which also guarantees the cancelation of the surface terms in the perturbation framework .
we will introduce the equality : _
_ d__+^_+ = _ _ d__- ^_- . 22.36
the solution of eq.([22.36 ] ) requires that the fields @xmath505 and @xmath506 ( and their first derivatives @xmath507 ) coincide on the boundary hypersurface @xmath65 : _(_)=(_)0 , 22.37 where , by definition , @xmath508 is the arbitrary  turning - point " field .    in the absence of the surface terms
, the existence of a nontrivial field @xmath509 has the influence only on the structure of the green functions g_++=<t_+_+ > ,  g_+-=<_+_- > , + g_-+=<_-_+ > ,  g_=<_-_- > , 22.38 where @xmath510 is the antitemporal time ordering operator .
these green functions must obey the equations : & ( ^2 + m^2)_x g_+- ( x - y)= ( ^2 + m^2)_x g_-+ ( x - y)=0 , + & ( ^2 + m^2)_x g_++ ( x - y)= ( ^2 + m^2)_x^ * g_ ( x - y)=(x - y ) , 22.39 and the general solution of these equations : g_ii = d_ii+g_ii , + g_ij = g_ij ,   ij 22.40 contain the arbitrary terms @xmath511 which are the solutions of homogenous equations : ( ^2 + m^2)_x g_ij ( x - y)=0 ,  i , j=+,- .
22.41 the general solution of these equations ( they are distinguished by the choice of the time contours @xmath110 ) g_ij(x - x)=d_1(q ) e^iq(x - x ) n_ij ( q ) 22.42 are defined through the functions @xmath512 which are the functionals of the ` turning - point ' field @xmath508 : if @xmath513 we must have @xmath514 .
our aim is to define @xmath512 .
we can suppose that @xmath515 the simplest supposition gives : n_ij~<_i_j>~<^2(_)>. 22.43 we will find the exact definition of @xmath512 starting from the @xmath0-matrix interpretation of the theory .
it was noted previously that the turning - point field @xmath509 may be arbitrary .
we will suppose that on the remote @xmath65 there are only free , on the mass - shell , particles .
formally it follows from ( [ 22.40 ] - [ 22.42 ] ) .
this assumption is natural also in the @xmath0-matrix framework . in other respects
the choice of boundary condition is arbitrary .
therefore , we wish to describe the evolution of the system in a background field of mass - shell particles . the restrictions connected with energy - momentum conservation laws will be taken into account and in other respects background particles are free . then our derivation is the same as in . here
we restrict ourselves mentioning only the main quantitative points .    calculating the product @xmath516
we describe a time ordered processes of particle creation and absorption described by @xmath482 and @xmath483 . in presence of the background particles ,
this time - ordered picture is slurred over because of the possibility to absorb particles before their creation appears .
the processes of creation and absorption are described in vacuum by the product of operators @xmath517 and @xmath518 .
we can derive ( see also ) the generalizations of ( [ 22.22 ] ) .
the presence of the background particles will lead to the following generating functional : r_cp = e^-i*n*(_i^*_j)r_0 ( _ ) , 22.44 where @xmath519 is the generating functional for the vacuum case , see ( [ 22.35 ] ) . the operator @xmath520 describes the external particles environment .
the operator @xmath521 can be considered as the creation and @xmath522 as the annihilation operator and the product @xmath523 acts as the activity operator .
so , in the expansion of @xmath524 we can leave only the first nontrivial term : ( _ i^*_j)= d(q ) ^*_i ( q ) n_ij _ j ( q ) , 22.45since no special correlation among background particles should be expected .
if the external ( nondynamical ) correlations are present then the higher powers of @xmath525 will appear in expansion ( [ 22.45 ] ) . following the interpretation of @xmath525
, we conclude that @xmath512 is the mean multiplicity of background particles .
computing @xmath526 we must conserve the translation invariance of amplitudes in the background field . then , to take into account the energy - momentum conservation laws one should adjust to each vertex of in - going @xmath125 particles the factor @xmath527 and for each out - going particle we have correspondingly @xmath528 .
so , the product @xmath529 can be interpreted as the probability factor of the one - particle ( creation+annihilation ) process .
the @xmath4-particles @xmath530 process probability is the simple product of these factors if there is no special correlations among background particles .
this interpretation is evident in the cm frame @xmath531 .    after these preliminaries , it is not too hard to find that in the cm frame we have : & n_++(q_0)=n_(q_0)= = + & = = ( |q_0| )
. 22.47 computing @xmath512 for @xmath532 we must take into account the presence of one more particle : n_+-(q_0)= ( q_0 ) + ( -q_0 ) = + = ( q_0)(1+(q_0 _ 1))+ ( -q_0)(-q_0
_ 1 ) 22.48 and n_-+(q_0)= ( q_0)(q_0 _ 2)+ ( -q_0)(1 + ( -q_0 _ 2 ) ) .
22.49 using ( [ 22.47 ] ) , ( [ 22.48 ] ) and ( [ 22.49 ] ) , and the definition ( [ 22.40 ] ) we find the green functions : g_i , j(x - x,())= e^iq(x - x ) _ ij ( q , ( ) ) 22.50 where i_ij ( q , ( ) ) = ( ) + +   + + 2(q^2 -m^2 ) ( ) 22.51 and a_()=-e^(|q_0|q_0 ) .
22.52 the corresponding generating functional has the standard form : r_p(j_)=e^-iv(-i_+)+iv(-i_- ) e^dx dx j_i ( x)g_ij(x - x,())j_j(x ) 22.53where the summation over repeated indexes is assumed .    inserting ( [ 22.53 ] ) in the equation of state ( [ 3 ] ) we can find that @xmath533 . if @xmath534 is a ` good ' parameter then @xmath535 coincides with the green functions of the real - time finite - temperature field theory and the kms boundary condition : g_+-(t - t)=g_-+(t - t-i ) , g_-+(t - t)=g_+-(t - t+i ) , 22.54 is restored .
the eq.([22.54 ] ) can be deduced from ( [ 22.51 ] ) by direct calculations .
we start this consideration from the assumption that the temperature fluctuations are large scale .
we can assume that the temperature is a ` good ' parameter in a cell the dimension of which is much smaller then the fluctuation scale of temperature .
( the ` good ' parameter means that the corresponding fluctuations are gaussian . )
let us divide the remote hypersurface @xmath536 on a @xmath537 and let us propose that we can measure the energy and momentum of groups of in- and out - going particles in each cell .
the 4-dimension of cells can not be arbitrary small because of the quantum uncertainty principle .
to describe this situation we decompose the @xmath68-function of the initial state constraint ( ) on the product of @xmath538 @xmath68-functions : @xmath539 where @xmath540 are the momentum of @xmath84-th in - going particle in the @xmath239-th cell and @xmath541 is the total 4-momenta of @xmath542 in - going particles in this cell , @xmath543 .
therefore , @xmath544    the same decomposition will be used for the second @xmath68-function of outgoing particle constraints .
we must take into account the multinomial character of particle decomposition on @xmath545 groups
. this will give the coefficient : @xmath546 where @xmath547 is the kronecker s symbol .
the summation over @xmath548 is assumed .    in result ,
the quantity & r_n_c(p , q)= _ \{n , m}_n_c |a_nm|^2 + & ^n_c_=1\ { ^(4)(q_-^m__k=1q_k , ) ^(4)(p_-^n__k=1p_k , ) } 34defines the probability to find in the @xmath239-th cell the fluxes of in - going particles with total 4-momentum @xmath541 and of out - going particles with the total 4-momentum @xmath549 . the sequence of these two measurements is not fixed .
the fourier transformation of @xmath68-functions in ( [ 34 ] ) gives : @xmath550 where @xmath551 has the form : _ n_c()=_=1^n_c\{^m__k=1 e^-i_1,q_k , ^n__k=1 e^-i_2,p_k , } |a_nm|^2 .
35inserting @xmath552 into ( [ 35 ] ) we find : & _ n_c()= \ { i_=1^n_c dx dx [ _ + ( x)d_+-(x - x;_2,)_-(x)- .
+ & .-_-(x)d_-+(x - x;_1,)_+(x ) ] } _ 0 ( ) , 36 where @xmath553 , and @xmath554 are the positive and negative frequency correlation functions
.    we must integrate over sets @xmath555 and @xmath556 if the distribution of momenta over cells is not fixed . in the result ,
r(p)=d^4_1(p ) d^4_2 ( p)_n_c ( ) , 37where the differential measure @xmath557 takes into account the energy - momentum conservation laws : @xmath558 the explicit integration gives that @xmath559 where @xmath560 is 3-vector of the cm frame .
choosing cm frame , @xmath561 , @xmath562    in this frame @xmath563 where @xmath564 and @xmath565 was defined in ( [ 36 ] ) with @xmath566 .
we will calculate integrals over @xmath498 using the stationary phase method .
the equations for the most probable values of @xmath498 : -k(e,\{}_n_c)= _ n_c ( ) ,  k=1,2 , 38always have unique positive solutions @xmath567 .
we propose that the fluctuations of @xmath568 near @xmath569 are small , i.e. are gaussian .
this is the basis of the local - equilibrium hypothesis . in this case
@xmath570 is the temperature in the initial state in the measurement cell @xmath239 and @xmath571 is the temperature of the final state in the @xmath239-th measurement cell .
the last formulation ( [ 37 ] ) implies that the 4-momenta @xmath555 and @xmath556 can not be measured .
it is possible to consider another formulation also .
for instance , we can suppose that the initial set @xmath555 is fixed ( measured ) but @xmath556 is not .
in this case we will have a mixed experiment : @xmath572 is defined by the equation : @xmath573 and @xmath574 is defined by the second equation in ( [ 38 ] ) .
considering the continuum limit , @xmath575 , the dimension of the cells tends to zero . in this case
we are forced by quantum uncertainty principle to assume that the 4-momenta sets @xmath576 and @xmath577 are not fixed .
this formulation becomes pure thermodynamical : we must assume that just @xmath578 and @xmath579 are measurable quantities .
for instance , we can fix @xmath578 and try to find @xmath579 as a function of the total energy @xmath10 and the functional of @xmath578 . in this case eqs.([38 ] ) become the functional equations .    in the considered microcanonical description , the finiteness of temperature does not touch the quantization mechanism .
indeed , one can see from ( [ 36 ] ) that all thermodynamical information is confined in the operator exponent @xmath580 the expansion of which describes the environment , and the ` mechanical ' perturbations are described by the functional @xmath581 .
this factorization was achieved by the introduction of the auxiliary field @xmath472 and is independent from the choice of boundary conditions , i.e. from the choice of the systems environment .
0.2 cm    * wigner functions *    0.2 cm    we will adopt the wigner functions formalism in the carruthers - zachariazen formulation . for the sake of generality , the @xmath52 into @xmath4 particles transition will be considered .
this will allow the inclusion of the heavy ion - ion collisions .    in the previous section ,
the generating functional @xmath565 was calculated by means of dividing the ` measuring device ' on the remote hypersurface @xmath536 into @xmath537 cells _ n_c()=e^-i*n*(;,z)_0 ( ) , 39where & * n*(;,z)= \{_=1^n_c dx dx. + & .(_+(x)d_+-(x - x;_2,,z_2)_-(x).- + & .-_-(x)d_-+(x - x;_1,,z_1)_+(x ) ) } 40is the particle number operator .
the frequency correlation functions @xmath482 and @xmath482 are defined by equalities : d_+-(x - x;_2,,z_2)=-id_1(q)e^iq_2,(x - x ) e^-_2,(q_2 , ) z_2(q_2 , ) a3.5 + d_-+(x - x,_1,,z_1)=id_1(q)e^-iq_1,(x - x ) e^-_1,(q_1,)z_1(q_1 , ) b3.5it was assumed that the dimension of the device cells tends to zero ( @xmath575 ) .
now we wish to specify the cells coordinates . in the result
we will get to the wigner function formalism .
let us introduce wigner variables : x - x=r ,  x+x=2y :  x = y+r/2 ,  x=y - r/2 .
3.6then & * n*(;,z)=-i_=1^n_c d(q)dr + & ( _ + ( y+r/2)_-(y - r/2)z_2(q_2,)e^iq_2,r e^-_2,(q_2,).+ + & + ._-(y+r/2)_+(y - r/2)z_1(q_1 , ) e^-iq_1,re^-_1,(q_1,))dy 3.10the boltzmann factor , @xmath582 , can be interpreted as the probability to find a particle with the energy @xmath583 in the final ( @xmath584 ) or initial ( @xmath585 ) state .
the total probability , i.e. the process of creation and further absorption of @xmath4 particles , is defined by multiplication of these factors . besides @xmath586 is the out - going particle momentum measured in the @xmath239-th cell .
generally it is impossible to adjust the 4-index of cell @xmath239 with coordinate @xmath587 .
for this reason the summation over @xmath239 and the integration over @xmath588 are performed in ( ) independently . but let us assume that the 4-dimension of the cell @xmath589 is higher then the scale of the characteristic quantum fluctuations @xmath590 , l>>l_q .
o6one can divide the 4-dimensional @xmath587 space into the @xmath589 dimensional cells . then , because of ( ) , the quantum fluctuations can not take away particles from this cell .
then we can adjust the index of the measurements cell with the index of the @xmath587 space cell .    as a result , & * n*(;,z)=-idyd_i(q)dr + & ( _ + ( y+r/2)_-(y - r/2)z_2(q_2,y)e^iq_2r e^-_2(y)(q_2)+ .
+ & .+_-(y+r/2)_+(y - r/2)z_1(q_1,y ) ) e^-iq_1re^-_1(y)(q_1 ) , 3.10where dy=__c()dy , o7and @xmath591 is the dimension @xmath589 of the @xmath587 space cell with index @xmath239 .
notice that the momentum @xmath592 did not carry the index @xmath239 ( or the index @xmath587 of the space cell ) .
our formalism allows the introduction of more general ` closed - path ' boundary conditions .
the presence of external black - body radiation will only reorganize the differential operator @xmath593 and a new generating functional @xmath526 has the same form :    @xmath594 the calculation of operator @xmath595 is strictly the same as in appendix b. introducing the cells we will find that @xmath596 where the occupation number @xmath597 carries the cell index @xmath587 : @xmath598 and ( @xmath599 ) @xmath600 @xmath601 @xmath602 for simplicity the cm system was used .
other calculations are the same as the constant temperature case .
first of all , two light - like 4-momenta @xmath603 are introduced . here
@xmath604 are momenta of colliding particles
. the final state particles momenta have the following representation : & p_1=_1p_2+_1p_1+p_1 ,   p_2=_2p_2+_2p_1+p_2 ,   + & k_i=_ip_2+_ip_1+k_i .
m1sudakov s parameters , @xmath605 , are not independent .
the mass shell conditions and the energy - momentum conservation laws give : & s_1_1=m^2+(p_1)^2=e_1 ^ 2 ,   s_2_2=e_1 ^ 2 ,  s_i_i = e_i^2 , + & _1+_2+_i=1 ,   _1+_2+_i=1 , m2where @xmath606 is the transverse energy .
we have for the multiperipheral kinematics : 1_1>>_1>> ... >>_n>>_2~ , + < < _1<<_1<< ... <<_n<<_2~1 m3and the transverse momenta are restricted :    considered cm frame : _
i= ,  |_i|>>|_i| , m5if the particle moves along @xmath607 , and a similar expression exists for particles moving in the opposite direction , where @xmath608 . in the ` central region ' of the cm frame @xmath609
the angles of produced particles are large and energies are small .
it should be underlined that all this excludes the ( mini)jets formation .
the final - state particles phase space volume element is & d_22+n = c^n_v + & ( _ 2-_1 ) ^n+1_i=1()^2(q_i^2)= dz_n , where @xmath610 and the 4-momentum of produced particle @xmath611 the square of pairs invariant mass : @xmath612 the energy conservation law takes the form : @xmath613 the trajectory of reggeized gluon is @xmath614 where @xmath615 is the gluon ` mass ' . if this virtuality is large , @xmath616 , then the gluon decay creating a pqcd jet .
but the constraint on the multiperipheral kinematics prevent this possibility .
0.2 cm    * deep inelastic reactions *    0.2 cm    for the pure deep inelastic case , when one of the initial hadrons is scattered at the angle @xmath617 have the energy @xmath618 in the cms of beams whereas the another is scattered at small angle and the large transfer momentum @xmath619 is distributed to the some number of the emitted particles due to evolution mechanism we have ( @xmath617 is small ) : & d_n^dis= dd_n ded , + & dd_n=()^n_m^2^q^2 _
m^2^k_n^2_m^2^k_2_x^1 d_n^(1)_n__n^1 d_n-1^(1)_n-1 ... + & _ _ 2 ^ 1d _ 1^(1)_1 p() ... p(_1 ) , p(z)=2 , where @xmath620 and the emission angle @xmath621 .
0.2 cm    * large angle production *    0.2 cm    for the large - angle particle production process the differential cross section ( as well as the total one ) fall down with cm energy @xmath208 . let us consider for definiteness the process of annihilation of electron - positron pairs to @xmath4 photons : & d_n^dl = df_n , + & df_n=()^n ^n_i=1d x_idy_i(x_i - y_i)(y_i - y_i-1 ) ( y_i)(x_i ) + & ( -x_n)(-y_n),y_i=,= , the similar formulae takes place for subprocess of quark - antiquark annihilation into @xmath4 large - angle moving gluons .    at the end , one can consider the following possibilities : + a ) pomeron regime ( p ) ; + b ) evolution regime ( dis ) ; + c ) double logarithmic regime ( dl ) ; + d ) dis+p regime ; + e ) p+dl+p regime .
+ the description of every regime may be performed in terms of effective ladder - type feynman diagrams .
this can be done using the blocks @xmath622 , @xmath623 and @xmath624 .
we will consider , see ( ) ( q , ;z)=_0^de^-p(q,;z)= ,  = ( s / m^2 ) .
s9as the ` propagator of the cut pomeron ' .
it will be assumed also that @xmath625 so , the resonance short range correlations will be ignored in this definition or propagator .
it was assumed also that the ` bare ' slope @xmath229 is @xmath129 independent .
it should be underlined that the ` propagator ' ( ) is written phenomenologically .
it absorbs the assumptions that ( i ) the diffraction cone shrinks with energy and ( ii ) the inclusive cross sections are universal , see ( )    the set of principal rules concerning multiperipheral kinematics of feynman diagrams is given in appendix d. the reggeon calculus supposes that the virtuality of each line of the feynman diagram is restricted .
this ignores ` hard jets ' , later known as the pqcd jets .
then the @xmath239 pomeron exchange eikonal diagram has only @xmath626 ways of being cut .
if the cut line goes through @xmath627 pomerons , then the corresponding contributions are : ^_( , q)= d_(m^_(q_1, ... ,q_))^2 y^__l=1^(q_i , _i;z)_i=+1^(q_i , _i;z=1 ) , s1where @xmath628 is the ` vertex function ' , the combinatorial coefficient is @xmath629 and the phase space element is d_=_i=1^ ( -_l=1^_l)^2(q-_i=1^q_l ) .
s2    as usual , the contribution ( ) leads to the following mean multiplicity of produced particles : ^_=.d_(s / m^2)^^_( , q=0)|_z=1 ~. s3therefore , ~s4is essential in the vhm region , where @xmath630 is assumed .    the impact parameter representation : ^_(s , q)=d()^ e^i^_( , q ) s5would be useful also .
the contribution ( ) describes interactions with impact parameter 4(s / m^2)/. s6notice @xmath631 is the number of cut pomerons independent of @xmath627 .
but , remembering that @xmath632 and that the regge model is only able to describe large distance interactions , @xmath633 , one can conclude that the regge pole description is valid only for n^2 .
s7thus is why the vhm region is defined by @xmath634 .
then the cut pomeron propagator in the impact parameter representation ( , ; z)=g(,)e^(z-1 ) , l3where g(,)=e^e^-^2/4 l4is the uncut pomeron profile function . using this definition
one can find that the contribution of the eikonal diagrams gives a contribution : _ 0
( , ; z)= ( 1-e^-^2g(,))- ( 1-e^2 ^ 2g(,)(e^(z-1)-1 ) ) , l5where @xmath615 is a constant .
first bracket is essential for @xmath636 .
so , with exponential accuracy , the first term is equal to @xmath637 let us now consider the second bracket . for @xmath638 ^24^2\{1+(1-e^(z-1))}= 4^2(,z ) l6are essential .
it is not hard to see that @xmath639 decreases if @xmath640 and @xmath639 is equal to zero for @xmath106 . in this case , by the definition of the generating function , the integral over @xmath230 of @xmath641 defines the contribution to the total cross section .
so , the model predicts the production of particles in the ring : 4^2(,z)^24^2 .
l7if @xmath638 , i.e. if @xmath642 .
notice also that @xmath643
. then , _ 0 ( , ; z=0)=_0 ( , ; z=1 ) .
l8so , the elastic part of the total cross section is half of the total cross section . this means , using the optical analogy , that the scattering on the absolutely black disk is well described .
the last conclusion means that the interaction radii should increase with @xmath4 in the vhm region . indeed , as follows from ( ) , _ 0 ( , ; z)(e^2 ^ 2g(,)(e^(z-1)-1)-1 ) , l9at @xmath132 and , therefore , 0 ^ 2b^2=4(+(z-1 ) ) l10are essential .
our purpose is to investigate the role of the exponential spectrum ( ) in the asymptotic region over multiplicity @xmath4 . in this case one can validate heavy resonance creation and such formulation of the problem have definite advantages .
\(i ) if creation of heavy resonances at @xmath644 is expected , then one can neglect the dependence on the resonance momentum @xmath645 .
so , the ` low - temperature ' expansion is valid in the vhm region .
\(ii ) having the big parameter @xmath4 , one can construct the perturbations expanding over @xmath142 .
\(iii ) we will be able to show at the end the range of applicability of these assumptions .
+ for this purpose , the following formal phenomena will be used .
the grand partition function t(z , s)=_n z^n _ n(s ) ,  t(1,s)=_tot(s ) ,   n / m_0n_max(s ) , gpfwill be introduced , see ( ) . then the inverse mellin transformation , see ( ) _ n(s)=t(z , s ) .
metrwill be performed expanding it in the vicinity of the solution @xmath646 of the equation of state , see ( ): n = zt(z , s ) .
eqst1it is assumed , and this should be confirmed at the end , that the fluctuations in the vicinity of @xmath130 are gaussian .
it is natural at first glance to consider @xmath647 as an increasing function of @xmath4 . indeed
, this immediately follows from the positivity of @xmath19 and the finiteness of @xmath648 at finite @xmath7 .
but one can consider the ` thermodynamical limit ' , see sec.2.3.1 , or the limit @xmath649 .
theoretically , the last one is right because of the pcac hypotheses and nothing should happen if the pion mass @xmath649 . in this sense @xmath137
may be considered as the whole function of @xmath129 .
then , @xmath647 would be an increasing function of @xmath4 if and only if @xmath137 is a regular function at @xmath106 .
the proof of this statement is as follows .
we should conclude , as follows from eq .
( ) , that z_c(n , s)z_s  at  n ,  and  at  s = const , asyi.e .
the singularity point @xmath138 @xmath650 @xmath130 in asymptotics over @xmath4 .
if @xmath651 , then @xmath652 , when @xmath4 tends to infinity .
the concrete realization of this possibility is shown in sec.3.3 .
but if @xmath653 , then @xmath654 in vhm region , see secs.3.1 , 3.5 .    on may use the estimation , see also ( ):
-= z_c(n , s)+o(1/n ) , est0where @xmath130 is the @xmath655 solution of ( ) .
it should be underlined that this estimation is @xmath656 of the character of singularity , i.e. the position @xmath138 is only important with @xmath657 accuracy .
0.2 cm    * partition function *    0.2 cm    introducing the ` grand partition function ' ( ) the ` two - level ' description means that & = _ k^k_i=1\{d_k(q)dm_i(q_i , z)e^-_i } + & n_k(q_1,q_2, ... ,q_k;)-(z , s ) , vdwhere @xmath658 .
this is our group decomposition .
the quantity @xmath659 may be considered as the local activity .
so , .|_=1~_totb_1(q ) r1cfif the resonance decay forms a group of particles with total 4-momentum @xmath592 , then @xmath660 is the mean number of such groups .
the second derivative gives : .|_=1~_tot \{b_2(q_1,q_2)-b_1(q_1)b_1(q_2)}_totk_2(q_1,q_2 ) r2cfwhere @xmath661 is the two groups correlation function , and so on .
one can consider @xmath662 as mayer s group coefficients , see sec.2.3.3 .
the lagrange multiplier @xmath62 was introduced in ( ) to each resonance : the boltzmann exponent @xmath663 takes into account the energy conservation law @xmath664 where @xmath10 is the total energy of colliding particles , @xmath665 in the cm frame .
this conservation law means that @xmath62 is defined by equation : = t(z , ) .
eqst2so , to define the state one should solve two equations of state ( ) and ( ) .    the solution @xmath81 of the eq .
( ) has the meaning of inverse temperature of the gas of resonances if and only if the fluctuations in vicinity of @xmath81 are gaussian , see sec.2.2.2 .    on the second level , we should describe the resonance decay into hadrons . using ( )
we can write in the vicinity of @xmath106 : ( q , z)=_n z^n ^r_n(q)=g^r()e^(z-1)|n(m ) . rgpfthe assumptions * b * and * d * , see ( ) , were used here
so , -(z , s)=_k^k_i=1\{dm_i(m_i , z)}_k(m ; ) , vd1where @xmath666 @xmath667 was defined in ( ) and _ k(m;)=_i=1^k\ { d_k(q)e^-_i(q_i ) } b_k(m;q ) .
1assuming now that @xmath668 are essential , _
k(m;)b_k(m ) _ i=1^k\{e^-m_i } est    following the duality assumption , one may write : b_k(m)=|b_k(m ) _ i=1^k
\{m_i^e^_0 m_i } stringand @xmath669 is a slowly varying function of @xmath666 : @xmath670 in the result , the low - temperature expansion looks as follows : -(z , s)=_k \{_m_0^dm m^+3/2 e^(z-1)|n^r(m)-(-_0)m}^k .
fe1we should assume that @xmath671 .
in this sense one may consider @xmath672 as the limiting temperature and the above mentioned constraint means that the resonance energies should be high enough .
0.2 cm    * thermodynamical parameters *    0.2 cm    remembering that the position of the singularity over @xmath129 is essential , let us assume that the resonance interactions can not renormalize it , i.e. that the sum ( ) is convergent . then , leaving the first term in the sum ( ) , -(z , s)= _ m_0^dm ( m / m_0)^+3/2 e^(z-1)|n^r(m)-(-_0)m .
fe2we expect that this assumption is satisfied if & _
m_0^dm m^+3/2 e^(z-1)|n^r(m)-(-_0)m > > + & \{_m_0^dm m^+3/2 e^(z-1)|n^r(m)-(-_0)m}^2 res1for n ,  s ,   < < 1 restr
so , we would solve our equations of state with the following free energy : -(z , s)= _ m_0^d ( ) ( ) ^-1 e^-(m / m_0 ) , fe3where , using ( ) , =+2(z-1)|n^r_0 + 5/2=2(z-1)|n^r_0 ,   = m_0(-_0)0 ,  = const .
nw    we have in terms of these new variables the following equation for @xmath129 , n = z .
eq1cthe equation for @xmath62 takes the form : n_max= , eq2cwhere @xmath673 and @xmath674 is the incomplete @xmath675-function : @xmath676    0.2 cm    * asymptotic solutions *    0.2 cm    following physical intuition , one should expect the cooling of the system when @xmath644 , for fixed @xmath208 , and heating when @xmath677 , for fixed @xmath4 .
but , as was mentioned above , since the solution of eq .
( ) @xmath81 is defined by the value of the total energy , one should expect that @xmath81 decreases in both cases .
so , the solution _
c0 ,   < 0  at  n ,   < 0  at  ssol2is natural for our consideration .
the physical meaning of @xmath129 is activity .
it defines at @xmath678 the work needed for one particle creation .
then , if the system is stable and @xmath137 may be singular at @xmath132 only , > 0  at  n ,   < 0  at  s. sol3    one should assume solving equations ( ) and ( ) that z_c^_c+1 < < ( _c+1,_c ) .
ineq1this condition contains the physical requirement that @xmath133 . in the opposite case
, the finiteness of the phase space for @xmath679 should be taken into account .    as was mentioned above
, the singularity @xmath138 attracts @xmath130 at @xmath644 .
for this reason one may consider the following solutions .
a. @xmath680 : @xmath681 , @xmath682 .    in this case ^-(,)~e^(/ ) .
a1this estimation gives the following equations : n = c_1(/)e^(/ ) , = c_2 ()<<1 , a2where @xmath683 are the unimportant constants .
the inequality is a consequence of ( ) .
these equations have the following solutions : _ c<<1 ,   _c~n>>1 . a3using them
one can see from ( ) that it gives _
n < o(e^-n ) . a4
b. @xmath684 : @xmath685 , @xmath686 .
one should estimate @xmath687 near the singularity at @xmath106 and in the vicinity of @xmath246 to consider the consequence of this solution .
expanding @xmath687 over @xmath248 at @xmath688 , ( ,)=()-^e^-+o(^+1 ) + o(1 ) .
b1this gives the following equations for @xmath689 : n = c_1e^(1/ ) .
b2the equation for @xmath248 has the form : n_max = c_2e^(+1)(1/ ) .
b3where @xmath690 are unimportant constants .    at 0<(1/)-1<<1 ,
i.e.  at  ( 1/)<<n<<^2(1/ ) , b4we find : _c~. b5inserting this solution into ( ): _ c~. b6it is remarkable that @xmath691 in the leading approximation is @xmath4 independent . by this reason @xmath692 becomes
@xmath4 independent also : _c~ :   z_c=1 + .
b7this means that _
n = o(e^-n ) b8and obeys the kno scaling with mean multiplicity @xmath693 .
considering particle creation in the dis processes , one should distinguish correlation of particles in the ( mini)jets and the correlations between ( mini)jets .
we will start from the description of the jet correlations .
one should introduce the inclusive cross section for the @xmath239 jets creation @xmath694 where @xmath695 , @xmath696 are the jets 4-momentum in the dis kinetics , @xmath697 .
having @xmath698 we can find the correlation functions @xmath699 where @xmath700 and @xmath701 defines the sort of created colour particle .
it is useful to introduce the generating functional f^ab(q^2,x;w)=_nd_n(k)_i=1^n w^r_i(k_i )    the amplitude , @xmath702 is the phase space volume and @xmath703 are the arbitrary functions .
it is evident that f^ab(q^2,x;w).|_w=1=d^ab(q^2,x ) .
11athe inclusive cross sections _^(r)(k_1,k_2, .... ,k_;q^2,x)=_i=1^f^ab(q^2,x;w).|_w=1 .
12the correlation function n_^(r)(k;q^2,x ) = _ 1=1^f^ab(q^2,x;w).|_w=1 .
13awe can find the partial structure functions @xmath704 , where @xmath4 is the number of produced ( time - like ) gluons , using there definitions .
it will be useful to introduce the laplace transform over the variable @xmath705 : f^ab(q^2,x;w)=_j<0()^j f^ab(q^2,j;w ) 14a    the expansion parameter of our problem @xmath706 .
for this reason one should take into account all possible cuts of the ladder diagrams .
so , calculating @xmath707 in the lla all possible cuts of the skeleton ladder diagrams are defined by the factor : \{^ab_rg_r_r^ab } , 15ai.e
. the cut line may not only get through the exact green function @xmath708 but through the exact vertex functions @xmath709 also ( @xmath710 are negative ) .
we have in the lla ( see appendix d ) @xmath711 and @xmath712 following our approximation , see previous section , we would not distinguish the way in which the cut line goes through the born amplitude @xmath713 we will simply associate @xmath714 to each rung of the ladder .    considering the asymptotics over @xmath4 , the time - like partons virtuality @xmath715 should be maximal . here
@xmath716 is the fraction of the longitudinal momentum of the jet .
then , slightly limiting the jets phase space , k_i^2=|q_i+1|^2(1+o((1/x)/|q_i+1|^2 ) ) .
16    as a result , introducing @xmath717 , where @xmath718 , @xmath719 in the lla variable , we can find the following set of equations : f_ab(q^2,j;w)=_c , r^r_ac(j)w^r ( ) f_ab(q^2,x;w ) , 17awhere
^r_ac(j)=_ac(j)=^1_0x^jp_ac(x ) 18aand @xmath720 is the regular kernel of the bethe - salpeter equation . at @xmath721
this equation is the ordinary one for @xmath707 .
we will search the correlation functions from eq .
( ) in terms of the laplace transform @xmath722 let us write : f_ab(q^2,j;w)=d_ab(q^2,j)\ { _ _ i=1^ ( ( w^r_i(_i)-1))n_ab^(r)_(k;q^2,j ) } 19ainserting ( ) into ( ) and expanding over @xmath723 we find the sequence of coupled equations .    omitting the cumbersome calculations , we write in the lla that & _ ab^(r)_(_1,_2, .... ,_;q^2,j)= + & = d_ac_1(j,_1)^r_1_c_1c_2(j)d_c_2c_3(j,_2)^r__c_c_+1(j)d_c_+1b(j,_+1 ) .
20one should take into account the conservation laws : _
1_2_+1= ,   _ 1<_2<  <_+1<.
21computing the laplace transform of this expression we find @xmath724    the kernel @xmath725 was introduced in ( ) .
let us write it in the form : d_ab(j,)=_=^_(j ) , 22where d^_qq=_s-_gg ,  d^_qg=_s-_qq ,   d^_qg=_gq ,  d^_gq=_qg 23and _
= \{_qq+_gg+^1/2}. 24if @xmath726 , then @xmath727 are essential . in this case , _ gg~_gq>>_qg~_qq = o(1 ) .
25this means the gluon jets dominance and n^g_gg=_gg+o(1 ) .
26one can find the following estimation of the two - jet correlation function : n_ab^r_1r_2(_1,_2;,)= o(max\{(_1/)^_gg,(_2/)^_gg , ( _ 1/_2)^_gg}}. 27this correlation function is small since in the lla @xmath728 .
this means that the jet correlation becomes high if and only if the mass of the correlated jets are comparable . but this condition shrinks the range of integration over @xmath729 and for this reason one may neglect the ` short - range ' correlations among jets . therefore , as follows from ( ) , f_ab(q^2,j;w)=d_gg(,j)\ { _
gg^__0 w^g( ) } 28we will use this expression to find the multiplicity distribution in the dis domain
.    0.2 cm    * generating function *    0.2 cm    to describe particle production , one should replace : @xmath730 where @xmath731 is the @xmath173 of @xmath4 particle production , _ nw^r_n=1 .
29having @xmath239 jets , one should take into account the conservation condition @xmath732 .
for this reason , the generating functions formalism is useful . in the result
, one can find that if we take ( ) w^g = w^g(,z ) ,  w^g(,z).|_z=1=1 , 30then @xmath733 defined by ( ) is the generating functional of the multiplicity distribution in the ` @xmath492 representation ' . in this expression @xmath734
is the generating function of the multiplicity distribution in the jet of mass @xmath735 .    in the result , see ( ) , f^ab(q^2,x;w)_j<0(1/x)^j e^_gg(,z ) 31where (,z)=^__0w^g(,z ) .
32noting the normalization condition ( ) , (,z=1)=. 33    the integral ( ) may be calculated by the steepest descent method .
it is not hard to see that jj_c=1+\{4n(,z)/(1/x)}^1/2 34is essential .
notice that @xmath736 should be essential but we find , instead of the constraint ( ) , that (,z)<<(1/x ) .
35 in the frame of this constraint , f^ab(q^2,x;w)\{4}. 36    generally speaking , there exist such values of @xmath129 that @xmath737 .
this is possible if @xmath738 is a regular function of @xmath129 at @xmath106 .
then @xmath130 should be an increasing function of @xmath4 and consequently @xmath739 would be an increasing function of @xmath4 .
therefore , one may expect that in the vhm domain @xmath737
.    then @xmath740 would be essential in the integral ( ) .
this leads to the following estimation : @xmath741 but this is impossible since @xmath742 should be an increasing function of @xmath129 .
this shows that the estimation ( ) has a finite range of validity .
solution of this problem with unitarity is evident .
one should take into account correlations among jets considering the expansion ( ) . indeed , smallness of @xmath743 may be compensated by large values of @xmath744 in the vhm domain .
one may neglect quark jets in the vhm region since the gluons mean multiplicity @xmath745 - quarks multiplicity and in the vhm region the leftmost singularities are important
. then we can write : t_j(,z)=t_j(,z ) ^__0d(t_j(,z)-1 ) , j1where @xmath746 and @xmath747 is the generating function of the distribution over the number of gluons @xmath748 : t_j(,z)=_n
z^n w_n ( ) ,  t_j(,z=1)=1 . j2we search a solution in the vhm region , where n>>n_j\ { } ,  a=12/11 .
j3let us consider the following solution : w_n=()^e^-n/|n_j .
j4it is useful to introduce _ k()=_n=1^n^k-1w_n ( ) j5for this solution . inserting ( ) in this expression , _
k()=|n^k()_k , j6where @xmath498 ( i ) should be positive and ( ii ) @xmath729 independent .
these conditions are satisfied for the following values of @xmath84 . indeed , at @xmath749 , _
k=_n=1^ ( ) ^k+e^-n/|n_j^-(k+ ) .
j7    the generating function @xmath750 has the following form in terms of @xmath751 : t_j(,z)=_n=1^(z)^k_k ( ) . j8inserting ( ) into ( ) and assuming that @xmath498 is a @xmath729 independent quantity , we find the following recurrent equation for @xmath498 : _ k=^k_k_1=1_k_1_k - k_1- 2_k .
j9therefore , if k > > j10then we can neglect last term in the right hand side of ( ) and in this case @xmath498 are positive and @xmath729 independent .
noting that in ( ) @xmath752 are essential the inequality ( ) means that the solution ( ) is correct if n>>|n_j|n_j , j11i.e . , only for this value of @xmath4 @xmath753 has the form ( ) and the corresponding generating functional has singularity at z_s=1 + . j12
it is important for the vhm experiment to have an upper restriction on the asymptotics .
we wish to show that @xmath70 decreases faster than any power of @xmath142 : _
n < o(1/n ) . 2.1to
prove this estimation , one should know the type of singularity at @xmath106 .
one can imagine that the points , where the external particles are created , form the system .
here we assume that this system is in equilibrium , i.e. there is not in this system macroscopical flows of energy , particles , charges and so on .
the lattice gas approximation is used to describe such a system .
this description is quite general and does not depend on details .
motion of the gas particles leads to the necessity to sum over all distributions of the particles on cells .
for simplicity we will assume that only one particle can occupy the cell .
so , we will introduce the occupation number @xmath754 in the @xmath174-th cell : @xmath755 means that we have no particle in the cell and @xmath756 means that a particle exists in a cell .
assuming that the system is in equilibrium , we may use the ergodic hypothesis and sum over all ` spin ' configurations of @xmath757 , with the restriction : @xmath758 .
it is evident that this restriction introduces the interactions .
the corresponding partition function in temperature representation ( , h)=de^-s ( ) 2.2awhere integration is performed over @xmath759 and , considering the continuum limit , @xmath760 . the action s_()=dx\{()^2 -^2+g^4- } 2.3where ~(1- ) ,  g~ ,   ~()^1/2h .
2.4and @xmath761 is the critical temperature .
0.2 cm    * unstable vacuum *    0.2 cm    we start this consideration from the case @xmath762 , i.e. assuming that @xmath763 . in this case
the ground state is degenerate if @xmath764 . the extra term @xmath765 in ( ) can be interpreted as the interaction with external magnetic field @xmath766 .
this term regulates the number of ` down ' spins with @xmath767 and is related to the activity : z^1/2=e^h , 2.5i.e .
@xmath766 coincides with the chemical potential .    the potential v()=-^2+g^4 ,  >0 , 2.6has two minima at @xmath768 if the dimension @xmath769 , no tunnelling phenomena exist .
but choosing @xmath770 the system in the correct minimum ( it corresponds to the state without particles ) becomes unstable .
the system tunneling into the state with an absolute minimum of energy .
the partition function @xmath185 becomes singular at @xmath764 because of this instability .
the square root branch point gives ( b , z)=e^-a_2()/h^2 ,  a_i>0 .
2.7note , @xmath771 at @xmath764 . deforming the contour in the mellin integral over @xmath129 on the branch line , _
n()=^_1 e^-4a_2 ^ 2/^2z .
2.8 in this integral z_c\{}^1/3 2.9is essential .
this leads to the following estimation : _ ne^-3(a_2 ^ 2)^1/3n^2/3<o(1/n ) . 2.10    it is useful to note at the end of this section that + ( i ) the value of @xmath69 is defined by @xmath772 and the metastable states , the decay of which gives a contribution into @xmath773 , are not important .
+ ( ii ) it follows from ( ) that in the vhm domain h~h_c~z_c~(1/n)^1/30 .
2.11so , the calculations are performed for the ` weak ' external field case , when the degeneracy is weakly broken .
it is evident that the life time of the unstable ( without particles ) state is large in this case and for this reason the semiclassical approximation is correct .
this is an important consequence of ( ) .
0.2 cm    * stable vacuum *    0.2 cm    let us consider now @xmath774 , i.e. @xmath775 .
the potential ( ) has only one minimum at @xmath776 in this case . the inclusion of an external field shifts the minimum to the point @xmath777 . in this case
the expansion in the vicinity of @xmath778 should be useful . as a result , ( , z)=\{dx _ c - w(_c ) } , 2.12where @xmath779
can be expanded over @xmath778 : w(_c)=_l_k\{dx_k_c(x_k;h)}_l ( x_1, ... ,x_l ) .
2.13 in this expression , @xmath780 is the one - particle irreducible green function , i.e. @xmath781 is the virial coefficient .
then @xmath778 can be considered as the effective activity of the correlated @xmath409-particle group .
the sum in ( ) should be convergent and , therefore , @xmath782 if @xmath783 .
but in this case the virial decomposition is equivalent to the expansion over the inverse density of particles . in the vhm region
it is high and the mean field approximation becomes correct . in result is _
c-()^1/3 :  |s_c|  if     to calculate @xmath69 . in this case z_ce^4gn^3  at  n , 2.17is essential and in the vhm domain _ ne^-4gn^4<o(e^-n ) .
2.18this result is an evident consequence of vacuum stability .
it should be noted once more that the conditions ( ) considerably simplify calculations .
* @xmath0-matrix unitarity constraints *      to explain our idea , let us consider the spectral representation of the one - particle amplitude : a_1(x_1 , x_2 ; e)= ,  + 0 , 3ait describes the transition of a particle with energy @xmath10 from point @xmath784 to @xmath785 .
according to our general idea , see introduction to sec.2.1 , we will calculate r_1(e)=dx_1 dx_2 a_1(x_1 , x_2 ; e)a^*_1 ( x_1 , x_2 ; e ) .
4 the integration over the end points @xmath784 and @xmath785 is performed only for the sake of simplicity .    inserting ( ) into ( ) and using ortho - normalizability of the wave functions
@xmath786 we find that & r_1(e)=_n||^2 = _
n\{-}= + & = _ n=_n(e - e_n ) . 5on other hand , the closed - path amplitude , offered for calculation in , & c_1(e)=_ndx=
_ n= + & = _ n\{p+ i(e - e_n)}=_np+ir_1(e ) .
4so , we wish to calculate only the imaginary part of the closed - path contribution : @xmath787 notice the extra factor @xmath788 in the left hand side .
the reason for this choice is evident : the real part of @xmath789 is equal to zero at @xmath790 , i.e. did not contribute to the measurable . to calculate the bound states energy spectrum ,
it is enough to know the only imaginary part of the closed - path amplitude .
this property is not accidental .
it is known as the optical theorem and is the consequence of the total probability conservation principles
. formal realization of it is the unitarity condition for the @xmath791-matrix : @xmath792 . in terms of the amplitudes @xmath793 , @xmath794
, the unitarity condition presents an infinite set of nonlinear operator equalities : i*a**a*^ * = * a * - * a*^*. 2notice that expressing the amplitude by the path integral one can see that the left hand side of this equality offers the double integral and , at the same time , the right hand side is the linear combination of integrals .
thus , the continuum contributions into the amplitudes should be canceled to provide the conservation of total probability . in this sense
it is a necessary condition .    indeed , to see the integral form of our approach ,
let us use the proper - time representation : a_1(x_1 , x_2 ; e)=_n _ n ( x_1)^*_n ( x_2)i ^_0dte^i(e - e_n+i)t 6 and insert it into ( [ 4 ] ) : r_1(e)=_n ^_0 dt_+dt_- e^-(t_++t_- ) e^i(e - e_n)(t_+-t_- ) .
7    we will introduce new time variables instead of @xmath795 : t_=t , 8 where , as follows from the jacobian of transformation , @xmath796 .
but we can put @xmath797 since @xmath798 is essential in the integral over @xmath799 . as a result , _
1(e)=2_n^_0 dt e^-2 t ^+_- e^2i(e - e_n ) .
9 in the last integral , the continuum of contributions with @xmath800 are canceled .
note that the product of amplitudes @xmath801 was ` linearized ' after the introduction of ` virtual ' time @xmath802 .
we wish to calculate the density matrix @xmath185 including the consequence of the unitarity condition cancelation of unnecessary contributions . here
we demonstrate the result and the intermediate steps we will formulate , without proof , as the statements offered in , where the formalities are described .          the statement , see and references cited therein , 0.1 cm _ s1 .
the unitarity condition unambiguously determines contributions in the path integrals for @xmath803 _ 0.1 cm looks like a tautology since @xmath804 , where @xmath805 is the action , is the unitary operator which shifts a system along the trajectory .
so , it seems evident that the unitarity condition is already included in the path integrals .
the rule as the path integrals should be calculated is weel known , see e.g. .
nevertheless the general path - integral solution contains unnecessary degrees of freedom ( unobservable states with @xmath806 in the above example ) .
we would define the path integrals in such a way that the condition of absence of unnecessary contributions in the final ( measurable ) result be loaded from the very beginning . just in this sense ,
the unitarity looks like the necessary and sufficient condition unambiguously determining the complete set of contributions .    _
s2 . the @xmath52- into @xmath4-particles transition ( unnormalized ) @xmath173 @xmath807 would have on the dirac measure the following symmetrical form _ : & r_nm(p_1, ... ,p_n , q_1, ... ,q_m)= < ^m_k=1|(q_k;u)|^2 ^n_k=1|(p_k;u)|^2>_u= + & = e^-i(j , e)dm(u)e^is_o(u)-iu(u , e ) ^m_k=1|(q_k;u)|^2 ^n_k=1|(p_k;u)|^2 + & ( u)^m_k=1|(q_k;u)|^2 ^n_k=1|(p_k;u)|^2
. 1.1ehere @xmath808 are the in(out)-going particle momenta .
it should be underlined that this representation is strict and is valid for arbitrary lagrange theory of arbitrary dimensions .
the eikonal approximation for inelastic amplitudes was considered in        the functionals @xmath813 and @xmath814 are defined by the equalities : s_o(u)=(s_0(u+e)-s_0(u - e ) ) + 2_c_+dx dte(x , t)(^2+m^2)u(x , t ) , u(u , e)= v(u+e)-v(u - e)-2_c_+ dx dte(x , t)v(u ) , 2.3awhere @xmath815 is the free part of the lagrangian and @xmath816 describes interactions .
the quantity @xmath814 is not equal to zero if @xmath817 have nontrivial topological charge ( see also ) .    according to s1 , considering motion in the phase space @xmath818 the measure @xmath819 has the dirac form : dm(u , p)= _ x , t du(x , t ) dp(x , t ) ( - ) ( + ) 2.4awith the total hamiltonian h_j(u , p)=dx \ { p^2 + ( u)^2 + v(u)-ju } . 2.5athis
last one includes the energy @xmath820 of quantum fluctuations .
* a. * only @xmath821 solutions of equations -=0 ,   + = 0 equwith @xmath466 should be taken into account .
this ` rigidness ' of the formalism means the absence of pseudo - solutions ( similar to multi - instanton , or multi - kink ) contribution .
the parameter @xmath826 is connected directly with @xmath827 particle energy , momentum , spin , polarization , charge , etc . , and is sensitive to the symmetry properties of the interacting fields system.@xmath828 , if @xmath829 is the magnetics order parameter , and the phase transition means that @xmath830 .
@xmath831 has just the same meaning as @xmath832 . ] for the sake of simplicity , @xmath833 is the real scalar field .
the generalization would be evident .    as a consequence of ( ) , @xmath826 is the function of the external particle momentum @xmath592 and is a @xmath834 functional of @xmath833 : ( q;u)=-dx e^iqx = dx e^iqx(^2 + m^2)u(x ) ,  q^2=m^2 , 1.2efor the mass @xmath52 field .
this parameter presents the momentum distribution of the interacting field @xmath833 on the remote hypersurface @xmath536 if @xmath833 is the regular function .
notice , the operator @xmath835 cancels the mass - shell states of @xmath833 .
the construction ( ) means , because of the klein - gordon operator and the external states being mass - shell by definition , the solution @xmath836 is possible for a particular topology ( compactness and analytic properties ) of @xmath837 field @xmath833 .
so , @xmath826 carries the following remarkable properties :          notice , the space - time topology of @xmath838 becomes important in calculating integral ( ) by parts .
this procedure is available if and only if @xmath838 is the regular function .
but the @xmath837 fields are always singular .
therefore , the solution @xmath839 is valid if and only if the quasiclassical approximation is exact .
just this situation is realized in the soliton sector of sin - gordon model .    despite evident ambiguity
@xmath826 carries the definite properties of the order parameter since the opposite solution @xmath836 can only be the dynamical display of an @xmath840 symmetry - matrix was introduced ` phenomenologically ' , see also the example considered in , postulating the lsz reduction formulae , see eq .
so , the formal constraints , e.g. the haag theorem , would not be taken into account on the chosen level of accuracy .
] , i.e. of the nontrivial topology of interacting fields , as the consequence of unbroken symmetry .
if ( ) have nontrivial solution @xmath41 , then this ` extended objects ' quantization problem arises .
we solve it by introducing convenient dynamical variables .
the main formal difficulty , see e.g. , of this program consists of transformation of the path - integral measure which was solved in
.    then + _ s3 .
the measure ( ) admits the transformation : u_c :  ( u , p)(,)w = g / g_c . o11and
the transformed measure has the form : dm(u , p)= _ x , t__c _ _
d(t ) d(t ) ( - ) ( + ) , o12where @xmath841 is the transformed hamiltonian .
_ 0.1 cm    it is evident that @xmath842 are parameters of integration of eqs .
( ) and they form the factor space @xmath843 . for instance , if one particle dynamics is considered , then one may choose @xmath844 and @xmath845 .
one may consider also the following possibility : @xmath846 and @xmath847 in this terms @xmath848 and new hamilton equations have the form : = 1-j ,   = j. o13so , we have at @xmath466 : @xmath849 and @xmath850 . by this reason
the ( action , angle)-type variables are mostly useful_.    according to ( ) there exists transformation of the perturbation generating operator : + _ s5 .
the operator @xmath851 has following transformed form : 2=dt \{__+ _ _ } , o14 in the factor space , where @xmath852 , @xmath853 , are new auxiliary variables . _    as a result of mapping of the perturbation generating operator @xmath851 on the manifold @xmath854 the equations of motion became linearized : dm=_t(--j _ ) ( -j _ ) .
o21then + _ s6 .
if the feynman s @xmath469-prescription is adopted , then the green function of eq .
( ) _ g(t - t)=(t - t ) o22    later on we will consider the soliton sector of sin - gordon model . in this case
@xmath855 is the coordinate and @xmath856 is the momentum of @xmath174-th soliton and @xmath545 is the number of solitons .
expansion of @xmath857 gives the ` strong coupling ' perturbation series .
its analysis shows that + _ s7 .
action of the integro - differential operator @xmath809 leads to the following representation _ : & r_nm(p , q)=_w\ { d(0)r^_nm(p , q).+ + & .
d(0)r^_nm(p , q)}. o20this means that the contributions into @xmath98 are accumulated strictly on the boundary ` bifurcation manifold ' @xmath858 , i.e. depends directly on the topology of @xmath854 .
performing the shifts in ( ): _
i ( t ) _ i ( t ) + dt g(t - t ) j_,i(t ) _ i ( t ) + _i ( t ) , + _
i ( t ) _ i ( t ) + dt g(t - t ) j_,i(t ) _ i ( t ) + _i ( t ) , 4.16we can get the green function @xmath859 into the operator exponent : ( ej ) = dt dt ( t - t)\{(t)_(t ) + (t)_(t)}. note the lorentz noncovariantness of our perturbation theory with green function ( ) .      using the definition :
@xmath860 the functional integrals on the measure ( ) are reduced to the ordinary integrals over initial data @xmath861 .
these integrals define zero modes volume .
notice that the zero - modes measure was defined without the faddeev - popov @xmath186 .
this strategy is necessary since it seems to be important to show the role of quantum corrections noting that for all physically acceptable field theories @xmath862 in the quasiclassical approximation .
the @xmath545-soliton solution @xmath863 depends upon @xmath864 parameters .
half of them , @xmath545 , can be considered as the position of solitons and the other @xmath545 as the solitons momentum . generally at @xmath865
the @xmath863 solution decomposed on the single solitons @xmath866 and on the double soliton bound states @xmath867 : @xmath868 note that this asymptotic is achieved if @xmath869 or / and @xmath870 .
this last one defines the bifurcation line of our model .
so , the one soliton @xmath866 and two - soliton bound state @xmath867 would be the main elements of our formalism .
its @xmath842 parametrizations have the form : u_s(x;,)=-\{(mx- ) } ,   = so and u_b(x;,)= -\ { } .
bo      in the quasiclassical approximation @xmath871 we have : u_n = u_n ( x;_0 , _ 0 + (_0)t ) .
5.7notice that the surface term dx^_(e^iqxu_n)=0 .
5.8then d^2x e^iqx(^2 + m^2)u_n ( x , t ) = -(q^2-m^2)d^2x e^iqxu_n ( x , t ) = 0 5.9since @xmath872 belongs to the mass shell by definition .
the condition ( ) is satisfied for all @xmath873 since @xmath863 belongs to schwarz space ( the periodic boundary condition for @xmath838 do not alter this conclusion ) .
therefore , in the quasiclassical approximation @xmath874 .
expanding the operator exponent in ( [ 5.1 ] ) we find that action of operators @xmath875 , @xmath876 creates terms ~d^2x e^iqx ( t - t ) ( ^2 + m^2)u_n ( x , t ) 0 .
5.12so , generally , if the quantum corrections are included , @xmath877 .
now we will show that the quasiclassical approximation is exact in the soliton sector of sin - gordon model .
the structure of the perturbation theory is readily seen in the ` normal - product ' form : r_2(q)=_n ^n_i=1 \{d_0 d_0}_i : e^-iu(u_n;/2i)e^is_0(u_n)|(q;u_n)|^2 : , 6.1where = _ - _ = _x 6.2and _ x = dt ( t - t)(t ) 6.3with the @xmath864-dimensional vector @xmath878 . in eq .
( ) @xmath879 is the ordinary symplectic matrix .
the colons in ( ) mean that the operator @xmath880 should stay to the left of all functions in the perturbation theory expansion over it .
the structure ( ) shows that each order over @xmath881 is proportional at least to the first order derivative of @xmath863 over the variable conjugate to @xmath882 .
the expansion of ( ) over @xmath883 can be written in the form ( omitting the quasiclassical approximation ) : r_2(q)=_n ^n_i=1 \{d_0 d_0}_i \{^2n_i=1p_x_i(u_n ) } , 6.4where @xmath884 is the infinite sum of ` time - ordered ' polynomials ( see ) over @xmath863 and its derivatives .
the explicit form of @xmath884 is complicated since the interaction potential is nonpolynomial .
but it is enough to know , see ( ) , that p_x_i(u_n)~_ij .
6.5      authors would like to take the opportunity to thank j. allaby , a. m. baldin , j. a. budagov , r. cashmore , g. a. chelkov , m. della negra , d. denegri , a. t. filippov , d. froidevaux , f. gianotti , i. a. golutvin , p. jenni , v. g. kadyshevski , o. v. kancheli , z. krumstein , e. a. kuraev , e. m. levin , l. n. lipatov , l. maiani , v. a. matveev , v. a. nikitin , a. nikitenko , a. g. olchevski , g. s. pogosyan , n. a. russakovich , v. i. savrin , m. v. saveliev , ju . schukraft , n. b. skachkov , a. n. tavkhelidze , v. m. ter - antonyan , h. t. torosyan , n.e .
tyurin , d. v. shirkov , e. sarkisian , s. tapprogge , o. i. zavialov , g. zinovyev for valuable discussions and important comments .
the scientific communications with them on the various stages were important .
we acknowledged also to members of the seminar ` symmetries and integrable systems ' of the n.n.bogolyubov laboratory of theoretical physics ( jinr ) , to the scientific seminars of the v. p. dzelepov laboratory of nuclear problems , d. v. scobeltsyn inp of the moscow state university , of the atlas experiment community for a number of important discussions .
we would like also to thank m. gostkin and n. shubitidze for the help in the formulation of the monte carlo simulation programs and for the discussions .
one of us ( j.m . ) was supported in part by georgian academy of sciences .
a.n.sissakian and l.a.slepchenko , preprint jinr , p2 - 10651 , 1977 ; + a.n.sissakian and l.a.slepchenko , fizika , 10 ( 1978 ) 21 ; + s. ch .
mavrodiev , v. k. mitryushkin , a. n. sissakian and h. t. torosyan , sov .
phys . , 30 ( 1979 ) 245      r.feynman , phys .
lett . , 23 ( 1969 ) 1415 ; + a.a.logunov and m.a.mestvirishvili , cern th/1659 ; + a.a.logunov , m.a.mestvirishvili , nguen van hieu , phys .
, 25b ( 1967 ) 611 ; + a.a.logunov and m.a.mestvirishvili , cern , th-1707 , 1973              g.giacomelli , nucl .
phys . 25b ( 1992 ) 30 ; + g.giacomelli , boll .
fis . , nuova ser . 6
( 1990 ) 20 , int
a5 ( 1990 ) 223 ; + h.d.wahl , acta phys .
, xxy ( 1983 ) 399 ; + a.b.kaidalov , survey in high energy physics , 7 ( 1994 ) 143 ; + p.g.ratcliff , proc .
8th int . symp . in high - en .
spin phys . , minneapolis , vol.1 ( 1988 ) 145 ; + d.grosh et al . ,
rev . , d49 ( 1994 ) 3113 ; + proc .
int . symp . in hadron interaction , eds . j.fischer , p.kolar , v.kundrat , czech , acad .
( 1988 ) 431                    e.kuraev , l.lipatov and v.fadin , sov .
jetp , 44 ( 1976 ) 443 ; zh .
fiz . , 71 ( 1976 ) 840
; + l.lipatov , sov.j .
20 ( 1975 ) 94 ; + v.gribov and l.lipatov , sov.j .
phys . , 15 ( 1972 ) 438 , 675 ; + g.altarelli and g.parisi , nucl .
phys . , b126 ( 1977 ) 298 ; + i.v.andreev , chromodynamics and hard processes at high energies ( nauka , moscow , 1981 )                                                                                i.p.kadanoff and g.baym , qunantum statistical mechanics ( benjamin , london , 1962 ) ; + g.niklasson and a.sjolander , ann.phys . , 49 ( 1968 ) 249 ; + p.danielevich , ann.phys . , 152 ( 1984 ) 239 + p.kadanoff and p.c.martin , ann.phys.(ny ) , 24 ( 1963 ) 419            e.wigner , phys .
, 40 , ( 1932 ) 749 ; + k.hisimi , proc .
jap . , 23 ( 1940 ) 264
; + r.j.glauber , phys .
10 ( 1963 ) 84 ; + e.c.g.sudarshan , phys . rev .
10 ( 1963 ) 177 ; + r.e.cahill and r.g.glauber , phys .
, 177 ( 1969 ) 1882 ; + s.mancini , v.i.manko and p.tombesi , quant .
, 7 ( 1995 ) 615 ; + v.i.manko , l.rose and p.vitale , hep - th/9806164