[논문 리뷰] FourPhonon: An extension module to ShengBTE for computing four-phonon scattering rates and thermal conductivity

Zherui Han, Xiaolong Yang외 3명|arXiv (Cornell University)|2021. 04. 10.

Thermal properties of materials참고 문헌 43인용 수 433

한 줄 요약

FourPhonon은 ShengBTE 프레임워크에 대한 오픈소스 확장으로, 적응형 넓힘 기반 기법을 사용하여 첫 번째 원리적 접근으로 4음파 산란률과 격자 열전도도를 계산할 수 있도록 한다. 기존 워크플로우와 통합되며, 대칭을 고려한 스크립트를 통해 4계 상호작용 힘 상수를 계산하고, 정확한 열전도도 예측을 제공한다 — Si, BAs, LiCoO2에서의 성능을 입증하였으며, 특히 BAs와 LiCoO2와 같은 물질에서 4음파 산란이 열전도에 중대한 영향을 미친다.

ABSTRACT

FourPhonon is a computational package that can calculate four-phonon scattering rates in crystals. It is built within ShengBTE framework, which is a well-recognized lattice thermal conductivity solver based on Boltzmann transport equation. An adaptive energy broadening scheme is implemented for the calculation of four-phonon scattering rates. In analogy with $thirdorder.py$ in ShengBTE, we also provide a separate python script, $Fourthorder.py$, to calculate fourth-order interatomic force-constants. The extension module preserves all the nice features of the well-recognized lattice thermal conductivity solver ShengBTE, including good parallelism and straightforward workflow. In this paper, we discuss the general theory, program design, and example calculations on Si, BAs and $\mathrm{LiCoO_2}$.

연구 동기 및 목표

• 결정에서 4음파 산란률과 열전도도를 계산하는 데 있어 첫 번째 오픈소스 도구를 개발하는 것.
• ShengBTE 프레임워크에 4음파 기능을 추가하면서도 병렬 처리와 워크플로우의 단순성을 유지하는 것.
• 정확한 4음파 산란률 계산을 위한 적응형 넓힘 기반 기법을 구현하는 것.
• 점군 대칭성과 이동 대칭성을 활용한 4계 원자간 힘 상수를 효율적으로 계산하기 위한 파이썬 스크립트( Fourthorder.py )를 제공하는 것.
• Si, BAs, LiCoO2와 같은 대표적 물질에서 방법을 검증하여 열전도에서 4음파 효과의 중요성을 입증하는 것.

제안 방법

• Boltzmann 운동방정식(BTE) 프레임워크 내에서 4음파 산란률을 계산하기 위해 ShengBTE에 새로운 모듈인 FourPhonon을 통합한다.
• 산란률의 수렴성과 정확도를 향상시키기 위해 군 속도 차이를 기반으로 한 적응형 에너지 넓힘 기법을 구현한다.
• 첫 번째 원리 계산과 대칭성 감소를 활용하여 4계 원자간 힘 상수(IfCs)를 계산하는 파이썬 스크립트인 Fourthorder.py 를 도입한다.
• 4음파 산란을 단일 모드 리프레시 타임 근사(SMRTA) 하에서 계산하고, 반복적 3음파 BTE 해법과 조합하는 수정된 워크플로우를 사용한다.
• CONTROL 파일에서 four_phonon=.TRUE.로 설정하여 4음파 계산을 활성화할 수 있도록 하며, 기존 ShengBTE 입력 및 출력 구조를 그대로 재사용한다.
• 재결합, 재분배, 분열 채널으로의 산란 과정 분해를 지원하며, Umklapp 및 정상 과정 기여도를 별도로 출력한다.

실험 결과

연구 질문

• RQ1 Si, BAs, LiCoO2 와 같은 물질에서 4음파 산란이 격자 열전도도에 얼마나 중요한가?
• RQ2 적응형 넓힘 기법이 4음파 산란률 계산의 정확도와 수렴성 향상에 효과적인가?
• RQ3 BAs와 LiCoO2 와 같이 강한 비선형성 특성을 가진 물질에서 4음파 과정이 3음파 과정에 비해 얼마나 지배적인가?
• RQ4 기존의 3음파 과정에만 기반한 접근 방식과 비교했을 때, 4음파 산란을 포함시키면 예측된 열전도도에 어떤 영향을 미치는가?
• RQ5 4음파 계산의 계산 비용과 확장성은 어떠한가? 그리고 오픈소스 도구를 통해 실용적으로 활성화할 수 있는가?

주요 결과

• 실리콘에서는 실온에서 4음파 산란률이 무시할 만큼 작으며, 허용 가능한 과정 수가 약 10^3개에 불과하여 이 시스템에서의 역할이 미미함을 확인한다.
• 붕소 arsenide(BAs)에서는 4음파 산란이 매우 강력하여 약 10^10개의 산란 과정이 가능하며, 특정 주파수 범위에서 3음파 산란을 초월하여 지배적임을 확인하였고, 이는 이전 예측 및 실험 결과와 일치한다.
• LiCoO2에서는 전체 영역의 열역학적 스펙트럼에서 4음파 산란률이 3음파 산란률과 유사한 수준이며, 특히 음향 모드 및 저주파수 옵티컬 모드에서 두드러진다.
• 재분배 채널(λ + λ′ → λ′′ + λ′′′) 이 4음파 산란에서 가장 큰 기여를 하며, LiCoO2 에서는 Umklapp 과정이 이 채널을 지배한다.
• 300 K에서 FourPhonon 은 LiCoO2 의 κ∥ = 9.35 W/mK 및 κ⊥ = 1.39 W/mK 를 예측하였으며, 문헌 값인 9.7 및 1.4 W/mK 와 양호한 일치를 보였다.
• 프로그램은 상세한 산란률, 상태 밀도, 다양한 산란 채널의 기여도를 성공적으로 계산하고 출력하여 4음파 효과의 체계적 분석을 가능하게 하였다.