Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

International Journal of Modern Physics B

A CONSTRAINED PATH MONTE CARLO METHOD FOR NUCLEON SYSTEMS

Kevin E. Schmidt, A. Sarsa, Stefano Fantoni

Published May 10, 2001• 1 citation
Publisher

Abstract

By combining diffusion Monte Carlo for the spatial degrees of freedom and auxiliary field Monte Carlo to separate the spin-isospin operators, we can solve for the ground state of many-nucleon systems. We use a path constraint to control the fermion sign problem and apply the method to neutron systems interacting with the Argonne v′8 two nucleon potential and the Urbana IX three-nucleon potential. We compare our results with fermion hypernetted chain calculations.

DOI
10.1142/s0217979201005994

Related papers

Fermion ground state methods

James E. Gubernatis, Naoki Kawashima, Philipp Werner

In this chapter, we discuss the fixed-node and constrained-path Monte Carlo methods for computing the ground state properties of systems of interacting electrons. These methods are arguably the two most powerful ones presently available for doing such calculations, but they are approximate. By sacrificing exactness, they avoid the exponential scaling of the Monte Carlo errors with system size that typically accompanies the simulation of systems of interacting electrons. This exponential scaling is called the Fermion sign problem . After a few general comments about the sign problem, we outline both methods, noting points of similarity and difference, plus points of strength and weakness. We also discuss the constrained-phase method, an extension of the constrained-path method, which controls the phase problem that develops when the ground state wave function cannot be real. Sign problem The “sign problem” refers to the exponential increase of the Monte Carlo errors with increasing system size or decreasing temperature (e.g., Loh et al., 1990, 2005; Gubernatis and Zhang, 1994) that often accompanies a Markov chain simulation whose limiting distribution is not everywhere positive. Such a case generally arises in simulations of Fermion and frustrated quantum-spin systems. It seems so inherent to Monte Carlo simulations of Fermion systems that the phrase “the sign problem” to many seems almost synonymous with the phrase “the Fermion sign problem.” Explanations for the cause of the sign problem vary and are still debated. In this chapter, we choose to summarize two explanations that seem to connect the causes in ground state Fermion simulations in the continuum and on the lattice. The sign problem, of course, is not limited to ground state calculations or even to Fermion simulations. While the cause we discuss in Sections 11.2 and 11.3 focuses on the low-lying states of diffusion-like operators, several topological pictures have been proposed (Muramatsu et al., 1992; Samson, 1993; Gubernatis and Zhang, 1994; Samson, 1995). Some of these discussions are done in the context of the zero- and finite-temperature determinant methods (Muramatsu et al., 1992; Gubernatis and Zhang, 1994). Others are done more analytically from a Feynman path-integral point of view (Samson, 1993, 1995). Some are for particles with statistics other than Fermions. The presentation of the sign problem in this chapter is appropriate for the Monte Carlo methods discussed in this chapter.

Cambridge University Press eBooks•Jan 31, 2016•0 citations

MONTE CARLO DIAGONALIZATION OF VERY LARGE MATRICES: APPLICATION TO FERMION SYSTEMS

Hans A. De Raedt, Wolfgang von der Linden

All known Quantum-Monte-Carlo algorithms for fermions suffer from the so-called “minus-sign-problem” which is detrimental to the application of these simulation methods to fermion systems at very low temperatures and/or of very many lattice sites. We identify the origin of this fundamental problem, demonstrate that it is a very general feature, not necessarily related to the presence of fermionic degrees of freedom. We describe an novel algorithm which does not suffer from the minus-sign problem. Illustrative results for the two-dimensional Hubbard model are presented

International Journal of Modern Physics C•Feb 1, 1992•6 citations

Green's function Monte Carlo calculations of few nucleon systems

Brian Pudliner

Recently, a realistic nucleon-nucleon potential has been developed with the goal of accurately reproducing isospin symmetry breaking effects in the nuclear interaction. We briefly discuss the motivation for this potential, called Argonne $v\sb{18},$ as well as the Urbana models of three-nucleon interactions. Our approach for modeling the ground states of nuclear systems begins with the development of a variational wave function. We discuss the basis for our variational wave functions and the Monte Carlo techniques used to evaluate expectation values with these wave functions. To determine the exact ground state, Green's function Monte Carlo (GFMC) techniques are used. We review the method of GFMC in general and then discuss its application to the problem of nuclear ground states in some detail.

