# Applications of Jarzynski’s relation in lattice gauge theories

###### Abstract

Preprint: CP3-Origins-2016-045 DNRF90

Jarzynski’s equality is a well-known result in statistical mechanics, relating free-energy differences between equilibrium ensembles with fluctuations in the work performed during non-equilibrium transformations from one ensemble to the other. In this work, an extension of this relation to lattice gauge theory will be presented, along with numerical results for the gauge model in three dimensions and for the equation of state in Yang-Mills theory in four dimensions. Then, further applications will be discussed, in particular for the Schrödinger functional and for the study of QCD in strong magnetic fields.

Applications of Jarzynski’s relation in lattice gauge theories

Alessandro Nada^{†}^{†}thanks: Speaker. , Michele Caselle, Gianluca Costagliola, Marco Panero and Arianna Toniato

Department of Physics, University of Turin & INFN, Turin

Via Pietro Giuria 1, I-10125 Turin, Italy

CP-Origins & Danish IAS, University of Southern Denmark

Campusvej 55, 5230 Odense M., Denmark

E-mail: anada@to.infn.it, caselle@to.infn.it, costagli@to.infn.it, marco.panero@unito.it, toniato@cp3.sdu.dk

\abstract@cs

## 1 Introduction

The computation of free energy differences in lattice gauge theories (LGTs) plays a crucial role in the study of a large set of physically interesting quantities and objects, such as ’t Hooft loops, interfaces between center domains at non-zero temperature and many others. The determination of the equation of state in QCD and QCD-like theories represents one of the most important examples: the pressure is indeed naturally related to the computation of the free energy density , through the relation (valid in the thermodynamic limit) .

The numerical determination of free-energy differences in LGTs is a nontrivial computational challenge that motivates the search for new algorithms that can be easily implemented to study a large and diverse set of physical quantities. Here we summarize our recent work [1], in which we presented a novel method for the computation of free energy differences based on the application to LGTs of Jarzynski’s relation, a very well known result in statistical mechanics that was proved almost 20 years ago by Jarzynski [2, 3]: it has been used since in a plethora of numerical studies in statistical mechanics and even verified experimentally (see, for instance, ref. [4]). We applied this method (briefly reviewed in section 2) for the calculation of two very different quantities: the free energy of an interface in lattice gauge theory and the pressure of the Yang-Mills gauge theory at non-zero temperature . Results for these observables are presented in sections 3 and 4 and are discussed along with possible future applications in section 5.

## 2 Jarzynski’s relation

Jarzynski’s equality relates the exponential statistical average of the work done on a system during a non-equilibrium process with the difference in free energy between the initial () and the final () state of the system. In the following, we denote the microscopic degrees of freedom of the system with , so that the partition function of a system with Hamiltonian (which depends also on a certain number of parameters and couplings) can be written as

(2.0) |

We introduce a set of parameters (such as the couplings and/or the temperature itself), which are varied either continuously or discretely during a transformation from an initial value , at which the system is in thermal equilibrium, to a final value . The implementation for a Markov chain in a Monte Carlo simulation is straightforward: the process from the initial to the final state of the system is discretized into steps, each of which is characterized by a certain (with ). After the -th step the set of parameters takes the new value which is used to update the old configuration of the system (obtained using ) with the appropriate algorithm, thus driving it out of equilibrium and never letting it thermalize.

The quantity of interest is the total work done on the system when performing the transformation from to : it can be written as the sum over all the sub-intervals of the difference in the energy, i.e. in the Hamiltonian:

(2.0) |

where and ; the energies are evaluated using the same configuration , before the system is updated with the new parameter . Jarzynski’s relation states the equality of the exponential average of the work over all possible isothermal transformations between the initial and final states and the exponential of the difference in free energy. It can be written as

(2.0) |

and a generalization for non-isothermal transformations has been introduced in ref. [5]. We remark that, in practice in a Monte Carlo simulation, the average is realized over a sufficiently large number of trajectories (which we denote with ) from the initial to the final state.

If we carry out a mapping of this relation from statistical mechanics to lattice gauge theory, we associate with the Euclidean action and transform the quantity appearing in eq. 2 into

(2.0) |

and we can reexpress Jarzynski’s relation using the ratio of the partition functions

(2.0) |

There is an important remark to be made concerning the convergence of this method to the correct result, since a systematic uncertainty appears between the result obtained when performing a trajectory in a certain direction (“direct” transformation) and in the opposite one (“reverse” transformation). This discrepancy depends both on the discretization of the transformation into a finite number of steps and on the finite number of realizations: thus, the combination of and is chosen in order to meet the desired level of uncertainty (which can be set, for example, to be negligible with respect to the statistical error) while also minimizing the computational cost. The determination of the optimal choice of and is a nontrivial problem which depends on the details of the system and/or the simulation; for a thorough discussion on the best practices to deal with this systematic uncertainty we refer to ref. [6].

