[논문 리뷰] Direct Calculation of Modal Contributions to Thermal Conductivity via Green-Kubo Modal Analysis: Crystalline and Amorphous Silicon

Wei Lv, Asegun Henry|arXiv (Cornell University)|2015. 03. 27.

Thermal properties of materials참고 문헌 53인용 수 155

한 줄 요약

이 논문은 Green-Kubo 모드 분석(GKMA)을 소개한다. 이는 음향파 속도나 음향가스 모델에 의존하지 않고 직접적으로 열전도도에 대한 모드 기여를 계산하는 제1원리 방법이다. 격자역학과 Green-Kubo 형식을 결합함으로써 GKMA는 완전한 비조화성과 모드-모드 상관관계를 포착하여, 결정성과 비정질 실리콘 모두에서 온도 의존성 열전도도를 정확하게 예측할 수 있다. 특히 양자 보정을 적용했을 때 비정질 실리콘의 실험 데이터와 가장 우수한 일치를 보이며, 지금까지의 보고된 결과들 중에서 가장 뛰어난 성능을 보인다.

ABSTRACT

In this letter we derive a new method for direct calculation of the modal contributions to thermal conductivity, which is termed Green-Kubo modal analysis (GKMA). The GKMA method combines the lattice dynamics formalism with the Green-Kubo formula for thermal conductivity, such that the thermal conductivity becomes a direct summation of modal contributions, where one need not define the phonon velocity. As a result the GKMA method can be applied to any material/group of atoms where the atoms vibrate around stable equilibrium positions, which includes not only crystalline line compounds, but also random alloys, amorphous materials and even molecules. The GKMA method provides new insight into the nature of phonon transport, as it casts the problem in terms of mode-mode correlation instead of scattering, and provides a general unified formalism that can be used to understand phonon-phonon interactions in essentially any class of materials or structures where the atoms vibrate around stable equilibrium sites.

연구 동기 및 목표

• 음향파 속도나 음향가스 모델에 의존하지 않는 모드 기여를 계산하기 위한 일반적인 형식을 개발하는 것.
• 비정기적 물질인 비정질 실리콘과 같이 음향모드가 잘 정의된 평면파가 아닌 시스템에서 열전도를 기술하는 데에 한계가 있는 기존 모델의 문제점을 해결하는 것.
• 기존 접근 방식이 실패하는 시스템에서도 완전한 비조화성과 모드-모드 상관관계를 열전도도 계산에 통합하는 것.
• 고전적 분자 동역학와 모드별 열용량에 대한 양자 보정을 조합하여 온도 의존성 열전도도를 정확하게 예측할 수 있도록 하는 것.

제안 방법

• 해당 방법은 변환 행렬을 통해 원자 이동과 속도를 정규 모드 기여의 초위상으로 표현하기 위해 격자역학을 사용한다.
• 모드 속도 성분에 기반하여 하디 열류 연산자를 사용하여 각 모드의 시간에 따른 열류 기여를 유도한다.
• 모드 해상도 열류에 대해 Green-Kubo 공식을 적용함으로써 각 모드의 열전도도에 대한 직접 계산이 가능해진다.
• 모드-모드 상관관계는 시간 상관 함수를 통해 계산되며, 서로 다른 주파수의 모드 간의 상호작용을 드러낸다.
• 각 모드의 기여를 양자와 고전적 열용량의 비율을 사용하여 스케일링함으로써 고전적 MD 결과에 대한 양자 보정이 가능하다.
• 이 방법은 결정성 및 비정질 실리콘에 적용하여 실험 데이터 및 이전 이론 모델과의 비교를 통해 검증되었다.

실험 결과

연구 질문

• RQ1 음향파 속도나 군속도가 잘 정의되어 있지 않은 상황에서도 열전도도에 대한 모드 기여를 직접 계산할 수 있는가?
• RQ2 비정질 물질인 실리콘과 같은 시스템에서, 특히 비대칭 항목(비대칭 항)이 포함된 모드-모드 상관관계는 열전도도에 어떻게 기여하는가?
• RQ3 분자 동역학를 통해 포착된 완전한 비조화성은 조화 근사와 비교해 볼 때 열전도도 예측을 얼마나 향상시키는가?
• RQ4 GKMA 프레임워크에서 고전적 MD 결과에 양자 보정을 적용하면 비정질 실리콘의 온도 의존성 열전도도를 정확하게 예측할 수 있는가?
• RQ5 GKMA 형식은 국소화된(locon) 모드와 확산되는(propagon/diffuson) 모드의 열전도도 기여에 대한 새로운 물리적 통찰을 제공하는가?

주요 결과

• GKMA 방법은 음향파 속도가 필요 없이도 모드 기여를 계산할 수 있어, 비정기적 시스템인 비정질 물질에까지 적용 가능하다.
• 비정질 실리콘의 경우, 양자 보정을 적용한 GKMA 방법이 지금까지 보고된 바 중에서 실험 데이터와 가장 우수한 일치를 보였다.
• 실온에서 모드-모드 상관관계는 총 열전도도의 약 30% 기여를 하며, 비조화성이 중요한 역할을 함을 시사한다.
• 상관관계 행렬은 약 16 THz에서 모드 상호작용의 급격한 전이가 나타나며, 이는 국소화된 모드(locons)의 출현과 일치한다. 이는 모드 성질에 대한 사전 지식 없이도 관측 가능하다.
• locon 모드의 자기상관관계는 무시할 만큼 작으며, 이는 열전도도에 미치는 기여가 미미함을 일관되게 반영한다.
• 저온 영역에서 모드 기여의 이동을 정확히 포착하며, GKMA로 확보한 주파수 해상도 기여에 대해 온도 의존성을 정확히 기대하기 위해서는 양자 보정이 반드시 필요하다는 점을 확인했다.