보존적 플러킹 모델에서 재진입 현상

초록

이 논문은 스핀‑속도 결합을 포함한 보존적 해밀토니안 플러킹 모델(HFM)에서, 자기 구동 강도 K를 증가시킬 때 균일 가스 → 클러스터 → 다시 균일 가스로 전이하는 재진입 현상을 수치적으로 확인한다. 재진입은 강한 스핀‑속도 결합이 횡방향 확산을 억제해 입자 이동성을 제한하고, 결과적으로 상분리 창을 닫는 ‘동역학적 좌절’ 메커니즘에 기인한다.

상세 분석

본 연구는 기존의 비보존적 활성 물질 모델과 달리, 전통적인 라그랑지안/해밀토니안 체계 안에서 스핀‑속도 결합(K)이라는 내부 구동원을 도입한 점이 핵심이다. 해밀토니안(1)은 입자 운동량 p_i에 스핀 S_i가 선형으로 결합되는 항 K S_i를 포함하고, 이는 전자기학에서 벡터 퍼텐셜에 대한 최소 결합과 수학적으로 동등하다. 이 구조는 스핀을 내부 게이지 장으로 해석하게 하며, 스핀‑속도 결합이 강해질수록 입자들의 실질적 질량이 증가하고, 이동성 매트릭스 M_i(θ_i)에서 비대각 성분 K n⊥_i가 크게 작용한다. 결과적으로 유효 마찰 Δ = γ_t γ_r + K²가 커져, 이동성 M_rr와 M_θθ이 감소하고, 특히 횡방향 확산을 억제한다.

동역학은 과잉 감쇠(over‑damped) Langevin 방정식(5)–(9)으로 기술되며, 잡음은 플럭투에이션‑디소시페이션 관계를 만족한다. 수치 시뮬레이션은 2차원 주기적 박스와 슬래브(길이 비율 L_x = 3 L_y) 설정에서 수행되었으며, 입자 수 N = 250–5000, 밀도 η ≈ 0.3, 온도 T = 0.6, 스핀‑스핀 결합 J = 1을 기본값으로 두었다.

주요 관찰은 다음과 같다. (1) K가 0에서 시작해 증가하면, 초기 균일 가스 상태에서 스핀 정렬에 의해 유도된 유효적인 인력이 발생, 이는 밀도 차이 Δρ = ρ_high − ρ_low을 크게 만든다. (2) K가 특정 값 K_peak ≈ √(3 γ_t γ_r) 근처에서 Δρ가 최대에 도달하고, 이후 K가 더 커지면 Δρ가 급격히 감소해 다시 균일 가스로 돌아간다. 이는 ‘재진입’ 현상이다. (3) K와 J의 비율 K/√J가 재진입 곡선을 결정하며, J/T가 임계값 J* ≈ 1.2 이하이면 스핀 정렬이 열 잡음에 의해 파괴돼 상분리가 일어나지 않는다. (4) 높은 K에서는 입자 속도 분포 P(v)가 Maxwell‑Boltzmann 형태에서 저속 중심의 비대칭 형태로 변하며, 이는 ‘동역학적 좌절’에 해당한다. (5) 국부 자기화 m_local와 부호 자기화 m_x의 확률 분포를 조사하면, 중간 K 구간에서는 이중 피크(또는 삼중 피크) 구조가 나타나 클러스터 내부에 반대 극성 도메인이 공존함을 보여준다. 반면 큰 K에서는 밀도는 균일하지만 m_x는 양극성 이중 피크를 유지해, 긴 시간 동안 도메인 벽이 거의 정지한 상태가 유지된다. 이는 슬래브 기하학에서 공간 대칭이 깨질 때만 관찰되는 현상이다.

이러한 결과는 기존 MIPS(모티리티‑유도 상분리)에서 ‘활성도(Pe)’가 너무 커지면 재진입이 일어난다는 보고와 정량적으로 일치한다. 다만, 여기서는 비보존적 에너지 주입이 아니라 보존적 스핀‑속도 결합이라는 내부 자유도가 동일한 비단조적 거동을 만들어낸다. 따라서 재진입 현상은 ‘활성’ 자체가 아니라, 구동과 이동성 사이의 경쟁, 즉 ‘동역학적 좌절’에 의해 보편적으로 발생할 수 있음을 시사한다.