•Jan 1, 1996•1 citation

Quantum Monte Carlo: What can random numbers tell us about correlated quantum systems?

S. Rombouts, K. Heyde

Abstract Over the past decades, quantum Monte Carlo methods have proven to be very valuable for the study of quantum‐mechanical many‐body systems. The basic ingredient of theses methods is a decomposition of the imaginary‐time evolution operator into path integrals. These paths can then be sampled using Monte Carlo techniques. Though very successful for bosons, these methods often have troubles in dealing with fermions: the fermionic antisymmetry leads to a ‘sign problem’. One approach, the auxiliary‐field method, has yielded slightly better results in dealing with the sign problem. Some recent developments on this method are highlighted. As a case study, it is applied to a pairing model for Fe nuclei, where indications are found for a phase transition related to nucleon–nucleon pairing.

physica status solidi (b)•Apr 24, 2003•5 citations

Quantum Monte Carlo calculations of neutron matter

Joseph Carlson, J. Morales, Vijay R. Pandharipande +1 more

Uniform neutron matter is approximated by a cubical box containing a finite number of neutrons, with periodic boundary conditions. We report variational and Green's function Monte Carlo calculations of the ground state of fourteen neutrons in a periodic box using the Argonne ${v8}^{\ensuremath{'}}$ two-nucleon interaction at densities up to one and half times the nuclear matter density. The effects of the finite box size are estimated using variational wave functions together with cluster expansion and chain summation techniques. They are small at subnuclear densities. We discuss the expansion of the energy of low-density neutron gas in powers of its Fermi momentum. This expansion is strongly modified by the large $\mathrm{nn}$ scattering length, and does not begin with the Fermi-gas kinetic energy, as assumed in both Skyrme and relativistic mean field theories. The leading term of neutron gas energy is approximately half the Fermi-gas kinetic energy. The quantum Monte Carlo results are also used to calibrate the accuracy of variational calculations employing Fermi hypernetted and single operator chain summation methods to study nucleon matter over a larger density range, with more realistic Hamiltonians including three-nucleon interactions.

arXiv:nucl-th/0302041•Physical Review C•Aug 15, 2003•176 citations

Monte Carlo methods for the nuclear shell model

Calvin W. Johnson, S. E. Koonin, G. H. Lang +1 more

We present novel Monte Carlo methods for treating the interacting shell model that allow exact calculations much larger than those heretofore possible. The two-body interaction is linearized by an auxiliary field; Monte Carlo evaluation of the resulting functional integral gives ground-state or thermal expectation values of few-body operators. The ``sign problem'' generic to quantum Monte Carlo calculations is absent in a number of cases. We discuss the favorable scaling of these methods with nucleon number and basis size and their suitability to parallel computation.

