아폴로니안 네트워크 위의 q 상태 포츠 모델

초록

본 논문은 아폴로니안 네트워크에 q‑상태 포츠 모델을 적용하여 몬테카를로 시뮬레이션과 전이 행렬(Transfer Matrix) 방법으로 자발적 자기화, 상관 길이, 엔트로피, 비열을 온도와 q의 함수로 조사한다. 다양한 q값에 대해 온도 스케일링 함수를 제시하고, 두 방법 간 정량적 일치를 확인하였다. 열역학적 극한에서도 임계 현상이 나타나지 않음을 보고한다.

상세 분석

아폴로니안 네트워크는 초기 삼각형에서 시작해 매 단계마다 기존의 삼각형 내부에 새로운 정점을 삽입하고, 그 정점을 기존의 세 정점과 연결하는 재귀적 과정을 통해 생성되는 스케일‑프리 및 작은 세계 특성을 동시에 갖는 복합 네트워크이다. 이러한 구조는 차수가 높은 정점이 소수의 중심에 집중되는 강한 이질성을 제공하므로, 전통적인 격자계에서 관찰되는 임계 현상이 어떻게 변형되는지를 탐구하기에 이상적인 시험대가 된다.

논문은 먼저 q‑상태 포츠 모델을 정의한다. 각 정점 i는 σ_i∈{1,…,q}의 스핀 변수를 가지며, 인접 정점 i와 j 사이의 상호작용 에너지는 −Jδ(σ_i,σ_j) 로 표현된다(δ는 크로네커 델타, J>0는 강자성 결합). 외부 자기장이 없으며, 자유 에너지와 열역학적 양은 전통적인 통계역학적 정의에 따라 계산된다.

연구진은 두 가지 독립적인 수치적 접근법을 채택하였다. 첫 번째는 표준 메트로폴리스 알고리즘을 이용한 몬테카를로 시뮬레이션으로, 네트워크 규모 N을 10^4에서 10^6까지 확장하면서 온도 T를 0.1부터 5.0까지 0.05 간격으로 샘플링하였다. 각 온도에서 10^5개의 마코프 체인 스텝을 열화 단계로 사용하고, 이후 5×10^5 스텝을 측정 단계로 활용해 평균 자기화 M, 상관 길이 ξ, 엔트로피 S, 비열 C를 추정하였다. 두 번째는 전이 행렬 기법이다. 아폴로니안 네트워크는 재귀적 생성 규칙을 갖기 때문에, 한 단계에서 다음 단계로의 전이 행렬 T_n을 정확히 구성할 수 있다. T_n은 각 가능한 스핀 배치의 가중치를 포함하며, 행렬의 최대 고유값 λ_max를 통해 자유 에너지 f=−k_BT ln λ_max/N을 얻는다. 이로부터 M, ξ, S, C를 미분적으로 계산한다.

두 방법 모두 q가 증가함에 따라 전이 온도(T_c와 유사한 특성)가 점차 상승하지만, ξ와 C의 피크가 네트워크 규모 N→∞에서 점점 얇아지며, 결국 유한한 값에 머무른다. 이는 전통적인 2차 상전이에서 관찰되는 발산형 상관 길이와는 근본적으로 다르다. 특히, 전이 행렬 분석을 통해 ξ∝N^0 (즉, 규모에 무관)임을 확인했으며, 이는 무한 네트워크에서도 장거리 순서가 형성되지 않음을 의미한다.

또한, 엔트로피 S는 온도가 낮을 때 거의 0에 수렴하고, 고온에서는 ln q에 가까운 값을 취한다. 이는 q개의 가능한 스핀 상태가 모두 동등하게 접근 가능한 무작위 상태로 전이함을 반영한다. 비열 C는 온도에 대해 넓은 단일 피크를 보이지만, 피크 높이는 N에 비례하지 않으며, q가 커질수록 피크가 약간 이동하고 폭이 넓어진다. 이러한 결과는 네트워크의 높은 차수 정점이 국부적인 에너지 장벽을 완화시켜, 전통적인 격자계에서 나타나는 급격한 에너지 변화를 완화시키는 역할을 함을 시사한다.

스케일링 분석에서는 온도와 q를 결합한 복합 스케일링 변수 x = (T−T_0)/q^α를 도입해 M, C, ξ 데이터를 모두 하나의 곡선으로 collapse 시켰다. 여기서 T_0는 q→∞에서의 한계 온도이며, α≈0.5 정도의 지수가 최적화되었다. 이는 q가 클수록 시스템이 더 높은 온도에서도 유사한 거동을 보이며, q에 대한 선형적인 스케일링이 아니라 제곱근 형태의 비선형 스케일링이 적용된다는 점을 강조한다.

결론적으로, 아폴로니안 네트워크 위의 q‑상태 포츠 모델은 전이 행렬과 몬테카를로 시뮬레이션 모두에서 일관된 열역학적 거동을 보이며, 무한 규모에서도 임계점이 존재하지 않는다. 이는 네트워크의 복합적인 차수 분포와 작은 세계 특성이 전이 현상을 억제하고, 전통적인 격자계와는 다른 보편적 거동을 만든다는 중요한 물리적 통찰을 제공한다.