스핀 1 BEC 밝은 솔리톤의 섭동 이론과 안정성 분석

본 논문은 스핀‑1 보스‑아인슈타인 응축(BEC)에서 나타나는 밝은 솔리톤의 동역학을 적분 가능한 모델을 기반으로 섭동 이론을 구축한다. 서론에서는 1차원 BEC에서 관측된 밝고 어두운 솔리톤의 실험적 배경과, 스핀 자유도가 추가된 스핀‑1 BEC가 벡터 솔리톤을 지원한다는 점을 소개한다. 기존 연구에서 적분 가능한 Yajima‑Oikawa 및 행렬 NLS 모델이 스핀‑1 BEC에 적용 가능함을 언급하고, 실제 실험에서는 적분 가능성 조건(c₀=c₂)이 완벽히 만족되지 않으므로, 이를 작은 섭동으로 취급하는 필요성을 제시한다. II절에서는 기본 모델을 제시한다. 세 개의 스핀 성분(Φ₊, Φ₀, Φ₋)으로 구성된 파동함수를 정의하고, Gross‑Pitaevskii 방정식 형태의 연동 방정식(2.1)을 제시한다. 스핀‑독립 상수 c₀와 스핀‑의존 상수 c₂를 도입하고, 적분 가능성을 보장하는 제약 c₀=c₂≡−c<0을 적용한다. 좌표와 시간의 정규화를 통해 차원 없는 형태(2.3)로 변환하고, 이를 2×2 행렬 Q로 재구성하여 행렬 NLS 방정식(2.5)로 정리한다. III절에서는 스펙트럼 문제를 위한 Jost 해와 스캐터링 행렬을 정의하고, 복소 평면에서의 해석적 성질을 상세히 논의한다. Jost 해는 Volterra 적분 방정식 형태로 전개되며, 상반부(C⁺)와 하반부(C⁻)에서 각각 다른 열/행이 해석적임을 보인다. 스캐터링 행렬은 상하 삼각 형태로 분해되어, RH 문제의 핵심 데이터가 된다. IV절에서는 두 해 ψ⁺와 ψ⁻¹을 이용해 리만‑히루베르트(RH) 문제를 구성한다. 실선에서의 점프 조건 ψ⁻¹ψ⁺=EGE⁻¹(4.1)와 정규화 조건 ψ⁺,ψ⁻¹→I(4.3)를 제시한다. RH 문제의 해는 스펙트럼 데이터(영점)의 존재 여부에 따라 정규와 비정규(영점 포함)로 나뉜다. 영점은 솔리톤을 나타내며, 영점의 기하학적·대수적 중복도에 따라 랭크‑1(단일 고유벡터)와 랭크‑2(두 개의 고유벡터) 솔리톤이 구분된다. V절에서는 실제 솔리톤 해를 도출한다. 랭크‑1 경우에는 전통적인 하이퍼볼릭 시컨트 형태의 단일 피크 솔리톤이 얻어지며, 이는 강자성(ferromagnetic) 상태와 동일함을 증명한다. 랭크‑2 경우에는 두 종류가 존재한다. 첫 번째는 두 개의 피크가 겹쳐 강자성 솔리톤과 동등함을 보이며, 두 번째는 피크가 분리된 극성(polar) 솔리톤이다. 극성 솔리톤은 총 스핀이 0이며 에너지가 더 높고, 프로파일이 이중 피크 형태를 가진다. 수치 실험에서 극성 솔리톤은 섭동이 가해지면 두 개의 강자성 솔리톤으로 분열함을 확인한다. 따라서 전체 섭동 이론은 강자성(랭크‑1) 솔리톤에만 집중하면 충분하다. VI절에서는 섭동 이론을 전개한다. 적분 가능성 조건의 작은 위반을 ε·δc 형태의 섭동으로 모델링하고, RH 프레임워크에서 파라미터(속도, 진폭, 위상, 스핀 성분 비율)의 진화 방정식을 도출한다. adiabatic 근사하에 시간에 따라 변하는 파라미터를 구하고, 섭동이 주로 솔리톤의 위상(주파수)만을 미세하게 이동시킴을 확인한다. 속도와 진폭, 스핀 비율은 1차 섭동에서 보존된다. VII절에서는 수치 시뮬레이션을 통해 이론을 검증한다. 초기 조건으로 강자성 및 극성 솔리톤을 설정하고, c₀≠c₂인 경우를 시뮬레이션한다. 결과는 강자성 솔리톤이 거의 변형 없이 이동하고, 주파수만 약간 이동함을 보여준다. 극성 솔리톤은 두 개의 강자성 솔리톤으로 분열하고, 이후 각각 위와 동일한 동역학을 보인다. 방출되는 복사와 프로파일 변형은 매우 작아, 섭동에 대한 솔리톤의 강인성을 입증한다. 결론에서는 본 섭동 이론이 행렬 NLS와 스핀‑1 BEC 시스템에 대해 강력한 분석 도구임을 강조한다. RH 문제 기반 접근법은 기존 Gel'fand‑Levitan 방식보다 대수적으로 간결하며, 랭크‑1과 랭크‑2 솔리톤을 통합적으로 다룰 수 있다. 향후 연구에서는 다차원 스핀‑1 BEC, 외부 포텐셜, 그리고 비선형 파라미터 변조에 대한 확장 가능성을 제시한다.