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

arXiv:nucl-th/9512027 · Physical Review Letters

Monte Carlo Computation of Pair Correlations in Excited Nuclei

Nicolas J. Cerf

Published Apr 1, 1996• 7 citations
arXivHTMLPDF

Abstract

We present a novel quantum Monte Carlo method based on a path integral in Fock space, which allows the computation of finite-temperature properties of a many-body nuclear system with a monopole pairing interaction in the canonical ensemble. It enables an exact calculation of the thermodynamic variables such as the internal energy, the entropy, or the specific heat, from the measured moments of the number of hops in a path of nuclear configurations. Monte Carlo calculations for a single-shell $({h}_{11/2}{)}^{6}$ model are consistent with an exact calculation from the many-body spectrum in the seniority model.

DOI
10.1103/physrevlett.76.2420

Related papers

Monte Carlo evaluation of path integrals for the nuclear shell model

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

We present in detail a formulation of the shell model as a path integral and Monte Carlo techniques for its evaluation. The formulation, which linearizes the two-body interaction by an auxiliary field, is quite general, both in the form of the effective ``one-body'' Hamiltonian and in the choice of ensemble. In particular, we derive formulas for the use of general (beyond monopole) pairing operators, as well as a novel extraction of the canonical (fixed-particle-number) ensemble via an activity expansion. We discuss the advantages and disadvantages of the various formulations and ensembles and give several illustrative examples. We also discuss and illustrate calculation of the imaginary-time response function and the extraction, by maximum entropy methods, of the corresponding strength function. Finally, we discuss the ``sign problem'' generic to fermion Monte Carlo calculations, and prove that a wide class of interactions are free of this limitation.

arXiv:nucl-th/9305009•Physical Review C•Oct 1, 1993•239 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

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

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

Finite-size effects in canonical and grand-canonical quantum Monte Carlo simulations for fermions

Zhenjiu Wang, Fakher F. Assaad, Francesco Parisen Toldin

We introduce a quantum Monte Carlo method at finite temperature for interacting fermionic models in the canonical ensemble, where the conservation of the particle number is enforced. Although general thermodynamic arguments ensure the equivalence of the canonical and the grand-canonical ensembles in the thermodynamic limit, their approach to the infinite-volume limit is distinctively different. Observables computed in the canonical ensemble generically display a finite-size correction proportional to the inverse volume, whereas in the grand-canonical ensemble the approach is exponential in the ratio of the linear size over the correlation length. We verify these predictions by quantum Monte Carlo simulations of the Hubbard model in one and two dimensions in the grand-canonical and the canonical ensemble. We prove an exact formula for the finite-size part of the free energy density, energy density and other observables in the canonical ensemble and relate this correction to a susceptibility computed in the corresponding grand-canonical ensemble. This result is confirmed by an exact computation of the one-dimensional classical Ising model in the canonical ensemble, which for classical models corresponds to the so-called fixed-magnetization ensemble. Our method is useful for simulating finite systems which are not coupled to a particle bath, such as in nuclear or cold atom physics.

arXiv:1706.01874•Physical review. E•Oct 16, 2017•24 citations

Quantum Monte Carlo simulation in the canonical ensemble at finite temperature

Kris Van Houcke, S. Rombouts, Lode Pollet

A quantum Monte Carlo method with a nonlocal update scheme is presented. The method is based on a path-integral decomposition and a worm operator which is local in imaginary time. It generates states with a fixed number of particles and respects other exact symmetries. Observables like the equal-time Green's function can be evaluated in an efficient way. To demonstrate the versatility of the method, results for the one-dimensional Bose-Hubbard model and a nuclear pairing model are presented. Within the context of the Bose-Hubbard model the efficiency of the algorithm is discussed.

arXiv:cond-mat/0508431•Physical Review E•May 16, 2006•39 citations
We present a novel quantum Monte Carlo method based on a path integral in Fock space, which allows the computation of finite-temperature properties of a many-body nuclear system with a monopole pairing interaction in the canonical ensemble. It enables an exact calculation of the thermodynamic variables such as the internal energy, the entropy, or the specific heat, from the measured moments of the number of hops in a path of nuclear configurations. Monte Carlo calculations for a single-shell (h11/2)6({h}_{11/2}{)}^{6}(h11/2​)6 model are consistent with an exact calculation from the many-body spectrum in the seniority model.
We present in detail a formulation of the shell model as a path integral and Monte Carlo techniques for its evaluation. The formulation, which linearizes the two-body interaction by an auxiliary field, is quite general, both in the form of the effective ``one-body'' Hamiltonian and in the choice of ensemble. In particular, we derive formulas for the use of general (beyond monopole) pairing operators, as well as a novel extraction of the canonical (fixed-particle-number) ensemble via an activity expansion. We discuss the advantages and disadvantages of the various formulations and ensembles and give several illustrative examples. We also discuss and illustrate calculation of the imaginary-time response function and the extraction, by maximum entropy methods, of the corresponding strength function. Finally, we discuss the ``sign problem'' generic to fermion Monte Carlo calculations, and prove that a wide class of interactions are free of this limitation.
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.
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.
We introduce a quantum Monte Carlo method at finite temperature for interacting fermionic models in the canonical ensemble, where the conservation of the particle number is enforced. Although general thermodynamic arguments ensure the equivalence of the canonical and the grand-canonical ensembles in the thermodynamic limit, their approach to the infinite-volume limit is distinctively different. Observables computed in the canonical ensemble generically display a finite-size correction proportional to the inverse volume, whereas in the grand-canonical ensemble the approach is exponential in the ratio of the linear size over the correlation length. We verify these predictions by quantum Monte Carlo simulations of the Hubbard model in one and two dimensions in the grand-canonical and the canonical ensemble. We prove an exact formula for the finite-size part of the free energy density, energy density and other observables in the canonical ensemble and relate this correction to a susceptibility computed in the corresponding grand-canonical ensemble. This result is confirmed by an exact computation of the one-dimensional classical Ising model in the canonical ensemble, which for classical models corresponds to the so-called fixed-magnetization ensemble. Our method is useful for simulating finite systems which are not coupled to a particle bath, such as in nuclear or cold atom physics.
A quantum Monte Carlo method with a nonlocal update scheme is presented. The method is based on a path-integral decomposition and a worm operator which is local in imaginary time. It generates states with a fixed number of particles and respects other exact symmetries. Observables like the equal-time Green's function can be evaluated in an efficient way. To demonstrate the versatility of the method, results for the one-dimensional Bose-Hubbard model and a nuclear pairing model are presented. Within the context of the Bose-Hubbard model the efficiency of the algorithm is discussed.