[논문 리뷰] Ab initio simulation of amorphous graphite

Rajendra Thapa, Chinonso Ugwumadu외 3명|PubMed|2022. 02. 22.

Graphene research and applications인용 수 4

한 줄 요약

이 연구는 밀도함수이론기반 분자역학 시뮬레이션을 사용하여, 불규칙한 단일층 그래핀 평면(오각형, 육각형, 칠각형 포함)으로 구성된 다층 탄소 물질인 비정질 그래파이트(a-G)의 형성을 입증한다. 전이 온도는 약 3000 K이며 밀도는 2.2–2.8 g/cm³ 범위에서 발생하며, 이질적 π-전자 분포에 의한 이질적 이중결합 공유결합이 반데르발스 힘보다 주요 상호작용으로 작용하고, 구조적 고리 결함으로 인해 전자적 전도도가著 감소한다.

ABSTRACT

An amorphous graphite material has been predicted from molecular dynamics simulation using ab initio methods. Carbon materials reveal a strong proclivity to convert into a sp^{2} network and then layer at temperatures near 3000 K within a density range of ca. 2.2-2.8 g/cm^{3}. Each layer of amorphous graphite is a monolayer of amorphous graphene including pentagons and heptagons in addition to hexagons, and the planes are separated by about 3.1 Å. The layering transition has been studied using various structural and dynamical analyses. The transition is unique as one of partial ordering (long range order of planes and galleries, but topological disorder in the planes). The planes are quite flat, even though monolayer amorphous graphene puckers near pentagonal sites. Interplane cohesion is due partly to non-Van der Waals interactions. The structural disorder has been studied closely, especially the consequences of disorder to electronic transport. It is expected that the transition elucidated here may be salient to other layered materials.

연구 동기 및 목표

• 고온에서의 열처리 조건에서 비정질 탄소 또는 무작위 구조에서 새로운 비정질 그래파이트(a-G) 상의 형성 메커니즘을 규명하는 것.
• 특히 다층 구조와 계면 간 응집력에 초점을 맞춰 a-G의 구조적, 에너제틱적, 전자적 성질을 규명하는 것.
• 오각형, 칠각형 등 구조적 고리 결함이 평면 내 전자적 전도도 및 전하 분포에 미치는 영향을 분석하는 것.
• 반데르발스 힘을 초월한 전자 상호작용이 a-G의 다층 구조 안정화에 미치는 역할을 밝히는 것.
• 비정질 그래파이트와 결정성 그래파이트, 결함이 있는 그래핀의 전자 구조 및 전도도를 비교하는 것.

제안 방법

• 3000 K에서 Nosé-Hoover 열역학장치를 사용한 PBE 함수를 적용한 ab initio 분자역학(AIMD) 시뮬레이션을 수행하였으며, 탄소 시스템의 밀도는 2.2–2.8 g/cm³ 범위였다.
• 다양한 시스템 크기, 初기 상태(비정질 또는 무작위), 스핀-오르빗 상호작용 함수(예: PBE, PBE+vdW, GAP-ML)를 갖는 M1–M6 모델에 대해 NVT 시뮬레이션을 실시하였다.
• 시뮬레이션 후 구조적 안정성 확인 및 에너지 최소화를 위해 공액 기울기 최적화를 적용하였다.
• 전자 구조 분석에는 밀도 상태(Dos), 전하 밀도 맵핑, 그리고 밴드 분해 전하 밀도 투영을 통해 π 및 π* 오비탈 기여도를 시각화하였다.
• Kubo-Greenwood 공식에 기반한 공간 투영 전도도(SPC)를 계산하여 전도 경로를 평가하고 전자 전도도 이방성을 정량화하였다.
• 구조 분석에는 상관 함수(RDF), 원자 위치 군집 기반 층 식별, 시간에 따른 원자 구조 시각화를 포함하였다.

실험 결과

연구 질문

• RQ1 비정질 탄소 또는 무작위 탄소 구조가 고온에서 열처리를 거치면 자발적으로 그래파이트 유사 다층 구조를 형성할 수 있는가?
• RQ2 형성된 비정질 그래파이트의 구조적 및 전자적 성질에 대해 오각형, 칠각형 등 고리 결함이 어떤 역할을 하는가?
• RQ3 a-G에서 계면 간 응집력은 반데르발스 힘과 어떻게 다를 수 있으며, 어떤 전자적 메커니즘이 다층 구조를 안정화시키는가?
• RQ4 고리 결함으로 인해 비정질 그래파이트의 평면 내 전자 전도도가 결정성 그래파이트에 비해 얼마나 감소하는가?
• RQ5 a-G의 전자 구조는 계면 간 간극 내에 확산된 π-전자에 의해 특징지어질 수 있으며, 이는 전하 분포에 어떤 영향을 미치는가?

주요 결과

• 비정질 그래파이트(a-G)는 약 3000 K와 2.2–2.8 g/cm³의 온도/밀도 창에서 안정적으로 형성되며, 불규칙한 단일층 그래핀 평면이 약 3.1 Å 간격으로 다층으로 배열된다.
• a-G의 총 에너지는 결정성 그래파이트보다 0.32 eV/atom 높지만, 구조적 고리 불규칙성에도 불구하고 열역학적으로 안정된 상태임을 시사한다.
• 계면 간 응집력은 반데르발스 힘보다는 확산된 π-전자들이 간극 내로 연장되는 결합 오비탈에 의해 주로 결정된다.
• 전하 밀도 맵은 간극 내에 낮은 밀도의 확산 전자 기체를 보이며, 고립된 평면에서 최대 전하 밀도의 2% 이상이 관찰되어 상당한 전자 결합이 존재함을 나타낸다.
• a-G의 전자 밀도 상태(Dos)는 패러미 레벨 근처에 넓은 피크를 보이며, 결정성 그래파이트의 반금속 특성과는 다릅니다.
• 공간 투영 전도도(SPC) 분석 결과, a-G의 전자 전도도는 그래파이트 대비 약 10⁻² 배 감소하였으며, 전도 경로는 오각형 및 칠각형 등의 고리 결함을 피하는 경향을 보였다.