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

A phase-free quantum Monte Carlo method for the nuclear shell model

J. Bonnard, Olivier Juillet

Published Jan 1, 2012• 0 citations
OpenAlex

Abstract

The shell model provides a powerful framework for nuclear structure calculations. The nucleons beyond an inert magic core are confined in a valence shell and interact through an effective two-body potential generally determined from the G-matrix method. However, the applicability of the shell model remains limited by the exponential growth of the many-body space with the number of valence nucleons and with the size of the one-body space. The quantum Monte Carlo (QMC) methods may thus be considered as a potentially attractive alternative to the direct diagonalization of the Hamiltonian matrix. The many-body problem is indeed reduced to a set of stochastic one-body problems describing independent particles in a fluctuating external field. Up to date, the QMC approaches of the shell model lead to the ground-state and the thermodynamic properties, but a detailed spectroscopy of nuclei is not reachable. Furthermore, they suffer from a common pathology of fermionic QMC methods, namely the sign/phase problem, which strongly contaminates the convergence of the calculation. In this context, we will present a new QMC method, providing the "yrast spectroscopy" of nuclei, and in which the sign/phase problem is managed through the constrained path approximation that is widely used in the ab initio nuclear calculations and in condensed matter physics.

Related papers

QUANTUM MONTE CARLO METHODS FOR NUCLEI AT FINITE TEMPERATURE

Y. Alhassid

We discuss finite temperature quantum Monte Carlo methods in the framework of the interacting nuclear shell model. The methods are based on a representation of the imaginary-time many-body propagator as a superposition of one-body propagators describing non-interacting fermions moving in fluctuating auxiliary fields. Fermionic Monte Carlo calculations have been limited by a "sign" problem. A practical solution in the nuclear case enables realistic calculations in much larger configuration spaces than can be solved by conventional methods. Good-sign interactions can be constructed for realistic estimates of certain nuclear properties. We present various applications of the methods for calculating collective properties and level densities.

arXiv:nucl-th/0009005•International Journal of Modern Physics B•May 10, 2001•29 citations

QUANTUM MONTE CARLO METHODS FOR NUCLEI AT FINITE TEMPERATURE

Y. Alhassid

We discuss finite temperature quantum Monte Carlo methods in the framework of the interacting nuclear shell model. The methods are based on a representation of the imaginary-time many-body propagator as a superposition of one-body propagators describing non-interacting fermions moving in fluctuating auxiliary fields. Fermionic Monte Carlo calculations have been limited by a ``sign'' problem. A practical solution in the nuclear case enables realistic calculations in much larger configuration spaces than can be solved by conventional methods. Good-sign interactions can be constructed for realistic estimates of certain nuclear properties. We present various applications of the methods for calculating collective properties and level densities.

•Jan 1, 2000•16 citations

MONTE CARLO SHELL MODEL CALCULATIONS FOR ATOMIC NUCLEI

Takaharu Otsuka, T. Mizusaki, Y. Utsuno +2 more

The formulation and recent applications of the Quantum Monte Carlo diagonalization (QMCD) method are reported. The QMCD has been proposed for solving the quantum many-body inteacting systems, providing us with energy eigenvalues, transition matrix elements and wave functions. Its application to the nuclear shell model is referred to as the Monte Carlo Shell Model. By the Monte Carlo Shell Model calculations, the level structure of low-lying states can be studied with realistic interactions, providing a useful tool for nuclear spectroscopy. The Monte Carlo Shell Model has been applied to the study of a variety of nuclei, and can be characterized as the importance truncation scheme to the full diagonalization which is infeasible in many cases due to extremely large dimensions.

International Journal of Modern Physics B•May 10, 2001•0 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

Electronic structure quantum Monte Carlo

Michal Bajdich, Luboš Mitáš

Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.

arXiv:1008.2369•Acta Physica Slovaca Reviews and Tutorials•Apr 1, 2009•33 citations

Auxiliary-field Monte Carlo methods for interacting fermions : application to the nuclear shell model

G. H. Lang

