Monte-Carlo approach to particle-field interactions and the kinetics of the chiral phase transition

Monte-Carlo approach to particle-field interactions and the kinetics of the chiral phase transition


The kinetics of the chiral phase transition is studied within a linear quark-meson- model, using a Monte-Carlo approach to semiclassical particle-field dynamics. The meson fields are described on the mean-field level and quarks and antiquarks as ensembles of test particles. Collisions between quarks and antiquarks as well as the annihilation to mesons and the decay of mesons is treated, using the corresponding transition-matrix elements from the underlying quantum field theory, obeying strictly the rule of detailed balance and energy-momentum conservation. The approach allows to study fluctuations without making ad hoc assumptions concerning the statistical nature of the random process as necessary in Langevin-Fokker-Planck frameworks.

1 Introduction

One of the motivations for the study of ultrarelativistic heavy-ion collisions is to gain a detailed understanding of the phase diagram of strongly interacting matter[1]. At the largest energies as achieved at the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) a hot and dense fireball is formed which can be described to a surprising accuracy as a nearly perfect fluid of strongly coupled quarks and gluons (QGP) undergoing a transition to a hot hadron-resonance gas. In these situations, where the net-baryon density or the baryon-chemical potential are small, lattice-QCD (lQCD) calculations indicate a crossover transition from confined to deconfined matter as well as from a phase where chiral symmetry is spontaneously broken to one where it is restored at a (common) transition temperature [2].

At lower collision energies, as studied in the RHIC beam-energy scan (BES) program, at the CERN SPS, and the future FAIR and NICA experiments, the produced medium starts at lower temperatures and larger net-baryon densities. Since in this situation the application of lQCD is challenging due to the “sign problem” at finite , one relies on effective chiral models, which predict the existence of a first-order chiral-phase-transition line ending in a critical point of a second-order phase transition.

For theory the challenge is to provide possible observables for this phase structure in heavy-ion collisions, like the (“grand-canonical”) fluctuations of conserved charges like net-baryon number or electric charge. Not only the question, how to effectively model the phase transition (e.g., with the Nambu-Jona-Lasinio (NJL) or the (linear) model with extensions taking into account gluonic degrees of freedom implementing Polyakov loops) arises but also, which of the features of the phase structure predicted for such models applying thermal quantum field theory (describing a medium in thermal and chemical equilibrium) like (critical) fluctuations of conserved charges survive for a rapidly expanding and cooling fireball as created in heavy-ion collisions.

To address the latter question, one relies on transport simulations to describe the off-equilibrium dynamics of the fireball. One approach is the use of ideal or viscous hydrodynamics to describe the bulk evolution of the fireball (assuming a state close to local thermal equilibrium), which successfully describes key phenomena of heavy-ion collisions, and adding the fluctuations by hand in a Langevin approach[3, 4, 5]. On the other hand this implies that the statistics of the random process has to be put in as an ad hoc assumption. Usually a Gaussian Markovian (“white noise”) is assumed, but the simulation of non-Markovian (“colored noise”) processes is feasible in principle [6].

On the other hand, one would like to study the implication of different transition scenarios (cross-over, 1 order, 2 order) as realized in the various quantum-field theoretical models on the nature of the hopefully observable fluctuations [7, 8, 9, 10, 11, 12, 13] on the statistics of the probably observable fluctuations.

In this work we present a novel Monte-Carlo approach to address this challenging problem using the most simple quark-meson linear model [14]. The meson fields are treated on the mean-field level and the quarks and antiquarks are realized in terms of a test-particle ensemble. Here the challenge is to implement “discrete” local interaction processes like elastic collisions and reactions like the annihilation to a meson and the decay of mesons to a pair admitting not only kinetic but also chemical equilibration starting from an off-equilibrium situation, using the transition-probability matrix elements of the underlying quantum field theory.

2 Linear quark-meson model

To investigate the feasibility of a kinetic description of the off-equilibrium dynamics of the chiral phase transition the most simple two-flavor chiral model, based on the chiral group is considered, using a flavor doublet of Dirac fields , describing and quarks and antiquarks and a four-dimensional real-valued set of scalar fields transforming under the SO(4) representation of the chiral group [15]. The Lagrangian reads


