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

Physical Review C

Circumventing the odd-particle-number sign problem in the shell model Monte Carlo method

Y. Alhassid, P. Fanto, C. Özen

Published Dec 23, 2024• 1 citation
Publisher

Abstract

The shell model Monte Carlo method is a powerful method for calculating exactly (up to statistical errors) thermal observables and statistical properties of atomic nuclei. However, its application has been limited by a sign problem at low temperatures that arises from the projection onto odd number of particles even for good-sign interactions. Here, we develop a technique---the partition function extrapolation method (PFEM)---to extract the ground-state energy of an odd-mass nucleus from the excitation partition function calculated at temperatures at which this sign problem is moderate. We validate the PFEM in heavy even-mass nuclei and systematically calculate ground-state energies for isotopic chains of heavy odd-mass nuclei. The PFEM can be extended to other finite-size quantum many-body systems.

DOI
10.1103/physrevc.110.l061303

Related papers

Circumventing the odd particle-number sign problem in the shell model Monte Carlo

Y. Alhassid, P. Fanto, C. Özen

The shell model Monte Carlo (SMMC) method is a powerful method for calculating exactly (up to statistical errors) thermal observables and statistical properties of atomic nuclei. However, its application has been limited by a sign problem at low temperatures that arises from the projection onto odd particle number even for good-sign interactions. Here, we develop a technique - the partition function extrapolation method (PFEM) - to extract the ground-state energy of an odd-mass nucleus from the excitation partition function calculated at temperatures at which this sign problem is moderate. We validate the PFEM in heavy even-mass nuclei and systematically calculate ground-state energies for isotopic chains of heavy odd-mass nuclei. The PFEM can be extended to other finite-size quantum many-body systems.

arXiv:2308.14698•Aug 28, 2023•0 citations

Odd-Particle Systems in the Shell Model Monte Carlo Method: Circumventing a Sign Problem

Abhishek Mukherjee, Y. Alhassid

The shell model Monte Carlo method is a powerful technique to calculate thermal and ground-state properties of strongly correlated finite-size systems. However, its application to odd-particle-number systems has been hampered by the sign problem that originates from the projection on an odd number of particles. We circumvent this sign problem for the ground-state energy by extracting the ground-state energy of the odd-particle-number system from the asymptotic behavior of the imaginary-time single-particle Green's function of the even-particle-number system. We apply this method to calculate pairing gaps of nuclei in the iron region. Our results are in good agreement with experimental pairing gaps.

arXiv:1201.3341•Physical Review Letters•Jul 17, 2012•14 citations

Recent Advances in the Microscopic Calculations of Level Den- sities by the Shell Model Monte Carlo Method

Y. Alhassid, M. Bonett-Matiz, S. Liu +2 more

The shell model Monte Carlo (SMMC) method enables calculations in model spaces that are many orders of magnitude larger than those that can be treated by conventional methods, and is particularly suitable for the calculation of level densities in the presence of correlations. We review recent advances and applications of SMMC for the microscopic calculation of level densities. Recent developments include (i) a method to calculate accurately the ground-state energy of an odd-mass nucleus, circumventing a sign problem that originates in the projection on an odd number of particles, and (ii) a method to calculate directly level densities, which, unlike state densities, do not include the spin degeneracy of the levels. We calculated the level densities of a family of nickel isotopes 59−64Ni and of a heavy deformed rare-earth nucleus 162Dy and found them to be in close agreement with various experimental data sets.

•Apr 1, 2014•3 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

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

J. Bonnard, Olivier Juillet

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.

•Jan 1, 2012•0 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.

arXiv:nucl-th/0009005•International Journal of Modern Physics B•May 10, 2001•29 citations
The shell model Monte Carlo method is a powerful method for calculating exactly (up to statistical errors) thermal observables and statistical properties of atomic nuclei. However, its application has been limited by a sign problem at low temperatures that arises from the projection onto odd number of particles even for good-sign interactions. Here, we develop a technique---the partition function extrapolation method (PFEM)---to extract the ground-state energy of an odd-mass nucleus from the excitation partition function calculated at temperatures at which this sign problem is moderate. We validate the PFEM in heavy even-mass nuclei and systematically calculate ground-state energies for isotopic chains of heavy odd-mass nuclei. The PFEM can be extended to other finite-size quantum many-body systems.
The shell model Monte Carlo (SMMC) method is a powerful method for calculating exactly (up to statistical errors) thermal observables and statistical properties of atomic nuclei. However, its application has been limited by a sign problem at low temperatures that arises from the projection onto odd particle number even for good-sign interactions. Here, we develop a technique - the partition function extrapolation method (PFEM) - to extract the ground-state energy of an odd-mass nucleus from the excitation partition function calculated at temperatures at which this sign problem is moderate. We validate the PFEM in heavy even-mass nuclei and systematically calculate ground-state energies for isotopic chains of heavy odd-mass nuclei. The PFEM can be extended to other finite-size quantum many-body systems.
The shell model Monte Carlo method is a powerful technique to calculate thermal and ground-state properties of strongly correlated finite-size systems. However, its application to odd-particle-number systems has been hampered by the sign problem that originates from the projection on an odd number of particles. We circumvent this sign problem for the ground-state energy by extracting the ground-state energy of the odd-particle-number system from the asymptotic behavior of the imaginary-time single-particle Green's function of the even-particle-number system. We apply this method to calculate pairing gaps of nuclei in the iron region. Our results are in good agreement with experimental pairing gaps.
The shell model Monte Carlo (SMMC) method enables calculations in model spaces that are many orders of magnitude larger than those that can be treated by conventional methods, and is particularly suitable for the calculation of level densities in the presence of correlations. We review recent advances and applications of SMMC for the microscopic calculation of level densities. Recent developments include (i) a method to calculate accurately the ground-state energy of an odd-mass nucleus, circumventing a sign problem that originates in the projection on an odd number of particles, and (ii) a method to calculate directly level densities, which, unlike state densities, do not include the spin degeneracy of the levels. We calculated the level densities of a family of nickel isotopes 59−64Ni and of a heavy deformed rare-earth nucleus 162Dy and found them to be in close agreement with various experimental data sets.
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.
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.