Universal Dynamic Scaling for Sub-Yield Aging and Over-Yield Rejuvenation in A Cyclically Driven Metallic Glass

기술 요약: 주기적 구동 하의 금속 유리에서 나타나는 하부 항복 노화 및 상부 항복 재활성화에 대한 보편적 동적 스케일링

문제 제기
주기적인 기계적 변형 하에서의 비정질 고체, 특히 금속 유리(MG)의 구조적 진화는 근본적인 이중성을 나타낸다. 변형 진폭($\gamma_A$)에 따라, 주기적 구동은 시스템을 더 깊고 안정적인 에너지 최솟값으로 이끄는 기계적 노화(하부 항복 영역) 또는 시스템을 더 얕고 높은 에너지 상태로 몰아넣는 구조적 재활성화(상부 항복 영역)를 유도한다. 이와 같이 상반되는 궤적이 서로 다른 미시적 메커니즘에서 기인하는지, 아니면 공통된 물리적 기원을 공유하는지에 대해서는 오랜 논쟁이 지속되어 왔다. 또한, 하부 항복 완화 역학(sub-yield relaxation kinetics)과 상부 항복 비선형 유동 유변학(over-yield non-linear flow rheology)을 연결하는 통일된 정량적 프레임워크, 특히 비단조적 에너지 진화가 발생하는 중간 변형 진폭에서의 전이 역학에 관한 통합적 틀은 현저히 부재한 상태였다.

방법론
저자들은 16,000개의 원자로 구성된 전형적인 이성분 금속 유리 시스템인 Cu$_{64}$Zr$_{36}$를 대상으로 LAMMPS 패키지를 이용한 광범위한 고전 분자 동역학(MD) 시뮬레이션을 수행하였다. 원자 간 상호작용은 검증된 임베디드 원자 모델(EAM) 포텐셜을 사용하여 모델링되었다.

• 시뮬레이션 프로토콜: 초기 구성은 액상(liquid phase)을 2000 K에서 300 K까지 1 K/ps의 속도로 냉각(quenching)한 후 열적 완화를 거쳐 준비되었다.
• 변형: x-방향을 따라 사인 함수 형태의 변형 프로토콜, $\gamma(t) = \gamma_A \sin(2\pi t/t_p)$를 사용하여 주기적 단축 인장-압축을 적용하였으며, 이는 넓은 범위의 변형 진폭($\gamma_A$)과 로딩 주기($t_p$)를 포괄한다.
• 데이터 분석: 포텐셜 에너지 변화($\Delta U$)를 탄성 기여분을 제거하기 위해 변형률이 0인 상태에서의 로딩 사이클 수($N_{cycle}$)에 대한 함수로 추적하였다. 통계적 평균은 이소구성법(isoconfigurational method)을 통해 생성된 10개의 독립적인 구성에 대해 계산되었다.

주요 기여 및 결과

1. 구별된 스케일링 영역의 식별:
   본 연구는 ($\gamma_A$, $t_p$) 매개변수 공간 전체의 동적 응답을 매핑하여 네 가지 작동 영역을 식별하였다:

   • 영역 1 (하부 항복): 단조적인 기계적 가속 노화.
   • 영역 2 & 3 (교차 영역): 비단조적 궤적(초기 재활성화 후 노화 발생)을 특징으로 하는 동적 전이대.
   • 영역 4 (상부 항복): 정상 상태 에너지 고점에 도달하는 지속적인 기계적 재활성화.

2. 보편적 스케일링 법칙:
   저자들은 두 가지 구별된 보편적 스케일링 지수를 발견하였다:

   • 상부 항복 재활성화 (영역 4): 정상 상태 에너지 고점($\Delta U_p$)은 허셜-벌클리(Herschel-Bulkley, H-B) 유변 모델에서 유도된 거시적 소성 작업량과 스케일링된다. 스케일링 관계식은 $\Delta U_p \propto \gamma_A \dot{\gamma}^\alpha$이며, 여기서 전단 박화 지수는 $\alpha \approx 0.17$이다.
   • 하부 항복 노화 (영역 1): 에너지 진화는 $\dot{\gamma}^\beta / \gamma_A$에 의해 스케일링될 때 마스터 커브를 따른다. 이때 노화 지수는 $\beta \approx 0.34$이다. 에너지 진화는 스트레치드 지수(stretched-exponential) 감쇠 함수로 피팅된다.