with the meson potential


where denotes the Yukawa coupling between quarks and mesons, the meson coupling constant (corresponding to a mass of ), the pion-decay constant, and . The potential (2) contains the explicit breaking of the chiral symmetry due to the finite current quark masses resulting in a non-zero pion mass, , of the pseudo-Goldstone modes . The constituent quark masses are given by .

The grand-canonical potential in mean-field approximation reads




where and the quark-degeneracy factor . The mean fields have to be evaluated self-consistently from the equilibrium condition


In the following we restrict ourselves to vanishing pion mean fields. A nice feature of this model is that by varying the Yukawa coupling , one finds different kinds of phase transition as illustrated in Fig. 1.

Figure 1: The phase diagram for the linear model in mean-field approximation at , . Left: the order parameter and the effective mass .

3 Semiclassical particle-field dynamics

The challenge in applying the above model to an off-equilibrium dynamical simulation of a system of particles (here quarks and antiquarks) and mean fields (representing the mesons) is that in order to reproduce the equilibrium-phase structure as depicted in 1 as the stationary limit, one has to ensure that both kinetic and chemical equilibration is possible through the introduction of the appropriate elastic collision terms for and scattering as well as quark-number changing processes such as . In a full kinetic approach this is achieved by a set of coupled Boltzmann-Vlasov equations, which read in our case schematically (again restricting ourselves to the case of vanishing pion-mean fields)


Here, and denote collision integrals contributing to the meson-mean-field and quark-phase-space distribution functions respectively.

In the following a novel scheme to Monte-Carlo simulate such a system of kinetic equations is defined, where one describes the mesons solely with a mean field and the quarks and antiquarks in terms of test particles. While the elastic-collision term is realized in a straight-forward way using the corresponding cross section from the underlying linear model, one has to find a way to realize the interactions in such a scheme, while still fulfilling energy-momentum conservation and the principle of detailed balance, which are the fundamental principles constraining the off-equilibrium dynamics and ensuring the proper (Maxwell-Boltzmann) equilibrium limit.

In our recently developed model (Dynamical Simulation of a Linear Sigma Model, DSLAM) this challenge is solved as follows: In order to properly simulate the collision terms on the right-hand sides of Eqs. (6) and (7) we define a space-time grid. In each time step at each spatial cell the cross section for the annihilation process is used to stochastically determine an energy-momentum transfer from the initial pair, located in the cell. This energy-momentum change is transferred to the mean field in terms of an appropriate relativistic Gaussian wave packet (which simulates the gain term in ), and the pair is taken out of the test-particle ensemble (which simulates the corresponding loss term in ).

Figure 2: While the energies of the quarks and the meson field in the thermal-box simulation show anti-correlated thermal fluctuations of the order for the quarks and for the field the numerical fluctuations of the total energy amount to only .
Figure 3: The total quark number in the thermal-box simulation shows fluctuations due to the dynamical pair-creation and annihilation processes, leading to an exchange of energy between the quarks and the mean field as shown in the left Fig. 3.

To simulate also the appropriate decay process, , we have to “particlize” the mean field locally in each spatial cell. This is done in the spirit of a course-graining procedure: First, the total energy-momentum content of the field within the cell in terms of the corresponding -field energy-momentum tensor is determined. Then one assumes a local thermal equilibrium phase-space distribution, equivalent to this energy-momentum tensor. In order to fulfill detailed balance, the temperature has to be the same as that for the corresponding procedure for the quarks and antiquarks. The temperature is related to the mean-field value which depends on the scalar quark-antiquark density. In this way a temperature can be determined. It is important to note that it is defined in the local rest frame of the heat bath and thus the -phase-space distribution is given by a Maxwell-Jüttner distribution , where is the four-velocity, given by the total field-four-momentum in the spacial cell under consideration, . Now in each time-step within each spatial cell one can choose an ensemble of particles according to this local Maxwell-Jüttner distribution and using the corresponding decay rate to determine the gain term to . The loss term for the mean-field equation in the collision term is again achieved by taking the appropriate amount of energy and momentum out of the mean field in terms of a Gaussian wave packet.

