pith. machine review for the scientific record. sign in

arxiv: 2604.24883 · v1 · submitted 2026-04-27 · ❄️ cond-mat.quant-gas · cond-mat.supr-con· cs.LG· physics.comp-ph

Recognition: unknown

Uncovering Exotic Paired States in the 2D Spin-Imbalanced Fermi Gas with Neural Wave Functions

Andres Perez Fadon, David Pfau, Gino Cassella, Halvard Sutterud, James S. Spencer, Johannes Knolle, Wan Tong Lou, W.M.C. Foulkes

Authors on Pith no claims yet

Pith reviewed 2026-05-07 17:10 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas cond-mat.supr-concs.LGphysics.comp-ph
keywords 2D Fermi gasspin imbalanceCooper pairsFulde-Ferrell-Larkin-Ovchinnikovneural wave functionsvariational Monte Carlotranslational symmetry breakingphase separation
0
0 comments X

The pith

Neural simulations of the 2D spin-imbalanced Fermi gas reveal an exotic crystal of Cooper pairs at intermediate interaction strengths.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper maps the zero-temperature phase diagram of the two-dimensional Fermi gas with unequal numbers of spin-up and spin-down particles under short-range attraction. It recovers the Fulde-Ferrell-Larkin-Ovchinnikov state at weak coupling, a polarized superfluid at strong coupling, and macroscopic phase separation at very strong coupling. At intermediate coupling the calculations instead find translational symmetry breaking: Cooper pairs arrange into a periodic crystal that coexists with a fluid of unpaired majority-spin particles. A sympathetic reader cares because this identifies a concrete microscopic realization of a long-conjectured paired crystal whose stability had remained uncertain.

Core claim

Using variational Monte Carlo with the AGPs FermiNet ansatz, the calculations observe translational symmetry breaking at intermediate interaction strengths, where the system forms an exotic crystal of Cooper pairs in a Fermi fluid of unpaired majority-spin particles. The same method recovers the Fulde-Ferrell-Larkin-Ovchinnikov phase in the weakly interacting BCS limit, a polarized superfluid in the BEC limit, and phase separation when interactions become very strong. The minority-spin momentum density is strongly suppressed inside the momentum region occupied by the unpaired majority particles once pairs are tightly bound.

What carries the argument

The AGPs FermiNet Ansatz, a neural-network variational wave function that is optimized by Monte Carlo sampling to approximate the ground state of the many-body Hamiltonian.

If this is right

  • The Fulde-Ferrell-Larkin-Ovchinnikov phase appears in the weakly interacting BCS regime.
  • A polarized superfluid is stable in the strongly interacting BEC regime.
  • Macroscopic phase separation into regions of tightly bound bosonic pairs and regions of unpaired majority particles occurs at very strong attraction.
  • The minority-spin momentum distribution develops a pronounced hole in the momentum range occupied by the unpaired majority fermions once pairs become tightly bound.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same neural-wave-function approach could be applied to three-dimensional or quasi-two-dimensional imbalanced gases to test whether the crystalline phase survives dimensional crossover.
  • Experiments that image pair correlations in two-dimensional Fermi gases at intermediate coupling could search for the predicted density modulation directly.
  • The method opens a route to study other conjectured exotic orders, such as supersolid or paired-density-wave states, in quantum gases where conventional ansatzes struggle.

Load-bearing premise

The finite-system AGPs FermiNet wave function optimized by variational Monte Carlo faithfully represents the true ground state without substantial bias from ansatz incompleteness or finite-size effects.

What would settle it

Direct experimental detection of periodic spatial modulation in the pair correlation function or in the local density at the interaction strength and population imbalance where the simulations predict the crystal phase.

Figures

Figures reproduced from arXiv: 2604.24883 by Andres Perez Fadon, David Pfau, Gino Cassella, Halvard Sutterud, James S. Spencer, Johannes Knolle, Wan Tong Lou, W.M.C. Foulkes.

Figure 1
Figure 1. Figure 1: FIG. 1: Mean-field (BCS) momentum density of the view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Momentum density of the view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Condensate fraction of the view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: One-particle density of the view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Occupied momentum states of the SIFG in our view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The spin-resolved momentum density and the view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Momentum density of the view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Momentum density of the view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Momentum density of the view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Momentum density of the view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: Momentum density of the view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Momentum density of the view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: Condensate fractions of the view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: Condensate fractions of the view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15: One-particle density of the 2D spin-imbalanced Fermi gas with view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16: One-particle density of the 2D spin-imbalanced Fermi gas with view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17: One-particle densities and their 1D projections for the view at source ↗
read the original abstract

