Monte Carlo evaluation of path integrals for the nuclear shell model

Stable recursive auxiliary field quantum Monte Carlo algorithm in the canonical ensemble: Applications to thermometry and the Hubbard model

Many experimentally accessible, finite-sized interacting quantum systems are most appropriately described by the canonical ensemble of statistical mechanics. Conventional numerical simulation methods either approximate them as being coupled to a particle bath or use projective algorithms which may suffer from nonoptimal scaling with system size or large algorithmic prefactors. In this paper, we introduce a highly stable, recursive auxiliary field quantum Monte Carlo approach that can directly simulate systems in the canonical ensemble. We apply the method to the fermion Hubbard model in one and two spatial dimensions in a regime known to exhibit a significant "sign" problem and find improved performance over existing approaches including rapid convergence to ground-state expectation values. The effects of excitations above the ground state are quantified using an estimator-agnostic approach including studying the temperature dependence of the purity and overlap fidelity of the canonical and grand canonical density matrices. As an important application, we show that thermometry approaches often exploited in ultracold atoms that employ an analysis of the velocity distribution in the grand canonical ensemble may be subject to errors leading to an underestimation of extracted temperatures with respect to the Fermi temperature.

Pairing in excited nuclei: a review

The present review summarizes the recent studies on the thermodynamic properties of pairing in many-body systems including superconductors, metallic nanosized clusters and/or grains, solid-state materials, focusing on the excited nuclei, that is nuclei at finite temperature and/or angular momentum formed via heavy-ion fusion, [Formula: see text]-induced fusion reactions, or inelastic scattering of light particles on heavy targets. Because of the finiteness of the systems, several interesting effects of pairing such as nonvanishing pairing gap, smoothing of superfluid-normal phase transition, first and second order phase transitions, pairing reentrance, etc, will be discussed in detail. Influences of exact and approximate thermal pairing on some nuclear properties such as temperature-dependent width of the giant dipole resonance, total level density, and radiative strength function of the [Formula: see text]-rays emission will be also analyzed. Finally, the first experimental evidence of the pairing reentrance phenomenon in a ¹⁰⁴Pd nucleus as well as its solid-state counterpart of ferromagnets under strong magnetic field will be presented.

Nuclear level density and related physics

The knowledge of the level density as a function of excitation energy and nuclear spin is necessary for the description of nuclear reactions and in many applied areas. We discuss the level density problem as a part of the general understanding of mesoscopic systems with strong interactions. The underlying physics is that of quantum chaos and thermalization which allows one to use statistical methods avoiding full diagonalization. The resulting level density is well described by the constant temperature model in agreement with experimental data. We discuss the effective temperature parameter and show that it is not related to the pairing phase transition being analogous to the limiting temperature in particle physics. Other aspects of underlying physics include the collective enhancement of the level density, random coupling of individual spins and the role of incoherent collision-like interactions.

One-Shot Decoupling and Page Curves from a Dynamical Model for Black Hole Evaporation

One-shot decoupling is a powerful primitive in quantum information theory and was hypothesized to play a role in the black hole information paradox. We study black hole dynamics modeled by a trilinear Hamiltonian whose semiclassical limit gives rise to Hawking radiation. An explicit numerical calculation of the discretized path integral of the S matrix shows that decoupling is exact in the continuous limit, implying that quantum information is perfectly transferred from the black hole to radiation. A striking consequence of decoupling is the emergence of an output radiation entropy profile that follows Page's prediction. We argue that information transfer and the emergence of Page curves is a robust feature of any multilinear interaction Hamiltonian with a bounded spectrum.

Split Orthogonal Group: A Guiding Principle for Sign-Problem-Free Fermionic Simulations

We present a guiding principle for designing fermionic Hamiltonians and quantum Monte Carlo (QMC) methods that are free from the infamous sign problem by exploiting the Lie groups and Lie algebras that appear naturally in the Monte Carlo weight of fermionic QMC simulations. Specifically, rigorous mathematical constraints on the determinants involving matrices that lie in the split orthogonal group provide a guideline for sign-free simulations of fermionic models on bipartite lattices. This guiding principle not only unifies the recent solutions of the sign problem based on the continuous-time quantum Monte Carlo methods and the Majorana representation, but also suggests new efficient algorithms to simulate physical systems that were previously prohibitive because of the sign problem.

Nuclear deformation at finite temperature

Deformation, a key concept in our understanding of heavy nuclei, is based on a mean-field description that breaks the rotational invariance of the nuclear many-body Hamiltonian. We present a method to analyze nuclear deformations at finite temperature in a framework that preserves rotational invariance. The auxiliary-field Monte Carlo method is used to generate a statistical ensemble and calculate the probability distribution associated with the quadrupole operator. Applying the technique to nuclei in the rare-earth region, we identify model-independent signatures of deformation and find that deformation effects persist to temperatures higher than the spherical-to-deformed shape phase-transition temperature of mean-field theory.

Crossover from vibrational to rotational collectivity in heavy nuclei in the shell-model Monte Carlo approach

