[논문 리뷰] ElasTool: An automated toolkit for elastic constants calculation

Zhongli Liu|arXiv (Cornell University)|2020. 02. 16.

Machine Learning in Materials Science참고 문헌 29인용 수 10

한 줄 요약

ElasTool는 DFT를 사용하여 2차원 및 3차원 결정에서 두 번째 차수 탄성 상수(SOECs)를 계산하기 위한 자동화된 파이썬 툴킷입니다. OHESS, ASESS 및 ULICS 세 가지 변형 매트릭스 세트를 사용하여 응력-변형률 방법을 통해 계산하며, 고온 및 영온도에서의 탄성 성질 계산이 효율적이고 정확하게 가능합니다. 이는 탄성 모odulus, 비대칭성, 디바이 온도, 안정성 분석 등 다양한 물리적 성질을 포함하며, 고성능 재료 설계 및 탄성 연구를 크게 향상시킵니다.

ABSTRACT

We here present, the ElasTool package, an automated toolkit for calculating the second-order elastic constants (SOECs) of any two- (2D) and three-dimensional (3D) crystal systems. It can utilize three kinds of strain-matrix sets, the high-efficiency strain-matrix sets (OHESS), the universal linear-independent coupling strains (ULICS) and the all-single-element strain-matrix sets (ASESS) to calculate the SOECs automatically. In an automatic manner, ElasTool can deal with both zero- and high-temperature elastic constants. The theoretical background and computational method of elastic constants, the package structure, the installation and run, the input/output files, the controlling parameters, and two representative examples of ElasTool are described detailedly. ElasTool is useful for either the exploration of materials' elastic properties or high-throughput new materials design. ElasTool is also available at our website: www.matdesign.cn.

연구 동기 및 목표

• 2차원 및 3차원 결정계에서 두 번째 차수 탄성 상수(SOECs)를 계산하기 위한 자동화되고 효율적인 툴킷 개발.
• 최적화된 변형 매트릭스 세트(OHESS, ASESS, ULICS)를 활용하여 DFT 기반 탄성 상수 계산의 계산 비용을 감소.
• 정적 응력(영온도) 및 동적 응력(분자 동역학를 통한 고온) 계산을 모두 지원.
• 탄성 모odulus, 음속, 디바이 온도, 비대칭성, 기계적 안정성 분석 등 포괄적인 출력 제공.
• 탄성 성질 계산의 간소화를 통해 고성능 재료 선별 및 이론적 설계를 가능하게 함.

제안 방법

• 후크의 법칙에 기반한 응력-변형률 방법을 사용: σ_i = Σ C_ij ε_j, 여기서 응력은 다양한 변형된 구조에서 DFT로 계산됨.
• 세 가지 변형 매트릭스 세트 통합: OHESS(최적화된 고성능 변형 세트), ASESS(단일 원소 변형 매트릭스), ULICS(일반적인 선형 독립 결합 변형).
• 전체 워크플로우 자동화: 입력 파싱, 사용자 정의 명령어(예: VASP)를 통한 DFT 응력 계산, 응력-변형률 피팅, 탄성 상수 추출.
• VASP 또는 기타 코드에서 계산된 DFT 응력 텐서를 사용하고, 선형 응력-변형률 반응에서 탄성 상수를 결정하기 위해 최소 제곱 피팅 적용.
• 기본 세포 및 원자 세포 입력 모두 지원하며, 정적 또는 동적 응력 평균화(예: MD 시뮬레이션에서 유도) 옵션 제공.
• C_ij, 체적 및 비틀림 모odulus(Voigt, Reuss, Voigt-Reuss-Hill 평균), 엥겔모듈러스, 푸아송 비, 음속, 디바이 온도, 비대칭성, 보른 기준에 따른 안정성 분석 등 결과 출력.

실험 결과

연구 질문

• RQ1 최적화된 변형 매트릭스 세트를 사용하여 DFT 기반 탄성 상수 계산의 계산 효율성을 어떻게 향상시킬 수 있는가?
• RQ2 다양한 2D 및 3D 결정계에서 영온도 및 고온 탄성 성질을 모두 처리할 수 있는 자동화된 파이프라인을 개발할 수 있는가?
• RQ3 OHESS 변형 매트릭스 세트가 DFT 계산 시간을 최소화하면서도 정확도를 유지하는 데 있어 ULICS 및 ASESS와 비교해 성능과 정확도가 어떻게 되는가?
• RQ4 ElasTool는 DFT 응력 데이터에서 기계적 안정성, 탄성 비대칭성, 디바이 온도를 얼마나 효과적으로 예측할 수 있는가?
• RQ5 ElasTool는 탄성 성질 계산의 수동 조작을 줄임으로써 고성능 재료 스크리닝에 얼마나 기여하는가?

주요 결과

• ElasTool는 OHESS, ASESS 및 ULICS 세 가지 변형 매트릭스 세트를 사용하여 어떤 2D 또는 3D 결정계에서도 두 번째 차수 탄성 상수(SOECs)를 성공적으로 계산함.
• 다이아몬드(입방계)의 경우 계산된 SOECs는 C11 = 1055.04 GPa, C12 = 136.56 GPa, C44 = 567.76 GPa이며, 평균 압력은 -0.08 GPa임.
• 체적 모odulus의 Voigt 평균(B_V)과 Reuss 평균(B_R)은 모두 442.72 GPa이며, Voigt-Reuss-Hill 평균(B_VRH) 역시 442.72 GPa임.
• 비틀림 모odulus 값은 G_V = 524.36 GPa, G_R = 518.73 GPa, G_VRH = 521.54 GPa이며, 엥겔모듈러스 E = 1123.47 GPa, 푸아송 비 ν = 0.0771임.
• 음속 값은 V_S = 12.20 km/s, V_B = 11.24 km/s, V_P = 18.03 km/s, V_M = 13.31 km/s이며, 디바이 온도 T_D = 1761.03 K임.
• 보른 탄성 안정성 기준에 따라 기계적 안정성이 확인되었으며, 탄성 비대칭성 A_U = 0.0542 및 A_C = 0.0054임.