{"id":"f6645204-c40d-4e97-abad-55b6f049e936","arxiv_id":"2501.16856","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In disordered bosonic systems, the 2D spectroscopy echo peak is controlled by the ratio of elastic to inelastic self-energies, and interaction-induced quantum fluctuations add broadening that cannot be rephased.","lead":"This paper builds a quantum field theory framework for reading two-dimensional terahertz spectroscopy maps of disordered materials, and it shows how the echo peak can reveal whether scattering is elastic or inelastic. It also predicts that interactions make perfect echo rephasing impossible, a limit ignored by the standard two-level-atom picture.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main plots (Figs. 8-9) use V0^2/m^2=3 in d=2, violating the stated ladder validity condition (mvd/2/V0)^2 >> 1 of Sec. V A; central vertex relation Eq. (39) and non-rephasing prediction Eq. (40) may be uncontrolled at shown parameters.","rationale":"The reader's weakest_assumption identifies the ladder (non-crossing) resummation validity, and the paper's own validity condition in Sec. V A is indeed violated by the parameters used in the main plots. I checked the algebra leading to Eq. (39): given the momentum-independent self-energy assumption and the self-consistent Born solution, the derivation of Eq. (39) is internally consistent. I also verified that Eq. (40) can be obtained from the generalized self-consistent relation by identifying Q = 4V0^2 sum_k |D_k(omega)|^2 ImSigma_{g-V}(k,omega), so the formal result is not algebraically flawed. The concern is therefore not about the internal consistency of the ladder calculation but about whether the ladder approximation itself is controlled at V0^2/m^2 = 3 in d=2. The paper states that crossing diagrams can be disregarded when E tau_e >> 1, with tau_e^{-1} ~ V0^2/(v d m); for the plotted parameters E tau_e ~ (m/V0)^2 <= 2/3, so the central approximation is applied outside its stated regime. This is an internal inconsistency between the stated validity condition and the numerical demonstration. It does not necessarily invalidate the qualitative predictions, but it does undermine the quantitative reliability of the central quantitative relation at the shown parameters. The proposed test recomputes the same observable inside the stated validity regime and directly estimates the size of the neglected crossing contribution. If the qualitative almond shape and the ratio scaling persist at V0^2/m^2 = 0.01, and if the crossing contribution at V0^2/m^2 = 3 is small, the concern is resolved. In the absence of such checks, a CONDITIONAL verdict is appropriate, and my stress-test does not change the reader's verdict.","tokens_in":27519,"tokens_out":11668,"duration_ms":96846,"concrete_test":"Recompute the echo contribution to the 2D map at a disorder strength satisfying the stated ladder condition, e.g. V0^2/m^2 = 0.01 in d=2 (so (m/V0)^2 = 100), with the same small g and gamma/m ratios as in Fig. 8, and check whether the almond-like echo peak and the Eq. (39) scaling of the diagonal and cross-diagonal slices survive. Additionally, evaluate the first crossing (non-ladder) disorder diagram contribution to the vertex Gamma_dis at V0^2/m^2 = 3 and compare its magnitude on the echo diagonal with the ladder result; if it exceeds roughly 10% of the ladder vertex, the plotted regime is not covered by the paper's stated validity condition.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central analytical results, Eqs. (39)-(40), are derived within the self-consistent Born and ladder (non-crossing) resummation of disorder diagrams. In Sec. V A the authors state this resummation is valid only when (mvd/2/V0)^2 >> 1, which for d=2 and v=1 reduces to V0^2/m^2 << 1. However, the main 2D maps in Fig. 8 use V0^2/m^2 = 1.5 and 3, and Fig. 9 uses V0^2/m^2 = 3, giving (m/V0)^2 = 2/3 and 1/3, respectively. These values are not >> 1, so the elastic collision time is comparable to or shorter than the oscillation period. In this regime crossing disorder diagrams and weak-localization-type corrections, which the paper itself acknowledges in Sec. V A and lists as future work in Sec. VIII, contribute at order one. If such diagrams modify the disorder vertex on the echo diagonal, the exact-looking Eq. (39) and the quantum-fluctuation cutoff of the divergence in Eq. (40) become quantitatively unreliable for the plotted parameters. The qualitative almond shape may survive, but the claimed quantitative mapping of echo peak slices to the self-energy ratio ImSigma_V/ImSigma_bath is not controlled.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a Keldysh field-theoretic framework for computing the third-order nonlinear response, and hence two-dimensional terahertz spectroscopy (2DTS) maps, of bosonic collective excitations in the presence of static disorder. The authors use a Luttinger-Ward functional combined with a ladder (non-crossing) resummation of disorder lines to derive the disorder-dressed four-point vertex, and they specialize to the weak-nonlinearity regime where the echo diagonal vertex takes the closed form T_g-V(ω,-ω) = 2ig ImΣ_V(ω)/ImΣ_bath(ω) (Eq. (39)). They further include one class of interaction-disorder self-energy corrections to obtain Eq. (40), showing that interaction-induced quantum fluctuations produce a nonrephasable broadening that prevents perfect echo rephasing. Numerical 2D maps illustrate the almond-like echo peak, its evolution with disorder strength, and its suppression or enhancement depending on the sign of the interaction.","tokens_in":27834,"tokens_out":5528,"duration_ms":49384,"significance":"The paper's main strength is that Eq. (39) is a derived, parameter-free relation connecting the echo-diagonal vertex to the elastic-to-inelastic self-energy ratio; this goes beyond a fit and provides a concrete, potentially falsifiable prediction for 2DTS experiments such as those on Josephson plasmons in cuprates. The framework also correctly reproduces the ITLS echo phenomenology in a collective-mode language and identifies a genuinely many-body effect (quantum-fluctuation-induced nonrephasable broadening). If the central results are robust, this is a valuable contribution to the theory of nonlinear spectroscopy of disordered quantum materials. The main caveat is that the quantitative predictions are obtained within a specific non-crossing approximation, and the numerical demonstration is carried out at parameter values that lie outside the stated validity regime of that approximation.","major_comments":[{"comment":"The ladder resummation is stated to be valid for (mvd/2/V0)^2 >> 1, which for d = 2 and v = 1 reduces to (m/V0)^2 >> 1. The main echo maps use V0^2/m^2 = 1.5 and 3 in Fig. 8(a)-(b) and V0^2/m^2 = 3 in Fig. 9, corresponding to (m/V0)^2 = 2/3 and 1/3, respectively. These parameters are not in the stated validity regime, so the quantitative echo-peak slices in Fig. 8(c) and the quantum-fluctuation broadening shown in Fig. 9 are not controlled by the non-crossing approximation. The authors should either rerun the main figures at parameters satisfying (m/V0)^2 >> 1, provide a quantitative estimate of the neglected crossing-disorder-diagram corrections at the plotted parameters, or explicitly justify that the ladder result is quantitatively accurate beyond its nominal validity range.","section":"Sec. V A; Figs. 8 and 9"},{"comment":"The conclusion that perfect rephasing is 'fundamentally unattainable in interacting many-body systems' is stronger than what is demonstrated. Equation (40) is derived within a specific Luttinger-Ward functional (Fig. 7) that retains one mixed disorder-interaction self-energy term (Eq. (35)) while neglecting the disorder bubble (Eq. (37)) and the inelastic vertex corrections of Appendix D. The formula is therefore established for this class of non-crossing diagrams, not as a general theorem for all interacting many-body systems. Please either prove that Q(omega) remains nonzero for all conserving approximations that include momentum-dependent self-energy corrections, or explicitly scope the conclusion to the non-crossing approximation used in this work.","section":"Sec. VII, Eq. (40)"},{"comment":"The statement that Eq. (39) is 'exact within perturbation theory in g in all spatial dimensions' should be qualified. Within the chosen self-energy structure (Eqs. (25)-(33)) the algebra leading to Eq. (39) is indeed exact, but the vertex Gamma_dis itself comes from the ladder resummation and excludes crossing disorder diagrams. As written, the sentence could be read as a model-independent exactness claim, which is not supported by the derivation.","section":"Sec. VI, text after Eq. (39)"}],"minor_comments":[{"comment":"'non-pertubative' should be 'non-perturbative'.","section":"Abstract and Sec. I"},{"comment":"The frequency convention for the four-point vertex Gamma(omega3, omega2, omega1) is not specified at the point of introduction; please state explicitly that the arguments follow the shifted-frequency convention of Eq. (12) to avoid confusion.","section":"Sec. III B, Eq. (15)"},{"comment":"The caption uses V0^2/m^2 v^2 = 2 while the main text sets v = 1; a single dimensionless notation (e.g., V0^2/m^2) should be used consistently.","section":"Fig. 2 caption"},{"comment":"'caclulation' should be 'calculation'.","section":"Appendix C, first paragraph"},{"comment":"The regulator eta is described as a finite broadening, while later inelastic damping is denoted gamma; please clarify the relationship between eta and gamma to avoid notational confusion.","section":"Sec. IV B, Eq. (22)"},{"comment":"Reference [99] is listed as 'A. Gomez Salvador, unpublished'; if this is intended as a companion work, it should be replaced by a proper citation or removed.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The framework and the derived relation Eq. (39) are sound within the stated approximation, but the main numerical figures are at disorder strengths outside the ladder validity condition stated by the authors themselves. The revision should focus on re-running or justifying the parameter regime, and on scoping the universality claim in Sec. VII."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version. This is a real contribution, not a repackaging: it takes the authors' earlier disorder-free 2DTS formalism and properly adds static disorder within the Keldysh/Baym-Kadanoff framework, producing one clean analytical result (Eq. 39) and a physically interesting prediction (Eq. 40) that interaction-induced quantum fluctuations make perfect rephasing impossible. The algebra behind Eq. (39) is transparent and, given the self-energy structure, exact. The comparison with the experimental nonrephasing peak from their earlier paper has no fitting parameters.\n\nThe main soft spot is the one the stress-test flags. The ladder (non-crossing) resummation is justified only when (mvd/2/V0)^2 >> 1, which for v=1, d=2 means V0^2/m^2 << 1. The main 2D maps in Figs. 8 and 9 use V0^2/m^2 = 1.5 and 3, so the ratio is 2/3 or 1/3, not >> 1. In that regime crossing diagrams and weak-localization corrections, which the paper itself consigns to future work, enter at order one. So Eq. (39) and the nonrephasing broadening in Eq. (40) are exact within the ladder approximation but are not quantitatively controlled at the plotted parameters. The qualitative almond shape and its asymmetry may well survive—the physics is not obviously wrong—but the paper's claim that the echo peak slices quantitatively map to ImSigma_V/ImSigma_bath is overstated for those figures.\n\nA second, milder issue is the scope of the 'no perfect rephasing' conclusion. It follows from their specific Luttinger-Ward truncation and the definition of Q(omega) within that truncation. The statement in the abstract that perfect rephasing is 'fundamentally unattainable in interacting many-body systems' reaches further than the derivation sustains.\n\nOn the plus side: the derivation is careful, the paper is honest about what is neglected (crossing diagrams, localization, fifth-order terms), and the connection to ITLS intuition is genuinely illuminating. Self-citation to Refs. [90] and [62] is fine; they are extending their own prior work, and the new vertex relation is not in those papers. No code or data is shipped, which is a minor omission for a mostly formal paper.\n\nWho gets value: theorists and experimentalists working on 2D THz spectroscopy of collective modes, especially in cuprates. It deserves a serious referee. I'd recommend acceptance after major revision: either soften the claims to match the regime actually plotted, or add a discussion of why crossing diagrams are suppressed at those parameters. Send it out.","headline":"Careful extension of the authors' 2DTS formalism to static disorder, with a clean vertex relation (Eq. 39) and an interesting quantum-fluctuation prediction, but the main plots sit outside the stated ladder-validity regime.","tokens_in":28411,"tokens_out":3264,"would_cite":true,"duration_ms":28767,"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 Keldysh-based theory of two-dimensional spectroscopy shows that for disordered bosonic collective excitations the echo peak is governed by the elastic-to-inelastic self-energy ratio, and that interaction-induced quantum fluctuations…","keywords":["two-dimensional terahertz spectroscopy","disordered bosonic systems","collective excitations","Keldysh formalism","third-order nonlinear response","echo peak","elastic and inelastic scattering","quantum fluctuations"],"falsifier":"Measure the echo peak of a disordered bosonic collective mode and extract its cross-diagonal width as a function of the independently determined inelastic linewidth; if Eq. (39) holds, that width should follow the total self-energy while the echo sharpness tracks $V_0^2/\\omega\\gamma$ in the elastic-dominated regime. Alternatively, compute the disorder-dressed vertex including crossing diagrams at the parameters used in the paper's plots, $(mvd/2/V_0)^2=1/3$; a sizable deviation from $2ig\\,\\mathrm{Im}\\Sigma_V/\\mathrm{Im}\\Sigma_{\\mathrm{bath}}$ would show that the ladder assumption underpinning the central formula fails in that regime.","tokens_in":27320,"feed_emoji":"⚛️","tokens_out":15614,"duration_ms":119140,"temperature":0.7,"pith_summary":"This paper develops a theoretical framework for two-dimensional terahertz spectroscopy of bosonic collective excitations in a disordered medium, based on a Keldysh path-integral treatment that sums the infinite series of non-crossing disorder diagrams. In the weak nonlinearity limit it derives a closed formula for the disorder-dressed four-wave-mixing vertex on the echo diagonal, $T_{g-V}(\\omega,-\\omega)=2ig\\,\\mathrm{Im}\\Sigma_V(\\omega)/\\mathrm{Im}\\Sigma_{\\mathrm{bath}}(\\omega)$, so the echo peak directly encodes the competition between elastic and inelastic scattering. When interaction-induced quantum fluctuations are included, the same vertex obeys $T_{g-V}=2ig[-\\mathrm{Im}\\Sigma_V(\\omega)+Q(\\omega)]/[-\\mathrm{Im}\\Sigma_{\\mathrm{bath}}(\\omega)-Q(\\omega)]$, and the paper argues that even with no inelastic bath there is no divergence, meaning perfect rephasing cannot be achieved in interacting many-body systems. A reader would care because this gives a concrete experimental diagnostic for separating static disorder from intrinsic many-body dephasing in materials such as cuprates where that separation has been difficult.","feed_headline":"Echo peak pins elastic-to-inelastic ratio in disordered bosons","feed_subtitle":"Ties the echo's shape to the elastic-inelastic self-energy ratio and says quantum fluctuations block perfect rephasing.","key_machinery":"The load-bearing object is the disorder-dressed four-wave-mixing vertex $T_{g-V}(\\omega_a,\\omega_b)$, obtained by solving a Bethe-Salpeter-type ladder equation for the infinite series of non-crossing impurity lines, $\\Gamma_{\\mathrm{dis}}(0;\\omega_a,\\omega_b)=2V_0^2/[1-V_0^2\\lambda(0;\\omega_a,\\omega_b)]$, and then contracting it with the $\\phi^4$ interaction. The disorder bubble $\\lambda(0;\\omega_a,\\omega_b)$ is the crucial identity carrier: evaluated at $\\omega_a=-\\omega_b=\\omega$ it equals $V_0^{-2}\\mathrm{Im}\\Sigma_V(\\omega)/\\mathrm{Im}\\Sigma(\\omega)$, which converts the vertex into the self-energy ratio of Eq. (39). The Keldysh contour supplies the causal time-ordering of the two-pulse protocol, and the overall self-consistent construction ensures the self-energy and vertex come from one generating functional, so the approximation is conserving.","core_discovery":"The central claim is that the rephasing (echo) nonlinearity of a disordered bosonic collective mode is controlled by a disorder-dressed vertex whose value on the echo diagonal is, within the paper's conserving ladder approximation, the ratio of the elastic and inelastic parts of the single-particle self-energy: $T_{g-V}(\\omega,-\\omega)=2ig\\,\\mathrm{Im}\\Sigma_V(\\omega)/\\mathrm{Im}\\Sigma_{\\mathrm{bath}}(\\omega)$ with $\\mathrm{Im}\\Sigma_{\\mathrm{bath}}=-2\\omega\\gamma$. This recovers the isolated-two-level-system phenomenology, but with a many-body difference: the almond-shaped peak is asymmetric and elongated toward higher energies, since a zero-momentum collective mode cannot scatter into finite-momentum states below its mass. Including interaction-induced self-energy corrections, the paper obtains $T_{g-V}(\\omega,-\\omega)=2ig[-\\mathrm{Im}\\Sigma_V(\\omega)+Q(\\omega)]/[-\\mathrm{Im}\\Sigma_{\\mathrm{bath}}(\\omega)-Q(\\omega)]$, where $Q(\\omega)$ collects momentum-dependent quantum fluctuations; because $Q(\\omega)$ acts like an effective inelastic process, the vertex no longer diverges even when $\\Sigma_{\\mathrm{bath}}=0$. The paper concludes that perfect rephasing is fundamentally unattainable in interacting many-body systems, and that repulsive interactions suppress while attractive interactions enhance the echo.","pith_inferences":["If Eq. (39) holds beyond the specific $\\phi^4$ model, experimental 2DTS maps could be used to extract the elastic-to-inelastic scattering ratio directly from the echo peak shape, turning the protocol into a quantitative disorder probe for quantum materials.","The prediction that no perfect rephasing occurs when $Q(\\omega)\\neq 0$ suggests that the residual cross-diagonal width of an echo, measured with the inelastic bath tuned away, is a direct measure of interaction-induced dephasing in a disordered bosonic system.","A numerical exact-diagonalization study of a disordered Bose-Hubbard model across the ladder-validity boundary $(mvd/2/V_0)^2\\sim 1$ could test where crossing diagrams begin to modify the vertex and map the quantitative limits of the central formula.","Applying the same ladder machinery to fifth-order responses could reveal whether bound states of collective excitations rephase sharply or inherit the nonrephasable broadening, a question the paper leaves open."],"forward_implications":["The echo peak of a disordered bosonic collective mode carries a direct signature of the ratio $\\mathrm{Im}\\Sigma_V/\\mathrm{Im}\\Sigma_{\\mathrm{bath}}$: in the elastic-dominated regime the echo sharpens and takes an almond shape, while inelastic scattering broadens it.","When vertex corrections are significant, the diagonal and cross-diagonal cuts of the echo no longer coincide with the linear-response spectral function, so 2DTS provides information beyond linear response.","In the absence of an inelastic bath and of quantum fluctuations, the dressed vertex diverges on the echo diagonal, giving perfect rephasing analogous to a zero-momentum diffuson pole for the bosonic excitation.","Interaction-induced quantum fluctuations, encoded in $Q(\\omega)$, add a nonrephasable broadening even when $\\Sigma_{\\mathrm{bath}}=0$; hence perfect rephasing is fundamentally unattainable in interacting many-body systems.","The sign of the interaction matters: repulsive ($g>0$) interactions suppress the echo, while attractive ($g<0$) interactions enhance it."],"supporting_citations":[{"why":"supplies the experimental Josephson-plasmon echo spectra that motivate the model and serve as the qualitative comparison target for the almond-shaped echo peak.","marker":"[62]"},{"why":"provides the isolated-two-level-system echo lineshape and the diagonal/cross-diagonal slicing procedure that the weak-nonlinearity result is shown to recover.","marker":"[84]"},{"why":"gives the earlier mean-field third-order response of collective excitations that this paper's RPA calculation generalizes to the quantum limit.","marker":"[90]"},{"why":"documents the nonperturbative character of disorder diagrams and the conditions under which non-crossing ladder diagrams dominate.","marker":"[96]"},{"why":"supplies the diffuson-pole and weak-disorder diagrammatic picture used to interpret the divergent rephasing vertex.","marker":"[97]"},{"why":"sets the Keldysh path-integral conventions, causality structure, and Green's function definitions used throughout the response calculation.","marker":"[86]"},{"why":"introduces the conserving approximation scheme that keeps the self-energy and vertex derived from a single functional.","marker":"[91, 92]"}],"fun_headline_variants":["Echo peak exposes elastic-to-inelastic split in disordered bosons","Quantum fluctuations block perfect rephasing in bosonic echo","Rephasing impossible: quantum fluctuations break echo in disordered bosons","2D spectroscopy reveals why perfect rephasing fails in bosons","Echo asymmetry signals quantum many-body effects in disordered bosons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the ladder (non-crossing) resummation of disorder diagrams captures the rephasing physics, an assumption the paper states holds only for $(mvd/2/V_0)^2 \\gg 1$, while the illustrative two-dimensional plots use $V_0^2/m^2=3$, where that combination equals $1/3$.","fun_headline_variants_meta":{"raw":{"variants":["Echo peak exposes elastic-to-inelastic split in disordered bosons","Quantum fluctuations block perfect rephasing in bosonic echo","Rephasing impossible: quantum fluctuations break echo in disordered bosons","2D spectroscopy reveals why perfect rephasing fails in bosons","Echo asymmetry signals quantum many-body effects in disordered bosons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001074,"raw_usage":{"total_tokens":4510,"prompt_tokens":972,"completion_tokens":3538,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":3460}},"tokens_in":588,"tokens_out":3538,"duration_ms":22721,"temperature":1.0,"reasoning_tokens":3460,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T10:09:15.388124+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the echo peak of a disordered bosonic collective mode and extract its cross-diagonal width as a function of the independently determined inelastic linewidth; if Eq. (39) holds, that width should follow the total self-energy while the echo sharpness tracks $V_0^2/\\omega\\gamma$ in the elastic-dominated regime. Alternatively, compute the disorder-dressed vertex including crossing diagrams at the parameters used in the paper's plots, $(mvd/2/V_0)^2=1/3$; a sizable deviation from $2ig\\,\\mathrm{Im}\\Sigma_V/\\mathrm{Im}\\Sigma_{\\mathrm{bath}}$ would show that the ladder assumption underpinning the central formula fails in that regime.","supporting_citations":[{"cited_title":"Houver, L","cited_arxiv_id":null,"evidence_quote":"supplies the experimental Josephson-plasmon echo spectra that motivate the model and serve as the qualitative comparison target for the almond-shaped echo peak."},{"cited_title":"Zhang, F","cited_arxiv_id":null,"evidence_quote":"provides the isolated-two-level-system echo lineshape and the diagonal/cross-diagonal slicing procedure that the weak-nonlinearity result is shown to recover."},{"cited_title":"McGinley, M","cited_arxiv_id":null,"evidence_quote":"gives the earlier mean-field third-order response of collective excitations that this paper's RPA calculation generalizes to the quantum limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the diffuson-pole and weak-disorder diagrammatic picture used to interpret the divergent rephasing vertex."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"sets the Keldysh path-integral conventions, causality structure, and Green's function definitions used throughout the response calculation."}],"review_version":1}