Heavy nuclei exhibit a crossover from vibrational to rotational collectivity as the number of neutrons or protons increases from shell closure towards midshell, but the microscopic description of this crossover has been a major challenge. We apply the shell model Monte Carlo approach to families of even-even samarium and neodymium isotopes and identify a microscopic signature of the crossover from vibrational to rotational collectivity in the low-temperature behavior of ⟨J(2)⟩(T), where J is the total spin and T is the temperature. This signature agrees well with its values extracted from experimental data. We also calculate the state densities of these nuclei and find them to be in very good agreement with experimental data. Finally, we define a collective enhancement factor from the ratio of the total state density to the intrinsic state density as calculated in the finite-temperature Hartree-Fock-Bogoliubov approximation. The decay of this enhancement factor with excitation energy is found to correlate with the pairing and shape phase transitions in these nuclei.

Odd-particle systems in the shell model Monte Carlo method: circumventing a sign problem

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.

Pairing reentrance phenomenon in heated rotating nuclei in the shell-model Monte Carlo approach

Rotational motion of heated 72Ge is studied within the microscopic shell-model Monte Carlo approach. We investigate the angular momentum alignment and nuclear pairing correlations associated with J^{π} Cooper pairs as a function of the rotational frequency and temperature. The reentrance of pairing correlations with temperature is predicted at high rotational frequencies. It manifests itself through the anomalous behavior of specific heat and level density.

Heavy deformed nuclei in the shell model Monte Carlo method

We extend the shell model Monte Carlo approach to heavy deformed nuclei using a new proton-neutron formalism. The low excitation energies of such nuclei necessitate low-temperature calculations, for which a stabilization method is implemented in the canonical ensemble. We apply the method to study a well-deformed rare-earth nucleus, 162Dy. The single-particle model space includes the 50-82 shell plus 1f_{7/2} orbital for protons and the 82-126 shell plus 0h_{11/2}, 1g_{9/2} orbitals for neutrons. We show that the spherical shell model reproduces well the rotational character of 162Dy within this model space. We also calculate the level density of 162Dy and find it to be in excellent agreement with the experimental level density, which we extract from several experiments.

Spin projection in the shell model monte carlo method and the spin distribution of nuclear level densities

We introduce spin projection methods in the shell model Monte Carlo approach and apply them to calculate the spin distribution of level densities for iron-region nuclei using the complete (pf + g9/2) shell. We compare the calculated distributions with the spin-cutoff model and extract an energy-dependent moment of inertia. For even-even nuclei and at low excitation energies, we observe a significant suppression of the moment of inertia and odd-even staggering in the spin dependence of level densities.

Variational grand-canonical electronic structure method for open systems

An ab initio method is developed for variational grand-canonical molecular electronic structure of open systems based on the Gibbs-Peierls-Boguliobov inequality. We describe the theory and a practical method for performing the calculations within standard quantum chemistry codes using Gaussian basis sets. The computational effort scales similarly to the ground-state Hartree-Fock method. The quality of the approximation is studied on a hydrogen molecule by comparing to the exact Gibbs free energy, computed using full configuration-interaction calculations. We find the approximation quite accurate, with errors similar to those of the Hartree-Fock method for ground-state (zero-temperature) calculations. A further demonstration is given of the temperature effects on the bending potential curve for water. Some future directions and applications of the method are discussed. Several appendices give the mathematical and algorithmic details of the method.

The well-tempered auxiliary-field Monte Carlo

The auxiliary-field Monte Carlo (AFMC) is a method for computing ground-state and excited-state energies and other properties of electrons in molecules. For a given basis set, AFMC is an approximation to full-configuration interaction and the accuracy is determined predominantly by an inverse temperature "beta" parameter. A considerable amount of the dynamical correlation energy is recovered even at small values of beta. Yet, nondynamical correlation energy is inefficiently treated by AFMC. This is because the statistical error grows with beta, warranting increasing amount of Monte Carlo sampling. A recently introduced multi-determinant variant of AFMC is studied, and the method can be tuned by balancing the sizes of the determinantal space and the beta-parameter with respect to a predefined target accuracy. The well-tempered AFMC is considerably more efficient than a naive AFMC. As a welcome "byproduct" low lying excitation energies of the molecule are supplied as well. We demonstrate the principles on dissociating hydrogen molecule and torsion of ethylene where we calculate the (unoptimized) torsional barrier and the vertical singlet-triplet splitting.

Behavior of shell effects with the excitation energy in atomic nuclei

We study the behavior of shell effects, like pairing correlations and shape deformations, with the excitation energy in atomic nuclei. The analysis is carried out with the finite temperature Hartree-Fock-Bogoliubov method and a finite range density dependent force. For the first time, properties associated with the octupole and hexadecupole deformation and with the superdeformation as a function of the excitation energy are studied. Calculations for the well quadrupole deformed 164Er and 162Dy, superdeformed 152Dy, octupole deformed 224Ra, and spherical 118Sn nuclei are shown. We find, in particular, the level density of superdeformed states to be 4 orders of magnitude smaller than for normal deformed ones.

Parity dependence of nuclear level densities

A simple formula for the ratio of the number of odd- and even-parity states as a function of temperature is derived. This formula is used to calculate the ratio of level densities of opposite parities as a function of excitation energy. We test the formula with quantum Monte Carlo shell model calculations in the (pf+g(9/2)) shell. The formula describes well the transition from low excitation energies where a single parity dominates to high excitations where the two densities are equal.