초고속 레이저 물질 상호작용을 위한 다중스케일 맥스웰 TTM MD 모델

초록

**
본 논문은 유한 차분 시간 영역(FDTD) 방식으로 전자기 파동을 전파하고, 두 온도 모델(TTM)과 분자 동역학(MD)을 결합한 다중스케일 시뮬레이션 프레임워크(M‑TTM‑MD)를 제시한다. 전자기장에 의한 에너지 흡수를 공간·시간적으로 해석하고, 전자‑격자 간 에너지 교환을 TTM으로, 원자 수준의 구조 변화를 MD으로 동시에 계산한다. 물질의 광학 상수는 보조 미분 방정식(ADE)으로 동적으로 업데이트되며, 표면 형상 변화가 FDTD 격자에 반영되어 반복적인 피드백 루프를 형성한다. MPI 기반 병렬 구현을 통해 대규모 시뮬레이션이 가능하며, 금속 표면에 대한 레이저 유도 주기 구조(LIPSS)·증발·스퍼터링 등 비평형 현상을 정량적으로 예측한다.

**

상세 분석

**
이 연구는 초초단 펨토초 레이저가 금속 표면에 입사할 때 발생하는 복합 물리 현상을 하나의 수치 체계에 통합한 점에서 혁신적이다. 기존의 두 온도 모델(TTM)은 전자와 격자 사이의 열 교환을 매크로 스케일에서 기술했지만, 전자기장의 공간적 비균일성이나 표면 플라스몬, 구조적 변형을 반영하지 못했다. 반면 전통적인 MD는 원자 간 상호작용을 정확히 다루지만, 외부 광원으로부터의 에너지 주입을 단순히 균일 가열 형태로 가정한다는 한계가 있었다. 저자들은 세 가지 수치 기법을 서로 다른 격자(FDTD, TTM, MD) 위에 배치하고, 격자 간 매핑 규칙을 정밀하게 정의함으로써 에너지 보존과 시간 동기화를 확보하였다.

특히 전자기 해석에서 ADE를 이용해 금속의 복소 유전율을 전자 온도와 레이저 주파수에 의존하도록 동적으로 업데이트한다는 점은, 레이저 펄스 진행 중에 발생하는 광학 상수 변화를 실시간으로 반영할 수 있게 한다. 이는 표면 형상이 변하면서 발생하는 산란 및 플라스몬 공명 변화를 정확히 포착하고, 결과적으로 Qₗₐₛ = –∇·S 형태의 에너지 흡수 항을 공간적으로 세분화된 형태로 TTM에 전달한다.

시간적 스케일링 측면에서 FDTD는 10⁻³ fs 수준의 타임스텝을, TTM은 10⁻² fs, MD는 수십 femtosecond 정도의 타임스텝을 사용한다. 저자들은 ‘카운터 변수’를 도입해 서로 다른 타임스텝을 동기화하고, FDTD 서브루프를 레이저 펄스가 전파될 때까지 진행한 뒤, 평균화된 전자 에너지 흐름을 TTM 격자에 매핑한다. 이 과정에서 전자 온도 구배에 따른 열 전도(Kₑ)와 전자‑격자 결합 상수(G)를 이용해 전자 온도 방정식을 풀고, 얻어진 전자 온도 분포를 MD에 전달한다. MD 단계에서는 임베디드 원자 모델(EAM) 포텐셜을 사용해 원자 간 힘을 계산하고, 온도와 압력은 각각 equipartition theorem과 virial theorem을 통해 추출한다.

병렬 구현은 MPI 기반 도메인 디컴포지션을 활용한다. 각 프로세스는 자체 서브도메인을 담당하고, ‘스킨 레이어’를 통해 경계 셀의 전자기장·원자 정보를 교환한다. 이는 통신 오버헤드를 최소화하면서도 대규모(수백만 원자) 시뮬레이션을 가능하게 한다. 또한, 물질 밀도가 일정 이하가 되면 해당 셀을 ‘비활성’으로 전환해 CPML 경계 조건과 결합함으로써 수치적 안정성을 유지한다.

이러한 설계는 레이저 펄스가 물질을 가열·용융·증발시키는 전 과정을 원자 수준에서 직접 관찰할 수 있게 하며, 특히 LIPSS 형성 메커니즘을 전자기 파동의 산란·플라스몬 공명, 전자‑격자 열 교환, 원자 재배열이라는 세 축으로 해석할 수 있게 만든다. 결과적으로 에너지 전달 효율, 표면 주기, 피크 온도 등 실험적 관측값과 정량적으로 비교·검증할 수 있는 강력한 도구를 제공한다.

**