[논문 리뷰] Fifth-degree elastic potential for predictive stress-strain relations and elastic instabilities under large strain and complex loading in Si

Hao Chen, Nikolai A. Zarkevich외 3명|arXiv (Cornell University)|2020. 02. 14.

Microstructure and mechanical properties참고 문헌 37인용 수 4

한 줄 요약

이 논문은 밀도함수이론(DFT) 데이터로부터 유도된 실리콘(Si I)에 대한 오목도 5차 탄성 포텐셜을 제안하여 복잡한 하중 조건 하에서 큰 변형률의 응력-변형률 거동과 탄성 불안정성을 정확하게 모델링한다. 불안정성 점을 포함한 큰 변형률 범위에서 오차를 최소화함으로써, 모델은 DFT 결과와 우수한 일치를 보이며, 카우치 응력-라그랑주 변형률 곡선, 상전이 조건, 비틀림 불안정성 등에 대해 우수한 예측 성능을 보이며, 이러한 현상을 신뢰성 있게 기록하지 못하는 저차수 포텐셜보다 뛰어난 성능을 발휘한다.

ABSTRACT

Materials under complex loading develop large strains and often transition via an elastic instability, as observed in both simple and complex systems. Here, we present Si I under large strain in terms of Lagrangian strain by an $5^{th}$-order elastic potential found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress-Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including elastic instability driving Si I$ ightarrow$II phase transformation (PT) and the shear instabilities. PT conditions for Si I$ ightarrow$II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such elastic potential permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum simulations of crystal behavior under extreme loading.

연구 동기 및 목표

• 복잡한 하중 조건 하에서 실리콘(Si I)의 큰 변형률 비선형 탄성 거동과 탄성 불안정성을 정확히 기술할 수 있는 예측 가능한 탄성 포텐셜을 개발하는 것.
• 3차 및 4차 탄성 상수의 한계를 극복하여 큰 변형률에서 격자 불안정성을 기술하지 못하는 문제를 해결하는 것.
• 극한의 정적 및 동적 하중 하에서 상전이, 미세결정의 슬립, 투과, 균열을 시뮬레이션하기 위한 기본 연속체 프레임워크를 제공하는 것.
• 임의의 응력 상태 하에서 Si I → Si II 등의 상전이 조건을 정확히 예측할 수 있도록 하는 것.
• 다차원 하중 조건 하에서 정규 응력과 비틀림 응력 간의 비선형 상호작용을 포괄하는 고차수 탄성 포텐셜을 구축하여, 다차원 하중 하에서의 재료 거동에 핵심적인 영향을 미치는 요소를 반영하는 것.

제안 방법

• 탄성 포텐셜은 라그랑주 변형률 성분(6개의 독립 성분)에 대한 5차 다항식으로 표현되며, 압축 표현을 위한 보이트 표기법을 사용한다.
• 5차 포텐셜의 계수는 큰 변형률 범위(불안정성 점 포함)에서 예측된 응력-변형률 곡선과 DFT로 계산된 카우치 응력-변형률 곡선 간 오차를 최소화하여 결정된다.
• 모델은 훈련을 위한 기준으로 DFT 데이터를 사용하여 탄성 및 불안정 영역 모두에서 정확성을 확보하며, Si I → Si II 상전이의 시작 시점까지도 포함된다.
• 다양한 복잡한 하중 조건(일축, 이축, 수압, 정규-비틀림 복합 하중) 하에서의 응력-변형률 거동을 계산하고 DFT 결과와 직접 비교한다.
• 상전이 조건은 해석적으로 유도되며 DFT 예측 결과와 검증되어, 입방축 하중 하에서 카우치 응력에 선형적으로 의존하는 것으로 나타났다.
• 이 방법을 통해 다양한 하중 경로가 비틀림 강도, 임계 변형률 및 안정성에 미치는 영향을 체계적으로 분석할 수 있으며, 원자적 연성의 지표로 볼 수 있는 플랫폼 형태의 영역이 나타나는 현상까지도 분석할 수 있다.

실험 결과

연구 질문

• RQ1 오목도 5차 탄성 포텐셜은 큰 변형률과 복잡한 하중 조건, 특히 탄성 불안정성을 포함한 실리콘의 응력-변형률 거동을 정확히 예측할 수 있는가?
• RQ2 저차수 탄성 상수에 비해 5차 항을 포함함으로써 상전이 조건(Si I → Si II 등)의 예측 성능이 어떻게 향상되는가?
• RQ3 정규 응력과 비틀림 응력의 초월된 조합(예: 정규 및 비틀림 변형률의 중첩)이 이론적 비틀림 강도와 불안정성에 대한 임계 변형률에 미치는 영향은 어떠한가?
• RQ4 모델은 정규 응력과 비틀림 응력 간의 비선형 상호작용을 얼마나 잘 반영하여 안정성 및 상전이 경로에 영향을 미치는가?
• RQ5 이 접근법은 균열 및 투과를 포함한 극한 하중 조건 하에서 결정의 거동을 연속체 시뮬레이션하기 위한 일반적이고 예측 가능한 프레임워크로 활용될 수 있는가?

주요 결과

• 오목도 5차 탄성 포텐셜은 임의의 복잡한 하중 조건 하에서 DFT로 계산된 카우치 응력-라그랑주 변형률 곡선을 매우 우수한 일치로 재현하며, 변형률이 0.35를 초과하는 영역에서도 성능을 유지를 한다.
• 입방축 응력 하에서 Si I → Si II 상전이 조건을 정확히 예측하며, 전이 조건은 카우치 응력에 선형적으로 의존하며 DFT 예측 결과와 정확히 일치한다.
• 비틀림 방향에 일축 또는 이축 압축을 중첩함으로써 최대 2 GPa까지 최대 비틀림 강도가 감소하고 임계 비틀림 변형률도 변화하며, 등방성 압축 하에서는 4.4 GPa(실제 DFT 결과는 4.31 GPa)로 감소한다.
• 모델은 특정 하중 경로(예: η₁ = η₂ = -0.5η₄)가 비틀림 응력-변형률 곡선에서 플랫폼 형태의 영역을 크게 증가시켜 원자적 연성 증가를 나타내는 것으로 밝혀졌다.
• 3차 및 4차 포텐셜은 큰 변형률에서 올바른 불안정성 거동과 응력-변형률 곡선을 재현하지 못하며, 이는 5차 항이 반드시 필요하다는 것을 확인한다.
• 이 모델은 방향 의존적인 극한 하중 조건 하에서 다수의 불안정성(상전이, 비정질화, 슬립, 투과, 균열 등)을 예측 가능한 시뮬레이션할 수 있으며, 연속체 역학의 새로운 기반을 제공한다.