Energy dissipation at the atomic scale explains how fracture energy depends on crack velocity in silica glass

분자 동역학 시뮬레이션과 기계 학습 전위를 활용하여 본 연구는 분기 임계값 이하에서 고유 표면 에너지 밀도의 증가와 나노 규모의 거칠기 증가가 결합되어 실리카 유리의 파단 에너지가 최대 33%까지 증가함을 규명함으로써, 동적 파열이 단순히 겉보기 표면적을 증가시키는 것이 아니라 근본적으로 다른 표면 구조를 생성함을 보여준다.

핵심

두꺼운 유리 조각을 부러뜨리려 한다고 상상해 보세요. 유리를 부러뜨리는 데 필요한 에너지는 단순히 유리를 구성하는 미세한 원자 결합을 끊는 데 드는 에너지, 마치 스파게티 한 가닥을 자르는 것과 같다고 생각할 수 있습니다. 하지만 과학자들은 오랫동안 유리를 부러뜨리는 데는 그 단순한 계산이 시사하는 것보다 훨씬 더 많은 에너지가 필요하다는 것을 알고 있었습니다. 마치 유리가 저항하며 부러뜨리려면 추가적인 노력이 필요하다고 요구하는 것처럼요.

오랫동안 연구자들은 이 "추가 비용"이 주로 균열이 빠르게 진행되면서 흔들리고 고르지 않게 되어 더 거친 표면적을 만들기 때문이라고 믿었습니다 (직선 대신 고르지 않은 띠로 종이를 찢는 것과 같습니다). 하지만 최신 컴퓨터 시뮬레이션을 활용한 새로운 연구는 더 복잡한 이야기를 밝혀냈습니다.

다음은 이 논문이 발견한 바를 간단히 설명한 것입니다:

1. "초고온" 균열 선단

균열이 유리를 매우 빠르게 통과할 때, 그 균열의 선단은 극도로 뜨거워집니다. 연구에 따르면 고속에서 균열 선단 바로 앞의 원자들은 8,000 켈빈의 온도에 도달합니다 (태양의 표면보다 더 뜨겁습니다!).

균열 선단을 단순히 파손 지점으로만 생각하지 말고, 미세한 규모의 토치처럼 생각하세요. 이 intense 한 열은 유리를 녹일 뿐만 아니라, 생성되는 표면의 본질을 근본적으로 변화시킵니다.

2. 유리를 부러뜨리는 데 더 많은 비용이 드는 두 가지 이유

연구자들은 개별 원자까지 볼 수 있는 디지털 현미경과 같은 초정밀 컴퓨터 모델을 사용하여 그 추가 에너지가 어디로 가는지 파악했습니다. 그들은 "추가 비용"이 다음 두 가지 요소에 대략 50 대 50으로 나뉜다는 것을 발견했습니다:

• "거칠기" 요인 (양): 균열이 가속화됨에 따라 그 뒤에 남는 표면은 완벽하게 매끄럽지 않습니다. 우주에서 본 산맥처럼 나노 단위로 거칠어집니다. 이는 균열이 외부에서 보이는 것보다 실제로 더 많은 표면적을 생성한다는 것을 의미합니다.
  • 비유: 천을 찢는다고 상상해 보세요. 천을 천천히 찢으면 가장자리가 곧습니다. 하지만 빠르게 찢으면 가장자리가 찢어지고 고르지 않게 됩니다. 그 고르지 않은 가장자리를 만드는 데 더 많은 천을 사용한 것입니다.
• "품질" 요인 (에너지 밀도): 이것이 새로운 발견입니다. 만약 그 고르지 않은 표면을 매끄럽게 다듬더라도, 차분하고 느리게 생성된 표면보다 여전히 더 많은 에너지가 필요합니다. 균열 선단의 극심한 열은 새로운 표면의 원자 구조를 변화시켜 더 높은 에너지를 갖거나 더 불안정하게 만듭니다.
  • 비유: 케이크를 굽는다고 상상해 보세요. 천천히 구운 케이크는 표준적인 질감을 가집니다. 하지만 토치로 강하게 구우면 바깥쪽이 타서 화학적으로 달라집니다. "탄" 표면은 근본적으로 다르며, 매끄럽고 천천히 구운 버전보다 생성하는 데 더 많은 에너지를 필요로 합니다.

3. "숨겨진" 거칠기

가장 흥미로운 점 중 하나는 컴퓨터가 발견한 "거칠기"가 원자 규모로 너무 작아서 공학자들이 깨진 유리를 측정하는 데 사용하는 표준 도구로는 완전히 놓쳐버린다는 것입니다.

