Abstract:We present a method for performing path integral molecular dynamics (PIMD) simulations for fermions and address its sign problem. PIMD simulations are widely used for studying many-body quantum systems at thermal equilibrium. However, they assume that the particles are distinguishable and neglect bosonic and fermionic exchange effects. Interacting fermions play a key role in many chemical and physical systems, such as electrons in quantum dots and ultracold trapped atoms. A direct sampling of the fermionic partition function is impossible using PIMD since its integrand is not positive definite. We show that PIMD simulations for fermions are feasible by employing our recently developed method for bosonic PIMD and reweighting the results to obtain fermionic expectation values. The approach is tested against path integral Monte Carlo (PIMC) simulations for up to 7 electrons in a two-dimensional quantum dot for a range of interaction strengths. However, like PIMC, the method suffers from the sign problem at low temperatures. We propose a simple approach for alleviating it by simulating an auxiliary system with a larger average sign and obtaining an upper bound to the energy of the original system using the Bogoliubov inequality. This allows fermions to be studied at temperatures lower than would otherwise have been feasible using PIMD, as demonstrated in the case of a three-electron quantum dot. Our results extend the boundaries of PIMD simulations of fermions and will hopefully stimulate the development of new approaches for tackling the sign problem.

Some bijective correspondences involving domino tableaux

Meron-Cluster Simulation of Quantum Spin Ladders in a Magnetic Field

Meron-Cluster Algorithms for Quantum Link Models

Solving the Fermion Sign Problem in Quantum Monte Carlo Simulations by Majorana Representation

Quantum Monte Carlo for gauge fields and matter without the fermion determinant

A Generic Solution of Fermion Sign Problem

Sign-Problem-Free Fermionic Quantum Monte Carlo: Developments and Applications

Majorana-Time-Reversal Symmetries: A Fundamental Principle for Sign-Problem-Free Quantum Monte Carlo Simulations

Preempting Fermion Sign Problem: Unveiling Quantum Criticality through Nonequilibrium Dynamics

From Spin Ladders to the 2-d O(3) Model at Non-Zero Density

Fermionic Sign Problem Minimization by Constant Path Integral Contour Shifts

Semigroup approach to the sign problem in quantum Monte Carlo simulations

Mitigating the fermion sign problem by automatic differentiation

Robust Fermi liquid instabilities in sign problem-free models

Path Integral Molecular Dynamics for Fermions: Alleviating the Sign Problem with the Bogoliubov Inequality

Fermion Bag Approach to Fermion Sign Problems

Asymptotic sign free in interacting fermion models

Quantum Monte Carlo simulations of highly frustrated magnets in a cluster basis: The two-dimensional Shastry-Sutherland model

Cluster-Algorithm-Amenable Models of Gauge Fields and Matter

Sign-problem-free effective models of triangular lattice quantum antiferromagnets