실리콘 레이저 절삭의 새로운 비열역학적 기작

초록

본 논문은 초고속 단일 펄스 레이저에 의해 강하게 여기된 실리콘의 비열역학적 거동을 조사한다. 열스파이크 모델(TSM)과 분자동역학(MD)을 결합하고, 여기 수준에 따라 변하는 결합 강도를 반영한 수정된 Tersoff 포텐셜(MOD*)을 적용해 대규모 병렬 시뮬레이션을 수행하였다. 결과적으로 비열융해, 비열증발, 전자‑구멍 구속에 의한 열전도 변화, 전자‑양성자 폭발 등 새로운 비열절삭 메커니즘을 제시하고, 전자 온도에 따른 실리콘의 상도표를 열역학적 적분법으로 구축하였다.

상세 분석

이 연구는 실리콘과 같은 공유 결합 반도체가 금속과 달리 레이저 광자에 의해 전자-정공쌍이 급격히 생성되고, 이들 캐리어가 결합 전자를 반결합(anti‑bonding) 상태로 전이시키면서 원자간 결합 강도가 실시간으로 약화된다는 점에 주목한다. 이를 정량화하기 위해 저자들은 기존의 Tersoff 포텐셜을 전자 온도(또는 여기 정도) 의존적으로 재파라미터화한 MOD*를 도입했으며, 이는 고전 MD 시뮬레이션에 비열 효과를 직접 삽입할 수 있게 한다. 열스파이크 모델(TSM)은 전자·양자 서브시스템의 에너지 전달을 2‑step(전자‑포논 결합)으로 기술하고, FD 격자와 MD 격자를 다중 CPU에 분산시켜 수십 억 원자 규모의 시뮬레이션을 가능하게 한다.

시뮬레이션 결과는 네 가지 주요 비열절삭 메커니즘을 밝혀낸다. 첫째, 비열증발은 전자 여기가 반결합 상태를 지배해 원자간 포텐셜이 순전히 반발적으로 변하면서 피코초 이하 시간에 물질이 급격히 기화한다. 둘째, 전‑충격파 비열융해는 전자 여기가 국소 압력을 급격히 상승시켜 광속에 근접하는 융해 전선이 전파되는 현상으로, 전통적인 음속 제한을 초월한다. 셋째, 비열공극 형성은 비열융해 과정에서 다이아몬드 격자가 고밀도 비정질 상태로 전이하면서 내부 압력이 급감하고, 물질 흐름이 이를 따라가지 못해 진공이 형성되는 현상이다. 넷째, 공극‑유도 액체 스팔레이션은 형성된 공극이 새로운 자유 표면을 제공하고, 반사된 충격파가 액체 혹은 기체 상태의 물질을 다시 표면으로 밀어내어 액적을 방출한다.

이러한 메커니즘은 기존 금속에서 관찰되는 **상폭발(phase explosion)**과는 근본적으로 다르며, 실리콘의 경우 비열 효과가 절삭 깊이와 크레이터 형태에 결정적인 영향을 미친다. MOD*를 적용한 시뮬레이션은 실험적 절삭 깊이와 거의 1:1 일치함을 보이며, MOD(전통 Tersoff)와 비교해 비열 효과를 무시했을 때 과도한 상폭발과 얕은 절삭을 예측한다. 또한, 저자들은 대규모 2D 시뮬레이션(수 마이크로미터 규모)과 작은 1D 시뮬레이션을 조합(composition)함으로써 계산 비용을 크게 절감하면서도 크레이터 형태와 비열 구역 직경을 정확히 재현한다는 실용적 방법론을 제시한다.

마지막으로, 전자 온도 의존적인 MOD*를 이용해 전자 온도‑압력 상도표를 열역학적 적분법으로 구축하였다. 전자 온도가 상승할수록 결합 약화가 일어나 융해·전이 온도가 낮아지는 현상이 확인되었으며, 이는 고여기 상태에서 실리콘이 β‑Sn 고압 상으로 전이하거나 비열적으로 액체‑비정질 혼합 상태에 머무를 가능성을 시사한다. 전체적으로 이 연구는 비열 효과를 정량화하고, 초고속 레이저 가공에서 실리콘의 물질 변형을 예측·제어하기 위한 이론·수치 기반을 제공한다.