실리콘 나노입자 레이저 전진 전사(LIFT) 중 결정화 메커니즘의 원자 규모 시뮬레이션

초록

본 연구는 분자동역학(MD) 시뮬레이션을 이용해 레이저‑유도 전진 전사(LIFT) 과정에서 실리콘 나노입자의 용융‑응고와 결정화 현상을 조사한다. Stillinger‑Weber 포텐셜을 적용해 2‑20 nm 직경 입자를 공기 중에서 냉각했으며, 입자 크기가 작아질수록 결정화 확률이 급격히 감소하고, 핵생성은 주로 표면 아래 5 Å 이내에서 일어난다는 것을 발견했다. 냉각 속도와 열전도도에 따라 단결정, 다결정, 비정질 구조가 전이되며, 느린 냉각(낮은 열전도도)에서는 100 % 결정화가 가능함을 제시한다. 이러한 결과는 LIFT 기반 단결정 실리콘 제조에 중요한 설계 지침을 제공한다.

상세 분석

본 논문은 레이저‑유도 전진 전사(LIFT) 프린팅 과정에서 실리콘(Si) 나노입자가 용융된 뒤 공기 중을 비행하면서 어떻게 고체화되는지를 원자 수준에서 규명하고자 한다. 이를 위해 저자들은 LAMMPS 기반의 분자동역학(MD) 시뮬레이션을 수행했으며, 실리콘의 결합 특성을 잘 재현하는 Stillinger‑Weber(SW) 포텐셜을 선택하였다. 시뮬레이션은 2 nm에서 20 nm까지 직경이 다른 구형 입자를 생성하고, 300 K에서 2000 K까지 10 ns 가열 후 5 ns 동안 2000 K에서 유지해 완전한 비정질 상태를 확보한다. 이후 1700 K에서 시작해 NVE ensemble으로 냉각시키며, 열전도도(6.38 × 10⁵ W m⁻² K⁻¹ 등)를 조절해 다양한 냉각 속도를 구현한다.

핵심 결과는 세 가지 축에서 도출된다. 첫째, 입자 크기가 감소할수록 용융‑응고 과정에서 관찰되는 잠열 방출이 사라지고, 결정화 온도와 융점이 모두 크게 낮아진다. 2 nm 입자는 어떠한 잠열 급증도 보이지 않아 사실상 비정질 상태를 유지한다. 이는 표면/부피 비율이 증가하면서 라플라스 압력이 상승하고, 원자 배열의 자유도가 제한돼 핵생성 장벽이 높아지는 현상과 일치한다.

둘째, 냉각 속도(열전도도)와 결정화 형태 사이의 상관관계가 명확히 드러난다. 높은 열전도도(빠른 냉각, ≈10⁸ K ns⁻¹)에서는 2 nm·4 nm 입자 모두 잠열 방출 없이 비정질 상태가 유지되지만, 8 nm·12 nm 입자는 잠열 급증과 함께 다결정(poly‑Si) 혹은 단결정(single‑Si) 구조가 형성된다. 열전도도를 1.99 × 10⁵ W m⁻² K⁻¹ 이하로 낮추면 냉각 속도가 약 13 K ns⁻¹로 감소하고, 8 nm 입자 10 회 실험 모두 결정화가 일어나 100 % 성공률을 보인다. 이는 고전적 핵생성 이론에서 제시하는 ‘느린 냉각 → 낮은 과냉각 → 높은 핵생성률’과 반대로, 충분히 낮은 과냉각이 핵의 성장 속도를 가속화해 큰 단결정이 형성된다는 점을 시사한다.

셋째, 핵생성 위치는 거의 항상 표면 아래 5 Å 이내의 서브‑서피스 영역에서 시작한다. BOP(결합 차수 파라미터) 분석을 통해 q₃ 값을 기준으로 결정 원자를 식별했으며, 가장 큰 결정 클러스터가 초기에는 얇은 구형 껍질 형태로 존재하다가 온도 급증 시점에 중심부로 이동하면서 급격히 성장한다. 단결정 형성 경우(예: case 10) 성장 속도는 1575 atoms ns⁻¹에 달했으며, 다결정 경우는 142 atoms ns⁻¹, 비정질 경우는 15 atoms ns⁻¹에 불과했다. 이는 핵이 서브‑서피스에서 형성된 뒤 전체 입자 내부로 퍼지는 과정이 결정 구조의 품질을 좌우함을 보여준다.

이러한 결과는 LIFT 프린팅에서 입자 크기와 냉각 환경을 정밀하게 제어하면, 원하는 단결정 Si 나노입자를 직접 제조할 수 있음을 의미한다. 특히, 열전도도(또는 주변 가스 압력·열전달 조건)를 조절해 냉각 속도를 10 K ns⁻¹ 이하로 낮추면 8 nm 이상 입자에서 거의 확실히 단결정이 형성된다. 반면, 작은 입자(≤4 nm)에서는 물리적 한계(라플라스 압력, 표면 에너지) 때문에 비정질 상태가 지속될 가능성이 크다. 따라서 LIFT 기반 고성능 Si 디바이스 제작 시, 입자 직경을 8 nm 이상으로 설계하고, 프린팅 후 급격한 열손실을 방지하는 환경(예: 저압 가스, 열전도도 낮은 매트릭스) 조성이 필요하다.

마지막으로, 시뮬레이션 설정의 제한점도 언급한다. 진공 환경 가정, 고정된 SW 포텐셜, 그리고 입자와 주변 가스 사이의 상호작용을 무시한 점은 실제 LIFT 공정에서 발생할 수 있는 열전달·충돌 효과를 완전히 반영하지 못한다. 향후 연구에서는 다중 입자 상호작용, 기판 충돌, 그리고 레이저 펄스에 의한 급격한 에너지 주입을 포함한 전동역학 모델링이 필요하다.