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

Diagrammatic Monte Carlo study of polaron systems

Jonas Vlietinck

Published Jan 1, 2015• 0 citations
OpenAlex

Abstract

We report on a study of a spin-down impurity strongly coupled to a spin-up Fermi sea (a so-called Fermi polaron) with the diagrammatic Monte-Carlo (DiagMC) technique.Conditions of zero temperature and three dimensions are considered for an ultracold atomic gas with resonant interactions in the zero-range limit.A Feynman diagrammatic series is developed for the one-body and two-body propagators providing information about the polaron and molecule channel respectively.The DiagMC technique allows us to reach diagram orders that are high enough for extrapolation to infinite order.The robustness of the extracted results is examined by checking various resummation techniques and by running the simulations with various choices for the propagators and vertex functions.It turns out that dressing the lines in the diagrams as much as possible is not always the optimal choice.We also identify classes of dominant diagrams for the one-body and two-body self-energy in the region of strong interaction.These dominant diagrams turn out to be the leading processes of the strong-coupling limit.The quasiparticle energies and Z-factor are obtained as a function of the interaction strength.We find that the DiagMC results for the molecule and polaron properties are very similar to those obtained with a variational ansatz.Surprisingly, this variational ansatz gives very good predictions for the quasiparticle residue even when this residue is significantly smaller than one.

Related papers

Quasiparticle properties of an impurity in a Fermi gas

Jonas Vlietinck, J. Ryckebusch, Kris Van Houcke

We report a study of a spin-down impurity strongly coupled to a spin-up Fermi sea (a so-called Fermi polaron) with the diagrammatic Monte Carlo (DiagMC) technique. Conditions of zero temperature and three dimensions are considered for an ultracold atomic gas with resonant interactions in the zero-range limit. A Feynman diagrammatic series is developed for the one-body and two-body propagators providing information about the polaron and molecule channel, respectively. The DiagMC technique allows us to reach diagram orders that are high enough for extrapolation to infinite order. The robustness of the extracted results is examined by checking various resummation techniques and by running the simulations with various choices for the propagators and vertex functions. It turns out that dressing the lines in the diagrams as much as possible is not always the optimal choice. We also identify classes of dominant diagrams for the one-body and two-body self-energy in the region of strong interaction. These dominant diagrams turn out to be the leading processes of the strong-coupling limit. The quasiparticle energies and $Z$ factor are obtained as a function of the interaction strength. We find that the DiagMC results for the molecule and polaron properties are very similar to those obtained with a variational ansatz. Surprisingly, this variational ansatz gives very good predictions for the quasiparticle residue even when this residue is significantly less than 1.

arXiv:1302.3505•Physical Review B•Mar 25, 2013•73 citations

Diagrammatic Monte Carlo study of the Fermi polaron in two dimensions

Jonas Vlietinck, J. Ryckebusch, Kris Van Houcke

We study the properties of the two-dimensional Fermi polaron model in which an impurity attractively interacts with a Fermi sea of particles in the zero-range limit. We use a diagrammatic Monte Carlo (DiagMC) method which allows us to sample a Feynman diagrammatic series to very high order. The convergence properties of the series and the role of multiple particle-hole excitations are discussed. We study the polaron and molecule energy as a function of the coupling strength, revealing a transition from a polaron to a molecule in the ground state. We find a value for the critical interaction strength which complies with the experimentally measured one and predictions from variational methods. For all considered interaction strengths, the polaron $Z$ factor from the full diagrammatic series almost coincides with the one-particle-hole result. We also formally link the DiagMC and the variational approaches for the polaron problem at hand.

arXiv:1402.4015•Physical Review B•Feb 19, 2014•50 citations

Diagrammatic Monte Carlo study of a mass-imbalanced Fermi-polaron system

Peter Kroiß, Lode Pollet

