{"id":"116afa09-2f56-443c-9238-a810b5d1a092","arxiv_id":"1908.02849","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Energies, densities, and correlations of 1D trapped Fermi polarons are obtained at any interaction strength by diagonalizing a variational basis built from g=0 and g=∞ states.","lead":"A method for trapped one-dimensional quantum gases builds a wavefunction basis from states at zero and infinite interaction strength, then diagonalizes to get energies and densities at any repulsion. It is checked against the exact two-particle solution and matrix product state simulations, agreeing well for small particle numbers.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The g=∞ input states are under-specified: for N>1 and the double well, the coefficients a_μ satisfying Eq. (7) are never stated, so every reported energy/density depends on an unverified external construction.","rationale":"I agree with the reader's weakest-assumption identification. The method's selling point is exactness at both limits; if the g=∞ basis is not the true one, the interpolation is built on sand. The reader's other concerns (no code/data, 6+1 MPS discrepancy, no error bars) are real but secondary: the 6+1 discrepancy is a convergence issue that could be cured by a larger basis, and missing code is a reproducibility issue. The a_μ specification is the one place where a wrong input would invalidate the central construction itself. The paper does give explicit a_μ for N=1 and the 1+1 agreement with the Busch solution is strong independent support, but that does not extend to N>1. The 2+1 harmonic-well agreement with MPS suggests that the authors' a_μ are probably correct for that case, but the double-well results depend on a different set of exchange integrals, and no derivation is shown. A focused recomputation of the double-well cases would settle whether this is a documentation lapse or a substantive error. Since the concern is concrete and addressable, and since the published evidence is otherwise consistent, the conditional verdict stands unchanged.","tokens_in":24577,"tokens_out":22950,"duration_ms":279948,"concrete_test":"Independently compute the a_μ coefficients for the double-well potential of Appendix C using the spin-chain construction of Refs [37,38]: build the single-particle orbitals, evaluate the exchange integrals, diagonalize the effective spin Hamiltonian, and insert the resulting a_μ into the Gram-Schmidt basis of Section 3.4. Then repeat the 2+1 and 6+1 ground-state calculations at g=1 and g=10 and compare with the published MPS-based densities (Figures 9, 12, 13) and energies (Figures 7-8). Reproduction of the published curves would reduce the concern to a documentation gap; failure to reproduce them would show that the variational input is incorrect and the central claim is unsupported for the double-well geometry.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the premise that the variational basis reproduces the exact limits g=0 and g=∞. Section 3.2 defines the g=∞ states as ψμ=a_l^{(μ)}Φ_q in each ordering sector, with coefficients constrained only by the normalization/orthogonality condition (7), and refers to Refs [37,38] for their exactness. But for N>1 and for the double-well geometry used in Sections 5.2 and 5.3, the paper never states how the vectors a_μ are computed. For N=1 the coefficients are given explicitly (the [1,-1] state), but the 2+1 and 6+1 calculations require sets of orthogonal a_μ for 3 and 7 sectors; these are not tabulated, derived, or even described. The construction in Refs [37,38] is trap-dependent because the effective spin-chain exchange integrals depend on the single-particle orbitals of the confining potential, so the a_μ for the harmonic well and for the smooth double well of Appendix C are not the same object. If the authors' unpublished a_μ are not the correct strong-coupling eigenstates for those traps, the basis does not actually contain the g=∞ limits, and every finite-g energy, density, and correlation function inherits the error. This is therefore a correctness risk for the double-well examples, not merely a missing code/data statement. The 1+1 benchmark cannot settle it, since that is the one case where the coefficients are fully specified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a variational method for a single impurity interacting through a contact potential with N identical fermions in a one-dimensional trap. The variational basis is constructed from low-energy eigenstates of the noninteracting Hamiltonian and low-energy exact eigenstates at infinite repulsion; the two sets are Gram-Schmidt orthonormalized, and energies, density matrices, momentum distributions, and correlation functions are obtained by diagonalizing the Hamiltonian in this basis. The method is benchmarked against the exact two-particle solution in a harmonic trap and against matrix-product-state calculations for 2+1 and 6+1 systems in both harmonic and double-well geometries. The central claim is that this truncated interpolatory basis accurately describes the system at arbitrary finite repulsion while being exact at the zero- and infinite-coupling limits.","tokens_in":24815,"tokens_out":17795,"duration_ms":204766,"significance":"If the central claim holds, the method is a valuable and computationally cheap tool for trapped one-dimensional Fermi polarons: it requires a single construction of the basis for all interaction strengths, and it is exact at both limits. The paper's strengths are the detailed derivations of overlaps and density-matrix elements in Section 3.3 and Appendices A and E, the excellent 1+1 benchmark with relative energy errors near 10^-10 at roughly 70 basis states, and the independent MPS comparisons for 2+1 and 6+1. The main weakness is that the infinite-coupling input states are not fully specified for N>1 and for the double-well geometry, which undermines the reproducibility of the strong-coupling results and the claimed exactness at g=∞. With the missing coefficients supplied, the method would constitute a useful contribution to few-body quantum-gas techniques.","major_comments":[{"comment":"The infinite-interaction basis states ψ_μ of Eq. (8) are under-specified for N>1. The coefficient vectors a_μ are constrained only by the orthonormality conditions in Eq. (7), but for the basis to contain the exact strong-coupling eigenstates these vectors must be eigenstates of the trap-dependent spin-chain Hamiltonian. The paper never states how the a_μ are obtained for the harmonic 2+1 and 6+1 calculations, and in particular for the double-well potential of Appendix C. For N=1 the coefficients are given explicitly ([1,−1] in §5.1.1), but for N=2 and N=6 the basis sizes reported in §5.2.1 and §5.3 presuppose definite a_μ that are neither tabulated nor described. The energy cutoff for the infinite-interaction states is also not defined for N>1, since it is not stated whether the spin-chain energy is included. Because every reported energy, density, and correlation function for the 2+1 and 6+1 systems inherits any error in this unspecified input, the central claim of exactness at the infinite-coupling limit is not verifiable for those systems. The authors should provide the a_μ, or a complete algorithm for computing them, for each geometry and particle number, and specify exactly how the basis cutoff is applied.","section":"§3.2 and §5.2–5.3"},{"comment":"The abstract claims the method describes the system accurately at arbitrary finite repulsion, but for N>1 the independent MPS comparisons are limited to g=1 for the 2+1 system and to g=1 and g=10 for the 6+1 system (Figures 7–13). No independent check is provided at weak coupling or at intermediate coupling where the interpolation error is typically largest, and the 6+1 agreement is visibly worse without a quantified discrepancy. At minimum, the paper should add one or two additional coupling points for the 2+1 or 6+1 system, or explicitly qualify the 'arbitrary finite repulsion' claim as an expectation from the variational construction rather than a demonstrated benchmark.","section":"§5.1.2 and §5.3"}],"minor_comments":[{"comment":"There is a typo: 'Gram-Schmidth ortonormalization' should be 'Gram-Schmidt orthonormalization'.","section":"§1"},{"comment":"The text says the extrapolated MPS value is represented by a black cross, while the caption of Figure 7 says a black star; the figure and caption should be made consistent.","section":"§5.2.1"},{"comment":"The word 'discrepeancy' should be 'discrepancy'.","section":"§5.2.3"},{"comment":"The sentence containing 'he factor 1/k' contains a typo: 'he' should be 'the'.","section":"Appendix A.1"},{"comment":"The notation for the energy cutoff, e.g., 'an energy not greater than 2 (above the lowest energy state)', should be defined once in the main text rather than only in the 1+1 example, because the cutoff is the single convergence parameter of the method.","section":"§5.1.1"}],"recommendation":"major_revision","confidential_remarks":"The central issue is the missing specification of the strong-coupling coefficient vectors a_μ for N>1 and for the double-well trap. This is fixable by adding an appendix or supplementary material, but without it the 2+1 and 6+1 results are not reproducible. The reliance on Refs. [37,38] is especially problematic for the double-well case, which is not covered in those references, and the overlapping authorship with Ref. [38] makes the omission more consequential. I recommend major revision rather than rejection because the 1+1 benchmark and the independent MPS comparisons indicate the underlying method is sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this is a solid method paper that does what its title says—turns the two-state interpolatory ansatz of Ref. [45] into a systematic truncated-basis method for 1D Fermi polarons, with clean derivations and good benchmarks. The main flaw is not the physics but the reproducibility: the infinite-interaction basis states are under-specified for N>1 and for the double well.\n\nWhat is genuinely new: the construction of a Gram-Schmidt orthonormalized basis from energy-truncated g=0 and g=∞ states, the recursive formulas for Hamiltonian matrix elements (Section 3.5), and the density-matrix, momentum-distribution, and correlation-function formulas in Section 4. Those go beyond the two-state paper and make the method usable. The 1+1 harmonic-trap benchmark is excellent—relative energy errors around 1e-10 with ~70 basis states—and the 2+1 comparisons with independent MPS simulations are convincing for both the harmonic and double-well geometries. The method is lightweight, gives all g from one basis construction, and gets the limits exactly right. That is a real contribution.\n\nSoft spots: the biggest one is exactly what your stress-test flags. Equation (7) only sets normalization/orthogonality constraints on the g=∞ coefficients a_μ; the actual values are never given for N>1 or for the double well. The 1+1 case is fine because [1,-1] is explicit, and the 2+1 MPS agreement suggests the authors did compute them correctly. But a reader cannot reproduce the basis for the 6+1 or double-well results without going back through Refs. [37,38] and re-deriving the trap-dependent spin-chain eigenstates. That is a reproducibility gap, not evidence of wrong numbers, but it should be fixed before publication: state how a_μ are obtained, and ideally tabulate them or ship code.\n\nTwo minor issues: no code or data is shipped, so convergence plots can't be independently checked; and the 6+1 density agreement is visibly worse, which the paper honestly acknowledges. The abstract's 'accurately describe' is a bit strong for larger N, but the method still captures the physics.\n\nThis is a paper I would send to a serious referee. The derivations are checkable, the benchmarks are meaningful, and the method is useful for people working on few-body 1D systems. The fixes are concrete: specify the a_μ construction, add reproducibility material, and calibrate the claims for large N. None of that shakes the central construction.","headline":"A genuinely useful generalization of the two-state interpolatory ansatz with solid benchmarks; the main gap is that the g=∞ basis coefficients are underspecified for N>1 and the double well.","tokens_in":25440,"tokens_out":4474,"would_cite":true,"duration_ms":45414,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A truncated basis built from the zero- and infinite-repulsion limits of a trapped one-dimensional Fermi polaron reproduces the exact energies and densities at all repulsion strengths when the Hamiltonian is diagonalized in that basis.","keywords":["Fermi polaron","one-dimensional quantum gas","contact interaction","variational method","interpolatory ansatz","trapped ultracold gases","strong interactions","momentum distribution"],"falsifier":"Independently compute the double-well $g=1$ and $g=10$ ground-state energies and densities with a method that does not use the paper's infinite-coupling states—for example, a lattice or tensor-network simulation with extrapolation—and compare; deviations larger than the matrix-product-state discrepancies shown in the paper would falsify the claimed accuracy for that geometry. A more direct check is to solve the strong-coupling eigenstate conditions (Eq. (7)) for the double-well trap numerically and compare the resulting coefficients $a_\\mu$ with the ones implicit in the paper's basis: a mismatch would show the reported results rest on an unverified input.","tokens_in":24305,"feed_emoji":"⚛️","tokens_out":8485,"duration_ms":86525,"temperature":0.7,"pith_summary":"The paper proposes a variational method for a single impurity moving among identical fermions in a one-dimensional trap with a contact repulsion. It builds a truncated basis from the low-energy exact eigenstates at zero coupling and at infinite coupling, orthonormalizes the union, and diagonalizes the Hamiltonian once to obtain energies and observables at any coupling $g$; the diagonalization is exact at the two limits by construction. The central claim is that this interpolation stays accurate across the whole range of repulsions, and the paper supports it with the analytic two-particle harmonic solution and with matrix-product-state numerics for larger systems. If the claim holds, strongly interacting few-body problems that usually need a fresh expensive simulation at every interaction strength can instead be solved from one small, fixed basis.","feed_headline":"Two-limit basis handles 1D polaron at any repulsion","feed_subtitle":"Building one basis from zero- and infinite-coupling states yields exact endpoints and accurate in-between energies.","key_machinery":"The central object is the Gram–Schmidt-orthonormalized set $\\{|\\varphi_i\\rangle,\\,|\\chi_\\mu\\rangle\\}$ built from zero-coupling states $|\\varphi_i\\rangle$ and infinite-coupling eigenstates $|\\psi_\\mu\\rangle$, whose coordinate form is given by Eq. (8) with coefficients $a_\\mu$ satisfying the orthogonality conditions of Eq. (7). The overlaps between the two families are evaluated through determinant identities (Appendix A) that convert integrals over ordered spatial regions into derivatives of determinants, which lets the matrix elements of the non-interacting Hamiltonian and the contact interaction be computed recursively. The machinery's job is to feed exact endpoint information into the finite-coupling regime with one basis construction, after which a single numerical diagonalization produces the spectrum and all observables.","core_discovery":"On its own terms, the paper establishes that the static properties of the trapped one-dimensional Fermi polaron—energies, spatial and momentum densities, and the minority–majority correlation function—can be recovered by diagonalizing the Hamiltonian in a finite basis obtained by Gram–Schmidt orthonormalizing the union of low-energy zero-interaction states and low-energy infinite-interaction eigenstates. The construction is exact at $g\\to 0$ and $g\\to+\\infty$ by design, and the paper shows numerically that it remains accurate at intermediate repulsion: relative energy errors near $10^{-10}$ for one majority fermion plus impurity in a harmonic trap, close agreement with matrix-product-state results for two majority fermions in both harmonic and double-well traps, and clearly larger but still reasonable agreement for six majority fermions.","pith_inferences":["If the infinite-coupling coefficients for an arbitrary trap could be supplied by a separate numerical solver, this construction would become a general few-body variational solver whose accuracy is set entirely by the endpoint inputs; the paper leaves that automation implicit.","The same interpolation logic should transfer to attractive contact interactions, provided suitable strong-coupling endpoint states are used, so a testable extension is to check the convergence of the method for negative $g$ across the crossover.","Because the basis is built once, entire observables such as momentum distributions can be produced pointwise across the crossover at almost no extra cost, which suggests the method could serve as a fast engine for checking sum rules or polaron crossover curves.","The worse agreement at 6+1 indicates the selection of which endpoint states to include matters as the particle number grows; an error estimate based on the truncation energy, rather than basis count, would be a concrete way to test and improve that choice."],"forward_implications":["Energies, spatial and momentum densities, and correlation functions at any repulsion $g$ follow from one basis construction and one diagonalization, so the cost of scanning interaction strengths is negligible once the basis exists.","The approach is exact at $g\\to 0$ and $g\\to+\\infty$, so the strong-coupling regime—the hardest for perturbative and many numerical methods—is where the interpolation is most reliable.","The same machinery applies to different trap geometries, demonstrated here for a harmonic well and a smooth double well, as long as the infinite-coupling states for that geometry are available.","Convergence is very fast for small particle numbers and improves with basis size, but it slows as the number of majority fermions grows, as seen in the 6+1 case."],"supporting_citations":[{"why":"It supplies the exact infinite-repulsion eigenstates used as the second half of the variational basis.","marker":"[37, 38]"},{"why":"It introduces the two-state interpolatory ansatz that this paper generalizes.","marker":"[45]"},{"why":"It gives the exact two-particle harmonic-trap solution used as a benchmark.","marker":"[46]"},{"why":"It supplies the matrix-product-state comparison data for the larger systems.","marker":"[47]"}],"fun_headline_variants":["Interpolating basis spans polaron extremes exactly","Truncated basis from two limits solves 1D polaron","Zero-plus-infinite basis captures polaron at all couplings","One ansatz, two limits: exact polaron energies","Exact endpoints, accurate middle: variational polaron basis"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes that correct infinite-coupling eigenstates—specifically the coefficients $a_\\mu$—are known for the trap being treated; for the double-well calculations the paper does not state how these coefficients were obtained, so those results inherit an unspecified external computation.","fun_headline_variants_meta":{"raw":{"variants":["Interpolating basis spans polaron extremes exactly","Truncated basis from two limits solves 1D polaron","Zero-plus-infinite basis captures polaron at all couplings","One ansatz, two limits: exact polaron energies","Exact endpoints, accurate middle: variational polaron basis"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000649,"raw_usage":{"total_tokens":2911,"prompt_tokens":809,"completion_tokens":2102,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":425,"completion_tokens_details":{"reasoning_tokens":2021}},"tokens_in":425,"tokens_out":2102,"duration_ms":16701,"temperature":1.0,"reasoning_tokens":2021,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:32:54.750726+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Independently compute the double-well $g=1$ and $g=10$ ground-state energies and densities with a method that does not use the paper's infinite-coupling states—for example, a lattice or tensor-network simulation with extrapolation—and compare; deviations larger than the matrix-product-state discrepancies shown in the paper would falsify the claimed accuracy for that geometry. A more direct check is to solve the strong-coupling eigenstate conditions (Eq. (7)) for the double-well trap numerically and compare the resulting coefficients $a_\\mu$ with the ones implicit in the paper's basis: a mismatch would show the reported results rest on an unverified input.","supporting_citations":[{"cited_title":"An interpolatory ansatz captures the physics of one- dimensional confined fermi systems,","cited_arxiv_id":null,"evidence_quote":"It introduces the two-state interpolatory ansatz that this paper generalizes."},{"cited_title":"Two cold atoms in a harmonic trap,","cited_arxiv_id":null,"evidence_quote":"It gives the exact two-particle harmonic-trap solution used as a benchmark."},{"cited_title":"Out-of-equilibrium dynamics with matrix product states,","cited_arxiv_id":null,"evidence_quote":"It supplies the matrix-product-state comparison data for the larger systems."}],"review_version":1}