Source author record

C. N. Gilbreth

C. N. Gilbreth appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2020arXiv

Pairing correlations across the superfluid phase transition in the unitary Fermi gas

In the two-component Fermi gas with a contact interaction, a pseudogap regime can exist at temperatures between the superfluid critical temperature $T_c$ and a temperature $T^* > T_c$. This regime is characterized by pairing correlations without superfluidity. However, in the unitary limit of infinite scattering length, the existence of this regime is still debated. To help address this, we have applied finite-temperature auxiliary-field quantum Monte Carlo (AFMC) to study the thermodynamics of the superfluid phase transition and signatures of the pseudogap in the spin-balanced homogeneous unitary Fermi gas. We present results at finite filling factor $ν\simeq 0.06$ for the condensate fraction, an energy-staggering pairing gap, the spin susceptibility, and the heat capacity, and compare them to experimental data when available. In contrast to previous AFMC simulations, our model space consists of the complete first Brillouin zone of the lattice, and our calculations are performed in the canonical ensemble of fixed particle number. The canonical ensemble AFMC framework enables the calculation of a model-independent gap, providing direct information on pairing correlations without the need for numerical analytic continuation. We use finite-size scaling to estimate $T_c$ at the corresponding filling factor. We find that the energy-staggering pairing gap vanishes above $T_c$, showing no pseudogap effects, and that the spin susceptibility shows a substantially reduced signature of a spin gap compared to previously reported AFMC simulations.

preprint2019arXiv

The contact in the unitary Fermi gas across the superfluid phase transition

A quantity known as the contact plays a fundamental role in quantum many-body systems with short-range interactions. The determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our results are in excellent agreement with the most recent precision ultracold atomic gas experiments. We also present results for the energy of the unitary gas as a function of temperature in the continuum limit.

preprint2016arXiv

Microsopic nuclear level densities by the shell model Monte Carlo method

The configuration-interaction shell model approach provides an attractive framework for the calculation of nuclear level densities in the presence of correlations, but the large dimensionality of the model space has hindered its application in mid-mass and heavy nuclei. The shell model Monte Carlo (SMMC) method permits calculations in model spaces that are many orders of magnitude larger than spaces that can be treated by conventional diagonalization methods. We discuss recent progress in the SMMC approach to level densities, and in particular the calculation of level densities in heavy nuclei. We calculate the distribution of the axial quadrupole operator in the laboratory frame at finite temperature and demonstrate that it is a model-independent signature of deformation in the rotational invariant framework of the shell model. We propose a method to use these distributions for calculating level densities as a function of intrinsic deformation.

preprint2015arXiv

Benchmarking mean-field approximations to level densities

We assess the accuracy of finite-temperature mean-field theory using as a standard the Hamiltonian and model space of the shell model Monte Carlo calculations. Two examples are considered: the nucleus $^{162}$Dy, representing a heavy deformed nucleus, and $^{148}$Sm, representing a nearby heavy spherical nucleus with strong pairing correlations. The errors inherent in the finite-temperature Hartree-Fock and Hartree-Fock-Bogoliubov approximations are analyzed by comparing the entropies of the grand canonical and canonical ensembles, as well as the level density at the neutron resonance threshold, with shell model Monte Carlo (SMMC) calculations, which are accurate up to well-controlled statistical errors. The main weak points in the mean-field treatments are seen to be: (i) the extraction of number-projected densities from the grand canonical ensembles, and (ii) the symmetry breaking by deformation or by the pairing condensate. In the absence of a pairing condensate, we confirm that the usual saddle-point approximation to extract the number-projected densities is not a significant source of error compared to other errors inherent to the mean-field theory. We also present an alternative formulation of the saddle-point approximation that makes direct use of an approximate particle-number projection and avoids computing the usual three-dimensional Jacobian of the saddle-point integration. We find that the pairing condensate is less amenable to approximate particle-number projection methods due to the explicit violation of particle-number conservation in the pairing condensate. Nevertheless, the Hartree-Fock-Bogoliubov theory is accurate to less than one unit of entropy for $^{148}$Sm at the neutron threshold energy, which is above the pairing phase transition.

preprint2014arXiv

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 the 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 higher temperatures than the spherical-to-deformed shape phase-transition temperature of mean-field theory.

preprint2014arXiv

Stabilizing Canonical-Ensemble Calculations in the Auxiliary-Field Monte Carlo Method

Quantum Monte Carlo methods are powerful techniques for studying strongly interacting Fermi systems. However, implementing these methods on computers with finite-precision arithmetic requires careful attention to numerical stability. In the auxiliary-field Monte Carlo (AFMC) method, low-temperature or large-model-space calculations require numerically stabilized matrix multiplication. When adapting methods used in the grand-canonical ensemble to the canonical ensemble of fixed particle number, the numerical stabilization increases the number of required floating-point operations for computing observables by a factor of the size of the single-particle model space, and thus can greatly limit the systems that can be studied. We describe an improved method for stabilizing canonical-ensemble calculations in AFMC that exhibits better scaling, and present numerical tests that demonstrate the accuracy and improved performance of the method.

preprint2013arXiv

Pair condensation in a Finite Trapped Fermi Gas

Superfluidity in the cold atomic two-species Fermi gas system in the unitary limit of infinite scattering length remains incompletely understood. In particular, a pseudogap phase has been proposed to exist above the superfluid critical temperature. Here we apply the auxiliary-field quantum Monte Carlo method to perform the first ab initio calculations of the temperature dependence of three quantities -- the energy-staggering pairing gap, the condensate fraction and the heat capacity -- in a trapped finite-size cold atom system. As the calculations of the energy-staggering pairing gap require the use of the canonical ensemble, we employ a novel algorithm for the stabilization of particle-number projection that is essential for reaching convergence in the size of the model space. We observe clear signatures of the superfluid phase transition in all three quantities, including a signature of the recently measured lambda peak in the heat capacity, but find no evidence of a pseudogap effect in the energy-staggering pairing gap.

preprint2011arXiv

A new effective interaction for the trapped fermi gas: the BEC-BCS crossover

We extend a recently introduced separable interaction for the unitary trapped Fermi gas to all values of the scattering length. We derive closed expressions for the interaction matrix elements and the two-particle eigenvectors and analytically demonstrate the convergence of this interaction to the zero-range two-body pseudopotential for s-wave scattering. We apply this effective interaction to the three- and four-particle systems along the BEC-BCS crossover, and find that their low-lying energies exhibit convergence in the regularization parameter that is much faster than for the conventional renormalized contact interaction. We find similar convergence properties of the three-particle free energy at unitarity.