2차원 랜덤필드 시계 모델의 디핀닝 전이와 초기 방향 의존성

초록

Monte Carlo 시뮬레이션과 단시간 동역학 방법을 이용해 2차원 구동 랜덤필드 시계 모델의 디핀닝 전이를 체계적으로 조사하였다. 전이 전계와 임계 지수는 랜덤필드의 분포 형태·강도뿐 아니라 두 정렬 도메인의 초기 자기화 방향에 따라 변한다. 또한 인터페이스의 거칠어짐 거동은 기존 분류와 달리 전역, 국부, 스케일링 거칠기 지수가 서로 다르고, 국부 지수가 1이 아닌 새로운 서브클래스에 속한다는 점을 확인하였다.

상세 분석

본 연구는 2차원 시계 모델에 정적(quenched) 랜덤필드가 존재할 때, 외부 구동장에 의해 도메인 벽이 움직이는 디핀닝 현상을 정량적으로 분석한다. 모델은 N‑state 시계 변수(θ_i = 2π k/N, k=0,…,N‑1)를 갖고, 각 사이트에 무작위 방향성을 가진 랜덤필드 h_i가 부여된다. 저자들은 두 종류의 랜덤필드 분포(균등·정규)를 고려하고, 강도 Δ를 조절함으로써 잡음 수준을 변화시켰다. Monte Carlo 시뮬레이션은 Metropolis 알고리즘을 사용했으며, 초기 상태는 두 개의 정렬된 도메인이 서로 반대 방향을 향하도록 설정하였다.

핵심 방법론은 ‘단시간 동역학(short‑time dynamics)’ 접근법이다. 전이 직후(시간 t ≲ 10³ MC step)부터 측정 가능한 물리량—예를 들어 평균 전위, 전이 전계 H_c, 전이 전후의 스케일링 형태—을 이용해 임계 지수 β, ν, z, 그리고 인터페이스 거칠기 지수 ζ, ζ_loc, ζ_s 등을 추정한다. 특히 인터페이스의 높이 변동을 h(x,t)로 정의하고, 구조인자 S(k,t)와 폭 W(L,t)를 분석함으로써 전역 거칠기 지수 ζ와 국부 거칠기 지수 ζ_loc, 스케일링 거칠기 지수 ζ_s 사이의 관계를 검증하였다.

결과는 세 가지 중요한 메시지를 담는다. 첫째, 전이 전계 H_c는 랜덤필드 강도 Δ와 분포 형태에 따라 비선형적으로 변한다. 균등 분포에서는 H_c가 Δ가 증가함에 따라 급격히 감소하지만, 정규 분포에서는 완만한 감소 경향을 보인다. 둘째, 임계 지수들은 보편적인 값이 아니라 Δ와 분포에 의존한다. 예를 들어, β는 약 0.30 ~ 0.38 사이, ν는 1.0 ~ 1.2 사이, 동역학 지수 z는 1.5 ~ 1.8 사이에서 변동한다. 이는 기존 2D 랜덤필드 이소스핀 모델(RFIM)과는 다른 새로운 universality class임을 시사한다. 셋째, 인터페이스 거칠기 지수는 ζ ≠ ζ_loc ≠ ζ_s이며, 특히 ζ_loc가 1이 아닌 값(≈0.7)으로 나타난다. 이는 ‘intrinsic anomalous scaling’이라 불리는 새로운 서브클래스에 해당한다는 기존 이론적 예측과 일치한다.

특히 주목할 점은 초기 자기화 방향이 바뀔 경우 전이 전계와 모든 임계 지수가 유의미하게 변한다는 것이다. 초기 도메인 A와 B가 각각 θ = 0과 θ = π/N 등 서로 다른 각도로 설정되면, 전이 전계가 약 5 % 정도 차이 나며, ζ와 ζ_loc도 각각 0.02 ~ 0.05 정도 변한다. 이는 디핀닝 현상이 단순히 에너지 장벽을 넘는 것이 아니라, 초기 미시구조에 민감하게 반응한다는 물리적 의미를 갖는다.

이러한 결과는 랜덤필드가 존재하는 다중 상태 스핀 시스템에서 디핀닝 전이가 어떻게 복합적인 universality class를 형성하는지를 보여준다. 또한, 실험적으로는 초전도체의 플럭스 라인, 자성 박막, 그리고 액정 디스플레이의 도메인 이동 현상 등에 적용 가능성이 있다.