반대 분산·비선형성에서의 매몰 솔리톤 변분 해석

본 논문은 기본파(Fundamental Harmonic, FH)와 2배파(Second Harmonic, SH) 사이에 존재하는 반대 분산(opposing dispersion)과 경쟁적 비선형성(competing nonlinearities)을 갖는 광학 시스템에서 매몰 솔리톤(Embedded Soliton, ES)과 일반 솔리톤(ordinary soliton)의 존재와 특성을 변분법을 통해 체계적으로 분석한다. 먼저, 저자들은 기존의 복합 비선형 편미분 방정식(1)·(2)를 제시한다. 여기에는 자기위상변조(self‑phase modulation, SPM), 교차위상변조(cross‑phase modulation, XPM), 그리고 χ^(2) 비선형성 항이 모두 포함된다. 이러한 방정식은 일반적으로 수치적 방법에 의존해 해를 구하지만, 저자들은 라그랑지안 L을 정의하고 시간에 대한 적분을 수행해 효과 라그랑지안 h_Li를 도출한다. 효과 라그랑지안은 전 모델(full model)과 단순화 모델(truncated model) 모두에 적용 가능하도록 구성된다. 전 모델은 라그랑지안 형태를 그대로 유지하고, 단순화 모델은 V⁴ 항을 제외하는 등 물리적 가정을 통해 비라그랑지안 형태가 된다. 그러나 Kaup와 Malomed가 제안한 바와 같이, 전 모델의 라그랑지안을 기반으로 V⁴ 항을 제거하고 시험함수에 대한 변분을 수행하면, 두 단계 근사법 없이도 동일한 결과를 얻을 수 있음을 저자들은 증명한다. 시험함수는 FH와 SH의 시간 의존성을 각각 U(t)=A sech(√2k t)와 V(t)=B sech²(√2k t) 형태로 가정한다. 여기서 A와 B는 변분 파라미터이며, 폭 파라미터 √2k는 고정한다. 이 시험함수를 효과 라그랑지안에 대입하고 시간 적분을 수행하면, 파라미터만 남는 식(14)이 얻어진다. 이를 A와 B에 대해 변분하면 두 비선형 대수식(15)·(16)이 도출된다. 전 모델에 대한 해는 (15)·(16)을 직접 풀어 A와 B를 구한다. 저자들은 구체적인 파라미터(k=0.25, γ₁=−0.05, γ₂=−0.025, δ=1, q=1)를 설정하고, 전 모델에서 B₁≈0.48, B₂≈5.69, B₃≈19.44를 얻는다. 물리적으로 의미 있는 해는 B₂에 대응하는 A≈7.21이며, 이는 전 모델에서의 솔리톤 진폭과 폭을 나타낸다. 반면, 단순화 모델에서는 V≪U라는 가정을 먼저 적용해 B²·B³ 항을 무시하고 (17)·(18)을 얻는다. 동일한 파라미터 조건에서 B≈6.72, A≈11.16을 얻으며, 이는 전 모델보다 더 큰 진폭과 더 좁은 폭을 보인다. 매몰 솔리톤에 대해서도 동일한 절차를 적용한다. 매몰 솔리톤은 연속 스펙트럼 안에 존재하는 특수한 해로, 파라미터 영역(11)에서 k와 q의 관계가 반대가 된다. 저자들은 k=0.6963, γ₁=−0.05, γ₂=−0.025, δ=1, q=1을 사용해 전 모델에서 B≈4.945, A≈5.65, 단순화 모델에서 B≈6.38, A≈9.99를 얻는다. 두 경우 모두 단순화 모델의 솔리톤이 더 뾰족하고 빠르게 감소하는 형태를 띤다. 에너지 정의 E=∫_{-∞}^{∞}(|u|²+2|v|²)dt에 따라 전 모델과 단순화 모델의 솔리톤 에너지를 계산하면, 항상 단순화 모델(E_T)이 전 모델(E_F)보다 크다. 이는 V가 약해질수록 FH의 기여가 상대적으로 커져 전체 에너지가 증가한다는 물리적 직관과 일치한다. 또한 매몰 솔리톤의 경우에도 동일한 경향이 관찰되어, V≪U 조건이 솔리톤 에너지 증폭에 기여한다는 점을 강조한다. 논문의 마지막에서는 이러한 변분 접근법이 전 모델과 단순화 모델 사이의 정량적 차이를 명확히 드러내며, 실험적 설계에 있어 V≪U 조건이 솔리톤 에너지와 전파 안정성을 향상시킬 수 있음을 제시한다. 또한, 라그랑지안 기반 변분법이 복잡한 비선형 파동 방정식의 해를 근사적으로라도 분석적으로 다룰 수 있는 강력한 도구임을 재확인한다.