We apply the diagrammatic Monte Carlo approach to three-dimensional Fermi-polaron systems with mass imbalance, where an impurity interacts resonantly with a noninteracting Fermi sea whose atoms have a different mass. This method allows us to go beyond frequently used variational techniques by stochastically summing all relevant impurity Feynman diagrams up to a maximum expansion order limited by the sign problem. Polaron energy and quasiparticle residue can be accurately determined over a broad range of impurity masses. Furthermore, the spectral function of an imbalanced polaron demonstrates the stability of the quasiparticle and allows us to locate in addition also the repulsive polaron as an excited state. The quantitative exactness of two-particle-hole wave functions is investigated, resulting in a relative lowering of polaronic energies in the mass-imbalance phase diagram. Tan's contact coefficient for the mass-balanced polaron system is found in good agreement with variational methods. Mass-imbalanced systems can be studied experimentally by ultracold atom mixtures such as $^{6}\mathrm{Li}\ensuremath{-}^{40}\mathrm{K}$.

arXiv:1411.1562•Physical Review B•Apr 20, 2015•32 citations

Bold diagrammatic Monte Carlo: A generic technique for polaron (and many-body?) problems

Nikolay V. Prokof'ev, Boris V. Svistunov

We develop a Monte Carlo scheme for sampling series of Feynman diagrams for the proper self-energy which are self-consistently expressed in terms of renormalized particle propagators. This approach is used to solve the problem of a single spin-down fermion resonantly interacting with the Fermi gas of spin-up particles. Though the original series based on bare propagators are sign-alternating and divergent one can still determine the answer behind them by using two strategies (separately or together): (i) using proper series re-summation techniques, and (ii) introducing renormalized propagators which are defined in terms of the simulated proper self-energy, i.e. making the entire scheme self-consistent. Our solution is important for understanding the phase diagram and properties of the BCS-BEC crossover in the strongly imbalanced regime. On the technical side, we develop a generic sign-problem tolerant method for exact numerical solution of polaron-type models, and, possibly, of the interacting many-body Hamiltonians.

arXiv:0801.0911•Jan 7, 2008•2 citations

Polarons and dimerons in the two-dimensional attractive Hubbard model

Gerard Pascual, Jordi Boronat, Kris Van Houcke

A two-dimensional (2D) spin-up ideal Fermi gas interacting attractively with a spin-down impurity in the continuum undergoes, at zero temperature, a first-order phase transition from a polaron to a dimeron state. Here, we study a similar system on a square lattice, by considering the attractive 2D Fermi-Hubbard model with a single spin-down and a finite filling fraction of spin-up fermions. We study polaron and dimeron quasiparticle properties via variational ansatz up to one particle-hole excitation. Moreover, we develop a determinant diagrammatic Monte Carlo algorithm for this problem based on expansion in bare on-site coupling U . This algorithm turns out to be sign-problem-free at any filling of spin-up fermions, allowing one to sample very high diagram order (larger than 200 in our study) and to do simulations for large U / t (we go up to U / t = − 20 with t the hopping strength). Both methods give qualitatively consistent results. With variational ansatz we go to even larger on-site attraction. At very low spin-up filling fraction, we observe the polaron-to-dimeron transition, in agreement with the continuum case. Upon increasing the filling fraction, however, the transition shifts to higher values of | U | / t and the transition disappears beyond a filling fraction of about 20 % . In this region, the polaron state always gives a lower energy and has a finite quasiparticle residue. Our findings are directly relevant to cold atom experiments with 2D optical lattices: A small and finite density of spin-down impurities in the ground state will form a superfluid at strong coupling at low spin-up filling fraction. Above some critical spin-up filling fraction, the system is expected to remain a normal Fermi liquid, even in the strong-coupling limit.

arXiv:2411.19725•Physical Review Research•Sep 11, 2025•2 citations

Polaron-to-molecule transition in a strongly imbalanced Fermi gas

Matthias Punk, Philipp T. Dumitrescu, W. Zwerger

