[논문 리뷰] Non-equilibrium dynamics in Ising like models with biased initial condition

Reshmi Roy, Parongama Sen|arXiv (Cornell University)|2021. 06. 24.

Opinion Dynamics and Social Influence참고 문헌 31인용 수 2

한 줄 요약

이 논문은 편향된 초기 조건을 가진 이징 유사 모델에서 비평형 상태의 영온도 글라우버 동역학을 조사하며, 고정점 안정성과 탈출 확률 E(x₀)에 초점을 맞춘다. 평균장 이론과 격자 시뮬레이션을 사용하여, z=2일 경우 자화 보존으로 인해 E(x₀) = x₀임을 보이며, z > 2인 경우 비선형적이고 체계 크기 의존적인 탈출 확률이 나타나며, 스케일링 변수 y = (x₀ − x_c)/x_c L^{1/ν}를 사용해 유니버설 스케일링 함수 f(y) = [1 + tanh(λy)]/2로 수렴함을 보여, 다양한 모델과 격자에서 유니버설 임계 지수를 규명한다.

ABSTRACT

We investigate the dynamical fixed points of the zero temperature Glauber dynamics in Ising-like models. The stability analysis of the fixed points in the mean field calculation shows the existence of an exponent that depends on the coordination number $z$ in the Ising model. For the generalised voter model, a phase diagram is obtained based on this study. Numerical results for the Ising model for both the mean field case and short ranged models on lattices with different values of $z$ are also obtained. A related study is the behaviour of the exit probability $E(x_0)$, defined as the probability that a configuration ends up with all spins up starting with $x_0$ fraction of up spins. An interesting result is $E(x_0) = x_0$ in the mean field approximation when $z=2$, which is consistent with the conserved magnetisation in the system. For larger values of $z$, $E(x_0)$ shows the usual finite size dependent non linear behaviour both in the mean field model and in Ising model with nearest neighbour interaction on different two dimensional lattices. For such a behaviour, a data collapse of $E(x_0)$ is obtained using $y = \frac{(x_0 - x_c)}{x_c}L^{1/ u}$ as the scaling variable and $f(y)=\frac{1+ anh(\lambda y)}{2}$ appears as the scaling function. The universality of the exponent and the scaling factor is investigated.

연구 동기 및 목표

• 편향된 스핀 분율을 가진 이징 유사 모델에서 영온도 글라우버 동역학의 동적 고정점 안정성을 분석하기.
• 연결 수 z가 탈출 확률 E(x₀)에 미치는 영향을 조사하고, z=2를 초월해 자화 보존이 유지되는지 여부를 탐구하기.
• 스케일링 변수 y = (x₀ − x_c)/x_c L^{1/ν}와 유니버설 스케일링 함수 f(y) = [1 + tanh(λy)]/2를 사용하여 E(x₀)의 유한 체계 스케일링 형태를 유도하고 검증하기.
• 평균장 및 국한된 상호작용 모델에서 2차원 격자에 대해 임계 지수 ν와 λ의 유니버설 성질을 확립하기.
• 평균장 고정점 안정성 분 析를 통해 일반화된 보어 모델의 상도도를 그려내기.

제안 방법

• 연결 수 z에 기반한 스핀 전환 확률을 사용하여 평균장 근사에서 상향 스핀 비율 x(t)에 대한 마스터 방정식 유도하기.
• 마스터 방정식에서 시간 진동 dx/dt를 계산하고, z = 4 이상일 경우 고정점 x* = 0, 0.5, 1을 식별하기.
• x = x* + δ(t)로 고정점 주변의 선형 안정성 분석을 수행하고, 이동 항의 도함수의 부호를 평가하기.
• 정사각형 및 삼각 격자, 평균장 케이스에서 이징-글라우버 모델을 수치적으로 시뮬레이션하여 다양한 체계 크기에서 E(x₀) 계산하기.
• 유한 체계 스케일링을 적용하여 스케일링 변수 y = (x₀ − x_c)/x_c L^{1/ν}를 도입하고, 함수 형태 f(y) = [1 + tanh(λy)]/2를 사용해 데이터 수렴 검증하기.
• z = 2, 4, 6 및 다양한 격자 유형 간에 추출된 지수 ν와 λ를 비교하여 유니버설 성 평가하기.

실험 결과

연구 질문

• RQ1 탈출 확률 E(x₀)가 모든 연결 수 z에 대해 선형적(즉, E(x₀) = x₀)으로 유지되는가, 아니면 z > 2일 경우 비선형적으로 변하는가?
• RQ2 다양한 z에 대해 평균장 이징 모델에서 x* = 0.5 고정점의 성격은 어떠한가? 이 고정점의 안정성은 z에 따라 어떻게 달라지는가?
• RQ3 스케일링 변수 y = (x₀ − x_c)/x_c L^{1/ν}를 사용하여, 다양한 z와 격자 유형 간에 E(x₀)의 유한 체계 스케일링이 동일한 함수로 수렴하는가?
• RQ4 스케일링 수렴을 통해 추출된 임계 지수 ν와 λ는 평균장 및 국한된 상호작용 모델 간에 유니버설한가?
• RQ5 일반화된 보어 모델의 상도도는 평균장 프레임워크에서 고정점의 안정성 분석을 통해 어떻게 도출되는가?

주요 결과

• z = 2일 경우 평균장 이징 모델은 자화를 보존하여 E(x₀) = x₀를 보이며, 이는 1차원 이징 및 보어 모델과 일치한다.
• z > 2일 경우 x* = 0.5 고정점은 불안정해지며, 이는 체계 크기 의존성과 함께 비선형 탈출 확률 행동로의 전이를 나타낸다.
• z = 4 및 z = 6인 정사각형 및 삼각 격자에서의 E(x₀) 수치 결과는 강한 체계 크기 의존성을 보이며, 데이터가 유니버설 스케일링 함수 f(y) = [1 + tanh(λy)]/2로 수렴한다.
• 스케일링 변수 y = (x₀ − x_c)/x_c L^{1/ν}에서의 지수 ν′은 평균장 케이스에서는 z에 독립적이지만, 2차원 격자에서 구한 값보다 작다.
• 함수 f(y) = [1 + tanh(λy)]/2는 다양한 z와 격자 유형 간에 데이터를 성공적으로 수렴시켜 스케일링 형태의 유니버설 성을 나타낸다.
• 결과는 열역학적 극한에서 z > 2일 경우 E(x₀)가 단계 함수로 수렴할 것임을 시사하며, 이는 x* = 0.5 고정점의 불안정성과 일치한다.