{"id":"f9661368-193c-45c0-ae21-d81808716888","arxiv_id":"2506.20891","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Numerically computed spectra, momentum distributions, and pair correlations for 2D trapped Fermi gases with up to six atoms, extended to finite effective range.","lead":"The paper computes numerically converged energy spectra and structural properties for two-dimensional Fermi gases with two, four, and six atoms in a harmonic trap, including the effect of a finite effective range. It provides reference predictions for mesoscopic cold-atom experiments and checks how far such small systems are from the regime of tightly bound molecules.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quasi-2D mapping (Eq. 5) is applied at r2D/lr^2=-0.2, where lz≈0.54lr; at the trap momenta probed (k~1/lr), klz is not small, so higher-order confinement corrections could be comparable to the effective-range term and are not quantified.","rationale":"The reader's weakest assumption was that the two-body physics is fully captured by a2D and r2D with higher-order terms in the effective-range expansion negligible. My concern is the same assumption, sharpened to the quasi-2D mapping Eq. (5), which is the one place where the paper invokes a specific external relation rather than a direct numerical fit. The numbers make the issue concrete: at r2D/lr^2=-0.2, Eq. (5) gives lz=0.54lr, so at the highest energies shown (eps_b=2 hbar omega_r, k~1.4/lr) the condition klz<<1 is violated by a factor near unity. The paper's own text in Section 2 states the mapping applies 'when lz is small (such that klz << 1)', yet the plotted range includes values where klz is not small. This is an internal consistency issue, not merely a disagreement with consensus, and it directly affects the quantitative strength of the headline claim that pairing is enhanced at fixed binding energy by finite negative effective range. I considered whether numerical convergence of the 3+3 excited-state manifold is a more load-bearing concern, but the paper provides a benchmark against the Busch spectrum for 1+1, an analytic non-interacting check for 2+2, and transparent statements about basis-size limitations; those are standard and do not single out the central claim. Similarly, the ambiguous condensate fraction is explicitly flagged by the authors. The effective-range truncation is the least secured external input, and it is testable. Because the reader already assigned CONDITIONAL and my concern reinforces that condition rather than overturning the paper's central 2D results, I recommend no change to the verdict.","tokens_in":19237,"tokens_out":10241,"duration_ms":117002,"concrete_test":"Solve the two-body problem in a true quasi-2D geometry: take the same short-range model potential in three dimensions with an axial harmonic confinement of length lz = 0.54 lr, and compute the trapped relative-energy spectrum and ground-state pair correlation for eps_b up to 2 hbar omega_r. Compare these against the 2D model using r2D = -0.2 lr^2 and Eq. (5). If the level energies or the small-r peak of the pair distribution function differ by more than a few percent from the 2D effective-range prediction, the finite-negative-range results in Section 3.5 are not a faithful representation of the quasi-2D system in the regime claimed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that a finite negative effective range enhances pairing at fixed binding energy depends on the quantitative validity of the mapping r2D=-lz^2 ln(2) (Eq. 5) and on the assertion in Section 2 that r0<=0.1lr makes higher-order terms in Eq. (4) negligible. The r0 condition controls only the short-range shape of the 2D model potential; it does not test the accuracy of the quasi-2D projection. For the largest effective range used, r2D/lr^2=-0.2, Eq. (5) gives lz/lr = sqrt(0.2/ln2) ≈ 0.54. The derivation of Eq. (5) assumes klz << 1, but the energy spectra and structural results extend to eps_b ≈ 2 hbar omega_r, corresponding to relative momenta k ~ sqrt(2)/lr and hence klz ≈ 0.76. Corrections of order (klz)^2 ≈ 0.6 can be comparable to, or larger than, the effective-range term (1/pi)k^2 r2D ≈ -0.06, and the next-order shape parameter is not computed. Thus the comparison between r2D≈0 and r2D=-0.2 at fixed eps_b may partly reflect truncation error in the quasi-2D mapping rather than a universal effective-range effect. The paper states that the mapping is used while remaining within its regime of validity, but it does not quantify the error at the highest energies shown.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies equal-mass spin-balanced two-component Fermi gases with two, four, and six atoms confined in a two-dimensional harmonic trap, interacting via a short-range Gaussian potential parameterized by the 2D scattering length and effective range. Using the explicitly correlated Gaussian method with stochastic variational optimization, the authors compute ground- and low-lying monopole energy spectra as functions of the two-body binding energy, natural-orbital occupation numbers from one- and reduced two-body density matrices, a molecular condensate fraction, momentum distributions of atoms and pairs, and radial and pair distribution functions. Benchmarks are provided against the Busch spectrum for 1+1 fermions and against analytical results for the non-interacting 2+2 system in Appendix A. The main physical conclusions are that for up to six atoms at binding energies eps_b <~ 2 hbar omega_r the gas remains outside the tightly bound molecular regime, and that at fixed binding energy a finite negative effective range, representing quasi-2D confinement, enhances pairing.","tokens_in":19564,"tokens_out":5278,"duration_ms":57194,"significance":"If the results hold, this paper provides a useful two-dimensional counterpart to the established three-dimensional few-fermion studies of Blume and co-workers, with direct relevance to ongoing mesoscopic Fermi-gas experiments. The strengths are the use of an established numerical method, explicit checks against exact analytic limits, self-contained calculations with no fitted observables, and an analytic Fourier transform of the density matrices that makes momentum-space observables cheap to evaluate. The paper is also candid about the limitations of the ECG approach for larger atom numbers and deeper binding. The main uncertainty is the quantitative validity of the quasi-2D effective-range mapping at the largest binding energies considered, which is load-bearing for the central claim about enhanced pairing at fixed binding energy; this needs to be addressed before the finite-effective-range results can be taken as universal rather than as artifacts of truncation of the mapping.","major_comments":[{"comment":"The central comparison between r2D/lr^2 = -0.2 and approximately zero is made without quantifying the error of the quasi-2D mapping in Eq. (5). This mapping is derived for klz << 1, but at the largest binding energies plotted, eps_b ~ 2 hbar omega_r, the relative momentum is k ~ sqrt(2)/lr, so with lz/lr = sqrt(0.2/ln 2) ~ 0.54 one obtains klz ~ 0.76. Higher-order confinement corrections of order (klz)^2 ~ 0.6 are then comparable to the effective-range term (1/pi)k^2 r2D ~ -0.06, and the next-order shape parameter is not computed. The condition r0 <= 0.1 lr controls only the short-range shape of the two-dimensional model potential and does not test the accuracy of the quasi-2D projection. I therefore request an explicit estimate of the omitted confinement corrections, or a restriction of the finite-range comparison to binding energies where klz is demonstrably small.","section":"Section 2, Eq. (5); Section 3.5, Figs. 6-9"},{"comment":"The manuscript reports no convergence residuals, basis-size dependence, or error bars for the energies or the derived structural observables. The text states that convergence cannot currently be achieved for six atoms at eps_b > 2 hbar omega_r, yet the spectra and structural results extend to eps_b ~ 2.1 hbar omega_r, and the central conclusion that eps_b <~ 2 hbar omega_r lies outside the strong-interaction regime relies on data in this region. Please add a quantitative convergence statement for the 3+3 states near the upper end of the range (for example, energy change with basis size or an extrapolated value) and indicate how the uncertainty propagates to the occupation numbers, momentum distributions, and distribution functions.","section":"Section 3.1 and Section 4"}],"minor_comments":[{"comment":"The displayed equality [n_up(k)]_AA' == [n_up(k)]_AA' is a typo; the left-hand side should presumably be defined as the matrix element and the right-hand side as the closed-form expression c1/(g1 g4) exp(k^2/(2 g4)).","section":"Eq. (31)"},{"comment":"The text states that V(rk - rl) = delta(r - rk - rl) is substituted into Eq. (35b), but the argument of the delta function should be r - (rk - rl); please correct the sign.","section":"Eq. (35b)"},{"comment":"The definition of the condensate fraction is not fully explicit about the set over which the maximum is taken: Eq. (22a) includes the n=1,m=0 term while the surrounding text says the sum applies for m > 0. Please state the complete summation rule.","section":"Section 3.2.3, Eqs. (21)-(22)"},{"comment":"The caption reads 'at very nearly zero effective' and appears to be missing the word 'range'; please correct.","section":"Fig. 1 caption"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this is a careful and useful numerical reference for 2D trapped few-fermion gases, but the finite-effective-range comparison that backs the 'pairing enhancement' claim is the least solid part of the paper.\n\nWhat's actually new: the 2D energy spectra for 2+2 and 3+3, the natural-orbital occupations, momentum distributions, and pair-correlation functions at near-zero effective range, plus the finite-range trends in Sec. 3.5 and the analytic non-interacting 2+2 results in Appendix A. The benchmarks are genuinely good: the 1+1 spectrum is checked against Busch, and the 2+2 non-interacting occupations are derived analytically. The method (ECG/SVM) is established, and the authors are honest about convergence limits and about the condensate fraction being unreliable at these particle numbers.\n\nThe soft spots are real but not fatal. There are no numerical error bars or convergence residuals anywhere, so the reader has to trust the claims of convergence. For a reference paper that's a bigger issue than the authors acknowledge. The larger concern is the quasi-2D mapping. Equation (5), r2D = -lz^2 ln2, is derived in the limit klz << 1. At r2D/lr^2 = -0.2, lz/lr ≈ 0.54, and at binding energies up to 2 hbar omega_r the relative momentum is k ~ sqrt(2)/lr, so klz ≈ 0.76. The paper says it stays within the regime of validity of the mapping but gives no error estimate. The central claim that a finite negative effective range enhances pairing at fixed binding energy rests on comparing r2D ≈ 0 with r2D = -0.2, so that claim is suggestive, not proven. This is not a reason to desk-reject, but the authors should either extend the mapping to next order or soften the conclusion.\n\nThe strictly-2D results (r2D ≈ 0) stand on their own and are the bulk of the paper. The paper is written honestly, the literature is handled fairly, and the analytic appendix is a nice touch. The lack of code or data is a practical drawback for a paper aimed to be a reference.\n\nWho this is for: cold-atom experimentalists working with few fermions in 2D, and theorists who want benchmarks for few-body calculations. It deserves a serious referee. I'd send it to review with a request for convergence data and a reworked discussion of the finite-range mapping.","headline":"Careful and useful 2D few-fermion reference data, but the finite-effective-range pairing claim rests on a quasi-2D mapping whose error is never quantified.","tokens_in":20146,"tokens_out":2996,"would_cite":true,"duration_ms":28834,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper provides numerically exact energy spectra and pairing correlations for two, four, and six trapped fermions in two dimensions, showing that binding energies up to $2\\hbar\\omega_r$ lie outside the regime of tightly bound…","keywords":["mesoscopic Fermi gases","two-dimensional Fermi gas","quasi-2D confinement","effective range","correlated Gaussian basis","stochastic variational method","pairing correlations","few-body physics"],"falsifier":"Compute the same few-body observables for two different short-range potentials that share the same 2D scattering length and effective range; if the energies or pair distributions differ systematically, the effective-range truncation used here is not sufficient.","tokens_in":19029,"feed_emoji":"⚛️","tokens_out":9693,"duration_ms":101578,"temperature":0.7,"pith_summary":"This paper asks what happens to small, spin-balanced Fermi gases in a two-dimensional harmonic trap as the attraction between opposite spins grows. Using a numerically exact stochastic variational calculation on a correlated Gaussian basis, it computes the ground- and low-lying excited-state energies, natural-orbital occupations, momentum distributions, and radial correlation functions for gases containing one, two, or three atoms per spin state, parameterised by the 2D scattering length and effective range. The central result is that for up to six particles at zero effective range, binding energies up to about $2\\hbar\\omega_r$ never reach the regime of tightly bound bosonic molecules. At a fixed binding energy, moving from strictly 2D to quasi-2D confinement, which corresponds to a finite negative effective range, enhances the pairing correlations. This gives a theoretical map of few-fermion physics that modern cold-atom experiments can realise and test.","feed_headline":"Six 2D trapped fermions stay outside the molecule regime","feed_subtitle":"Numerically exact spectra and pair correlations show quasi-2D confinement, not stronger attraction, promotes pairing.","key_machinery":"The machinery is the two-dimensional effective-range expansion, $\\cot[\\delta(k)] = (2/\\pi)[\\gamma + \\ln(ka_{2D}/2)] + (1/\\pi) k^2 r_{2D} + \\cdots$, which reduces all short-range interaction details to the scattering length $a_{2D}$ and effective range $r_{2D}$. A finite-range Gaussian potential is tuned so that its width stays small, $r_0 \\lesssim 0.1 l_r$, making higher-order terms negligible, and quasi-2D confinement is encoded through the mapping $r_{2D} = -l_z^2 \\ln(2)$. The $N$-body Schrödinger equation is then solved by expanding the relative wave function in explicitly correlated Gaussian basis states with stochastically optimised width parameters, which allows numerically exact extraction of density matrices and their analytical Fourier transforms.","core_discovery":"On the paper's own terms, the discovery is a complete numerical characterisation of trapped two-dimensional Fermi gases with $1+1$, $2+2$, and $3+3$ particles, including low-lying monopole spectra, one- and two-body density matrices, natural-orbital occupation numbers, atom and pair momentum distributions, and radial pair distribution functions. The key finding is that even at the strongest binding energies converged, $\\epsilon_b \\approx 2\\hbar\\omega_r$, the opposite-spin pair correlations have not developed the tight composite-boson structure of the deep BEC limit: the pair distribution function still shows a broad dimer-dimer peak, and the one-body density matrix can be reconstructed from only a handful of natural orbitals. Adding a finite, negative effective range, the signature of quasi-2D confinement, increases the short-distance molecular peak and lowers the occupation of the lowest natural orbitals at fixed $\\epsilon_b$, so pairing is stronger in quasi-2D geometry. The paper also reports that the molecular \"condensate fraction\" extracted from the reduced two-body density matrix is non-monotonic in $\\epsilon_b$ and has limited interpretability for such small systems.","pith_inferences":["A direct experimental test would be to measure the pair momentum distribution $n(K)$ and the real-space pair distribution in a quasi-2D microtrap with two to twelve atoms while varying the axial confinement; the paper's results predict that $n(K)$ changes very little with binding energy, but the short-distance pair peak grows as the effective range becomes more negative.","Because the correlated Gaussian approach fails to converge for $\\epsilon_b > 2\\hbar\\omega_r$ in the six-fermion case, the paper leaves open whether the non-monotonic condensate fraction turns over and approaches unity at stronger coupling; alternative few-body methods could probe this part of the crossover.","The strong effect of effective range at fixed binding energy suggests a practical tuning route for experiments: changing the trap aspect ratio may mimic changing the interaction strength even when the scattering length is held fixed."],"forward_implications":["The energy spectra for $1+1$, $2+2$, and $3+3$ fermions provide finite-system benchmarks that can be compared directly with deterministic few-atom experiments and with extrapolations toward many-body 2D Fermi gas theories.","For binding energies $\\epsilon_b \\lesssim 2\\hbar\\omega_r$, the one-body density matrix is accurately decomposed into at most the six lowest natural orbitals, so few-body pairing in this regime is describable by a small effective single-particle Hilbert space.","At fixed $\\epsilon_b$, a more negative effective range, corresponding to stronger quasi-2D confinement, shifts energies upward, strengthens the short-distance peak of the pair distribution function, and lowers the lowest natural-orbital occupations, meaning confinement promotes molecule formation without changing the momentum distributions.","The non-monotonic \"condensate fraction\" defined from the reduced two-body density matrix is not a reliable order parameter for these small systems, so molecular condensation in 2D few-fermion gases should be diagnosed through pair distribution functions and pair momentum distributions instead.","Because the method cannot converge for $3+3$ fermions at binding energies above about $2\\hbar\\omega_r$, the deep BEC side of the crossover remains an open target for other few-body techniques."],"supporting_citations":[{"why":"Supplies the correlated-Gaussian stochastic-variational method, the effective-range convention, and the earlier spectra and pair-correlation results this paper extends.","marker":"[11]"},{"why":"Its three-dimensional treatment of up to six trapped fermions provides the reduced two-body density matrix, natural-orbital occupation numbers, condensate fraction, and momentum-distribution framework adapted here to two dimensions.","marker":"[25]"},{"why":"The exact analytic two-atom harmonic-trap spectrum is used to validate the $1+1$ energy spectrum shown in Fig. 1(a).","marker":"[35]"},{"why":"Gives the two-dimensional effective-range expansion (Eq. 4) that defines the scattering length $a_{2D}$ and effective range $r_{2D}$ used to parameterise the model.","marker":"[18-20]"},{"why":"Establishes the quasi-2D confinement mapping $r_{2D} = -l_z^2 \\ln(2)$ and supports the finite-range Gaussian model with a tunable shape resonance.","marker":"[21-24]"},{"why":"The same finite-range Gaussian potential and effective-range framework was previously used to model breathing modes and shape resonances in two-dimensional Fermi gases.","marker":"[24]"},{"why":"Supplies the structural-property framework and the pair-distribution integrations giving the approximate $1/2$ and $1/3$ molecular fractions for the $2+2$ and $3+3$ systems.","marker":"[26]"},{"why":"Provides the two-dimensional BCS-BEC crossover classification, locating the strong-coupling regime around $\\ln(k_F a_{2D})=0$, which is used to interpret the few-body results.","marker":"[32,33]"}],"fun_headline_variants":["Quasi-2D geometry enhances pairing more than attraction","No deep molecules: 2D fermions pair without BEC limit","Even strong binding keeps 2D fermion pairs far from molecules","Six trapped fermions: pairing dominated by geometry, not attraction","Trapped 2D fermions: pairing without deep molecules"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculations assume that the short-range interaction is fully described by the 2D scattering length and effective range, so that the detailed shape of the interatomic potential never matters at the trap energies studied.","fun_headline_variants_meta":{"raw":{"variants":["Quasi-2D geometry enhances pairing more than attraction","No deep molecules: 2D fermions pair without BEC limit","Even strong binding keeps 2D fermion pairs far from molecules","Six trapped fermions: pairing dominated by geometry, not attraction","Trapped 2D fermions: pairing without deep molecules"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000905,"raw_usage":{"total_tokens":3894,"prompt_tokens":949,"completion_tokens":2945,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":2857}},"tokens_in":565,"tokens_out":2945,"duration_ms":22188,"temperature":1.0,"reasoning_tokens":2857,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:39:09.944777+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same few-body observables for two different short-range potentials that share the same 2D scattering length and effective range; if the energies or pair distributions differ systematically, the effective-range truncation used here is not sufficient.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the correlated-Gaussian stochastic-variational method, the effective-range convention, and the earlier spectra and pair-correlation results this paper extends."},{"cited_title":"Blume and K","cited_arxiv_id":null,"evidence_quote":"Its three-dimensional treatment of up to six trapped fermions provides the reduced two-body density matrix, natural-orbital occupation numbers, condensate fraction, and momentum-distribution framework adapted here to two dimensions."},{"cited_title":"Busch, B.-G","cited_arxiv_id":null,"evidence_quote":"The exact analytic two-atom harmonic-trap spectrum is used to validate the $1+1$ energy spectrum shown in Fig. 1(a)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The same finite-range Gaussian potential and effective-range framework was previously used to model breathing modes and shape resonances in two-dimensional Fermi gases."},{"cited_title":"von Stecher, C","cited_arxiv_id":null,"evidence_quote":"Supplies the structural-property framework and the pair-distribution integrations giving the approximate $1/2$ and $1/3$ molecular fractions for the $2+2$ and $3+3$ systems."}],"review_version":1}