A single down-spin fermion with an attractive zero-range interaction with a Fermi sea of up-spin fermions forms a polaronic quasiparticle. The associated quasiparticle weight vanishes beyond a critical strength of the attractive interaction, where a many-body bound state is formed. From a variational wave function in the molecular limit, we determine the critical value for the polaron-to-molecule transition. The value agrees well with the diagrammatic Monte Carlo results of Prokof'ev and Svistunov and is consistent with recent rf-spectroscopy measurements of the quasiparticle weight by Schirotzek et al. [Phys. Rev. Lett. 102, 230402 (2009)]. In addition, we calculate the contact coefficient of the strongly imbalanced gas, using the adiabatic theorem of Tan and discuss the implications of the polaron-to-molecule transition for the phase diagram of the attractive Fermi gas at finite imbalance.

arXiv:0908.1343•Physical Review A•Nov 6, 2009•190 citations
We report on a study of a spin-down impurity strongly coupled to a spin-up Fermi sea (a so-called Fermi polaron) with the diagrammatic Monte-Carlo (DiagMC) technique.Conditions of zero temperature and three dimensions are considered for an ultracold atomic gas with resonant interactions in the zero-range limit.A Feynman diagrammatic series is developed for the one-body and two-body propagators providing information about the polaron and molecule channel respectively.The DiagMC technique allows us to reach diagram orders that are high enough for extrapolation to infinite order.The robustness of the extracted results is examined by checking various resummation techniques and by running the simulations with various choices for the propagators and vertex functions.It turns out that dressing the lines in the diagrams as much as possible is not always the optimal choice.We also identify classes of dominant diagrams for the one-body and two-body self-energy in the region of strong interaction.These dominant diagrams turn out to be the leading processes of the strong-coupling limit.The quasiparticle energies and Z-factor are obtained as a function of the interaction strength.We find that the DiagMC results for the molecule and polaron properties are very similar to those obtained with a variational ansatz.Surprisingly, this variational ansatz gives very good predictions for the quasiparticle residue even when this residue is significantly smaller than one.
We report a study of a spin-down impurity strongly coupled to a spin-up Fermi sea (a so-called Fermi polaron) with the diagrammatic Monte Carlo (DiagMC) technique. Conditions of zero temperature and three dimensions are considered for an ultracold atomic gas with resonant interactions in the zero-range limit. A Feynman diagrammatic series is developed for the one-body and two-body propagators providing information about the polaron and molecule channel, respectively. The DiagMC technique allows us to reach diagram orders that are high enough for extrapolation to infinite order. The robustness of the extracted results is examined by checking various resummation techniques and by running the simulations with various choices for the propagators and vertex functions. It turns out that dressing the lines in the diagrams as much as possible is not always the optimal choice. We also identify classes of dominant diagrams for the one-body and two-body self-energy in the region of strong interaction. These dominant diagrams turn out to be the leading processes of the strong-coupling limit. The quasiparticle energies and ZZZ factor are obtained as a function of the interaction strength. We find that the DiagMC results for the molecule and polaron properties are very similar to those obtained with a variational ansatz. Surprisingly, this variational ansatz gives very good predictions for the quasiparticle residue even when this residue is significantly less than 1.
We study the properties of the two-dimensional Fermi polaron model in which an impurity attractively interacts with a Fermi sea of particles in the zero-range limit. We use a diagrammatic Monte Carlo (DiagMC) method which allows us to sample a Feynman diagrammatic series to very high order. The convergence properties of the series and the role of multiple particle-hole excitations are discussed. We study the polaron and molecule energy as a function of the coupling strength, revealing a transition from a polaron to a molecule in the ground state. We find a value for the critical interaction strength which complies with the experimentally measured one and predictions from variational methods. For all considered interaction strengths, the polaron ZZZ factor from the full diagrammatic series almost coincides with the one-particle-hole result. We also formally link the DiagMC and the variational approaches for the polaron problem at hand.
We apply the diagrammatic Monte Carlo approach to three-dimensional Fermi-polaron systems with mass imbalance, where an impurity interacts resonantly with a noninteracting Fermi sea whose atoms have a different mass. This method allows us to go beyond frequently used variational techniques by stochastically summing all relevant impurity Feynman diagrams up to a maximum expansion order limited by the sign problem. Polaron energy and quasiparticle residue can be accurately determined over a broad range of impurity masses. Furthermore, the spectral function of an imbalanced polaron demonstrates the stability of the quasiparticle and allows us to locate in addition also the repulsive polaron as an excited state. The quantitative exactness of two-particle-hole wave functions is investigated, resulting in a relative lowering of polaronic energies in the mass-imbalance phase diagram. Tan's contact coefficient for the mass-balanced polaron system is found in good agreement with variational methods. Mass-imbalanced systems can be studied experimentally by ultracold atom mixtures such as $^{6}\mathrm{Li}\ensuremath{-}^{40}\mathrm{K}$.
We develop a Monte Carlo scheme for sampling series of Feynman diagrams for the proper self-energy which are self-consistently expressed in terms of renormalized particle propagators. This approach is used to solve the problem of a single spin-down fermion resonantly interacting with the Fermi gas of spin-up particles. Though the original series based on bare propagators are sign-alternating and divergent one can still determine the answer behind them by using two strategies (separately or together): (i) using proper series re-summation techniques, and (ii) introducing renormalized propagators which are defined in terms of the simulated proper self-energy, i.e. making the entire scheme self-consistent. Our solution is important for understanding the phase diagram and properties of the BCS-BEC crossover in the strongly imbalanced regime. On the technical side, we develop a generic sign-problem tolerant method for exact numerical solution of polaron-type models, and, possibly, of the interacting many-body Hamiltonians.
A two-dimensional (2D) spin-up ideal Fermi gas interacting attractively with a spin-down impurity in the continuum undergoes, at zero temperature, a first-order phase transition from a polaron to a dimeron state. Here, we study a similar system on a square lattice, by considering the attractive 2D Fermi-Hubbard model with a single spin-down and a finite filling fraction of spin-up fermions. We study polaron and dimeron quasiparticle properties via variational ansatz up to one particle-hole excitation. Moreover, we develop a determinant diagrammatic Monte Carlo algorithm for this problem based on expansion in bare on-site coupling U . This algorithm turns out to be sign-problem-free at any filling of spin-up fermions, allowing one to sample very high diagram order (larger than 200 in our study) and to do simulations for large U / t (we go up to U / t = − 20 with t the hopping strength). Both methods give qualitatively consistent results. With variational ansatz we go to even larger on-site attraction. At very low spin-up filling fraction, we observe the polaron-to-dimeron transition, in agreement with the continuum case. Upon increasing the filling fraction, however, the transition shifts to higher values of | U | / t and the transition disappears beyond a filling fraction of about 20 % . In this region, the polaron state always gives a lower energy and has a finite quasiparticle residue. Our findings are directly relevant to cold atom experiments with 2D optical lattices: A small and finite density of spin-down impurities in the ground state will form a superfluid at strong coupling at low spin-up filling fraction. Above some critical spin-up filling fraction, the system is expected to remain a normal Fermi liquid, even in the strong-coupling limit.
A single down-spin fermion with an attractive zero-range interaction with a Fermi sea of up-spin fermions forms a polaronic quasiparticle. The associated quasiparticle weight vanishes beyond a critical strength of the attractive interaction, where a many-body bound state is formed. From a variational wave function in the molecular limit, we determine the critical value for the polaron-to-molecule transition. The value agrees well with the diagrammatic Monte Carlo results of Prokof'ev and Svistunov and is consistent with recent rf-spectroscopy measurements of the quasiparticle weight by Schirotzek et al. [Phys. Rev. Lett. 102, 230402 (2009)]. In addition, we calculate the contact coefficient of the strongly imbalanced gas, using the adiabatic theorem of Tan and discuss the implications of the polaron-to-molecule transition for the phase diagram of the attractive Fermi gas at finite imbalance.