[논문 리뷰] Large temperature-up-jump simulations of a binary Lennard-Jones system

Aude Amari, Lorenzo Costigliola외 1명|arXiv (Cornell University)|2026. 01. 25.

Material Dynamics and Properties인용 수 0

한 줄 요약

이 논문은 대규모 온도 상승에 노출된 이진 Kob–Andersen–유형 Lennard-Jones 액체에 대해 도구-나라야나스와미 물질-시간 노화 설명을 시험하고, 시간-상관 함수가 물질-시간 차이에 대하여 collapse 되는지 시뮬레이션으로 평가한다.

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)?

연구 동기 및 목표

• 대규모 온도 상승 후 노화를 선형화할 수 있는 단일 물질 시간은 퍼텐셜 에너지로부터 도출되는가를 평가한다.
• 노화 중 퍼텐셜에너지 자기상관에 대한 삼각 관계(CK 시간 재매개변수화)를 평가한다.
• 다섯 가지 관측량에 대해 물질-시간 차이에 대한 Collapse가 가능한지 여부를 평가한다.
• 최종 온도로의 소-또는 대규모 상승에 따른 노화 특성을 비교하고 TN 노화 설명의 한계를 평가한다.

제안 방법

• 샘플링된 이진 Lennard-Jones 혼합물(80% A, 20% B)을 사용하고, 시프트-포스 컷오프를 적용한 GPU 최적화 코드 RUMD로 시뮬레이션한다.
• Cuu(t1,t2)=phi_eq(xi(t2)-xi(t1))를 통해 퍼텐셜에너지 자기상관으로부터 물질 시간 xi(t)를 정의한다.
• 온도 상승 후 수백만 개의 시간 삼중항을 순회하며 CK 삼각 관계 C13=F_A(C12,C23)을 검증한다.
• aging 중 다섯 가지 양을 모니터링한다: Cuu, Fs, <Delta r^2>, chi4, alpha2를 T=0.43 및 T=0.37에서 T=0.48로의 상승에서 관찰한다.
• xi2−xi1에 대해 그래프를 그렸을 때 이 관찰량들이 t2−t1이 아닌 물질 시간 차이에 대해 collapse하는지 평가한다.
• 소-상승과 대-상승에 대한 결과를 비교한다.

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 액체의 대규모 온도 상승 후 노화를 적절히 설명하는가?
• RQ2 오름 상승 후 노화 중 퍼텐셜에너지 자기상관에 대해 CK 삼각 관계가 얼마나 잘 성립하는가?
• RQ3 다수의 동적 관찰량이 물질-시간 차이의 보편적 함수로 collapse하는가, 그리고 이 collapse가 상승 크기에 의존하는가?
• RQ4 상승이 크면 노화를 설명하기 위해 지역적이거나 다중 물질 시간을 필요로 하는가?

주요 결과

• 삼각 관계는 퍼텐셜에너지 자기상관에 대해 잘 준수되며, 퍼텐셜 에너지에 기반한 물질 시간을 가능하게 한다.
• 작은 상승(0.43→0.48)에서 공통 물질 시간을 사용할 때 관찰된 노화 observables의 collapse가 더 잘 나타나며, 큰 상승(0.37→0.48)에서는 그렇지 않다.
• 작은 상승의 경우 Cuu와 Fs는 강한 collapse를 보이고, MSD, chi4, alpha2는 약하지만 부분적인 collapse를 보인다.
• 더 큰 상승의 경우 다섯 관찰량 중 어느 것도 단일 물질 시간으로 설득력 있는 collapse를 보이지 않아 단일 시계 TN 노화의 한계를 시사한다.
• 물질 시간은 최종 온도에서의 평형에 접근함에 따라 장기간에는 실제 시간에 비례하게 되지만, 짧은 시간대에는 편차가 나타난다.
• 저자들은 각 관찰량이 자체 물질 시간을 가질 수 있거나, 동적 이질성이 존재하는 경우 지역 물질 시간이 필요할 가능성을 시사한다.

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