일반 현미경으로 깨진 유리 조각을 보면 매끄러운 표면이 보일 것입니다. 모든 추가 에너지가 표면을 "더 뜨겁게" 또는 더 에너지가 많게 만드는 데 들어갔다고 가정할 것입니다. 하지만 이 연구는 그 에너지의 상당 부분이 실제로 표면을 물리적으로 더 크고 거칠게 만드는 데 들어갔음을 보여줍니다. 다만 그 규모가 우리 눈에는 너무 작을 뿐입니다.

4. 균열 이동 속도에 대한 수학 공식 수정

이 논문은 가해진 힘에 따라 균열이 얼마나 빠르게 이동하는지 예측하는 데 오랫동안 사용되어 온 공식을 바로잡았습니다. 기존 공식 ("프룬드 모델") 은 고속에서 약간 흐릿해진 지도와 같았습니다. 새로운 연구는 데이터와 완벽하게 부합하는 더 나은 공식 ("제곱근 관계") 을 발견했습니다.

이 수정은 중요합니다. 왜냐하면 이전의 깨진 유리의 열을 측정했던 실험들 (균열에서 방출되는 빛, 즉 파열발광을 사용) 이 속도 예측과 완전히 일치하지 않았던 이유를 설명하는 데 도움이 되기 때문입니다. 새로운 공식을 사용하면 예측된 속도와 온도가 마침내 컴퓨터 시뮬레이션이 보여준 결과와 일치하게 됩니다.

결론

유리를 부러뜨리는 것은 단순히 결합을 끊는 것만이 아닙니다. 균열이 빠르게 이동할 때, 그것은 다음과 같은 역할을 하는 작은 초고온 레이저처럼 행동합니다:

1. 표면을 물리적으로 더 거칠게 만들어 (더 많은 면적을 생성합니다).
2. 표면을 화학적으로 변형시켜 더 에너지가 많게 만듭니다.

이 연구는 유리를 부러뜨리는 데 필요한 에너지가 고정된 숫자가 아니라, 부러뜨리는 속도에 따라 변하며, 파손의 형태와 선단에서의 극심한 열 모두에 의해 주도된다는 것을 증명합니다.

기술 요약

문제 제기
실리카 유리 같은 취성 재료의 파괴 에너지는 균열 속도가 증가함에 따라 증가하는 것으로 관찰되며, 이는 결합 파괴의 열역학적 비용을 초과하는 현상입니다. 이러한 속도 의존성은 종종 동적 불안정성 (예: 미세 가지 형성) 으로 인한 표면 거칠기 증가로 귀결되지만, 가지 형성 임계값 이하에서의 거동은 여전히 잘 이해되지 않고 있습니다. 고분자 (예: PMMA) 에 대한 이전 연구들은 가지 형성 시작 전 파괴 에너지가 약 30% 증가하는 것을 관찰했으며, 이를 해석되지 않은 소산 과정으로 귀인했습니다. 핵심적인 질문은 여전히 남아 있습니다: 가지 형성 전 에너지 증가는 (표준 측정으로는 보이지 않는 나노 규모의 거칠기를 통해) 표면적의 양이 증가했기 때문인지, 아니면 균열 선단에서의 극한 조건으로 인해 생성된 표면의 질 (고유 에너지 밀도) 이 증가했기 때문인지? 또한, 이전의 분자 동역학 (MD) 연구에서 사용된 경험적 원자 간 퍼텐셜은 고에너지 결합 파괴 사건을 정확하게 포착하지 못해 가지 형성 속도를 과소평가하고 과정 영역 (process zone) 설명을 부정확하게 만들었습니다.

방법론
이러한 한계를 극복하기 위해 저자들은 융합 실리카 유리에 대한 광범위한 1 차 원리 분자 동역학 시뮬레이션을 수행했습니다.

• 원자 간 퍼텐셜: 경험적 퍼텐셜 대신, Allegro 아키텍처 기반의 머신 러닝 원자 간 퍼텐셜 (MLIP) 을 사용했습니다. 이 모델은 용융 실리카, 유리질 구조, 그리고 15,000 K 까지 도달하는 기계적 변형을 포함한 광범위한 상태를 포괄하는 r2SCAN 함수를 사용한 밀도 범함수 이론 (DFT) 데이터로 훈련되었습니다.
• 시뮬레이션 설정: 시스템은 사전에 존재하는 노치를 가진 1,230 만 개의 원자로 구성된 실리카 유리 스트립 (300 × 8.5 × 75 nm³) 으로 구성되었습니다. 시뮬레이션은 300 K 와 1 bar 에서 수행된 후, 균열 전파를 유도하기 위해 급격한 인장 변위가 가해졌습니다.
• 매개변수: 변위 속도를 달리하여 14 개의 시뮬레이션을 수행했으며, 그 결과 초기 인장 응력은 3.9 GPa 에서 7 GPa 사이였습니다.
• 분석: 저자들은 "구조적 파괴 에너지" (냉각 후 손상된 구조에 저장된 에너지를 투영 면적으로 정규화한 값) 와 "고유 표면 에너지 밀도" (델로네 삼각분할을 통해 측정한 실제 거친 표면적으로 정규화한 에너지) 를 구분했습니다. 또한 0.5 ps 간격으로 운동 에너지를 평균화하여 균열 선단 온도를 분석했습니다.

