화학적 패턴 기판 위 이중유체의 표면유도 스핀노달 분해: 분자동역학 시뮬레이션

초록

본 연구는 체커보드 형태의 화학적 패턴을 가진 비정질 기판 위에서 이진 유체 혼합물의 표면유도 스핀노달 분해(SDSD)를 분자동역학(MD)으로 조사한다. 패턴 주기 λ와 혼합물의 스핀노달 길이 λ_c의 관계에 따라 패턴이 유체에 전이되는 ‘레지스트리(registry)’ 현상이 나타나며, L⊥(t)∝t¹⁄³, L∥(z,t)∝t^α와 같은 확산적·동적 스케일링 법칙이 확인된다. 또한 패치 크기에 따른 레지스트리 형성·소멸 시간의 스케일링도 제시한다.

상세 분석

본 논문은 기존의 Cahn‑Hilliard‑Cook(CHC) 기반 연구와 달리 전통적인 유체역학을 자연스럽게 포함하는 전통적인 분자동역학(MD) 접근법을 채택했다는 점에서 의미가 크다. Lennard‑Jones(LJ) 포텐셜을 이용해 대칭적인 A‑B 혼합물을 모델링하고, 기판은 비정질 구조로 설계해 평면성에 의한 층화 효과를 최소화하였다. 화학적 패턴은 체커보드 형태의 정사각형 패치(M×M)로 구성했으며, 각 패치는 A 또는 B 입자를 선호하도록 ε_AS와 ε_BS 값을 조절하였다.

핵심 변수는 패턴 주기 λ(≈2Mσ)와 혼합물의 스핀노달 길이 λ_c(≈2π/ξ_B)이다. λ>λ_c인 경우, 즉 패턴이 스핀노달 길이보다 크게 될 때, 기판의 화학적 정보를 유체가 ‘레지스트리’ 형태로 전이한다는 현상이 뚜렷이 관찰된다. 이는 표면‑법선 방향 상관함수 C⊥(z,t)와 표면 필드와의 상관분석을 통해 정량화되었으며, 연관된 퍼짐 길이 L⊥(t)는 λ이 작아질수록(즉, 패턴이 촘촘해질수록) 증가한다.

시간에 따른 L⊥(t)의 성장률은 t^{1/3}에 근접하는 확산적 스케일링을 보였으며, 이는 균일한 기판에서 보고된 습윤층 성장(R₁(t)∝t^{1/3})과 일치한다. 이는 표면‑유도 스핀노달 분해가 패턴 유무에 관계없이 기본적인 확산 메커니즘에 의해 지배된다는 점을 시사한다.

또한, 표면에 평행한 x‑y 평면에서의 조성 파동이 형성되었으며, 그 파동의 파장은 L∥(z,t)∝t^{α}로 시간에 따라 확대된다. α는 위치(z)에 따라 변동하는데, 초기에는 ‘표면‑레지스트리’ 영역에서 α≈0.20.3 정도의 느린 성장률을 보이다가, 레지스트리가 소멸하고 순수한 상분리 단계에 진입하면 α≈1/31/2 수준으로 전이한다. 이는 레지스트리 영역이 존재할 때는 기판의 화학적 강제조건이 파동 전파를 억제하지만, 시간이 흐르면 bulk의 상분리 역학이 우세해지는 과정을 반영한다.

패치 크기(M)에 대한 스케일링 분석에서는 레지스트리 형성 시간 τ_form∝M^{z_f}, 소멸 시간 τ_melt∝M^{z_m}와 같은 거듭제곱 관계가 확인되었다. 구체적인 지수값은 시뮬레이션 파라미터에 따라 다소 차이가 있지만, 일반적으로 τ_form>τ_melt이며, 큰 패치일수록 레지스트리 유지가 오래 지속된다. 이는 실험적 패턴 설계 시, 원하는 도메인 구조를 얻기 위해 패치 크기를 조절함으로써 레지스트리의 수명과 전이 정도를 제어할 수 있음을 의미한다.

마지막으로, MD 시뮬레이션이 제공하는 전통적인 유체 흐름(수소역학) 효과는 기존 CHC 기반 연구에서 관찰되지 않았던 빠른 도메인 성장과 복잡한 파동 전파를 가능하게 한다. 이는 특히 나노·마이크로 스케일 디바이스에서 패턴 전이 효율을 높이는 데 중요한 역할을 할 것으로 기대된다.