실제 원자 모델로 본 다결정 그래핀의 기계적 특성

초록

본 연구는 화학기상증착(CVD)으로 합성된 다결정 그래핀의 원자 규모 모델을 구축하고, 이를 이용해 파단 메커니즘과 기계적 강도를 조사한다. 결정립 경계가 교차하는 지점에서 균열이 시작되어 주로 지그재그(zigzag) 또는 팔머(armchair) 방향으로 전파되며, 내재 강도는 단결정 그래핀의 약 50% 수준으로 정규 분포를 보인다. 입자 크기가 15 nm 이상일 경우 파단 변형률은 약 0.09, 영률은 약 600 GPa로 예측된다. 이러한 강도 감소는 경계 교차점을 그리피스 균열로 보는 연속체 모델로 설명된다.

상세 분석

이 논문은 CVD 공정으로 생산되는 다결정 그래핀의 실제 미세구조를 재현하기 위해, 무작위 회전 및 위치 변화를 적용한 다중 결정립을 원자 수준에서 결합하는 알고리즘을 제시한다. 기존 연구에서는 주로 단일 결정립이나 이상적인 경계 모델만을 사용했으나, 여기서는 입자 크기(5–30 nm)와 결정립 사이의 미스오리엔테이션(0–30°)을 다양하게 설정해 30여 개의 시뮬레이션 샘플을 생성하였다. 분자동역학(MD) 시뮬레이션은 REBO(반경-결합 차수) 포텐셜을 기반으로 수행했으며, 변형률을 0.5 %/ps의 속도로 가하면서 응력‑변형률 곡선을 기록하였다.

주요 발견은 다음과 같다. 첫째, 균열은 결정립 경계가 교차하는 ‘점’에서 시작한다. 이러한 점은 원자 결함이 집중된 영역으로, 응력 집중도가 높아 파단의 시발점이 된다. 둘째, 균열 전파는 주로 그래핀 격자 내의 지그재그(zigzag) 혹은 팔머(armchair) 방향을 따라 진행한다. 이는 실험적 전자현미경 관찰과 일치하며, 결정립 내부의 강도 이방성이 아니라 경계 교차점의 구조적 취약성이 파손을 좌우함을 시사한다. 셋째, 내재 강도와 파단 변형률은 미스오리엔테이션에 거의 의존하지 않는다. 오히려 입자 크기가 작아질수록 강도가 감소하고 변형률이 증가하는 경향을 보였으며, 15 nm 이상에서는 강도가 약 50 GPa, 변형률이 0.09에 수렴한다. 넷째, 영률은 입자 크기에 따라 약 600 GPa로 일정하게 유지되며, 이는 단결정 그래핀(≈1 TPa)보다 약 40 % 낮은 값이다.

이러한 결과를 기존의 연속체 파괴 이론에 적용하면, 경계 교차점을 유효한 ‘그리피스 균열’ 길이로 모델링할 수 있다. 균열 길이가 약 2 nm 정도로 추정되며, 이때의 임계 응력은 실험값과 일치한다. 따라서 복잡한 원자 구조를 가진 다결정 그래핀도, 적절히 정의된 마이크로크랙 파라미터만으로 거시적 강도 예측이 가능함을 보여준다.