We study the zero-temperature phase diagram of the 2D spin-imbalanced Fermi gas with short-ranged attractive interactions using the recently developed neural network variational Monte Carlo method with the AGPs FermiNet Ansatz. The Fulde-Ferrell-Larkin-Ovchinnikov phase is observed in the weakly interacting BCS limit and a polarised superfluid is seen in the strongly interacting BEC limit. When the interactions are strong, the minority-spin momentum density is reduced almost to zero in the momentum-space region occupied by the unpaired majority-spin electrons. When the interactions are very strong, phase separation occurs, with regions containing bosonic pairs and unpaired regions occupied by the remaining majority-spin particles. In addition, we observe translational symmetry breaking at intermediate interaction strengths, where the system forms an exotic crystal of Cooper pairs in a Fermi fluid of unpaired majority-spin particles. We provide a possible explanation for the formation of the crystalline phase, explain the origins of the k-space momentum-density hole when the pairs are tightly bound, and discuss how our approach opens new directions for future work.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript applies neural-network variational Monte Carlo with the AGPs FermiNet ansatz to map the zero-temperature phase diagram of the two-dimensional spin-imbalanced Fermi gas with short-range attraction. It reports the Fulde-Ferrell-Larkin-Ovchinnikov phase at weak coupling, a polarized superfluid at strong coupling, near-complete suppression of minority-spin momentum density inside the majority Fermi sea at intermediate-to-strong coupling, phase separation into bosonic-pair and unpaired-majority domains at very strong coupling, and, most notably, spontaneous translational symmetry breaking at intermediate interaction strengths that produces a crystalline lattice of Cooper pairs embedded in a Fermi sea of excess majority particles.

Significance. If the reported translational-symmetry-broken phase survives finite-size extrapolation and ansatz validation, the work would constitute a significant advance: it supplies the first microscopic evidence for an exotic pair crystal in a continuum 2D Fermi gas and demonstrates that unconstrained neural wave functions can discover unanticipated broken-symmetry states without prior bias. The method also yields concrete, falsifiable predictions for momentum distributions and pair-density modulations that could be tested in ultracold-atom experiments.

major comments (3)
  1. [results on intermediate-coupling regime] The central claim of translational symmetry breaking (abstract and results section) is obtained from a single, unconstrained AGPs FermiNet ansatz optimized at finite particle number. No comparison is presented to a translationally invariant variant of the same ansatz, nor is a finite-size extrapolation of the pair-structure-factor peak height or wave-vector reported; without these controls the observed modulation could be an optimization artifact rather than a thermodynamic-limit property.
  2. [abstract and methods] The abstract states that the variational Monte Carlo energies supply the phase boundaries, yet no quantitative statistical or systematic error bars, convergence diagnostics with respect to Monte Carlo steps or network parameters, or direct benchmarks against exact diagonalization or diffusion Monte Carlo on small systems are provided. Because VMC furnishes only an upper bound, the location of the reported phase transitions remains sensitive to possible bias in the neural ansatz.
  3. [strong-coupling momentum distributions] The reported suppression of minority-spin spectral weight inside the majority Fermi sea (strong-coupling regime) is presented as a robust feature, but the manuscript does not quantify how this feature evolves with system size or with the number of hidden units; a residual finite-size effect or ansatz incompleteness could artificially deepen the momentum-space hole.
minor comments (2)
  1. [figures] Figure captions should explicitly state the particle numbers, twist angles, and network hyperparameters used for each data point so that reproducibility is immediate.
  2. [methods] The definition of the pair-density operator and the precise procedure for extracting the modulation wave vector should be moved from the supplemental material into the main text.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive criticism of our manuscript. We address each major comment in detail below. Where the concerns identify genuine gaps, we have revised the manuscript to incorporate additional controls, diagnostics, and analysis.

