-
Spectrum of density, spin and pairing fluctuations of an attractive two-dimensional Fermi gas
Authors:
Christian Apostoli,
Patrick Kelly,
Annette Lopez,
Kaelyn Dauer,
Gianluca Bertaina,
Davide Emilio Galli,
Ettore Vitali
Abstract:
We leverage random phase approximation and unbiased auxiliary-field quantum Monte Carlo methods to compute dynamical correlations for a dilute homogeneous two-dimensional attractive Fermi gas. Our main purpose is to quantitatively study the collective excitations of the system to generate robust benchmark results and to shed light into fermionic superfluidity in the strongly correlated regime. In…
▽ More
We leverage random phase approximation and unbiased auxiliary-field quantum Monte Carlo methods to compute dynamical correlations for a dilute homogeneous two-dimensional attractive Fermi gas. Our main purpose is to quantitatively study the collective excitations of the system to generate robust benchmark results and to shed light into fermionic superfluidity in the strongly correlated regime. In particular we are motivated by a recent paper suggesting that the Higgs mode can be detected in the spectrum of spin fluctuations. Despite the fact that we are somewhat limited by finite-size effects, our study pinpoints indeed a clear peak in the spin channel at low momentum, but a detailed analysis suggests that such a peak, though certainly interesting, does not correspond to the Higgs mode. We propose a different explanation for the shape of the spin structure factor. On the other hand, our results clearly show that the Higgs mode can be detected in the density channel at very small wave vectors, although very good resolution is necessary.
△ Less
Submitted 5 July, 2024; v1 submitted 20 March, 2024;
originally announced March 2024.
-
Quantum Monte Carlo and perturbative study of two-dimensional Bose-Fermi mixtures
Authors:
Jacopo D'Alberto,
Lorenzo Cardarelli,
Davide Emilio Galli,
Gianluca Bertaina,
Pierbiagio Pieri
Abstract:
We derive analytically the leading beyond-mean field contributions to the zero-temperature equation of state and to the fermionic quasi-particle residue and effective mass of a dilute Bose-Fermi mixture in two dimensions. In the repulsive case, we perform quantum Monte Carlo simulations for two representative bosonic concentrations and equal masses, extending a method for correcting finite-size ef…
▽ More
We derive analytically the leading beyond-mean field contributions to the zero-temperature equation of state and to the fermionic quasi-particle residue and effective mass of a dilute Bose-Fermi mixture in two dimensions. In the repulsive case, we perform quantum Monte Carlo simulations for two representative bosonic concentrations and equal masses, extending a method for correcting finite-size effects in fermionic gases to Bose-Fermi mixtures. We find good agreement between analytic expressions and numerical results for weak interactions, while significant discrepancies appear in the regime close to mechanical instability, above which we provide evidence of phase separation of the bosonic component.
△ Less
Submitted 8 May, 2024; v1 submitted 21 February, 2024;
originally announced February 2024.
-
Analysis of spin-squeezing generation in cavity-coupled atomic ensembles with continuous measurements
Authors:
A. Caprotti,
M. Barbiero,
M. G. Tarallo,
M. G. Genoni,
G. Bertaina
Abstract:
We analyze the generation of spin-squeezed states via coupling of three-level atoms to an optical cavity and continuous quantum measurement of the transmitted cavity field in order to monitor the evolution of the atomic ensemble. Using analytical treatment and microscopic simulations of the dynamics, we show that one can achieve significant spin squeezing, favorably scaling with the number of atom…
▽ More
We analyze the generation of spin-squeezed states via coupling of three-level atoms to an optical cavity and continuous quantum measurement of the transmitted cavity field in order to monitor the evolution of the atomic ensemble. Using analytical treatment and microscopic simulations of the dynamics, we show that one can achieve significant spin squeezing, favorably scaling with the number of atoms $N$. However, contrary to some previous literature, we clarify that it is not possible to obtain Heisenberg scaling without the continuous feedback that is proposed in optimal approaches. In fact, in the adiabatic cavity removal approximation and large $N$ limit, we find the scaling behavior $N^{-2/3}$ for spin squeezing and $N^{-1/3}$ for the corresponding protocol duration. These results can be obtained only by considering the curvature of the Bloch sphere, since linearizing the collective spin operators tangentially to its equator yields inaccurate predictions. With full simulations, we characterize how spin-squeezing generation depends on the system parameters and departs from the bad cavity regime, by gradually mixing with cavity-filling dynamics until metrological advantage is lost. Finally, we discuss the relevance of this spin-squeezing protocol to state-of-the-art optical clocks.
△ Less
Submitted 29 May, 2024; v1 submitted 27 November, 2023;
originally announced November 2023.
-
Phase Noise in Real-World Twin-Field Quantum Key Distribution
Authors:
Gianluca Bertaina,
Cecilia Clivati,
Simone Donadello,
Carlo Liorni,
Alice Meda,
Salvatore Virzì,
Marco Gramegna,
Marco Genovese,
Filippo Levi,
Davide Calonico,
Massimiliano Dispenza,
Ivo Pietro Degiovanni
Abstract:
The impact of noise sources in real-world implementations of Twin-Field Quantum Key Distribution (TF-QKD) protocols is investigated, focusing on phase noise from photon sources and connecting fibers. This work emphasizes the role of laser quality, network topology, fiber length, arm balance, and detector performance in determining key rates. Remarkably, it reveals that the leading TF-QKD protocols…
▽ More
The impact of noise sources in real-world implementations of Twin-Field Quantum Key Distribution (TF-QKD) protocols is investigated, focusing on phase noise from photon sources and connecting fibers. This work emphasizes the role of laser quality, network topology, fiber length, arm balance, and detector performance in determining key rates. Remarkably, it reveals that the leading TF-QKD protocols are similarly affected by phase noise despite different mechanisms. This study demonstrates duty cycle improvements of over a factor of two through narrow-linewidth lasers and phase-control techniques, highlighting the potential synergy with high-precision time/frequency distribution services. Ultrastable lasers, evolving toward integration and miniaturization, offer promise for agile TF-QKD implementations on existing networks. Properly addressing phase noise and practical constraints allows for consistent key rate predictions, protocol selection, and layout design, crucial for establishing secure long-haul links for the Quantum Communication Infrastructures under development in several countries.
△ Less
Submitted 29 May, 2024; v1 submitted 12 October, 2023;
originally announced October 2023.
-
Quantum Monte Carlo study of the role of p-wave interactions in ultracold repulsive Fermi gases
Authors:
Gianluca Bertaina,
Marco G. Tarallo,
Sebastiano Pilati
Abstract:
Single-component ultracold atomic Fermi gases are usually described using noninteracting many-fermion models. However, recent experiments reached a regime where $p$-wave interactions among identical fermionic atoms are important. In this paper, we employ variational and fixed-node diffusion Monte Carlo simulations to investigate the ground-state properties of single-component Fermi gases with shor…
▽ More
Single-component ultracold atomic Fermi gases are usually described using noninteracting many-fermion models. However, recent experiments reached a regime where $p$-wave interactions among identical fermionic atoms are important. In this paper, we employ variational and fixed-node diffusion Monte Carlo simulations to investigate the ground-state properties of single-component Fermi gases with short-range repulsive interactions. We determine the zero-temperature equation of state, and elucidate the roles played by the $p$-wave scattering volume and the $p$-wave effective range. A comparison against recently derived second-order perturbative results shows good agreement in a broad range of interaction strength. We also compute the quasiparticle effective mass, and we confirm the perturbative prediction of a linear contribution in the $p$-wave scattering volume, while we find significant deviations from the beyond-mean-field perturbative result, already for moderate interaction strengths. Finally, we determine ground-state energies for two-component unpolarized Fermi gases with both interspecies and intraspecies hard-sphere interactions, finding remarkable agreement with a recently derived fourth-order expansion that includes $p$-wave contributions.
△ Less
Submitted 8 May, 2023; v1 submitted 18 December, 2022;
originally announced December 2022.
-
Toward precision Fermi liquid theory in two dimensions
Authors:
Silas R. Beane,
Gianluca Bertaina,
Roland C. Farrell,
William R. Marshall
Abstract:
The ultra-cold and weakly-coupled Fermi gas in two spatial dimensions is studied in an effective field theory framework. It has long been observed that universal corrections to the energy density to two orders in the interaction strength do not agree with Monte Carlo simulations in the weak-coupling regime. Here, universal corrections to three orders in the interaction strength are obtained for th…
▽ More
The ultra-cold and weakly-coupled Fermi gas in two spatial dimensions is studied in an effective field theory framework. It has long been observed that universal corrections to the energy density to two orders in the interaction strength do not agree with Monte Carlo simulations in the weak-coupling regime. Here, universal corrections to three orders in the interaction strength are obtained for the first time, and are shown to provide agreement between theory and simulation. Special consideration is given to the scale ambiguity associated with the non-trivial renormalization of the singular contact interactions. The isotropic superfluid gap is obtained to next-to-leading order, and nonuniversal contributions to the energy density due to effective range effects, p-wave interactions and three-body forces are computed. Results are compared with precise Monte Carlo simulations of the energy density and the contact in the weakly-coupled attractive and repulsive Fermi liquid regimes. In addition, the known all-orders sum of ladder and ring diagrams is compared with Monte Carlo simulations at weak coupling and beyond.
△ Less
Submitted 14 April, 2023; v1 submitted 9 December, 2022;
originally announced December 2022.
-
Evolution of static and dynamical density correlations in a one-dimensional soft-core gas from the Tonks-Girardeau limit to a clustering fluid
Authors:
Martina Teruzzi,
Christian Apostoli,
Davide Pini,
Davide Emilio Galli,
Gianluca Bertaina
Abstract:
Repulsive soft-core atomic systems may undergo clustering if their density is high enough that core overlap is unavoidable. In one-dimensional quantum systems, it has been shown that this instability triggers a transition from a Luttinger liquid to various cluster Luttinger liquids. Here, we focus on the Luttinger liquid regime and theoretically study the evolution of key observables related to de…
▽ More
Repulsive soft-core atomic systems may undergo clustering if their density is high enough that core overlap is unavoidable. In one-dimensional quantum systems, it has been shown that this instability triggers a transition from a Luttinger liquid to various cluster Luttinger liquids. Here, we focus on the Luttinger liquid regime and theoretically study the evolution of key observables related to density fluctuations, that manifest a striking dependence on density. We tune the interaction so that the low-density regime corresponds to a Tonks-Girardeau gas, and show that as the density is increased the system departs more and more from Tonks-Girardeau behavior, displaying a much larger compressibility as well as rotonic excitations that finally drive the clustering transition. We compare various theoretical approaches, which are accurate in different regimes. Using quantum Monte Carlo methods and analytic continuation as a benchmark, we investigate the regime of validity of the mean-field Bogoliubov and the real-time multiconfiguration time-dependent Hartree-Fock approaches. Part of the behavior that we describe should be observable in ultracold Rydberg-dressed gases, provided that system losses are prevented.
△ Less
Submitted 22 June, 2021;
originally announced June 2021.
-
Quantum Monte Carlo simulations of two-dimensional repulsive Fermi gases with population imbalance
Authors:
S. Pilati,
G. Orso,
G. Bertaina
Abstract:
The ground-state properties of two-component repulsive Fermi gases in two dimensions are investigated by means of fixed-node diffusion Monte Carlo simulations. The energy per particle is determined as a function of the intercomponent interaction strength and of the population imbalance. The regime of universality in terms of the s-wave scattering length is identified by comparing results for hard-…
▽ More
The ground-state properties of two-component repulsive Fermi gases in two dimensions are investigated by means of fixed-node diffusion Monte Carlo simulations. The energy per particle is determined as a function of the intercomponent interaction strength and of the population imbalance. The regime of universality in terms of the s-wave scattering length is identified by comparing results for hard-disk and for soft-disk potentials. In the large imbalance regime, the equation of state turns out to be well described by a Landau-Pomeranchuk functional for two-dimensional polarons. To fully characterize this expansion, we determine the polarons' effective mass and their coupling parameter, complementing previous studies on their chemical potential. Furthermore, we extract the magnetic susceptibility from low-imbalance data, finding only small deviations from the mean-field prediction. While the mean-field theory predicts a direct transition from a paramagnetic to a fully ferromagnetic phase, our diffusion Monte Carlo results suggest that the partially ferromagnetic phase is stable in a narrow interval of the interaction parameter. This finding calls for further analyses on the effects due to the fixed-node constraint.
△ Less
Submitted 24 March, 2021;
originally announced March 2021.
-
Representing Molecular Ground and Excited Vibrational Eigenstates with Nuclear Densities obtained from Semiclassical Initial Value Representation Molecular Dynamics
Authors:
Chiara Aieta,
Gianluca Bertaina,
Marco Micciarelli,
Michele Ceotto
Abstract:
We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in protonated glycine molecule [C. Aieta et. al. Nat. Commun. 11, 4348 (2020)]. Here, we first validate and di…
▽ More
We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in protonated glycine molecule [C. Aieta et. al. Nat. Commun. 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, and to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond-lengths, angles and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching the typical harmonic nodal pattern is absent in the anharmonic distribution.
△ Less
Submitted 7 December, 2020;
originally announced December 2020.
-
Emergence of an Ising critical regime in the clustering of 1D soft matter revealed through string variables
Authors:
Francesco Mambretti,
Simone Molinelli,
Davide Pini,
Gianluca Bertaina,
Davide Emilio Galli
Abstract:
Soft matter systems are renowned for being able to display complex emerging phenomena such as clustering phases. Recently, a surprising quantum phase transition has been revealed in a one-dimensional (1D) system composed of bosons interacting via a pairwise soft potential in the continuum. It was shown that the spatial coordinates undergoing two-particle clustering could be mapped into quantum spi…
▽ More
Soft matter systems are renowned for being able to display complex emerging phenomena such as clustering phases. Recently, a surprising quantum phase transition has been revealed in a one-dimensional (1D) system composed of bosons interacting via a pairwise soft potential in the continuum. It was shown that the spatial coordinates undergoing two-particle clustering could be mapped into quantum spin variables of a 1D transverse Ising model. In this work we investigate the manifestation of an analogous critical phenomenon in 1D classical fluids of soft particles in the continuum. In particular, we study the low-temperature behavior of three different classical models of 1D soft matter, whose inter-particle interactions allow for clustering. The same string variables highlight that, at the commensurate density for the two-particle cluster phase, the peculiar pairing of neighboring soft particles can be nontrivially mapped onto a 1D discrete classical Ising model. We also observe a related phenomenon, namely the presence of an anomalous peak in the low-temperature specific heat, thus indicating the emergence of Schottky phenomenology in a non-magnetic fluid.
△ Less
Submitted 14 October, 2020; v1 submitted 31 July, 2020;
originally announced July 2020.
-
Dynamical structure factor of a fermionic supersolid on an optical lattice
Authors:
Ettore Vitali,
Patrick Kelly,
Annette Lopez,
Gianluca Bertaina,
Davide Emilio Galli
Abstract:
Interfacing unbiased quantum Monte Carlo simulations with state-of-art analytic continuation techniques, we obtain exact numerical results for dynamical density and spin correlations in the attractive Hubbard model, describing a spin-balanced two-dimensional cold Fermi gas on an optical lattice. We focus on half-filling, where on average one fermion occupies each lattice site, and the system displ…
▽ More
Interfacing unbiased quantum Monte Carlo simulations with state-of-art analytic continuation techniques, we obtain exact numerical results for dynamical density and spin correlations in the attractive Hubbard model, describing a spin-balanced two-dimensional cold Fermi gas on an optical lattice. We focus on half-filling, where on average one fermion occupies each lattice site, and the system displays an intriguing supersolid phase: a superfluid with a checkerboard density modulation. The coexistence of $U(1)$ broken symmetry and the density modulations makes this regime very challenging and interesting for the calculation of dynamical properties. We compare our unbiased results with state-of-the-art Generalized Random Phase Approximation calculations: both approaches agree on a well-defined low-energy Nambu-Goldstone collective mode in the density correlations, while the higher energy structures appear to differ significantly. We also observe an interesting high-energy spin mode. We argue that our results provide a robust benchmark for Generalized Random Phase Approximation techniques, which are widely considered to be the method of choice for dynamical correlations in Fermi gases. Also, our calculations yield new physical insight in the high-energy behavior of the dynamical structure factor of the attractive Hubbard model, which is a well known prototype lattice model for superconductors and is a fertile field to target the observation of collective modes in strongly correlated systems.
△ Less
Submitted 4 December, 2020; v1 submitted 12 December, 2019;
originally announced December 2019.
-
Reduced Rovibrational Coupling Cartesian Dynamics for Semiclassical Calculations: Application to the Spectrum of the Zundel Cation
Authors:
Gianluca Bertaina,
Giovanni Di Liberto,
Michele Ceotto
Abstract:
We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from [Huang et al., J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesia…
▽ More
We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from [Huang et al., J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesian dynamics and carefully reduce the coupling to global rotations in the definition of normal modes. We apply the time-averaging filter and obtain clean power spectra relative to suitable reference states, that highlight the spectral peaks corresponding to the fundamental excitations of the system. Our trajectory-based approach allows us for physical interpretation of the very challenging proton transfer modes. We find that it is important, for such a floppy molecule, to selectively avoid to initially excite lower energy modes, in order to obtain cleaner spectra. The estimated vibrational energies display a mean absolute error (MAE) of $\sim29\text{cm}^{-1}$ with respect to available Multi-Configuration time-dependent Hartree calculations and $\text{MAE}\sim14\text{cm}^{-1}$ when compared to the optically active experimental excitations of the Ne-tagged Zundel cation. The reasonable scaling in the number of trajectories for Monte Carlo convergence is promising for applications to higher dimensional protonated cluster systems.
△ Less
Submitted 29 July, 2019;
originally announced July 2019.
-
Two-dimensional repulsive Fermi polarons with short and long-range interactions
Authors:
Raúl Bombín,
Tommaso Comparin,
Gianluca Bertaina,
Ferran Mazzanti,
Stefano Giorgini,
Jordi Boronat
Abstract:
We study the repulsive polaron problem in a two-component two-dimensional system of fermionic atoms. We use two different interaction models: a short-range (hard-disk) potential and a dipolar potential. In our approach, all the atoms have the same mass and we consider the system to be composed of a uniform bath of a single species and a single atomic impurity. We use the diffusion Monte Carlo meth…
▽ More
We study the repulsive polaron problem in a two-component two-dimensional system of fermionic atoms. We use two different interaction models: a short-range (hard-disk) potential and a dipolar potential. In our approach, all the atoms have the same mass and we consider the system to be composed of a uniform bath of a single species and a single atomic impurity. We use the diffusion Monte Carlo method to evaluate polaron properties such as its chemical potential and pair distribution functions, together with a discussion on the deficit of volume induced by the impurity. We also evaluate observables that allow us to determine the validity of the quasi-particle picture: the quasi-particle residue and the effective mass of the polaron. Employing two different potentials allows us to identify the universality regime, where the properties depend only on the gas parameter $n a_s^2$ fixed by the bath density and the two-dimensional scattering length.
△ Less
Submitted 18 July, 2019; v1 submitted 24 May, 2019;
originally announced May 2019.
-
Static density response of one-dimensional soft bosons across the clustering transition
Authors:
Martina Teruzzi,
Davide Pini,
Stefano Rossotti,
Davide Emilio Galli,
Gianluca Bertaina
Abstract:
One-dimensional bosons interacting via a soft-shoulder potential are investigated at zero temperature. The flatness of the potential at short distances introduces a typical length, such that, at relatively high densities and sufficiently strong interactions, clusters are formed, even in the presence of a completely repulsive potential. We evaluate the static density response function of this syste…
▽ More
One-dimensional bosons interacting via a soft-shoulder potential are investigated at zero temperature. The flatness of the potential at short distances introduces a typical length, such that, at relatively high densities and sufficiently strong interactions, clusters are formed, even in the presence of a completely repulsive potential. We evaluate the static density response function of this system across the transition from the liquid to the cluster liquid phases. Such quantity reveals the density modulations induced by a weak periodic external potential, and is maximal at the clustering wavevector. It is known that this response function is proportional to the static structure factor in the classical regime at high temperature, while for this zero-temperature quantum systems, we extract it from the dynamical structure factor evaluated with quantum Monte Carlo methods.
△ Less
Submitted 9 April, 2018;
originally announced April 2018.
-
Quantum Critical Behavior of One-Dimensional Soft Bosons in the Continuum
Authors:
Stefano Rossotti,
Martina Teruzzi,
Davide Pini,
Davide Emilio Galli,
Gianluca Bertaina
Abstract:
We consider a zero-temperature one-dimensional system of bosons interacting via the soft-shoulder potential in the continuum, typical of dressed Rydberg gases. We employ quantum Monte Carlo simulations, which allow for the exact calculation of imaginary-time correlations, and a stochastic analytic continuation method, to extract the dynamical structure factor. At finite densities, in the weakly-in…
▽ More
We consider a zero-temperature one-dimensional system of bosons interacting via the soft-shoulder potential in the continuum, typical of dressed Rydberg gases. We employ quantum Monte Carlo simulations, which allow for the exact calculation of imaginary-time correlations, and a stochastic analytic continuation method, to extract the dynamical structure factor. At finite densities, in the weakly-interacting homogeneous regime, a rotonic spectrum marks the tendency to clustering. With strong interactions, we indeed observe cluster liquid phases emerging, characterized by the spectrum of a composite harmonic chain. Luttinger theory has to be adapted by changing the reference lattice density field. In both the liquid and cluster liquid phases, we find convincing evidence of a secondary mode, which becomes gapless only at the transition. In that region, we also measure the central charge and observe its increase towards c = 3/2, as recently evaluated in a related extended Bose-Hubbard model, and we note a fast reduction of the Luttinger parameter. For 2-particle clusters, we then interpret such observations in terms of the compresence of a Luttinger liquid and a critical transverse Ising model, related to the instability of the reference lattice density field towards coalescence of sites, typical of potentials which are flat at short distances. Even in the absence of a true lattice, we are able to evaluate the spatial correlation function of a suitable pseudo-spin operator, which manifests ferromagnetic order in the cluster liquid phase, exponential decay in the liquid phase, and algebraic order at criticality.
△ Less
Submitted 25 November, 2017; v1 submitted 27 July, 2017;
originally announced July 2017.
-
Statistical and computational intelligence approach to analytic continuation in Quantum Monte Carlo
Authors:
G. Bertaina,
D. E. Galli,
E. Vitali
Abstract:
The term analytic continuation emerges in many branches of Mathematics, Physics, and, more generally, applied Science. Generally speaking, in many situations, given some amount of information that could arise from experimental or numerical measurements, one is interested in extending the domain of such information, to infer the values of some variables which are central for the study of a given pr…
▽ More
The term analytic continuation emerges in many branches of Mathematics, Physics, and, more generally, applied Science. Generally speaking, in many situations, given some amount of information that could arise from experimental or numerical measurements, one is interested in extending the domain of such information, to infer the values of some variables which are central for the study of a given problem. For example, focusing on Condensed Matter Physics, state-of-the-art methodologies to study strongly correlated quantum physical systems are able to yield accurate estimations of dynamical correlations in imaginary time. Those functions have to be extended to the whole complex plane, via analytic continuation, in order to infer real-time properties of those physical systems. In this Review, we will present the Genetic Inversion via Falsification of Theories method, which allowed us to compute dynamical properties of strongly interacting quantum many--body systems with very high accuracy. Even though the method arose in the realm of Condensed Matter Physics, it provides a very general framework to face analytic continuation problems that could emerge in several areas of applied Science. Here we provide a pedagogical review that elucidates the approach we have developed.
△ Less
Submitted 6 March, 2017; v1 submitted 25 November, 2016;
originally announced November 2016.
-
Dynamical structure factor of one-dimensional hard rods
Authors:
M. Motta,
E. Vitali,
M. Rossi,
D. E. Galli,
G. Bertaina
Abstract:
The zero-temperature dynamical structure factor $S(q,ω)$ of one-dimensional hard rods is computed using state-of-the-art quantum Monte Carlo and analytic continuation techniques, complemented by a Bethe Ansatz analysis. As the density increases, $S(q,ω)$ reveals a crossover from the Tonks-Girardeau gas to a quasi-solid regime, along which the low-energy properties are found in agreement with the n…
▽ More
The zero-temperature dynamical structure factor $S(q,ω)$ of one-dimensional hard rods is computed using state-of-the-art quantum Monte Carlo and analytic continuation techniques, complemented by a Bethe Ansatz analysis. As the density increases, $S(q,ω)$ reveals a crossover from the Tonks-Girardeau gas to a quasi-solid regime, along which the low-energy properties are found in agreement with the nonlinear Luttinger liquid theory. Our quantitative estimate of $S(q,ω)$ extends beyond the low-energy limit and confirms a theoretical prediction regarding the behavior of $S(q,ω)$ at specific wavevectors $\mathcal{Q}_n=n 2 π/a$, where $a$ is the core radius, resulting from the interplay of the particle-hole boundaries of suitably rescaled ideal Fermi gases. We observe significant similarities between hard rods and one-dimensional $^4$He at high density, suggesting that the hard-rods model may provide an accurate description of dense one-dimensional liquids of quantum particles interacting through a strongly repulsive, finite-range potential.
△ Less
Submitted 27 August, 2016;
originally announced August 2016.
-
Microscopic study of static and dynamical properties of dilute one-dimensional soft bosons
Authors:
M. Teruzzi,
D. E. Galli,
G. Bertaina
Abstract:
We study static properties and the dynamical structure factor of zero-temperature dilute bosons interacting via a soft-shoulder potential in one dimension. Our approach is fully microscopic and employs state-of-the-art quantum Monte Carlo and analytic continuation techniques. By increasing the interaction strength, our model reproduces the Lieb-Liniger gas, the Tonks-Girardeau and the Hard-Rods mo…
▽ More
We study static properties and the dynamical structure factor of zero-temperature dilute bosons interacting via a soft-shoulder potential in one dimension. Our approach is fully microscopic and employs state-of-the-art quantum Monte Carlo and analytic continuation techniques. By increasing the interaction strength, our model reproduces the Lieb-Liniger gas, the Tonks-Girardeau and the Hard-Rods models.
△ Less
Submitted 11 January, 2017; v1 submitted 18 July, 2016;
originally announced July 2016.
-
Linear response of one-dimensional liquid $^4$He to external perturbations
Authors:
M. Motta,
G. Bertaina,
E. Vitali,
D. E. Galli,
M. Rossi
Abstract:
We study the response of one-dimensional liquid $^4$He to weak perturbations relying on the dynamical structure factor, $S(q,ω)$, recently obtained via ab-initio techniques [Phys. Rev. Lett. 116, 135302 (2016)]. We evaluate the drag force, $F_v$, experienced by an impurity moving along the system with velocity $v$ and the static response function, $χ(q)$, describing the density modulations induced…
▽ More
We study the response of one-dimensional liquid $^4$He to weak perturbations relying on the dynamical structure factor, $S(q,ω)$, recently obtained via ab-initio techniques [Phys. Rev. Lett. 116, 135302 (2016)]. We evaluate the drag force, $F_v$, experienced by an impurity moving along the system with velocity $v$ and the static response function, $χ(q)$, describing the density modulations induced by a periodic perturbation with wave vector $q$.
△ Less
Submitted 25 November, 2016; v1 submitted 18 July, 2016;
originally announced July 2016.
-
One-dimensional liquid $^4$He: dynamical properties beyond Luttinger liquid theory
Authors:
G. Bertaina,
M. Motta,
M. Rossi,
E. Vitali,
D. E. Galli
Abstract:
We compute the zero-temperature dynamical structure factor of one-dimensional liquid $^4$He by means of state-of-the-art Quantum Monte Carlo and analytic continuation techniques. By increasing the density, the dynamical structure factor reveals a transition from a highly compressible critical liquid to a quasi-solid regime. In the low-energy limit, the dynamical structure factor can be described b…
▽ More
We compute the zero-temperature dynamical structure factor of one-dimensional liquid $^4$He by means of state-of-the-art Quantum Monte Carlo and analytic continuation techniques. By increasing the density, the dynamical structure factor reveals a transition from a highly compressible critical liquid to a quasi-solid regime. In the low-energy limit, the dynamical structure factor can be described by the quantum hydrodynamic Luttinger liquid theory, with a Luttinger parameter spanning all possible values by increasing the density. At higher energies, our approach provides quantitative results beyond the Luttinger liquid theory. In particular, as the density increases, the interplay between dimensionality and interaction makes the dynamical structure factor manifest a pseudo {\it{particle-hole}} continuum typical of fermionic systems. At the low-energy boundary of such region and moderate densities, we find consistency, within statistical uncertainties, with predictions of a power-law structure by the recently-developed non-linear Luttinger liquid theory. In the quasi-solid regime we observe a novel behavior at intermediate momenta, which can be described by new analytical relations that we derive for the hard-rods model.
△ Less
Submitted 10 March, 2016; v1 submitted 22 December, 2014;
originally announced December 2014.
-
Condensed phase of Bose-Fermi mixtures with a pairing interaction
Authors:
Andrea Guidini,
Gianluca Bertaina,
Davide Emilio Galli,
Pierbiagio Pieri
Abstract:
We study the condensed phase of a Bose-Fermi mixture with a tunable pairing interaction between bosons and fermions with many-body diagrammatic methods and fixed-node diffusion Quantum Monte Carlo simulations. A universal behavior of the condensate fraction and bosonic momentum distribution with respect to the boson concentration is found to hold in an extended range of boson-fermion couplings and…
▽ More
We study the condensed phase of a Bose-Fermi mixture with a tunable pairing interaction between bosons and fermions with many-body diagrammatic methods and fixed-node diffusion Quantum Monte Carlo simulations. A universal behavior of the condensate fraction and bosonic momentum distribution with respect to the boson concentration is found to hold in an extended range of boson-fermion couplings and concentrations. For vanishing boson density, we prove that the bosonic condensate fraction reduces to the quasiparticle weight Z of the Fermi polaron studied in the context of polarized Fermi gases, unifying in this way two apparently unrelated quantities.
△ Less
Submitted 9 February, 2015; v1 submitted 8 December, 2014;
originally announced December 2014.
-
Energy of strongly attractive Bose-Fermi mixtures
Authors:
Andrea Guidini,
Elisa Fratini,
Gianluca Bertaina,
Pierbiagio Pieri
Abstract:
We discuss how approximate theories for Bose-Fermi mixtures recover in the molecular limit the expected expression for Fermi-Fermi mixtures of molecules and unpaired fermions. In particular, we compare the energy of the system resulting from a T-matrix diagrammatic approach with that obtained with Variational and Fixed-Node Diffusion Quantum Monte Carlo methods.
We discuss how approximate theories for Bose-Fermi mixtures recover in the molecular limit the expected expression for Fermi-Fermi mixtures of molecules and unpaired fermions. In particular, we compare the energy of the system resulting from a T-matrix diagrammatic approach with that obtained with Variational and Fixed-Node Diffusion Quantum Monte Carlo methods.
△ Less
Submitted 4 December, 2014;
originally announced December 2014.
-
Implementation of the Linear Method for the optimization of Jastrow-Feenberg and Backflow Correlations
Authors:
M. Motta,
G. Bertaina,
D. E. Galli,
E. Vitali
Abstract:
We present a fully detailed and highly performing implementation of the Linear Method [J. Toulouse and C. J. Umrigar (2007)] to optimize Jastrow-Feenberg and Backflow Correlations in many-body wave-functions, which are widely used in condensed matter physics. We show that it is possible to implement such optimization scheme performing analytical derivatives of the wave-function with respect to the…
▽ More
We present a fully detailed and highly performing implementation of the Linear Method [J. Toulouse and C. J. Umrigar (2007)] to optimize Jastrow-Feenberg and Backflow Correlations in many-body wave-functions, which are widely used in condensed matter physics. We show that it is possible to implement such optimization scheme performing analytical derivatives of the wave-function with respect to the variational parameters achieving the best possible complexity O(N^3) in the number of particles N.
△ Less
Submitted 10 December, 2014; v1 submitted 1 December, 2014;
originally announced December 2014.
-
Quantum Monte Carlo study of the indirect Pauli exclusion effect in Bose-Fermi mixtures
Authors:
G. Bertaina,
A. Guidini,
P. Pieri
Abstract:
We study the momentum distributions of a three-dimensional resonant Bose-Fermi mixture in the molecular limit at zero temperature. For concentration of the bosons with respect to the fermions less or equal to one, each boson is bound to a fermion and the system is composed of fermionic molecules plus excess fermions. Not only the bosonic condensate fraction goes to zero, signaling a quantum phase…
▽ More
We study the momentum distributions of a three-dimensional resonant Bose-Fermi mixture in the molecular limit at zero temperature. For concentration of the bosons with respect to the fermions less or equal to one, each boson is bound to a fermion and the system is composed of fermionic molecules plus excess fermions. Not only the bosonic condensate fraction goes to zero, signaling a quantum phase transition towards a normal phase, but a finite region of low momenta is depleted, depending on the concentration. This phenomenon is named indirect Pauli exclusion effect and is demonstrated via Fixed-Node Diffusion Monte Carlo simulations and T-matrix calculations.
△ Less
Submitted 12 May, 2015; v1 submitted 27 November, 2014;
originally announced November 2014.
-
Bose-Fermi mixtures in the molecular limit
Authors:
A. Guidini,
G. Bertaina,
E. Fratini,
P. Pieri
Abstract:
We consider a Bose-Fermi mixture in the molecular limit of the attractive interaction between fermions and bosons. For a boson density smaller or equal to the fermion density, we show analytically how a T-matrix approach for the constituent bosons and fermions recovers the expected physical limit of a Fermi-Fermi mixture of molecules and atoms. In this limit, we derive simple expressions for the s…
▽ More
We consider a Bose-Fermi mixture in the molecular limit of the attractive interaction between fermions and bosons. For a boson density smaller or equal to the fermion density, we show analytically how a T-matrix approach for the constituent bosons and fermions recovers the expected physical limit of a Fermi-Fermi mixture of molecules and atoms. In this limit, we derive simple expressions for the self-energies, the momentum distribution function, and the chemical potentials. By extending these equations to a trapped system, we determine how to tailor the experimental parameters of a Bose-Fermi mixture in order to enhance the 'indirect Pauli exclusion effect' on the boson momentum distribution function. For the homogeneous system, we present finally a Diffusion Monte Carlo simulation which confirms the occurrence of such a peculiar effect.
△ Less
Submitted 8 March, 2014; v1 submitted 3 January, 2014;
originally announced January 2014.
-
Two-dimensional short-range interacting attractive and repulsive Fermi gases at zero temperature
Authors:
Gianluca Bertaina
Abstract:
We study a two-dimensional two-component Fermi gas with attractive or repulsive short-range interactions at zero temperature. We use Diffusion Monte Carlo with Fixed Node approximation in order to calculate the energy per particle and the opposite-spin pair distribution functions. We show the relevance of beyond mean field effects and verify the consistency of our approach by using Tan's Contact r…
▽ More
We study a two-dimensional two-component Fermi gas with attractive or repulsive short-range interactions at zero temperature. We use Diffusion Monte Carlo with Fixed Node approximation in order to calculate the energy per particle and the opposite-spin pair distribution functions. We show the relevance of beyond mean field effects and verify the consistency of our approach by using Tan's Contact relations.
△ Less
Submitted 15 March, 2013;
originally announced March 2013.
-
Quantum Monte Carlo Study of a Resonant Bose-Fermi Mixture
Authors:
G. Bertaina,
E. Fratini,
S. Giorgini,
P. Pieri
Abstract:
We study a resonant Bose-Fermi mixture at zero temperature by using the fixed-node diffusion Monte Carlo method. We explore the system from weak to strong boson-fermion interaction, for different concentrations of the bosons relative to the fermion component. We focus on the case where the boson density $n_B$ is smaller than the fermion density $n_F$, for which a first-order quantum phase transiti…
▽ More
We study a resonant Bose-Fermi mixture at zero temperature by using the fixed-node diffusion Monte Carlo method. We explore the system from weak to strong boson-fermion interaction, for different concentrations of the bosons relative to the fermion component. We focus on the case where the boson density $n_B$ is smaller than the fermion density $n_F$, for which a first-order quantum phase transition is found from a state with condensed bosons immersed in a Fermi sea, to a Fermi-Fermi mixture of composite fermions and unpaired fermions. We obtain the equation of state and the phase diagram, and we find that the region of phase separation shrinks to zero for vanishing $n_B$.
△ Less
Submitted 15 March, 2013; v1 submitted 26 November, 2012;
originally announced November 2012.
-
BCS-BEC crossover in a two-dimensional Fermi gas
Authors:
G. Bertaina,
S. Giorgini
Abstract:
We investigate the crossover from Bardeen-Cooper-Schrieffer (BCS) superfluidity to Bose-Einstein condensation (BEC) in a two-dimensional Fermi gas at T=0 using the fixed-node diffusion Monte Carlo method. We calculate the equation of state and the gap parameter as a function of the interaction strength, observing large deviations compared to mean-field predictions. In the BEC regime our results sh…
▽ More
We investigate the crossover from Bardeen-Cooper-Schrieffer (BCS) superfluidity to Bose-Einstein condensation (BEC) in a two-dimensional Fermi gas at T=0 using the fixed-node diffusion Monte Carlo method. We calculate the equation of state and the gap parameter as a function of the interaction strength, observing large deviations compared to mean-field predictions. In the BEC regime our results show the important role of dimer-dimer and atom-dimer interaction effects that are completely neglected in the mean-field picture. Results on Tan's contact parameter associated with short-range physics are also reported along the BCS-BEC crossover.
△ Less
Submitted 21 March, 2011; v1 submitted 16 November, 2010;
originally announced November 2010.
-
Itinerant ferromagnetism of a repulsive atomic Fermi gas: a quantum Monte Carlo study
Authors:
S. Pilati,
G. Bertaina,
S. Giorgini,
M. Troyer
Abstract:
We investigate the phase diagram of a two-component repulsive Fermi gas at T=0 by means of quantum Monte Carlo simulations. For a given value of the positive s-wave scattering length, both purely repulsive and purely attractive model potentials are considered in order to analyze the limits of the universal regime where the details of interatomic forces can be neglected. The equation of state of bo…
▽ More
We investigate the phase diagram of a two-component repulsive Fermi gas at T=0 by means of quantum Monte Carlo simulations. For a given value of the positive s-wave scattering length, both purely repulsive and purely attractive model potentials are considered in order to analyze the limits of the universal regime where the details of interatomic forces can be neglected. The equation of state of both balanced and unbalanced systems is calculated as a function of the interaction strength and the critical density for the onset of ferromagnetism is determined. The energy per particle of the strongly polarized gas is calculated and parametrized in terms of the physical properties of repulsive polarons, which are relevant for the stability of the fully magnetized ferromagnetic state. Finally, we analyze the phase diagram in the polarization/interaction plane under the assumption that only phases with homogeneous magnetization can be produced.
△ Less
Submitted 7 April, 2010;
originally announced April 2010.
-
First and second sound in cylindrically trapped gases
Authors:
G. Bertaina,
L. Pitaevskii,
S. Stringari
Abstract:
We investigate the propagation of density and temperature waves in a cylindrically trapped gas with radial harmonic confinement. Starting from two-fluid hydrodynamic theory we derive effective 1D equations for the chemical potential and the temperature which explicitly account for the effects of viscosity and thermal conductivity. Differently from quantum fluids confined by rigid walls, the harmon…
▽ More
We investigate the propagation of density and temperature waves in a cylindrically trapped gas with radial harmonic confinement. Starting from two-fluid hydrodynamic theory we derive effective 1D equations for the chemical potential and the temperature which explicitly account for the effects of viscosity and thermal conductivity. Differently from quantum fluids confined by rigid walls, the harmonic confinement allows for the propagation of both first and second sound in the long wave length limit. We provide quantitative predictions for the two sound velocities of a superfluid Fermi gas at unitarity. For shorter wave-lengths we discover a new surprising class of excitations continuously spread over a finite interval of frequencies. This results in a non-dissipative damping in the response function which is analytically calculated in the limiting case of a classical ideal gas.
△ Less
Submitted 14 October, 2010; v1 submitted 1 February, 2010;
originally announced February 2010.
-
Density profiles of polarized Fermi gases confined in harmonic traps
Authors:
G. Bertaina,
S. Giorgini
Abstract:
On the basis of the phase diagram of the uniform system we calculate the density profiles of a trapped polarized Fermi gas at zero temperature using the local density approximation. By varying the overall polarization and the interaction strength we analyze the appearance of a discontinuity in the profile, signalling a first order phase transition from a superfluid inner core to a normal outer s…
▽ More
On the basis of the phase diagram of the uniform system we calculate the density profiles of a trapped polarized Fermi gas at zero temperature using the local density approximation. By varying the overall polarization and the interaction strength we analyze the appearance of a discontinuity in the profile, signalling a first order phase transition from a superfluid inner core to a normal outer shell. The local population imbalance between the two components and the size of the various regions of the cloud corresponding to different phases are also discussed. The calculated profiles are quantitatively compared with the ones recently measured by Shin {\it et al.}, Phys. Rev. Lett. {\bf 101}, 070404 (2008).
△ Less
Submitted 23 January, 2009; v1 submitted 12 November, 2008;
originally announced November 2008.