3. 통합 스케일링 항등식 ($\beta \approx 2\alpha$):
   중심적인 발견은 $\beta \approx 2\alpha$ (0.34 $\approx$ 2 $\times$ 0.17)라는 정량적 관계이다. 저자들은 이 상반된 두 영역을 설명하기 위해 다음과 같은 통합된 물리적 메커니즘을 제안한다:

   • 1차 응답 (상부 항복): 재활성화는 거시적 소성 전단 작업($\propto \dot{\gamma}^\alpha$)에 의해 직접적으로 구동된다.
   • 2차 응답 (하부 항복): 전역적인 소성 흐름이 없는 상태에서, 구조적 노화는 2차 국부 탄성 응력 변동($\langle \delta\sigma^2_{local} \rangle$)에 의해 매개된다. 이러한 변동은 구동되는 매트릭스의 응력 척도($\propto \dot{\gamma}^\alpha$)를 상속받는 기계적 노이즈원 역할을 한다. 여기서 교란 에너지(agitation energy)는 $\langle \delta\sigma^2_{local} \rangle \propto (\dot{\gamma}^\alpha)^2 = \dot{\gamma}^{2\alpha}$로 스케일링된다.
   • 메커니즘: 강렬한 고율 교란은 깊은 포텐셜 에너지 최솟값을 탈출할 확률을 높여, 구조적 완화의 유지 효율을 억제한다. 저자들은 $\beta \approx 2\alpha$임을 입증함으로써, 하부 항복 완화와 상부 항복 유변학을 연결하는 타당한 물리적 기원을 제시하였다.

4. 교차 역학 ($N_c$):
   중간 진폭(영역 2 및 3)의 경우, 시스템은 특성 사이클 수 $N_c$에서 재활성화에서 노화로의 과도기적 교차를 보인다. 저자들은 $N_c$를 1차 소성 작업 주입과 2차 탄성 변동 소산 사이의 동적 균형으로 정의한다. 또한 보편적 스케일링 관계식 $N_c \dot{\gamma}^{2\alpha} = G(\dot{\gamma}^\alpha / \gamma_A)$를 도출하였으며, 여기서 $G$는 보편적 스케일링 함수이다. 이 "동적 시계(dynamic clock)"는 단기적 흥분(재활성화)에서 장기적 집합적 완화(노화)로의 전이를 결정한다.

의의 및 주장
본 논문은 하부 항복 노화와 상부 항복 재활성화 사이의 이중성을, 에너지 랜드스케이프 탐색을 규정하는 통합된 동적 시계를 확립함으로써 해결했다고 주장한다.

• 이론적 가교: 이는 비선형 유변학(Herschel-Bulkley 프레임워크)과 응력 변동 역학을 직접 연결하여, 구조적 진화를 예측할 수 있는 단일 프레임워크를 제공한다.
• 메커니즘적 통찰: 하부 항복 노화를 단순히 수동적인 완화가 아니라, 2차 국부 탄성 응력 변동이 "기계적 촉매"이자 노이즈원으로 작용하는 과정으로 재정의한다.
• 예측 능력: 스케일링 항등식 $\beta \approx 2\alpha$와 교차 경계의 마스터 콜랩스(master collapse)를 식별함으로써, 본 연구는 주기적 기계적 훈련을 통해 비정질 재료의 구조적 상태를 설계할 수 있는 예측 가능한 프로토콜을 제공한다.

저자들은 특정 지수 값($\alpha \approx 0.17, \beta \approx 0.34$)이 Cu$_{64}$Zr$_{36}$ 모델의 단축 NPT 조건에 특화된 것이지만, 변동 기반 프레임워크와 $\beta \approx 2\alpha$라는 스케일링 관계는 응력 변동과 완화 역학 사이의 근본적인 연결을 나타내며, 이는 다양한 로딩 모드에서 추가적인 검증이 필요함을 명시하였다.