[논문 리뷰] Understanding Phonon Transport Properties Using Classical Molecular Dynamics Simulations

Murali Gopal Muraleedharan, Kiarash Gordiz외 5명|arXiv (Cornell University)|2020. 11. 02.

Thermal properties of materials참고 문헌 373인용 수 6

한 줄 요약

이 논문은 다양한 재료에서 열전도도 및 계면 열전도도와 같은 진동수 전도성 특성을 예측하기 위한 고전적 분자역학(MD) 시뮬레이션에 대한 종합적인 리뷰를 제시한다. MD와 격자역학을 융합함으로써 고전적 근사와 부정확한 원자 간 잠재에너지 함수의 한계에도 불구하고 진동 모드 기여도를 세밀하게 분석할 수 있으며, 잠재에너지 함수의 정확성 향상 및 양자 보정 전략을 강조한다.

ABSTRACT

Predictive modeling of the phonon/thermal transport properties of materials is vital to rational design for a diverse spectrum of engineering applications. Classical Molecular Dynamics (MD) simulations serve as a tool to simulate the time evolution of the atomic level system dynamics and enable calculation of thermal transport properties for a wide range of materials, from perfect periodic crystals to systems with strong structural and compositional disorder, as well as their interfaces. Although MD does not intrinsically rely on a plane wave-like phonon description, when coupled with lattice dynamics calculations, it can give insights to the vibrational mode level contributions to thermal transport, which includes plane-wave like modes as well as others, rendering the approach versatile and powerful. On the other hand, several deficiencies including the lack of vibrationally accurate interatomic potentials and the inability to rigorously include the quantum nature of phonons prohibit the widespread applicability and reliability of Molecular Dynamics simulations. This article provides a comprehensive review of classical Molecular Dynamics based formalisms for extracting thermal transport properties: thermal conductivity and thermal interfacial conductance and the effects of various structural, compositional, and chemical parameters on these properties. Here, we highlight unusual property predictions, and emphasize on the needs and strategies for developing accurate interatomic potentials and rigorous quantum correction schemes.

연구 동기 및 목표

• 구조적 및 조성적 불순물 정도가 다른 재료에서 열전도성 특성을 예측하는 데 있어 고전적 분자역학(MD)을 체계적으로 활용하는 방법을 개괄하는 것.
• 고전적 MD의 핵심적 한계, 즉 진동에 기반한 정확한 원자 간 잠재에너지 함수의 부재와 진동수의 고전적 처리로 인한 예측 정확도 저하 문제를 규명하고 해결하는 것.
• MD 기반 열전도성 예측의 신뢰성을 향상시키기 위해 정확한 원자 간 잠재에너지 함수 개발과 철저한 양자 보정 기법의 중요성을 강조하는 것.
• MD와 격자역학을 융합하여 평면파 유사 모드와 비평면파 모드를 포함한 다양한 진동 모드가 열전도성에 기여하는 방식을 분석하는 것.
• 공학적 응용을 위한 열전도성 특성이 조절된 재료 설계에 효과적으로 MD 시뮬레이션을 적용할 수 있도록 연구자들을 안내하는 것.

제안 방법

• 원자 시스템의 시간적 진화를 계산하고 열전도도 및 계면 열전도도와 같은 열전도성 특성을 추출하기 위해 고전적 분자역학(MD) 시뮬레이션을 활용한다.
• MD와 격자역학 계산을 융합하여 열전도성의 기여도를 평면파 유사 모드와 비평면파 모드를 포함한 개별 진동 모드로 분해한다.
• 균형 상태 MD 시뮬레이션을 기반으로 열전도도를 계산하기 위해 그린-쿠보 공식을 적용한다. 이는 열류 자기상관함수를 통해 수행된다.
• 계면에 온도 기울기를 도입함으로써 계면 열전도도를 계산하기 위해 역방향 비균형 분자역학(RNMD) 방법을 사용한다.
• MD 시뮬레이션을 통해 구조적 결함, 조성 불순물, 화학적 이질성이 진동수 산산각과 열전도성에 미치는 영향을 분석한다.
• 기계학습 기반 접근법과 제1원리 데이터에 적합한 방법을 통해 원자 간 잠재에너지 함수의 정확성을 향상시키는 전략을 제안하고, 제로점 운동과 양자화된 진동수 통계를 고려하기 위한 양자 보정 기법을 논의한다.

실험 결과

연구 질문

• RQ1 완벽한 결정에서 불순물이 있는 시스템에 이르기까지 다양한 재료의 열전도도를 효과적으로 예측하기 위해 고전적 분자역학 시뮬레이션을 어떻게 활용할 수 있는가?
• RQ2 고전적 MD는 진동수의 양자적 성격을 어떻게 반영하지 못하며, 이러한 한계는 어떻게 보완할 수 있는가?
• RQ3 구조적 및 조성적 불순물은 진동수 전도성에 어떤 영향을 미치며, MD 시뮬레이션은 원자 수준에서 이를 어떻게 밝혀낼 수 있는가?
• RQ4 MD와 격자역학을 융합할 경우, 특정 진동 모드의 기여도가 열전도성에 얼마나 정밀하게 분리될 수 있는가?
• RQ5 열전도성 예측 정확도를 향상시키기 위해 가장 효과적인 정확한 원자 간 잠재에너지 함수 및 양자 보정 전략은 무엇인가?

주요 결과

• 고전적 MD 시뮬레이션은 격자역학과 융합함으로써 비평면파 모드를 포함한 개별 진동 모드의 기여도를 해석할 수 있으며, 기존의 평면파 근사 이론을 뛰어넘는 세밀한 그림을 제공한다.
• 균형 상태 MD에서 적용된 그린-쿠보 방법은 주기적이고 결함이 없는 시스템에서 열전도도를 정확하게 계산할 수 있으며, 제1원리 계산 결과와 일치하는 결과를 도출한다.
• 역방향 비균형 분자역학(RNMD)은 특히 초합성 구조나 비정질층이 포함된 복잡한 계면에서 계면 열전도도 계산에 효과적이다.
• 구조적 및 조성적 불순물은 진동수 산산각을 강화시켜 열전도도를 크게 감소시키며, MD 시뮬레이션은 이러한 영향을 원자 척도에서 정량화할 수 있다.
• MD 예측의 정확도는 원자 간 잠재에너지 함수의 품질에 매우 민감하며, 부정확한 잠재에너지 함수는 물리적으로 비합리적인 진동수 분산 및 열전도도 값을 초래할 수 있다.
• 제로점 운동을 고려하고 실험 데이터와의 일치도를 향상시키기 위해, 드바이 모델 기반 또는 제1원리에서 유도된 보정 기법과 같은 양자 보정 기법이 필수적이다.