주요 결과

1. 속도 의존적 파괴 에너지: 가지 형성 임계값 (레이리 파 속도 $v_R$의 0.72 로 정의됨) 이하에서도 균열 속도가 증가함에 따라 구조적 파괴 에너지는 최대 33% 까지 상승합니다.
2. 에너지 증가의 분해: 이 상승은 두 가지 메커니즘 사이에 대략 균등하게 분배됩니다:
   • 나노 규모 거칠기: 투영 면적에 비해 실제 파괴 표면적 ($A_{real}$) 이 약 16% 증가합니다. 이 거칠기는 균열 주변에 더 두꺼운 손상 영역을 생성하지만, 표준 사후 광학 분석으로는 보이지 않는 길이 척도에서 발생합니다.
   • 고유 표면 에너지 밀도: 고유 표면 에너지 밀도 ($\gamma_s$) 가 약 15% 증가합니다. 이는 과정 영역이 불과 3~6 nm 일 때 균열 선단에서 온도가 3,000 K 를 초과하는 극한의 열적 조건에 기인합니다.
3. 균열 속도와 불안정성: 시뮬레이션은 균열 속도가 응력 강도 인자와 제곱근 관계 ($v_c \propto \sqrt{1 - (K_{Ic}/K_I)^2}$) 를 따르며, 표준 Freund 모델과 다르다는 것을 보여줍니다. 가지 형성은 약 2.4 km/s (0.72 $v_R$) 에서 발생하며, 이는 융합 실리카의 실험적 관찰과 일치합니다.
4. 과정 영역 크기: 과정 영역 반지름은 응력 강도에 따라 증가하며, 평면 균열의 경우 약 3.1 nm 에서 가지 형성 균열의 경우 약 6.8 nm 까지 다양합니다. 이 크기는 Dugdale–Barenblatt 예측과 일치하며, 이전 경험적 퍼텐셜 시뮬레이션에서 유도된 값 (~10 nm) 보다 작습니다.
5. 열적 불일치 해결: 이 연구는 시뮬레이션 온도와 프랙토루미네선스 기반 실험 예측 사이의 불일치를 해결합니다. 새로 확인된 제곱근 속도–응력 관계를 적용함으로써, 주어진 속도에 대한 예측된 균열 선단 온도는 더 낮아져 분석 모델이 시뮬레이션 데이터와 더 잘 일치하게 됩니다.

의의 및 주장
본 논문은 실리카 유리의 동적 파괴가 단순히 더 많은 겉보기 표면적을 생성함으로써 파괴 에너지를 증가시키는 것이 아니라, 나노 규모에서 표면을 근본적으로 변화시킴으로써 파괴 에너지를 증가시킨다고 주장합니다. 저자들은 다음과 같은 점을 입증합니다:

• 파괴 에너지는 일정한 재료 특성이 아닙니다. 거시적 가지 형성이 없더라도 속도에 의존적입니다.
• 표준 사후 분석은 나노 규모 거칠기의 기여와 증가된 고유 에너지 밀도의 기여를 혼동하여, 전체 증가분을 표면 에너지 밀도의 상승으로 해석할 것입니다.
• MLIP 를 사용한 1 차 원리 MD 시뮬레이션은 이제 경험적 튜닝 없이 동적 파괴를 정량적으로 정확하게 모델링할 수 있으며, 실험 데이터와 일치하는 탄성 특성, 파괴 인성, 그리고 가지 형성 임계값을 성공적으로 재현합니다.
• 수정된 속도–응력 관계는 프랙토루미네선스 실험에서 추론된 균열 선단 온도를 해석하기 위한 더 정확한 틀을 제공합니다.

이 연구는 PMMA 와 같은 재료에서 이전에 관찰되었으나 설명되지 않았던 가지 형성 전 파괴 에너지의 증가가 나노 규모 거칠기와 극한 국부 온도에 의해 주도되는 상승된 고유 표면 에너지 밀도의 조합에서 비롯된다고 결론지었습니다.