read point-by-point responses
  1. Referee: [results on intermediate-coupling regime] The central claim of translational symmetry breaking (abstract and results section) is obtained from a single, unconstrained AGPs FermiNet ansatz optimized at finite particle number. No comparison is presented to a translationally invariant variant of the same ansatz, nor is a finite-size extrapolation of the pair-structure-factor peak height or wave-vector reported; without these controls the observed modulation could be an optimization artifact rather than a thermodynamic-limit property.

    Authors: We agree that these controls are necessary to establish the robustness of the translational symmetry breaking. In the revised manuscript we add a direct comparison to a translationally invariant variant of the AGPs FermiNet (constructed by enforcing uniform phase factors across particles and averaging over twisted boundary conditions). The symmetry-broken ansatz yields a lower variational energy than the invariant one at the same interaction strength. We also include finite-size extrapolations of the pair-structure-factor peak height and ordering wave-vector for system sizes from N=20 to N=100 particles. The peak height extrapolates to a nonzero value in the thermodynamic limit while the wave-vector remains stable, supporting the crystalline phase as a genuine feature rather than an artifact. revision: yes

  2. Referee: [abstract and methods] The abstract states that the variational Monte Carlo energies supply the phase boundaries, yet no quantitative statistical or systematic error bars, convergence diagnostics with respect to Monte Carlo steps or network parameters, or direct benchmarks against exact diagonalization or diffusion Monte Carlo on small systems are provided. Because VMC furnishes only an upper bound, the location of the reported phase transitions remains sensitive to possible bias in the neural ansatz.

    Authors: The referee is correct that these quantitative controls were missing from the original submission. We have revised the methods and results sections to report statistical error bars obtained from Monte Carlo sampling variance together with systematic convergence tests versus the number of Monte Carlo steps and network hyperparameters (depth, width, and hidden units). For small systems (N≤8) we now include direct benchmarks against exact diagonalization; where diffusion Monte Carlo data are available we also compare. While VMC energies remain upper bounds, the phase boundaries are determined from crossings of relative energies and from order-parameter diagnostics; we explicitly discuss the residual bias and its effect on transition locations in the revised text. revision: yes

  3. Referee: [strong-coupling momentum distributions] The reported suppression of minority-spin spectral weight inside the majority Fermi sea (strong-coupling regime) is presented as a robust feature, but the manuscript does not quantify how this feature evolves with system size or with the number of hidden units; a residual finite-size effect or ansatz incompleteness could artificially deepen the momentum-space hole.

    Authors: We acknowledge the need to demonstrate size and ansatz independence. In the revised manuscript we present the minority-spin momentum distribution for system sizes ranging from N=16 to N=64 and for several values of the number of hidden units. The suppression inside the majority Fermi sea becomes sharper with increasing system size and remains stable once the network capacity exceeds a modest threshold, indicating that the feature is physical and tied to the formation of tightly bound pairs. We have also expanded the discussion of the physical origin of the k-space hole. revision: yes

Circularity Check

0 steps flagged

No circularity: variational Monte Carlo results are independent numerical outputs

full rationale

The paper computes the phase diagram by variationally minimizing the energy of the AGPs FermiNet ansatz via neural VMC for the 2D spin-imbalanced Fermi gas. Observed features (FFLO phase, polarized superfluid, pair crystal at intermediate coupling, phase separation) are extracted from the converged wavefunction properties such as pair density modulations and momentum distributions. No equation or claim reduces a derived quantity to a fitted parameter by construction, nor imports a uniqueness theorem or ansatz choice from self-citation that forces the result. Self-references to prior FermiNet/AGP development supply the computational method but leave the physical observations as falsifiable outputs of the optimization, not tautological redefinitions. The chain is a direct numerical simulation under the variational principle.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on the variational principle, the representational power of the chosen neural ansatz, and standard assumptions of zero-temperature quantum many-body physics for short-range attractive fermions.

free parameters (2)
  • interaction strength parameter
    Varied continuously across BCS to BEC regimes to locate phase boundaries
  • spin-imbalance ratio
    Fixed or scanned to control polarization
axioms (2)
  • standard math Variational principle: the trial wave function yields an upper bound to the true ground-state energy
    Foundation of all VMC calculations
  • domain assumption The AGPs FermiNet ansatz can represent the relevant pairing and symmetry-breaking correlations with sufficient accuracy
    Core modeling choice whose validity determines the phase diagram

