Researcher profile

Joaquín E. Drut

Joaquín E. Drut contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2021arXiv

Pairing patterns in polarized unitary Fermi gases above the superfluid transition

We non-perturbatively study pairing in the high-temperature regime of polarized unitary two-component Fermi gases by extracting the pair-momentum distribution and shot-noise correlations. Whereas the pair-momentum distribution allows us to analyze the propagation of pairs composed of one spin-up and one spin-down fermion, shot-noise correlations provide us with a tomographic insight into pairing correlations around the Fermi surfaces associated with the two species. Assuming that the dominant pairing patterns right above the superfluid transition also govern the formation of condensates in the low-temperature regime, our analysis suggests that the superfluid ground state is homogeneous and of the Bardeen-Cooper-Schrieffer-type over a wide range of polarizations.

preprint2020arXiv

Emergence of a pseudogap in the BCS-BEC crossover

Strongly correlated Fermi systems with pairing interactions become superfluid below a critical temperature $T_c$. The extent to which such pairing correlations alter the behavior of the liquid at temperatures $T > T_c$ is a subtle issue that remains an area of debate, in particular regarding the appearance of the so-called pseudogap in the BCS-BEC crossover of unpolarized spin-$1/2$ nonrelativistic matter. To shed light on this, we extract several quantities of crucial importance at and around the unitary limit, namely: the odd-even staggering of the total energy, the spin susceptibility, the pairing correlation function, the condensate fraction, and the critical temperature $T_c$, using a non-perturbative, constrained-ensemble quantum Monte Carlo algorithm.

preprint2020arXiv

Pairing patterns in one-dimensional spin- and mass-imbalanced Fermi gases

We study spin- and mass-imbalanced mixtures of spin-$\tfrac{1}{2}$ fermions interacting via an attractive contact potential in one spatial dimension. Specifically, we address the influence of unequal particle masses on the pair formation by means of the complex Langevin method. By computing the pair-correlation function and the associated pair-momentum distribution we find that inhomogeneous pairing is present for all studied spin polarizations and mass imbalances. To further characterize the pairing behavior, we analyze the density-density correlations in momentum space, the so-called shot noise, which is experimentally accessible through time-of-flight imaging. At finite spin polarization, the latter is known to show distinct maxima at momentum configurations associated with the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) instability. Besides those maxima, we find that additional features emerge in the noise correlations when mass imbalance is increased, revealing the stability of FFLO-type correlations against mass imbalance and furnishing an experimentally accessible signature to probe this type of pairing.

preprint2019arXiv

Thermodynamics of spin-orbit coupled bosons in two dimensions from complex Langevin

We investigate the thermal properties of interacting spin-orbit coupled bosons with contact interactions in two spatial dimensions. To that end, we implement the complex Langevin method, motivated by the appearance of a sign problem, on a square lattice with periodic boundary conditions. We calculate the density equation of state non-perturbatively in a range of spin-orbit couplings and chemical potentials. Our results show that mean-field solutions tend to underestimate the average density, especially for stronger values of the spin-orbit coupling. Additionally, the finite nature of the simulation volume induces the formation of pseudo-condensates. These have been observed to be destroyed by the spin-orbit interactions.

preprint2013arXiv

Cooper pairing above the critical temperature in a unitary Fermi gas

We present an ab initio determination of spin responses of the unitary Fermi gas, based on finite temperature quantum Monte Carlo calculations and the Kubo linear-response formalism. We determine the temperature dependence of the spin susceptibility and the spin conductivity. We show that both quantities exhibit suppression above the critical temperature of the superfluid to normal phase transition due to presence of the Cooper pairs. The spin diffusion transport coefficient does not display the existence of a minimum in the vicinity of the critical temperature and it drops to very low values D_s approx 0.8hbar/m in the superfluid phase. All these spin observables show a smooth and monotonic behavior with temperature when crossing the critical temperature T_c, until the Fermi liquid regime is attained at the temperature T*, where the pseudogap regime disappears.

preprint2013arXiv

Imaginary polarization as a way to surmount the sign problem in ab initio calculations of spin-imbalanced Fermi gases

From ultracold atoms to quantum chromodynamics, reliable ab initio studies of strongly interacting fermions require numerical methods, typically in some form of quantum Monte Carlo calculation. Unfortunately, (non)relativistic systems at finite density (spin polarization) generally have a sign problem, such that those ab initio calculations are impractical. It is well-known, however, that in the relativistic case imaginary chemical potentials solve this problem, assuming the data can be analytically continued to the real axis. Is this feasible for nonrelativistic systems? Are the interesting features of the phase diagram accessible in this manner? By introducing complex chemical potentials, for real total particle number and imaginary polarization, the sign problem is avoided in the nonrelativistic case. To give a first answer to the above questions, we perform a mean-field study of the finite-temperature phase diagram of spin-1/2 fermions with imaginary polarization.

