{"id":"4125d13a-2075-4c48-b76f-49f75455b58f","arxiv_id":"2411.17895","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Exact Bethe-ansatz calculation of the 1D Hubbard Fermi polaron spectral function reveals anomalous Fermi singularities at the Brillouin-zone boundary and collective polaron quasiparticle peaks near quarter filling.","lead":"Using the Bethe ansatz, the authors compute the impurity spectral function of a single mobile fermion in a one-dimensional lattice Fermi bath, finding that at high momentum the usual Fermi edge singularities become one-sided 'anomalous' singularities, and near quarter filling broad polaron quasiparticle peaks appear. The results provide an exact benchmark for approximate polaron theories and make concrete predictions for cold-atom experiments in optical lattices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The polaron-quasiparticle identification rests on an uncontrolled extrapolation of the finite-size exponent alpha_II to zero; the data show alpha_II = 0.31 and 0.17 for two broadening choices, so the coexistence claim is not yet established.","rationale":"The reader's conditional verdict identifies the truncation of the Bethe-state sum and the finite-size scaling of the polaron peak as the weakest assumptions. I agree with that assessment and sharpen it: the single most load-bearing assumption is that alpha_II vanishes in the thermodynamic limit. The paper's headline result is the coexistence of two distinct spectral features, and the polaron peak is the only feature whose classification is not backed by an independent analytical prediction or a robust numerical trend. The reported exponents (0.31 and 0.17) are clearly nonzero at the accessible system sizes, and the change with broadening choice shows that the extrapolation is sensitive to the fitting procedure. The sum rule, while reassuring for global completeness, does not distinguish a broad quasiparticle from a dense set of small-residue states forming a singularity; both can have small total weight. The integrated weight test is the decisive check because a quasiparticle carries a finite spectral weight in the thermodynamic limit, whereas a power-law singularity's weight vanishes as delta^{1-alpha}. The paper makes a strong technical contribution — the exact Bethe form factors, including irregular states, and the successful reproduction of the Q = 0 exponent are real achievements — so a conditional verdict is appropriate pending this quantitative test. No change to the verdict is needed.","tokens_in":13711,"tokens_out":4514,"duration_ms":42792,"concrete_test":"Compute the integrated spectral weight W_II(L) = integral of A(pi, omega) over a window around peak II that excludes the singularity I, for L = 60, 80, 100, 120, 160 with delta = t/L and delta = t/(2L). For a genuine polaron quasiparticle, W_II(L) should approach a nonzero constant as L -> infinity; for a power-law singularity with exponent alpha, W_II(L) ~ delta^{1-alpha} -> 0. If W_II(L) decreases systematically with L, the 'polaron' is actually a Fermi singularity. Additionally, refit alpha_II including a subleading correction A_max = c1 delta^{-alpha} + c2 and report confidence intervals; if the best-fit alpha_II is consistent with zero within errors, the quasiparticle interpretation is supported. This test directly settles whether the central coexistence claim survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most novel claim is the coexistence of anomalous Fermi singularities with polaron quasiparticles near quarter filling at large momentum. The singularity classification is supported by finite-size exponents (alpha_I ~ 0.89, alpha_III ~ 0.50), but the polaron classification of peak II rests entirely on the finite-size scaling of its maximum value A_max. In Appendix B, the authors fit A_max ~ delta^{-alpha_II} and obtain alpha_II = 0.31 using delta = 4t/L and alpha_II = 0.17 using delta = t/L with L = 60, 80, 100. They then write that it is 'reasonable to believe' that alpha_II -> 0 with more numerical effort, but this is an extrapolation from three points with no error bars and no model for the overlap with the nearby singularity I. A nonzero alpha_II would mean that peak II is itself a power-law singularity, not a quasiparticle, invalidating the central 'coexistence' claim. The large sum-rule saturation (rho_s > 99.8%) does not resolve this: a singularity with small total weight can still produce a diverging peak height. The distinction between a quasiparticle and a weak singularity is quantitative and must be tested by a quantity that distinguishes them, such as the integrated spectral weight under the peak.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Using the Bethe-ansatz solution of the one-dimensional Hubbard model, the authors derive explicit form factors for the overlap between the non-interacting impurity-plus-Fermi-sea product state and the regular Bethe states, as well as the irregular spin-flip and eta-pairing states. They then evaluate the impurity spectral function A(Q,omega) on finite lattices by summing over states with up to three pseudo particle-hole pairs, and analyze the resulting spectra for various momenta and fillings. The main findings are: (i) at small momentum the spectrum shows two power-law Fermi-edge singularities, with the Q=0 exponent matching the known Castella-Zotos value; (ii) at large momentum the singularities become two-sided and eventually, at Q=pi, develop into 'anomalous' Fermi singularities with low-energy-oriented power-law tails; and (iii) near quarter filling, two broad peaks at large momentum are interpreted as polaron quasiparticles collectively generated by many states. The classification of the singularities and the purported polaron peaks rests on finite-size scaling of peak heights A_max ~ delta^{-alpha}, with exponents alpha=0.89 and 0.50 for the singularities and alpha=0.31 (for delta=4t/L) or 0.17 (for delta=t/L) for the claimed polaron peak.","tokens_in":13863,"tokens_out":8049,"duration_ms":70669,"significance":"If established, this work provides the first exact Bethe-ansatz spectral function for lattice Fermi polarons including the contributions of irregular spin-flip and eta-pairing states, and it uncovers a qualitatively new spectral structure at large momentum that is absent in higher dimensions. The paper has clear strengths: the form-factor derivation is nontrivial, the sum rule exceeds 99.8%, the Q=0 singularity exponent (alpha_fit=0.871) agrees with the independent analytical value (0.875), and the Chevy-ansatz comparison gives a useful benchmark. However, the central novelty—the coexistence of anomalous Fermi singularities with polaron quasiparticles—is only as convincing as the evidence that the second peak is not itself a power-law singularity. The finite-size scaling data show a nonzero exponent for both broadening choices, and the extrapolation to zero is not controlled. Thus the significance would be greatly enhanced by a more rigorous identification of the polaron peak, either through a controlled scaling that isolates the overlap with the adjacent singularity or through a quantity that directly distinguishes a true quasiparticle from a power-law singularity.","major_comments":[{"comment":"The classification of peak II as a polaron quasiparticle, rather than as a power-law singularity, is load-bearing for the central claim of coexistence. The finite-size scaling in Figs. 6(b) and 7 yields alpha_II = 0.31 for delta = 4t/L and alpha_II = 0.17 for delta = t/L, both clearly nonzero. The sentence 'It is reasonable to believe that, by taking the limit delta -> 0 with more numerical efforts, we might eventually confirm alpha_II = 0' is an unsupported extrapolation from three data points without error bars and without a model for the overlap with singularity I that the authors themselves invoke to explain the L-dependence. A nonzero alpha_II would imply that peak II is itself a power-law singularity with vanishing integrated weight in the thermodynamic limit, which would invalidate the coexistence claim. The authors should either provide a controlled scaling that explicitly removes the overlap contribution and shows alpha_II -> 0, or measure a quantity that distinguishes a quasiparticle from a singularity, such as the integrated spectral weight under peak II as a function of L (a delta-function peak has finite weight, whereas a power law with 0<alpha<1 has weight scaling as delta^{1-alpha} -> 0).","section":"Appendix B, Figs. 6-7"},{"comment":"The paper repeatedly calls the result 'exact', but the sum in Eq. (4) is truncated to three pseudo particle-hole pairs. The stated sum rule rho_s > 99.8% bounds the total omitted residue, but it does not by itself guarantee that the omitted states contribute negligibly to a specific feature such as peak II, whose total weight may be only a few percent of the total spectral weight. A convergence check (e.g., comparing the spectral function and the peak height with truncations at two and three pairs, or reporting the integrated weight under peak II as a function of L) is needed to support the statement that 'the omitted states have negligible effect on the convergence of our numerics' and to secure the polaron classification.","section":"Spectral function section (page 3)"}],"minor_comments":[{"comment":"The phrase 'the leftest hole state' should be corrected to 'the leftmost hole state'.","section":"Page 3, near Fig. 3"},{"comment":"The phrase 'eight-spoked asterisks' is awkward; consider 'eight-pointed asterisks' or simply 'asterisks'.","section":"Caption of Fig. 2"},{"comment":"The broadening factor delta = 4t/L is introduced without justification; a brief comment on the sensitivity of the results to this choice (partly addressed in Appendix B) would be helpful for reproducibility.","section":"Eq. (4) and spectral function section"},{"comment":"The term 'exact spectral function' is used despite the truncation to three pseudo particle-hole pairs; consider using 'numerically exact' or 'highly accurate' to align the wording with the stated limitation in the spectral function section.","section":"Abstract and Introduction"}],"recommendation":"major_revision","confidential_remarks":"The main technical contribution—the form-factor derivation for regular and irregular Bethe states—appears sound and will be useful to the community. However, the central claim of coexistence depends on the polaron-peak classification, which is under-supported by the scaling analysis. I would recommend asking the authors to either strengthen this evidence with a more controlled finite-size analysis or to soften the claim. The paper does not include a data/code availability statement; requesting one would be reasonable given the computational nature of the results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is worth a serious look: it computes the impurity spectral function of the 1D Hubbard model from Bethe-ansatz form factors for all regular excited states plus the irregular spin-flip and eta-pairing states. That's a genuine extension of Castella-Zotos to the lattice at finite U, and it's the first time these irregular states are included. The Q=0 check is the strongest evidence: the numerical exponent 0.871 tracks the analytical 0.875 to three digits. That tells me the form-factor machinery is sound.\n\nThe new physics is in the large-momentum features: the two-sided Fermi singularities that turn into 'anomalous' singularities at Q=pi, with power-law tails to low energy. These are well supported by the residue distributions and by exponents alpha_I=0.89 and alpha_III=0.50, which behave like singularities.\n\nThe soft spot is the second peak at Q=pi, claimed to be a polaron quasiparticle. The argument there is a finite-size scaling of its maximum height: alpha_II=0.31 with delta=4t/L, dropping to 0.17 with delta=t/L. That is not zero, and the extrapolation to alpha_II=0 is a hope, not a result. The authors attribute the residual exponent to overlap with the neighboring singularity, but they don't show that separation works, and they don't report error bars on the fits. A weak singularity with small total weight could look exactly like this. The proper diagnostic is the integrated spectral weight under the peak, which should saturate to a constant for a quasiparticle and shrink for a singularity; they don't show it.\n\nThe truncation to three particle-hole pairs is well covered by the sum rule (rho_s>99.8%), and the omitted weight is unlikely to change the qualitative picture. Still, 'exact' in the title is slightly stronger than what is computed. No code or data is released, which doesn't help.\n\nOverall: the anomalous Fermi singularities are on solid ground; the polaron coexistence is plausible but not established. The paper deserves peer review, and the referee should push on the alpha_II question. I'd read it for the form-factor technique alone.","headline":"Solid exact Bethe-ansatz spectral function for the 1D Hubbard polaron with a convincing Q=0 benchmark, but the central 'coexistence' claim rests on a finite-size exponent that doesn't yet extrapolate to zero.","tokens_in":14516,"tokens_out":3201,"would_cite":true,"duration_ms":27199,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.38.-k","71.10.Fd","03.75.Ss","05.30.Fk"],"model":"deepseek-v4-flash","headline":"A single mobile impurity in a one-dimensional lattice Fermi gas, at large momentum, develops anomalous Fermi singularities with low-energy power-law tails and broad polaron quasiparticle peaks that are collectively generated by many…","keywords":["Fermi polaron","one-dimensional Hubbard model","Bethe ansatz","spectral function","Fermi edge singularity","polaron quasiparticle","form factor","optical lattice"],"falsifier":"An independent, untruncated calculation of A(Q=π,ω) at quarter filling with U=4t — for instance, by time-dependent density-matrix renormalization group or by including four or more particle-hole pairs in the Bethe sum — should reproduce the anomalous low-energy tail and a polaron peak whose height stays finite as the broadening δ→0. Alternatively, a cold-atom Ramsey interferometry measurement of an impurity accelerated to momentum Q=π in a one-dimensional optical lattice at quarter filling should show the predicted power-law tail and the collective peak; if the tail vanishes or the peak height scales like a singularity exponent, the central claim fails.","tokens_in":13331,"feed_emoji":"⚛️","tokens_out":8694,"duration_ms":68592,"temperature":0.7,"pith_summary":"This paper establishes that the spectral function of a single mobile impurity in a one-dimensional lattice Fermi gas contains two previously unrecognized features at large impurity momentum: anomalous Fermi singularities, whose power-law tails extend to low energy at the Brillouin-zone boundary, and broad polaron quasiparticle peaks that are generated collectively by many excited many-body states rather than by a single dominant state. The results are derived exactly from the Bethe-ansatz solution of the one-dimensional Hubbard model, with explicit form factors computed for regular Bethe states and for the irregular spin-flip and eta-pairing states. Near quarter filling and momentum Q=π, the anomalous singularities and the polaron peaks coexist in the spectrum, a situation with no analogue in two or three dimensions. If the results are right, they provide benchmark spectra for approximate polaron theories and a concrete prediction for cold-atom Ramsey interferometry in optical lattices.","feed_headline":"Bethe ansatz reveals anomalous Fermi singularities in 1D polarons","feed_subtitle":"Power-law tails and collective quasiparticle peaks at momentum π are testable in cold-atom lattices.","key_machinery":"The central object is the Bethe-ansatz form factor: the overlap of an eigenstate of the one-dimensional Hubbard model with the non-interacting state in which the impurity is a plane wave and the bath fermions form a Fermi sea. The authors write this overlap as an (N+1)×(N+1) determinant, using a Slater-determinant identity, and obtain the norm of the Bethe wavefunction from the same determinant structure. The machinery also includes a full classification of the eigenstates — real-k, k−Λ, spin-flip, and η-pairing — so that the sum over states in the spectral function is nearly complete, and a sum rule ρ_s > 99.8% certifies the truncation at three pseudo particle-hole pairs.","core_discovery":"Using the Lieb-Wu Bethe ansatz for the Hubbard model with one spin-down impurity and N spin-up fermions, the authors derive determinant expressions for the form factor of every eigenstate: the overlap of a Bethe state with the non-interacting product of a plane-wave impurity and a Fermi sea. This includes regular states with real quasi-momenta, k−Λ states with a complex momentum pair, and the irregular spin-flip and η-pairing states. The impurity spectral function A(Q,ω) is then obtained as a sum over roughly half a million states (up to three pseudo particle-hole pairs, with sum rule above 99.8 percent) with a Lorentzian broadening δ=4t/L. At zero momentum the spectrum shows the two conventional power-law Fermi edge singularities; as Q increases each singularity becomes two-sided, and at Q=π the high-energy side disappears, leaving anomalous singularities with low-energy power-law tails. Near quarter filling a broad polaron peak appears at ω≈3t for U=4t, built from many states with residue around $10^{-3}$ (a second bundle near ω≈6t is masked by the upper singularity); finite-size scaling of the peak height distinguishes this collective polaron peak (exponent α≈0.31, decreasing toward 0) from the singularities (α≈0.50 and α≈0.89).","pith_inferences":["If the collective polaron peak is genuine, the conventional notion of quasiparticle weight (the largest single residue) may need to be replaced by an integrated 'bundle weight' over a small energy window; its scaling with system size could define a new diagnostic for quasiparticles in integrable systems.","The same determinant form-factor technique, applied to two impurities, could reveal whether the anomalous singularities persist and how the collective polaron peaks hybridize; the paper's hint at spin-charge separation suggests spinon contributions would appear in the two-impurity spectrum near Q=0 or Q=π.","A concrete experimental test: measure the Ramsey overlap S(t) after rapidly transferring an impurity to momentum Q=π; the power-law exponent of the decay tail should match the anomalous singularity exponent α≈0.89, while the polaron peak should appear as a slowly decaying component at a different frequency."],"forward_implications":["At Q=π and quarter filling, the impurity spectrum simultaneously shows an anomalous Fermi singularity and a polaron quasiparticle peak; both are exact predictions of the Bethe ansatz.","The polaron peak is collective: no single eigenstate carries it, so conventional single-pole quasiparticle language (one dominant residue) does not describe it.","Finite-size scaling distinguishes the two kinds of features: peak heights of singularities diverge as δ^{-α} with α≈0.89 and 0.50, while the polaron peak's exponent is much smaller (≈0.31, decreasing to ≈0.17 with smaller broadening), consistent with a constant height in the thermodynamic limit.","The Chevy variational ansatz with two-particle-hole excitations reproduces the polaron peak, giving a benchmark for approximate many-body theories.","The same spectral structure appears for attractive interactions U=-4t under a particle-hole transformation, indicating the coexistence is a generic feature of the one-dimensional lattice polaron, independent of the sign of interaction."],"supporting_citations":[{"why":"Supplies the exactly solvable one-dimensional Hubbard model (Lieb-Wu) on which the entire calculation is based.","marker":"[35]"},{"why":"Provides the classification of Bethe eigenstates and excitations of the 1D Hubbard model used to enumerate the states in the spectral sum.","marker":"[36]"},{"why":"Prior exact form-factor calculation for a particle in a 1D Fermi gas that established the Fermi edge singularity and whose determinant method is generalized here.","marker":"[38]"},{"why":"Constructs the irregular spin-flip and eta-pairing eigenstates that are required for completeness of the spectral sum.","marker":"[41]"},{"why":"Introduces the determinant representation of the Bethe wavefunction that underlies the form-factor computation.","marker":"[50]"},{"why":"Gives the Slater-determinant overlap identity used to evaluate the form factors and normalization factors.","marker":"[52]"},{"why":"Provides the Chevy variational spectral function with two-particle-hole excitations used as a quantitative benchmark for the exact results.","marker":"[26]"},{"why":"Demonstrates the use of the determinant form factor and the sum rule for spectral functions of impurities in 1D quantum liquids, the technical route adopted here.","marker":"[42]"}],"fun_headline_variants":["Exact 1D polaron spectrum reveals anomalous Fermi singularities","Bethe ansatz exposes exotic singularities in 1D Fermi polarons","1D polarons: exact spectrum shows anomalous edge singularities","Collective polaron peak predicted from exact Bethe ansatz","Anomalous Fermi singularities emerge in exact 1D polaron spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation's load-bearing premise is that the spectral function is converged after summing only up to three pseudo particle-hole pairs (about half a million states) at finite sizes L=20–100, so that the thermodynamic-limit singularities and the polaron peak identified by finite-size scaling with a chosen broadening δ=4t/L are not artifacts of the truncation.","fun_headline_variants_meta":{"raw":{"variants":["Exact 1D polaron spectrum reveals anomalous Fermi singularities","Bethe ansatz exposes exotic singularities in 1D Fermi polarons","1D polarons: exact spectrum shows anomalous edge singularities","Collective polaron peak predicted from exact Bethe ansatz","Anomalous Fermi singularities emerge in exact 1D polaron spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000193,"raw_usage":{"total_tokens":1401,"prompt_tokens":1045,"completion_tokens":356,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":661,"completion_tokens_details":{"reasoning_tokens":263}},"tokens_in":661,"tokens_out":356,"duration_ms":3362,"temperature":1.0,"reasoning_tokens":263,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:43:26.338610+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent, untruncated calculation of A(Q=π,ω) at quarter filling with U=4t — for instance, by time-dependent density-matrix renormalization group or by including four or more particle-hole pairs in the Bethe sum — should reproduce the anomalous low-energy tail and a polaron peak whose height stays finite as the broadening δ→0. Alternatively, a cold-atom Ramsey interferometry measurement of an impurity accelerated to momentum Q=π in a one-dimensional optical lattice at quarter filling should show the predicted power-law tail and the collective peak; if the tail vanishes or the peak height scales like a singularity exponent, the central claim fails.","supporting_citations":[{"cited_title":"Deguchi, F","cited_arxiv_id":null,"evidence_quote":"Provides the classification of Bethe eigenstates and excitations of the 1D Hubbard model used to enumerate the states in the spectral sum."},{"cited_title":"Exact spectral properties of Fermi polarons in one-dimensional lattices: Anomalous Fermi singularities and polaron quasiparticles","cited_arxiv_id":"2411.17895","evidence_quote":"Prior exact form-factor calculation for a particle in a 1D Fermi gas that established the Fermi edge singularity and whose determinant method is generalized here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Constructs the irregular spin-flip and eta-pairing eigenstates that are required for completeness of the spectral sum."},{"cited_title":"Our integer quantum number sj is related to Ij by sj = Ij + 1/2","cited_arxiv_id":null,"evidence_quote":"Introduces the determinant representation of the Bethe wavefunction that underlies the form-factor computation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Slater-determinant overlap identity used to evaluate the form factors and normalization factors."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Chevy variational spectral function with two-particle-hole excitations used as a quantitative benchmark for the exact results."},{"cited_title":"Here, we adopt a single- particle wavefunction χj(y) that is more suitable to han- dle the k − Λ solutions with a pair of complex-valued quasi-momenta [51]","cited_arxiv_id":null,"evidence_quote":"Demonstrates the use of the determinant form factor and the sum rule for spectral functions of impurities in 1D quantum liquids, the technical route adopted here."}],"review_version":1}