Transition to thermal equilibrium in a deformed crystal

Andrey S. Murachev; Anton M. Krivtsov

REVIEW 3 major objections 6 minor 39 references

Transition to thermal equilibrium in a deformed crystal

T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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…

desk verdict A small, honest harmonic-chain calculation with one genuinely new formula and one unquantified linearization; worth refereeing, not a breakthrough. read the letter →

arxiv 1908.05464 v2 pith:FFAESECY submitted 2019-08-15 cond-mat.stat-mech

Anton M. Krivtsov , Andrey S. Murachev This is my paper · ORCID

classification cond-mat.stat-mech

keywords kinetictemperaturethermalequilibrationanharmonicchainBesselfunctioninstantaneousdeformationvirialtheoremfemtoseconddynamicsgeneralizedLagrangian

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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$.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague—

This 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.

The 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.

Minor 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.

Verdict: 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.

Referee Report

3 major / 6 minor

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.

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 (3)

[Section IV, Eq. (14)] 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 IV, Eqs. (18)-(19)] 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 IV, Eq. (14)] 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.

minor comments (6)

[Section IV] 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 II and Section V] 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.
[Abstract and Introduction] 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.
[Fig. 1] 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 IV, Eq. (16)] 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 V] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the transient temperature formula follows from stated initial conditions and a parameter-free Bessel solution, with no fitted inputs used as predictions.

full rationale

The central result, Eq. (21), is obtained by linearizing the equations of motion to a harmonic chain with renormalized stiffness C+2αε (Eq. (15)), prescribing independent random initial velocities and displacements through Eq. (16), and then solving the covariance initial-value problem (18). The solution (19) is imported from the authors' earlier work [18], but that is a parameter-free published mathematical result for a similar linear problem, not a fit calibrated to the present data; under the stated harmonic approximation the subsequent algebra from (19)-(20) to (21) is a direct derivation. The equilibrium limit (13) is likewise derived independently from the virial relation (A4) and the initial variances, and the long-time limit of (21) coincides with it automatically because J0(4ωt)→0. Numerical simulation with αε/C=-0.1 confirms the analytical expression, providing an external check. The main limitation, namely the unquantified neglect of αΔε_n^2 at Eq. (14), concerns the validity domain of the harmonic approximation rather than any circularity between the inputs and the claimed prediction. No fitted parameter is renamed as a prediction, no result is defined in terms of the quantity it claims to derive, and the self-citations used are not load-bearing in a circular sense.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central derivation contains no fitted free parameters: C, α, ε, T0, and m are physical inputs, and the temperature formula (21) is an exact consequence of the linearized equations and the imported Bessel solution. The main burden is the small-nonlinearity linearization and the equilibrium initial conditions. No new entities are introduced.

assumptions (4)

standard math The covariance initial-value problem (18) has solution L_k = -(T0 kB αε/C) J_{2k}(4ωt), imported from [18].
Eq. (19) is cited to the earlier paper; the present paper does not rederive it, so the Bessel solution is an external mathematical input.
domain assumption Nonlinearity is small enough that αΔε_n² is negligible in Eq. (14), keeping the dynamics linear except for the effective stiffness C + 2αε.
Section IV; this linearization is load-bearing because the Bessel solution describes a linear system only.
domain assumption Before loading the chain is in thermodynamic equilibrium with kBT0 = C<ε_n²> and <v_n²> = kBT0/m, with velocities and deformations independent, uncorrelated, zero-mean random quantities.
Eqs. (9) and (16); these initial conditions determine L0 = -T0 kB αε/C in Eq. (18).
domain assumption The chain is infinite (or periodic in numerics) with nearest-neighbor interactions, so translational invariance holds and covariance functions depend only on the index difference.
Section II and Fig. 1 setup; needed for the generalized-energy formulation and to avoid boundary effects.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Transition to thermal equilibrium in a deformed crystal." pith.science (2026). https://pith.science/paper/FFAESECY

@misc{pith2026190805464,
  author       = {Pith},
  title        = {Pith review of: Transition to thermal equilibrium in a deformed crystal},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FFAESECY}},
  note         = {Machine review of arXiv:1908.05464}
}

read the original abstract

