{"id":"4262e804-5f1b-4741-9747-1b933d3715f9","arxiv_id":"1908.09225","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"No on-site finite-difference velocity can give time-step-independent kinetic temperature for the GJF scheme, while an infinite one-parameter family of two-point velocities does.","lead":"This paper shows that the GJF thermostat, which samples positions correctly, cannot have an on-site velocity with correct time-step-independent kinetic temperature, but does have an infinite family of two-point (leap-frog) velocity definitions that give correct kinetic energy for linear systems. The result gives simulation practitioners a one-parameter menu of velocities for measuring temperature and diffusion in Langevin molecular dynamics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The no-go for on-site velocities in §II.B is narrower than the abstract: it holds only if an on-site velocity must reduce to central difference as αdt→0.","rationale":"The paper does what it sets out to do for two-point velocities: the derivation in §II.C is explicit, the one-parameter family follows from Eqs. (31)–(33), and the Lennard-Jones data support the kinetic-temperature behavior of Cases A, B, and C. My only substantive objection targets the first advertised conclusion. The no-on-site proof in §II.B is not a no-go for all on-site finite-difference velocities; it is a no-go for on-site velocities whose zero-friction limit is the central difference. That limit is imposed by hand as part of the definition of 'reasonable on-site velocity'. Relaxing it to first-order consistency admits explicit finite-coefficient solutions to the same mean-square condition, so the abstract's 'not possible' is an overstatement. Because the overstatement is confined to the no-go headline and does not affect the constructive family, the reader's CONDITIONAL verdict stands.","tokens_in":13367,"tokens_out":26755,"duration_ms":248873,"concrete_test":"Independently re-derive Eq. (27) from Eq. (22) (being careful with the cross-term factors) and substitute the backward stencil γ1=0, γ2=1, γ3=−1, γ5=0 with γ4=(−b+√(b^2+b))/2. If the expression equals 1 for b∈(0,1), the accepted no-go statement must be narrowed to the central-difference limiting definition; if it does not, report the correct cross-term factor and re-test with the root of the corrected quadratic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The main construction—the one-parameter two-point family of Section II.C—is internally consistent and simulation-supported. The load-bearing weakness is the first conclusion advertised in the abstract: that no on-site velocity can give time-step-independent kinetic temperature. Section II.B proves this only under the added requirement that a 'reasonable' on-site velocity reduce to the central-difference stencil γ1=1/2, γ3=−1/2 in the limit αdt→0. That requirement is what makes γ4=bγ3/[4(1−b)] diverge. If the central-difference limit is relaxed to the weaker first-order condition γ1−γ3→1 (with γ3→0), Eq. (27) admits solutions with finite coefficients. For instance, on the backward-looking three-point stencil γ1=0, γ2=1, γ3=−1, γ5=0, Eq. (27) reduces to b+4(1−b)γ4+4(1−b)γ4^2/b=1, which has the finite root γ4=(−b+√(b^2+b))/2 for every b∈(0,1). This velocity is a time-shifted member of the paper's own two-point family, so whether the abstract's 'impossible' survives is a matter of how 'on-site' is defined. The Discussion does state the central-difference qualifier, but the abstract and the opening of Section IV do not. The two-point family itself, and the simulation comparisons of Cases A, B, and C, are not threatened by this concern.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies possible velocity definitions that can accompany the Grønbech-Jensen-Farago (GJF) thermostat. Starting from the GJF configurational trajectory, the authors write a general finite-difference velocity involving three consecutive positions and the two adjacent noise increments [Eq. (18)], impose the requirement that the mean-square velocity equal the equipartition value k_B T/m for every stable time step in a harmonic oscillator, and solve the resulting constraints. They conclude that no on-site velocity in this class can have time-step-independent kinetic temperature, while a one-parameter family of two-point velocities can; they give the explicit parametrization [Eqs. (31)-(34)], express the resulting leap-frog algorithms [Eqs. (35)-(41)], identify Cases A, B, and C, and test them on Lennard-Jones solids and liquids. The central derivation for the two-point family is explicit and internally consistent, and the simulations corroborate the predicted time-step independence for the highlighted choices.","tokens_in":13669,"tokens_out":18210,"duration_ms":174923,"significance":"The central constructive result is a useful unification: one parameter γ1 generates infinitely many two-point velocities with correct, time-step-independent mean kinetic energy for linear systems, including the previously known 2GJ half-step velocity (Case A) and the Farago velocity (Case B). The analytic derivation is not fitted to data, the algorithmic forms in Eqs. (35)-(41) are directly implementable, and the paper is careful to note that correct kinetic energy alone is not sufficient (Case D) and to give Green-Kubo consistency criteria. If the construction holds, the paper provides practical guidance for choosing a velocity estimator in GJF-based Langevin simulations. The main weakness is the advertised no-go statement for on-site velocities, which is narrower than the abstract suggests and needs qualification.","major_comments":[{"comment":"The no-go result for on-site velocities is proved only under the convention, stated in the text after Eq. (29), that a 'reasonable' on-site velocity must approach the central-difference stencil γ1=1/2, γ3=−1/2 as αdt→0. The abstract and the first sentence of Section IV state the impossibility without this qualifier. The convention is load-bearing: with only first-order consistency required, the factor bγ1 in Eq. (29) can vanish. For example, γ1=0, γ2=1, γ3=−1, γ5=0 makes Eq. (27) reduce to b + 4(1−b)γ4^2/b + 4(1−b)γ4 = 1, which has the finite root γ4=(−b+√(b^2+b))/2 for every b∈(0,1). This is a backward two-point velocity, a time-shifted member of the family used for the paper's positive result, with finite coefficients and no diverging noise term. The conclusion should therefore be reworded as a conditional statement, or proved without the central-difference limiting assumption.","section":"Abstract and §II.B (after Eq. (29))"},{"comment":"The claim of a 'complete investigation of all possible finite difference approximations' exceeds what is shown. The analysis covers only the three-position, two-noise ansatz of Eq. (18); velocities built from more than three consecutive positions or from additional noise increments (for example β_{n−1}) are not treated. Both main conclusions should be stated as applying within the family Eq. (18), unless a supplementary argument establishes that no wider ansatz can work.","section":"Abstract and §II.C (Eq. (18))"}],"minor_comments":[{"comment":"The notation for the square root is ambiguous; the expression should be written explicitly as sqrt{b(1−b^2 γ1^2)/(1−b)} so that the restriction (bγ1)^2≤1 follows immediately.","section":"Eq. (33)"},{"comment":"Refs. [27]-[29] are cited with DOIs and arXiv identifiers instead of complete bibliographic entries; please complete them in the journal's reference format.","section":"References"},{"comment":"The text contains typographical artifacts, for example an extra space in 'G JF' in the title line and an awkward line break inside 'kinetic averages' in Section I; a careful proofread would improve readability.","section":"Global"},{"comment":"The sentence describing the Lennard-Jones kinetic results as 'as expected from the analysis above' extrapolates an exact harmonic-oscillator result to a nonlinear system; the agreement is numerical validation rather than a consequence of the analytic derivation, and the text should say so explicitly.","section":"Section III"}],"recommendation":"major_revision","confidential_remarks":"The two-point family and its algorithmic implementation are solid, and the simulation evidence supports the main constructive claim. The revision should focus on the overreaching no-go statement for on-site velocities and on the 'complete investigation' wording; with those qualified, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper is worth reading if you care about defining velocities in the GJF thermostat. The main result is a one-parameter family of two-point velocities that give time-step-independent mean kinetic energy for linear systems, with explicit formulas and a leap-frog algorithm. It subsumes the two previously known velocities (Jensen-Grønbech-Jensen and Farago) as special cases, so it really is a classification rather than just another example. The derivation for linear systems is straightforward and the key condition, Eq. (32), is consistent with the special cases. The Lennard-Jones simulations show the family holds up in nonlinear systems, and configurational statistics are unaffected, as expected.\n\nThe soft spots are all about overstatement. The abstract says 'all possible finite difference approximations' and concludes no on-site velocity can have correct, time-step-independent kinetic temperature. The actual proof in Sec. II.B assumes the on-site velocity must reduce to the central-difference stencil (γ1→1/2, γ3→−1/2) as αdt→0. That is a reasonable physical requirement, but it is an assumption. The stress-test note gives a concrete backward-looking stencil with γ3→−1 that satisfies Eq. (27) with finite coefficients; it is a time-shifted member of the paper's own two-point family. So the no-go holds only for on-site definitions that approach central difference in the zero-friction limit. The Discussion does state this qualifier, but the abstract and the opening of Sec. IV do not. Similarly, 'all possible' really means 'all within the three-point/two-noise ansatz of Eq. (18).' Both are easy to fix with rewording.\n\nThe extension to nonlinear systems is empirical rather than proven—the LJ simulations are convincing but not a proof for general potentials. No code is shipped, which is a minor inconvenience but not a serious flaw. The citation pattern is mostly self-referential, but the earlier results are independently established, so this is not a red flag.\n\nOverall, the central construction is solid and the paper is honest about its assumptions once you read past the abstract. It is a within-subfield methodological contribution; its audience is molecular dynamics practitioners using Langevin thermostats. I would send it to peer review with a request to correct the overstatements in the abstract and conclusion. It should be conditionally accepted.\n\nBest","headline":"A solid extension of the GJF kinetic-velocity program; the one-parameter family is the real result, and the no-on-site claim needs a qualifier.","tokens_in":14183,"tokens_out":3389,"would_cite":true,"duration_ms":27197,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.10.-a","02.70.Ns"],"model":"deepseek-v4-flash","headline":"For the GJF thermostat, no on-site velocity can give time-step-independent kinetic temperature, but a one-parameter family of two-point velocities can.","keywords":["GJF thermostat","Langevin dynamics","discrete-time velocity","kinetic temperature","leap-frog velocity","finite-difference approximation","equipartition","molecular dynamics"],"falsifier":"Take the on-site ansatz $w=(r_{n+1}-r_{n-1})/(2dt)+\\lambda\\beta_n/m$ with a finite, fixed $\\lambda$, and insert the GJF trajectory correlations into the mean-square velocity; if any finite $\\lambda$ makes $\\langle w^2\\rangle=k_B T/m$ for all stable $\\Omega_0 dt$, the no-go claim for on-site velocities is wrong. Conversely, testing any two-point member of the family on a strongly anharmonic oscillator and seeing a systematic time-step trend in $\\langle w^2\\rangle$ larger than statistical error would bound how far the linear-system proof extends.","tokens_in":13145,"feed_emoji":"⚛️","tokens_out":6526,"duration_ms":57007,"temperature":0.7,"pith_summary":"The paper asks whether discrete-time Langevin simulations can measure kinetic temperature as reliably as they measure positions. Working from the GJF thermostat, an algorithm already known to produce time-step-independent configurational statistics, it examines every finite-difference approximation to velocity. It concludes that no on-site (same-time) velocity can give correct, time-step-independent kinetic energy, but that infinitely many two-point (half-step) velocities can, parameterized by one free constant. Concrete members of the family are written into practical leap-frog algorithms. If correct, the result tells molecular-dynamics users exactly which velocity formulas to use so that kinetic and configurational thermometers agree.","feed_headline":"Infinite velocity choices give exact GJF kinetic energy","feed_subtitle":"A one-parameter family of leap-frog velocities makes kinetic and configurational thermometers agree for GJF.","key_machinery":"The engine of the argument is the general finite-difference velocity ansatz, Eq. (18), with five unitless coefficients ($\\gamma_1$ through $\\gamma_5$) spanning all three-point position differences plus the two noise increments that bracket a time step. Plugging the GJF trajectory correlations into the mean-square velocity produces Eq. (27), whose requirement of time-step independence imposes algebraic constraints. For the two-point subclass the condition collapses to a single equation whose solution parameterizes the whole family by $\\gamma_1$. The analysis of on-site velocities turns on the limiting behavior of $\\gamma_3$ as $\\alpha dt\\to 0$, where the divergence in the noise coefficient $\\gamma_4$ eliminates the possibility.","core_discovery":"The central claim is that kinetic statistics can be made exactly consistent with the GJF trajectory by choosing the right velocity definition. For the linear harmonic oscillator $f=-\\kappa r$, the paper derives the general three-point finite-difference velocity $w=(\\gamma_1 r_{n+1}+\\gamma_2 r_n+\\gamma_3 r_{n-1})/dt+(\\gamma_4 \\beta_n+\\gamma_5 \\beta_{n+1})/m$ and imposes that $\\langle w w\\rangle=k_B T/m$ for every stable time step. This forces the coefficients of $(\\Omega_0 dt)^{-2}$ and $(\\Omega_0 dt)^2$ in the mean-square velocity to vanish. The inverse-squared term rules out on-site velocities: in the frictionless limit any on-site velocity must approach the central difference $\\gamma_1=1/2$, $\\gamma_3=-1/2$, which drives the required noise coefficient to infinity. For two-point velocities ($\\gamma_3=\\gamma_4=0$), the remaining condition is $b\\gamma_1^2+4(1-b)\\gamma_5^2/b+4(1-b)\\gamma_1\\gamma_5=1$, yielding a one-parameter family $\\gamma_5(\\gamma_1)$ with $|b\\gamma_1|\\le 1$. Three highlighted members are Case A ($\\gamma_1=1/\\sqrt{b}$, $\\gamma_5=0$, the 2GJ velocity), Case B ($b\\gamma_1=1$, $\\gamma_5=-1/2$), and Case C ($\\gamma_1=1$). The paper also provides a leap-frog algorithm update that realizes any member and notes that configurational sampling remains unchanged.","pith_inferences":["The one-parameter freedom is likely a resource for tuning higher-order statistics: members of the family share the same mean kinetic energy but differ in velocity autocorrelation and fluctuations, so future work could pick $\\gamma_1$ to minimize variance of kinetic temperature estimates.","A similar no-go argument may apply to any integrator whose natural velocity is on-site and must match the central-difference limit, not just GJF; the divergence mechanism is generic.","The linear-system result rigorously covers harmonic potentials, but the paper's Lennard-Jones tests suggest the family may work for nonlinear systems too; a systematic nonlinear proof or a counterexample remains open.","For non-equilibrium simulations with a drift velocity, Case C's property of giving the correct ballistic velocity may make it preferable even though it is not a true half-step velocity; this distinction is a practically testable prediction."],"forward_implications":["Users of the GJF thermostat can now choose a velocity definition matched to their measurement: Case A for equilibrium diffusive transport, Case C for drift or ballistic motion, and Case B as a maximal-amplitude option.","All members of the family give exactly the same configurational sampling as GJF, since they are built on the same position trajectory; kinetic and configurational temperatures agree across the full stability range for linear systems.","The family subsumes the two previously reported kinetically correct velocities, showing they are special cases of one continuous parameter.","Lennard-Jones simulations across two thermodynamic states and three friction values confirm that the highlighted definitions measure kinetic temperature with time-step-independent averages, while the on-site GJF velocity departs increasingly with time step.","Green-Kubo diffusion evaluation changes with the chosen velocity: Case A requires a right-Riemann sum and Case B a trapezoidal sum to reproduce the Einstein diffusion coefficient."],"supporting_citations":[{"why":"Introduces the GJF thermostat whose position trajectory all velocity definitions in this paper are built on.","marker":"[15]"},{"why":"Identifies the first statistically correct half-step velocity (2GJ), which becomes Case A of the family.","marker":"[19]"},{"why":"Provides a related GJF velocity formulation that becomes Case B of the family.","marker":"[26]"},{"why":"Establishes the configurational sampling robustness of the GJF method that the new velocities must preserve.","marker":"[24]"},{"why":"Defines the standard Verlet central-difference and leap-frog velocities that the paper's ansatz generalizes and must reduce to in the zero-friction limit.","marker":"[4]"},{"why":"Supports the claim that Case A is the true half-step velocity correctly cross-correlated with the spatial degree of freedom.","marker":"[29]"}],"fun_headline_variants":["One-parameter velocity family fixes GJF kinetic stats","GJF kinetics exact via leap-frog, not on-site velocities","Infinite GJF velocities make kinetic energy exact","GJF kinetic statistics solved by leap-frog velocity choice","No on-site velocity can preserve GJF kinetic temperature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results assume that a kinetically correct velocity must match the standard central-difference or leap-frog form when friction vanishes, and that matching equipartition for the harmonic oscillator at every stable time step is the defining test of correct kinetic statistics.","fun_headline_variants_meta":{"raw":{"variants":["One-parameter velocity family fixes GJF kinetic stats","GJF kinetics exact via leap-frog, not on-site velocities","Infinite GJF velocities make kinetic energy exact","GJF kinetic statistics solved by leap-frog velocity choice","No on-site velocity can preserve GJF kinetic temperature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00026,"raw_usage":{"total_tokens":1661,"prompt_tokens":1090,"completion_tokens":571,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":706,"completion_tokens_details":{"reasoning_tokens":492}},"tokens_in":706,"tokens_out":571,"duration_ms":5572,"temperature":1.0,"reasoning_tokens":492,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:19:08.564611+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the on-site ansatz $w=(r_{n+1}-r_{n-1})/(2dt)+\\lambda\\beta_n/m$ with a finite, fixed $\\lambda$, and insert the GJF trajectory correlations into the mean-square velocity; if any finite $\\lambda$ makes $\\langle w^2\\rangle=k_B T/m$ for all stable $\\Omega_0 dt$, the no-go claim for on-site velocities is wrong. Conversely, testing any two-point member of the family on a strongly anharmonic oscillator and seeing a systematic time-step trend in $\\langle w^2\\rangle$ larger than statistical error would bound how far the linear-system proof extends.","supporting_citations":[{"cited_title":"Vanden-Eijnden, G","cited_arxiv_id":null,"evidence_quote":"Introduces the GJF thermostat whose position trajectory all velocity definitions in this paper are built on."},{"cited_title":"Paquet, H","cited_arxiv_id":null,"evidence_quote":"Identifies the first statistically correct half-step velocity (2GJ), which becomes Case A of the family."},{"cited_title":"Jackson, J","cited_arxiv_id":null,"evidence_quote":"Provides a related GJF velocity formulation that becomes Case B of the family."},{"cited_title":"Rickayzena and J","cited_arxiv_id":null,"evidence_quote":"Establishes the configurational sampling robustness of the GJF method that the new velocities must preserve."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the standard Verlet central-difference and leap-frog velocities that the paper's ansatz generalizes and must reduce to in the zero-friction limit."},{"cited_title":"Farago, Physica A 534, 122210 (2019)","cited_arxiv_id":null,"evidence_quote":"Supports the claim that Case A is the true half-step velocity correctly cross-correlated with the spatial degree of freedom."}],"review_version":1}