parameter estimation is a central issue in system modeling : a typical problem is to start from a certain amount of information on a given system ( e.g. its correlation functions ) and then extract the parameters of a model which is supposed to describe the its properties @xcite .
the parameter estimation procedure gives insight on the validity of the model and can suggest the introduction of more appropriate and efficient models .
a usual approach is to extract the parameters from an instance of the problem in certain conditions and subsequently testing the model in other instances . from this point of view
is useful to deal with systems in conditions where the relation between observables and model parameters is more transparent : e.g. for a statistical mechanics system this corresponds to high / low temperature or field .
once the parameters have been estimated , one moves to more interesting parameter regions , where the full complexity of the system shows up .
such an approach , when translated into the wide arena of complex systems , generally can not be carried out since no knob such as temperature or field is available , so that we may be faced with the inverse problem in the hardest region .
a huge interest in obtaining accurate parameter estimation stems from the current availability of large datasets in several areas of biology , economy and social sciences , to name a few examples dna sequences , stock market time series and internet traffic data ( see more references in @xcite ) .
this great amount of data has made even more pressing the quest for efficient models , allowing us to extract and encode the relevant information .
various techniques have been developed in order to solve this problem : two general approaches which can be flexibly adapted to the specific problems are bayesian model comparison @xcite and boltzmann - machine learning @xcite .    in the past decade a significant contribution to topic of parameter estimation came from the application of typical statistical mechanics techniques which turned out to be very useful in the modeling and study of different fields ranging from neurobiology @xcite to the economy @xcite .
the description of a system using statistical models ( and in particular ising - like models ) appears natural in many contexts : e.g. , effective ising models generally arise when the space of states is intrinsically discrete ( e.g. , for dna and proteins ) and , even when this is not the case , some ising variables may be lurking behind the continuous ones . in the statistical physics realm
such emergence of effective ising models could occur near a critical point when the microscopic model is in the ising universality class @xcite - but one can also find more subtle examples where discrete ising - like spin degrees of freedom describe some hidden order , e.g. the chiral ordering in frustrated continuous spin models @xcite .    a paradigmatic example considered by the statistical physics community in the context of parameter estimation is , of course , the inverse ising problem , i.e. the problem of inferring the values of the coupling constants of a general ising model from the correlation functions .
the inverse ising problem has been tackled by numerical and analytical methods , often adapting old techniques to the problem at hand . among these attempts
we mention monte carlo optimization @xcite , message passing based algorithms @xcite and thouless - anderson - palmer equations approaches @xcite ( see @xcite for a review ) .
field theoretical techniques has been used by sessak and monasson @xcite who perturbatively calculated , in terms of the correlations , expressions for the interaction parameters of a general ( heterogeneous ) ising model with two - body interaction and external field .
most of the available results on the inverse ising problem concerns ising models having two - spin interactions : in this context exact methods , solving the inverse ising problem with general multispin interactions , are welcome .
we decided to concentrate in this paper on the inverse ising problem in one dimension .
the motivation is three - fold : firstly one - dimensionality allows for exact solutions . in this manuscript
we indeed present explicit analytical formulas to exactly perform the inversion for one - dimensional ising systems having general multispin interactions .
our results therefore provide a theoretical laboratory where different approximate inverse ising techniques @xcite can be benchmarked against the exact results obtained using our method : in the following we compare some other approximate methods with exact results .
the possibility of testing approximate methods against exact results in one - dimensional systems is not our only motivation : indeed one - dimensional classical models are often employed to describe the conformational transition of systems , as proteins or dna , naturally possessing an underlying one - dimensional structure .
such simple models are found to capture some of the global properties of these complex systems as long as conformational properties are concerned .
the existence of exact methods would then help to determine the parameters and the important interactions of effective models describing the properties of such systems .
more in detail , the use of one - dimensional statistical mechanics models applied to systems like proteins or dna is usually based on the individuation of a reduced set of states representing the conformational state of a given elementary unit : e.g. in protein systems the states could be chosen as helix , coil and sheet ( for aminoacids belonging to an @xmath0-helix , to a coil and to a @xmath1-sheet , respectively ) .
the task is then , given this reduced set of states , to estimate the probabilities to have the consecutive elements in different states @xcite and then our method ( working for ising and potts models ) would then allow for the determination of the parameters of effective discrete models .
we notice that in our method we can consider also longer ranged couplings ( _ i_.e . , longer than nearest - neighbour ) emulating interactions among aminoacids distant along the chain , but near in physical space @xcite .
another motivation for our work is based on the fact that one - dimensional ising - like models can also be used to deal with stationary time - series of correlated data : as we will later discuss in the conclusions , it is possible to connect stationary time - series of data by using a mapping onto an equilibrium discrete markov chain having finite memory . for this application , the inversion task ( to which we can refer as an inverse markov problem )
consists in extracting from the data the transition probabilities of the associated guessed markov chain : therefore , given the similarity of the two inversion ( ising and markov ) problems , the existence of exact techniques can provide in perspective a way to effectively attack the inverse markov problem .
we observe that the method of using markov chains to describe sequences of data may prove useful even in biological realms when statistical properties of e.g. dna sequences are concerned @xcite .    in the following
we study the one - dimensional inverse ising problem with general finite - range multispin interactions : by finite range @xmath2 we mean that that two spins exceeding the distance @xmath2 do not interact ( this implies that at maximum @xmath2-spin couplings can be present ) .
we will then consider the reconstruction of ising models having exponential or power - law two - spin couplings ( and three- or four - spin interactions ) , approximating them with a finite range @xmath2 and checking the validity of the reconstructed couplings .
a remark about dimensionality is due : being the dimension set to one , the system can not order at finite temperature .
however we show that mean - field like interactions can be included in our formalism , so that one can treat systems having finite - range multispin couplings and long - range mean - field interactions giving rise to finite - temperature transitions
. another possibility would be to extend the range of the interaction and perform a so - called finite - range scaling .
such a technique has been employed @xcite to the ising model with power - law @xmath3 decaying interactions , a model exhibiting a rich behaviour including a berezinskii - kosterlitz - thouless transition ( for @xmath4 ) @xcite and gaussian and non - gaussian rg fixed points ( in the range of @xmath0 between @xmath5 and @xmath6 ) as the decay exponent @xmath0 is varied @xcite .    in this paper
we present the solution of the inverse problem for a one - dimensional ising model with finite - range arbitrary interactions , i.e. not restricted to the one- and two - body type .
the main result of our paper is formula which expresses the entropy of a one - dimensional translational invariant system ( in equilibrium ) in terms of a sufficiently large number of correlation functions , from which the inversion formula immediately follow .
we observe that ising chains are usually treated via transfer matrix method , but when longer range or multispin types of interaction are included the search for the parameters reproducing the observables might become very onerous .
our method provides a direct method of estimating the parameters when a sufficiently large number of correlation function is known .
the inclusion of many - body interactions may prove useful for the description of complex systems where the two - body assumption is not justified or in more traditional many - body systems with long - range interaction , where the construction of low - energy effective theories quite naturally leads to the appearance of multispin interactions @xcite .
the paper is structured as follows : in section [ notation_and_statement ] we introduce our notations and we state the mathematical problem .
section [ inversion_formula ] contains our main result on the entropy in terms of the correlation functions and the resulting inversion formula .
the obtained result is illustrated on simple problems in section [ simple_examples ] . in section [ realistic_examples ]
we examine more complicated examples where the usefulness of our result is shown .
we analyze models formally not have finite - range interactions and having exponential or power - law two - spin interactions plus multispin interactions .
the data generated by monte carlo simulations are analyzed with our technique which correctly detects the structure of interactions . in section [ mf_lr ]
we briefly discuss how the developed formalism may be modified to allow mean - field interactions .
finally we draw our conclusions in [ conclusions_and_perspectives ] .
the appendix present checks of the obtained findings for small values of the range @xmath2 using the transfer matrix method , and as well as supplementary material on the @xmath7 ising model .
we consider a general one - dimensional ising model with multispin interactions defined by the hamiltonian @xmath8 @xmath9 where @xmath10 is the configuration of the @xmath11 ising spins ( @xmath12 ) ; periodic boundary conditions will be assumed so that @xmath13 for @xmath14 .
the sums runs over distinct couples , triples and so on ; the temperature dependence is absorbed in the coupling constants : explicitly , @xmath15 , @xmath16 , and so on ( where _ e_.g .
@xmath17 is the two - body coupling among a spin in @xmath18 and a spin in @xmath19 - as usual @xmath20 ) .
the couplings @xmath21 are assumed to be invariant under translation by @xmath22 spins ( for simplicity we will assume @xmath23 is an integer , but since we are interested in the @xmath24 limit this is not strictly necessary ) : this condition reads @xmath25 ( if the indices on the right hand side exceed @xmath11 , they have to be replaced by the indices equivalent modulo @xmath11 contained in the set @xmath26 ) .
finally we assume that the couplings are zero if their indices can not be brought by a translation of a multiple of @xmath22 to a subset of @xmath27 .
since the use of the form of the hamiltonian may be cumbersome , it is convenient introduce a more compact notation , rewriting hamiltonian as @xmath28 where @xmath29 is a subset of @xmath30 ( this is encoded in the writing @xmath31 , which stands for `` the range of the interaction is less or equal than @xmath2 '' ) .
@xmath22 is the periodicity of the interaction and @xmath32 is an operator associated to the subset @xmath33 ( @xmath34 is the number of elements of @xmath29 ) translated by @xmath35 which acts on the spins as @xmath36 for the null subset @xmath37 we define @xmath38 . the prime in the sum over @xmath29 in indicates that the null subset ( which would contribute just to a constant in the hamiltonian ) is not included and that the terms related by a translation of a multiple of @xmath22 are counted only once , in order to avoid the presence of equivalent operators in the hamiltonian .
once the hamiltonian is specified we proceed in the usual calculation of the thermodynamic quantities , defining the partition function @xmath39 the free energy per elementary unit cell in the infinite volume limit ( i.e. @xmath22 spins ) @xmath40 and the correlation functions associated to the operator @xmath29 @xmath41 ( by definition , @xmath42 ) .
the inverse problem for the system introduced in the previous section is stated as follows : given the set of correlations @xmath43 determine the couplings @xmath44 .
the hamiltonian is the one specified in equation , i.e. the most general finite - range multispin hamiltonian ; in section [ mf_lr ] we will extend this treatment to include long - range mean - field interactions .
the calculation is based on the evaluation of the entropy per unit cell @xmath45 characterized by the set of correlation functions @xmath43 .
once @xmath45 is known we may compute the free energy @xmath46 where @xmath47 is the energy of a unit cell which is readily evaluated using directly and on a state specified by the set of correlations @xmath43 .
the minimization of the above expression yields the inversion formulas : @xmath48    we may state now our main result for the entropy @xmath45 , which is given by @xmath49 the entropy is written in terms of the functions @xmath50 ( to which we may refer as `` entropy at range @xmath51 '' ) , given by @xmath52 where @xmath53 is the configuration of @xmath51 auxiliary ising spins .
notice that the sum over @xmath29 now includes every subset , including the null one .
the entropy can be shown to be convex in the variables @xmath43 , thus the equation admits a solution , unless some of the @xmath54 s used in the calculation become negative , signaling a set of `` nonphysical '' correlations .
we now discuss the derivation of the formula .
let us think how the measurement of a correlation @xmath55 is operatively defined : we look at @xmath2 consecutive spins and we perform the measurement .
each of the microscopic configuration @xmath56 will occur with a given probability @xmath57 which would give rise to a mean value of @xmath55 given by @xmath58 since we know all of the correlations within the subsets of the @xmath2 spins , the system of the equations above may be inverted giving rise to with @xmath59 . then the boltzmann formula @xmath60 is applied to this set of probabilities obtaining the expression @xmath61 ( always for @xmath59 ) . to derive we calculate the entropy of the unit cell of size @xmath22 , regardless of the state of the remaining @xmath62 spins ; in terms of number of states it is @xmath63 where @xmath64 denotes the number of microstates of a set of spins @xmath65 ( subject to the constraints imposed by the correlations ) .
it should be noted that the @xmath62 spins to be traced out can not be chosen at will : by inspection it turns out that the first @xmath62 spins is a good choice . thus taking the logarithm of we obtain our expression for the entropy of a state characterized by the set of correlations @xmath66 .
the number of correlations required to specify the state can be shown by simple counting to be equal to @xmath67 .
the above procedure can be formally applied also to a finite system of size @xmath2 : this is achieved by letting @xmath68 . in this case
the system looks like a set of @xmath69 disjoint assemblies of @xmath2 spins for which the entropy is given by @xmath70 being @xmath71 .
this result refers to a general finite system of ising spins with arbitrary heterogeneous couplings without translational invariance and it should be used if one wants to treat datasets obtained from finite heterogeneous systems and extract ising couplings @xcite . in general , the number of needed correlations to be known grows exponentially with the system size : in the case of the present paper , exponentially with @xmath2 .
therefore , our result might be of practical importance if @xmath2 is small , @xmath11 is very large ( _ i.e. _ , near the thermodynamical limit ) and the underlying ising model is supposed to be one - dimensional .
although the number of required correlation grows exponentially with the range @xmath2 , some simplifications may occur .
for example if the system is known to be invariant under reflection @xmath72 , then odd couplings vanish and we do not have to measure the corresponding correlation function which are trivially zero . more generally , if we know that a coupling is equal to a given value @xmath73 then one has the additional ( nonlinear ) equation : @xmath74 thus reducing the number of independent correlation functions .
the technique described may be easily adapted to other discrete spin systems as potts or blume - emery - griffith models by using two or more ising spins to encode the state of the discrete spin , although such a mapping may obscure the symmetries of the original model .    to conclude this section
, we finally observe that it is possibler to show in general the equivalence between the problem we have solved with the problem of finding , from a set of known correlations , the transition rates of an equilibrium - i.e. satisfying detailed balance - markov chain with finite memory : we postpone such discussion to the conclusions .
legs may be mapped onto a spin chain with @xmath75 and @xmath76 .
the filled circles represent spins which are present in the different operators considered ( and associated to the couplings @xmath77 ) . in the middle part of the figure
we represent them on the ladder system , while on the bottom part we show how they look like on the chain .
explicitly these operators correspond to the following terms in the hamiltonian : @xmath78 , @xmath79 , @xmath80 , @xmath81 , @xmath82 and @xmath83     legs may be mapped onto a spin chain with @xmath75 and @xmath76 .
the filled circles represent spins which are present in the different operators considered ( and associated to the couplings @xmath77 ) . in the middle part of the figure we represent them on the ladder system , while on the bottom part we show how they look like on the chain .
explicitly these operators correspond to the following terms in the hamiltonian : @xmath78 , @xmath79 , @xmath80 , @xmath81 , @xmath82 and @xmath83    as a first application of the results presented in the previous section , we consider a model with @xmath84 and @xmath85 ( i.e. the hamiltonian is @xmath86 where @xmath87 is the magnetic field and @xmath88 is the coupling ) .
the only independent correlations are the one - body correlator , i.e. the magnetization @xmath89 , and the nearest - neighbour correlator @xmath90 . using the entropy
is calculated as @xmath91 which agrees with the expression obtained in @xcite by combinatorial means . in figure [ entropy_g_m ]
we plot the entropy : the convexity of @xmath92 guarantees to obtain the field and nearest - neighbour interaction in terms of @xmath93 and @xmath94 .
the system we have just described presents no phase transitions , apart from the zero temperature ones which occur at the border of the surface depicted in figure [ entropy_g_m ] ; in section [ mf_lr ] we will see how the addition of a mean - field type interaction is easily included , making phase transitions possible .
interestingly on the lines @xmath95 the system is frustrated and our approach readily provides an expression for the ground state degeneracy .
differentiation of the entropy allows us to obtain the couplings , field @xmath96 and nearest - neighbour interaction @xmath97 conjugated to @xmath93 and @xmath94 respectively : @xmath98 in appendix [ transfer_matrices ] we examine this example ( @xmath84 and @xmath85 ) and higher range ones ( @xmath99 and @xmath85 ) , explicitly checking the validity of the inversion formula using the transfer matrix method .
we will consider now a translationally invariant spin ladder with interaction among the nearest - neighbors of the same and other chain .
for simplicity we will restrict ourselves to even interactions , i.e. in the hamiltonian only terms containing an even number of spins enter .
as shown explicitly in figure [ ladder_to_chain ] this system may be mapped onto a chain system with @xmath75 and @xmath76 , where the original interactions ( allowed by symmetry ) and the new ones are shown . in terms of our subset notation used in the interaction parameters are defined as @xmath100 .
this is easily generalized to ladders made up of more than two chains and higher interaction range , and thus our method is suited to treat general finite - range ladder systems .
we point out that the inversion formula allows to explicitly write the relation among the @xmath88 s and @xmath94 s while the transfer matrix approach , e.g. in the simple ladder system described above , already entails the solution of a fourth order algebraic equation .      in this section
we consider another simple example , in which nearest - neighbour and next - to - nearest - neighbour interactions are present together with a four - spin interaction : denoting the coefficients @xmath101 , @xmath102 and @xmath103 by @xmath104 , @xmath105 and @xmath106 respectively , the hamiltonian reads @xmath107 this hamiltonian can be exactly treated in our framework ; the case @xmath108 ( the @xmath7 model ) is discussed in the appendix . here
we aim at comparing approximate inverse ising methods against exact results , focusing in particular on the low - correlation expansion ( lce ) @xcite which is in the following compared with exact findings .
we will use the lce discussed in @xcite using as input a finite number of correlations , coherently with what is done in this work ( notice that for the present case our method needs just four correlation functions in order to recover exactly the couplings ) ; the maximal range of two - body correlators for the lce is denoted by @xmath109 .
the lce is discussed in @xcite to present an approximate technique for inverse ising models having at most two - spin interactions : since hamiltonian has only two - spin interactions for @xmath110 , we present lce results for the case @xmath110 in the appendix where we discuss the @xmath7 model in detail , presenting the explicit solution using the transfer matrix approach .
as expected , for low temperatures ( i.e. large couplings ) , the lce breaks down and , as it can be seen in the right panel of figure [ j1j2sm ] , for moderate temperatures the expansion may settle to an incorrect value of the coupling as the range @xmath109 is increased . in order to further test the performance of the lce against exact results we present in figure [ cv_chi ] results for @xmath110 and @xmath111 : although the lce is developed in @xcite for two - spin interactions ( and the extension to treat multispin interactions is expected to be cumbersome ) , the lce reconstructed two - spin couplings with @xmath112 may partially take into account the effect of the four - spin interaction . to test to
what extent this occurs , we consider two observables , susceptibility and specific heat , calculated using the reconstructed couplings : the comparison with the exact results is in figure [ cv_chi ] .
as we can see in the left panel the specific heat is more sensitive than the susceptibility to reconstruction errors , even without four - spin interaction ( _ i_.e .
@xmath110 ) .
this may be traced back to this type of lce inversion procedure which aims at reproducing two body correlators which in this model are required to calculate the susceptibility while the specific heat already contains averages of four body operators .
obviously the lce , being not designed to infer models with multispin interaction , gives no hint on the value @xmath106 but in the right panel of figure [ cv_chi ] we apply it for @xmath113 in order to test how it can reproduce the considered observables anyway : we can see that the addition of such an operator reduce the temperature range where the observables are correctly reproduced . as noted above the specific heat
is more subjected to errors than the susceptibility at contains higher - body operators . from numerical inspection
we saw that the lce rather well performs in the high temperature regime even for relatively large values of @xmath106 , but deviates from exact results at lower temperature even for small values of @xmath106 as shown in figure [ cv_chi ] .
is varied .
the value @xmath114 is held fixed at the value @xmath5 .
the figure on the right includes the four - spin interaction whose coupling is fixed at the value @xmath115 .
the inset shows the ratio between the predicted values of specific heat and inverse susceptibility and the exact ones .
both figures are obtained keeping @xmath116 correlation functions and using third order loop resummed expansion to reconstruct the couplings , as presented in @xcite.,title="fig : " ]   is varied .
the value @xmath114 is held fixed at the value @xmath5 .
the figure on the right includes the four - spin interaction whose coupling is fixed at the value @xmath115 .
the inset shows the ratio between the predicted values of specific heat and inverse susceptibility and the exact ones .
both figures are obtained keeping @xmath116 correlation functions and using third order loop resummed expansion to reconstruct the couplings , as presented in @xcite.,title="fig : " ]
[ reconstruction ]     ( empty symbols ) and the values of the couplings really used to generate the correlations ( filled circles ) - the full line is a guide to the eye .
the figures refers to hamiltonian with parameters @xmath117 , @xmath118 , @xmath119 ( left ) and @xmath120 , @xmath121 , @xmath122 ( right).,title="fig : " ]   ( empty symbols ) and the values of the couplings really used to generate the correlations ( filled circles ) - the full line is a guide to the eye .
the figures refers to hamiltonian with parameters @xmath117 , @xmath118 , @xmath119 ( left ) and @xmath120 , @xmath121 , @xmath122 ( right).,title="fig : " ]     ( top ) and the inferred couplings @xmath123 ( bottom ) for the @xmath29 s allowed by translational symmetry ( filled blue circles ) and exact values ( empty purple triangles ) .
the black diamonds refer to lce results obtained with the perturbative expansion up to third order @xcite .
the figure at the bottom is a logarithmic plot of the absolute values of the @xmath123 s , while the inset gives for comparison the linear plot of the same @xmath123 s .
the figures refer to hamiltonian with parameters @xmath117 , @xmath118 .
the reconstruction range is @xmath124 . in abscissa
are reported the @xmath29 s
denoting the various coupling and correlators with the subset notation introduced in section [ notation_and_statement ] , e.g. @xmath125 is the correlation of the subset @xmath126 , @xmath127 is the correlation of the subset @xmath128 and so on . ]    in this section we consider examples where we can not access the full knowledge of our system : our inversion procedure will therefore yield approximate results .
first , we consider an ising model with an exponentially decaying two - body interaction @xmath129 since the interaction now is not formally of `` finite range '' i.e. it does not vanish for distances beyond a given value of @xmath2 , the transfer matrix method is not viable ( although we still may perform a finite - range scaling in the size of the transfer matrix @xcite ) .
the set of synthetic correlation functions is generated by a monte carlo method .
of course we will not record all of the correlation function , but we will fix a maximal range @xmath130 , thus we will have to measure on the order of @xmath131 correlation functions .
the results for such a reconstruction are shown in figure [ reconstruction ] for two values of the parameters .
we see that the agreement improves as the value of @xmath130 is increased ( at the expense of calculating a larger number of correlation functions ) . in figure [ gj1 ] a full set of the correlations and inferred couplings
is shown ; if we look at the lower panel , the nonzero couplings are clearly singled out ( even for a value of @xmath130 as low as @xmath132 ) , thus our reconstruction procedure gives useful hints to build a faithful model of an unknown system . in figure [ gj1 ] results obtained with the lce
are also reported : one sees that there is a good agreement for the considered value of the temperature between the lce results and the findings obtained using our reconstruction procedure .    in order to test the procedure on a system with more - body couplings we consider the hamiltonian @xmath133 which includes three- and four - body interactions .
as we see in figure [ gj2 ] , even in this case the reconstruction procedure gives useful hints on the couplings present in the system , although some of the inferred couplings , which were zero in the starting model , are predicted to be of comparable size to the nonzero ones ( especially the @xmath134 coupling ) .
this is due to the finite reconstruction range and to the , albeit small ( indeed smaller than the symbols in the figures [ gj1 ] , [ gj2 ] , [ gj3 ] ) , errors in the determination of the correlation .
this implies that in order to clearly distinguish the contribution of the different couplings , the correlation should be known with high accuracy . as for modeling purposes , this is not a problem since the values of the coupling obtained give rise to a set of correlations not distinguishable from the original one . for reference we also
plot in figure [ gj2 ] the results of the tap equation approach developed in @xcite , which of course provides no information on the multispin couplings , but as far as one - body operators and two - body operators are concerned this approach at the temperature considered in figure [ gj2 ] performs very well .
finally we examine a model with power - law decay of the interaction @xmath135 as can be seen in figure [ gj3 ] ( with @xmath136 ) , the results are good also in this case : it can be generally observed that the reconstructed interactions are higher than the exact ones , due to the fact that the interaction within the reconstruction range have to account for the interactions lying outside this range .     except for the black diamonds which are obtained with the tap approach @xcite .
the figures refer to hamiltonian with parameters @xmath137 , @xmath118 , @xmath138 , @xmath139 .
the reconstruction range is @xmath124 .
the arrows mark the multispin interactions . ]    .
measured values of the correlations @xmath55 ( top ) and inferred couplings @xmath123 ( bottom ) for the power - law decaying hamiltonian with parameters @xmath140 , @xmath136 . the reconstruction range is @xmath124 . ]
in this section we briefly discuss how our previous results can be used in presence of mean - field long - range interactions , showing that our inversion approach may be used on this class of systems .
we consider a system with energy @xmath141 of general form , i.e. a non - linear function of the correlators .
the number of couplings entering the energy @xmath141 should still equal the number of independent correlation functions in order to perform , at least in principle , the inversion procedure .
such an energy will be denoted by @xmath142 where the index @xmath93 runs over the mean - field couplings .    by requiring the free energy @xmath143 to have a minimum when the @xmath144 are set to the known values will give the equations implicitly determining all the couplings , including the mean - filed ones @xmath145 . when the energy is differentiable these equations read @xmath146 this set of equations will in general have multiple solutions or possibly no solutions at all .
if more solutions are present the one ( or ones ) rendering the free energy minimal should be chosen .
the points where the absolute minimum of the free energy branches or it changes discontinuously will signal a phase transition . as described in the section [ inversion_formula ] some values of the couplings may be known in advance thus reducing the number of equations to be solved .
another possibility is that a function of the coupling is fixed ; a notable example being the energy itself , corresponding to the microcanonical description of the system .
we remark that in the class of models we consider all of these steps can be carried out exactly since the explicit form of entropy is known .    as an example we may consider the first model examined in section [ simple_examples ] ( @xmath84 , @xmath85 ) by adding mean - field two - body couplings . instead of the energy @xmath147 ( where @xmath148 and @xmath149 ) we will set @xmath150 .
the appearance of the nonlinear term @xmath151 in the energy is due to the presence of non - local mean - field operators in the spin hamiltonian like @xmath152 .
it should be noted that such an operator is not uniquely defined , e.g. the operator @xmath153 and other kac - rescaled nonextensive potentials give rise to the term @xmath154 in the energy density at the thermodynamic limit when evaluated on a state with magnetization @xmath93 and nearest - neighbour correlation @xmath94 @xcite .
if we set @xmath155 we get the usual mean - field model , otherwise we obtain a model with competing mean - field and short - range coupling introduced by kardar @xcite , exhibiting a complex phase diagram which shows nonequivalence between the canonical and microcanonical description @xcite ( for a review on inequivalence between ensembles and other issues concerning nonextensive systems see @xcite ) . as already observed the entropy of such a model within our approach
is easily computed [ see ] and it could be generalized thus allowing one to treat one - dimensional models possessing multiple competing finite - range and mean - field interactions .
in this paper we presented the explicit solutions of the inverse ising problem for a one - dimensional translational invariant model with arbitrary finite - range multispin interactions once a number @xmath156 ( where @xmath2 is the range of the interactions ) of independent correlations is known . when applied to unknown systems
this method correctly detects arbitrary interactions ; our results are then applied to systems with a range extending beyond the one set by maximum distance of the spins of the recorded correlation functions , giving useful hints on the interactions that should be kept in an effective model .    as an application
, we reconstructed the couplings of chain ising hamiltonians having exponential or power - law two - spin plus three- or four - spin couplings .
we also discussed the generalization of the method to ladders .
mean - field interactions can be also included in the framework , allowing us to describe systems exhibiting phase transitions .
the presence of both finite - range ( local ) and mean - field ( nonlocal ) interactions can give rise to interesting competition effects greatly enhancing the descriptive power of the models we can exactly solve with our techniques .
our results provide then a theoretical laboratory where different approximate inverse ising techniques can be benchmarked against the exact results obtained using our method : in the paper we performed such comparison is some illustrative examples .
the one - dimensional inverse ising problem we have solved in the present paper is analogous to what may be called the inverse markov chain problem : given a specific set of correlations at equilibrium , find the corresponding transition rates .
the two inversion ( ising and markov ) problems are related since it is possible to associate an ising model to an equilibrium markov ising chain with finite memory in full generality : to show it , for definiteness let us consider ising variables ( although extensions to other discrete state spaces is straightforward ) .
the finite - range @xmath2 in our solution of the inverse one - dimensional ising problem is the counterpart of the finite memory in the inverse markov problem : let the state of the next @xmath22 spins be ruled by the state of the preceding @xmath62 spins .
we put the new spins on the right side of the old ones : the time of the markov process is increasing from left to right .
the correlations impose constraints on the transition rates : it turns out that the number of independent correlations required to solve the inverse markov problem is the same as the number needed to solve the ( related ) inverse ising problem in one dimension of range @xmath2 and period @xmath22 , i.e. the number of independent correlations is @xmath67 .
adapting the procedure discussed in section [ inversion_formula ] which lead to it is possible to compute the transition rates from the state @xmath157 of @xmath62 spins to the state @xmath158 of @xmath22 spins .
these transitions are given by @xmath159 where the @xmath160 s are calculated according to formula where the input correlations @xmath66 s are plugged in , and @xmath161 is the configuration of @xmath2 spins obtained by juxtaposing the states @xmath157 ( on the left ) and @xmath158 ( on the right ) . the transition matrix obtained this way is a @xmath162 by @xmath163 matrix ; in order to have a square matrix we have to fold more steps of the the transition matrix until we obtain the probability to go from a set of @xmath164 to the next @xmath164 spins @xcite .
the transition matrix satisfies detailed balance by construction , therefore this is a reversible markov chain .
this mapping has already been worked out in a discrete different form in @xcite where the connection to discrete statistical models is also discussed .    according to the previous discussion , we may thus associate an ising model to the an equilibrium markov ising chain with finite memory in full generality , allowing us to treat systems where one direction ( typically time ) is singled out . as a possible
example deserving future investigation one could consider time - series of financial data and try to estimate with the procedure discussed before the transition probabilities of the associated guessed markov chain in order to test the validity of such a description . it is intended - in order to apply the previous results - that the analyzed data should be discretized on a timescale such that nontrivial correlations occur , and that the whole time of observation is such that the system can be reliably considered at equilibrium .
obviously a way to encode significant information in an ising variable has to be devised , being this in general a nontrivial task ; for example , we may think of `` up '' and `` down '' spin corresponding to a price raise or lowering respectively . the next step would be to analyze the correlations among different time - series of data to determine if and how correlations among different stocks occur .
for example the ising ladder system depicted in figure [ ladder_to_chain ] may reproduce the correlations among two stocks whose state depends on the value of the other stock at the same timestep and on the value of the same or other stock at the previous timestep ( for simplicity , in figure [ ladder_to_chain ] odd interactions are not depicted ) .
, it should be noted that the inclusion of the interaction dubbed @xmath165 in figure [ ladder_to_chain ] may reproduce some kind of nontrivial many - body interaction among the stocks .
extending the number of chains in the ladder system and/or the range allows us to treat larger sets of stocks with longer correlations in time .
we think that studying stationary time - series of correlated data using the techniques presented in this paper ( and as well the mapping onto markov chains discussed in this section ) will be an interesting subject of future research . in perspective
, one could apply the method here discussed to datasets and/or statistical mechanics models which are supposed to be described by effective one - dimensional ising chains near the thermodynamical limit . to this aim
, one should address in the future a treatment of the case where large errors in the measured correlation functions are present and/or some of the correlations are missing ; our exact result could be a good starting point to move in that direction .
next , our results could be extended in higher dimensional cases ( where of course one does not expect to find closed formulas ) , hierarchical or tree - like models .
some preliminary results in the two dimensional case seem to indicate that this approach leads to equations resembling the dobrushin - lanford - ruelle ones @xcite .
another interesting direction could be to use the a renormalization group approach on the correlation functions , in order to study how the couplings determined by correlations at some scale @xmath2 are related to the ones computed at a larger scale @xmath166 .
_ acknowlegments : _ we wish to thank i. mastromatteo and m. marsili for many very useful discussions .
this work has been supported by the grants instans ( from esf ) and 2007jhlpez ( from miur ) .
in this appendix we work out the transfer matrix for the general translational invariant ( @xmath85 ) ising model with range @xmath167 and we check that the results obtained analytically and numerically with the transfer matrix formalism are fully consistent with the predictions of the formula .
we start by briefly recalling the method ( see e.g. @xcite ) . in general the transfer matrix @xmath168
is built by identifying the @xmath169 states of a block of spins @xmath170 with independent and orthogonal vector of a space of dimension @xmath169 such that the matrix elements of @xmath168 are @xmath171 where @xmath172 is the interaction energy among two consecutive blocks of spins @xmath173 and @xmath174 ( @xmath173 is placed to the left of @xmath174 ) , @xmath175 is the interaction energy among the spins belonging to the same block .
the vector corresponding to the spin state @xmath173 is denoted , using the ket notation , by @xmath176 .
the size of the blocks @xmath170 has to be chosen according to the range and size of the unit cell @xmath22 of the system in order to have all of the interaction terms contained in @xmath177 or @xmath178 .
the partition function of the system of size @xmath11 is simply given by @xmath179 so that in the infinite size limit the free energy per unit cell may be written in terms of the largest eigenvalue @xmath180 of @xmath168 : @xmath181 the existence and unicity of @xmath180 is guaranteed by perron - frobenius theorem , being all elements of @xmath168 strictly positive ( for nonvanishing temperature ) .
the correlation functions may be obtained just by differentiation of the free energy @xmath182 .
we start with the simple example of @xmath84 which has been worked out in section [ simple_examples ] .
the two independent couplings will be denoted as usual by @xmath183 ( magnetic field ) and @xmath184 ( nearest - neighbour coupling ) .
the @xmath6 by @xmath6 transfer matrix @xmath168 , with the states identification @xmath185 and @xmath186 , reads @xmath187 the free energy is then given by @xmath188.\ ] ] differentiating @xmath189 with respect to @xmath87 and @xmath88 gives the magnetization @xmath93 and the nearest - neighbour correlation @xmath94 respectively : @xmath190 the inversion of the above formulas with respect to @xmath87 and @xmath88 yields the expressions .
now we consider a model with only even interactions and range @xmath191 , i.e. containing the coupling @xmath192 and @xmath193 which we dub @xmath7 model : the corresponding hamiltonian is given by with @xmath110 .
this model may be mapped onto the previous example ( @xmath84 , @xmath85 ) by introducing `` kink '' variables @xmath194 with the identifications @xmath195 , @xmath196 ( we do not report the corresponding results ) . in figure [ j1j2sm ]
we plot the couplings reconstructed according to the low - coupling expansion ( lce ) introduced in @xcite against the exact results which can be found by the method discussed in this paper or by the transfer matrix approach .
the lce allows to infer the magnetic fields and the two - body couplings from the two - body correlators and magnetizations . in @xcite
the lce has been carried out , in the zero field case , up to seventh order in the correlations with loop resummation .
we have used the lce as both the order of the expansion and the number of correlation function we assume to know are increased .
the maximal range of two - body correlations which are used as input is denoted by @xmath109 . in figure [ j1j2sm ]
we depict the reconstructed couplings @xmath197 , @xmath198 and the exact ones for different values of @xmath109 .
the ratio @xmath199 is kept fixed , thus @xmath104 serves as inverse temperature . as expected as the temperature is lowered the agreement gets worse , and it may be noticed that the inclusion of higher order terms in this case does not significantly improve the performance of the inversion : as one can see , the lower order results depicted in the left panel of figure [ j1j2sm ] is more reliable at lower temperatures ( this is the reason why in figure [ cv_chi ] we employed the third order lce ) . in the right panel of figure [ j1j2sm ]
we can see how the the increase of @xmath109 improves the quality of the inversion but beyond a given range the reconstructed couplings @xmath200 settle to a value which , in the lower temperature case examined , deviates from the exact one .    .
the various lines denote how many correlation functions are kept in the reconstruction formulas as indicated .
the continuous and dotted line refer respectively to the third and seventh order in the lce of @xcite ( including the loop contributions ) . on the right the ratio of reconstructed @xmath201 ( above ) ,
@xmath202 ( below ) and the exact value of @xmath104 in terms of the reconstruction range @xmath109 are shown .
only the result for the seventh order lce is shown . in both figures
the green straight lines are the exact values and the nearest - neighbour coupling @xmath105 is set to the half of @xmath104.,title="fig : " ] .
the various lines denote how many correlation functions are kept in the reconstruction formulas as indicated .
the continuous and dotted line refer respectively to the third and seventh order in the lce of @xcite ( including the loop contributions ) . on the right the ratio of reconstructed @xmath201 ( above ) ,
@xmath202 ( below ) and the exact value of @xmath104 in terms of the reconstruction range @xmath109 are shown .
only the result for the seventh order lce is shown . in both figures
the green straight lines are the exact values and the nearest - neighbour coupling @xmath105 is set to the half of @xmath104.,title="fig : " ]    we move on to the next example , the @xmath191 case with no restriction on the symmetry of the couplings . identifying the states as @xmath203 , @xmath204 , @xmath205 , @xmath206 and the couplings as @xmath207 , @xmath208 , @xmath209 , @xmath210 , the transfer matrix reads @xmath211    proceeding as before , we can obtain the entropy and correlation functions .
the comparison between the resulting obtained findings and our inversion formula are shown in the left panel of figure [ tm4x48x8 ] .
the plot shows the numerically calculated entropy and the analytical one where the numerical correlations are plugged in , for some values of the coupling constants .
the entropy in this case is given by the expression ( the subscripts in the @xmath94 s just indicate to which coupling they are conjugated ) : @xmath212    as a last case we consider a translation invariant chain with range @xmath75 : identifying the states as @xmath213 , @xmath214 , @xmath215 , @xmath216 , @xmath217 , @xmath218 , @xmath219 , @xmath220 and the couplings as @xmath207 , @xmath208 , @xmath209 , @xmath221 , @xmath222 , @xmath223 , @xmath224 , @xmath225 , the @xmath226 by @xmath226 transfer matrix reads : @xmath227 where the @xmath228 matrices @xmath229 are @xmath230 results are shown in figure [ tm4x48x8 ] ( we do not write down the entropy for this case ) . by inspecting figure [ tm4x48x8 ] we find that the agreement between the two approaches is complete .
we do not report here the other checks we performed for higher values of @xmath2 .
( left ) for the transfer matrix and @xmath231 ( right ) for the transfer matrix .
the values of the remaining couplings are specified on the @xmath232-axes of the plots .
the arrows point to the relevant @xmath232-axis .
the continuous curves are obtained by numerical calculation , while the dots are calculated by means of .,title="fig : " ]   ( left ) for the transfer matrix and @xmath231 ( right ) for the transfer matrix .
the values of the remaining couplings are specified on the @xmath232-axes of the plots .
the arrows point to the relevant @xmath232-axis .
the continuous curves are obtained by numerical calculation , while the dots are calculated by means of .,title="fig : " ]


