{"id":"3b60a655-9bf6-4c13-a471-45b289bbc67c","arxiv_id":"1908.05464","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"After a sudden stiffness change in a harmonic chain, the kinetic temperature follows a decaying Bessel-function oscillation toward a new equilibrium value.","lead":"This paper derives a formula for how the temperature of a one-dimensional chain of particles oscillates and slowly settles after its stiffness is suddenly changed. The result connects fast mechanical deformation to thermal behavior on a femtosecond timescale.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Linearization at Eq. (14) is unquantified: the Bessel transient (21) is verified only for αε/C=-0.1, and no bound guarantees neglect of αΔε_n²; a parameter sweep is needed.","rationale":"The reader identified the same load-bearing assumption, and I agree; hence agreement_with_reader = agree. My read does not move the verdict: the harmonic derivation is coherent, the virial limit is consistent, and the single numerical case supports Eq. (21) at αε/C=-0.1. The absence of a quantitative small-nonlinearity bound justifies the existing CONDITIONAL verdict rather than rejection. I would not classify the 'unexpected' framing as a correctness issue because the cited [17,18] already reported Bessel oscillations for thermal perturbations; the novelty here is the deformation-initiated setting, not the absence of all prior knowledge. The proposed sweep is a direct test of whether the neglected αΔε_n² term changes the prediction.","tokens_in":7116,"tokens_out":8824,"duration_ms":87573,"concrete_test":"Run a symplectic MD sweep of the full equations (3)–(5), N=5·10^4, 10^3 realizations, same T0 and same ω=sqrt((C+2αε)/m) (adjust C and α accordingly), with αε/C in {-0.5,-0.3,-0.2,-0.1,0.1,0.2,0.3,0.5}. For each run, extract the first-peak amplitude A_1=max|T(t)-T_∞| and the late-time envelope exponent γ from |T(t)-T_∞|~t^{-γ}. Compare A_1 to T0|αε/C|·max|1-J0(4ωt)| and γ to 1/2. If deviations exceed statistical error as |αε/C| grows, Eq. (21) is restricted to a range that the paper must state; if deviations stay within error across all values, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result (21) is obtained by the linearization at Eq. (14)→(15): the term αΔε_n² is discarded, leaving a harmonic chain with renormalized stiffness C+2αε. Everything after Eq. (15)—the covariance equations, the Bessel solution (19), and Eq. (21)—is exact for that linear system, so the paper is internally consistent within its stated model. The load-bearing question is whether the discarded term is actually negligible in the regime the paper claims to treat. The introduction says the nonlinearity must be 'large enough to cause adiabatic heating' yet 'sufficiently small' for harmonic analysis, and Eq. (14) itself says the cubic term is neglected 'in the case of small deformations,' but no quantitative bound on αε/C or on α√(kBT0/C)/C is supplied. The only numerical test, Fig. 1, uses αε/C=-0.1 with 10^3 realizations and no error bars; it demonstrates agreement at one operating point, not a validity domain. Since the same cubic nonlinearity is the mechanism that makes the temperature change (via C+2αε), the prediction's physical reach is exactly as wide as this unquantified linearization. If αε/C is not small, mode coupling through αΔε_n² can alter the oscillation amplitude and the t^{-1/2} envelope of (21). This is a significant limitation, but not a mathematical flaw within the stated harmonic approximation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies an infinite one-dimensional chain with cubic (alpha-FPU) nearest-neighbor interactions, initially at thermal equilibrium at temperature T0. At t=0 the chain is subjected to an instantaneous homogeneous deformation epsilon, which is modeled as a sudden change of the harmonic stiffness from C to C+2*alpha*epsilon. Using covariance analysis of the linearized equations of motion, the authors derive Eq. (21): T(t)=T0+T0*(alpha*epsilon/C)*(1-J0(4*omega*t)). This predicts that the kinetic temperature oscillates with a Bessel-function envelope decaying as t^(-1/2) and approaches the virial-theorem value T0*(1+alpha*epsilon/C) at long times. The analytic result is compared with molecular dynamics for N=5e4 particles at alpha*epsilon/C=-0.1, showing good agreement over several tens of oscillation periods. The paper also gives an estimate for a carbon ring.","tokens_in":7356,"tokens_out":14867,"duration_ms":139761,"significance":"If Eq. (21) is correct, the result is a clean, parameter-free prediction for the thermal transient after an ultrafast mechanical load: the temperature does not relax monotonically but oscillates and decays algebraically. The derivation is internally consistent: no fitting parameters are introduced, the long-time limit reproduces the virial theorem (13), and the numerical comparison at one parameter point is quantitative. The main value is as a benchmark for ultrafast thermomechanics and as a building block for higher-dimensional generalizations. The paper is less novel mathematically, since the Bessel-function solution is taken from [18]; its contribution is the formulation of the mechanical-loading quench and the closed-form result (21).","major_comments":[{"comment":"The linearization that drops alpha*Delta(epsilon_n^2) is not quantitatively justified. The retained nonlinear term 2*alpha*epsilon*Delta(epsilon_n) is first order in the applied strain times the thermal fluctuation, whereas the dropped term is second order in the thermal fluctuation; the validity condition is roughly |epsilon| >> sqrt(k_B*T0/C) (or alpha*sqrt(k_B*T0/C)/C << 1). The paper states only that deformations are small, but the same coupling alpha is responsible for the stiffness renormalization C+2*alpha*epsilon, so the domain of validity of Eq. (21) is exactly the unquantified part of the approximation. Please state the small parameter(s), estimate the size of the neglected term, and supplement Fig. 1 with a sweep over alpha*epsilon/C (and over the thermal fluctuation amplitude) to show where the Bessel prediction breaks down.","section":"Section IV, Eq. (14)"},{"comment":"The derivation of the initial-value problem for L_k is compressed into a single sentence, and the solution (19) is imported from [18]. Because Eq. (21) is the central claim, the covariance calculation should be given in an appendix: differentiating (17) with the linearized equations (15), evaluating the initial values L_k(0) and dL_k/dt(0), and solving the resulting discrete wave equation by Fourier transform to obtain J_{2k}(4*omega*t). This would also make the paper self-contained and allow the reader to check the factor 4 in the argument and the sign of the initial condition.","section":"Section IV, Eqs. (18)-(19)"},{"comment":"The notation alpha*Delta(epsilon_n^2) is ambiguous. If Delta(epsilon_n) = epsilon_{n+1}-epsilon_n, the nonlinear force difference is alpha*(epsilon_{n+1}^2-epsilon_n^2) = alpha*(epsilon_{n+1}-epsilon_n)*(epsilon_{n+1}+epsilon_n), not alpha*(epsilon_{n+1}-epsilon_n)^2. If instead alpha*Delta(epsilon_n^2) is intended, this should be written explicitly. As it stands, a reader following the displayed equation would conclude there is an algebraic error in the equation of motion.","section":"Section IV, Eq. (14)"}],"minor_comments":[{"comment":"The sentence \"Differentiating of (17) and equations of motion (14) lead to...\" should refer to the linearized equations (15), not the full nonlinear equations (14).","section":"Section IV"},{"comment":"The units of epsilon are not stated: in Eq. (2) and (16), epsilon is a displacement difference with units of length because C is a force constant, while in the text it is called a deformation, which suggests a dimensionless quantity. The same applies to the example in Section V, where epsilon=0.01*a. Please state explicitly that epsilon is a length, and give the units of alpha.","section":"Section II and Section V"},{"comment":"The term \"adiabatic transition\" is potentially misleading; the protocol is an instantaneous quench, not a slow adiabatic change. \"Thermally isolated\" or \"sudden parameter change\" would be clearer.","section":"Abstract and Introduction"},{"comment":"The figure does not include error bars or a discussion of numerical accuracy. A short statement on the dependence of the result on the integration step, the number of particles, and the number of realizations would strengthen the comparison.","section":"Fig. 1"},{"comment":"The initial conditions specify variances but not the distribution of rho_n and varrho_n; for the linear system only second moments matter, but stating that the variables are Gaussian (or that the result is independent of the distribution) would be helpful.","section":"Section IV, Eq. (16)"},{"comment":"The estimate should note that alpha from Eq. (23) is an anharmonic parameter whose accuracy limits the predicted temperature change; a one-sentence discussion of uncertainty would be appropriate.","section":"Section V"}],"recommendation":"major_revision","confidential_remarks":"The mathematical content is largely an application of the covariance method of [18] to a new physical setup and initial condition; the authors should make the incremental contribution clearer, possibly by citing [22] for the derivation. The paper fits the scope of the journal. The revision should focus on the validity of the linearization and on a self-contained derivation of the central result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nThis is a short paper with one clean result: after an instantaneous homogeneous deformation of an infinite 1D chain with a small cubic term, the kinetic temperature evolves as T0 + T0(αε/C)(1 − J0(4ωt)), then settles to the virial value T0(1 + αε/C). That formula is new in application—earlier work from the same school dealt with instantaneous heating, not mechanical loading—and the derivation inside the linearized model is standard and, as far as I can tell, correct. No fitting, no free parameters, and the limiting cases check out: Eq. (21) reduces to Eq. (13) as t→∞ and starts at T0 at t=0. The one MD comparison in Fig. 1 covers many oscillation periods and looks quantitatively good.\n\nThe soft spot is exactly what the stress-test says. Eq. (14) drops αΔε_n² \"for small deformations\" without a criterion. Since the same cubic term is what shifts the stiffness from C to C+2αε, the physical reach of the result is only as wide as that linearization, and the paper does not bound αε/C, α√(kBT0/C)/C, or the thermal fluctuation amplitude. Fig. 1 uses one value, αε/C = −0.1, with no parameter sweep and no error bars. That is not a mathematical flaw—within the linear chain everything follows cleanly—but it is a real gap between the claimed \"transition to thermal equilibrium\" and what is actually demonstrated. I would also call the \"usually unexpected\" statement in the abstract overstated: Klein–Prigogine and the group's own earlier papers already show oscillatory energy after a quench; what is new here is the mechanical trigger, not the ringing itself.\n\nMinor points: Eq. (18) is sketched rather than derived, and the Bessel solution is imported from [18]. That is legitimate—[18] is a published, parameter-free solution—but I would ask the authors to show enough of the covariance step to make the paper self-contained. Section V's carbon estimate is a nice sanity check, though it is just one number.\n\nVerdict: useful analytic benchmark for ultrafast thermomechanics and for people testing MD codes against harmonic transients. It deserves serious peer review; the referee request should be a quantitative smallness criterion and at least one more parameter value. It does not deserve a desk reject.","headline":"A small, honest harmonic-chain calculation with one genuinely new formula and one unquantified linearization; worth refereeing, not a breakthrough.","tokens_in":7925,"tokens_out":4566,"would_cite":true,"duration_ms":45538,"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":"This paper shows that an instantaneous homogeneous deformation of an anharmonic crystal does not heat it monotonically: its kinetic temperature oscillates with a Bessel-function envelope and relaxes as the inverse square root of time to…","keywords":["kinetic temperature","thermal equilibration","anharmonic chain","Bessel function","instantaneous deformation","virial theorem","femtosecond dynamics","generalized Lagrangian"],"falsifier":"Integrate the full equations of motion (14) without dropping $\\alpha\\Delta\\varepsilon_n^2$, for example with $\\alpha\\varepsilon/C = -0.3$, or run the existing $-0.1$ case for many more oscillation periods. If the temperature envelope departs from the $t^{-1/2}$ Bessel decay while the oscillation amplitude is still visible, the closed formula (21) holds only in a narrower range than the paper's small-nonlinearity premise states.","tokens_in":6863,"feed_emoji":"🌡️","tokens_out":6005,"duration_ms":56476,"temperature":0.7,"pith_summary":"This paper claims that when an infinite one-dimensional crystal is deformed instantaneously, its temperature does not relax smoothly toward equilibrium: it oscillates, with an amplitude that decays only as the inverse square root of time. The setting is a chain of point masses with a small cubic anharmonicity, so the sudden deformation acts as an instantaneous change of the effective bond stiffness. From the equations of motion and a covariance analysis the authors derive an explicit time-dependent formula for kinetic temperature, $T = T_0 + T_0\\frac{\\alpha\\varepsilon}{C}\\bigl(1 - J_0(4\\omega t)\\bigr)$, and confirm it by numerical simulation. This matters because ultrafast mechanical loads on solids are now realizable on femtosecond timescales, where thermal and mechanical processes happen on the same clock.","feed_headline":"A sudden squeeze sends crystal temperature into long-lived oscillations","feed_subtitle":"One-dimensional model predicts ringing that fades as the inverse square root of time.","key_machinery":"The load-bearing object is the generalized Lagrangian $L_k = K_k - U_k$, the difference between the generalized kinetic and potential energies of particle pairs separated by $k$ lattice spacings. Its equations of motion form a discrete wave equation, and the instantaneous load enters only through the initial condition $L_k = -T_0 k_B \\frac{\\alpha\\varepsilon}{C}\\delta_k$. The solution is $L_k = -T_0 k_B \\frac{\\alpha\\varepsilon}{C} J_{2k}(4\\omega t)$, and the $k=0$ component $L_0$ is exactly the difference between kinetic temperature and half the total energy. Bessel functions are the natural Green's functions of this discrete wave equation, so this object is what turns the mechanical shock into ringing, decaying temperature oscillations.","core_discovery":"The central discovery is that an adiabatic transition between two equilibrium states with different stiffness is oscillatory rather than monotonic. After a homogeneous strain $\\varepsilon$ is suddenly applied, the kinetic temperature is $T(t)=T_0 + T_0\\frac{\\alpha\\varepsilon}{C}\\bigl(1-J_0(4\\omega t)\\bigr)$, where $T_0$ is the initial temperature, $\\alpha$ and $C$ are the cubic and quadratic stiffness coefficients, and $\\omega=\\sqrt{(C+2\\alpha\\varepsilon)/m}$ is the modified phonon frequency. The asymptotic form of the Bessel function makes the oscillation amplitude decay as $t^{-1/2}$, so the transient is long-lived. At long times the formula reduces to $T_0(1+\\alpha\\varepsilon/C)$, exactly the value obtained independently from the virial theorem in Eq. (13). Numerical simulations of a chain of $5\\cdot10^4$ particles, averaged over $10^3$ realizations, match the analytical curve over tens of oscillation periods.","pith_inferences":["If Eq. (21) is right, a pump-probe experiment on a femtosecond-shocked film should see lattice-temperature ringing at the phonon timescale; observing faster damping would indicate that nonlinear mode coupling matters even at small deformation.","The oscillation frequency $4\\omega$ depends on $C+2\\alpha\\varepsilon$, so the transient itself offers a way to measure the cubic stiffness coefficient from a non-equilibrium measurement, independent of the equation of state.","The same covariance method, applied in two or three dimensions, will likely replace the Bessel function by a product or integral of Bessel functions; the qualitative $t^{-d/2}$ amplitude decay is a testable prediction.","Because the linearization discards $\\alpha\\Delta\\varepsilon_n^2$, the formula cannot capture energy exchange between phonon modes; at larger strains the decay should cross over from algebraic to exponential, and locating that crossover would sharpen the model's validity domain."],"forward_implications":["A fast mechanical load produces a measurable oscillatory thermal transient, not a monotonic approach to the new equilibrium temperature.","The transient decays only as $t^{-1/2}$, so appreciable temperature oscillations can persist for many vibrational periods after an ultrashort load.","At long times the kinetic and potential energies equalize and the temperature reaches the virial-theorem value $T_0(1+\\alpha\\varepsilon/C)$.","For a carbon ring with 1% bond deformation the model gives a temperature change of about 0.6% and an oscillation period near 10 femtoseconds.","The covariance approach extends to two- and three-dimensional crystals, so the same Bessel-type ringing should appear in more realistic lattices."],"supporting_citations":[{"why":"Supplies the virial relation $2K = \\langle\\varepsilon_n F_n\\rangle$ used to write the initial kinetic temperature before and after loading.","marker":"[8]"},{"why":"Provides the pioneering direct solution showing that energy oscillations after an instantaneous thermal perturbation are described by a Bessel function.","marker":"[17]"},{"why":"Gives the solution of the analogous initial-value problem for the generalized Lagrangian, which the paper adapts to obtain Eq. (19).","marker":"[18]"},{"why":"Justifies that thermal processes in crystals with sufficiently low nonlinearity are described accurately by the harmonic approximation on the timescales considered.","marker":"[20]"},{"why":"Provides the covariance-analysis approach used to derive the generalized Lagrangian equations and initial conditions.","marker":"[22]"},{"why":"Supplies the Bessel-function asymptotic expansion that yields the $t^{-1/2}$ decay of the temperature oscillation amplitude.","marker":"[35]"}],"fun_headline_variants":["Crystal temperature rings down after sudden squeeze","Sudden stiffness jump triggers long-lived temperature oscillations","Temperature ringing: crystal equilibration is non-monotonic","Crystal temperature rings like a bell after sudden strain"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes the cubic nonlinearity is small enough that the equation of motion can be linearized, keeping the deformation only through the modified stiffness $C+2\\alpha\\varepsilon$ and dropping the quadratic deformation term $\\alpha\\Delta\\varepsilon_n^2$.","fun_headline_variants_meta":{"raw":{"variants":["Crystal temperature rings down after sudden squeeze","Sudden stiffness jump triggers long-lived temperature oscillations","Temperature ringing: crystal equilibration is non-monotonic","Crystal temperature rings like a bell after sudden strain"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000561,"raw_usage":{"total_tokens":2606,"prompt_tokens":829,"completion_tokens":1777,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":445,"completion_tokens_details":{"reasoning_tokens":1716}},"tokens_in":445,"tokens_out":1777,"duration_ms":15028,"temperature":1.0,"reasoning_tokens":1716,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:13:53.429099+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Integrate the full equations of motion (14) without dropping $\\alpha\\Delta\\varepsilon_n^2$, for example with $\\alpha\\varepsilon/C = -0.3$, or run the existing $-0.1$ case for many more oscillation periods. If the temperature envelope departs from the $t^{-1/2}$ Bessel decay while the oscillation amplitude is still visible, the closed formula (21) holds only in a narrower range than the paper's small-nonlinearity premise states.","supporting_citations":[{"cited_title":"A., Krivtsov, A","cited_arxiv_id":null,"evidence_quote":"Supplies the virial relation $2K = \\langle\\varepsilon_n F_n\\rangle$ used to write the initial kinetic temperature before and after loading."},{"cited_title":"Physica 19 (1953)","cited_arxiv_id":null,"evidence_quote":"Provides the pioneering direct solution showing that energy oscillations after an instantaneous thermal perturbation are described by a Bessel function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the solution of the analogous initial-value problem for the generalized Lagrangian, which the paper adapts to obtain Eq. (19)."},{"cited_title":"A., Krivtsov, A","cited_arxiv_id":null,"evidence_quote":"Justifies that thermal processes in crystals with sufficiently low nonlinearity are described accurately by the harmonic approximation on the timescales considered."},{"cited_title":"N., Krivtsov, A","cited_arxiv_id":null,"evidence_quote":"Provides the covariance-analysis approach used to derive the generalized Lagrangian equations and initial conditions."},{"cited_title":"and Stegun, I","cited_arxiv_id":null,"evidence_quote":"Supplies the Bessel-function asymptotic expansion that yields the $t^{-1/2}$ decay of the temperature oscillation amplitude."}],"review_version":1}