Coarse-graining amorphous plasticity: impact of rejuvenation and disorder

기술 요약: 비정질 소성(Amorphous Plasticity)의 거친 입자화(Coarse-graining): 재활성화(Rejuvenation) 및 무질서도의 영향

문제 정의
에이징(aging), 재활성화(rejuvenation), 전단 밴딩(shear banding), 바우싱거 효과(Bauschinger effect)와 같은 현상으로 특징지어지는 비정질 고체의 기계적 거동은 탄성 매질을 통해 상호작용하는 국부적인 구조적 재배열(전단 변형, Shear Transformations)에서 발생한다. 원자 단위 시뮬레이션은 이러한 특징들을 포착할 수 있지만, 계산 비용이 매우 높다. 반면, 중간 규모의 탄소성 모델(Elasto-Plastic Models, EPMs)은 계산 효율적인 프레임워크를 제공하지만, 종종 원자 단위의 결과를 정량적으로 재현하는 데 실패할 수 있는 지나치게 단순화된 가정에 의존한다. 본 연구가 다루는 핵심 문제는 "거친 입자화(coarse-graining)" 과제이다. 즉, 비정질 소성의 물리학을 재현하기 위해 원자 규모에서 중간 규모의 격자 모델로 파라미터를 어떻게 전달할 것인가 하는 문제이며, 특히 기저에 깔린 무질서한 퍼텐셜 에너지 지형(Potential Energy Landscape, PEL), 하중 시 무질서도의 진화, 그리고 유리 상태의 재활성화 역할에 초점을 맞춘다.

방법론
본 연구는 비등방성 준정적(Athermal Quasistatic, AQS) 원자 단위 시뮬레이션과 단순화된 스칼라 중간 규모 EPM 사이의 비교 접근법을 사용한다.

1. 원자 단위 시뮬레이션:

   • 세 가지 서로 다른 열 이력을 가진 2D 이성분 레너드-존스(Lennard-Jones) 유리가 준비되었다:
     • MCT: 모드 결합 전이(Mode Coupling Transition, $T_p \approx 0.38$) 근처의 온도에서 즉각적으로 퀜칭(quench)된 샘플.
     • ESL: 더 높은 온도($T_p = 0.351$)에서 즉각적으로 퀜칭되어, 결함이 많은(연한) 유리를 나타내는 샘플.
     • GQ: $T_p = 0.351$에서 점진적으로 퀜칭되어, 잘 어닐링된(단단한) 유리를 나타내는 샘플.
   • 샘플들은 전역 변형률 $\gamma = 5$까지 AQS 전단 하중을 받았다.
   • 국부적 특성은 "동결 행렬 방법(frozen matrix method)"을 사용하여 측정되었으며, 이는 원형 영역을 격리하여 국부 항복 응력($\sigma_c$), 내부 응력 변동($\delta\sigma$), 그리고 18개의 전단 방향에 따른 응력 강하($\Delta\sigma$)를 측정한다.

2. 중간 규모 모델:

   • 스칼라 격자 모델이 이중 주기적 경계 조건을 가진 정사각형 그리드 상에 구현되었다.
   • 역학(Dynamics): 이 모델은 국부 응력이 임계값 $\sigma_c$를 초과할 때 셀이 항복하는 국부 임계값 역학을 특징으로 한다. 플라스틱 사건(plastic events)은 장범위 탄성 응력 재분배(Eshelby kernel)와 국부 임계값의 갱신을 유발한다.
   • 주요 파라미터: 연구는 두 가지 무차원 파라미터에 집중한다:
     • 무질서도 강도 ($s$): 항복 응력 변동과 플라스틱 사건에 의한 응력 강하 사이의 비율.
     • 에이징/재활성화 파라미터 ($r$): $r = (\sigma_{c}^{initial} - \sigma_{c}^{renewed}) / \sigma_{c}^{renewed}$로 정의된다. 이는 초기 유리 상태와 플라스틱 재배열 이후의 상태 사이의 차이를 정량화한다.
   • 업스케일링 전략: 중간 규모 모델의 파라미터(항복 응력 분포, 슬립 증분, 갱신 분포)는 원자 단위 데이터로부터 직접 추출되었다. 연구는 다음 두 가지 가설을 테스트했다:
     • 무질서 시나리오(Disordered Scenario): 항복 응력의 갱신 분포가 초기 분포와 동일함 ($r=0$).
     • 재활성화 시나리오(Rejuvenation Scenario): 갱신 분포가 MCT 유리(더 "신선한" 상태)에 해당하며, 이는 $r > 0$임을 의미함.
     • 상관관계 가설(Correlation Hypothesis): 국부 항복 응력과 후속 슬립 증분 크기 사이에 직접적인 상관관계($\varepsilon_0 \propto \sigma_c^2$)를 도입하였으며, 이는 PEL의 계층적 구조를 반영한다.

