유리 파손 메커니즘: 손상 영역이 표면·파괴 에너지 차이를 만든다

초록

**
본 연구는 다중 스케일 시뮬레이션을 통해 비정질 실리카 유리의 자유표면에너지와 파괴에너지 사이의 차이를 규명한다. 원자 규모 MD 모델을 연속체 위상장 모델에 동화시켜 손상 확산 영역(≈16–23 Å)을 정량화하고, 표면 에너지와 손상 에너지의 기여를 명확히 구분한다. 플라스틱 변형은 무시할 수 있을 정도로 작으며, 손상은 주로 고리 구조 변화와 실리콘 원자 배위수 변동에 기인한다. 결과는 실리카 파괴인성 측정값과 일치한다.

**

상세 분석

**
이 논문은 비정질 실리카의 파괴 거동을 원자 수준과 연속체 수준을 연결하는 ‘멀티스케일’ 접근법으로 분석한다. 원자 간 상호작용을 BKS 포텐셜로 기술한 분자동역학(MD) 시뮬레이션을 수행하고, 얻어진 원자 배열을 가우시안 컨볼루션(코스 그레이닝)으로 연속체 변수(변위, 변형률, 응력, 포텐셜 에너지 등)로 변환한다. 이렇게 얻은 연속체 필드에 대해 손상 변수 d를 정의하고, 손상 에너지 ψ_el 과 무손상 탄성 에너지 ψ₀ 의 차이로 손상 분포를 계산한다. 손상은 압축에서는 발생하지 않으며, d=0은 무손상, d=1은 완전 파손을 의미한다.

핵심 결과는 두 가지 에너지 기여가 명확히 구분된다는 점이다. 첫 번째는 ‘자유표면에너지’(γ)로, 원자 결합이 끊어지고 표면 원자들이 재배열될 때 발생하는 잠재에너지 증가이다. 이는 직접적인 표면 절단 실험(에너지 최소화 후 2γ = ΔΨ_pot / (Lx Lz))을 통해 정량화되었다. 두 번째는 ‘손상 확산 에너지’로, 균열 경로 주변 16–23 Å 범위에 걸쳐 구조적 재배열(고리 구조 변형, Si 배위수 변화)이 일어나면서 저장된 탄성 에너지가 소멸한다. 손상 영역의 폭(l)은 가우시안 피팅을 통해 추정되었으며, 위상장 모델의 내부 길이 스케일 l_c와 일치한다.

플라스틱 변형은 전통적인 금속이나 고분자와 달리 거의 관찰되지 않았으며, 손상 에너지의 대부분이 비탄성적인 구조 재배열에 기인한다. 이는 실리카가 ‘극한 취성(brittle)’ 재료임을 재확인하면서도, 기존의 ‘표면에너지 = 파괴에너지/2’ 가정이 실리카에서는 성립하지 않음을 설명한다.

또한, 위상장 파괴 모델에 필요한 파괴인성 g_c와 내부 길이 l_c를 원자 시뮬레이션 결과와 최소오차 기준으로 피팅하였다. FEMU(Finite Element Update) 절차를 통해 20개의 단면에서 손상 필드를 재현하고, 단일 스냅샷만으로 균열 길이, g_c, l_c를 추정할 수 있음을 보여준다. 이는 실험적으로 측정하기 어려운 내부 손상 길이를 간접적으로 파악하는 새로운 방법론으로 평가된다.

결론적으로, 이 연구는 비정질 실리카 파괴에서 표면 에너지와 손상 확산 에너지의 비율이 약 1:4 정도이며, 이는 실험적으로 보고된 파괴인성(g_c ≈ 5γ)과 일치한다는 점을 입증한다. 손상 영역은 고리 구조와 배위수 변화를 통해 정의되며, 이는 향후 유리 재료 설계와 파괴 예측 모델에 중요한 물리적 인사이트를 제공한다.

**