{"id":"c4c9bac9-40dd-44ba-9b8d-83678d19ffd5","arxiv_id":"2608.11004","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the interacting SSH model, QMC data show that the low-frequency response of chiral density to an electric field remains equal to the electrical conductivity, so local Hubbard interactions do not renormalize the 1D chiral effect.","lead":"This paper uses quantum Monte Carlo simulations to test whether electron-electron interactions change the chiral effect in a one-dimensional Dirac semimetal. It finds that, for the studied parameters, the relation between induced chiral charge and electric current survives unchanged, even when interactions open a gap.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"SAC reconstruction is the load-bearing step: no error bars or validation show that it preserves the low-frequency ratio, so the observed λ≈σ could be a reconstruction artifact.","rationale":"The reader's conditional verdict is appropriately cautious. The central numerical result is entirely mediated by SAC; without uncertainty quantification or validation on an exactly solvable case, the equality at ω→0 is not established at the claimed confidence. The paper's own Section VI admits that in the gapped regime the equilibrium linear-response anomaly formula no longer applies, so the broad conclusion that interactions do not renormalize the chiral effect is stronger than the evidence presented. A synthetic-data benchmark or a direct τ-space ratio would settle the SAC issue and is cheap given the public code. The reader already identified this same weakest assumption, so no adjustment to the CONDITIONAL verdict is needed.","tokens_in":46,"tokens_out":13900,"duration_ms":206421,"concrete_test":"Benchmark the SAC pipeline on an exactly known case: using the exact noninteracting spectral functions for finite L, generate synthetic imaginary-time correlators with Gaussian noise at the same statistics as the QMC runs, reconstruct λ and σ with the identical SAC algorithm and default model, and report the distribution of the ratio λ(0)/σ(0). If the median deviates from 1 by more than the claimed numerical accuracy, the interacting result is unreliable; if it is unbiased, the SAC concern is resolved. A complementary direct check is to compute R(τ)=⟨ρ5(0)J(τ)⟩/⟨J(0)J(τ)⟩ from the raw QMC data and test whether R(τ)≈1 at large τ for the nontrivial U cases, e.g., Model 2 at U=1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section V obtains the central equality λ(ω)≈σ(ω) by stochastic analytic continuation of two different imaginary-time correlators, Eqs. (13)–(14). SAC is an ill-posed inversion; the paper gives no error bars and no benchmark on an exactly known case. In the free case the exact spectra contain δ(ω) peaks (Fig. 1), so the reconstruction necessarily broadens them; for interacting U, the same smoothing could push independent λ and σ reconstructions toward similar low-frequency shapes even if the true zero-frequency weights differ. The high-frequency parts of λ and σ are clearly different (Figs. 2, 3, 9, 11), yet no test shows that spectral-weight leakage from high frequencies does not bias the ω→0 ratio. Additionally, in the gapped U=4 Model 2 case both reconstructed functions vanish, so the claim there is the trivial 0=0 and provides no constraint on renormalization. The bare operator ρ5=ψ†σ2ψ is also assumed to remain the physical chiral density after interactions; if renormalization mixes it with other local operators, the measured equality may not represent the chiral effect.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates whether the relation j = e v_F rho_5 between the electric current and the axial charge density, previously derived analytically for the noninteracting Su-Schrieffer-Heeger model with finite dissipation [34], survives the addition of local Hubbard interactions. Two variants are studied by finite-temperature determinant Quantum Monte Carlo: spinful electrons with the on-site Hubbard term (Model 1) and spinless electrons with an intracell density-density term (Model 2). The axial-charge response lambda(omega) and the electrical conductivity sigma(omega) are obtained from imaginary-time density-current and current-current correlators [Eqs. (13)-(14)] via stochastic analytic continuation, and the single-particle spectral function A(omega,k) is reconstructed to monitor the interaction-induced gap. The authors report that lambda(omega) is approximately equal to sigma(omega) as omega tends to zero for all studied interactions, including Model 2 in the gapped U=4 case, and conclude that the chiral effect is not renormalized by local density-density interactions, suggesting a possible topological origin of the relation.","tokens_in":28,"tokens_out":17375,"duration_ms":221947,"significance":"If the equality lambda(0)=sigma(0) is quantitatively established, the paper would provide the first numerical confirmation of the non-renormalization of the 1D chiral effect in a tight-binding model with interactions, extending Ref. [34] beyond the noninteracting limit and lending support to the dimensional-reduction picture of the chiral magnetic effect. Strengths include the exact QMC treatment of the model, the independence of the lambda and sigma computations from the prior analytic result (no fitting parameters), the free-case consistency check of the delta-function weights, and the availability of the open-source BSS-QMC code [56] used for the simulations. The Model 2 results, for which no exact solution exists, are potentially the most valuable contribution. The main weakness is that the central comparison currently rests on unquantified analytic-continuation reconstructions; the evidence is suggestive but not yet at a standard that would support the strength of the stated conclusion.","major_comments":[{"comment":"The central claim rests on the visual coincidence of two stochastic analytic continuation (SAC) reconstructions, with no error bars, no propagated uncertainties, and no quantitative measure such as lambda(0)/sigma(0) with confidence intervals reported anywhere in the paper; the U=0 benchmark in Fig. 1 is the exact free-fermion result rather than a validation of the SAC pipeline itself. Because the kernels in Eqs. (13)-(14) become nearly tau-independent for frequencies below about 1/beta, the behavior of lambda(omega) and sigma(omega) in the immediate vicinity of omega=0 at beta=4 is an extrapolation governed by the SAC entropy prior rather than a direct observable, and two independent ill-posed inversions could be smoothed toward similar low-frequency shapes even if the true zero-frequency weights differ. To make the equality load-bearing, the authors should provide (i) statistical error bars on the low-frequency parts of lambda(omega) and sigma(omega), for instance from bootstrap resampling of the QMC correlators and multiple SAC runs with different seeds and default models; (ii) an explicit validation of the same SAC pipeline on a case with a known answer, such as the U=0 model with analytically known delta-function weights or the exactly solvable Model 1 at small U; and (iii) a quantitative statement of the ratio lambda(omega_min)/sigma(omega_min) or an equivalent low-frequency integral ratio with its uncertainty.","section":"§V, Eqs. (13)–(14), Figs. 2–12"},{"comment":"In the Model 2, U=4 case the reconstructed single-particle gap is approximately 1.8, well above the temperature T=0.25, and both lambda(omega) and sigma(omega) vanish at small omega; the equality lambda(0)=sigma(0) reduces to the trivial 0=0 and provides no constraint on the possible renormalization of the chiral effect. The paper itself acknowledges in §VI that the linear-response anomaly formula no longer applies in this regime, yet the abstract states that 'even in the regime where sufficiently strong interactions drive the system into a Mott insulating phase, the proportionality ... remains unchanged,' which overstates what the computation can show. The nontrivial content lies in the cases with nonzero low-frequency weight (Model 2 at U less than or similar to 1 and the thermally activated Model 1 cases); the authors should either drop the insulator case from the central claim or explicitly label it as a null test that is consistent with, but does not test, the chiral effect.","section":"§V.B, §VI, Figs. 11–12; Abstract"},{"comment":"The chiral charge density is taken as the bare lattice operator rho_5 = psi^dagger sigma_2 psi and is assumed to remain the physical chiral density in the interacting model; possible interaction-induced operator renormalization (mixing with other local operators with the same quantum numbers) is not discussed. Since the conclusion 'the chiral effect is not renormalized' is a statement about this specific operator, the paper should state explicitly that the bare lattice operator is the operative definition and, ideally, discuss whether lattice Ward identities fix its normalization in the interacting case; otherwise the measured equality could reflect a particular choice of operator rather than a property of the physical chiral charge.","section":"§IV.A, Eq. (12), §VI"},{"comment":"The spectral function lambda(omega) of the axial-density-current correlator is sign-changing (it is visibly negative at intermediate frequencies in Figs. 3, 9, and 11), while the standard stochastic analytic continuation method cited in Refs. [57-59] assumes a positive-definite spectral function. The paper does not describe how its SAC implementation handles a non-positive spectral function (for example, decomposition into positive and negative parts, or a modified prior), nor does it state the default model and annealing parameters. Without this documentation the negative excursions and, more importantly, the reconstructed low-frequency behavior of lambda(omega) are not auditable; the authors should specify the implementation and justify its validity for sign-changing spectral functions.","section":"§IV.A, Eqs. (13)–(14), Figs. 3, 9, 11"}],"minor_comments":[{"comment":"The caption of Fig. 1 should state explicitly that the U=0 curves are the exact binned delta-function spectra and not SAC reconstructions, since this is the only point of comparison with an exact result.","section":"Fig. 1 caption"},{"comment":"The manuscript does not state the boundary conditions (periodic versus open) for the L=100 lattice or the treatment of the boundary links in the Peierls-substituted current operator; this information is needed to reproduce the free-fermion delta-function weights and the interacting results.","section":"§IV.A"},{"comment":"The simulation parameters (number of QMC sweeps, imaginary-time discretization Delta_tau, number of SAC sampling steps, and annealing parameters) are not given; a short methods paragraph or table is needed for reproducibility of the central result.","section":"§IV.A"},{"comment":"The phrase 'numerically exact Quantum Monte Carlo simulations' could be misread as applying to the reconstructed spectra; the QMC data are statistically exact, but the spectral functions are obtained from an approximate inversion, and the text should qualify this distinction.","section":"§VI"},{"comment":"Reference [2] and Reference [5] are the same article (Z. Qin et al., Phys. Rev. B 108, 195103 (2023)) and are listed twice with different citation numbers; Reference [7] is missing its author list.","section":"References"},{"comment":"Typos and formatting issues include 'thechiral effect' (page 2), '1DDirac semimetal' in the running title, 'chirality matrix' in §VI, and the ungrammatical 'for simplicity of expressions units' in §III.","section":"Throughout"},{"comment":"The normalization convention in Eq. (15) (cosh/cosh kernel without an explicit factor of 1/2) should be checked against the standard fermionic spectral representation, since the reported gap values in Figs. 4-6 and 13 depend on the absolute scale of A(omega,k).","section":"Eq. (15)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of cond-mat.str-el and the underlying QMC computation appears sound, but the evidentiary standard for the central claim is currently below what a quantitatively-minded journal would require, mainly because of the unquantified analytic-continuation step and the triviality of the deep-insulating comparison. The citation pattern is concentrated on the authors' own prior works (Refs. [14], [15], [34], and others), and the paper does not clearly delineate which part of the result is new relative to the exactly solvable Model 1 case; these are presentation issues rather than correctness concerns. Should the authors provide error bars and a validated SAC benchmark, I would be inclined to support acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plain take: this is a modest but legitimate numerical extension of the same group's earlier analytical result. The genuinely new thing is the first QMC check that the 1D chiral relation j = ev_F rho_5 survives local Hubbard interactions for the studied parameters. It does not open a new direction, but it is an honest, mostly careful piece of work.\n\nWhat it does well: The QMC is numerically exact for the model, the code is public, and the free-case spectra recover the expected delta-peak structure. The paper considers both spinful and spinless models, and it is candid that in the strongly gapped U=4 Model 2 case both response functions vanish, so that point is not evidence. The abstract appropriately limits the claim to \"within numerical resolution and for the parameters studied.\" The topological-invariant discussion is also honest that the invariant vanishes once the gap opens.\n\nThe soft spots are real but not disqualifying. The central comparison lambda(omega) ~ sigma(omega) comes from stochastic analytic continuation of two different imaginary-time correlators, and the paper gives no error bars and no validation on an exactly known interacting case. SAC broadens delta functions in the free case; if the two reconstructions are smoothed differently, the apparent low-frequency agreement could be an artifact. That concern deserves a direct test, for instance comparing reconstruction against exact diagonalization on a small system or using a known sum rule. The U=4 Model 2 case is indeed 0=0 and should be presented as such. There is also an overreach in the conclusions: the statement that the chiral effect is \"not renormalized by any kinds of interactions\" goes beyond the local density-density interactions actually studied. The bare operator rho_5 is assumed to remain the physical chiral density; interaction-induced mixing with other local operators is not discussed. Heavy self-citation is present, but the QMC calculation is independent of [34], so this is not circular.\n\nWho is this for? People working on the chiral magnetic effect in one-dimensional systems, dimensional reduction of anomalies, and QMC spectral-function methods. It is not going to change the field, but it is a useful numerical confirmation. I would send it to peer review, but ask for a revision that adds quantitative uncertainty estimates, a validation benchmark for the analytic continuation, and a more restrained conclusion.","headline":"A solid, modest numerical check that the 1D chiral relation survives local Hubbard interactions, undermined mainly by unquantified analytic-continuation errors and an overbroad concluding claim.","tokens_in":11749,"tokens_out":2338,"would_cite":true,"duration_ms":21580,"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":"Interacting 1D Dirac semimetals preserve the chiral effect, quantum Monte Carlo shows.","keywords":["chiral effect","chiral magnetic effect","1D Dirac semimetal","Su-Schrieffer-Heeger model","Hubbard interaction","quantum Monte Carlo","axial charge density","electric conductivity"],"falsifier":"Compute the exact spectral functions of a small interacting SSH chain by exact diagonalization, then apply the same stochastic analytic continuation to the exact correlators; if the reconstructed $\\lambda(\\omega)$ and $\\sigma(\\omega)$ deviate from the exact ones differently at low frequency, the observed equality in the paper could be an artifact of the continuation method.","tokens_in":10837,"feed_emoji":"⚡","tokens_out":8715,"duration_ms":71269,"temperature":0.7,"pith_summary":"The paper asks whether the chiral effect—the proportional relation between axial charge density and electric current in a 1D Dirac semimetal—survives when electrons interact. Using exact quantum Monte Carlo simulations of the Su-Schrieffer-Heeger model with two versions of local Hubbard interactions, the authors find that the proportionality remains the same as in the noninteracting case, within numerical accuracy and for the parameters studied. This holds for both a spinful and a spinless model, including interaction strengths that open a Mott gap. The result indicates that the chiral effect is not renormalized by local density-density interactions, suggesting a possible topological origin.","feed_headline":"Interacting 1D Dirac semimetals preserve the chiral effect","feed_subtitle":"Exact simulations show the axial-charge/current relation survives Hubbard interactions, even in a Mott-gapped phase.","key_machinery":"The key identity is $j = ev_F \\rho_5$, where $j$ is the electric current density and $\\rho_5 = \\rho_R - \\rho_L$ is the chiral (axial) charge density. In the noninteracting SSH model, this identity follows from the chiral anomaly combined with a relaxation-time approximation. The paper tests the identity in the interacting case by evaluating, from quantum Monte Carlo data, the frequency-dependent spectral functions $\\lambda(\\omega)$ and $\\sigma(\\omega)$ associated with the axial-charge–current and current–current correlators, and checking their equality at low frequency. The relation is expected to be governed by the topological invariant $N$ of Eq. (17) in the gapless regime.","core_discovery":"The central claim is that the relation $j = ev_F \\rho_5$, which equates the electric current with the axial charge density times the Fermi velocity and charge, holds in the interacting 1D Dirac semimetal. The authors verify this by computing the static linear-response coefficients for axial density and electric current via Kubo formulas from imaginary-time correlators, extracting the spectral functions $\\lambda(\\omega)$ and $\\sigma(\\omega)$ through stochastic analytic continuation, and showing that they coincide in the $\\omega \\to 0$ limit. This holds for the SSH model with both spinful and spinless Hubbard-type interactions, at temperatures where the interaction-induced gap is invisible or small, and even when the gap exceeds the temperature in the spinless model. Thus, within numerical accuracy, the chiral effect is unrenormalized by local interactions.","pith_inferences":["One could test whether $\\lambda(\\omega) = \\sigma(\\omega)$ holds at all frequencies, not just in the $\\omega \\to 0$ limit; if it does, the identity may stem from a general Kubo-type relation analogous to the non-renormalization of the chiral anomaly.","Because stochastic analytic continuation is an ill-posed inversion, a validation against exactly solvable small systems or a different continuation algorithm would make the equality more convincing; the paper does not include such a check.","The bare chiral density operator $\\rho_5$ is assumed to remain the correct observable in the interacting model; if interactions renormalize it, the measured equality may correspond to a different effective operator.","The study is restricted to local density-density interactions; whether nonlocal or spin-exchange interactions also preserve the relation is open."],"forward_implications":["If the relation $j = ev_F \\rho_5$ is universal, the chiral effect in 1D Dirac semimetals is a general feature that does not depend on the strength of local interactions.","The equality also holds in a gapped Mott phase, implying that the low-frequency linear response of both axial density and electric current vanish in a coordinated way; this could be probed in cold-atom realizations of the SSH model.","The result provides a numerical benchmark against which nonperturbative theories of the chiral magnetic effect in interacting systems can be tested.","It strengthens the case that the relation has a topological origin, as conjectured from the noninteracting analysis."],"supporting_citations":[{"why":"Defines the SSH model that is the subject of the study.","marker":"[1]"},{"why":"Establishes the $j = ev_F \\rho_5$ relation for the continuum 3D Dirac semimetal via Keldysh technique, forming the theoretical basis that this paper tests at the lattice level.","marker":"[14]"},{"why":"Proves the relation analytically for the noninteracting SSH model; this paper seeks to verify the same relation when interactions are added.","marker":"[34]"},{"why":"Supplies the auxiliary-field method at the heart of the quantum Monte Carlo simulations of the interacting model.","marker":"[53]"},{"why":"Provides the specific BSS-QMC algorithm used to sample the interacting Hamiltonian.","marker":"[55]"},{"why":"Provides the stochastic analytic continuation method used to extract the spectral functions $\\lambda(\\omega)$ and $\\sigma(\\omega)$ from imaginary-time correlators.","marker":"[58]"}],"fun_headline_variants":["Hubbard interactions don't alter chiral effect in 1D Dirac","Chiral effect persists under Hubbard interactions in 1D Dirac","Interactions leave chiral effect intact in 1D Dirac semimetal","1D Dirac chiral effect unchanged by Hubbard interactions","Chiral effect unrenormalized by Hubbard in Dirac semimetal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The equality of the two spectral functions at low frequency rests on the reconstruction of noisy imaginary-time Monte Carlo data by stochastic analytic continuation, an ill-posed inverse problem, and on the assumption that the bare chiral charge density operator remains the correct definition after interactions are turned on.","fun_headline_variants_meta":{"raw":{"variants":["Hubbard interactions don't alter chiral effect in 1D Dirac","Chiral effect persists under Hubbard interactions in 1D Dirac","Interactions leave chiral effect intact in 1D Dirac semimetal","1D Dirac chiral effect unchanged by Hubbard interactions","Chiral effect unrenormalized by Hubbard in Dirac semimetal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000935,"raw_usage":{"total_tokens":3967,"prompt_tokens":881,"completion_tokens":3086,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":2997}},"tokens_in":497,"tokens_out":3086,"duration_ms":18876,"temperature":1.0,"reasoning_tokens":2997,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:13:43.371874+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact spectral functions of a small interacting SSH chain by exact diagonalization, then apply the same stochastic analytic continuation to the exact correlators; if the reconstructed $\\lambda(\\omega)$ and $\\sigma(\\omega)$ deviate from the exact ones differently at low frequency, the observed equality in the paper could be an artifact of the continuation method.","supporting_citations":[{"cited_title":"A version of the calculation is represented in Ref","cited_arxiv_id":null,"evidence_quote":"Defines the SSH model that is the subject of the study."},{"cited_title":"Selch, M","cited_arxiv_id":null,"evidence_quote":"Proves the relation analytically for the noninteracting SSH model; this paper seeks to verify the same relation when interactions are added."},{"cited_title":"Impurity and dispersion effects on the linear magnetoresistance in the quantum limit","cited_arxiv_id":"2212.00383","evidence_quote":"Provides the specific BSS-QMC algorithm used to sample the interacting Hamiltonian."},{"cited_title":"White, D","cited_arxiv_id":null,"evidence_quote":"Provides the stochastic analytic continuation method used to extract the spectral functions $\\lambda(\\omega)$ and $\\sigma(\\omega)$ from imaginary-time correlators."}],"review_version":1}