[논문 리뷰] Thermal properties of materials from ab-initio quasi-harmonic phonons

Stefano Baroni, Paolo Giannozzi외 1명|arXiv (Cornell University)|2011. 12. 21.

High-pressure geophysics and materials참고 문헌 135인용 수 158

한 줄 요약

이 논문은 열역학적 성질인 열용량, 엔트로피, 격자 진동 기여도를 압력과 온도의 광범위한 범위에서 정확하게 예측하기 위해 밀도함수이론(DFT)과 준고전진동모형(QHA)을 조합한 최초의 원리적 접근을 제시한다. 이 방법은 영점 진동 효과를 잘 포착하며, 규소산염, 산화물 및 고압 상에서의 실험적 경향을 잘 재현하여 지구물리학적 조건에서 경험적 잠재력보다 뛰어난 성능을 발휘한다.

ABSTRACT

This paper gives a short overview of the calculation of thermal properties of materials from first principles, using the Quasi-Harmonic Approximation (QHA). We first introduce some of the thermal properties of interest and describe how they can be calculated in the framework of the QHA; then we briefly recall Density-Functional Perturbation Theory as a tool for calculating phonons from first principles, and present some codes that implement it; finally we review recent applications of first-principle QHA.

연구 동기 및 목표

• 실험이 불가능한 극한의 압력과 온도 조건에서 물질의 열역학적 성질을 예측할 수 있는 프레임워크를 개발하기 위해.
• 반경험적 상호작용 잠재력의 한계를 극복하기 위해, 진동 기여도를 준고전진동모형을 통해 원리적 DFT를 사용하기 위해.
• 특히 융해점 근처와 고압 조건에서 고체 내 양자적 영점 진동과 비선형 효과를 정확히 고려하기 위해.
• 지구 내부에서 중요한 역할을 하는 규소산염, 산화물, 퍼보스카이트와 같은 지구물리학적 조건에서의 물질을 모델링할 수 있는 신뢰할 수 있는 계산 도구를 제공하기 위해.
• 그리드 컴퓨팅과 같은 확장 가능한 계산 인프라를 활용해 열적 및 구조적 거동에 대한 대규모 시뮬레이션을 가능하게 하기 위해.

제안 방법

• 온도에 따른 격자 역학을 모델링하기 위해 격자 상수나 압력에 따라 진동 주파수를 변화시킬 수 있는 준고전진동모형(QHA)을 사용한다.
• 조화 진동 주파수와 고유벡터를 원리적 전자 구조 계산에서 밀도함수퍼투urbation이론(DFPT)을 사용해 계산한다.
• 영점 에너지 보정을 포함한 진동 모드에 대한 통합을 통해 자유 에너지를 부피의 기능으로 계산한다.
• 표준 통계역학 관계를 사용해 자유 에너지에서 열용량, 엔트로피, 열팽창 계수와 같은 열역학적 성질을 유도한다.
• 진동 스펙트럼과 열역학 반응 함수의 온도 및 압력 의존성을 계산하기 위해 이 방법을 적용한다.
• 대규모 병렬 및 그리드 컴퓨팅 아키텍처에서 확장 가능한 DFT+DFPT 계산을 수행하기 위해 Quantum ESPRESSO 소프트웨어 스위트를 사용한다.

실험 결과

연구 질문

• RQ1 준고전진동모형과 함께 원리적 DFT를 조합하여 지구물리학적 조건에서의 압력과 온도에서 물질의 열역학적 성질을 정확히 예측할 수 있는가?
• RQ2 코룬덤과 스티쇼바이트와 같은 고압 상의 안정성과 열역학에 영점 진동 효과가 얼마나 중요한가?
• RQ3 경험적 원자 간 잠재력은 칼세이트와 아라곤사이트와 같은 산화물의 올바른 열역학적 거동을 얼마나 잘못 포착하는가?
• RQ4 저온에서 α-Al₂O₃ 상이 진동으로 안정화되는 방식은 고압 다형체와 비교해 어떻게 다른가?
• RQ5 QHA-DFT 방법은 SiO₂와 MgSiO₃ 다형체의 실험적 열용량과 열팽창 계수를 얼마나 잘 재현하는가?

주요 결과

• QHA-DFT 방법은 α-쿼츠와 스티쇼바이트 SiO₂의 광범위한 온도 범위에서 실험적 열용량을 양적 정확도로 잘 재현한다.
• 영점 진동 에너지는 스티쇼바이트의 열역학적 성질에 큰 영향을 미치며, 특히 저온에서 열팽창보다 더 강한 영향을 미친다.
• QHA는 저온에서 영점 진동이 α-Al₂O₃ 상을 안정화시키지만, 실온에서는 진동만으로는 안정화되지 않음을 예측한다.
• 칼세이트(CaCO₃)의 경우, 고온에서 영점 진동이 평형 부피에 기여하는 정도가 열팽창보다 더 크며, 이는 그 중요성을 강조한다.
• 이 방법은 α-Al₂O₃ 상이 고압에서 CaIrO₃ 및 U₂S₃ 유형의 다형체로 전이됨을 밝혀 실험적 상도와 일치한다.
• QHA-DFT 방법은 초고압 조건에서 SiO₂의 페라이트-코탄나이트 상전이를 규명하였으며, 이는 외계행성 내부에 관련된 것으로 보고된 지질학적 불연속성과는 대조된다.