An adiabatic transition between two equilibrium states corresponding to different stiffnesses in an infinite chain of particles is studied. Initially, the chain particles have random displacements and random velocities corresponding to a uniform initial temperature. An instant change of parameters of chain initiates a transient process. Analytical expressions for the chain temperature as a function of time are obtained from the statistical analysis of the dynamics equations. It is shown that the transition process is oscillatory and it converges non-monotonically to a new equilibrium state. Such behavior is usually unexpected for thermal processes. The analytical results are supplemented by numerical simulations.

Figures

Figures reproduced from arXiv: 1908.05464 by the authors.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

[18]

Krivtsov, A. M. Doklady Physics, 59, 9 (2014)

work page 2014
[1]

V., Morozov, N

Goldstein, R. V., Morozov, N. F. Physical Mesomechan- ics, 10, 5-6 (2007)

work page 2007
[2]

F., Golovneva, E

Golovnev, I. F., Golovneva, E. I., Fomin, V. M. Comp. Mat. Sci. 36, 176 (2006)

work page 2006
[3]

N., Alyokhin, V

Korobeynikov, S. N., Alyokhin, V. V., Annin, B. D., Babichev, A. V. Archives of Mechanics. 64, 367 (2012)

work page 2012
[4]

M., Morozov, N

Krivtsov, A. M., Morozov, N. F. Physics of the Solid State, 44, 12 (2002)

work page 2002
[5]

A., Murzaev, R

Baimova, Y. A., Murzaev, R. T., Dmitriev, S. V. Physics of the Solid State, 56, 10 (2014)

work page 2014
[6]

Vikesland, P. J. Nature Nanotechnology, 13, 8, (2018)

work page 2018
[7]

Simulation and Control of Chaotic Nonequilibrium Systems: Advanced Series in Nonlinear Dynamics: V

Hoover, W.G., Hoover, C.G. Simulation and Control of Chaotic Nonequilibrium Systems: Advanced Series in Nonlinear Dynamics: V. 27. World Scientiﬁc, Singapore. (2015). 324 p

work page 2015

Show all 39 references

[8]

A., Krivtsov, A

Kuzkin, V. A., Krivtsov, A. M. Discrete and continuum thermomechanics, Encyclopedia of Continuum Mechan- ics, Springer-Verlag, (2018)

2018
[9]

Kuzkin, V. A. Continuum Mech. Thermodyn., 31 (2019)

2019
[10]

Bass, R. Phys. Rev. B, 32, 4 (1985)

1985
[11]

A., Krivtsov, A

Kuzkin, V. A., Krivtsov, A. M. Phys. Stat. Sol. b, 252, 7 (2015)

2015
[12]

Lepri, S., Livi, R, Politi, A. Phys. Rep. 377, 1 (2003)

2003
[13]

L., Lieb, E., J

Rieder, Z., Lebowitz, J. L., Lieb, E., J. Math. Phys. 8, 5 (1967)

1967
[14]

Kannan, V., Dhar, A., Lebowitz, J. L. Phys. Rev. E 85, 4 (2012)

2012
[15]

Xiong, D., Zhang, Y., Zhao, H. Phys. Rev. E,88 5 (2013)

2013
[16]

P., Tildesley, A

Allen, M. P., Tildesley, A. K. Clarendon Press, Oxford. (1987). 390 p

1987
[17]

Physica 19 (1953)

Klein, G., Prigogine, I. Physica 19 (1953)

1953
[19]

A., Krivtsov, A

Kuzkin, V. A., Krivtsov, A. M. Journal of Physics: Con- densed Matter, 29, 50 (2017)

2017
[20]

A., Krivtsov, A

Kuzkin, V. A., Krivtsov, A. M. Physics of the Solid State, 59, 5 (2017)

2017
[21]

Murachev, A. S., A.M. Krivtsov, A. M., Tsvetkov, D.V. Journal of Physics: Condensed Matter, 31, 9 (2019)

2019
[22]

N., Krivtsov, A

Gavrilov, S. N., Krivtsov, A. M. Phys. Rev. E, 100, 2 (2019)

2019
[23]

A., Krivtsov, A

Sokolov, A. A., Krivtsov, A. M., M¨ uller, W.,H., Vilchevskaya, E. N. Phys. Rev. E, 99, 4 (2019)

2019
[24]

E., Kuzkin V

Berinskii I. E., Kuzkin V. A. Philosophical Transactions of the Royal Society A: Mathematical, Physical and En- gineering Sciences, 378, 2162 (2019)

2019
[25]

Sederberg, S., Kong, F., Corkum, P. B. . Physical Review X, 10(1) (2020)

2020
[26]

