[논문 리뷰] ElATools: A tool for analyzing anisotropic elastic properties of the 2D and 3D materials

Shahram Yalameha, Zahra Nourbakhsh외 1명|arXiv (Cornell University)|2021. 05. 15.

Chemical and Physical Properties of Materials참고 문헌 3인용 수 5

한 줄 요약

ElATools는 2차 에너지 강성 텐서를 사용하여 2차원 및 3차원 물질의 이방성 탄성 특성을 분석하는 사용자 友好的한 Fortran 90으로 작성된 계산 도구이다. Voigt, Reuss 및 Hill 평균화 방법을 통해 부피율, 전단율, 엘라스틱 모odulus, 푸아송 비율 및 이방성 지수와 같은 주요 기계적 성질을 계산하고, 통합된 플로팅 라이브러리를 통해 음향 속도 및 탄성 이방성을 시각화한다. 이 도구는 음수 푸아송 비율 및 음수 선형 압축성과 같은 이례적인 거동를 식별하여 나노스케일 장치 설계를 위한 예측 분석을 가능하게 한다.

ABSTRACT

We introduce a computational method and a user-friendly code with a terminal-based graphical user interface (GUI), named ElATools, developed to analyze mechanical and anisotropic elastic properties. ElATools enables facile analysis of the second-order elastic stiffness tensor of two-dimensional (2D) and three-dimensional (3D) crystal systems. It computes and displays the main mechanical properties including the bulk modulus, Young's modulus, shear modulus, hardness, p-wave modulus, universal anisotropy index, Chung-Buessem anisotropy index, log-Euclidean anisotropy parameter, Cauchy pressures, Poisson's ratio, and Pugh's ratio, using three averaging schemes of Voigt, Reuss, and Hill. It includes an online and offline database from the Materials Project with more than 13,000 elastic stiffness constants for 3D materials. The program supports output files of the well-known computational codes IRelast, IRelast2D, ElaStic, and AELAS. Four types of plotting and visualization tools are integrated to conveniently interface with GNUPLOT, XMGRACE, view3dscene, and plotly libraries, offering immediate post-processing of the results. It can efficiently identify anomalous mechanical properties, such as negative linear compressibility, negative Poisson's ratio, and highly-anisotropic elastic modulus in 2D and 3D materials, which are central properties to design and develop high-performance nanoscale electromechanical devices. Six case studies on selected material systems, namely, ZnAu$_2$ (CN)$_4$, CrB$_2$, $δ$-phosphorene, Pd$_2$O$_6$Se$_2$ monolayer, and GaAs, and a hypothetical set of systems with cubic symmetry are presented to demonstrate the descriptive and predictive capabilities of ElATools.

연구 동기 및 목표

• 결정성 물질에서의 실험적 탄성 상수 부족 문제를 해결하기 위해 고속 계산 분석을 가능하게 하기 위해.
• 다양한 2차원 및 3차원 물질을 대상으로 기계적 및 이방성 탄성 특성을 통합된 플랫폼에서 분석하기 위해.
• 고급 나노스케일 장치에 핵심적인 음수 푸아송 비율 및 음수 선형 압축성과 같은 이례적인 기계적 거동을 탐지하기 위해.
• 탄성 텐서의 고유값 분석을 통해 신뢰할 수 있는 기계적 안정성 평가를 지원하기 위해.
• 주요 DFT 코드 및 시각화 도구와 통합되어 결과를 즉각적인 후처리에 활용할 수 있도록 하기 위해.

제안 방법

• VASP, WIEN2k, Quantum Espresso 및 ElaStic, IRelast 등의 외부 패키지를 포함한 DFT 기반 코드를 통해 계산된 2차 탄성 강성 텐서를 입력으로 사용한다.
• Voigt, Reuss 및 Hill 평균화 방법을 적용하여 부피율, 전단율, 엘라스틱 모odulus 등 효과적인 기계적 성질을 계산한다.
• 유니버설 이방성(A_U), 로그-유클리드(A_L), Chung-Buessem(A_CB), Reuss(A_R)와 같은 이방성 지수를 계산하여 탄성 이방성을 정량화한다.
• 행렬 대각화를 통해 크리스토펠 방정식을 해결하여 음향 위상 속도 및 군 속도(파워 플로우 각도 포함)를 계산한다.
• G NUPLOT, X MGRACE, view3dscene, Plotly와 통합하여 탄성 성질의 실시간 2차원 및 3차원 시각화를 지원한다.
• IRelast, IRelast2D, ElaStic 및 AELAS 출력 형식의 입력을 지원하며, Materials Project에서 유래한 13,000개 이상의 탄성 상수를 포함한 오프라인 데이터베이스에 액세스할 수 있다.

실험 결과

연구 질문

• RQ1 2차 탄성 강성 텐서에서 2차원 및 3차원 물질의 이방성 탄성 특성을 체계적으로 계산하고 시각화하는 방법은 무엇인가?
• RQ2 탄성 텐서의 고유값 기반으로 2차원 및 3차원 물질의 기계적 안정성 기준은 무엇인가?
• RQ3 음수 푸아송 비율 또는 음수 선형 압축성을 포함한 이례적인 기계적 거동을 가장 신뢰성 있게 식별하는 데 가장 적합한 이방성 지수는 무엇인가?
• RQ4 복잡한 결정 대칭성을 가진 물질에서 크리스토펠 방정식을 사용하여 음향 속도 이방성을 얼마나 정확하게 예측할 수 있는가?
• RQ5 다양한 대칭성을 가진 다양한 물질에서 탄성 성질을 효율적으로 처리하고 시각화할 수 있는 통합 계산 도구는 존재하는가?

주요 결과

• ElATools는 ZnAu₂(CN)₄가 무한한 이방성(A_M = ∞)과 음수 선형 압축성을 나타내어 이례적인 기계적 반응을 보임을 성공적으로 식별하였다.
• CrB₂의 경우 계산된 A_M = 2.227는 중간 수준의 이방성을 나타내며, 부피율은 최소 53.51 GPa에서 최대 119.22 GPa까지, 전단율은 최소 125.04 GPa에서 최대 197.40 GPa까지 변동한다.
• δ-화이트 포스포레닌의 경우 푸아송 비율은 최소 -0.267에 도달하고, 엘라스틱 모odulus는 최소 62.277 N/m에서 최대 142.877 N/m까지 변동하여 강한 이방성을 나타낸다.
• Pd₂O₆Se₂ 이중층은 최대 푸아송 비율 1.315와 큰 A_R = 4.257을 보이며 강한 탄성 이방성을 나타낸다.
• GaAs의 경우 SS 모드에서 위상 속도 이방성(A_M)은 1.35에 도달하고 군 속도 이방성은 최대 1.41까지 증가하여 음향 전파의 방향적 의존성이 뚜렷하게 나타난다.
• 이 도구는 ZnAu₂(CN)₄의 -51.689 TPa⁻¹와 같이 음수 선형 압축성을 보이는 물질과 δ-P의 -0.267과 같이 음수 푸아송 비율을 보이는 물질을 정확히 식별하여 이례적인 기계적 거동에 대한 예측 능력을 확인하였다.