Ab initio study of the neutron and Fermi polarons on the lattice
Pith reviewed 2026-05-18 08:55 UTC · model grok-4.3
The pith
Lattice quantum Monte Carlo calculations give benchmark energies for Fermi and neutron polarons.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors establish that auxiliary-field quantum Monte Carlo on the lattice, with parameters fixed to two-body physics through the parametric matrix model and Luscher's formula, produces reliable energies for the polaron in both cold-atom and nuclear contexts that can serve as stringent benchmarks for future theoretical and experimental work.
What carries the argument
Auxiliary-field quantum Monte Carlo applied to a lattice Hamiltonian for the fermionic polaron, with parameters tuned by the parametric matrix model to reproduce two-body energies in a periodic box according to Luscher's formula.
If this is right
- The computed polaron energies provide direct reference points for cold-atom experiments across scattering lengths.
- The neutron polaron results constrain models of neutron-rich matter in nuclear physics.
- The same lattice framework can be applied to other impurity problems in fermionic systems.
- Future many-body calculations in both atomic and nuclear regimes can be tested against these lattice results.
Where Pith is reading between the lines
- These benchmarks could help clarify the equation of state in low-density neutron matter by offering an independent lattice anchor.
- The method's generality suggests it could be extended to study polarons in mixed atomic-nuclear hybrid systems or other tunable interactions.
- Varying the lattice volume or spacing in follow-up runs would test how well the current tuning suppresses finite-size effects.
Load-bearing premise
Tuning the lattice Hamiltonian parameters to two-body energies via Luscher's formula and the parametric matrix model introduces no uncontrolled discretization or finite-volume errors that propagate into the many-body polaron energies.
What would settle it
A precise experimental measurement or independent exact calculation of the polaron energy at unitarity or for the neutron case that differs significantly from the lattice Monte Carlo value without an identified source of error.
Figures
read the original abstract
We have used the auxiliary-field quantum Monte Carlo (AFQMC) many-body approach on the lattice to study the equation of state for a fermionic impurity interacting with a background sea of spin-polarized fermions. The impurity, or polaron, is an interesting system in both cold atomic and nuclear physics. Our approach is general, and we are able to straightforwardly study the polaron across these regimes. We first study the Fermi polaron at unitarity and for a wide range of scattering lengths, comparing against previous theoretical and experimental studies. We then explore the neutron polaron which has been shown to be an important constraint for nuclear physics. We have also employed the recently developed parametric matrix model to emulate AFQMC solutions to the two-body problem on the lattice, to accelerate the tuning of our lattice Hamiltonian parameters directly to two-body energies in a periodic box, following Luscher's formula. Our lattice quantum Monte Carlo results for the polaron in both a cold atomic and nuclear physics context can serve as stringent benchmarks for future theoretical and experimental research.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports an auxiliary-field quantum Monte Carlo (AFQMC) lattice study of the equation of state for a fermionic impurity (polaron) interacting with a spin-polarized Fermi sea. It examines the Fermi polaron at unitarity and across a range of scattering lengths, with comparisons to prior theory and experiment, and extends the method to the neutron polaron. Lattice Hamiltonian parameters are tuned via a parametric matrix model to match two-body energies in a periodic box using Luscher's formula; the resulting many-body energies are presented as benchmarks for cold-atom and nuclear-physics applications.
Significance. If discretization and finite-volume errors are shown to be controlled at the level of the reported precision, the numerical results would constitute useful ab initio benchmarks bridging ultracold atoms and nuclear structure. The generality of the AFQMC lattice framework across regimes and the practical acceleration provided by the parametric matrix model for two-body tuning are clear strengths. The work supplies concrete numerical data rather than a new analytic derivation, which is appropriate for a benchmark-oriented study.
major comments (2)
- [Section describing the parametric matrix model and two-body tuning] The central benchmark claim rests on the assumption that tuning the lattice parameters to two-body energies via the parametric matrix model and Luscher's formula introduces no uncontrolled errors that propagate into the many-body polaron energies. No explicit sensitivity analysis or residual-mismatch propagation test is provided for the extracted polaron energies.
- [Results section on Fermi and neutron polaron energies] Comparisons to previous theoretical and experimental studies for the Fermi polaron are stated, yet the manuscript does not supply full error budgets, complete data tables, or volume-extrapolation checks. This limits the ability to judge whether the new results achieve the precision needed to serve as stringent benchmarks.
minor comments (2)
- [Hamiltonian and parameter definitions] Clarify the notation distinguishing the lattice interaction parameters from the physical scattering length and effective range throughout the text and equations.
- [Figures and captions] Ensure all figures displaying polaron energies include explicit error bars, lattice-spacing labels, and volume sizes consistent with the accompanying text.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments, which help clarify the presentation of our benchmark results. We address each major comment below and have revised the manuscript to incorporate the requested analyses and supporting data.
read point-by-point responses
-
Referee: [Section describing the parametric matrix model and two-body tuning] The central benchmark claim rests on the assumption that tuning the lattice parameters to two-body energies via the parametric matrix model and Luscher's formula introduces no uncontrolled errors that propagate into the many-body polaron energies. No explicit sensitivity analysis or residual-mismatch propagation test is provided for the extracted polaron energies.
Authors: We agree that an explicit sensitivity analysis strengthens the central benchmark claim. In the revised manuscript we have added a dedicated subsection that quantifies the propagation of two-body tuning uncertainties. We vary the lattice interaction parameters within the residual mismatch allowed by the parametric matrix model and Luscher's formula (at the level of the two-body energy precision), recompute the many-body polaron energies on representative volumes, and demonstrate that the induced shifts remain well below the statistical Monte Carlo uncertainties reported for the final results. This test is now included for both the unitary Fermi polaron and the neutron polaron cases. revision: yes
-
Referee: [Results section on Fermi and neutron polaron energies] Comparisons to previous theoretical and experimental studies for the Fermi polaron are stated, yet the manuscript does not supply full error budgets, complete data tables, or volume-extrapolation checks. This limits the ability to judge whether the new results achieve the precision needed to serve as stringent benchmarks.
Authors: We acknowledge that the original manuscript presented only selected comparisons and summary figures. The revised version now includes (i) complete tabulated data for the polaron energy at all studied scattering lengths and lattice volumes, (ii) a detailed error budget separating statistical, systematic (including discretization and tuning), and extrapolation uncertainties, and (iii) explicit volume-extrapolation plots and fits for the Fermi polaron at unitarity and selected scattering lengths. The extrapolated infinite-volume values with combined uncertainties are reported and used in the comparisons to prior theory and experiment. The same level of documentation is provided for the neutron polaron results. revision: yes
Circularity Check
No significant circularity; self-contained numerical benchmarks
full rationale
The paper calibrates lattice Hamiltonian parameters to two-body energies via Luscher's formula using the parametric matrix model, then computes many-body polaron energies directly with AFQMC. This is explicit external calibration to a standard result (Luscher), not a reduction of the central many-body claim to the fit by construction. No self-citation load-bearing for the benchmark utility, no self-definitional steps, and no renaming or ansatz smuggling. The results are presented as numerical data rather than derived predictions equivalent to inputs.
Axiom & Free-Parameter Ledger
free parameters (1)
- lattice interaction parameters
axioms (1)
- domain assumption Auxiliary-field quantum Monte Carlo accurately samples the fermionic many-body problem without prohibitive sign issues
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Our lattice Hamiltonian is written in second quantization... tuned... via Lüscher’s formula and the parametric matrix model
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
AFQMC... constrained path approximation... to deal with the fermion sign problem
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
L. D. Landau, Phys. Z. Sowjet.3, 664 (1933)
work page 1933
-
[2]
P. Massignan, M. Zaccanti, and G. M. Bruun, Rep. Prog. Phys.77, 034401 (2014)
work page 2014
-
[3]
C. Franchini, M. Reticcioli, M. Setvin, and U. Diebold, Nat. Rev. Mater.6, 560 (2021)
work page 2021
- [4]
-
[5]
M. M. Forbes, A. Gezerlis, K. Hebeler, T. Lesinski, and A. Schwenk, Phys. Rev. C89, 041301 (2014)
work page 2014
- [6]
- [7]
-
[8]
S.-G. Zhou, J. Meng, P. Ring, and E.-G. Zhao, Phys. Rev. C82, 011301 (2010)
work page 2010
-
[9]
J. Tanaka, Z. Yang, S. Typel, S. Adachi, S. Bai, P. van Beek, D. Beaumel, Y. Fujikawa, J. Han, S. Heil, S. Huang, A. Inoue, Y. Jiang, M. Kn¨ osel, N. Kobayashi, Y. Kubota, W. Liu, J. Lou, Y. Maeda, Y. Matsuda, K. Miki, S. Nakamura, K. Ogata, V. Panin, H. Scheit, F. Schindler, P. Schrock, D. Symochko, A. Tamii, T. Ue- saka, V. Wagner, K. Yoshida, J. Zenihi...
work page 2021
-
[10]
C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. Mod. Phys.82, 1225 (2010)
work page 2010
-
[11]
A. Schirotzek, C.-H. Wu, A. Sommer, and M. W. Zwier- lein, Phys. Rev. Lett.102, 230402 (2009)
work page 2009
-
[12]
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
-
[13]
I. Fritsche, C. Baroni, E. Dobler, E. Kirilov, B. Huang, R. Grimm, G. M. Bruun, and P. Massignan, Phys. Rev. A103, 053314 (2021)
work page 2021
- [14]
-
[15]
C. Peng, R. Liu, W. Zhang, and X. Cui, Phys. Rev. A 103, 063312 (2021)
work page 2021
- [16]
- [17]
-
[18]
M. Punk, P. T. Dumitrescu, and W. Zwerger, Phys. Rev. A80, 053605 (2009)
work page 2009
-
[19]
K. Van Houcke, F. Werner, and R. Rossi, Phys. Rev. B 101, 045134 (2020)
work page 2020
-
[20]
C. Lobo, A. Recati, S. Giorgini, and S. Stringari, Phys. Rev. Lett.97, 200403 (2006)
work page 2006
-
[21]
S. Bour, D. Lee, H.-W. Hammer, and U.-G. Meißner, Phys. Rev. Lett.115, 185301 (2015)
work page 2015
- [22]
-
[23]
M. H. Kalos, Phys. Rev.128, 1791 (1962)
work page 1962
-
[24]
J. B. Anderson, J. Chem. Phys.63, 1499 (1975)
work page 1975
- [25]
-
[26]
K. E. Schmidt and S. Fantoni, Phys. Lett. B446, 99 (1999)
work page 1999
- [27]
- [28]
- [29]
-
[30]
A. Roggero, A. Mukherjee, and F. Pederiva, Phys. Rev. Lett.112, 221103 (2014)
work page 2014
- [31]
-
[32]
S. Gandolfi, A. Lovato, J. Carlson, and K. E. Schmidt, Phys. Rev. C90, 061306 (2014)
work page 2014
-
[33]
D. Lonardoni, I. Tews, S. Gandolfi, and J. Carlson, Phys. Rev. Res.2, 022033 (2020)
work page 2020
- [34]
-
[35]
A. Gezerlis, S. Gandolfi, K. E. Schmidt, and J. Carlson, Phys. Rev. Lett.103, 060403 (2009)
work page 2009
- [36]
-
[37]
R. Blankenbecler, D. J. Scalapino, and R. L. Sugar, Phys. Rev. D24, 2278 (1981)
work page 1981
- [38]
- [39]
- [40]
- [41]
-
[42]
M. Qin, H. Shi, and S. Zhang, Phys. Rev. B94, 235119 (2016)
work page 2016
-
[43]
J. Carlson, S. Gandolfi, K. E. Schmidt, and S. Zhang, Phys. Rev. A84, 061602 (2011)
work page 2011
-
[44]
H. Shi, S. Chiesa, and S. Zhang, Phys. Rev. A92, 033603 (2015)
work page 2015
-
[45]
P. Magierski, G. Wlaz lowski, A. Bulgac, and J. E. Drut, Phys. Rev. Lett.103, 210403 (2009)
work page 2009
- [46]
- [47]
-
[48]
B.-N. Lu, N. Li, S. Elhatisari, D. Lee, J. E. Drut, T. A. L¨ ahde, E. Epelbaum, and U.-G. Meißner, Phys. Rev. Lett.125, 192502 (2020)
work page 2020
-
[49]
Z. Ren, S. Elhatisari, T. A. L¨ ahde, D. Lee, and U.-G. Meißner, Phys. Lett. B850, 138463 (2024)
work page 2024
- [50]
- [51]
-
[52]
P. Cook, D. Jammooa, M. Hjorth-Jensen, D. D. Lee, and D. Lee, Nat. Commun.16, 5929 (2025)
work page 2025
-
[53]
R. Somasundaram, C. L. Armstrong, P. Giuliani, K. Godbey, S. Gandolfi, and I. Tews, Phys. Lett. B866, 139558 (2025). 9
work page 2025
- [54]
- [55]
-
[56]
D. R. Hamann and S. B. Fahy, Phys. Rev. B41, 11352 (1990)
work page 1990
-
[57]
J. E. Hirsch, Phys. Rev. B28, 4059 (1983)
work page 1983
-
[58]
P. J. Reynolds, D. M. Ceperley, B. J. Alder, and W. A. Lester, Jr., J. Chem. Phys.77, 5593 (1982)
work page 1982
-
[59]
F. Assaad, inLecture Notes of the Winter School on Quantum Simulations of Complex Many-Body Systems: From Theory to Algorithms, Vol. 10 (2002)
work page 2002
- [60]
-
[61]
F. Werner and Y. Castin, “Exact relations for quantum- mechanical few-body and many-body problems with short-range interactions in two and three dimensions,” (2010), arXiv:1001.0774 [cond-mat]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[62]
S. R. Beane, P. F. Bedaque, A. Parreno, and M. J. Sav- age, Phys. Lett. B585, 106 (2004)
work page 2004
-
[63]
P. Klos, J. E. Lynn, I. Tews, S. Gandolfi, A. Gezerlis, H.-W. Hammer, M. Hoferichter, and A. Schwenk, Phys. Rev. C94, 054005 (2016)
work page 2016
- [64]
- [65]
-
[66]
S. K¨ onig, A. Ekstr¨ om, K. Hebeler, D. Lee, and A. Schwenk, Phys. Lett. B810, 135814 (2020)
work page 2020
-
[67]
J. A. Melendez, C. Drischler, A. J. Garcia, R. J. Furn- stahl, and X. Zhang, Phys. Lett. B821, 136608 (2021)
work page 2021
-
[68]
S. Wesolowski, I. Svensson, A. Ekstr¨ om, C. Forss´ en, R. J. Furnstahl, J. A. Melendez, and D. R. Phillips, Phys. Rev. C104, 064001 (2021)
work page 2021
-
[69]
T. Dj¨ arv, A. Ekstr¨ om, C. Forss´ en, and H. T. Johansson, Phys. Rev. C105, 014005 (2022)
work page 2022
-
[70]
M. Companys Franzke, A. Tichai, K. Hebeler, and A. Schwenk, Phys. Rev. C109, 024311 (2024)
work page 2024
- [71]
-
[72]
W. G. Jiang, C. Forss´ en, T. Dj¨ arv, and G. Hagen, Phys. Rev. C109, 064314 (2024)
work page 2024
-
[73]
W. G. Jiang, C. Forss´ en, T. Dj¨ arv, and G. Hagen, Phys. Rev. C109, L061302 (2024)
work page 2024
- [74]
-
[75]
B. T. Reed, R. Somasundaram, S. De, C. L. Armstrong, P. Giuliani, C. Capano, D. A. Brown, and I. Tews, ApJ 974, 285 (2024)
work page 2024
-
[76]
Inferring three-nucleon couplings from multi-messenger neutron- star observations,
R. Somasundaram, I. Svensson, S. De, A. E. Deneris, Y. Dietz, P. Landry, A. Schwenk, and I. Tews, “Inferring three-nucleon couplings from multi-messenger neutron- star observations,” (2024), arXiv:2410.00247 [nucl-th]
-
[77]
C. L. Armstrong, P. Giuliani, K. Godbey, R. Soma- sundaram, and I. Tews, Phys. Rev. Lett.135, 142501 (2025)
work page 2025
-
[78]
E. Bonilla, P. Giuliani, K. Godbey, and D. Lee, Phys. Rev. C106, 054322 (2022)
work page 2022
-
[79]
M. J. H. Ku, A. T. Sommer, L. W. Cheuk, and M. W. Zwierlein, Science335, 563 (2012)
work page 2012
-
[80]
C. Kohstall, M. Zaccanti, M. Jag, A. Trenkwalder, P. Massignan, G. M. Bruun, F. Schreck, and R. Grimm, Nature485, 615 (2012)
work page 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.