The thermodynamic model for nuclear multifragmentation

A great many observables seen in intermediate energy heavy ion collisions can be explained on the basis of statistical equilibrium. Calculations based on statistical equilibrium can be implemented in microcanonical ensemble (energy and number of particles in the system are kept fixed), canonical ensemble (temperature and number of particles are kept fixed) or grand canonical ensemble (fixed temperature and a variable number of particles but with an assigned average). This paper deals with calculations with canonical ensembles. A recursive relation developed recently allows calculations with arbitrary precision for many nuclear problems. Calculations are done to study the nature of phase transition in intermediate energy heavy ion collision, to study the caloric curves for nuclei and to explore the possibility of negative specific heat because of the finiteness of nuclear systems. The model can also be used for detailed calculations of other observables not connected with phase transitions, such as populations of selected isotopes in a heavy ion collision.The model also serves a pedagogical purpose. For the problems at hand, both the canonical and grand canonical solutions are obtainable with arbitrary accuracy hence we can compare the values of observables obtained from the canonical calculations with those from the grand canonical. Sometimes, very interesting discrepancies are found.To illustrate the predictive power of the model, calculated observables are compared with data from the central collisions of Sn isotopes.     

Baryon number projection for an interacting hadron gas

Calculations of hadronic matter usually enforce conservation of the average baryon number density using the grand canonical ensemble. We have performed calculations for an interacting system in the canonical ensemble with fixed baryon numberN _b , as appropriate for a finite fireball of the type produced in ultra relativistic heavy ion collisions. These results are compared with those obtained from calculations in the grand canonical ensemble. For an interacting nucleon gas the two ensembles yield free energies which differ by approximately 5%.     

A new simulation method for the determination of phase equilibria in mixtures in the grand canonical ensemble

A grand canonical Monte Carlo (GCMC) simulation method is presented for the determination of the phase equilibria of mixtures. The coexistence is derived by expanding the pressure into a Taylor series as a function of the temperature and the chemical potentials that are the independent intensive variables of the grand canonical ensemble. The coefficients of the Taylor series can be calculated from ensemble averages and fluctuation formulae that are obtained from GCMC simulations in both phases. The method is able to produce the equilibrium data in a certain domain of the (T, p) plane from two GCMC simulations. The vapour-liquid equilibrium results obtained for a Lennard-Jones mixture agree well with the corresponding Gibbs ensemble Monte Carlo data.     

On the meaning of the canonical ensemble

Thermal equilibrium between (quantum) systems is taken to mean stability for the combined system. Necessary and sufficient conditions for such stability are found and used to show that any system in equilibrium with suitably complex second system (heat bath) will be characterized by a canonical ensemble. Thus the notion of temperature is derived directly from that of equilibrium, without, for example, recourse to microcanonical ensembles or information theory. Discussed briefly are the generalization of these results to grand canonical ensembles and their application to the equilibrium between a black hole and the surrounding radiation field.     

The grand canonical ensemble Monte Carlo method applied to electrolyte solutions

A grand canonical ensemble Monte Carlo method is developed for application to electrolyte solutions.

The method proves to be of comparable accuracy and speed to the conventional NVT Monte Carlo method for electrolytes but has the added advantage of being able to fix the chemical potential. This latter point is vital for the study of surface phenomena. Application of the method to 1:1 and 2:2 primitive model electrolytes is made as is a comparison with the results of approximate statistical mechanical treatments.     

Linear response,persistent currents and related problems

Based on a general linear response approach we provide a systematic and unified survey of existing theories on persistent currents. The central notions in this context are equilibrium and dynamic persistent currents which are analyzed with respect to their similarities and differences in the canonical and grand canonical ensemble. We present criteria which relate the existence of persistent currents to the equipartition law and ergodicity for current correlators. We find that in additive Fermion systems at low temperatures both kinds of persistent currents coincide in the canonical ensemble whereas they differ in the grand canonical ensemble. Comparing different works on averaged persistent currents in diffusive mesoscopic rings within our framework and discussing several methods of calculating canonical currents with the help of grand canonical ensembles, we clarify some misunderstandings which have arisen in methodologically different approaches to the phenomenon of persistent currents. Finally, we relate the presence of dynamic persistent currents to the Hall conductivity on a finite cylinder and the center coordinate Kubo formula for the Hall conductivity.     

Connection between the local Maxwell and local grand canonical distributions

The one-particle distribution corresponding to the local grand canonical ensemble is calculated rigorously. It is shown to coincide with the local Maxwell distribution provided the macroscopic parameters characterizing the ensembles are chosen properly. Their physical meaning is discussed.     

Novel approach to deriving the canonical generating functional in lattice QCD at a finite chemical potential

A novel approach to the problem of deriving the generating functional for the canonical ensemble in lattice QCD at a nonzero chemical potential is proposed. The derivation proceeds in several steps. First, the baryon density for imaginary values of the chemical potential is obtained. Then, again for imaginary values of the chemical potential, the generating functional of the grand canonical ensemble is derived. In this analysis, a fit of baryon density is employed toward simplifying the procedure of numerical integration. Finally, the generating potential for the canonical ensemble is derived using a high-precision numerical Fourier transform. The generating functional for the canonical ensemble is also derived using the known hopping-parameter expansion, and the results obtained with the two methods are compared for the deconfinement phase in the lattice QCD with two flavors.     

Density-functional theory of inhomogeneous fluids in the canonical ensemble

