고차 푸시안 방정식으로 본 격자 이징 모델의 자기 감수도 chi6

1. 서론 및 배경 이징 모델의 자기 감수도 χ는 고온 χ⁺와 저온 χ⁻ 두 부분으로 나뉘며, 각각은 자기 변수 w=½·s/(1+s²) 로 표현된다. χ는 n‑입자 기여 χ⁽ⁿ⁾(w) 로 전개되며, χ⁽ⁿ⁾는 (n‑1) 차원 적분식으로 정의된다. n이 커질수록 적분 차원이 증가해 급수 전개가 급격히 복잡해진다. 기존 연구에서는 χ⁽⁴⁾까지는 정확한 정수 연산으로, χ⁽⁵⁾는 소수 모듈러 연산으로 ODE를 찾아냈다. 본 연구는 이 방법을 n=6에 확대한다. 2. 급수 생성 및 ODE 공식 χ⁽⁶⁾(w)의 급수를 소수 p=32749 모듈러 하에서 w 변수에 대해 6,500 차까지 계산하였다. 급수 길이와 ODE 차수(Q)·다항식 차수(D) 사이의 관계를 N=(Q+1)(D+1)−f 로 나타내는 “ODE formula”를 이용해, 최소 항수 N₀를 구하고 최적 ODE를 설계했다. χ⁽⁶⁾에 대한 최적 ODE는 차수 Q≈84, D≈73 로, 약 5,300~5,600 항만 필요했다. 3. 감쇠된 급수 Φ⁽⁶⁾ 정의 및 ODE 도출 χ⁽⁶⁾에서 저차 입자 기여를 적절히 빼낸 감쇠된 급수 Φ⁽⁶⁾=χ⁽⁶⁾−(2/3)χ⁽⁴⁾+(2/45)χ⁽²⁾ 를 정의한다. 이는 χ⁽⁶⁾보다 낮은 차수의 ODE를 갖게 하며, 실제로 차수 46의 선형 푸시안 ODE를 얻었다. 4. 연산자 인수분해 Φ⁽⁶⁾를 만족하는 연산자 L₄₆을 소수 모듈러 환경에서 인수분해하였다. 좌측 6차 인자는 완전 타원 적분 E를 나타내는 2차 연산자 L_E의 대칭 5제곱과 동형이며, 이는 “대칭 전력” 구조가 고차 입자 기여에 보편적임을 시사한다. 오른쪽 인자는 직접합 형태로, (Z₂·N₁)⊕V₂⊕(F₃·F₂·L_{s1}) 등 여러 하위 연산자로 구성된다. 여기서 Z₂는 무게 1 모듈라 형식, V₂는 χ⁽²⁾와 동등한 2차 연산자, F₂·F₃는 전역적으로 영( nilpotent)인 연산자이다. 5. 다중 소수에서의 정확한 재구성 여러 소수(32749, 65537, 10007 등)에서 얻은 인자들을 이용해, 각 하위 연산자를 정수 계수 다항식 형태로 재구성하였다. 특히 V₂와 Z₂는 정확히 알려진 타원 적분 연산자와 동일함을 확인했으며, 직접합 구조의 각 성분도 정밀히 복원되었다. 6. 특이점 및 임계 거동 분석 ODE의 특이점 구조를 분석해, 물리적 임계점 w_c=1/2(√2−1) 근처에서 로그‑다항식 형태의 발산이 나타남을 확인하였다. 또한, 복소 평면상의 “자연 경계”(natural boundary)라 불리는 곡선이 실제 ODE의 정규 특이점이 아니라, 모듈러 연산에서 나타나는 가짜 특이점임을 증명하였다. 이는 χ⁽⁵⁾에서 제시된 가설을 n=6에서도 유지한다는 중요한 물리적 의미를 가진다. 7. 계산 복잡도와 효율성 χ⁽⁶⁾ 급수 생성은 O(N⁴ log N) 복잡도를 가지며, 6,500 차까지 계산하는 데 약 65,000 CPU‑hour가 소요되었다. 그러나 ODE formula를 활용해 최적 ODE를 찾는 데는 5,100~5,400 항만 필요했으며, 이는 최소 차수 ODE(≈49,000 항) 대비 90% 이상의 연산 절감 효과를 보여준다. 8. 결론 본 연구는 χ⁽⁶⁾에 대한 푸시안 선형 미분 방정식을 성공적으로 도출하고, 연산자를 구조적으로 인수분해함으로써 고차 입자 기여가 타원 적분과 직접합 구조를 갖는다는 일반적인 패턴을 확인하였다. 또한, 모듈러 소수 계산과 인수분해 기법을 통해 대규모 기호 연산 없이도 정확한 정수 형태의 연산자를 복원할 수 있음을 입증하였다. 이러한 방법론은 χ⁽ⁿ⁾(n>6)에도 확장 가능하며, 이징 모델의 복잡한 특이점 구조와 임계 현상을 이해하는 데 중요한 도구가 될 것이다.