다층 탄소 나노구조의 sp³ 결합이 비선형 변형·파괴에 미치는 영향: 다중스케일 보조노드(MAN) 모델

초록

본 논문은 원자 수준의 REBO 포텐셜을 고차 연속성의 연속 시트에 매핑하는 다중스케일 보조노드(MAN) 방법을 제안한다. 이동 최소제곱(MLS) 보간을 통해 원자 변위를 노드 변위로 변환하고, 비선형 강성 방정식을 반복적으로 풀어 2D 탄소 나노구조(이중층 그래핀, 다이아메인 및 중간 전이 구조)의 대변형·파괴 거동을 예측한다. 결과는 sp³ 결합의 밀도·분포가 인장 강도와 파괴 패턴을 크게 좌우함을 보여준다.

상세 분석

MAN 방법은 두 개의 가상 연속 시트를 구축하고, 각 원자에 영향 영역을 정의한 뒤 MLS 근사를 이용해 원자 위치와 노드 변위 사이의 고차 연속 보간함수를 만든다. 이렇게 얻은 보간함수는 원자 간 결합 거리와 각도를 REBO 포텐셜에 입력해 계산된 원자 포텐셜 에너지를 연속 시트의 변형 에너지 형태로 변환한다. 비선형 최소화 과정에서 전체 시스템의 자유도는 노드 수에 의해 조절 가능하므로, 계산 비용을 크게 절감하면서도 원자 수준의 정확도를 유지한다.

연구에서는 sp³ 결합이 전부 존재하는 완전 다이아메인부터, 일부만 sp³으로 전이된 구조까지 다양한 모델을 구축하였다. 시뮬레이션 결과, 다이아메인의 Young’s modulus와 전단계수는 각각 1.1 TPa, 0.44 TPa 수준으로 다이아몬드에 근접하고, 그래핀(≈1.0 TPa, 0.4 TPa)보다 현저히 높았다. 특히, sp³ 결합이 균일하게 분포될 경우 인장 파괴 전 응력‑변형 곡선이 급격히 상승하고 파괴 인장이 30 % 이상 증가한다. 반면, sp³ 결합이 국부적으로 집중되면 응력 집중이 발생해 파괴가 초기화되며, 전체 강도 향상 효과가 제한적이다.

MAN 모델은 전단·압축·전단-압축 복합 하중까지 확장 가능하나, 현재 연구는 인장 하중에만 초점을 맞추었다. 또한, 층간 반데르발스 상호작용을 무시했으며, REBO 포텐셜의 절단 함수가 1.7–2.0 Å 구간에서 비물리적 힘 증가를 일으킬 수 있다는 한계도 명시하였다. 이러한 제약에도 불구하고, MAN은 원자‑연속계 사이의 에너지 매핑을 체계적으로 구현함으로써, 대규모 2D 탄소 구조의 비선형 거동과 파괴 메커니즘을 효율적으로 탐구할 수 있는 강력한 도구임을 입증한다.