In summary we have achieved a scheme which enables us to simulate the set of Boltzmann-Vlasov Eqs. (6) and (7) using test particles for the quarks and antiquarks and restricting the description of the mesons strictly to the mean-field level. The scheme by construction fulfills energy-momentum conservation through the Gaussian wave-packet description for the exchange of energy and momentum between the mean field and the test particles. At the same time also the principle of detailed balance is fulfilled, using the coarse-graining approach to locally map the field-energy-momentum distribution to a local-equilibrium Maxwell-Jüttner distribution to reinterpret the mean field as a meson phase-space distribution and using the leading-order transition rates (cross sections) of the underlying QFT linear model fulfilling the detailed-balance principle.

4 Proof of principle: “Box” calculations

As a plausibility check of the simulation method the stability of an equilibrium situation has been tested in a finite cubic box with periodic boundary conditions[14]. The stability of the energy conservation is demonstrated in Fig. 3. While the kinetic energy of the quarks and the energy of the mean field show anti-correlated thermal fluctuations the total energy stays stable within the numerical accuracy of the simulation.

The annihilation and creation processes of quark-antiquark pairs lead to thermal fluctuations in the total particle number, as shown in Fig. 3.

Figure 4: The distribution of the quarks’ kinetic energy in a thermal-box calculation show an excellent agreement with the expected Maxwell-Boltzmann distribution.
Figure 5: The distribution of the kinetic field energy to its Fourier modes. The classical “UV catastrophe”, i.e., the equipartition theorem associating an average energy of for each Fourier mode, is avoided by the use of a finite finite width of the Gaussian wave packets for the transfer of energy and momentum between particles and field.

The energy distribution of the quarks shows the expected Maxwell-Boltzmann distribution at the expected temperature, as demonstrated in Fig. 5. Also the spectral analysis of the mean field shows that the thermal equipartition theorem for the kinetic field energy is fulfilled at small wave numbers (long wavelengths), i.e., each mode contains an average energy of . On the other hand the “UV catastrophe” must be avoided due to the finite total-energy content within a finite box. Indeed, the energy-momentum transfer between particles and the field due to the pair-annihilation and -creation processes is not strictly “local” but occur within a finite volume whose scale is fixed by the finite width, , of the Gaussian wave packets used as field increments to keep care of the correct energy-momentum transfer in each process. Thus the large wave numbers (short wave lengths) are effectively cut off at a scale , as shown by plotting the corresponding Gaussian on top of the -distribution of the kinetic field energy (Fig. 5).

Figure 6: Starting with a uniform mean field, a short time after the start of the simulation the field is disturbed by some Gaussian wave packets due to quark-pair annihilation processes (left panel). In the long-time limit the full equilibrium-thermal-field fluctuations have developed.
Figure 7: The field fluctuations are Gaussian distributed around an average value which can slowly drift with time due to the dynamical fluctuations of the field momentum.

The dynamical behavior of the field in our thermal-box simulation is illustrated in Fig. 7, showing the field distribution within the plane of the simulation. Starting the simulation with a uniform mean-field equilibrium value, after a short time some local blob-like disturbances have developed (left panel) due to quark-pair-annihilation processes, leading to the propagation of Gaussian wave packets on top of the still quite uniform average field value, as implemented by our concept of describing the exchange between field and particles in these processes. The right panel shows the fully developed equilibrium-thermal-field fluctuations in the long-time limit of the simulation. In this limit the field values are Gaussian distributed around a mean value (cf. Fig. 7), which is expected due to many random energy-momentum transfers between the field and particles. Thereby the average field value can slowly drift with time due to the thermal fluctuations of the field energy and momentum.

5 Conclusions and outlook

In this talk it was demonstrated that the dynamical description of the chiral phase transition in a simple quark-meson model is feasible with a novel Monte-Carlo-simulation technique for a corresponding coupled Boltzmann-transport equation for the meson mean-field and the quark and anti-quark phase-space distribution functions, implementing both elastic quark and anti-quark scattering as well as quark-antiquark-pair creation and annihilation processes, enabling both kinetic and chemical equilibration between particles and the mean field (mesons).