주요 기여 및 결과

1. 응력-변형률 거동 및 과도 현상:

   • 항복 응력과 슬립 증분 사이의 상관관계가 없는 단순 모델들은 원자 단위 시뮬레이션에서 관찰되는 과도 응력 피크를 재현하지 못하고 단조로운 곡선만을 생성했다.
   • 국부 항복 응력과 슬립 증분 진폭 사이의 상관관계를 도입함으로써, 중간 규모 모델은 연한(ESL) 유리와 단단한(GQ) 유리 모두에 대해 과도 피크와 정상 상태 유동 응력을 포함한 전체 응력-변형률 곡선을 재현할 수 있었다.
   • 정상 상태 유동 응력은 슬립 진폭의 평균값과 연결된 무질서도 강도 파라미터 $s$에 의해 주로 제어됨을 확인하였다.

2. 재활성화 및 국부 항복 응력 진화:

   • 연한 유리 (ESL): 원자 단위 시뮬레이션은 약간의 경화(평균 국부 항복 응력 증가) 후 플래토(plateau)를 보이는 것을 나타냈다. 재활성화($r \approx 0.14$)를 적용한 중간 규모 모델은 이러한 진화에 대해 준양적인 일치를 달성했다. 재활성화가 없으면 모델은 항복 응력을 과대평가했다.
   • 단단한 유리 (GQ): 원자 단위 시뮬레이션은 전단 밴딩과 관련된 초기 경화 후 연화(softening) 구간을 보였다. 재활성화가 없는 중간 규모 모델은 연화 단계를 재현하지 못하고 지속적인 경화만을 보여주었다. 재활성화 시나리오($r \approx 0.65$)를 포함하는 것이 평균 항복 응력의 감소를 질적으로 회복하는 데 결정적이었으나, 양적인 일치는 여전히 불완전했다.

3. 국부화 및 전단 밴딩:

   • 단단한 유리(GQ)의 경우, 원자 단위 시뮬레이션은 명확한 전단 밴딩을 나타냈다.
   • 재활성화와 상관관계가 적용된 중간 규모 모델은 국부화 지수(localization index)의 진화를 성공적으로 재현했다: 밴드 핵생성 동안 급격히 증가한 후 밴드가 넓어짐에 따라 점진적으로 감소함.
   • 재활성화가 없는 모델은 상당한 국부화를 보였으나, 전단 밴드의 강도와 지속성을 과소평가했다.

4. 변형 및 응력 변동:

   • 중간 규모 모델은 변형 분산의 진화 경향은 재현했으나, 원자 단위 결과에 비해 변동의 크기를 약 1.5배 체계적으로 과대평가했다.
   • 내부 응력장 변동은 훨씬 더 크게(2~3배) 과대평가되었다. 저자들은 이 불일치를 모델의 스칼라 성질(슬립 방향의 방위 무질서도를 무시함)과 중간 규모 공식에서의 탄성 불균질성 및 비정형(non-affine) 효과의 누락 때문이라고 설명한다.

의의 및 주장
본 논문은 특정 구조적 정보가 원자 단위 규모로부터 통합될 경우, 거친 입자화된 스칼라 모델이 비정질 소성의 주요 특징들을 재현할 수 있다고 주장한다. 주요 의의는 양적인 일치를 위한 두 가지 핵심 요소의 식별에 있다:

1. 상관관계: 국부적인 응력 강하를 국부 항복 응력과 상관시키는 것(더 깊은 우물은 더 큰 슬립에 대응한다는 PEL의 계층적 구조를 반영).
2. 재활성화: 플라스틱 재배열이 국부적인 유리 구조를 모드 결합 전이 근처에서 퀜칭된 유리 상태와 유사하게 변화시킨다는 재활성화 시나리오의 검증. 이 메커니즘은 잘 어닐링된 유리에 나타나는 연화 거동과 전단 밴딩을 재현하는 데 필수적이다.

저자들은 자신들의 모델이 완전한 양적 일치를 위해서는 "너무 단순하고 캐리커처 같다"고 겸허히 인정하며, 특히 변동 진폭에 관한 부분이다. 결론적으로, 스칼라 모델이 유동 응력과 국부화와 같은 핵심 물리학은 포착할 수 있지만, 변동 통계의 불일치를 완전히 해결하기 위해서는 향후 탄성 불균질성, 텐서 무질서도, 그리고 메모리 효과를 다루어야 한다고 밝히고 있다.