Efficient variational approach to the Fermi polaron problem in two dimensions, both in and out of equilibrium
Pith reviewed 2026-05-24 11:02 UTC · model grok-4.3
The pith
Non-Gaussian variational states track the two-dimensional Fermi polaron across its transition to molecules without truncating particle-hole excitations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The non-Gaussian variational states enable an unbiased analysis of the polaron-to-molecule phase transition without relying on truncations in the total number of particle-hole excitations, and they allow exploration of long-time polaron evolution and spectral properties.
What carries the argument
Non-Gaussian variational states whose functional form is chosen to capture impurity-fermion correlations beyond Gaussian level.
If this is right
- Ground-state energy and quasiparticle residue remain consistent with Monte Carlo data across the transition.
- Long-time dynamics after a quench become accessible without artificial cutoffs on excitations.
- Radiofrequency spectroscopy can be designed to extract both polaronic and molecular spectral functions.
- Fermionic density profiles near the impurity exhibit either persistent oscillations or rapid relaxation depending on interaction strength.
Where Pith is reading between the lines
- The same variational family might be tested on related impurity problems such as the Bose polaron in two dimensions.
- Time-dependent spectra obtained this way could be compared directly to existing ultracold-atom radio-frequency data.
- If the ansatz remains accurate at strong coupling, it offers a route to parameter-free predictions for quasiparticle lifetimes.
Load-bearing premise
The chosen non-Gaussian variational family is flexible enough to capture the essential physics of the 2D Fermi polaron across the transition without missing important correlations.
What would settle it
Exact diagonalization or quantum Monte Carlo results that differ quantitatively from the variational energies or residues near the polaron-molecule transition point would show the ansatz misses key correlations.
Figures
read the original abstract
We develop a non-Gaussian variational approach that enables us to study both equilibrium and far-from-equilibrium physics of the two-dimensional Fermi polaron. This method provides an unbiased analysis of the polaron-to-molecule phase transition without relying on truncations in the total number of particle-hole excitations. Our results -- which include the ground state energy and quasiparticle residue -- are in qualitative agreement with the known Monte Carlo calculations. The main advantage of the non-Gaussian states compared to conventional numerical methods is that they enable us to explore long-time polaron evolution and, in particular, study various spectral properties accessible to both solid-state and ultracold atom experiments. We design two types of radiofrequency spectroscopies to measure polaronic and molecular spectral functions. Depending on the parameter regime, we find that these spectral functions and fermionic density profiles near the impurity display either long-lived oscillations between the repulsive and attractive polaron branches or exhibit fast relaxational dynamics to the molecular state.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a non-Gaussian variational method for the two-dimensional Fermi polaron that avoids explicit truncations on the number of particle-hole excitations. It reports qualitative agreement with existing Monte Carlo results for the ground-state energy and quasiparticle residue, and uses the ansatz to compute long-time dynamics and radiofrequency spectra that distinguish polaronic versus molecular branches.
Significance. If the chosen non-Gaussian manifold proves sufficiently expressive, the approach would offer a useful route to far-from-equilibrium polaron physics that is difficult to access with truncated expansions or sign-problematic Monte Carlo. The absence of an explicit particle-hole cutoff is a genuine technical advantage over conventional variational treatments.
major comments (2)
- [Abstract] Abstract: only qualitative agreement with Monte Carlo is stated for ground-state energy and residue; no numerical values, error estimates, or direct comparison tables are supplied, leaving the accuracy of the variational family near the polaron-molecule transition unquantified.
- [Abstract] Abstract: the non-equilibrium claims (long-time evolution, spectral functions, and density profiles) rest on the untested assumption that the specific non-Gaussian functional form captures all relevant correlations; no convergence tests with respect to variational parameters or comparisons against exact benchmarks in the dynamical regime are provided.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. We address the two major comments point by point below. Revisions have been made to the abstract and additional material has been added to the manuscript.
read point-by-point responses
-
Referee: [Abstract] Abstract: only qualitative agreement with Monte Carlo is stated for ground-state energy and residue; no numerical values, error estimates, or direct comparison tables are supplied, leaving the accuracy of the variational family near the polaron-molecule transition unquantified.
Authors: We agree that the abstract would benefit from greater quantification. In the revised manuscript we have updated the abstract to report explicit numerical values for the ground-state energy and quasiparticle residue at representative interaction strengths, including near the polaron-molecule transition. We now explicitly reference the figures and tables that contain the direct Monte Carlo comparisons, and we note that the variational energies constitute rigorous upper bounds, which serves as a built-in error indicator. revision: yes
-
Referee: [Abstract] Abstract: the non-equilibrium claims (long-time evolution, spectral functions, and density profiles) rest on the untested assumption that the specific non-Gaussian functional form captures all relevant correlations; no convergence tests with respect to variational parameters or comparisons against exact benchmarks in the dynamical regime are provided.
Authors: We acknowledge that the dynamical results rely on the expressiveness of the chosen non-Gaussian manifold. To strengthen the presentation we have added an appendix containing convergence tests obtained by systematically increasing the number of variational parameters. While exact dynamical benchmarks for the 2D Fermi polaron do not exist in the literature, we have expanded the discussion to include recovery of known equilibrium limits, perturbative regimes, and internal consistency checks. A new paragraph has been inserted that explicitly states the assumptions of the ansatz and its expected range of validity. revision: partial
Circularity Check
No significant circularity in variational derivation
full rationale
The paper introduces a new non-Gaussian variational family as an ansatz for the 2D Fermi polaron, claiming it avoids explicit particle-hole truncations. Ground-state energies and residues are compared qualitatively to independent Monte Carlo results from the literature, with no indication that parameters are fitted to those data or that predictions reduce to the inputs by construction. No self-citation chains, uniqueness theorems, or renamings of known results are invoked as load-bearing steps in the provided text. The derivation chain is therefore self-contained: the variational states are chosen, the equations of motion or energy minimization are solved within that manifold, and external benchmarks are used only for validation rather than as hidden inputs.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
evidence for the polaron-to-molecule transition. No- tably, one finds that in the case of equal masses of impu- rity and bath particles, one can obtain a good description of many-body polaronic states by including only a sin- gle particle-hole excitation [33]. These so-called Chevy ansatz (CA) wave functions work surprisingly well even at unitarity when th...
-
[2]
P. W. Anderson, Phys. Rev. Lett. 18, 1049 (1967)
work page 1967
- [3]
-
[4]
G. D. Mahan, Many-Particle Physics (Springer US, Boston, MA, 2000)
work page 2000
-
[5]
P. W. Anderson and G. Yuval, Phys. Rev. Lett. 23, 89 (1969)
work page 1969
- [6]
- [7]
-
[8]
M. Hentschel, D. Ullmo, and H. U. Baranger, Phys. Rev. B 72, 035310 (2005)
work page 2005
-
[9]
Y. V. Nazarov and Y. M. Blanter, Quantum Trans- port: Introduction to Nanoscience (Cambridge Univer- 12 sity Press, Cambridge, 2009)
work page 2009
-
[10]
A. Schirotzek, C.-H. Wu, A. Sommer, and M. W. Zwier- lein, Phys. Rev. Lett. 102, 230402 (2009)
work page 2009
- [11]
-
[12]
R. Schmidt, H. R. Sadeghpour, and E. Demler, Phys. Rev. Lett. 116, 105302 (2016)
work page 2016
- [13]
-
[14]
M. P. Sarachik, E. Corenzwit, and L. D. Longinotti, Phys. Rev. 135, A1041 (1964)
work page 1964
-
[15]
A. C. Hewson, The Kondo Problem to Heavy Fermions , Cambridge Studies in Magnetism (Cambridge University Press, Cambridge, 1993)
work page 1993
- [16]
-
[17]
S. Kirchner, S. Paschen, Q. Chen, S. Wirth, D. Feng, J. D. Thompson, and Q. Si, Rev. Mod. Phys. 92, 011002 (2020)
work page 2020
-
[18]
V. G. Storchak and N. V. Prokof’ev, Rev. Mod. Phys. 70, 929 (1998)
work page 1998
-
[19]
Rosch, Advances in Physics 48, 295 (1999)
A. Rosch, Advances in Physics 48, 295 (1999)
work page 1999
- [20]
-
[21]
P. E. Dolgirev, Y.-F. Qu, M. B. Zvonarev, T. Shi, and E. Demler, Phys. Rev. X 11, 041015 (2021)
work page 2021
-
[22]
G. E. Astrakharchik and L. P. Pitaevskii, Phys. Rev. A 70, 013608 (2004)
work page 2004
-
[23]
A. Y. Cherny, J.-S. Caux, and J. Brand, Front. Phys. 7, 54 (2012)
work page 2012
- [24]
-
[25]
C. Kohstall, M. Zaccanti, M. Jag, A. Trenkwalder, P. Massignan, G. M. Bruun, F. Schreck, and R. Grimm, Nature 485, 615 (2012)
work page 2012
-
[26]
M. Koschorreck, D. Pertot, E. Vogt, B. Fr¨ ohlich, M. Feld, and M. K¨ ohl, Nature485, 619 (2012)
work page 2012
- [27]
- [28]
-
[29]
Z. Yan, P. B. Patel, B. Mukherjee, R. J. Fletcher, J. Struck, and M. W. Zwierlein, Phys. Rev. Lett. 122, 093401 (2019)
work page 2019
-
[30]
G. Ness, C. Shkedrov, Y. Florshaim, O. K. Diessel, J. von Milczewski, R. Schmidt, and Y. Sagi, Phys. Rev. X 10, 041019 (2020)
work page 2020
-
[31]
C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. Mod. Phys. 82, 1225 (2010)
work page 2010
- [32]
-
[33]
N. V. Prokof’ev and B. V. Svistunov, Phys. Rev. B 77, 125101 (2008)
work page 2008
- [34]
-
[35]
R. Combescot, A. Recati, C. Lobo, and F. Chevy, Phys. Rev. Lett. 98, 180402 (2007)
work page 2007
- [36]
-
[37]
S. Z¨ ollner, G. M. Bruun, and C. J. Pethick, Phys. Rev. A 83, 021603 (2011)
work page 2011
-
[38]
M. M. Parish, Phys. Rev. A 83, 051603 (2011)
work page 2011
-
[39]
M. M. Parish and J. Levinsen, Phys. Rev. A 87, 033616 (2013)
work page 2013
- [40]
-
[41]
C. Peng, R. Liu, W. Zhang, and X. Cui, Phys. Rev. A 103, 063312 (2021)
work page 2021
- [42]
-
[43]
D. K. Efimkin and A. H. MacDonald, Phys. Rev. B 95, 035417 (2017)
work page 2017
-
[44]
C. Fey, P. Schmelcher, A. Imamoglu, and R. Schmidt, Phys. Rev. B 101, 195417 (2020)
work page 2020
-
[45]
C. Kuhlenkamp, M. Knap, M. Wagner, R. Schmidt, and A. Imamo˘ g˘ glu, Phys. Rev. Lett.129, 037401 (2022)
work page 2022
-
[46]
T. Shi, E. Demler, and J. Ignacio Cirac, Annals of Physics 390, 245 (2018)
work page 2018
-
[47]
D. A. Abanin and L. S. Levitov, Phys. Rev. Lett. 94, 186803 (2005)
work page 2005
-
[48]
M. Knap, A. Shashi, Y. Nishida, A. Imambekov, D. A. Abanin, and E. Demler, Phys. Rev. X 2, 041020 (2012)
work page 2012
-
[49]
R. Schmidt, M. Knap, D. A. Ivanov, J.-S. You, M. Cetina, and E. Demler, Rep. Prog. Phys. 81, 024401 (2018)
work page 2018
-
[50]
J. B. McGuire, J. Math. Phys. 6, 432 (1965)
work page 1965
-
[51]
J. B. McGuire, J. Math. Phys. 7, 123 (1966)
work page 1966
-
[52]
C. J. M. Mathy, M. B. Zvonarev, and E. Demler, Nature Phys 8, 881 (2012)
work page 2012
-
[53]
M. Knap, C. J. M. Mathy, M. Ganahl, M. B. Zvonarev, and E. Demler, Phys. Rev. Lett. 112, 015302 (2014)
work page 2014
-
[54]
O. Gamayun, A. G. Pronko, and M. B. Zvonarev, New J. Phys. 18, 045005 (2016)
work page 2016
-
[55]
O. Gamayun, O. Lychkovskiy, E. Burovski, M. Malcom- son, V. V. Cheianov, and M. B. Zvonarev, Phys. Rev. Lett. 120, 220605 (2018)
work page 2018
-
[56]
O. Gamayun, O. Lychkovskiy, and M. B. Zvonarev, Sci- Post Phys. 8, 53 (2020)
work page 2020
-
[57]
R. Schmidt, T. Enss, V. Pietil¨ a, and E. Demler, Phys. Rev. A 85, 021602 (2012)
work page 2012
-
[58]
M. M. Parish and J. Levinsen, Phys. Rev. B 94, 184303 (2016)
work page 2016
-
[59]
W. E. Liu, J. Levinsen, and M. M. Parish, Phys. Rev. Lett. 122, 205301 (2019)
work page 2019
-
[60]
W. E. Liu, Z.-Y. Shi, M. M. Parish, and J. Levinsen, Phys. Rev. A 102, 023304 (2020)
work page 2020
-
[61]
H. S. Adlong, W. E. Liu, F. Scazza, M. Zaccanti, N. D. Oppong, S. F¨ olling, M. M. Parish, and J. Levinsen, Phys. Rev. Lett. 125, 133401 (2020)
work page 2020
- [62]
-
[63]
T. D. Lee, F. E. Low, and D. Pines, Phys. Rev. 90, 297 (1953)
work page 1953
-
[64]
E. Demler, S.-C. Zhang, N. Bulut, and D. J. Scalapino, Int. J. Mod. Phys. B 10, 2137 (1996)
work page 1996
- [65]
- [66]
-
[67]
J. Vlietinck, J. Ryckebusch, and K. Van Houcke, Phys. Rev. B 89, 085119 (2014)
work page 2014
- [68]
-
[69]
Y. Tang, J. Gu, S. Liu, K. Watanabe, T. Taniguchi, J. Hone, K. F. Mak, and J. Shan, Nat. Nanotechnol. 16, 52 (2021)
work page 2021
-
[70]
J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 (1966)
work page 1966
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.