pion multiplicity distribution and their associated fluctuation and correlations have been of interest for several reasons .
several models predict large fluctuations such as the disoriented chiral condensate @xcite and in density fluctuations from droplets arising in a first order phase transition @xcite .
a well known procedure for studying correlations uses the bose - einstein symmetries associated with pions in a hanbury brown - twiss analysis @xcite .
such an analysis gives information about the space time history of the collision through measurements of source parameters .
if the density of pions becomes large , bose - einstein correlation may also lead to a strongly emitting system which has been called a pion laser @xcite .
the pion laser model has been recently solved analytically by t. csrg and j. zimanyi @xcite .
the importance of bose - einstein correlations has also been illustrated in the observation of a condensation of atoms in a harmonic oscillator or laser trap @xcite .
previous interest in pionic distributions have centered around the possibility of intermittency behavior @xcite and fractal structure based on parallels with turbulent flow in fluids .
a distribution widely used to discuss such features has been the negative binomial ( nb ) distribution @xcite with its associated clan structure @xcite and kno scaling feature @xcite .
kno scaling properties have been interpreted in terms of a phase transition associated with a feynman - wilson gas @xcite .
various other issues associated with pions include evidence for thermalization @xcite , critical point fluctuations @xcite , fluctuations from a first order phase transition @xcite , charge particle ratios and question of chemical equilibrium @xcite , the behavior of fluctuations in net charge in a qg plasma for transition @xcite .    for lower energy heavy ion collisions ,
multifragmentation of nuclei takes place .
the fragment distribution can also be described statistically by considering all the possible partition of @xmath4 nucleons into smaller clusters @xcite .
this study gives a tool for the description of nuclear multifragmentation distributions @xcite , nuclear liquid - gas phase transition @xcite , critical exponent , intermittency , and chaotic behavior @xcite of nuclear multifragmentation .
the same model can describe pionic distribution .
this possibility arises in our approach which is based on feynman path integral methods where symmetrization of @xmath4 bosons or anti - symmetrization of @xmath4 fermions leads to a cycle class decomposition of the permutations associated with these symmetries .
the correspondance comes from the identification of clusters of size @xmath5 and the cycles of length @xmath5 in a permutation as discussed below . in nuclear multifragmentation
the number of nucleons , @xmath4 , is fixed and thus a canonical partition function approach is appropriate . by contrast , the number of pions , @xmath4 , in particle production is not fixed and thus we need to use a grand canonical partition . using this parallel , we also note that some results from cluster yields can be carried over into particle multiplicity distributions and associated properties .
namely , we show the importance of the fisher critical exponent @xmath3 and relate it to one of the parameters called @xmath0 in our approach which has three main parameters @xmath0 , @xmath1 , @xmath2 .
moreover , taking special values of @xmath0 , or equivalently @xmath3 , reduces our unified model to various specific cases that are frequently used in particle production phenomenology .
quantities that appear in this development can also be related combinants as will be discussed . in so doing we can give a physical significance to the combinants and show how the fisher exponent appears in them and how the resulting hierarchical structure and void scaling relations also depend on its value .    in next section a summary of a generalized statistical model based on a grand canonical partition function
will be given .
this generalized model can be used to study multiplicity distributions associated with particle production such as at rhic .
various models of particle multiplicity distribution using the grand canonical partition function will then be developed . in section ii - b
, we will summarize the multiplicity distribution and its various moments of all order for a general grand canonical ensemble which we will then use for particle production .
the physical meaning of the parameters used in this model is given here .
the differences and relations between the canonical nuclear multifragmentation and the grand canonical multiparticle production are sumarized in sect .
ii - c . in sect .
ii - d and e , moments ( especially the mean and the fluctuation ) in the general grand canonical model are related to variables used in other existing standard models in describing multiparticle production .
the connection of the general grand canonical model developed here to existing standard models are summarized in sect .
moreover , we derive a generalized model ( hga ) of the grand canonical partition which can further be reduced to a geometric , negative binomial , and lorentz / catalan model .
various moments and multiplicity distribution for a generalized hga model are summarized in sect .
ii - g and another generalized model ( grw1d ) in sect .
ii - h . in sect .
[ sectcomp ] , we compare various statistical properties between different models within the generalized hga model and with data for charged particle distribution in jets of l3 and h1 .
fluctuation and multiplicity distribution in hga are compared with other standard models of multiparticle production .
consider a system composed of @xmath6 different types of species or objects which could be the fragments in a fragmentation or in a cycle class description .
any event of such a system can be associated with a vector @xmath7 or @xmath8 where the non - negative integer @xmath9 is the number of individuals of species @xmath5 . for example @xmath9 can be the number of clusters of size @xmath5 or the number of cycles of length @xmath5 in a given permutation of @xmath10 particles .
the later is important for bose - einstein and fermi - dirac statistics and particle multiplicity distributions .
a general block picture of @xmath11 is shown in fig.[fig1]a .
fig.[fig1]b shows how the various partition can be developed as an evolution from successively smaller systems .
the number of species @xmath6 can be infinity in general .
( 25,8)(7.5,0.0 )    ( 0.5,0.5 )    ( 12.2,6.2)*a ) *    ( 8.0,0.25)*@xmath12 *    ( 8.01 , 1.51)(0,1)4.5 ( 8.0 , 5.5)(0.5,0)2 ( 8.0 , 5.0)(0.5,0)2 ( 8.0 , 4.5)(0.5,0)4 ( 8.0 , 4.0)(0.5,0)5 ( 8.0 , 3.5)(0.5,0)5 ( 8.0 , 3.0)(0.5,0)7 ( 8.0 , 2.5)(0.5,0)7 ( 8.0 , 2.0)(0.5,0)7 ( 8.0 , 1.5)(0.5,0)11 ( 8.99 , 5.0)(0,1)1.0 ( 9.99 , 4.5)(0,1)0.5 ( 10.49 , 3.5)(0,1)1.0 ( 11.49 , 2.0)(0,1)1.5 ( 13.49 , 1.5)(0,1)0.5 ( 8.01 , 1.51)(1,0)5.5 ( 11.51 , 2.01)(1,0)2.0 ( 10.51 , 3.51)(1,0)1.0 ( 10.01 , 4.51)(1,0)0.5 ( 9.01 , 5.01)(1,0)1.0 ( 8.01 , 6.01)(1,0)1.0 ( 7.5,6.5)(0,0.01)2(1,0)0.5 ( 9.5,6.5)(0,0.01)2(-1,0)0.5 ( 8.0,6.2)(1.0,0.6)*@xmath9 * ( 7.6,3.2)(0,0.01)2(0,-1)1.7 ( 7.6,4.3)(0,0.01)2(0,1)1.7 ( 7.4,3.5)(0.3,0.6)[r]*@xmath5 * ( 15.5,0 )    ( 10,8)(0,0 )    ( 0.5,0.5 )    ( 8,6.2)*b ) *    ( 3.5 , 7.0 ) ( 4.2,7.3)@xmath13 ( 3.6,6.9)(-3,-2)1.05 ( 3.9,6.9)(2,-3)0.85    ( 2.0 , 5)(0,0.5)2 ( 2.7,5.6)@xmath14 ( 2.05,4.9)(-2,-1)0.65 ( 2.4,4.9)(3,-2)1.1    ( 4.5 , 5.0)(0.5,0)2 ( 5.7,5.4)@xmath15 ( 4.52,4.9)(-1,-1)0.78 ( 5.45,4.9)(3,-4)0.95    ( 1 . , 3.0)(0,0.5)3 ( 1.7,4.0)@xmath16 ( 1.1,2.9)(-1,-1)0.75 ( 1.4,2.9)(1,-4)0.3    ( 3.5 , 3.0)(0,0.5)2 ( 4.0 , 3.0 ) ( 4.2,3.8)@xmath17 ( 3.6,2.9)(-4,-3)1.6 ( 3.9,2.9)(0,-1)1.75 ( 4.4,2.9)(3,-4)1.3    ( 6.0 , 3.0)(0.5,0)3 ( 7.1,3.7)@xmath18 ( 6.4,2.9)(-1,-3)0.6 ( 7.3,2.9)(1,-2)1.1    ( 0 . , 0.0)(0,0.5)4 ( 0.7,1.6)@xmath19    ( 1.5 , 0.0)(0,0.5)3 ( 2.0 , 0.0 ) ( 2.15,1.1)@xmath20    ( 3.5 , 0.0)(0,0.5)2 ( 4.0 , 0.0)(0,0.5)2 ( 4.7,0.7)@xmath21    ( 5.5 , 0.0)(0.5,0)3 ( 5.5 , 0.5 ) ( 6.2,0.8)@xmath22    ( 8.0 , 0.0)(0.5,0)4 ( 9.6,0.7)@xmath23    various probability distributions related with this system can be developed by assigning an appropriate weight @xmath24 to each type @xmath5 .
a weight @xmath25 is then given to each event @xmath11 and the type of weight that will be considered has the structure : @xmath26        \label{weight}\end{aligned}\ ] ] the @xmath27 are gibbs factorials .
such a weight structure appears in feynman s path integral approach and bose - einstein problems @xcite .
the @xmath24 will be given below and contains various physical quantities . in sect.[sectiib ] we will show that the @xmath24 s are also the combinants which in turn can be related to the factorial cumulants .
summing the weight @xmath25 over all the possible events of @xmath11 , the grand canonical partition function @xmath28 of the system is given as @xmath29     = \exp\left[\sum_{k=1}^n x_k\right ]          \label{grandptf }   \end{aligned}\ ] ] the last equation holds due to the form of eq.([weight ] ) of the weight @xmath25 , i.e. , the factor @xmath30 is the @xmath9-th order expansion term of @xmath31 .
introducing other quantities @xmath32 to each individual entity or group of type @xmath5 , @xmath33 , we can define a canonical partition function @xmath34 with a fixed @xmath4 as @xmath35              \label{canoptf }    \end{aligned}\ ] ] with @xmath36 . here
@xmath37 is the summation over all events with a fixed value of @xmath4 , i.e. , over a canonical ensemble of a fixed @xmath4 , and the @xmath38 is a summation over all the possible values of @xmath4 ; it becomes @xmath39 for the case of @xmath40 with positive integer @xmath5 .
the multiplicity of a partitions is @xmath41 the case @xmath40 is encountered in fragmentation problems and permutation problems .
various type of @xmath32 can be used depending on the physics related with the quantity @xmath4 as discussed in ref.@xcite .
the physics of the canonical ensemble depends on the choice of the quantity @xmath32 @xcite .
if we take @xmath42 , then the canonical partition function @xmath34 is the @xmath43 dependent term of the grand canonical partition function @xmath28 .
if we take @xmath44 , then @xmath1 counts the number of clusters @xmath45 explicitly and the @xmath46 dependent term of the grand canonical partition function @xmath28 is the partition function for events with multiplicity @xmath45 .
there always is at least one event having @xmath47 , i.e. , the event where all @xmath9 s are zero , @xmath48 .
thus @xmath49 if all @xmath32 are non zero positive since then there are no other possible events having @xmath47 .
due to the form of the weight @xmath25 given by eq.([weight ] ) the canonical partition function @xmath34 satisfies a recurrence relation @xcite @xmath50 this relation is nothing but the constraint eq.([asum ] ) in terms of the mean @xmath51 using eq.([factma ] ) . for non - zero positive @xmath32
, there is no case having @xmath52 , i.e. , @xmath53 for @xmath52 and thus the @xmath54 can be obtained by the recurence relation of eq.([recur ] ) starting from @xmath49 .      in a canonical ensemble of fixed @xmath4
, we can define a probability distribution of a specific partition @xmath11 as @xmath55 with this probability , various mean values , fluctuations and correlations of the number of species @xmath9 can be evaluated as a ratio of canonical partition functions for two different values of @xmath4 such as @xcite @xmath56 thus we have @xmath57 .
the recurrence relation of eq.([recur ] ) then follows simply using the fact that @xmath58 .
this distribution has been used in describing various fragment distributions in nuclear multifragmentation @xcite .    now knowing the partition functions @xmath28 and @xmath34 , we can associate a probability @xmath59 of the system to have a fixed value of @xmath4 in a grand canonical ensemble as @xmath60_{z=0 }                    \label{pax }   \\
z(\vec x , z ) & = & \sum_a z_a(\vec x ) z^a      = \exp\left[\sum_k x_k z^{\alpha_k}\right ]    \label{grandzxz}\end{aligned}\ ] ] the last step follows from the fact that the @xmath43 power term of @xmath61 is @xmath34 if we put @xmath62 . for particle production @xmath59
is the probability of having @xmath4 particles .
thus the generating function @xmath61 of @xmath63 can also be looked at as a grand canonical partition function with the weight @xmath24 replaced to be @xmath64 , where the variable @xmath2 counts @xmath4 explicitly and @xmath65 . if we consider @xmath66 then @xmath2 has two roles ; one as a weight which is assigned the same to each constituent and another as a generating parameter of the probability @xmath59 . for the case
that @xmath49 , the void probability @xmath67 @xcite is the inverse of the grand canonical partition function , i.e. , @xmath68 inversely the canonical partition function @xmath34 is the probability @xmath63 normalized or rescaled by the void probability @xmath67 , i.e. , @xmath69_{z=0 }     = \frac{p_a(\vec x)}{p_0(\vec x)}\end{aligned}\ ] ] another type of generating function of @xmath63 that is frequently used may be defined as @xmath70 \right ]             \label{genu }    \\
p_a(\vec x ) & = & \frac{1}{\gamma(a+1 ) } \left[\left(-\frac{d}{d u}\right)^a          g(\vec x , u)\right]_{u=1}\end{aligned}\ ] ] we see that @xmath71 , @xmath72 and @xmath73 .
once the probability @xmath59 is determined , various statistical quantities can be evaluated .
also the grand canonical partition function @xmath74 is given once we know the thermodynamic grand potential @xmath75 this result can be used to study the statistical properties of a system .
moreover various moments and cumulants , mean values and fluctuations may be obtained using the generating function @xcite .
this subsection gives general expression for various quantities that will be used later when we discuss specific models .
since the probability of a specific event @xmath11 in a grand canonical ensemble is given by @xmath76 the mean of a quantity @xmath77 in a grand canonical ensemble is related to the mean of @xmath77 in a canonical ensemble as @xmath78 we can easily show that @xmath79 this result shows that the weight factor @xmath24 in this model is the mean number @xmath80 in a grand canonical ensemble and that the mean multiplicity @xmath81 is related to the thermodynamic potential .
the @xmath82-th power moment of @xmath4 and its factorial moments are given simply by @xmath83_{z=1 }        \\
\left<\frac{\gamma(a+1)}{\gamma(a - m+1)}\right>(\vec x )    & \equiv & \sum_{\vec n } \frac{\gamma(a+1)}{\gamma(a - m+1 ) } p(\vec x , \vec n )     = \frac{1}{z(\vec x ) } \left[\left(\frac{d}{d z}\right)^m            z(\vec x , z)\right]_{z=1 }       \end{aligned}\ ] ] similarly the @xmath82-th cumulants , which is the power moments of @xmath32 , and the factoral cumulants are @xmath84_{z=1 }     = \sum_{k=1}^\infty \alpha_k^m x_k      \label{alpham }   \\
f_m(\vec x ) & \equiv & \left <
\sum_{k=1}^n            \frac{\gamma(\alpha_k+1)}{\gamma(\alpha_k - m+1 ) } n_k\right >     = \left[\left(\frac{d}{d z}\right)^m \ln z(\vec x , z)\right]_{z=1 }                \nonumber    \\    & = & \left[\left(- \frac{d}{d u}\right)^m \ln g(\vec x , u)\right]_{u=0 }     = \sum_{k=1}^n \frac{\gamma(\alpha_k+1)}{\gamma(\alpha_k - m+1 ) } x_k                  \label{factfm } \end{aligned}\ ] ] due to eqs.([grandzxz ] ) and ( [ grandpot ] ) we can see easily that , for the power moments of @xmath32 , @xmath85 the power moments of @xmath32 are directly related to the power moments of @xmath4 measured from the mean @xmath86 , i.e. , the cumulants @xmath87 are same with the central moments of @xmath4 .
this simple relation does not hold for @xmath88 but we can evaluate them starting from @xmath89 using the recurrence relation @xmath90 similarly the @xmath82-th factorial cumulants @xmath91 , which is the factorial moments of @xmath32 , can be found using recurrence relation @xmath92 starting from @xmath93 the reduced factorial cumulants @xmath94 defined in ref.@xcite corresponds to the factorial cumulants @xmath91 normalized with mean number @xmath95 as @xmath96 thus with @xmath97 and @xmath98 . the cummulants are directly related with the generating functions as @xmath99     = \exp\left[\sum_{m=1}^\infty \frac{(\ln[1-u])^m}{m ! }
\bra\alpha^m\ket(\vec x)\right ]                      \\
z(\vec x , z ) & = & \exp\left[\sum_{m=0}^\infty \frac{(z-1)^m}{m ! } f_m(\vec x)\right ]      = \exp\left[\sum_{m=0}^\infty \frac{(\ln z)^m}{m ! } \bra\alpha^m\ket(\vec x)\right]\end{aligned}\ ] ] the generating function @xmath61 differs from the generating function @xmath100 only by an extra term of @xmath101 in their exponent .
the relation between the @xmath24 s and the @xmath102 shows that the @xmath24 s are also the combinants of ref.@xcite . in the approach presented here ,
the combinants are given an underlying significance through the partition weight @xmath25 of eq.([weight ] ) .
in turn the combinants @xmath24 can be related to the factorial cumulants @xmath91 defined by @xmath103 the factorial cumulants @xmath91 are the @xmath82-th order factorial moments of @xmath32 of eq.([factfm ] ) .
thus @xmath104 for @xmath40 which is the eq.([factfm ] ) with @xmath105 .
the normalized factorial cumulant , i.e. , the reduced cumulant , is @xmath106 for a negative binomial ( nb ) distribution , this result of eq.([kappam ] ) shows that @xmath94 for nb has an hierarchical structure of a distribution at the reduced cumulant level which was realized for the nb distribution in ref.@xcite .
this result will be generalized later .    also using the above power moments and factorial moments we can study voids and void scaling relation , hierarchical structure , combinant and cummulant properties which will be discussed below .
since our approach was first used to discuss multifragmentation and then later extended to include multiparticle production , we briefly mention some of difference between multifragmentation and multiparticle production .
we also show how the fisher exponent @xmath3 , which initially appeared in cluster yields around a crititical point , may manifest itself in particle production yields .    in nuclear multifragmentation
, @xmath40 is the number of nucleons in a fragment and @xmath9 is the number of fragments of size @xmath5 .
the total number of fragments is @xmath107 and the total number of nucleons is @xmath108 . in nuclear multifragmentation ,
the total number of nucleons @xmath4 is usually fixed and we study the distribution in size @xmath5 of the mean multiplicity @xmath109 of fragments of size @xmath5 in a canonical ensemble . in multiparticle production
the @xmath4 is the total number ( multiplicity ) of produced particles and is not fixed and we study the multiplicity distribution of produced particles @xmath110 as a function of @xmath4 in a grand canonical ensemble .
the mean multiplicity @xmath111 and @xmath51 in a canonical ensemble are related with grand canonical ensemble as @xmath112 the weight @xmath24 can be determined experimentally from the mean multiplicity @xmath80 of cluster size @xmath5 in a grand canonical ensemble .
in our initial fragmentation studies @xcite we have used @xmath113 .
this choice gives @xmath114 and thus @xmath115 as @xmath116 where the @xmath117 with @xmath118 the chemical potential . from eq.([factma ] ) , @xmath119 , thus @xmath120 as @xmath116 in the case of nuclear fragmentation in the canonical ensemble of fixed @xmath4 for @xmath113 and the yields fall as a power law @xmath121 at a critical point ( @xmath122 ) for this choice of @xmath24 .
the yield in the grand canonical ensemble is @xmath123 and this becomes the same as @xmath120 in @xmath116 limit for @xmath124 with @xmath125 . in general ,
cluster yields fall for large @xmath5 as a power law @xmath126 at a critical point where @xmath3 is the fisher critical exponent .
the above choice of @xmath24 can be generalized to @xmath127 so that the cluster yields fall as @xmath126 .
the @xmath3 determines the grand canonical partition function @xmath128 = \exp[x \sum_k z^k /k^\tau]$ ] .
forms for @xmath24 whose asymptotic behavior is @xmath129 can also be used . in section ii.f
we will give various choice for @xmath24 for particle production .
for example , the above choice of @xmath24 with @xmath130 gives the frequently used negative binomial distribution with @xmath1 the negative binomial parameter which determines the degree of depature from poisson statistics .
other distributions can be generated and their associated @xmath3 dependence will also be given .
then we will show how the @xmath3 dependence of the combinants shows up in hierarchical structure and void scaling relations in sect ii.g and in other places .
the physical meaning of @xmath1 and @xmath2 will also be discussed for these other cases . since an exact description of particle multiplicity yields is not known , we consider @xmath3 as a free exponent to be determined by comparing various resulting distributions with experiment .
this procedure is what is done in cluster yields . in multi - fragmentation distribution ,
the mean number of clusters @xmath81 is determined by @xmath1 and @xmath2 is directly related with the number of constituent particles @xmath4 .
the large @xmath5 behavior of the mean number of cluster @xmath80 gives the fisher power @xmath3 .
in multiparticle distribution , @xmath1 and @xmath2 determine the mean particle number @xmath86 or the peak position of the particle number distribution @xmath63 and the fluctuation @xmath131 or the width of the particle distribution .
the fisher power @xmath3 determines the large @xmath5 behavior of the weight @xmath24 and offers a new parameter in particle production .
@xmath132 in grand canonical ensemble .
bose - einstein condensation of atoms in a box of sides @xmath133 of dimensions @xmath134 have @xmath135 with @xmath136 and have @xmath137 .
feynman used random walk arguments , the closing of a cycle parallels a closed random walk , to discuss his choice of @xmath24 in his discussion of a superfluid phase transition in liquid helium .
thermal emission of pions based on statistical mechanics and equilibrium ideas have been popular descriptions of pions coming from relativistic heavy ion collisions . for thermal models
@xcite , the @xmath138 for a cycle length @xmath5 or a cluster of size @xmath5 with @xmath139 a mac donald function . for low temperatures , @xmath140 and the boltzmann factor in mass , @xmath141 , suppresses large fluctuations . in the high temperature limit and/or zero pion mass limit @xmath142 .
thus the high temperature limit has @xmath143 and the low temperature limit has @xmath144 . the @xmath24 can be used to generate the pion probability distribution @xmath145 .
the thermal models can be combined with hydrodynamic descriptions and an application was given @xcite to @xmath146 gev pb+pb data measured by the cern / na44 and cern / na49 collaborations .
the results of ref.@xcite showed a gaussian distribution with a width about 20 % larger than the poisson result .
van hove and giovannini have introduced clan variables @xmath147 and @xmath148 to describe a general class of probability distributions , with most discussions of these variables centering around the negative binomial distribution @xcite .
these variables are defined as @xmath149 where the mean number of clans is @xmath147 and the @xmath148 is the mean number of members per clan .
the @xmath102 is the grand canonical generating function and thus @xmath150 where @xmath24 is the cycle class weight distribution @xmath151 .
the @xmath152 is the mean number of total members ( particles ) .
the @xmath9 here is the number of clans of size @xmath5 having @xmath5 members .
the clan variable @xmath147 is also related to the void probability @xmath153 ; thus @xmath154 . an important function in void
analysis is @xmath155 .
thus the void parameters , void probability @xmath67 and void function @xmath156 @xcite , are equivalent to the generalized clan parameters @xmath147 and @xmath148 with the equivalence given by : @xmath157 the @xmath158 is the normalized grand potential @xmath159 for the mean @xmath86 .    void analysis looks for scaling properties associated with @xmath160 ; specifically , @xmath160 is a function of the combination @xmath161 where @xmath162 is the coefficient of @xmath163 in the fluctuation @xmath164 .
since @xmath165 , the variance of @xmath4 in a grand canonical ensemble becomes @xmath166\end{aligned}\ ] ] with @xmath167 , i.e. , the normalized factorial cumulant .
since the variance for a poissonian distribution is the same as the mean , @xmath168 , the @xmath169 ; thus , the parameter @xmath170 represents a degree of departure from poissonian fluctuation normalized by mean @xmath86 of the distribution .
a well known non - poissonian example is a nb distribution which has @xmath171 and this becomes plank distribution with @xmath172 .
thus the variable @xmath1 in @xmath24 determines the strength of the non - poissonian fluctuation term for the negative binomial distribution . in general , @xmath1 and @xmath2
determine the mean particle number or the peak position of the particle multiplicity distribution and the fluctuation or the width of the particle distribution . using the recurrence relation eq.([recfm ] ) for @xmath91
, we can show that @xmath173          - 2 \frac{\xi(\vec x , z)}{\bra a\ket }   \end{aligned}\ ] ] a nb distribution has @xmath174 while the lorentz / catalan ( lc ) distribution discussed in ref .
@xcite and below has @xmath175 .
we will study in sect .
[ voidscal ] , @xmath160 _ vs _ @xmath176 , i.e. , the void or clan variable _ vs _ the fluctuation for various choices of @xmath24 summarized in table [ tabl1 ] .      the lc model , which has @xmath177 , was shown to be a useful model for discussing an underlying splitting dynamics when ancestral or evolutionary variables @xmath178 and @xmath179 were introduced into @xmath1 and @xmath2 as discussed in ref.@xcite .
specifically @xmath180 and @xmath181 giving @xmath182 where @xmath183 .
percolation or splitting dynamics with a branching probability @xmath178 and survival probability @xmath184 has a hierarchical topology as shown in fig.[lcfig ] .
weighting each diagram by @xmath185 , the evolutionary or ancestral variables are related to the clan variables @xmath186 and @xmath187 .
@xmath188 is the number of diagrams of size @xmath5 shown in fig.[lcfig ] . for this case
the evolutionary dynamics is just that of the lc model .
then with @xmath179 set equal to 1 , @xmath189 , @xmath190 , @xmath191 , @xmath192 , etc .
the interpretation of this set of @xmath24 s reads as follows : @xmath193 has 1 surviving line without a branch ( @xmath194 ) and one diagram ( @xmath195 ) , @xmath196 has 1 branch ( @xmath197 ) leading to 2 surviving lines ( @xmath198 ) and one diagram ( @xmath199 ) , @xmath200 has 2 branch points ( @xmath201 ) leading 3 surviving lines ( @xmath202 ) and two diagrams ( @xmath203 ) , @xmath204 has 3 branch points(@xmath205 ) , 4 surviving lines ( @xmath206 ) and 5 diagrams ( @xmath207 ) , etc . in these evolutionary / ancestral variables the @xmath208 which determines @xmath102 is @xmath209 for all @xmath210 . for @xmath211 ,
@xmath212 is no longer a constant and is @xmath213 . to keep @xmath214 a constant without changing @xmath215 ,
a @xmath216 was introduced in ref.@xcite .
the @xmath217 for @xmath210 and @xmath218 for @xmath211 . for @xmath211
there is a finite probability ( @xmath219 ) that the splitting will go on forever ( @xmath220 ) . in percolation above a certain @xmath178 an infinite cluster is formed and @xmath221 is similar to the strength of the infinite cluster .
moreover the sudden appearance of @xmath221 is similar to the behavior of an order parameter in a phase transition .
the appearance of @xmath221 can also be interpreted as the sudden appearance of a jet without pion ( @xmath222 ) , if we introduce @xmath223 instead of @xmath216 .    since the clan variables are @xmath224 and @xmath225 , then , for the lc model with evolutionary or ancestral varables @xmath179 and @xmath178 , @xmath226        \nonumber \\
\beta & = & n_c \frac{n_c - ( 1 \pm 1)/2}{n_c - 1 }                 \nonumber\end{aligned}\ ] ] these can be reduced to @xmath227 , @xmath228 , and @xmath229 for @xmath210 while @xmath230 , @xmath231 , and @xmath232 for @xmath233 .
since the branching probability @xmath178 varies in @xmath234 , the clan variable @xmath148 has @xmath235 with @xmath236 at @xmath237 and 1 , only one member per clan in average , and @xmath238 at @xmath239 , infinitly many members per clan .
the lc model thus connects the clan variable @xmath148 to the probability @xmath178 of branching in the evolutionary or ancestral picture of fig.[lcfig ] or in a percolation model . for poisson processes @xmath240 ( no branching ) ,
@xmath241 ( only unit cycles and no be correlations ) and @xmath236 ( one member in each clan in average ) .    (
35,7)(1.2,4.0 ) ( 3,7)(0,1)4 ( 2,6.0)@xmath242 ( 1,4.5)unit cycles ( 8.5,9)(0,1)2 ( 8.0,9)(1,0)1 ( 8.5,9 ) ( 8.7,9.4)@xmath178 ( 8.0,7)(0,1)2 ( 9.0,7)(0,1)2 ( 7.0,6.0)@xmath243 ( 5.6,4.5)period 2 cycles ( 13.2,10)(0,1)1 ( 12.5,10)(1,0)1.5 ( 13.2,10 ) ( 13.4,10.4)@xmath178 ( 12.5,7)(0,1)3 ( 14.0,8.5)(0,1)1.5 ( 13.5,8.5)(1,0)1 ( 14.0,8.5 ) ( 14.2,8.9)@xmath178 ( 13.5,7)(0,1)1.5 ( 14.5,7.0)(0,1)1.5 ( 16.8,10)(0,1)1 ( 16.0,10)(1,0)1.5 ( 16.8,10 ) ( 17.0,10.4)@xmath178 ( 17.5,7)(0,1)3 ( 16.0,8.5)(0,1)1.5 ( 15.5,8.5)(1,0)1 ( 16.0,8.5 ) ( 16.2,8.9)@xmath178 ( 15.5,7)(0,1)1.5 ( 16.5,7.0)(0,1)1.5 ( 13.0,6.0)@xmath244 ( 12.2,4.5)period 3 cycles ( 20.9,10)(0,1)1 ( 20.2,10)(1,0)1.5 ( 20.9,10 ) ( 21.1,10.4)@xmath178 ( 20.2,7)(0,1)3 ( 21.7,9.0)(0,1)1.0 ( 21.2,9.0)(1,0)1 ( 21.7,9.0 ) ( 21.9,9.4)@xmath178 ( 21.2,7.0)(0,1)2.0 ( 22.2,8.0)(0,1)1.0 ( 21.7,8.0)(1,0)1 ( 22.2,8.0 ) ( 22.4,8.4)@xmath178 ( 21.7,7.0)(0,1)1.0 ( 22.7,7.0)(0,1)1.0 ( 24.1,10)(0,1)1 ( 23.4,10)(1,0)1.5 ( 24.1,10 ) ( 24.3,10.4)@xmath178 ( 23.4,7)(0,1)3 ( 24.9,9.0)(0,1)1.0 ( 24.4,9.0)(1,0)1 ( 24.9,9.0 ) ( 25.1,9.4)@xmath178 ( 24.4,8.0)(0,1)1.0 ( 25.4,7.0)(0,1)2.0 ( 23.9,8.0)(1,0)1 ( 24.4,8.0 ) ( 24.6,8.4)@xmath178 ( 23.9,7.0)(0,1)1.0 ( 24.9,7.0)(0,1)1.0 ( 27.4,10)(0,1)1 ( 26.7,10)(1,0)1.5 ( 27.4,10 ) ( 27.6,10.4)@xmath178 ( 26.7,8.5)(0,1)1.5 ( 28.2,8.5)(0,1)1.5 ( 26.2,8.5)(1,0)1 ( 26.7,8.5 ) ( 26.9,8.9)@xmath178 ( 26.2,7)(0,1)1.5 ( 27.2,7.0)(0,1)1.5 ( 27.7,8.5)(1,0)1 ( 28.2,8.5 ) ( 28.5,8.9)@xmath178 ( 27.7,7)(0,1)1.5 ( 28.7,7.0)(0,1)1.5 ( 30.8,10)(0,1)1 ( 30.0,10)(1,0)1.5 ( 30.8,10 ) ( 31.0,10.4)@xmath178 ( 31.5,7)(0,1)3 ( 30.0,9.0)(0,1)1.0 ( 29.5,9.0)(1,0)1 ( 30.0,9.0 ) ( 30.2,9.4)@xmath178 ( 29.5,7.0)(0,1)2.0 ( 30.5,8.0)(0,1)1.0 ( 30.0,8.0)(1,0)1 ( 30.5,8.0 ) ( 30.7,8.4)@xmath178 ( 30.0,7.0)(0,1)1.0 ( 31.0,7.0)(0,1)1.0 ( 34.1,10)(0,1)1 ( 33.3,10)(1,0)1.5 ( 34.1,10 ) ( 34.3,10.4)@xmath178 ( 34.8,7)(0,1)3 ( 33.3,9.0)(0,1)1.0 ( 32.8,9.0)(1,0)1 ( 33.3,9.0 ) ( 33.5,9.4)@xmath178 ( 32.8,8.0)(0,1)1.0 ( 33.8,7.0)(0,1)2.0 ( 32.3,8.0)(1,0)1 ( 32.8,8.0 ) ( 33.0,8.4)@xmath178 ( 32.3,7.0)(0,1)1.0 ( 33.3,7.0)(0,1)1.0 ( 25.5,6.0)@xmath245 ( 24.7,4.5)period 4 cycles    a connection of the lc model can also be made with a ginzburg - landau theory of phase transitions and a feynman - wilson gas . these connections are discussed in ref .
@xcite . where @xmath1 and @xmath2 are related to coefficients in the ginzburg - landau approach .
once we identify an appropriate @xmath24 for a physical system , then we may use our general model to study the statistical behavior of the system .
the various models used in pion distribution can be related to our general model with @xmath40 by choosing @xmath24 as a term in a gauss hypergeometric series @xmath246 ; @xmath247_m [ b]_m}{[c]_m } \frac{z^m}{m ! }
\label{hyperf }      \\   { [ a]_m } = \frac{\gamma(a+m)}{\gamma(a ) }     & = & \frac{\gamma(a+n)}{\gamma(a ) } \frac{\gamma(a+n + m - n)}{\gamma(a+n ) }             \label{amfc }     = [ a]_n [ a+n]_{m - n }        \end{aligned}\ ] ] usefull values of @xmath248_n$ ] are @xmath249_n = ( 2n)!/(n ! 2^{2n})$ ] , @xmath250_n = n!$ ] , @xmath251_n = ( n+1)!$ ] . considering only positive @xmath5 , we choose @xmath252_{k-1 } [ b]_{k-1}}{[c]_{k-1 } } \frac{z^k}{(k-1 ) ! }        = x \frac{\gamma(a+k-1)}{\gamma(a ) } \frac{\gamma(c)}{\gamma(c+k-1 ) }           \frac{\gamma(b+k-1)}{\gamma(b ) } \frac{z^k}{(k-1)!}\end{aligned}\ ] ] for this case , the thermodynamic grand potential or the generating function is @xmath253 if we allow jets without a pion , then we may allow @xmath222 also .
for such a case , @xmath252_{k } [ b]_{k}}{[c]_{k } } \frac{z^k}{k ! }        = x \frac{\gamma(a+k)}{\gamma(a ) } \frac{\gamma(c)}{\gamma(c+k ) }           \frac{\gamma(b+k)}{\gamma(b ) } \frac{z^k}{k ! }    \\   f_0(\vec x ) & = & f_0(x , z ) = \log z(x , z )     = - \omega(x , z ) = \sum_{k=0}^\infty x_k        = x f(a , b ; c ; z)\end{aligned}\ ] ] we can see the only difference of the generating functions between the above two cases is the extra factor @xmath2 for the grand potential . we will mostly concentrate on the first case , i.e. , @xmath254 .    using eq .
( [ factfm ] ) or the recurrence relation eq.([recfm ] ) and eqs.([hyperf ] )  ( [ f0hg ] ) , @xmath255_{m } [ b]_{m}}{[c]_{m } } z^{m+1 }              f(a+m , b+m ; c+m ; z )                  \nonumber   \\    & &       + x m \frac{[a]_{m-1 } [ b]_{m-1}}{[c]_{m-1 } } z^{m }              f(a+m-1 , b+m-1 ; c+m-1 ; z)\end{aligned}\ ] ] the second order normalized factorial cumulant @xmath167 and the void variable @xmath160 are then @xmath256_2 [ b]_2}{[c]_2 } z^2 f(a+2 , b+2 ; c+2 ; z )            + 2 \frac{a b}{c } z f(a+1 , b+1 ; c+1 ; z ) }          { \frac{a b}{c } z f(a+1 , b+1 ; c+1 ; z ) + f(a , b ; c ; z ) }      \\
\chi(x , z ) & = & \frac{f_0(x , z)}{\bra a\ket }      = \frac{f_0(x , z)}{f_1(x , z ) }      \nonumber   \\     & = & \frac{f(a , b ; c ; z ) }           { \frac{a b}{c } z f(a+1 , b+1 ; c+1 ; z ) + f(a , b ; c ; z)}\end{aligned}\ ] ]    for some values of @xmath0 , @xmath257 , and @xmath258 , the hypergeometric function become a simple function ; @xmath259 various models of pion distributions can be related with these functions as listed in table [ tabl1 ] .
.various models with specific choice of @xmath40 and @xmath24 in hypergeometric series @xmath246 of eq.([hyperf ] ) . here
@xmath260 is not included , and thus @xmath261 . here
@xmath262 with @xmath105 except for poisson which has a finite @xmath263 .
fisher exponent @xmath3 for each @xmath24 as discussed in sect.[sectfrag ] are given too . [
cols="^,^,^,^,^,^,^",options="header " , ]     because of a unique role played by the poisson distribution and the form @xmath264 $ ] of the generating function , the cluster distribution @xmath265 is usually taken to be a poisson .
however , as noted before , other divisions are possible . using the same approach used above for @xmath266
, we can expand any distribution using a nb instead of poisson , i.e. , with @xmath267 using the form of the generating function @xmath268 for the nb . replacing @xmath2 by @xmath269 , the normalization factor of a new distribution
, we have @xmath270^{-x }       = \left[1 - { \cal n}(z ) { \cal g}(u=0)\right]^{-x}\end{aligned}\ ] ] the @xmath271 is @xmath272z)}{{\cal n}(z ) }      = \sum_k ( 1-u)^k { \cal p}_k\end{aligned}\ ] ] while the @xmath273 is @xmath274z)}{z(x , z ) }      = \frac{\left[1 - { \cal
n } { \cal g}(u)\right]^{-x } }             { \left[1 - { \cal n}\right]^{-x}}\end{aligned}\ ] ] thus we have @xmath275^m      = \sum_{m=0}^\infty p_m^{\rm nb}(x , { \cal n } ) [ { \cal g}(u)]^m           \nonumber   \\    & = & \sum_{m=0}^\infty p_m^{\rm nb }                \left[\sum_{j=0}^\infty ( 1-u)^j { \cal p}_j\right]^m      = \sum_{m=0}^\infty p_m^{\rm nb }          \sum_{\{n_k\}_m }   \prod_{k } \left[(1-u)^k{\cal p}_k\right]^{n_k }      \\
p_n(\vec x ) & = & \sum_m p_m^{\rm nb}(x , { \cal n } )           \sum_{\{n_k\}_m } \prod_k { \cal p}_k^{n_k}(\vec x )    \label{compnb}\end{aligned}\ ] ] the result of eq.([compnb ] ) shows that the distribution @xmath145 can be written as a compound probability distribution of a negative binomial @xmath276 with another probability @xmath277 distribution generated from @xmath271 . for the case of @xmath278 , which may be considered as the fugacity @xmath279 for a particle with chemical potential @xmath280
, @xmath281 can further be decomposed as @xmath282 i.e. , @xmath277 for this case is poisson @xmath283 .
if @xmath284 , then @xmath285 without @xmath222 and @xmath286^x = e^{f_0}$ ] . for @xmath214 given in table [ tabl1 ] , the @xmath277 becomes the same probability with @xmath24 replaced by @xmath287 .
as an example the hga with @xmath288\right)$ ] can be decomposed as a sequential process consisting of a nb distribution of clusters with @xmath289 which is then followed by a breakup of clusters distributed with a hga distribution with @xmath290 but without voids , i.e. , with @xmath291 .
thus this decomposition separates the parameter @xmath1 assigned to cluster from other parameters .
event - by - event studies from high energy collisions are being used to study the details of particle multiplicity distributions as , for example , those associated with pions .
such studies not only give information about the mean number of particles produced , but also information about fluctuations and higher order moments of the probability distribution which are important tools for studying the underlying processes and mechanisms that operate .
they are also useful in distinguishing various phenomenological models .
issues associated with fluctuations play an important role in many areas of physics and departures from poisson statistics are of current interest .
one purpose of this paper was an investigation of various models of particle multiplicity distributions that can be used in event - by - event analysis .
these various phenomenological models are developed using a general form of a unified model which is based on a grand canonical partition function and an underlying weight arising from feynman s path integral approach to statistical processes .
a resulting distribution has three control parameters called @xmath0 , @xmath1 , @xmath2 .
the relationships of these parameters to various physical quantities are discussed .
one important result is the connection of the parameter @xmath0 to the fisher exponent @xmath3 ; namely @xmath292 for hga .
this connection arises from a parallel we developed between the model for particle multiplicity distributions considered here and our previous approach to cluster yields .
since an exact description of particle multiplicity distributions is not known , we have considered several cases with different @xmath3 s or @xmath0 s which are contained in our unified description . moreover , many of the existing distributions currently used in particle phenomenology are shown to be special choices of @xmath3 or @xmath0 which appear in a generalized hypergeometric model called hga .
these include the poisson distribution coming from coherent emission , chaotic emission producing a negative binomial distribution , combinations of coherent and chaotic processes leading to signal / noise distributions and field emission from lorentzian line shapes producing the lorentz / catalan distribution . using the hga model combinants , cumulants and moments
are discussed and a physical significance is given to combinants in terms of the underlying partition weights of a feynman path integral approach to statistical processes .
the parameter @xmath0 or fisher exponent @xmath3 is shown to play an important role in the behavior of the combinants which manifest itself in various physical relationhips . the hga model and its associated special cases
are used to explore a wide variety of phenomena .
these include : linked pair approximations leading to hierarchical scaling relations on the reduced cumulant level , generalized void scaling relations , clan variable descriptions and their connections with stochastic variables and branching processes , kno scaling behavior , enhanced non - poissonian fluctuations .
models based on an underlying random walk description are also discussed .    in this paper
we compared various particle multiplicity distributions within the hypergeometric model hga .
our results show that even though various distributions have the same mean and fluctuation , the distribution itself or the underlying mechanism could be different .
comparisons within the hga model also show that just comparing void variables @xmath160 and @xmath162 or mean @xmath293 and fluctuation @xmath294 or @xmath295 is not enough to distinguish different models that descibe particle multiplicity data .
thus , to find the correct distribution and underlying mechanism from various data more information than just the mean and fluctuation are necessary and new variables should be found which are quite different between different models .
for example , higher order reduced factorial cumulants need also to be evaluated such as the third order cumulant @xmath296 .
applications of our approach to the charged multiplicity data of l3 and h1 are given .
the mean @xmath297 , fluctuation variable @xmath162 , void scaling variable @xmath160 , and third order reduced cumulant variable @xmath296 obtained from these experimnets are compared with various models discussed in this paper .
the value of the third cumulant @xmath296 shows that the grw1d model can fit the charge multiplicity data of l3 and h1 .
the hga model also gives a fit to the l3 data , but is not as successful as grw1d with respect to the third order reduced cumulant of the h1 data .
both have very similar multiplicity distributions as shown in fig.[jetdat ] which compares the two models with the charged particle multiplicity distributions of jets for l3 and h1 data .
differences appear in the high multiplicity events with reduced multiplicity @xmath298 .
the reduced third cumulant is negative for the h1 data of table [ tabldat ] .
negative values of a cumulant require negative values of the parameter @xmath0 which are allowed in grw1d if the parameters @xmath1 and @xmath2 are also taken as negative , but a negative value of @xmath0 is not acceptable in hga since this would require a value of @xmath2 greater than 1 and thus a complex @xmath1 .
we also discuss the kno scaling behavior for these data and conclude that the kno scaling behavior would not be clearly exhibited in the jets for l3 and h1 .
in this paper we have also generalized the compound distribution that arises from sequential process which may reveal the dynamical structure of the distribution .
specifically , the underlying sequential picture involves a two step process where the final distribution arises from the production of clusters followed by a subsequent decay of the clusters . for the hga model
, the final distribution is obtained from compounding a poisson distribution of clusters with a nb distribution coming from the decay of each of the clusters .
the hga may arise through a three step sequential process of poisson - poisson - logarithmic compound distribution .
it is also shown that the hga can arise from a two step sequential process of a nb distribution followed by a new hga with a different mean value .
k. rajagopal and f. wilczek , nucl .
* b404 * , 557 ( 1993 ) . s. gavin and b. mller , phys . lett . * b329 * , 486 ( 1994 ) .
g. bayer and h. heiselberg , phys
b469 * , 7 ( 1999 ) . v. heinz and b.v .
jacak , ann .
nucl . part .
, * 49 * , ( 1999 ) . s. pratt , phys . lett . *
b301 * , 159 ( 1993 ) .
csrg and j. zimanyi , phys .
lett . * 80 * , 916 ( 1998 ) .
andrews , c.g .
townsend , h.j .
meisner , d.s .
durfee , d.m .
korn , w. ketterle , science * 275 * , 637 ( 1997 ) .
a. bialas and r. peschanski , nucl . phys . *
b272 * , 703 ( 1986 ) ; * b306 * , 857 ( 1988 ) .
p. carruthers and c.c .
shih , phys .
* 127b * , 242 ( 1983 ) . l. van hove , phys . lett . *
b232 * , 509 ( 1989 ) .
l. van hove and a. giovannini , z. phys .
* c30 * , 391 ( 1988 ) .
z. koba , h.b .
nielson and p. olessen , nucl .
* b40 * , 317 ( 1972 ) .
p. carruthers and i. sarcevic , phys .
b189 * , 442 ( 1987 ) .
m. gadzicki and s. mrwczynski , z. phys .
* c54 * , 127 ( 1992 ) .
m. stephanov , k. rajagopal , e. shuryak , phys .
lett . * 81 * , 4816 ( 1998 ) .
haiasz , a.d .
jackson , r.e .
shrock , m.a .
stepanov , and j.j.m .
verbaarschot , phys . rev . *
d58 * , 96007 ( 1998 ) .
g. baym and h. heiselberg , phys .
b469 * , 7 ( 1999 ) . s. jeon and v. koch , phys . rev .
lett * 83 * , 5435 ( 1999 ) ; phys .
* 85 * , 2076 ( 2000 ) .
m. asakawa , u. heinz and b. mller , phys .
lett . * 85 * , 2072 ( 2000 ) .
h. heiselberg and a.d .
jackson , phys .
rev . * c63 * , 064904 ( 2001 ) .
mekjian and s.j .
lee , phys .
rev . * a44 * , 6294 ( 1991 ) .
s.j . lee and a.z .
mekjian , phys .
rev . * c45 * , 1284 ( 1992 ) .
s.j . lee and a.z .
mekjian , phys . rev . *
c45 * , 365 ( 1992 ) .
s.j . lee and a.z .
mekjian , phys . rev .
* c56 * , 2621 ( 1997 ) .
s.j . lee and a.z .
mekjian , phys . rev .
* c47 * , 2266 ( 1993 ) .
s.j . lee and a.z .
mekjian , phys . rev . *
c50 * , 3025 ( 1994 ) .
k.c . chase and a.z .
mekjian , phys . rev . * c 49 * , 2164 ( 1994 ) ; * 50 * , 2078 ( 1994 ) . s. lupia , a. giovannini , and r. ugoccioni , z. phys .
* c59 * , 427 ( 1993 ) . m. gyulassy and s. kaufmann , phys .
* 40 * , 298 ( 1978 ) ; j. phys .
* a11 * , 1715 ( 1978 ) .
mekjian , b.r .
schlei and d. strottman , phys . rev .
* c58 * , 3627 ( 1998 ) .
mekjian , phys .
lett . * 86 * , 220 ( 2001 ) .
a. mekjian , phys . rev . *
c65 * , 014907 ( 2002 ) .
p. carruthers and c.c .
shih , int . j. modern phys .
* a2 * , 1447 ( 1987 ) .
f. becattini , a. giovannini and s. lupia , z. phys . *
c72 * , 43 ( 1996 ) .
r. hwa , in _ proc . of the santa fe workshop on inbtermittency in high energy collisions _
( world scientific , singapore , 1990 , f. cooper , r. hwa and i. sarcevic , edit ) .
s. hegyi , phys . lett . *
b309 * , 443 ( 1993 ) ; * b318 * , 642 ( 1993 ) .
feynman , _ statistical physics _
( benjamin / cummings , new york , 1972 ) .
j. klauder and e.c.g .
sudarshan , _ fundamentals of quantum optics _ ( benjamin , n.y . , 1968 ) .
h. biyajima , prog .
physics * 69 * , 966 ( 1983 ) .
l3 collaboration , lanl preprint hep - ex/0110072 ( 2001 ) .
s. aid et al .
, h1 collaboration , z. phys . * c72 * , 573 ( 1996 ) ; lanl hep - ex/9608011 ( 1996 ) .