Field-theoretic simulation


A field-theoretic simulation is a numerical strategy to calculate structure and physical properties of a many-particle system within the framework of a statistical field theory, like e.g. a polymer field theory. A convenient possibility is to use Monte Carlo algorithms, to sample the full partition function integral expressed in field-theoretic representation. The
procedure is then called the auxiliary field Monte Carlo method. However, it is well known that MC sampling in conjunction with the basic field-theoretic representation of the partition function integral, directly obtained via the Hubbard-Stratonovich transformation, is impracticable, due to the so-called numerical sign problem. The
difficulty is related to the complex and oscillatory nature of the resulting distribution function, which causes a bad statistical convergence of the ensemble averages of the desired structural and thermodynamic quantities. In such cases special analytical and numerical techniques are required to accelerate the statistical convergence of the field-theoretic simulation.

Shifted-contour Monte Carlo technique

Mean field representation

To make the field-theoretic methodology amenable for computation, Baeurle proposed to shift the contour of integration of the partition function integral through the homogeneous mean field solution using Cauchy's integral theorem, which provides its so-called mean-field representation. This strategy was previously successfully employed in field-theoretic electronic structure calculations. Baeurle could demonstrate that this technique provides a significant acceleration of the statistical convergence of the ensemble averages in the MC sampling procedure.

Gaussian equivalent representation

In subsequent works Baeurle et al. applied the concept of tadpole renormalization, which originates from quantum field theory and leads to the Gaussian equivalent representation of the partition function integral, in conjunction with advanced MC techniques in the grand canonical ensemble. They could convincingly demonstrate that this strategy provides an additional boost in the statistical convergence of the desired ensemble averages.

Alternative techniques

Other promising field-theoretic simulation techniques have been developed recently, but they either still lack the proof of correct statistical convergence, like e.g. the Complex Langevin method, and/or still need to prove their effectiveness on systems, where multiple saddle points are important.