pith-pipeline@v0.9.0 · 5531 in / 1441 out tokens · 62124 ms · 2026-05-07T17:10:26.440651+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

62 extracted references · 10 canonical work pages

  1. [1]

    D. Pfau, J. S. Spencer, A. G. d. G. Matthews, and W. M. C. Foulkes, Physical Review Research2, 033429 (2020), arXiv:1909.02487 [physics]

  2. [2]

    W. T. Lou, H. Sutterud, G. Cassella, W. M. C. Foulkes, J. Knolle, D. Pfau, and J. S. Spencer, Physical Review X 14, 021030 (2024)

  3. [3]

    C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Re- views of Modern Physics82, 1225 (2010)

  4. [4]

    Q. Chen, J. Stajic, S. Tan, and K. Levin, Physics Reports 412, 1 (2005)

  5. [5]

    Zwerger (Springer Berlin Heidelberg, Berlin, Heidel- berg, 2012) pp

    M.Randeria, W.Zwerger,andM.Zwierlein,inThe BCS- BEC Crossover and the Unitary Fermi Gas, edited by W. Zwerger (Springer Berlin Heidelberg, Berlin, Heidel- berg, 2012) pp. 1–32

  6. [6]

    D. C. W. Foo and G. J. Conduit, Physical Review A100, 063602 (2019)

  7. [7]

    A. E. Feiguin and F. Heidrich-Meisner, Physical Review Letters102, 076403 (2009)

  8. [8]

    Moreo and D

    A. Moreo and D. J. Scalapino, Physical Review Letters 66, 946 (1991)

  9. [9]

    Y. Chen, Z. D. Wang, F. C. Zhang, and C. S. Ting, Exploring Exotic Superfluidity of Polarized Ultracold Fermions in Optical Lattices (2007), arXiv:0710.5484 [cond-mat]

  10. [10]

    He and P

    L. He and P. Zhuang, Physical Review A78, 033613 (2008)

  11. [11]

    T. K. Koponen, T. Paananen, J.-P. Martikainen, M. R. Bakhtiari, and P. Törmä, New Journal of Physics10, 045014 (2008)

  12. [12]

    M. J. Wolak, B. Grémaud, R. T. Scalettar, and G. G. Batrouni, Physical Review A86, 023630 (2012)

  13. [13]

    W. Ong, C. Cheng, I. Arakelyan, and J. E. Thomas, Physical Review Letters114, 110403 (2015)

  14. [14]

    Mitra, P

    D. Mitra, P. T. Brown, P. Schauß, S. S. Kondov, and W. S. Bakr, Physical Review Letters117, 093601 (2016)

  15. [15]

    Fulde and R

    P. Fulde and R. A. Ferrell, Physical Review135, A550 (1964)

  16. [16]

    A. I. Larkin and Y. N. Ovchinnikov, Zhurnal Eksperi- mental’noi i Teoreticheskoi Fiziki47, 1136 (1964)

  17. [17]

    Sarma, Journal of Physics and Chemistry of Solids 24, 1029 (1963)

    G. Sarma, Journal of Physics and Chemistry of Solids 24, 1029 (1963)

  18. [18]

    Randeria, J.-M

    M. Randeria, J.-M. Duan, and L.-Y. Shieh, Physical Re- view Letters62, 981 (1989)

  19. [19]

    Randeria, J.-M

    M. Randeria, J.-M. Duan, and L.-Y. Shieh, Physical Re- view B41, 327 (1990)

  20. [20]

    L. He, H. Lv, G. Cao, H. Hu, and X.-J. Liu, Physical Re- view A92, 023620 (2015), arXiv:1506.07156 [cond-mat]. 7

  21. [21]

    M. M. Parish, S. K. Baur, E. J. Mueller, and D. A. Huse, Physical Review Letters99, 250403 (2007)

  22. [22]

    G. G. Batrouni, M. H. Huntley, V. G. Rousseau, and R. T. Scalettar, Physical Review Letters100, 116405 (2008)

  23. [23]

    Taylor, A

    E. Taylor, A. Griffin, and Y. Ohashi, Physical Review A 76, 023614 (2007), arXiv:cond-mat/0703433

  24. [24]

    Potapova, I

    A. Potapova, I. Pilé, T.-C. Yi, R. Mondaini, and E. Burovski, Physical Review A107, 053301 (2023)

  25. [25]

    E. Y. Loh, J. E. Gubernatis, R. T. Scalettar, S. R. White, D. J. Scalapino, and R. L. Sugar, Physical Review B41, 9301 (1990)

  26. [26]

    Troyer and U.-J

    M. Troyer and U.-J. Wiese, Physical Review Letters94, 170201 (2005)

  27. [27]

    T.-C.Yi, S.Cheng, I.Pilé, E.Burovski,andR.Mondaini, Physical Review B110, 085131 (2024)

  28. [28]

    Carleo and M

    G. Carleo and M. Troyer, Science355, 602 (2017)

  29. [29]

    J. S. Spencer, D. Pfau, A. Botev, and W. M. C. Foulkes, Better, Faster Fermionic Neural Networks (2020), arXiv:2011.07125 [physics]

  30. [30]

    Vitali, P

    E. Vitali, P. Rosenberg, and S. Zhang, Physical Review Letters128, 203201 (2022)

  31. [31]

    Pieri and G

    P. Pieri and G. C. Strinati, Physical Review Letters96, 150404 (2006)

  32. [32]

    Y.-i.Shin, A.Schirotzek, C.H.Schunck,andW.Ketterle, Physical Review Letters101, 070404 (2008)

  33. [33]

    W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Ra- jagopal, Reviews of Modern Physics73, 33 (2001)

  34. [34]

    Gao and S

    N. Gao and S. Günnemann, inProceedings of the 40th International Conference on Machine Learning, Proceed- ings of Machine Learning Research, Vol. 202, edited by A. Krause, E. Brunskill, K. Cho, B. Engelhardt, S. Sabato, and J. Scarlett (PMLR, 2023) pp. 10708– 10726

  35. [35]

    K. Choo, A. Mezzacapo, and G. Carleo, Nature Commu- nications11, 2368 (2020)

  36. [36]

    Hermann, Z

    J. Hermann, Z. Schätzle, and F. Noé, Nature Chemistry 12, 891 (2020)

  37. [37]

    von Glehn, J

    I. von Glehn, J. S. Spencer, and D. Pfau, A Self- Attention Ansatz for Ab-initio Quantum Chemistry (2023), arXiv:2211.13672 [physics]

  38. [38]

    Cassella, H

    G. Cassella, H. Sutterud, S. Azadi, N. D. Drummond, D. Pfau, J. S. Spencer, and W. M. C. Foulkes, Physical Review Letters130, 036401 (2023)

  39. [39]

    Cassella, W

    G. Cassella, W. M. C. Foulkes, D. Pfau, and J. S. Spencer, Neural network variational Monte Carlo for positronic chemistry (2024), arXiv:2310.05607 [cond- mat, physics:physics]

  40. [40]

    Pescia, J

    G. Pescia, J. Han, A. Lovato, J. Lu, and G. Carleo, Phys- ical Review Research4, 023138 (2022)

  41. [41]

    X. Li, C. Fan, W. Ren, and J. Chen, Physical Review Research4, 013021 (2022)

  42. [42]

    Luo and B

    D. Luo and B. K. Clark, Physical Review Let- ters122, 226401 (2019), arXiv:1807.10770 [cond-mat, physics:physics]

  43. [43]

    Bajdich, L

    M. Bajdich, L. Mitas, G. Drobný, L. K. Wagner, and K. E. Schmidt, Physical Review Letters96, 130201 (2006)

  44. [44]

    Casula, C

    M. Casula, C. Attaccalite, and S. Sorella, J. Chem. Phys 121, 7110 (2004)

  45. [45]

    Casula and S

    M. Casula and S. Sorella, The Journal of Chemical Physics119, 6500 (2003), arXiv:cond-mat/0305169

  46. [46]

    Bouchaud, A

    J. Bouchaud, A. Georges, and C. Lhuillier, Journal de Physique49, 553 (1988)

  47. [47]

    Flügge,Practical Quantum Mechanics(Springer Berlin Heidelberg, Berlin, Heidelberg, 1971)

    S. Flügge,Practical Quantum Mechanics(Springer Berlin Heidelberg, Berlin, Heidelberg, 1971)

  48. [48]

    R. J. Needs, M. D. Towler, N. D. Drummond, P. L. Rios, and J. R. Trail, The Journal of Chemical Physics152,154106(2020),arXiv:2003.06506[cond-mat, physics:physics]

  49. [49]

    This gives us confidence that the results are qualitatively robust

    Although the system simulated is too small to allow us to reliably mitigate finite-size errors by twist aver- aging and/or size-extrapolation techniques, the form of the phase diagram remains the same when we reduce the (already small) system size or change the neural network architecture. This gives us confidence that the results are qualitatively robust

  50. [50]

    D. S. Petrov, M. A. Baranov, and G. V. Shlyapnikov, Physical Review A67, 031601 (2003)

  51. [51]

    Bertaina and S

    G. Bertaina and S. Giorgini, Physical Review Letters 106, 110403 (2011)

  52. [52]

    J. S. Spencer, D. Pfau, and contributors, FermiNet (2020)

  53. [53]

    Google/jax, Google (2024)

  54. [54]

    Martens and R

    J. Martens and R. Grosse, inProceedings of the 32nd In- ternational Conference on Machine Learning, Proceed- ings of Machine Learning Research, Vol. 37, edited by F. Bach and D. Blei (PMLR, Lille, France, 2015) pp. 2408–2417

  55. [55]

    Botev and J

    A. Botev and J. Martens, KFAC-JAX (2022)

  56. [56]

    Glorot and Y

    X. Glorot and Y. Bengio, inProceedings of the Thirteenth International Conference on Artificial Intelligence and Statistics(JMLRWorkshopandConferenceProceedings,

  57. [57]

    Fermi surfaces

    pp. 249–256. 8 END MA TTER Appendix IIn this section, we list the set of interac- tion strength parameter,v 0, and the corresponding val- ues ofη= ln(k F as)in Table (I). The Fermi energyϵF is v0 η|ϵ b|/ϵF Region 0.10 4.367 4.061×10 −4 (I)0.15 1.956 5.047×10 −2 0.20 0.747 5.661×10 −1 0.25 0.018 2.431 0.30 -0.470 6.457 (II)0.35 -0.821 13.04 0.40 -1.087 22....

  58. [58]

    ˆψα(r1) E (15) ρ(2) αβ(r′ 1,r ′ 2;r 1,r 2) = D ˆψ† α(r′

  59. [59]

    ˆψβ(r2) ˆψα(r1) E (16) 11 Parameter Value Batch size 4096 Training iterations 3e5 Pretraining iterations None Learning rate(2e4 +t) −1 Local energy clipping 5.0 KFAC Momentum 0 KFAC Covariance moving average decay 0.95 KFAC Norm constraint 1e-3 KFAC Damping 1e-3 MCMC Proposal std. dev. (per dimension) 0.02 MCMC Steps between parameter updates 10 TABLE II:...

  60. [60]

    drN R |Ψ(R)|2dR ,(23) ρ(2) αβ(r1,r 2;r ′ 1,r ′

    =N α R |Ψ(R)|2 Ψ(r′ 1) Ψ(r1) dr2 . . . drN R |Ψ(R)|2dR ,(23) ρ(2) αβ(r1,r 2;r ′ 1,r ′

  61. [61]

    drN R |Ψ(R)|2dR ,(24) whereαandβdenote the spin species

    =N α(Nβ −δ αβ) R |Ψ(R)|2 Ψ(r′ 1,r′ 2) Ψ(r1,r2) dr3 . . . drN R |Ψ(R)|2dR ,(24) whereαandβdenote the spin species. To reduce the effect of the one-body contribution to the TBDM, we use an improved estimator of the TBDM in Eq. (25) that removes the one-body contribution explicitly [48]: ρ(2) αβ(r1,r 2;r ′ 1,r ′

  62. [62]

    drN R |Ψ(R)|2dR .(25) All of these quantities can then be estimated using Monte Carlo sampling

    =N α(Nβ −δ αβ) R |Ψ(R)|2 h Ψ(r′ 1,r′ 2) Ψ(r1,r1) − Ψ(r′ 1,r2) Ψ(r1,r2) Ψ(r1,r′ 2) Ψ(r1,r2) i dr3 . . . drN R |Ψ(R)|2dR .(25) All of these quantities can then be estimated using Monte Carlo sampling. All Data In this section, we present all the results of our simulations, including the momentum densities, the condensate fraction and the one-particle real-s...