다중모드 스핀보존 모델의 대칭과 정확해

본 연구는 스핀‑보존 모델, 즉 두 수준 시스템(스핀)과 다중 보존 모드가 상호작용하는 일반적인 해밀토니언을 대상으로, 그 대칭 구조를 체계적으로 규명하고 정확한 에너지 스펙트럼을 도출하는 방법을 제시한다. 첫 번째 단계에서는 기본 해밀토니언 \(H_M=\sum_{j=1}^N\omega_j a_j^\dagger a_j+\sigma_x\sum_{k=1}^N g_k(a_k^\dagger+a_k)+\Delta\sigma_z\) 의 대칭을 분석한다. 여기서 알려진 Z₂ 짝수 대칭(스핀·보존 전체의 짝수성) 외에도, 시간반전 연산자 \(T\)와 그와 짝수 연산자 \(\Pi\)의 조합이 보존됨을 확인한다. 특히, 스핀 X축 회전을 보존수 연산자 \(P\)와 결합한 유니터리 변환 \(U=e^{i\frac{\pi}{2}\sigma_x P}\) 을 도입함으로써, 스핀 부분을 완전 대각화하고 보존 연산자 \(\Pi=\sigma_z P\)가 스핀‑보존 상호작용을 단순히 부호 전환으로만 작용하도록 만든다. 변환 후 해밀토니언은 \(U^\dagger H_M U=\sum_j\omega_j a_j^\dagger a_j-i\sum_k g_k(a_k^\dagger-a_k)+\Delta\Pi\) 의 형태가 되며, 이는 스핀과 보존이 분리된 구조를 제공한다. 다음으로, 보존 모드가 2배 광자(두 포톤)와 결합하는 확장형 \(H_{MS}=\sum_j\omega_j a_j^\dagger a_j+\sigma_x\sum_k g_k(a_k^{\dagger2}+a_k^2)+\Delta\sigma_z\) 을 고려한다. 여기서는 \(\theta=\pi/4\) 회전을 적용한 변환 \(\tilde U=e^{i\frac{\pi}{4}\sigma_x P}\) 를 사용해 \(\tilde U^\dagger H_{MS}\tilde U=\sum_j\omega_j a_j^\dagger a_j-i\sum_k g_k(a_k^{\dagger2}-a_k^2)+\Delta\Pi_q\) 을 얻는다. \(\Pi_q\)는 Z₄ 대칭을 구현하는 연산자로, \(\sqrt{P}\)와 \(\sigma_{y,z}\)의 조합으로 정의된다. 이 변환은 전체 힐베르트 공간을 네 개의 불변 부분공간(짝수·홀수 보존수와 스핀 양·음)의 직교합으로 분해시키며, 각 부분공간에서 스핀‑보존 결합이 완전히 해제된 형태의 유효 해밀토니언을 만든다. 정확해를 구하기 위해 다변량 Bargmann 공간을 도입한다. 보존 연산자를 복소 변수와 미분 연산자로 치환하면, 스핀 대각화 후 두 블록(±)은 다음과 같은 미분 연산자로 표현된다. \(\tilde H_\pm =\sum_j\omega_j z_j\partial_{z_j}+ \sum_k g_k(z_k+\partial_{z_k})\pm\Delta T\) 여기서 \(T\)는 모든 변수 \(z_j\)를 \(-z_j\)로 반전시키는 연산자이다. 스크루딩어 방정식은 \((\omega\cdot z+g)\cdot\partial_z\psi(z)+(g\cdot z-E)\psi(z)\pm\Delta\psi(-z)=0\) 이며, 이는 Braak이 제시한 G‑함수 방법을 그대로 적용할 수 있다. 다중 인덱스 \(\mathbf n=(n_1,\dots,n_N)\)에 대한 급수 전개와 재귀 관계식(10)를 통해 계수 \(A_{\mathbf n}\)를 구하고, 최종적으로 \(G_\pm^{(N)}(X)=\sum_{\mathbf n\in\mathbb N^N}\frac{A_{\mathbf n}}{X-\omega\cdot\mathbf n}\) 의 영점이 에너지 스펙트럼을 결정한다. 여기서 \(X=E+\sum_j g_j^2/\omega_j\)이며, \(\omega\cdot\mathbf n\)는 보존 모드의 자유 진동수와 정수 인덱스의 내적이다. 수치 검증으로 N=2인 두 모드 사례를 선택하였다. 파라미터 \(\omega_1=1,\ \omega_2=0.92,\ g_1=0.7,\ g_2=0.78,\ \Delta=0.68\)에 대해 G‑함수 \(G_\pm^{(2)}(X)\)를 계산하고, 그 극점(보존 모드가 결합이 없을 때의 에너지)과 영점을 플롯하였다. 영점은 각 극점 근처에 위치해 단일 모드 Rabi 모델과 동일한 구조를 보이며, 이는 제시된 해법이 다중 모드에서도 정확히 작동함을 입증한다. 또한 \(g_1,g_2\) 평면에서 첫 번째 에너지 레벨의 등고선을 그려, 두 모드 간 교환 대칭(\(g_1\leftrightarrow g_2\))과 강결합 영역에서의 스펙트럼 변화를 시각화하였다. 결론적으로, 본 논문은 (1) 시간반전·Z₂·보존 교환 대칭을 동시에 활용한 새로운 스핀‑보존 해밀토니언 대각화 기법, (2) Z₄ 대칭을 포함한 다중 보존 모드의 정확해 구조, (3) 다변량 Bargmann 공간과 일반화된 G‑함수를 결합한 전역적인 해석적 해법을 제시한다. 이는 기존 단일 모드 Rabi 모델의 정확해를 다중 모드로 확장하는 중요한 이론적 토대를 제공하며, 고차원 힐베르트 공간에서의 스핀‑보존 상호작용을 이해하고 설계하는 데 필수적인 도구가 될 것이다.