[논문 리뷰] Long-range interacting classical systems: universality in mixing weakening

Alessandro Campa, Andrea Giansanti외 2명|arXiv (Cornell University)|2000. 07. 06.

Statistical Mechanics and Entropy인용 수 30

한 줄 요약

이 논문은 $d=2$ 및 $d=3$ 차원의 관성적 $XY$ 모형에 $r^{-\alpha}$ 상호작용을 가진 장거리 상호작용을 갖는 고전적 시스템에서의 동적 혼돈을 분자역학 시뮬레이션을 통해 조사한다. 결과적으로, 스케일링된 최대 리아프노프 지수는 $N^{-\kappa}$로 감소하며, 이때 $\kappa$는 $\alpha/d$의 보편 함수로 나타나며, $\alpha/d$가 0에서 1로 증가함에 따라 $1/3$에서 0으로 감소함을 보여주며, 이는 비확장 통계역학과 일치하며 비확장 영역에서 혼합의 보편적인 약화를 나타낸다.

ABSTRACT

Through molecular dynamics, we study the $d=2,3$ classical model of $N$ coupled rotators (inertial XY model) assuming a coupling constant which decays with distance as $r_{ij}^{-α}$ ($α\ge 0$). The total energy $$ is asymptotically $\propto N { ilde N}$ with ${ ilde N} \equiv [N^{1-α/d}-(α/d)]/[1-α/d]$, hence the model is thermodynamically extensive if $α/d>1$ and nonextensive otherwise. We numerically show that, for energies above some threshold, the (appropriately scaled) maximum Lyapunov exponent is $\propto N^{-κ}$ where $κ$ is an {\it universal} (one and the same for $d=1,2$ and 3, and all energies) function of $α/d$, which monotonically decreases from 1/3 to zero when $α/d$ increases from 0 to 1, and identically vanishes above 1. These features are consistent with the nonextensive statistical mechanics scenario, where thermodynamic extensivity is associated with {\it exponential} mixing in phase space, whereas {\it weaker} (possibly {\it power-law} in the present case) mixing emerges at the $N o \infty$ limit whenever nonextensivity is observed.

연구 동기 및 목표

• 장거리 상호작용을 갖는 고전적 다체계에서 시스템 크기의 영향을 동적 혼돈에 대해 조사하는 것.
• 최대 리아프노프 지수의 스케일링 지수 $\kappa$가 차원과 에너지에 관계없이 보편적인가를 테스트하는 것.
• Hamiltonian 시스템에서 $r^{-\alpha}$ 상호작용을 갖는 시스템에서의 비확장성과 약화된 혼합 간의 관계를 조사하는 것.
• 모든 차원과 에너지 밀도에 대해 $\kappa$가 $\alpha/d$에만 의존하며, 임계 에너지 이하에서 성립하는 가설을 검증하는 것.

제안 방법

• $d=2$ 및 $d=3$ 격자에서 $N$개의 결합된 회전자에 대해 분자역학 시뮬레이션을 수행하며, $r_{ij}^{-\alpha}$ 상호작용을 고려한다.
• 비확장 영역에서 열역학적 확장성을 확보하기 위해 해밀토니안을 $\tilde{N} = [N^{1-\alpha/d} - (\alpha/d)] / [1 - \alpha/d]$ 로 스케일링한다.
• 최대 리아프노프 지수 $\lambda_{\text{max}}'$ 를 계산하고, $\lambda_{\text{max}}' \propto N^{-\kappa}$ 를 통해 $N$-의존성에 대해 스케일링한다.
• 에너지 보존 오차가 $10^{-4}$ 이하인 속도-버렛 적분을 사용하며, 총 운동량이 0이고 에너지가 고정된 초기 조건을 설정한다.
• $\lambda_{\text{max}}'$ 와 $N$ 의 로그-로그 플롯을 통해 $\kappa$ 를 추출하고, 비율 $\alpha/d$ 를 이용해 차원 간에 데이터를 캐릭터화한다.
• 보편적인 $\kappa(\alpha/d)$ 함수를 기술하기 위해 히우리스틱 피팅 $\kappa = [1 - (\alpha/d)^2] / [3 + (\alpha/d)^2/2]$ 을 사용한다.

실험 결과

연구 질문

• RQ1 스케일링 지수 $\kappa$ 가 스케일링된 최대 리아프노프 지수에 대해 $\alpha/d$ 에만 의존하며, 공간 차원 $d$ 와 독립적인가?
• RQ2 임계 임계 에너지 이상에서, 최대 리아프노프 지수의 $N$-의존성은 다양한 에너지 밀도에서 보편적인가?
• RQ3 장거리 상호작용 시스템에서 혼합의 약화는 통계역학에서의 비확장성과 어떻게 관련이 있는가?
• RQ4 관측된 $\kappa(\alpha/d)$ 함수는 $1/r^\alpha$ 상호작용을 갖는 다른 고전적 시스템으로 일반화될 수 있는가?

주요 결과

• 임계 에너지 이상에서는 스케일링된 최대 리아프노프 지수가 $\lambda_{\text{max}}' \propto N^{-\kappa}$ 로 스케일링되며, 이때 $\kappa$ 는 $\alpha/d$ 에만 의존한다.
• $\alpha/d < 1$ 인 경우, $\alpha/d = 0$ 에서 $\kappa = 1/3$ 으로 시작하여 $\alpha/d = 1$ 에서 0으로 감소하며, 이는 시스템 크대에 따라 혼돈이 약화됨을 나타낸다.
• $\kappa(\alpha/d)$ 함수는 $d=1,2,3$ 전역에서 보편적이며, $\alpha/d$ 를 축으로 하여 모든 데이터가 하나의 곡선에 수렴한다.
• $\alpha/d > 1$ 인 경우, $\kappa(\alpha/d)$ 함수는 정확히 0이 되며, 이는 지수적 혼합과 광범위한 열역학과 일치한다.
• 다양한 에너지 밀도에서도 스케일링이 성립함을 확인하였으며, $\alpha=0.8$, $d=2$ 인 경우 $U=5$ 와 $U=8$ 에서 동일한 $\kappa$ 값을 관측하였다.
• 관측된 $\kappa(\alpha/d)$ 함수는 비확장 통계역학의 예측과 일치하며, 약화된 혼합을 비확장성과 연결한다.