{"id":"5698ddf4-9f84-4b16-9694-0d04b0c10566","arxiv_id":"2505.00861","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A path-integral, coherent-state derivation yields an exact non-Markovian stochastic master equation for any system linearly coupled to a harmonic lattice, demonstrated on the Fröhlich model.","lead":"This paper derives a stochastic master equation for a quantum system coupled to a vibrating crystal lattice, treating lattice vibrations as coherent waves rather than particles. It applies the method to the Fröhlich model and shows that zero-point lattice motion causes temperature-independent electron scattering.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Exactness of Eq. (8) turns on the complex-noise construction in S2.B: the covariances (7) are admitted not positive semidefinite, and no check shows the FFT noise preserves the required influence functional.","rationale":"The reader's weakest_assumption (single-particle wavepacket vs Fermi sea) concerns the physical interpretation of the Fröhlich application; it is real but does not attack the methodological exactness claim. I focused on the step where the central claim is least secure: the passage from the Feynman–Vernon influence functional to a well-defined stochastic process. The main text and S1 derive the influence functional in a standard way; from there Eq. (S28) is formally the usual Hubbard–Stratonovich transformation. The fracture point is the complex-noise extension, which the paper explicitly says involves freely chosen unphysical covariances. Nothing in the manuscript verifies that the generated noise satisfies the moment conditions that make the identity exact. This is not a disagreement with community consensus; it is an unproven step inside the derivation. Since the step is likely fixable (e.g., by providing the positive-semidefinite covariance construction used in Ref. [20] and a numerical consistency check), the appropriate verdict remains CONDITIONAL, so I set UNCHANGED relative to the reader.","tokens_in":21452,"tokens_out":12694,"duration_ms":138579,"concrete_test":"Set up the smallest exact benchmark: one bath oscillator (or two) with a system such as a harmonic oscillator or two-level system, linear coupling g_q(x), initial bath coherent state, and time-evolve the full composite Hamiltonian exactly over several oscillator periods. Generate η, ν with the S2.B FFT method using the chosen 'extra' covariances, propagate the stochastic master equation (8) for ~10^4–10^5 realizations, and compare the averaged reduced density matrix ⟨ρ̃_S⟩_W element-wise with the exact reduced dynamics. If the error does not shrink toward numerical tolerance as the number of realizations grows, the complex-noise construction does not reproduce the influence functional. A more direct version is to estimate E_W exp(i∫η·v + ν·u) from the generated noise for fixed test functions v, u and compare it numerically to exp(-Ψ) computed from Eq.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central exactness claim depends entirely on the Hubbard–Stratonovich identity (S28): for the averaged density matrix ρ_S = ⟨ρ̃_S⟩_W to equal the true reduced density matrix, the Gaussian noise W must have exactly the second moments (7)/(S27). Supplement S2.B concedes that the specified covariances are not positive semidefinite ('thus yielding the ill-defined normal distribution') and remedies this by making η, ν complex and 'freely' choosing extra, unphysical covariances so that the overall covariance becomes positive semidefinite. What is not shown is that this complex extension leaves the characteristic functional E_W exp(i∫ η·v+ν·u) equal to the Feynman–Vernon influence phase. Since the system propagator is a nonlinear functional of the noise, unconstrained changes to second-order statistics can change the averaged density matrix. The FFT recipe (S50) is only as good as the covariance matrix fed into it; if that matrix is chosen merely for positive semidefiniteness rather than to satisfy (7), then Eq. (8) is an exact master equation for a different stochastic process, not necessarily for the original linear-coupling model. The Fröhlich results would then not benchmark Eq. (8). This is an internal-consistency gap in the proof of the paper's central claim, independent of the physical applicability issues.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a path-integral formalism for a 'quantum acoustics' program in which a harmonic lattice bath is represented in a coherent-state (wave) picture rather than a Fock-state (particle) picture. The central result is a non-Markovian stochastic master equation, Eq. (8), for a system coupled linearly to bath coordinates through arbitrary system-dependent functions g_q(r). The derivation follows the Feynman-Vernon influence-functional route and uses a Hubbard-Stratonovich transformation to replace the bath influence by Gaussian noises η_q, ν_q with prescribed covariances (Eq. 7). The authors claim the equation is exact beyond weak coupling and the Markovian approximation, and they demonstrate the method on the Fröhlich model, comparing momentum relaxation times and wavepacket spreads obtained from the stochastic master equation, from a pure mean-field treatment, and from perturbation theory for Copper (weak coupling) and Bi2212 (strong coupling). The Supplement contains the full derivation, simulation details, and the perturbation-theory formulas used as benchmarks.","tokens_in":21602,"tokens_out":2247,"duration_ms":25919,"significance":"If the central exactness claim holds, the paper provides a useful and nontrivial extension of the Stockburger–Grabert stochastic Liouvillian machinery to baths with non-bilinear coupling (linear in bath coordinates but arbitrary in system coordinates), including coherent and discrete-mode baths. This fills a genuine gap in the open-quantum-systems toolbox and offers a potentially powerful non-perturbative simulation method for electron-lattice dynamics. The application to the Fröhlich model is plausible and the mean-field and perturbation-theory comparisons provide independent consistency checks. The paper also makes explicit that the method reproduces the spontaneous-emission contribution absent from the mean-field treatment, which is a physically meaningful and checkable prediction. However, the exactness proof is incomplete at a load-bearing point: the complex-noise extension in Supplement S2.B, which is needed because the required noise covariances are not positive semidefinite, is not shown to preserve the influence functional that defines the original model. This gap, if unresolved, would leave Eq.","major_comments":[{"comment":"The comparison between the stochastic master-equation results and the perturbation-theory benchmark is not fully quantitative. The reader is told that τ is extracted from an exponential fit of the momentum decay and that the analytic formulas (S60)-(S61) serve as the perturbative reference, but the figure shows separate curves for the numerical simulation and the analytic approximation with no error bars, no noise realizations count, and no convergence study with respect to either the number of trajectories or the time step. The claims that the stochastic and mean-field methods 'align well' and that the difference is due solely to spontaneous emission would be far more convincing with a convergence analysis in the number of noise realizations and a quantitative goodness-of-fit measure between the simulation and the analytic perturbation-theory curves.","section":"Section III.A and Figure 2"}],"minor_comments":[{"comment":"The caption of Figure 3 says the dashed line indicates the Debye temperature, but the legend does not explicitly show the dashed line. Please make sure each line and symbol is identified in the caption and legend.","section":"Figure 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the proposed method is potentially valuable, but the central exactness claim is not proven because of the complex-noise issue in Supplement S2.B. This is a technical point that could be fixed with a rigorous argument or by clearly labeling the result as a stochastic unraveling with an approximate noise construction. The single-particle limitation for metals is honestly discussed by the authors, but the presentation in Section III.A may still overstate the applicability to Copper and Bi2212. With a substantial revision that addresses the exactness proof and reframes the physical application, the paper could merit publication in a strong quantum-physics journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: competent, honest application of Stockburger's stochastic Liouville-von Neumann method to the Fr\\\"ohlich model, but the central exactness claim currently rests on an unproven step in the complex-noise construction. Worth refereeing, not worth accepting as-is.\n\nWhat's good: the path-integral derivation follows Feynman-Vernon and the supplement is thorough. The Fr\\\"ohlich application is concrete, and the interpretation of zero-point fluctuations as spontaneous emission is physically clear and connects naturally to the mean-field deformation-potential picture. The authors are transparent about the low-temperature behavior being unphysical for metals, which is a point in their favor.\n\nThe main problem is the exactness proof. Equation (8) hinges on the Hubbard-Stratonovich identity (S28), which requires the Gaussian noise to have precisely the covariances (S27). The supplement concedes those covariances are not positive semidefinite, so the noise is made complex and the extra \"unphysical\" covariances are chosen freely to restore positive semidefiniteness. That is a legitimate strategy, but the paper never proves that the resulting noise still satisfies the characteristic functional required by (S28). The FFT recipe (S50) only propagates the covariance you feed it; if that covariance no longer matches (S27), the averaged density matrix is the solution for a different stochastic process. This is an internal-consistency gap in the proof of the central claim, independent of the applications. It may well be fixable, but it has to be shown.\n\nSecond, the numerical validation compares the single-particle master equation against perturbation theory that includes Fermi statistics. Those are different models. At low temperature they disagree, and the paper explains this as fermiology suppressing spontaneous emission, but that means the master equation is not actually benchmarked in that regime. The high-temperature agreement is fine, but it is the regime where mean-field and full dynamics nearly coincide, so it is a weak test.\n\nThird, the novelty claim that non-bilinear coupling with coherent baths \"has not been explored\" is overstated. Feynman-Vernon and SLN methods already accommodate arbitrary system coupling operators; the genuinely new piece is the coherent-state bath specialization and the Fr\\\"ohlich application, which is enough.\n\nWho it's for: people working on stochastic master equations, electron-phonon dynamics, and the quantum-acoustic program. A useful incremental contribution once the complex-noise proof is supplied. My recommendation: send to peer review, with the requirement that the authors close the exactness gap and either benchmark against a single-electron calculation without Pauli blocking or include the blocking in the master equation.","headline":"The master equation is a coherent-state specialization of Stockburger's SLN method; the paper is worth refereeing, but the exactness claim rests on an unproven complex-noise construction and the numerical benchmark compares two different models.","tokens_in":22234,"tokens_out":3908,"would_cite":true,"duration_ms":43777,"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 coherent-state path integral yields a stochastic master equation claimed to be exact for any system linearly coupled to a harmonic lattice, demonstrated on the Fröhlich model.","keywords":["quantum acoustics","coherent states","stochastic master equation","path integral","Hubbard-Stratonovich transformation","non-Markovian dynamics","electron-phonon interaction","Fröhlich model"],"falsifier":"Compare the noise-averaged density matrix from Eq. (8) against an exact solution of the same model in a case where one is available—for instance, a single system mode linearly coupled to one or two harmonic bath modes solved by direct numerical integration of the composite Schrödinger equation; if the averages disagree for any such benchmark, the exactness claim fails.","tokens_in":21134,"feed_emoji":"⚛️","tokens_out":9112,"duration_ms":84467,"temperature":0.7,"pith_summary":"This paper tries to establish a wave-based, rather than particle-based, description of lattice vibrations—quantum acoustics—in which phonons behave like photons and acoustic waves like electromagnetic waves. Its central claim is that a non-Markovian stochastic master equation derived from a coherent-state path integral captures the exact reduced dynamics of any system coupled linearly to a harmonic lattice bath, with no restriction to weak coupling, Markovian behavior, or continuous spectral densities. The authors demonstrate the procedure on the Fröhlich model of electron-phonon interaction, computing momentum relaxation and wavepacket spreading in weak- and strong-coupling regimes and comparing with mean-field and perturbation theory. If the claim holds, the framework supplies a first-principles, non-perturbative tool for electron-lattice dynamics and places quantum acoustics on the same footing as quantum optics.","feed_headline":"Exact master equation promises wave view of electron-lattice dynamics","feed_subtitle":"Coherent-state path integrals turn lattice vibrations into a stochastic potential beyond weak coupling.","key_machinery":"The machinery is the coherent-state path integral for the lattice bath, with the bath initialized as a multimode coherent state (or a thermal ensemble of such states). The bath's effect on the system is packaged in an influence functional, the phase factor that records how the bath responds to the system's history; this functional is split into a mean-field action $S_{\\mathrm{mf}}$ linear in the coupling and an influence phase $\\Psi$ quadratic in the coupling. A Hubbard–Stratonovich transformation introduces two Gaussian noise fields $\\eta_q(t)$ and $\\nu_q(t)$ whose covariances are fixed by the lattice propagators; substituting them turns the influence phase into averaged noise actions and yields the stochastic master equation. The operational object is the stochastic pseudopotential $H^{\\pm}_{B+I} = H_{\\mathrm{mf}} - \\sum_q(\\eta_q(t)\\cdot g_q \\mp \\frac{1}{2}\\nu_q(t)\\cdot g_q)$, which turns the problem into two non-Hermitian Schrödinger equations per noise realization; averaging over noise recovers the physical density matrix.","core_discovery":"The authors' central discovery is the stochastic Liouville–von Neumann equation (Eq. 8), $i\\partial_t \\tilde{\\rho}_S = [H_S + H_{\\mathrm{mf}}, \\tilde{\\rho}_S] - \\sum_q( [\\eta_q(t)\\cdot g_q, \\tilde{\\rho}_S] - \\frac{1}{2}\\{\\nu_q(t)\\cdot g_q, \\tilde{\\rho}_S\\})$, whose noise average $\\rho_S = \\langle \\tilde{\\rho}_S \\rangle_W$ reproduces the exact reduced density matrix of the system. The derivation starts from a multimode coherent-state density matrix for the lattice and the Feynman–Vernon influence functional, decomposes the latter into a linear mean-field action and a quadratic influence phase, and applies the Hubbard–Stratonovich transformation to convert the phase into Gaussian noise. The paper claims this master equation is exact for any system with coupling linear in the bath coordinates and nonlinear in the system coordinates, beyond weak coupling and Markovian approximation, and that it can be unraveled into two independent Schrödinger equations with non-Hermitian Hamiltonians, allowing wave-packet propagation in real space. Applied to the Fröhlich model, it yields relaxation times and spatial spreads that agree with mean-field and perturbation theory at high temperature, while exposing spontaneous emission from zero-point lattice fluctuations at low temperature.","pith_inferences":["The same master equation should apply to optical phonons and to molecular or cavity vibrational modes, since the derivation only requires coupling linear in bath coordinates; testing those settings would widen the quantum-acoustics analogy beyond acoustic deformation-potential coupling.","A natural next step is to replace the single Gaussian wavepacket with a multi-electron initial state that respects fermionic statistics, which would test whether the low-temperature spontaneous-emission saturation survives Pauli blocking and would turn the reduced-model results into material predictions.","Because the stochastic potential is additive over modes, the approach is naturally parallelizable and could be combined with existing wavefunction propagation codes; a practical speedup over quantum Monte Carlo or hierarchical-equation methods is a plausible but unproven consequence."],"forward_implications":["Any system with coupling linear in the bath coordinates and nonlinear in the system coordinates can be propagated non-perturbatively, without weak-coupling or Markovian assumptions, by sampling the noise and averaging the resulting wavefunctions.","In the weak-coupling limit the mean-field deformation-potential picture is recovered, because the mean-field action is linear in the coupling while the influence phase is quadratic; the master equation therefore reduces to the established quasi-classical description as the coupling shrinks.","Zero-point lattice fluctuations generate spontaneous emission that the mean-field approximation omits, so at low temperature the exact momentum relaxation time saturates to a finite value while the mean-field value vanishes.","Because the formalism accommodates arbitrary dispersion, discrete or coherent baths, and any system Hamiltonian, it extends beyond the Fröhlich model to optical modes, molecular vibrations, or other linear-coupling settings without a Lee–Low–Pines restriction to a vacuum bath.","Benchmarking on the Fröhlich model shows the method tracks perturbation theory and mean-field results where expected, giving a validated baseline for entering the non-perturbative strong-coupling regime where no perturbative benchmark exists."],"supporting_citations":[{"why":"supplies the stochastic Liouville–von Neumann method that this work extends to non-bilinear coupling with coherent and discrete baths.","marker":"[19]"},{"why":"provides the complex-noise covariance construction and fast-Fourier-transform sampling used to generate the Gaussian noise.","marker":"[20]"},{"why":"gives the Feynman–Vernon influence-functional formalism underlying the path-integral propagator for the reduced density matrix.","marker":"[27]"},{"why":"the Hubbard–Stratonovich transformation converts the quadratic influence phase into noise averages and yields the stochastic master equation.","marker":"[28, 29]"},{"why":"establishes the coherent-state deformation-potential picture of lattice vibrations and provides the mean-field baseline the paper compares against.","marker":"[10]"},{"why":"define the Fröhlich model of electron-phonon coupling used as the demonstration system.","marker":"[22, 23]"},{"why":"supplies the standard separable-initial-state and reduced-density-matrix framework on which the derivation builds.","marker":"[24]"}],"fun_headline_variants":["Wave paradigm yields exact non-Markovian master equation","Path integrals unlock wave view of quantum acoustics","Exact master equation from coherent-state path integrals","Beyond particles: exact dynamics for lattice-coupled systems","Quantum acoustics meets path integrals: exact stochastic dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that one Gaussian wavepacket launched at the Fermi momentum stands in for a real charge carrier, even though the model neglects the occupied electron states around it and the exclusion principle that would block many decays; the paper itself says the resulting low-temperature spontaneous emission is unphysical for metals.","fun_headline_variants_meta":{"raw":{"variants":["Wave paradigm yields exact non-Markovian master equation","Path integrals unlock wave view of quantum acoustics","Exact master equation from coherent-state path integrals","Beyond particles: exact dynamics for lattice-coupled systems","Quantum acoustics meets path integrals: exact stochastic dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000285,"raw_usage":{"total_tokens":1674,"prompt_tokens":937,"completion_tokens":737,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":663}},"tokens_in":553,"tokens_out":737,"duration_ms":7297,"temperature":1.0,"reasoning_tokens":663,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:34:19.320394+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the noise-averaged density matrix from Eq. (8) against an exact solution of the same model in a case where one is available—for instance, a single system mode linearly coupled to one or two harmonic bath modes solved by direct numerical integration of the composite Schrödinger equation; if the averages disagree for any such benchmark, the exactness claim fails.","supporting_citations":[{"cited_title":"Feynman and A","cited_arxiv_id":null,"evidence_quote":"supplies the stochastic Liouville–von Neumann method that this work extends to non-bilinear coupling with coherent and discrete baths."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the complex-noise covariance construction and fast-Fourier-transform sampling used to generate the Gaussian noise."},{"cited_title":"Weiss, Quantum dissipative systems (World Scientific, 2012)","cited_arxiv_id":null,"evidence_quote":"gives the Feynman–Vernon influence-functional formalism underlying the path-integral propagator for the reduced density matrix."},{"cited_title":"Fr¨ ohlich and N","cited_arxiv_id":null,"evidence_quote":"supplies the standard separable-initial-state and reduced-density-matrix framework on which the derivation builds."}],"review_version":1}