International Journal of Fracture, 142, 3-4, (2007)

Zhang, H., Ravi-Chandar, K. International Journal of Fracture, 142, 3-4, (2007)

2007
[27]

A., Petrov, Y

Morozov, V. A., Petrov, Y. V., Sukhov, V. D. Technical Physics, 64,5 (2019)

2019
[28]

V., Portnoi, M

Kruglyak, V. V., Portnoi, M. E. Generation of fem- tosecond current pulses using the inverse magneto-optical Faraday eﬀect. Technical Physics Letters, 31, 12 (2005)

2005
[29]

Lawrence and D.G

Toyserkani, E., Rasti, N., Ultrashort pulsed laser surface texturing, in Laser Surface Engineering, Edited by: J. Lawrence and D.G. Waugh, Woodhead Publishing, 718 6 p. (2015)

2015
[30]

Tzallas, P., Charalambidis, D., Papadogiannis, N., Witte, K., Tsakiris, G. D. Nature, 426, (267) (2003)

2003
[31]

APL Photonics, 4, (8) (2019)

Orfanos, I., Makos, I., Liontos, I., Skantzakis, E., F¨ org, B., Charalambidis, D., and Tzallas, P. APL Photonics, 4, (8) (2019)

2019
[32]

The European Physical Journal D, 72, 9 (2018)

Feng, L., Li, Y. . The European Physical Journal D, 72, 9 (2018)

2018
[33]

and Ulam S

Fermi E., Pasta J. and Ulam S. Los Alamos report LA-1940 (1955), published later in Collected Papers of Enrico Fermi, E. Segr´ e (Ed.), University of Chicago Press (1965)

1955
[34]

P., Izrailev, F

Berman, G. P., Izrailev, F. M. Chaos: An Interdisci- plinary Journal of Nonlinear Science, 15,1 (2005)

2005
[35]

and Stegun, I

Abramowitz, M. and Stegun, I. A. Handbook of Math- ematical Functions With Formulas, Graphs, and Mathe- matical Tables. New York: Dover. (1964) 1046 p

1964
[36]

Chemical Bonding in Solids

Burdett Jeremy K. Chemical Bonding in Solids. New York: Oxford University Press, 1995: 152

1995
[37]

International Journal of Solids and Structures

Berinskii, I.E., Krivtsov, A.M. International Journal of Solids and Structures. 96, 152 (2016)

2016
[38]

From nonlinear oscillations to equation of state in simple discrete systems Chaos, Solitons & Frac- tals 17(1), (2003), pp

Krivtsov A.M. From nonlinear oscillations to equation of state in simple discrete systems Chaos, Solitons & Frac- tals 17(1), (2003), pp. 79–87

2003
[39]

J., Johnston C., Chalker P

Moelle C., Klose S., Sz¨ ucs F., Fecht H. J., Johnston C., Chalker P. R., Werner M. Diamond and Related Mate- rials, 6(5-7) (1997)

1997

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [18]

Krivtsov, A. M. Doklady Physics, 59, 9 (2014)

work page 2014

[2] [1]

V., Morozov, N

Goldstein, R. V., Morozov, N. F. Physical Mesomechan- ics, 10, 5-6 (2007)

work page 2007

[3] [2]

F., Golovneva, E

Golovnev, I. F., Golovneva, E. I., Fomin, V. M. Comp. Mat. Sci. 36, 176 (2006)

work page 2006

[4] [3]

N., Alyokhin, V

Korobeynikov, S. N., Alyokhin, V. V., Annin, B. D., Babichev, A. V. Archives of Mechanics. 64, 367 (2012)

work page 2012

[5] [4]

M., Morozov, N

Krivtsov, A. M., Morozov, N. F. Physics of the Solid State, 44, 12 (2002)

work page 2002

[6] [5]

A., Murzaev, R

Baimova, Y. A., Murzaev, R. T., Dmitriev, S. V. Physics of the Solid State, 56, 10 (2014)

work page 2014

[7] [6]

Vikesland, P. J. Nature Nanotechnology, 13, 8, (2018)

work page 2018

[8] [7]

Simulation and Control of Chaotic Nonequilibrium Systems: Advanced Series in Nonlinear Dynamics: V

Hoover, W.G., Hoover, C.G. Simulation and Control of Chaotic Nonequilibrium Systems: Advanced Series in Nonlinear Dynamics: V. 27. World Scientiﬁc, Singapore. (2015). 324 p

work page 2015