We present a density-functional approach for dealing with inhomogeneous fluids in the canonical ensemble. A general relation is proposed between the free-energy functionals in the canonical and the grand canonical ensembles. The minimization of the canonical-ensemble free-energy functional gives rise to Euler-Lagrange equations which involve averaged Ornstein-Zernike equations of second and third order. The theory is especially appropriate for systems with a small, fixed number of particles. As an example of application we obtain accurate results for the density profile of a hard-sphere fluid in a closed spherical cavity that contains only a few particles.     

Relationship between statistical sums of canonical and grand canonical ensembles of interacting particles

A strict asymptotic relation between statistical sums of canonical and grand canonical ensembles is derived. Institute of Applied Physics at Irkutsk State University. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 6, pp. 43–46, June, 2000.     

The thermodynamic model for relativistic heavy ion collisions

The collision between two heavy ions at high energies is discussed within the framework of a statistical model which is developed step by step so that it can be followed by a reader unfamiliar with this approach. As we hope to show, the statistical model forms a comprehensive framework for investigating various aspects of such collisions. When statistical thermodynamics is coupled with a one and two fireball model or a firestreak model, single particle and multiparticle inclusive cross sections can be evaluated. Specifically, single particle proton, deuteron, triton and pion cross sections are calculated and compared with experiment. Two particle correlations are discussed using the microcanonical ensemble. Multiplicity distributions are evaluated using the canonical ensemble and composite particle formation is simply obtained from the grand canonical ensemble. The power law behaviour of composite particle cross sections is studied. The thermodynamic model is justified by calculating various reaction rates. Many comparisons with experimental data are made.     

A molecular dynamics method for simulations in the canonical ensemble

A molecular dynamics simulation method which can generate configurations belonging to the canonical (T, V, N) ensemble or the constant temperature constant pressure (T, P, N) ensemble, is proposed. The physical system of interest consists of N particles (f degrees of freedom), to which an external, macroscopic variable and its conjugate momentum are added. This device allows the total energy of the physical system to fluctuate. The equilibrium distribution of the energy coincides with the canonical distribution both in momentum and in coordinate space. The method is tested for an atomic fluid (Ar) and works well.     

The disassembly of nuclear matter

A statistical model is applied for multi-fragment final states in nuclear collisions with bombarding energies E/A ≈ 100 MeV. A portion of the intermediate system formed is assumed to decay according to the available classical non-relativistic phase space, calculated in a grand canonical ensemble. The model correlates and predicts many experimental observables in terms of three parameters: the available energy per nucleon, the isospin asymmetry, and the effective interaction volume.     

刚性多原子分子的正则系综分子动力学算法

通过将Nose所定义的扩展系统哈密顿量推广到刚性多原子分子系统中，严格地推导了刚性多原子分子在正则系综下的运动方程.在此基础上，证明了依靠所给出的运动方程能得到平衡态的正则分布，并给出了该方程所对应的可逆而守恒的积分格式. 关键词： 刚性多原子分子 正则系综 分子动力学算法     

Gibbs ensembles in relativistic classical and quantum mechanics

Classical and quantum Gibbs ensembles are constructed for equilibrium statistical mechanics in the framework of an extension to many-body theory of a relativistic mechanics proposed by Stueckelberg. In addition to the usual chemical potential in the grand canonical ensemble, there is a new potential corresponding to the mass degree of freedom of relativistic systems. It is shown that in the nonrelativistic limit the relativistic ensembles we have obtained reduce to the usual ones, and mass fluctuations for the free-particle gas approach the fluctuations in N. The ultrarelativistic limit of the canonical ensemble for the free-particle gas differs from the corresponding limit of the ensemble proposed by Jüttner and Pauli. Due to the mass degree of freedom, the quantum counting of states is different from that of the nonrelativistic theory. If the mass distribution is sufficiently sharp, the thermodynamical effects of this multiplicity will not be large. There may, however, be detectable effects such as a shift in the Fermi level and the critical temperature for Bose-Einstein condensation, and some change in specific heats.     

Improvement of the grand canonical quantum Monte Carlo method at low temperatures

We introduce an improvement of the algorithm for the simulation of correlated fermi systems in the grand canonical ensemble. Using this new method the computer time grows no more with the square but essentially linearly with the inverse temperature. At low temperatures the number of operations is diminished by a factor typically between 5 and 20. We present results for correlation functions of the one-dimensional Hubbard model at various band fillings.     

Linear scaling of computer time with the inverse temperature for the grand canonical quantum Monte Carlo method

Using a new method of storing and updating one particle Green's functions the computer time for the simulation of correlated fermi systems in the grand canonical ensemble scales linearly with the inverse temperature. This reduces the number of operations by a factor of typically 20 and more at low temperatures compared with the standard algorithm.     

Phase transitions in four-dimensional AdS black holes with a nonlinear electrodynamics source

In this work we consider black hole solutions to Einstein's theory coupled to a nonlinear power-law electromagnetic field with a fixed exponent value. We study the extended phase space thermodynamics in canonical and grand canonical ensembles, where the varying cosmological constant plays the role of an effective thermodynamic pressure. We examine thermodynamical phase transitions in such black holes and find that both first- and second-order phase transitions can occur in the canonical ensemble while, for the grand canonical ensemble, Hawking–Page and second-order phase transitions are allowed.     

Canonical hadron gas interpretation ofp-A E802 data

Hadron gas models have proved successful in predicting particle production in relativistic nucleus-nucleus collisions. The extension of these models to the smaller systems formed in proton-nucleus collisions requires that the finite size of the system be considered. We study two features introduced by the finite size: the need to conserve strangeness and baryon number exactly by performing calculations in the canonical ensemble, and the inclusion of a finite size geometrical correction term in the single particle density of states. We find significant differences between the grand canonical and canonical ensembles and a strong dependence on the baryon number of the system.