글래스 동역학의 수치적 재규격화

초록

글래스 시스템의 긴 시간 동역학 시뮬레이션은 계산 비용이 시간의 세제곱에 비례하는 난제였다. 본 연구는 2차원 그린 함수 보간을 기반으로 한 새로운 알고리즘을 제안하여 계산 비용의 서브리니어 스케일링을 달성했다. 이 방법을 spherical mixed p-spin 모델의 퀀치 동역학에 적용하여, 초기 상태의 유한 온도에서 강한 에르고딕성 붕괴와 약한 에르고딕성 붕괴 사이의 상전이가 존재함을 규명했다. 기존보다 천 배 긴 시간尺度에 도달하여 이 상전이의 임계 지수를 측정한 결과, 이들이 연속적으로 변하는 비보편적(non-universal) 성질을 가짐을 발견했다.

상세 분석

본 논문의 핵심 기술적 기여는 글래스와 같이 장기 메모리와 비정상 상태(non-stationary) 특성을 보이는 시스템의 동역학 방정식을 효율적으로 푸는 ‘동역학적 재규격화 알고리즘’이다. 시스템의 나이(t)가 증가함에 따라 발생하는 emergent 시간尺度를 재규격화 그룹(RG) 관점에서 해석한다. 즉, 거친 특징(scale ~ t)은 RG 관련(relevant) 변수로, 먼 과세의 세부사항은 RG 무관(irrelevant) 변수로 간주하여 후자를 체계적으로 제거하는 전략을 취한다.

구체적인 알고리즘은 두 시간 변수 (t, t’)로 표현되는 그린 함수(상관 함수 C, 응답 함수 R)를 (t, θ=t’/t)로 재매개변수화한다. θ에 대해 고정된 불규칙 격자(N개 점)를 정의하고, t에 대해서만 적응형(adaptive) 시간 단계를 사용하여 진화시킨다. 이 때 θ 격자는 미시적 시간尺度(τ_micro) 부근에서 조밀하도록 설계되어, 모든 RG 관련 메모리 효과를 포착한다. 시간 t가 b배 증가할 때 θ 격자 점의 수 N을 고정하는 것은 RG의 coarse-graining 단계에 해당한다.

가장 큰 계산 난제는 메모리 적분(예: ∫ds A(t,s)B(s,t’))의 평가이다. 알고리즘은 (t, θ) 공간에서 정의된 불규칙한 2차원 격자에서 고차(3차) 에르미트 보간법을 사용하여 이 적분을 상수 시간(O(N²))에 계산한다. 보간 계수는 사전 계산 가능하며, 적응형 t-격자에서의 보간은 이전 시간 단계 결과를 초기값으로 사용하여 효율적으로 수행된다. 결과적으로 각 시간 단계의 계산 비용이 고정되고, 적응형 t-격자와 결합되어 전체 계산 비용이 서브리니어(초선형) 스케일링을 달성한다.

이 방법을 spherical mixed p-spin 모델의 영온도 퀀치 동역학에 적용했다. 모델의 동역학은 카다노프-바움 방정식 형태의 적분-미분 방정식(7)으로 기술된다. 초기 온도 T에 따른 장시간 동역학을 시뮬레이션한 결과, 세 가지 상을 발견했다: 낮은 T에서는 초기 상태와 유한 중첩을 유지하는 ‘강한 글래스’ 상(강한 에르고딕성 붕괴), 중간 T에서는 초기 상태 기억을 완전히 잃는 ‘약한 글래스’ 상, 매우 높은 T에서는 에이징 없이 정상 상태로 수렴하는 ‘강자성’ 상. 강한 글래스와 약한 글래스 사이의 상전이 임계 온도 T_c^u를 최초로 정확히 규명했으며, 이 상전이의 임계 지수(η, ν)가 모델 파라미터(p, s, λ)에 따라 연속적으로 변화하는 비보편적 성질을 보인다는 점이 특히 중요하다. 이는 해당 상전이가 평형 상전이와는 근본적으로 다른 비보편적 동역학적 상전이임을 시사한다.