{"id":"5ef1c782-259c-44aa-8019-dd0f7202ba5b","arxiv_id":"2608.06883","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Van Hove singularities in rf spectra of the 3D attractive Fermi Hubbard model provide kinematics-based markers from which the pairing gap and chemical potential can be extracted without full spectral fitting.","lead":"This paper shows that sharp features in the radio-frequency spectrum of fermions in a 3D optical lattice, known as van Hove singularities, can be used to read off the pairing gap and chemical potential from simple algebraic formulas. The approach works even when only partial momentum resolution is available, which matters for current quantum simulation experiments.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Broadening-induced shifts of vHS peak positions are not quantified; the paper never extracts Δ and μ from the broadened spectra it computes, so the central two-marker inversion is unvalidated.","rationale":"The paper's kinematic derivation from the BCS-like Green's function is internally sound: the delta-function resonance condition (Eq. 4) implies the vHS detunings are exactly given by Eq. (6), and the singularity locations in momentum are fixed by ∇ε_k=0. The paper deserves credit for recognizing this structural feature and for testing visibility under phenomenological broadening. However, the numerical validation is incomplete in a way that directly bears on the central claim. The broadened spectral function (Eq. 9) replaces delta functions with Lorentzians, so the vHS features in I(ν) are rounded maxima rather than singularities. The position of a rounded maximum depends on the broadening, on the coherence-factor and Fermi-function prefactors, and on overlap with neighboring features; it is not automatically equal to the kinematic detuning ν_i. The paper asserts that the features remain 'centered close' to ν_i, but this assertion is never quantified, and the illustrative extraction in Sec. III.E plugs the known self-consistent values into Eq. (10) rather than reading positions off the broadened spectra. My concern is therefore that the proposed two-marker protocol may return biased (Δ, μ) even if the BCS-like ansatz is exact, because the measured feature positions are shifted. This is the most load-bearing soft spot for the practical claim, and it is testable within the paper's own framework by applying the inversion to the computed broadened spectra and reporting extraction errors. The reader's chosen weakest assumption (validity of the BCS-like spectral form) is related but distinct; it concerns external validity rather than internal consistency. I thus partially agree with the reader and keep the CONDITIONAL verdict, since the proposed test could either confirm or refute the central claim.","tokens_in":13680,"tokens_out":11139,"duration_ms":123147,"concrete_test":"For each (n, U) and broadening pair (Γ0, Γ1) used in Figs. 4 and 7, compute I(ν) from Eq. (9) (or the kz-integrated version) on a fine frequency grid with the self-consistent input (Δ, μ). Locate the vHS features using a generic peak-finding algorithm (e.g., local maxima after smoothing, or inflection points) for the recommended markers, e.g., M1 and M2 at k∥=(0,π). Insert the observed positions into the inversion equations (Eq. 10) to obtain (Δ_ext, μ_ext). Compare with the input (Δ, μ) and report the relative error as a function of broadening and coupling. If the error exceeds ~10–20% for the parameter range where the features are 'well resolved', the central claim of accurate two-marker extraction fails. This test settles the quantitative validity of the protocol within the paper's own model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the vHS features observed in the broadened rf spectra can be identified with the kinematic detunings ν_i of Eq. (6). However, the paper's numerical validation never actually performs the proposed extraction: Sec. III.E inverts Eq. (6) using the known self-consistent (Δ, μ), not positions read off from the broadened spectra of Figs. 4–7. With the phenomenological self-energy of Eq. (8), the spectral function (Eq. 9) is broadened, so the vHS are not singularities but rounded maxima on a smooth background. The observed maximum of a rounded, asymmetric feature (log divergence for saddles M1/M2, square-root onset for extrema M0/M3) can be shifted from ν_i by an amount of order the broadening Γ0, Γ1 = 0.05 (in units of 6t), which is comparable to the detuning separations at weak coupling and to Δ itself in the examples shown. Overlap between adjacent features (e.g., M0 and M1 merging in Fig. 4(c,d)) further biases peak positions. The paper asserts (Sec. III.B) that features 'remain centered close to the kinematic detunings' and are 'weakly sensitive' to broadening, but no bias or error bars are reported. If the shift between the measured feature position and ν_i is a significant fraction of Δ or of the inter-marker spacing, then the two-marker inversion (Eq. 10) will return biased (Δ, μ), undermining the claim that two well-resolved singularities determine both parameters. This is internal to the paper's own formalism: it does not depend on the external validity of the BCS-like ansatz, and it can be tested by applying the proposed extraction to the broadened spectra.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a spectroscopic diagnostic for the 3D attractive Fermi Hubbard model: the rf detunings of four van Hove singularities are claimed to depend only on the pairing gap Δ and the effective chemical potential μ through the algebraic relations in Eq. (6), ν_i = ε_i − μ + sqrt((ε_i − μ)^2 + Δ^2), with ε_i = 0, 2/3, 4/3, 2 in units of 6t. The authors solve the G0G T-matrix pair-fluctuation equations self-consistently for μ, Δ_sc, and Δ_pg, then compute momentum-integrated and k_z-integrated rf spectra with phenomenological broadening parameters Γ0 and Γ1. They argue that two well-resolved vHS features determine Δ and μ algebraically, and that at half filling a single feature suffices because particle-hole symmetry fixes μ = 1. The numerical section shows broadened spectra with vHS features at positions consistent with Eq. (6) and argues that the features remain identifiable under realistic broadening. The paper explicitly acknowledges that the results are obtained within a particular approximation and that other many-body schemes may give different quantitative relations.","tokens_in":13980,"tokens_out":6393,"duration_ms":64448,"significance":"If the central claim is established, this is a useful and experimentally practical diagnostic: vHS-based markers survive partial momentum integration, become more prominent at strong coupling, and avoid full spectral fitting. The derivation of Eq. (6) is transparent and the algebraic inversion is simple. The paper also gives credit for being self-consistent: it solves the G0G equations numerically, identifies concrete observable EDCs at k_∥ = (0, π) as the most favorable, and candidly states the approximation-dependent limitations in the conclusion. However, the numerical support is substantially internal to the same BCS-like/G0G framework used to derive the kinematic relations, and the central extraction step is not actually performed on the broadened spectra. This currently limits the strength of the 'robust route' claim and requires additional numerical demonstration before the protocol can be regarded as validated.","major_comments":[{"comment":"The proposed two-marker inversion is never actually applied to positions read from broadened spectra. The text states that 'using the parameters of Fig. 4(c) ... the M1 and M2 detunings extracted from the spectrum' yield μ ≈ 1.0 and Δ ≈ 0.5, but no algorithm (peak detection, local maximum, or fit) is described, and the values coincide with the kinematic detunings obtained from the known self-consistent (μ, Δ). For the central claim to be supported, the authors must demonstrate the full pipeline: compute the broadened I(ν) and I(k, ν), identify feature positions with a specified operational definition, invert Eq. (10), and compare the resulting (Δ, μ) with the self-consistent input values for representative (U, n) and for several Γ0, Γ1 values, reporting the bias and any associated uncertainty.","section":"Sec. III.E, Eq. (10)"},{"comment":"Broadening-induced shifts of the vHS feature positions are asserted but not quantified. The paper claims that features 'remain centered close to the kinematic detunings' and are 'only weakly sensitive' to Γ0 and Γ1, but no numerical measure of the shift is given. Because the un-broadened features are logarithmic divergences (M1, M2) or square-root onsets (M0, M3), the maxima of the broadened asymmetric features can shift by an amount of order Γ. With Δ = 0.207 at U_c, n = 0.1 (Fig. 1) and Γ0 = Γ1 = 0.05, such a shift is a substantial fraction of Δ and of the inter-marker spacing; for overlapping features such as M0 and M1 in Fig. 4(c,d) the bias is likely larger. This directly affects the accuracy of the inversion, so a quantitative analysis is needed.","section":"Sec. III.B and Sec. II.D"},{"comment":"The validation is circular with respect to the spectral-function ansatz. The broadened spectra are generated from Eq. (9), which is based on the same BCS-like two-branch quasiparticle form used to derive Eq. (6); the numerical results therefore cannot test whether a more general spectral function, such as one with a frequency-dependent pair-fluctuation self-energy or one obtained from the GG or G0G0 schemes, would place the vHS features at the same detunings. The conclusion explicitly acknowledges this limitation. The paper should state more prominently, including in the abstract, that the protocol is established within the G0G/BCS-like approximation, and should either provide independent benchmarks (e.g., spectra from a different approximation or from quantum Monte Carlo data where available) or clearly frame the results as a consistency check rather than a validation of the kinematic relations.","section":"Sec. IV and Eq. (9)"},{"comment":"The phenomenological self-energy in Eq. (8) has a nonzero real part, but its effect on the vHS detuning positions is not analyzed. The kinematic relations in Eq. (6) are derived from the resonance condition ξ_k − ν + E_k = 0 with E_k = sqrt(ξ_k^2 + Δ^2), which ignores the real part of Σ^R_pg. In the full spectral function of Eq. (9), the quasiparticle dispersion is modified by ReΣ^R_pg, so the actual peak locations in Figs. 4–7 may differ from Eq. (6) not only through broadening but also through a real energy shift. The paper should either show that ReΣ^R_pg produces a negligible shift at the vHS momenta, or include this shift in the inversion formula.","section":"Eq. (8) and Eq. (9)"}],"minor_comments":[{"comment":"The vertical lines in Fig. 4 are labeled collectively as 'M0–M3 vHSs'; since some of these detunings nearly coincide at n = 1, please distinguish the individual markers with different line styles or a legend.","section":"Fig. 4 caption"},{"comment":"The k_z-integrated expression contains an unexplained factor [1 − f(...)]^3 and a prefactor 1/π^3; please provide the full derivation or correct the expression, since Eq. (7) is used to generate the k_z-integrated spectra in Figs. 5–7.","section":"Eq. (7)"},{"comment":"The statement that Γ0 and Γ1 'become negligible at low temperatures due to the formation of stable pairing' is not supported by any temperature-dependent calculation; please clarify whether all presented results use fixed Γ0 = Γ1 = 0.05 and how the parameters are expected to scale with temperature.","section":"Sec. II.D"},{"comment":"After broadening, the features are no longer true van Hove singularities but rounded maxima; consider using 'van Hove features' or 'vHS remnants' when referring to the broadened spectra to avoid overstating the singularity.","section":"Sec. III (terminology)"},{"comment":"The Supplemental Material citation lacks an arXiv identifier or DOI, which makes it difficult for readers to locate the self-consistent equations and the additional numerical results for n = 0.2, 0.4, 0.6, and 0.8.","section":"Reference [48]"}],"recommendation":"major_revision","confidential_remarks":"The kinematic relations in Eq. (6) are elegant and potentially useful, and the paper is within the scope of cond-mat.quant-gas. The main reason for major revision is that the central extraction protocol is not actually demonstrated on the broadened spectra: the inversion in Sec. III.E is performed with the known self-consistent (Δ, μ) rather than with positions read off from the computed spectra, and the broadening-induced bias is not quantified. This is fixable within the manuscript's scope by adding the missing numerical analysis and by more carefully framing the approximate nature of the spectral-function ansatz. I would not reject the paper, but the current form overstates the robustness of the proposed diagnostic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know about this paper is that the central idea is new and useful: identify the four van Hove singularities in rf spectra of the 3D attractive Fermi Hubbard model and use their detuning positions as algebraic markers for the pairing gap Δ and chemical potential μ. Equation (6) is a clean consequence of the BCS-like spectral ansatz, and the Figure 8 map is a practical resource. The authors also make a genuinely helpful move by analyzing k_z-integrated spectra and recommending the (0,π) EDC as the most favorable experimental observable. They are honest about the approximation dependence of their scheme in the conclusion, flagging that GG and G0G0 give different gaps and that a fully self-consistent treatment would soften the features.\n\nThe soft spots are real but not fatal. The validation is substantially circular: the numerical spectra are produced by the same G0G pair-fluctuation theory that generates the Δ and μ used in Eq. (6), and there is no independent benchmark against QMC or experiment. More concretely, the stress-test concern lands: the paper never actually performs its proposed extraction on the broadened spectra it computes. Section III.E inverts Eq. (6) using the known self-consistent parameters, not positions read off from Figs. 4–7. The claim that features remain \"centered close\" to the kinematic detunings under broadening is asserted without error bars or bias estimates. With Γ0 = Γ1 = 0.05 in units of 6t, the broadening is comparable to the inter-marker spacing at weak coupling and to Δ itself; overlapping features (as in Fig. 4(c,d)) can shift apparent maxima. That is a testable concern internal to the paper's own formalism, and it should be addressed by applying the inversion to the broadened spectra and reporting the resulting bias.\n\nMinor issues: the paper leans heavily on the authors' prior work, but that is not a flaw given the formalism is theirs and the vHS relations are new. No code or data is provided, which is a modest reproducibility gap.\n\nWho is this for? People working on rf spectroscopy of lattice Fermi gases, especially quantum simulation groups looking for practical diagnostics. The proposal deserves a serious referee; it is not a desk reject. The authors should be required to demonstrate the two-marker inversion on their own broadened spectra and, if possible, compare against a QMC spectral function or an experiment. If that test passes, the method becomes genuinely practical.\n\nRecommendation: send to peer review, with the request for a direct extraction test on broadened spectra.","headline":"A genuinely new set of vHS-based markers for gap extraction in the lattice Fermi Hubbard model, but the paper never demonstrates the actual extraction on its own broadened spectra; worth refereeing, needs a real inversion test.","tokens_in":14589,"tokens_out":1477,"would_cite":true,"duration_ms":16384,"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":"Van Hove singularities in rf spectra pin down the pairing gap and chemical potential of the 3D Fermi Hubbard model","keywords":["Fermi Hubbard model","pairing gap","van Hove singularity","rf spectroscopy","pseudogap","BCS-BEC crossover","optical lattice","quantum simulation"],"falsifier":"Measure the rf spectrum of a 3D attractive Fermi Hubbard gas at known density and interaction strength, resolve at least two van Hove features, and solve the two corresponding equations of Eq. (6) for $\\Delta$ and $\\mu$. If the resulting pair disagrees between different choices of vHS markers, or fails to match an independent measurement of $\\mu$ (for example from band mapping), the single-gap assumption is falsified.","tokens_in":13427,"feed_emoji":"⚛️","tokens_out":9882,"duration_ms":86944,"temperature":0.7,"pith_summary":"The paper sets out to make van Hove singularities in radio-frequency spectra a practical diagnostic for pairing in the three-dimensional attractive Fermi Hubbard model. It claims that four types of singular features appear at rf detunings that depend only on the pairing gap $\\Delta$ and the chemical potential $\\mu$, through $\\nu_i = (\\epsilon_i-\\mu) + \\sqrt{(\\epsilon_i-\\mu)^2+\\Delta^2}$, where the $\\epsilon_i$ are the four lattice band critical-point energies. If this holds, resolving any two of these features gives both $\\Delta$ and $\\mu$ without full spectral fitting, and at half filling one feature suffices. The authors support the claim with numerical spectra from a pair-fluctuation theory, and they show the features remain visible when lifetime and scattering broadenings are added, becoming clearer at strong coupling where back-bending diagnostics lose sensitivity. This matters because quantum-simulated 3D Fermi Hubbard systems currently lack a direct spectroscopic probe of pairing in exactly the regime where theory predictions differ most.","feed_headline":"Two rf kinks reveal the pairing gap in 3D Fermi Hubbard gases","feed_subtitle":"Four lattice singularities fix gap and chemical potential algebraically, no full spectral fitting needed.","key_machinery":"The load-bearing object is the van Hove singularity of the cubic-lattice dispersion, a momentum point where $\\nabla_{\\mathbf{k}}\\xi_{\\mathbf{k}}=0$ and the density of states is non-analytic. The lattice has four such points (labelled $M_0,\\dots,M_3$, realised at $\\Gamma$, $X$, $M$, $R$), and the rf resonance condition $\\xi_{\\mathbf{k}}-\\nu+E_{\\mathbf{k}}=0$ maps each one to a sharp feature in the rf current. The argument is carried by the kinematic relation $\\nu_i=(\\epsilon_i-\\mu)+\\sqrt{(\\epsilon_i-\\mu)^2+\\Delta^2}$, which converts measured detunings into the two pairing parameters.","core_discovery":"Within the $G_0G$ T-matrix pair-fluctuation scheme, the spectral function of the attractive Fermi Hubbard model is taken to have the BCS two-branch form $A(\\mathbf{k},\\omega)=2\\pi[u_{\\mathbf{k}}^2\\delta(\\omega-E_{\\mathbf{k}})+v_{\\mathbf{k}}^2\\delta(\\omega+E_{\\mathbf{k}})]$ with $E_{\\mathbf{k}}=\\sqrt{\\xi_{\\mathbf{k}}^2+\\Delta^2}$. The paper shows that when this is fed into the rf response $I(\\nu)=\\frac{1}{2\\pi}\\sum_{\\mathbf{k}} A(\\mathbf{k},\\omega)f(\\omega)|_{\\omega=\\xi_{\\mathbf{k}}-\\nu}$, the lattice critical points—$\\Gamma(0,0,0)$, $X(\\pi,0,0)$, $M(\\pi,\\pi,0)$, and $R(\\pi,\\pi,\\pi)$—produce four singular features whose detunings obey $\\nu_i=(\\epsilon_i-\\mu)+\\sqrt{(\\epsilon_i-\\mu)^2+\\Delta^2}$ with $\\epsilon_i=0,\\,2/3,\\,4/3,\\,2$ in units of $6t$. Inverting any two of these relations yields $\\Delta$ and $\\mu$ algebraically, and numerical spectra with phenomenological broadenings confirm that the features stay near these kinematic positions. At half filling, particle-hole symmetry fixes $\\mu=1$, reducing the extraction to a single vHS measurement.","pith_inferences":["If a fully self-consistent treatment introduces frequency-dependent pair-fluctuation self-energy, the vHS positions may shift or soften; this can be tested by computing spectra with finite pair momentum included and comparing extracted $(\\Delta,\\mu)$ with the kinematic map.","The same reasoning should transfer to other lattice geometries whenever the band has critical points and pairing is nearly momentum-independent, though the band-energy offsets $\\epsilon_i$ would change.","A consistency check within one dataset—extracting $(\\Delta,\\mu)$ from two different vHS pairs—would reveal whether a single-gap description actually holds; disagreement would signal momentum-dependent pairing beyond the BCS-like ansatz."],"forward_implications":["Two resolved vHS features determine both $\\Delta$ and $\\mu$ from a single dataset, with no dispersion fitting.","The scheme remains usable in $k_z$-integrated spectra, where the strongest signal appears at $\\mathbf{k}_\\parallel=(0,\\pi)$.","At half filling, one well-resolved vHS fixes $\\Delta$ because particle-hole symmetry pins $\\mu=1$.","The features sharpen relative to the background at stronger coupling, exactly where back-bending loses sensitivity.","Because the detuning–gap map is kinematic, it can serve as a cross-method benchmark for comparing gap values from different many-body theories."],"supporting_citations":[{"why":"Supplies the rf current formula linking the spectral function to detuning, the starting point for the vHS analysis.","marker":"[39]"},{"why":"Provides the equation-of-motion derivation underlying the G0G T-matrix scheme used for the self-consistent solutions.","marker":"[22]"},{"why":"Gives the pair-fluctuation theory and T-matrix formalism used to compute $\\Delta$, $\\mu$, and the rf spectra.","marker":"[45]"},{"why":"Introduces the phenomenological lifetime and scattering broadenings ($\\Gamma_0$, $\\Gamma_1$) used to test vHS visibility.","marker":"[46]"},{"why":"Reports experimental rf spectra in strongly interacting Fermi gases that support the BCS-like two-branch spectral function.","marker":"[43]"},{"why":"Demonstrates momentum-resolved photoemission/rf spectroscopy in a 2D Fermi-Hubbard system, motivating the partial-momentum-integration analysis.","marker":"[31]"},{"why":"Provides the detailed self-consistent equations and additional numerical results for intermediate fillings.","marker":"[48]"}],"fun_headline_variants":["Two rf kinks yield pairing gap and chemical potential algebraically","No spectral fitting: rf van Hove singularities give pairing gap","Four lattice vHS solve gap and mu with two rf measurements","At half filling, one rf kink extracts pairing gap in Hubbard gases"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The extraction rests on the spectral function keeping its BCS-like two-branch form with one momentum-independent gap $\\Delta$; if pairing fluctuations distort that form, the kinematic relations for the vHS detunings no longer hold.","fun_headline_variants_meta":{"raw":{"variants":["Two rf kinks yield pairing gap and chemical potential algebraically","No spectral fitting: rf van Hove singularities give pairing gap","Four lattice vHS solve gap and mu with two rf measurements","At half filling, one rf kink extracts pairing gap in Hubbard gases"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000693,"raw_usage":{"total_tokens":3161,"prompt_tokens":998,"completion_tokens":2163,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":614,"completion_tokens_details":{"reasoning_tokens":2089}},"tokens_in":614,"tokens_out":2163,"duration_ms":17132,"temperature":1.0,"reasoning_tokens":2089,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:00:36.164965+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the rf spectrum of a 3D attractive Fermi Hubbard gas at known density and interaction strength, resolve at least two van Hove features, and solve the two corresponding equations of Eq. (6) for $\\Delta$ and $\\mu$. If the resulting pair disagrees between different choices of vHS markers, or fails to match an independent measurement of $\\mu$ (for example from band mapping), the single-gap assumption is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the rf current formula linking the spectral function to detuning, the starting point for the vHS analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the pair-fluctuation theory and T-matrix formalism used to compute $\\Delta$, $\\mu$, and the rf spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the phenomenological lifetime and scattering broadenings ($\\Gamma_0$, $\\Gamma_1$) used to test vHS visibility."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports experimental rf spectra in strongly interacting Fermi gases that support the BCS-like two-branch spectral function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates momentum-resolved photoemission/rf spectroscopy in a 2D Fermi-Hubbard system, motivating the partial-momentum-integration analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the detailed self-consistent equations and additional numerical results for intermediate fillings."}],"review_version":1}