arXiv:nucl-th/9210014•Physical Review Letters•Nov 30, 1992•133 citations
By combining diffusion Monte Carlo for the spatial degrees of freedom and auxiliary field Monte Carlo to separate the spin-isospin operators, we can solve for the ground state of many-nucleon systems. We use a path constraint to control the fermion sign problem and apply the method to neutron systems interacting with the Argonne v′8 two nucleon potential and the Urbana IX three-nucleon potential. We compare our results with fermion hypernetted chain calculations.
In this chapter, we discuss the fixed-node and constrained-path Monte Carlo methods for computing the ground state properties of systems of interacting electrons. These methods are arguably the two most powerful ones presently available for doing such calculations, but they are approximate. By sacrificing exactness, they avoid the exponential scaling of the Monte Carlo errors with system size that typically accompanies the simulation of systems of interacting electrons. This exponential scaling is called the Fermion sign problem . After a few general comments about the sign problem, we outline both methods, noting points of similarity and difference, plus points of strength and weakness. We also discuss the constrained-phase method, an extension of the constrained-path method, which controls the phase problem that develops when the ground state wave function cannot be real. Sign problem The “sign problem” refers to the exponential increase of the Monte Carlo errors with increasing system size or decreasing temperature (e.g., Loh et al., 1990, 2005; Gubernatis and Zhang, 1994) that often accompanies a Markov chain simulation whose limiting distribution is not everywhere positive. Such a case generally arises in simulations of Fermion and frustrated quantum-spin systems. It seems so inherent to Monte Carlo simulations of Fermion systems that the phrase “the sign problem” to many seems almost synonymous with the phrase “the Fermion sign problem.” Explanations for the cause of the sign problem vary and are still debated. In this chapter, we choose to summarize two explanations that seem to connect the causes in ground state Fermion simulations in the continuum and on the lattice. The sign problem, of course, is not limited to ground state calculations or even to Fermion simulations. While the cause we discuss in Sections 11.2 and 11.3 focuses on the low-lying states of diffusion-like operators, several topological pictures have been proposed (Muramatsu et al., 1992; Samson, 1993; Gubernatis and Zhang, 1994; Samson, 1995). Some of these discussions are done in the context of the zero- and finite-temperature determinant methods (Muramatsu et al., 1992; Gubernatis and Zhang, 1994). Others are done more analytically from a Feynman path-integral point of view (Samson, 1993, 1995). Some are for particles with statistics other than Fermions. The presentation of the sign problem in this chapter is appropriate for the Monte Carlo methods discussed in this chapter.
All known Quantum-Monte-Carlo algorithms for fermions suffer from the so-called “minus-sign-problem” which is detrimental to the application of these simulation methods to fermion systems at very low temperatures and/or of very many lattice sites. We identify the origin of this fundamental problem, demonstrate that it is a very general feature, not necessarily related to the presence of fermionic degrees of freedom. We describe an novel algorithm which does not suffer from the minus-sign problem. Illustrative results for the two-dimensional Hubbard model are presented
Recently, a realistic nucleon-nucleon potential has been developed with the goal of accurately reproducing isospin symmetry breaking effects in the nuclear interaction. We briefly discuss the motivation for this potential, called Argonne $v\sb{18},$ as well as the Urbana models of three-nucleon interactions. Our approach for modeling the ground states of nuclear systems begins with the development of a variational wave function. We discuss the basis for our variational wave functions and the Monte Carlo techniques used to evaluate expectation values with these wave functions. To determine the exact ground state, Green's function Monte Carlo (GFMC) techniques are used. We review the method of GFMC in general and then discuss its application to the problem of nuclear ground states in some detail.
Abstract Over the past decades, quantum Monte Carlo methods have proven to be very valuable for the study of quantum‐mechanical many‐body systems. The basic ingredient of theses methods is a decomposition of the imaginary‐time evolution operator into path integrals. These paths can then be sampled using Monte Carlo techniques. Though very successful for bosons, these methods often have troubles in dealing with fermions: the fermionic antisymmetry leads to a ‘sign problem’. One approach, the auxiliary‐field method, has yielded slightly better results in dealing with the sign problem. Some recent developments on this method are highlighted. As a case study, it is applied to a pairing model for Fe nuclei, where indications are found for a phase transition related to nucleon–nucleon pairing.
Uniform neutron matter is approximated by a cubical box containing a finite number of neutrons, with periodic boundary conditions. We report variational and Green's function Monte Carlo calculations of the ground state of fourteen neutrons in a periodic box using the Argonne ${v8}^{\ensuremath{'}}$ two-nucleon interaction at densities up to one and half times the nuclear matter density. The effects of the finite box size are estimated using variational wave functions together with cluster expansion and chain summation techniques. They are small at subnuclear densities. We discuss the expansion of the energy of low-density neutron gas in powers of its Fermi momentum. This expansion is strongly modified by the large nn\mathrm{nn}nn scattering length, and does not begin with the Fermi-gas kinetic energy, as assumed in both Skyrme and relativistic mean field theories. The leading term of neutron gas energy is approximately half the Fermi-gas kinetic energy. The quantum Monte Carlo results are also used to calibrate the accuracy of variational calculations employing Fermi hypernetted and single operator chain summation methods to study nucleon matter over a larger density range, with more realistic Hamiltonians including three-nucleon interactions.
We present novel Monte Carlo methods for treating the interacting shell model that allow exact calculations much larger than those heretofore possible. The two-body interaction is linearized by an auxiliary field; Monte Carlo evaluation of the resulting functional integral gives ground-state or thermal expectation values of few-body operators. The ``sign problem'' generic to quantum Monte Carlo calculations is absent in a number of cases. We discuss the favorable scaling of these methods with nucleon number and basis size and their suitability to parallel computation.