시간 의존 슈뢰딩거 방정식의 흡수 경계와 소스 항을 위한 간편 크랭크 니콜슨 구현

초록

본 논문은 크랭크‑니콜슨 방법과 카이리 형태의 유한 차분 연산자를 결합해, 위치표현에서 진화 연산자 행렬을 직접 해석적으로 역전시킨다. 이를 통해 매트리카 6.0 및 유사 과학 계산 환경에서 소스 항을 포함한 시간 의존 슈뢰딩거 방정식을 손쉽게 구현할 수 있는 알고리즘을 제시한다.

상세 요약

시간 의존 슈뢰딩거 방정식( TDSE )은 양자역학 시뮬레이션에서 핵심적인 역할을 하지만, 수치적 해법은 두 가지 큰 난관에 직면한다. 첫 번째는 행렬 연산을 최소화하면서도 높은 정확도를 유지하는 방법이며, 두 번째는 외부 구동(소스) 항을 자연스럽게 포함시킬 수 있는 구조이다. 기존의 분할-운동법(split‑operator)이나 피스톤-스펙트럼 방법은 행렬 역연산을 회피하지만, 비선형 소스 항을 도입할 경우 복잡도가 급증한다. 반면, 크랭크‑니콜슨(CN) 방식은 단위 연산자를 보존하면서도 행렬 방정식 형태로 전개되지만, 매트리카와 같은 고수준 언어에서는 N×N 행렬을 직접 역산하는 것이 비효율적이다.

저자들은 이러한 딜레마를 해소하기 위해, 1차원 격자에서 라플라시안 연산자를 3‑대각 행렬 형태로 이산화하고, 연산자 exp(−iĤΔt/ħ)를 카이리 변환 (1−iĤΔt/2ħ)⁻¹(1+iĤΔt/2ħ) 로 근사한다. 핵심 아이디어는 이 3‑대각 행렬의 역을 폐쇄형으로 구할 수 있다는 점이다. 저자는 연속적인 가우시안 소거와 역행렬의 재귀 관계를 이용해, N×N 행렬의 모든 원소를 O(N) 연산으로 표현한다. 결과적으로, 매트리카 내부에서 Inverse 함수를 호출하지 않아도, 직접 계산된 계수 배열만으로 다음 시간 단계의 파동함수를 업데이트할 수 있다.

또한, 소스 항 S(x,t) 를 RHS에 추가하면, CN 스킴은 선형 시스템 A ψ^{n+1}=B ψ^{n}+Δt·S^{n+½} 로 변형된다. 여기서 A와 B는 앞서 정의한 카이리 연산자의 분자·분모 행렬이며, S^{n+½}는 시간 중앙값으로 보간한다. 저자는 소스 항을 포함한 경우에도 동일한 폐쇄형 역행렬을 재사용할 수 있음을 증명한다. 흡수 경계 조건(ABC)은 복소수 포텐셜를 격자 가장자리 영역에 삽입함으로써 구현되며, 이는 파동이 인공적인 반사 없이 외부로 “흡수”되도록 만든다. 이러한 ABC는 행렬 구조에 영향을 주지 않으므로, 기존의 해석적 역행렬 공식에 그대로 적용된다.

수치 실험에서는 Gaussian 파킷이 자유 입자와 퍼텐셜 장벽을 통과하는 과정을 시뮬레이션하고, 소스 항을 이용해 지속적인 파킷 방출을 모델링한다. 결과는 에너지 보존, 반사 최소화, 그리고 시간 단계에 대한 2차 정확도를 만족한다. 특히, 매트리카 6.0에서 구현된 코드는 1000 격자점 기준으로도 0.1초 이내에 전체 시뮬레이션을 완료한다는 점에서, 기존의 LinearSolve 기반 구현보다 5배 이상 빠른 성능을 보인다.

이 논문의 주요 공헌은 (1) 카이리 형태와 3‑대각 행렬의 폐쇄형 역을 제공함으로써, CN 스킴을 “매트리카 친화적”으로 만든 점, (2) 소스 항과 흡수 경계 조건을 자연스럽게 통합한 일반화된 프레임워크, (3) 구현 코드가 간결하고, 다른 고수준 언어(예: Python/NumPy, MATLAB)에도 그대로 이식 가능하도록 설계한 점이다. 이러한 접근은 양자 동역학 시뮬레이션을 처음 접하는 비전문가에게도 실용적인 도구를 제공한다는 점에서 학문적·교육적 가치가 크다.