## 3 Benchmark study I: interface free energy in gauge model

In the first part of this study we applied Jarzynski’s relation to compute the free energy associated with the creation of an interface in the lattice gauge theory in 3 dimensions. The study of interfaces can be very insightful in high-energy physics: in particular they can be related to the world-sheet spanned by flux tubes in confining gauge theories and as such they can be analyzed both from the perspective of an effective theory and through numerical simulations on the lattice.

In the three-dimensional gauge model interfaces can be studied with extraordinary precision; here the degrees of freedom are variables defined on the links between nearest-neighbour sites of a cubic lattice. Remarkably, a confining phase exists for low values of the Wilson parameter . The Wilson action describing the dynamics of this model can be exactly rewritten using the Kramers-Wannier duality as the three-dimensional Ising model on the dual lattice, whose Hamiltonian reads

(3.0) |

where are variables defined on the sites of the lattice, are (anti)ferromagnetic couplings from the site in the direction , and .

An interface can be created by inducing a frustration on the system, i.e. by imposing the condition for the couplings in a chosen direction and for a specific slice only of the lattice, while setting all the remaining ones to . This is equivalent to imposing antiperiodic boundary conditions in one direction: the free energy of the interface (denoted as ) created this way can thus be defined as

(3.0) |

where and are the partition functions of the system with antiperiodic and periodic boundary conditions respectively; the factor accounts for the possibility for the interface to be located anywhere in the direction. An improved definition () which accounts for multiple interfaces has been introduced in [7].

The ratio can be evaluated using Jarzynski’s relation (eq. 2) by identifying the couplings used to create the frustration as the parameters which are varied during the non-equilibrium transformation. More specifically, we vary such couplings using

(3.0) |

which interpolates from to linearly; a similar implementation of Jarzynski’s relation was used in recent works on the 2-dimensional Ising model [8, 5, 9].

Results for the interface free energy obtained in Monte Carlo simulations using the “direct” (switching the couplings from to ) and the “reverse” (from to ) transformations are presented in fig. 1. They clearly converge on the same value of at fixed value of when the number of discretization steps is large enough; moreover, they show an excellent agreement with the results calculated in ref. [7] using thermodynamic integration.

After assessing the reliability of this method, we obtained high precision results at fixed for different values of the interface size using intermediate steps for each out-of-equilibrium transformation and carrying out the average of eq. 2 over different trajectories. These results can be compared with the analytical prediction of the effective string model that describes the transverse fluctuations of the interface at low energy. In particular one can look at the difference between numerical results and the Nambu-Gotō action prediction and examine its dependence on , in order to understand the nature of the terms that do not arise from the Nambu-Gotō low-energy expansion. For the details of this analysis we refer to section 3 of ref. [1].

## 4 Benchmark study II: pressure in gauge theory

In the second part of this study we focus on the determination of the equation of state in the confining phase of the Yang–Mills theory in four spacetime dimensions. We discretize it on a hypercubic lattice of spacing using Wilson’s action; the temperature is defined via , where is the extent of the periodic, compactified Euclidean-time dimension, while we take the lattice sizes in the three other directions to be equal () and sufficiently large to avoid finite-volume effects. Note that in order to control the temperature of the system, we used the relation between and determined in ref. [10], so that we were able to change the temperature simply by varying at fixed . In this work we focus on the computation of the pressure , which in the thermodynamic limit () equals minus the free-energy density ().

One of the most popular techniques use to calculate the pressure non-perturbatively on the lattice is the “integral method”, introduced in ref. [11]. In a nutshell, the pressure as a function of the temperature is determined by integrating over plaquette expectation values computed on lattices of size

(4.0) |

where the lower integration limit corresponds to a temperature low enough at which the pressure is negligible. Moreover, a quartic ultraviolet divergence has been removed by subtracting the value of the plaquette at (denoted as ) computed on a symmetric lattice with size .

Jarzynski’s relation can be naturally extended to the computation of pressure differences by performing non-equilibrium transformations in Monte Carlo simulations on a lattice in which the role of the parameter is taken by the Wilson parameter . The transformation starts at a certain value (which corresponds to a certain temperature ) which is changed linearly after each update of the lattice variables until it reaches the final value (corresponding to the desired temperature ). In this way eq. 2 can be rewritten in order to compute differences in pressure:

(4.0) |

where represents the total variation in the Wilson action (eq. 2) in a transformation from to . Like in eq. 4, the divergence has been removed: the contribution is calculated performing the same transformation on a lattice and then subtracting the corresponding total difference in the action; the exponent is the ratio between the sizes of the two lattices. Results for the confining phase have been calculated using eq. 4 and reported in fig. 2: very good convergence can be observed between “direct” and “reverse” transformations and moreover they show excellent agreement with older results obtained with the integral method from ref. [10]. The transformations were performed independently from one value of to the next and they were discretized using either or intermediate steps; different realizations were carried out in order to compute the exponential average.

