[論文レビュー] Large temperature-up-jump simulations of a binary Lennard-Jones system

Aude Amari, Lorenzo Costigliola+1 more|arXiv (Cornell University)|Jan 25, 2026

Material Dynamics and Properties被引用数 0

ひとこと要約

この論文は、Large temperature up-jumps に対してバイナリ Kob–Andersen 型 Lennard-Jones 液の Tool-Narayanaswamy 材料時刻 aging の記述を検証し、時刻-dependent の時間相関関数が材料時刻差に收束するかをシミュレーションで評価する。

ABSTRACT

This paper presents simulations of the physical aging of a binary Kob-Andersen-type Lennard-Jones liquid following large temperature up-jumps from equilibrated states of high relaxation time. The purpose is to investigate how well the Tool-Narayanaswamy (TN) material-time concept works for this rather extreme case of aging. First the triangular relation of the potential energy is studied. This is found to be well obeyed, making it possible to define a potential-energy-based material time $ξ$. We proceed to study aging toward equilibrium at the final temperature 0.48 for jumps from the two temperatures 0.43 and 0.37, monitoring the following five quantities: the potential energy, the self-intermediate scattering function, the mean-square displacement, the dynamic susceptibility $χ_4$, and the non-Gaussian parameter $α_2$. The TN material-time prediction is that all time-autocorrelation functions should collapse to only depend on the material-time difference $ξ_2-ξ_1$. This is found to work better for the $0.43 o 0.48$ temperature jump than for the $0.37 o 0.48$ jump, however. Our findings thus confirm the general understanding that the TN aging formalism works best for systems that are never very far from equilibrium. This raises two questions for future work: Is the collapse significantly improved if each aging quantity is allowed its own material time? Can better collapse be obtained if the material-time is generalized to be locally defined (in order to reflect dynamic heterogeneity)?

研究の動機と目的

• 大きな温度ジャンプ後の aging を線形化できる潜在エネルギーから導かれる単一の材料時刻が成立するかを評価する。
• aging 中の潜在エネルギー自己相関について CK 三角関係（CK time-reparametrization）を評価する。
• five observables の時間自己相関関数が材料時刻差に収束するかを検証する。
• final temperature へ向かう小さなジャンプと大きなジャンプの aging の挙動を比較し、TN aging の記述の限界を評価する。

提案手法

• 80% A, 20% B のバイナリ Lennard-Jones 系を使用し、シフトフォースカットオフを用いたGPU最適化コード RUMD でシミュレーションする。
• Cuu(t1,t2)=phi_eq(xi(t2)-xi(t1)) によってポテンシャルエネルギー自己相関から材料時刻 xi(t) を定義する。
• 温度ジャンプ後に多数の時刻三つ挙動をループして CK 三角関係 C13=F_A(C12,C23) の成り立ちを検証する。
• aging 中に5つの量をモニタする：Cuu, Fs, <Delta r^2>, chi4, alpha2。ジャンプは T=0.43 → 0.48 および T=0.37 → 0.48 から。
• これらの観測量を t2−t1 ではなく xi2−xi1 に対してプロットしたときの収束を評価する。
• 小さなジャンプと大きなジャンプの結果を比較する。

Figure 1 : Physical aging. (a) Schematic drawing of a temperature up jump starting and ending in thermal equilibrium (with time on the x-axis). Even though not all quantities relax in identical manner, the standard TN aging formalism operates with a single “global” clock controlling all relaxations

実験結果

リサーチクエスチョン

• RQ1 潜在エネルギーから導かれる単一の TN 材料時刻が Kob–Andersen 型 LJ 液における大きな温度ジャンプ後の aging を適切に記述するか。
• RQ2 aging 後の潜在エネルギー自己相関について CK 三角関係がどの程度成立するか。
• RQ3 複数の動的 observables が材料時刻差の普遍関数に収束するか、収束はジャンプサイズに依存するか。
• RQ4 ジャンプが大きい場合、 aging を記述するのに局所的または複数の材料時刻が必要か。

主な発見

• 潜在エネルギー自己相関に対して CK 三角関係は良好に成り立ち、潜在エネルギーを基にした材料時刻が可能である。
• 小さなジャンプ（0.43→0.48）に対して共通の材料時刻を用いると観測量はより良く収束するが、大きなジャンプ（0.37→0.48）では収束が悪い。
• 小さなジャンプで、Cuu と Fs は強く収束する。一方、MSD、chi4、alpha2 は弱いが部分的な収束を示す。
• 大きなジャンプでは、5つの観測量のいずれも単一材料時刻で納得の収束を示さず、単一時計 TN aging の限界を示唆する。
• 材料時刻は最終温度の平衡へ近づく長時間には実時間に比例するが、短時間では偏差が生じる。
• 著者らは、各観測量が独自の材料時刻を持つべきか、あるいはダイナミック・ヘテロジニティの存在下で局所的な材料時刻が必要かもしれないと示唆している。

Figure 2 : Parametric representation of the triangular relation, Eq. ( 2 ). (a) Illustration of the sampling of time triplets used to calculate correlation “triangles”. The values on the time axis illustrate the logarithmic increase in time intervals within a single simulation block. (b) Illustratio