preprint2013arXiv

Precise determination of the structure factor and contact in a unitary Fermi gas

We present a high-precision determination of the universal contact parameter in a strongly interacting Fermi gas. In a trapped gas at unitarity we find the contact to be $3.06 \pm 0.08$ at a temperature of 0.08 of the Fermi temperature in a harmonic trap. The contact governs the high-momentum (short-range) properties of these systems and this low temperature measurement provides a new benchmark for the zero temperature homogeneous contact. The experimental measurement utilises Bragg spectroscopy to obtain the dynamic and static structure factors of ultracold Fermi gases at high momentum in the unitarity and molecular Bose-Einstein condensate (BEC) regimes. We have also performed quantum Monte Carlo calculations of the static properties, extending from the weakly coupled Bardeen-Cooper-Schrieffer (BCS) regime to the strongly coupled BEC case, which show agreement with experiment at the level of a few percent.

preprint2012arXiv

Improved lattice operators for non-relativistic fermions

In this work I apply a recently proposed improvement procedure, originally conceived to reduce finite lattice spacing effects in transfer matrices for dilute Fermi systems, to tuning operators for the calculation of observables. I construct, in particular, highly improved representations for the energy and the contact, as a first step in an improvement program for finite-temperature calculations. I illustrate the effects of improvement on those quantities with a ground-state lattice calculation at unitarity.

preprint2012arXiv

Shear Viscosity of a Unitary Fermi Gas

We present the first ab initio determination of the shear viscosity eta of the Unitary Fermi Gas, based on finite temperature quantum Monte Carlo calculations and the Kubo linear-response formalism. We determine the temperature dependence of the shear viscosity to entropy density ratio eta/s. The minimum of eta/s appears to be located above the critical temperature for the superfluid-to-normal phase transition with the most probable value being eta/s approx 0.2 hbar/kB, which almost saturates the Kovtun-Son-Starinets universal value hbar/(4 pi kB).

preprint2012arXiv

The Equation of State of the Unitary Fermi Gas: An Update on Lattice Calculations

The thermodynamic properties of the unitary Fermi gas (UFG) have recently been measured to unprecedented accuracy at the MIT. In particular, these measurements provide an improved understanding of the regime below T/eF ~ 0.20, where a transition into a superfluid phase occurs. In light of this development, we present an overview of state-of-the-art auxiliary field quantum Monte Carlo (AFQMC) results for the UFG at finite temperature and compare them with the MIT data for the energy, chemical potential, and density. These AFQMC results have been obtained using methods based on the hybrid Monte Carlo (HMC) algorithm, which was first introduced within the context of lattice QCD.

preprint2011arXiv

Exact-exchange density functional theory for neutron drops

We compute the ground-state properties of finite systems of neutrons in an external harmonic trap, interacting via the Minnesota potential, using the "exact-exchange" form of orbital-dependent density functional theory. We compare our results with Hartree-Fock calculations and find very close agreement. Within the context of the interaction studied, we conclude that this simple orbital-dependent functional brings conventional nuclear density functional theory to the level of Hartree-Fock in an ab initio fashion. Our work is a first step towards higher-order ab initio nuclear functionals based on realistic nucleon-nucleon interactions.

preprint2011arXiv

Momentum Distribution and Contact of the Unitary Fermi gas

We calculate the momentum distribution n(k) of the Unitary Fermi Gas using Quantum Monte Carlo calculations at finite temperature T/ε_F as well as in the ground state. At large momenta k/k_F, we find that n(k) falls off as C/k^4, in agreement with the Tan relations. From the asymptotics of n(k), we determine the contact C as a function of T/ε_F and present a comparison with theory. At low T/ε_F, we find that C increases with temperature, and we tentatively identify a maximum around T/ε_F \simeq 0.4. Our calculations are performed on lattices of spatial extent up to N_x = 14 with a particle number per unit volume of \simeq 0.03 - 0.07.

preprint2011arXiv

Strongly Coupled Graphene on the Lattice