This thesis presents the path-integral formulation of the nuclear shell model using the Hubbard-Stratonovich transformation, which linearizes the two-body interaction by auxiliary fields. The path-integral was evaluated via Monte Carlo. The method scales favorably with valence-nucleon number and shell-model basis: full-basis calculations can be done up to the rare-earth region, which cannot be treated by other methods. Observables are calculated for the ground state and in a thermal ensemble. Dynamical correlations are obtained, from which strength functions are extracted through the Maximum Entropy method. Examples in the s-d shell, where exact diagonalization can be carried out, compare well with exact results. The "sign problem", which is generic to fermion Monte Carlo calculations, is proved to be absent in a wide class of interactions including the attractive pairing-plus-multipole interactions. The formulation is general for interacting fermion systems and is well suited for parallel computation. The method has been implemented on the Intel Touchstone Delta System, achieving better than 99% parallelization.

•Nov 18, 2009•0 citations
The shell model provides a powerful framework for nuclear structure calculations. The nucleons beyond an inert magic core are confined in a valence shell and interact through an effective two-body potential generally determined from the G-matrix method. However, the applicability of the shell model remains limited by the exponential growth of the many-body space with the number of valence nucleons and with the size of the one-body space. The quantum Monte Carlo (QMC) methods may thus be considered as a potentially attractive alternative to the direct diagonalization of the Hamiltonian matrix. The many-body problem is indeed reduced to a set of stochastic one-body problems describing independent particles in a fluctuating external field. Up to date, the QMC approaches of the shell model lead to the ground-state and the thermodynamic properties, but a detailed spectroscopy of nuclei is not reachable. Furthermore, they suffer from a common pathology of fermionic QMC methods, namely the sign/phase problem, which strongly contaminates the convergence of the calculation. In this context, we will present a new QMC method, providing the "yrast spectroscopy" of nuclei, and in which the sign/phase problem is managed through the constrained path approximation that is widely used in the ab initio nuclear calculations and in condensed matter physics.
We discuss finite temperature quantum Monte Carlo methods in the framework of the interacting nuclear shell model. The methods are based on a representation of the imaginary-time many-body propagator as a superposition of one-body propagators describing non-interacting fermions moving in fluctuating auxiliary fields. Fermionic Monte Carlo calculations have been limited by a "sign" problem. A practical solution in the nuclear case enables realistic calculations in much larger configuration spaces than can be solved by conventional methods. Good-sign interactions can be constructed for realistic estimates of certain nuclear properties. We present various applications of the methods for calculating collective properties and level densities.
We discuss finite temperature quantum Monte Carlo methods in the framework of the interacting nuclear shell model. The methods are based on a representation of the imaginary-time many-body propagator as a superposition of one-body propagators describing non-interacting fermions moving in fluctuating auxiliary fields. Fermionic Monte Carlo calculations have been limited by a ``sign'' problem. A practical solution in the nuclear case enables realistic calculations in much larger configuration spaces than can be solved by conventional methods. Good-sign interactions can be constructed for realistic estimates of certain nuclear properties. We present various applications of the methods for calculating collective properties and level densities.
The formulation and recent applications of the Quantum Monte Carlo diagonalization (QMCD) method are reported. The QMCD has been proposed for solving the quantum many-body inteacting systems, providing us with energy eigenvalues, transition matrix elements and wave functions. Its application to the nuclear shell model is referred to as the Monte Carlo Shell Model. By the Monte Carlo Shell Model calculations, the level structure of low-lying states can be studied with realistic interactions, providing a useful tool for nuclear spectroscopy. The Monte Carlo Shell Model has been applied to the study of a variety of nuclei, and can be characterized as the importance truncation scheme to the full diagonalization which is infeasible in many cases due to extremely large dimensions.
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.
Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.
This thesis presents the path-integral formulation of the nuclear shell model using the Hubbard-Stratonovich transformation, which linearizes the two-body interaction by auxiliary fields. The path-integral was evaluated via Monte Carlo. The method scales favorably with valence-nucleon number and shell-model basis: full-basis calculations can be done up to the rare-earth region, which cannot be treated by other methods. Observables are calculated for the ground state and in a thermal ensemble. Dynamical correlations are obtained, from which strength functions are extracted through the Maximum Entropy method. Examples in the s-d shell, where exact diagonalization can be carried out, compare well with exact results. The "sign problem", which is generic to fermion Monte Carlo calculations, is proved to be absent in a wide class of interactions including the attractive pairing-plus-multipole interactions. The formulation is general for interacting fermion systems and is well suited for parallel computation. The method has been implemented on the Intel Touchstone Delta System, achieving better than 99% parallelization.