The implementation of Jarzynski’s relation has proved to be more efficient than the integral method, since only a fraction of the configurations was needed to obtain results with comparable errors and to keep the systematic uncertainty under control.

## 5 Discussion and future applications

In this work we showed how an extension of Jarzynski’s relation can be used to compute free-energy differences in lattice gauge theories. This novel method successfully reproduced results obtained with other techniques both for the interface free energy in the gauge model and the pressure in the Yang-Mills theory. Convergence to the correct result is achieved by increasing the number of realizations of the non-equilibrium transformation and/or the number of intermediate steps between initial and final state; under such conditions the efficiency of this method proved to be very competitive (and in many cases clearly superior) to that of other algorithms.

Since this novel method is very general and does not require strong assumptions, we envision a number of future applications of this relation with a particular attention to lattice gauge theories with dynamical fermions.
Specifically, we want to emphasize the possibility of applying Jarzynski’s relation to studies involving the Schrödinger functional [12, 13] for the computation of the renormalized coupling : it could be used to compute changes in the effective action induced by changes in the parameters that specify the initial and final states on the boundaries of the lattice in the Euclidean time direction.
Another interesting application would be in the determination of the magnetic susceptibility of QCD in the presence of a strong background magnetic field , which requires the computation of a free energy density difference between different values of (see for example ref. [14]).
We leave these and other potential applications to future studies.

Acknowledgements.

The work of A. T. is partially supported by the Danish National Research Foundation grant DNRF90.

## References

- [1] M. Caselle, G. Costagliola, A. Nada, M. Panero, and A. Toniato, Jarzynski’s theorem for lattice gauge theory, Phys. Rev. D94 (2016), no. 3 034503, [arXiv:1604.0554].
- [2] C. Jarzynski, Nonequilibrium Equality for Free Energy Differences, Phys. Rev. Lett. 78 (1997) 2690–2693, [cond-mat/9610209].
- [3] C. Jarzynski, Equilibrium free-energy differences from nonequilibrium measurements: A master-equation approach, Phys. Rev. E56 (1997) 5018–5035, [cond-mat/9707325].
- [4] J. Liphardt, S. Dumont, S. B. Smith, I. Tinoco, and C. Bustamante, Equilibrium Information from Nonequilibrium Measurements in an Experimental Test of Jarzynski’s Equality, Science 296 (2002), no. 5574 1832–1835.
- [5] C. Chatelain, Temperature-extended Jarzynski relation: Application to the numerical calculation of the surface tension, J. Stat. Mech. 0704 (2007) P04011, [cond-mat/0702044].
- [6] A. Pohorille, C. Jarzynski, and C. Chipot, Good Practices in Free-Energy Calculations, J. Phys. Chem. B114 (2010) 10235–10253.
- [7] M. Caselle, M. Hasenbusch, and M. Panero, The Interface free energy: Comparison of accurate Monte Carlo results for the 3D Ising model with effective interface models, JHEP 0709 (2007) 117, [arXiv:0707.0055].
- [8] C. Chatelain and D. Karevski, Probability distributions of the work in the two-dimensional Ising model, Journal of Statistical Mechanics: Theory and Experiment 2006 (2006), no. 06 P06005, [cond-mat/0602580].
- [9] H. Híjar, J. Quintana-H., and G. Sutmann, Non-equilibrium work theorems for the two-dimensional Ising model, Journal of Statistical Mechanics: Theory and Experiment 2007 (2007), no. 04 P04010.
- [10] M. Caselle, A. Nada, and M. Panero, Hagedorn spectrum and thermodynamics of SU(2) and SU(3) Yang-Mills theories, JHEP 07 (2015) 143, [arXiv:1505.0110].
- [11] J. Engels, J. Fingberg, F. Karsch, D. Miller, and M. Weber, Nonperturbative thermodynamics of SU(N) gauge theories, Phys. Lett. B252 (1990) 625–630.
- [12] K. Symanzik, Schrödinger Representation and Casimir Effect in Renormalizable Quantum Field Theory, Nucl. Phys. B190 (1981) 1.
- [13] M. Lüscher, Schrödinger Representation in Quantum Field Theory, Nucl. Phys. B254 (1985) 52–57.
- [14] C. Bonati, M. D’Elia, M. Mariti, F. Negro, and F. Sanfilippo, Magnetic susceptibility and equation of state of QCD with physical quark masses, Phys. Rev. D89 (2014), no. 5 054506, [arXiv:1310.8656].