반무한 구간을 위한 라게르‑소보레프 다항식과 스펙트럴 경계값 해법

초록

본 논문은 (0,∞) 구간에서 정의되는 라게르‑소보레프 직교다항식 군을 구축하고, 이들을 이용해 1/x 형태의 특이 전위가 포함된 슈뢰딩거형 방정식의 디리클레 경계값 문제를 선형 시스템 없이 재귀적으로 해결할 수 있는 완전 대각화 스펙트럴 방법을 제시한다. 연결식, 재귀 관계, 비동등점 근사 및 Bessel 함수를 이용한 생성함수를 도출하고, 수치 실험을 통해 지수적 수렴을 확인한다.

상세 분석

논문은 먼저 고전 라게르 다항식 Lₙ^{(α)}(x)의 기본 성질(정규화, 3항 재귀, 미분·구조 관계, 복소평면에서의 외부 비동등점 근사)을 정리한다. 이후 Sobolev 내적 ⟨u,v⟩λ = λ∫₀^∞ u(x)v(x)/x dx + ∫₀^∞ u’(x)v’(x) dx 를 도입하고, 이 내적에 대해 (0,∞)에서 x e^{-x/2}·p(x) 형태의 함수들이 사라지는 공간 P x e^{-x/2}에 대한 직교 기저 {S_n(x)}를 Gram‑Schmidt 과정으로 정의한다. 핵심은 고전 라게르 다항식 Lₙ^{(1)}(x)와 새로운 다항식 S_n(x) 사이의 단순 연결식 Lₙ^{(1)}(x)=S_n(x)+a{n-1}S_{n-1}(x) (n≥1) 를 증명한 점이다. 여기서 계수 a_n 은 재귀식 a_n = (n+2)/(4λ+2(n+1)) - n(n+1)/(4λ+2(n+1))·a_{n-1} 로 정의되며, 초기값 a_0=2/(4λ+2) 를 갖는다. 이 재귀는 a_n이 0<a_n<1을 만족함을 보이며, 명시적으로 a_n = (n+2)n! L_n^{(1)}(-4λ)·L_{n+1}^{(1)}(-4λ) 로 표현될 수 있다. 비동등점 근사를 이용하면 a_n = 1 - √(4λ)/√n + O(1/n) 로 수렴한다는 중요한 비동등점 결과를 얻는다. 또한 S_n(x)와 L_n^{(1)}(-4λ) 사이의 관계식 S_n(x)L_n^{(1)}(-4λ)^{n+1}=∑{k=0}^n (-1)^{n-k}L_k^{(1)}(x)L_k^{(1)}(-4λ)^{k+1} 를 통해 S_n의 명시적 전개를 얻고, 이를 바탕으로 lim{n→∞}√{n+1} S_n(x)L_n^{(1)}(x)=0 (복소평면에서 양의 실축 제외) 를 증명한다. 생성함수 측면에서는 ∑{n=0}^∞ S_n(x)L_n^{(1)}(-4λ)^{n+1}ω^n = (1-ω)^{-1}e^{-(x-4λ)ω/(1-ω)}(1+ω)^{1/2}J{1/4}(2√{xλω})/(1-ω)^{1/2} 형태를 도출해 Bessel 함수와의 깊은 연관성을 밝힌다. 이러한 이론적 토대를 바탕으로, 슈뢰딩거형 BVP -u’’+λ u/x = f(x), u(0)=u(∞)=0 의 해를 u(x)=∑{n=0}^∞ b_n S_n(x)e^{-x/2} 로 전개한다. 내적 정의와 적분 부분적분을 이용하면 b_n·‖S_n e^{-x/2}‖λ^2 = ∫₀^∞ f(x)S_n(x)e^{-x/2}dx 가 되므로, 선형 시스템을 풀 필요 없이 f_n = ∫ f L_n^{(1)}e^{-x/2}dx 와 연결식 g_n = f_n + a{n-1}f{n-1} 로부터 b_n을 재귀적으로 구할 수 있다. Sobolev 노름 s_n 역시 동일한 재귀식으로 얻어져 전체 알고리즘은 O(N) 복잡도로 구현 가능하며, 전통적인 라게르 스펙트럴 방법과 동일한 연산량이면서 특이 전위에 대한 수치적 안정성을 크게 향상시킨다. 마지막으로 1/x 전위를 가진 비동질 슈뢰딩거 방정식에 대해 N=30~80 정도의 차원에서 실험을 수행했으며, 오류가 exp(-cN) 형태로 급격히 감소함을 확인, 이론적 지수 수렴을 실증하였다. 전체적으로 라게르‑소보레프 다항식의 구조적 특성을 활용해 반무한 구간의 특이 BVP를 효율적으로 해결하는 새로운 스펙트럴 프레임워크를 제시한다.