The two-dimensional carbon allotrope graphene has recently attracted a lot of attention from researchers in the disciplines of Lattice Field Theory, Lattice QCD and Monte Carlo calculations. This interest has been prompted by several remarkable properties of the conduction electrons in graphene. For instance, the conical band structure of graphene at low energies is strongly reminiscent of relativistic Dirac fermions. Also, due the low Fermi velocity of v_F = c/300, where c is the speed of light in vacuum, the physics of the conduction electrons in graphene is qualitatively similar to Quantum Electrodynamics in a strongly coupled regime. In turn, this opens up the prospect of the experimental realization of gapped, strongly correlated states in the electronic phase diagram of graphene. Here, we review the experimental and theoretical motivations for Lattice Field Theory studies of graphene, and describe the directions that such research is likely to progress in during the next few years. We also give a brief overview of the two main lattice theories of graphene, the hexagonal Hubbard theory and the low-energy Dirac theory. Finally, we describe the prospect of extracting response functions, such as the electric conductivity, using Lattice Field Theory calculations.

preprint2011arXiv

The unitary Fermi gas at finite temperature: momentum distribution and contact

The Unitary Fermi Gas (UFG) is one of the most strongly interacting systems known to date, as it saturates the unitarity bound on the quantum mechanical scattering cross section. The UFG corresponds to a two-component Fermi gas in the limit of short interaction range and large scattering length, and is currently realized in ultracold-atom experiments via Feshbach resonances. While easy to define, the UFG poses a challenging quantum many-body problem, as it lacks any characteristic scale other than the density. As a consequence, accurate quantitative predictions of the thermodynamic properties of the UFG require Monte Carlo calculations. However, significant progress has also been made with purely analytical methods. Notably, in 2005 Tan derived a set of exact thermodynamic relations in which a universal quantity known as the "contact" C plays a crucial role. Recently, C has also been found to determine the prefactor of the high- frequency power-law decay of correlators as well as the right-hand-sides of shear- and bulk viscosity sum rules. The contact is therefore a central piece of information on the UFG in equilibrium as well as away from equilibrium. In this talk we describe some of the known aspects of Fermi gases at and around unitarity, show our latest Monte Carlo results for the contact at finite temperature, and summarize the open questions in the field, some of which we are starting to answer using large-scale Monte Carlo calculations by adapting methods from Lattice QCD.

preprint2010arXiv

First-order chiral transition in the compact lattice theory of graphene and the case for improved actions

A comparison of the compact and non-compact lattice versions of the low-energy theory of graphene is presented. The compact theory is found to exhibit a chiral phase transition which appears to be of first order, at a critical coupling of beta_c = 0.42 +/- 0.01. We confirm that the non-compact theory exhibits a second-order transition at beta_c = 0.072 +/- 0.003, and determine the effects of UV-divergent tadpole contributions in both cases. Upon tadpole improvement of the non-compact theory we find beta_c^{TI} = 0.163 +/- 0.002, which strengthens the case for a semimetal-insulator transition in graphene at strong Coulomb coupling. Finally, we highlight the need for systematic studies using improved lattice actions.

preprint2010arXiv

Graphene: from materials science to particle physics

Since its discovery in 2004, graphene, a two-dimensional hexagonal carbon allotrope, has generated great interest and spurred research activity from materials science to particle physics and vice versa. In particular, graphene has been found to exhibit outstanding electronic and mechanical properties, as well as an unusual low-energy spectrum of Dirac quasiparticles giving rise to a fractional quantum Hall effect when freely suspended and immersed in a magnetic field. One of the most intriguing puzzles of graphene involves the low-temperature conductivity at zero density, a central issue in the design of graphene-based nanoelectronic components. While suspended graphene experiments have shown a trend reminiscent of semiconductors, with rising resistivity at low temperatures, most theories predict a constant or even decreasing resistivity. However, lattice field theory calculations have revealed that suspended graphene is at or near the critical coupling for excitonic gap formation due to strong Coulomb interactions, which suggests a simple and straightforward explanation for the experimental data. In this contribution we review the current status of the field with emphasis on the issue of gap formation, and outline recent progress and future points of contact between condensed matter physics and Lattice QCD.

preprint2010arXiv

Preconditioning the non-relativistic many-fermion problem

Preconditioning is at the core of modern many-fermion Monte Carlo algorithms, such as Hybrid Monte Carlo, where the repeated solution of a linear problem involving an ill-conditioned matrix is needed. We report on a performance comparison of three preconditioning strategies, namely Chebyshev polynomials, strong-coupling approximation and weak-coupling expansion. We use conjugate gradient (CG) on the normal equations as well as stabilized biconjugate gradient (BiCGStab) as solvers and focus on the fermion matrix of the unitary Fermi gas. Our results indicate that BiCGStab is by far the most efficient strategy, both in terms of the number of iterations and matrix-vector operations.