{"id":"18f96f93-fae4-4505-b901-8aaba60024b2","arxiv_id":"1908.05809","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A new set of stochastic Ehrenfest relations for the full stochastic projected Gross-Pitaevskii equation, including energy damping and all projector terms, is derived and tested on center-of-mass motion, matching simulations.","lead":"This paper derives exact stochastic equations for how the center of mass, momentum, and other collective observables of a hot Bose-Einstein condensate evolve when the condensate exchanges particles and energy with a thermal reservoir. The equations, called stochastic Ehrenfest relations, give theorists a way to compute damping and noise in hot condensates analytically instead of only by simulation.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The COM validation tests only the simplified, projector-free OU limit of the SERs, so it does not actually exercise the exact projector and noise-correlation terms that constitute the paper's central claim.","rationale":"The paper's formal derivation is careful and the internal-consistency comparison is a reasonable first step, but the validation target is narrower than the central claim. The analytic COM solution necessarily drops the projector terms that the SERs are designed to include, and the numerical check of Appendix B.3 monitors only two of those terms. Consequently, even perfect agreement in Figs. 2-3 would not certify the exactness of (50a)-(50e) or the noise correlations (46)-(49). This is a concern about evidence rather than an identified algebraic error: the exact SERs may well be correct, but the paper's tests do not yet establish that. The reader's TIRA/Kohn's-theorem concern is valid and related, but it is not the only reason the validation is weak; the under-testing of the projector terms is more direct. A term-by-term increment test against SPGPE trajectories would settle the question cleanly. Since this supports the reader's CONDITIONAL verdict rather than overturning it, the verdict is left unchanged.","tokens_in":23475,"tokens_out":18880,"duration_ms":173733,"concrete_test":"Perform a term-by-term check of Eqs. (50a)-(50b) on equilibrated SPGPE trajectories: for an ensemble of trajectories, compute the actual finite-difference increments ΔR_j and ΔP_j over a small time step dt, and compare their ensemble means and covariances with the sums of all drift, trace, and projector terms in those equations plus the noise covariances (46)-(49), evaluated at the same field configurations. If the full SERs reproduce the SPGPE increments to numerical tolerance, the exactness claim is supported; if not, the discrepancy identifies the specific projector or noise term that is incorrect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eqs. (50a)-(50e) are exact for the full SPGPE, retaining all projector terms and the noise correlations (46)-(49). The reported validation, however, never directly exercises those exact equations. The analytic COM model is obtained by (i) imposing the Thomas-Fermi ansatz Eq. (53), (ii) linearizing in z(t), and (iii) neglecting all projector terms, after which the SERs reduce to the Ornstein-Uhlenbeck system (57)-(58). The only numerical check on the neglected terms is Appendix B.3, which tracks just two of the many projector contributions, qH_z and dzε, and argues that the rest are smaller because they carry small damping rates. Thus an algebraic error confined to, for example, the unmonitored qγ_z, dγ_z, or the diffusion-projector term Dε_A in Eq. (70), would leave Figs. 2-3 essentially unchanged. The agreement in Figs. 2-3 therefore tests the simplified Thomas-Fermi/OU reduction against SPGPE simulations; it is not a test of the exact SERs. This under-determination is independent of, and compounded by, the acknowledged TIRA/Kohn's-theorem limitation: even if the harmonic trap were a fully legitimate SPGPE setting, the reported comparison would not certify the projector and noise-correlation structure that constitutes the paper's claimed novelty.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives stochastic Ehrenfest relations (SERs) for the stochastic projected Gross-Pitaevskii equation (SPGPE) including both number-damping and energy-damping reservoir terms, retaining all projector corrections that arise from the energy cutoff. The main result is the set of equations (50a)-(50e) together with the noise correlations (46)-(49). The authors then apply this formalism to the centre-of-mass motion of a quasi-1D harmonically trapped condensate, using a Thomas-Fermi ansatz to reduce the SERs, after linearizing in the centre-of-mass displacement and neglecting all projector terms, to an Ornstein-Uhlenbeck system (57)-(58). Analytic steady-state correlation functions and spectra, Eqs. (63a)-(64c), are compared with direct 1D SPGPE simulations in Figs. 2 and 3 and reported to be in close agreement. The paper explicitly acknowledges that the harmonic-trap test system violates Kohn's theorem under the time-independent reservoir approximation (TIRA), while arguing that the formalism remains relevant for non-harmonic and multicomponent systems.","tokens_in":23725,"tokens_out":3402,"duration_ms":35224,"significance":"If the formal derivation is correct, the SERs provide a genuinely useful analytic tool: they extend the earlier number-damping-only Ehrenfest relations of Bradley, Blakie and Gardiner to the full SPGPE with energy damping, and they convert, in special ansatz-based cases, multiplicative field noise into additive noise in collective variables. The derivation is a genuine derivation rather than a fit; the analytic centre-of-mass solution is parameter-free given the simulation parameters, and the comparison against direct SPGPE simulations is an appropriate internal consistency check of the reduction. The manuscript is clear about several limitations, including the TIRA/Kohn-theorem caveat. The central weakness is that the validation exercises only the simplified, projector-free limit of the SERs, so the exact projector and noise-correlation structure that constitutes the paper's headline novelty is not numerically certified.","major_comments":[{"comment":"The reported validation does not actually test the exact SERs that are the paper's central claim. The analytic centre-of-mass model is obtained by imposing the Thomas-Fermi ansatz, linearizing in z(t), and explicitly neglecting all projector terms, as stated before Eq. (55) and again in the passage leading to Eq. (57). Appendix B.3 numerically monitors only two of the many projector contributions, qH_z and dz_epsilon, and argues that the remaining terms are smaller because they carry small damping rates. Consequently an algebraic error in, for example, q_gamma_z, d_gamma_z, or the diffusion-projector term D_epsilon_A in Eq. (70) would leave Figs. 2 and 3 essentially unchanged. The agreement in those figures therefore certifies the Thomas-Fermi/Ornstein-Uhlenbeck reduction against 1D SPGPE simulations, not the exact projector and noise-correlation structure claimed in the abstract and in Sec. 3.5.","section":"Sec. 4.1-4.2 and Appendix B.3"},{"comment":"The chosen test system is one in which the underlying reservoir theory is acknowledged to be suspect: the paper states that the time-independent reservoir approximation violates Kohn's theorem for a harmonically trapped scalar BEC, and that c-field theory is currently best suited to systems with a time-independent high-energy reservoir. Since the 1D SPGPE simulation itself inherits the TIRA, the comparison in Figs. 2 and 3 is an internal consistency check between an approximate analytic reduction and the same approximate equation of motion, not a test of the SERs in a physically validated regime. This limitation is acknowledged in Sec. 5.2, but its consequences for the strength of the validation claim should be stated more prominently, and possible tests in a non-harmonic or two-component system where Kohn's theorem does not apply should be discussed.","section":"Sec. 4, opening paragraphs"},{"comment":"The numerical comparison lacks error bars and a quantitative measure of agreement, yet the abstract and conclusions describe the agreement as 'close' and 'excellent.' The visible deviations at larger tau (Fig. 2) and the absence of statistical uncertainties from the 5000-trajectory ensemble make it difficult to assess whether the residual differences are consistent with projector-term neglect, finite averaging time, or an actual discrepancy in the reduced SERs. Reporting a quantitative fit statistic or confidence interval, especially for the tails of the correlation functions, would materially strengthen the validation claim.","section":"Sec. 4.2, Figs. 2 and 3"}],"minor_comments":[{"comment":"There are numerous typographical errors that should be corrected: 'analtyic' in the Introduction, 'interction' in Sec. 2.2.2, 'linterature' in Sec. 1, 'avarages' in Sec. 3.5(i), 'caonical' in Sec. 2.2.1, 'Futher' in Sec. 2.2, and 'to find the an SDE' in Sec. 3.2.","section":"Throughout"},{"comment":"The sentence describing the Stratonovich correction says 'the first term in the second line' but the term in question appears in the final term of Eq. (31) and in the second line of Eq. (32); please clarify the wording so the reader can identify which term is being named.","section":"Sec. 3.1, Eq. (31)"},{"comment":"The ergodic-averaging description is ambiguous: 'the remaining time interval t = 5t_omega' should specify the starting and ending times of the averaging window explicitly, for example 'from t = 5 t_omega to the end of the trajectory.'","section":"Sec. 4.2, footnote 7"},{"comment":"The denominator |dz/dt| is said to be strictly non-zero for harmonic motion, but for the initial condition x(0)=p(0)=0 used in the simulations the centre-of-mass variable z(0)=0 and its initial rate of change is also zero; the reported relative magnitudes at early times should be interpreted with this in mind.","section":"Appendix B.3, Eq. (83)"},{"comment":"Reference [30] is missing the journal name and volume (it should be Phys. Rev. A 93, 063603), and reference [2] is incomplete as printed. Please check all references for completeness.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely publishable after revision, but I would ask the authors to address the mismatch between the advertised claim ('exact' SERs retaining all projector terms) and the validation, which exercises only a projector-free Ornstein-Uhlenbeck limit. If testing the full projector structure is impractical, the paper should be reframed as a derivation plus a simplified-limit consistency check, with the exact-equation claim confined to the formal section. This is a presentation-of-scope issue more than a correctness issue, but it is load-bearing for the paper's central assertion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers what its title promises: a detailed derivation of stochastic Ehrenfest relations for the SPGPE including energy damping and all projector terms. That derivation is the core contribution, and it is the best part of the paper. The projected functional calculus is handled cleanly, the noise correlations for one-body moments and multiple moments are explicit, and the reduction to an Ornstein-Uhlenbeck model for the centre of mass is a good example of how the formalism can be used. As far as I can tell from reading the algebra, the central equations (50a)-(50e) are correct. The authors have also been honest about the physical limitation of their test system: they clearly state that the time-independent reservoir approximation violates Kohn's theorem for a scalar BEC in a harmonic trap.\n\nThe soft spot is the validation, and it is exactly the soft spot the stress-test identifies. Figures 2 and 3 compare SPGPE simulations with the analytic solution of the Thomas-Fermi, linearized, projector-free OU limit of the SERs. That comparison does not exercise the exact projector terms and noise correlations that constitute the paper's claimed novelty. Appendix B.3 monitors two of the many projector contributions and argues that the others are small because they carry small damping rates, but that is not a direct test of the full equations. If there were an error in an unmonitored projector term, the shown agreement would be unchanged. The lack of error bars on the simulation points is a minor issue by comparison, but the visible degradation of agreement at larger tau is not addressed quantitatively.\n\nThe internal-consistency nature of the test is also worth stating more plainly. The SERs are derived from the same SPGPE used in the simulation, so the comparison is a check that the derived equations reduce correctly to the simulated model, not a test of the SPGPE itself against experiment. That is fine for a formal paper, but it is a different claim than 'experimentally accessible hot BECs'.\n\nNone of this breaks the paper's central argument. The derivation is the product, and it looks solid. The validation section under-sells the main result by testing only a special limit; it would be improved by either a test system where the reservoir model is clean, error bars, or at least a clear statement that the comparison verifies the reduction, not the exact projector terms.\n\nWho is this for? The c-field/SPGPE community. It will be a useful reference for anyone doing analytic work on damped collective modes. It deserves a serious referee and, in my view, publication after revision. The authors should be asked to make the scope of the validation explicit and to quantify the long-time discrepancies.","headline":"A careful and likely correct formal derivation of stochastic Ehrenfest relations for the full SPGPE, with a validation section that only tests a simplified limit of those equations.","tokens_in":24276,"tokens_out":2919,"would_cite":true,"duration_ms":27329,"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":"The paper derives exact stochastic Ehrenfest relations for the full damped SPGPE, giving closed stochastic equations for one-body observables of hot Bose–Einstein condensates.","keywords":["stochastic projected Gross-Pitaevskii equation","Ehrenfest relations","finite-temperature Bose-Einstein condensate","number damping","energy damping","centre-of-mass fluctuations","Ornstein-Uhlenbeck process","c-field theory"],"falsifier":"Compute the projector corrections given in Appendix B for a cutoff chosen so that the population of the highest coherent mode is not small: if those corrections do not stay below a few percent during equilibration, the analytic Ornstein–Uhlenbeck solution is not the theory's prediction. Alternatively, measure the centre-of-mass position spectrum of a hot condensate in a non-harmonic trap and compare its line shape with Eq. (64a); disagreement beyond the computed projector corrections would show the static-reservoir assumption fails.","tokens_in":23234,"feed_emoji":"🧊","tokens_out":7924,"duration_ms":73404,"temperature":0.7,"pith_summary":"This paper derives exact stochastic Ehrenfest relations for the stochastic projected Gross-Pitaevskii equation (SPGPE), the classical-field equation used to simulate finite-temperature Bose-Einstein condensates. The derivation keeps both reservoir channels—number damping and energy damping—and retains every projector term coming from the energy cutoff that separates the coherent low-energy system from the incoherent thermal reservoir. The resulting equations give closed stochastic motion for position, momentum, angular momentum, energy, and particle number, so collective excitations can be studied analytically instead of by simulating the full fluctuating field. As a test, the paper shows that the analytic centre-of-mass correlations and spectra for a quasi-one-dimensional harmonic trap agree closely with direct SPGPE simulations.","feed_headline":"Full damped SPGPE now has exact Ehrenfest relations","feed_subtitle":"Number and energy damping are both kept; analytic centre-of-mass correlations match simulation.","key_machinery":"The load-bearing machinery is the stochastic change-of-variables formula for functionals of the projected c-field, applied after rewriting the SPGPE so that fields and noises are independent at equal times. This converts the multiplicative noise of the field equation into explicit drift, diffusion, and projector-correction terms for any one-body observable; the projector terms describe mode mixing between the highest-energy coherent mode and the lowest-energy incoherent mode, and vanish as the cutoff is raised. For the centre-of-mass test, a Thomas–Fermi wavefunction ansatz evaluates the integrals, leaving a linear stochastic system with constant drift and diffusion matrices—an Ornstein–Uhlenbeck process—whose exact steady-state correlations and spectra are the objects compared with simulation.","core_discovery":"The central claim is that the set of equations (50a)–(50e) are exact stochastic Ehrenfest relations for the full SPGPE: each one-body observable obeys a stochastic differential equation whose drift combines the ordinary Hamiltonian Ehrenfest forces with explicit damping drifts from number and energy exchange with the reservoir, and whose noise has the explicit correlations (46)–(49). All terms generated by the projector—the Hilbert-space cutoff separating system from reservoir—are kept, and the paper shows they are small in the centre-of-mass application. There, a Thomas–Fermi ansatz reduces the position and momentum equations to a two-dimensional Ornstein–Uhlenbeck process, and the steady-state correlation functions and spectra obtained analytically match c-field simulations.","pith_inferences":["The same projection machinery could be used to derive cutoff-corrected equations for other collective modes, such as quadrupole or scissors modes, where the relevant projector overlaps may not be as small as they are for the centre of mass.","If the additive-noise reduction holds more broadly, soliton and vortex Brownian motion could be treated analytically by reading effective damping and diffusion constants off the general noise correlations, without solving the full SPGPE.","Because the paper's static-reservoir approximation is most questionable in exactly harmonic traps, the strongest physical tests of the theory are in non-harmonic, toroidal, or two-species systems, which the paper lists as future targets.","The spectral linewidth formulas (64a)–(64c) offer an experimental route to measuring the reservoir damping rates $\\Lambda_\\gamma$ and $\\Lambda_\\varepsilon$ from equilibrium noise spectra alone, without preparing a non-equilibrium initial state."],"forward_implications":["Equations (50a)–(50e) provide closed stochastic equations for one-body observables of the full SPGPE, so collective-mode dynamics can be studied without simulating the full fluctuating field.","For many one-body operators the multiplicative noise of the SPGPE becomes additive noise in the collective equations, opening these equations to linear-noise and Ornstein–Uhlenbeck techniques.","At equilibrium the SERs yield fluctuation–dissipation relations, such as $\\langle\\langle L-\\mu\\rangle\\rangle=k_B T N$ for pure number damping, that serve as ensemble-level consistency checks.","In the centre-of-mass application, the analytic steady-state correlations and spectra agree with direct SPGPE simulations, validating both the neglect of projector corrections and the overall formalism in that regime."],"supporting_citations":[{"why":"Derives the complete SPGPE with number- and energy-damping reservoir terms and the energy cutoff; it is the object of the paper's SER derivation.","marker":"[10]"},{"why":"Derived the earlier number-damping Ehrenfest relations that this work generalizes by retaining all noises and cutoff terms.","marker":"[21]"},{"why":"Provides the one-dimensional SPGPE and its damping/noise coefficients used in the centre-of-mass test.","marker":"[18]"},{"why":"Presents the SPGPE formulation whose energy-damping interaction and scattering kernel enter the SERs.","marker":"[13]"},{"why":"Introduced the projected Gross-Pitaevskii equation and the cutoff basis that underlies the SPGPE projector.","marker":"[16]"},{"why":"First used the present analytic approach for soliton decay in a toroidal trap, cited as precedent for the method applied here.","marker":"[19]"}],"fun_headline_variants":["Exact stochastic Ehrenfest relations for damped SPGPE","Hot BECs: full damped SPGPE now has exact Ehrenfest equations","Number and energy damping: exact Ehrenfest relations for hot BECs","Analytic Ehrenfest relations match simulations for hot BECs","Full SPGPE with damping now has exact stochastic Ehrenfest relations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the time-independent reservoir approximation—that the high-energy incoherent region can be treated as a static thermal reservoir—which the paper itself notes is not strictly valid for a scalar condensate in a purely harmonic trap, because the exact centre-of-mass mode should be protected by harmonicity.","fun_headline_variants_meta":{"raw":{"variants":["Exact stochastic Ehrenfest relations for damped SPGPE","Hot BECs: full damped SPGPE now has exact Ehrenfest equations","Number and energy damping: exact Ehrenfest relations for hot BECs","Analytic Ehrenfest relations match simulations for hot BECs","Full SPGPE with damping now has exact stochastic Ehrenfest relations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000165,"raw_usage":{"total_tokens":1164,"prompt_tokens":774,"completion_tokens":390,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":390,"completion_tokens_details":{"reasoning_tokens":299}},"tokens_in":390,"tokens_out":390,"duration_ms":3842,"temperature":1.0,"reasoning_tokens":299,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:04:11.465893+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the projector corrections given in Appendix B for a cutoff chosen so that the population of the highest coherent mode is not small: if those corrections do not stay below a few percent during equilibration, the analytic Ornstein–Uhlenbeck solution is not the theory's prediction. Alternatively, measure the centre-of-mass position spectrum of a hot condensate in a non-harmonic trap and compare its line shape with Eq. (64a); disagreement beyond the computed projector corrections would show the static-reservoir assumption fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the complete SPGPE with number- and energy-damping reservoir terms and the energy cutoff; it is the object of the paper's SER derivation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derived the earlier number-damping Ehrenfest relations that this work generalizes by retaining all noises and cutoff terms."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the one-dimensional SPGPE and its damping/noise coefficients used in the centre-of-mass test."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Presents the SPGPE formulation whose energy-damping interaction and scattering kernel enter the SERs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"First used the present analytic approach for soliton decay in a toroidal trap, cited as precedent for the method applied here."}],"review_version":1}