[논문 리뷰] The generalized quasiharmonic approximation via space group irreducible derivatives

Mark A. Mathis, Amey Khanolkar외 10명|arXiv (Cornell University)|2022. 02. 28.

Acoustic Wave Phenomena Research참고 문헌 93인용 수 18

한 줄 요약

이 논문은 공간군의 기저로 사용되는 기저 도함수를 활용하여 일반화된 준조화 근사(Quasiharmonic Approximation, QHA)를 수립함으로써, 결정에서 온도 및 응력에 의존하는 열역학적 성질을 정확하고 효율적으로 계산할 수 있도록 한다. 이중 두 가지 변형 파arametrization(격자 보간과 테일러 전개)을 사용하여 이차 이동 도함수의 비대칭성을 제거하고, 토리아(ThO2)에서 탄성 상수와 열팽창을 정밀하고 재현 가능한 방식으로 계산할 수 있다. 이는 제1원리 계산 및 중성자 산란 데이터와의 비교를 통해 검증되었다.

ABSTRACT

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. Here we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of the elastic constant tensor and the thermal expansion in thoria (ThO$_2$) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data driven applications.

연구 동기 및 목표

• 일반화된 준조화 근사(QHA)를 압력이 아닌 일정 응력 조건에서 작동하도록 개발한다.
• 임의의 변형에서 진동 모드를 평가하는 데 발생하는 계산 복잡성을 해결하기 위해 공간군 기반 기저 도함수를 사용한다.
• 응력에 따라 변하는 탄성 에너지 및 동역학 행렬을 계산하기 위한 실용적이고 수치적으로 안정된 알고리즘을 제공한다.
• 재료에서 온도 및 응력에 따라 변화하는 탄성 상수와 열팽창을 정확히 예측할 수 있도록 한다.
• 토리아(ThO2)에서 밀도함수 이론과 중성자 산란을 활용하여 실험 데이터와의 검증을 수행한다.

제안 방법

• 공간군 기반 기저 도함수를 사용하여 계산 비용을 줄이고 대칭성 일致성을 향상시키기 위해 일반화된 QHA를 수립한다.
• 유일한 기저 도함수 접근법을 활용하여 응력에 따라 변하는 이차 핵 이동 도함수를 효율적으로 계산한다.
• 두 가지 변형 파arametrization을 적용한다: (1) 이산 응력 격자 보간과 (2) 대칭화된 응력 변수에서의 테일러 급수 전개.
• 응력에 따라 변하는 동역학 행렬과 탄성 에너지를 제1원리 밀도함수 이론을 사용하여 다수의 상호작용-상관 기능을 활용해 계산한다.
• 응력 의존성에 따라 헬름홀츠 자유 에너지에서 기브스 자유 에너지로의 레지외르 변환을 수행한다.
• 시간 영역 브릴루앙 산란 및 비탄성 중성자 산란 측정을 통해 ThO2의 탄성 상수를 검증한다.

실험 결과

연구 질문

• RQ1 대칭 적응 도함수를 사용하여 일정 응력 조건에서 일반화된 QHA를 체계적으로 수립할 수 있는가?
• RQ2 응력에 따라 변하는 진동 모드를 파arametrization하기 위해 격자 보간과 테일러 급수 전개의 상대적 이점과 수치적 안정성은 무엇인가?
• RQ3 기저 도함수의 사용이 QHA 계산에서 계산 비용을 크게 줄이면서도 정확성을 유지하는 데 기여하는가?
• RQ4 일반화된 QHA는 토리아(ThO2)에서 탄성 상수와 열팽창의 온도 및 압력 의존성을 얼마나 정확하게 예측하는가?
• RQ5 시간 영역 브릴루앙 산란 및 비탄성 중성자 산란을 활용한 실험 측정과 비교했을 때 QHA 예측의 탄성 상수는 얼마나 잘 일치하는가?

주요 결과

• 기저 도함수 접근법을 통해 응력에 따라 변하는 동역학 행렬과 탄성 에너지를 효율적이고 정밀하게 계산할 수 있으며, 일반화된 QHA의 구현을 단순화한다.
• 테일러 급수와 격자 보간 양측 모두 정확한 결과를 도출하며, 테일러 급수는 응력에 대해 주어진 차수까지 해석적 제어를 가능하게 한다.
• 이 방법은 토리아(ThO2)에서 탄성 상수 텐서의 온도 및 압력 의존성을 성공적으로 예측하였으며, 시간 영역 브릴루앙 산란 및 비탄성 중성자 산란 측정 결과와 양호한 일치를 보였다.
• 다양한 상호작용-상관 기능을 사용한 계산에서 토리아(ThO2)의 열팽창 계수는 일관되게 유사한 값을 보이며, 이는 접근법의 강건성을 입증한다.
• 대칭화된 응력 변수와 기저 성분을 사용함으로써 파arametrization이 결정의 점군 대칭성을 유지하여 임의의 모드 혼합을 방지한다.
• 일반화된 QHA 프레임워크는 실험적으로 관련된 일정 응력 조건 하에서 재현 가능하고 데이터 기반의 열역학적 성질 예측을 가능하게 한다.