The simulation is set up in a way that both energy-momentum conservation and the principle of detailed balance are precisely realized (within the limits of achievable numerical accuracy). While the elastic quark-scattering processes are simulated with a straight-forward test-particle realization, the particle-field kinetics has demanded the development of a novel scheme.

The annihilation process is evaluated by the Monte-Carlo sampling according to the corresponding transition matrix elements from the underlying quantum-field theoretical interpretation of the linear model. The corresponding energy and momentum are precisely transferred to the -mean field in terms of an appropriate disturbance in form of a relativistic Gaussian wave packet. To obey the principle of detailed balance, also the inverse decay has to be simulated. To that purpose a coarse-graining procedure based on the local energy-momentum content of the field has been used for a “particlization” of the mean field in terms of a local Boltzmann-Jüttner equilibrium distribution, which in turn enables a Monte-Carlo sampling of the -meson decay according to the decay rate from the quantum-field theory picture of the model, which automatically implements detailed balance.

This scheme has now been applied to off-equilibrium situations as in a “thermal quench” in a box where the particles and fields are initialized with different temperatures and it was demonstrated that the distribution can describe a possible phase transition from the initial state to a late-time (equilibrium) state.

Last but not least also the case of “expanding fireballs”, mimicking the situation of the medium created in heavy-ion collisions, is under study.


This work was partially supported by the Bundesministerium für Bildung und Forschung (BMBF Förderkennzeichen 05P12RFFTS) and by the Helmholtz International Center for FAIR (HIC for FAIR) within the framework of the LOEWE program (Landesoffensive zur Entwicklung Wissenschaftlich-Ökonomischer Exzellenz) launched by the State of Hesse. C. W. and A. M. acknowledge support by the Helmholtz Graduate School for Hadron and Ion Research (HGS-HIRe), and the Helmholtz Research School for Quark Matter Studies in Heavy Ion Collisions (HQM). Numerical computations have been performed at the Center for Scientific Computing (CSC). H. v. H. has been supported by the Deutsche Forschungsgemeinschaft (DFG) under grant number GR 1536/8-1. C.W. has been supported by BMBF under grant number 0512RFFTS.



  1. Friman B, Hohne C, Knoll J, Leupold S, Randrup J et al. 2011 Lect.Notes Phys. 814 pp. 980 URL
  2. Philipsen O 2013 Prog. Part. Nucl. Phys. 70 55–107 URL
  3. Nahrgang M, Leupold S, Herold C and Bleicher M 2011 Phys. Rev. C 84 024912 URL
  4. Nahrgang M, Leupold S and Bleicher M 2012 Phys. Lett. B 711 109–116 URL
  5. Herold C, Nahrgang M, Mishustin I and Bleicher M 2013 Phys. Rev. C 87 014907 URL
  6. Schmidt J, Meistrenko A, van Hees H, Xu Z and Greiner C 2015 Phys. Rev. E 91 032125 URL
  7. Stephanov M A, Rajagopal K and Shuryak E V 1999 Phys. Rev. D 60 114028 URL
  8. Schäfer B J, Pawlowski J M and Wambach J 2007 Phys. Rev. D 76 074023 URL
  9. Skokov V, Friman B and Redlich K 2011 Phys. Rev. C 83 054904 URL
  10. Skokov V, Stokic B, Friman B and Redlich K 2010 Phys. Rev. C 82 015206 URL
  11. Braun-Munzinger P, Friman B, Karsch F, Redlich K and Skokov V 2012 Nucl. Phys. A 880 48–64 URL
  12. Schaefer B J 2012 Phys. Atom. Nucl. 75 741–743
  13. Morita K, Skokov V, Friman B and Redlich K 2014 Eur. Phys. J. C 74 2706 URL
  14. Wesp C, van Hees H, Meistrenko A and Greiner C 2015 Phys. Rev. E 91 043302 URL
  15. Gell-Mann M and Lévy M 1960 Nuovo Cim. 16 705 URL
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description