558 고리 갖는 새로운 2D 탄소 동소체 열전도도: 분자 동역학 시뮬레이션 통찰

요약

이 연구는 역비평형 및 평형 분자 동역학 시뮬레이션을 사용해 새롭게 제안된 두 2차원 탄소 동소체인 팔각형-오각형 그래핀-라인(OPG-L)과 팔각형-오각형 그래핀-지그재그(OPG-Z)의 열 수송을 조사합니다. 열전도도(TC)는 시료 크기와 함께 꾸준히 증가하며, 무한 길이 외삽 시 OPG-L 310-332 W m⁻¹ K⁻¹, OPG-Z 247-228 W m⁻¹ K⁻¹ 값을 산출합니다. 그래핀(≈ 3000-5000 W m⁻¹ K⁻¹)과 비교해 낮은 TC는 감소된 포논 군속도와 평균 자유 경로에 기인합니다. 온도와 인장 변형이 체계적으로 TC를 억제해 열관리 응용을 위한 조정 가능성을 제공합니다.

서론

다이아몬드, 탄소 나노튜브, 그래핀 같은 탄소 동소체는 뛰어난 열전도도로 유명하며, 단일벽 나노튜브와 그래핀은 수천 W m⁻¹ K⁻¹에 도달합니다. 최근 이론 연구는 그래다인, 옥타그래핀, T-그래핀, H-넷 등 독특한 고리 모티프와 유망한 전자적, 기계적 특성을 가진 다양한 2D 탄소 네트워크를 공개했습니다. Su 등은 팔각형과 오각형(558 고리)으로 구성된 에너지적으로 유리하고 운동론적으로 안정한 두 2D 시트 OPG-L과 OPG-Z를 소개했습니다. 이 구조들은 각각 0.31, 0.34 eV atom⁻¹ 형성 에너지를 가져 그래다인(0.76 eV atom⁻¹)보다 상당히 낮고 이방성 전자 거동을 보입니다. 그러나 열 방산 특성은 미개척 상태로 남아 본 조사를 촉진했습니다.

모델 및 방법

OPG-L(그림 1a)과 OPG-Z(그림 1b) 원자 모델은 반복되는 558 리본에서 구축되며, 카이랄리티는 그래핀(암체어, 지그재그)과 유사하게 정의됩니다. LAMMPS를 사용해 수정된 최적화된 테르소프 포텐셜로 C-C 상호작용을 정확히 포착했습니다. 얻어진 평형 격자 상수(OPG-L: OA = 3.63 Å, OB = 9.38 Å; OPG-Z: OA = 6.78 Å, OB = 5.04 Å)는 선행 제1원리 결과와 일치합니다.

역비평형 분자 동역학(rNEMD)으로 TC를 계산했습니다. 300 K에서 0.25 ns Nosé-Hoover 평형화 후 시료를 열전달 방향 따라 50 슬래브로 분할했습니다. 차가운 슬래브에서 가장 뜨거운 원자와 뜨거운 슬래브에서 가장 차가운 원자가 운동 에너지를 교환해 안정적 열 플럭스를 생성했습니다. 결과적 온도 구배로 푸리에 법칙을 통해 TC를 얻습니다: \[ \kappa = \frac{J}{2A\,\partial T/\partial L} \] 여기서 인자 2는 양방향 열 흐름을 고려하며, 두께는 그래핀 층간 간격인 0.34 nm로 취합니다.

평형 MD(EMD) 계산은 그린-쿠보 형식으로 수행했습니다. 열 플럭스 자기상관 함수에서 진행 TC를 얻고, 1-2 ns 상관 시간 평균으로 무한 크기 값을 추출했습니다.

결과 및 논의

크기 의존성: TC는 유효 시료 길이(L_eff = L/2)와 함께 단조롭게 증가합니다. 50 nm OPG-L(지그재그)과 OPG-Z(지그재그) TC 값은 각각 125, 94 W m⁻¹ K⁻¹이며, 1000 nm에서 296, 236 W m⁻¹ K⁻¹로 상승합니다. κ⁻¹ 대 L_eff⁻¹ 선형 피팅으로 외삽된 무한 크기 TC: OPG-L 310-332 W m⁻¹ K⁻¹, OPG-Z 247-228 W m⁻¹ K⁻¹로 EMD 결과(OPG-L 313-344, OPG-Z 261-233 W m⁻¹ K⁻¹)와 우수한 일치를 보입니다.

온도 효과: 200-500 K 범위에서 움클랍 산란 강화로 TC 감소합니다. 300 K에서 500 K로 TC가 36-42% 하락합니다. 상대 전도도는 그래핀에서 관찰된 추세를 따라 2D 격자에서 포논-포논 산란 보편성을 확인합니다.

변형 효과: 열전달 방향을 따라 단축 인장 변형 적용 시 두 동소체 모두 TC가 단조롭게 억제됩니다. 10% 변형에서 TC 감소는 OPG-L 지그재그 49%, OPG-Z 암체어 31%에 도달합니다. 진동 상태 밀도 분석은 면내 모드에서 현저한 포논 연화 드러내며 열 수송 감소 설명합니다.

기계적 특성: 단축 인장 하 OPG-L과 OPG-Z는 각각 538-648 GPa, 492-550 GPa 영률을 보이며 그래다인과 비슷하나 순수 그래핀보다 낮습니다. 최대 변형률은 OPG-L 지그재그 17.2%, OPG-L 암체어 8.7%로 견고하면서도 변형 민감 거동 나타냅니다.

결론

포괄적 rNEMD와 EMD 시뮬레이션으로 OPG-L과 OPG-Z의 크기, 온도, 변형 의존 TC를 정량화했습니다. 전도도는 포논 속도와 평균 자유 경로 감소로 그래핀보다 현저히 낮으나 변형 공학으로 조정 가능합니다. 이 통찰은 OPG-L과 OPG-Z를 마이크로 및 나노스케일 소자 맞춤형 열관리 유망 후보로 확립합니다.

약어

558

오각-오각-팔각 고리

OPG-L

팔각형-오각형 그래핀-라인

OPG-Z

팔각형-오각형 그래핀-지그재그

rNEMD

역비평형 분자 동역학

TC

열전도도

VDOS

진동 상태 밀도