반투명 경계조건을 이용한 포텐셜 분할 크랭크 니콜슨 스키마

초록

**
본 논문은 반무한 스트립 영역에서 가변 계수를 갖는 2차원 시간 의존 슈뢰딩거 방정식을 풀기 위해, 포텐셜에 대한 스트랭(Strang) 분할을 적용한 크랭크‑니콜슨 유한 차분 스키마와 이산 투명 경계조건(Discrete Transparent Boundary Conditions, DTBC)을 결합한다. 제안된 방법은 시간에 대해 무조건적인 $L^2$ 안정성을 보장하고, FFT 기반 직접 알고리즘을 통해 일반 포텐셜에 대해 효율적으로 구현한다. 수치 실험에서는 직사각형 장벽을 통한 터널 효과를 검증하고, 오류 분석을 통해 2차 정확도와 안정성을 확인한다.

**

상세 분석

**
논문은 먼저 반무한 스트립 $\Omega=(0,\infty)\times(0,Y)$ 상에서 정의된 일반화된 슈뢰딩거 방정식 $i\rho D_t\psi=(\mathcal H_0+V)\psi$ 를 소개한다. 여기서 $\mathcal H_0$ 은 비대칭이지만 양의 정칙성을 갖는 2차 미분 연산자이며, 계수 $\rho$, $B_{pq}$, $V$ 는 $x$와 $y$에 따라 변한다. 물리적 요구에 따라 경계조건은 $x=0$에서 디리클레, $x\to\infty$에서 방사형 소멸, $y$ 방향은 주기적 혹은 디리클레 조건을 적용한다.

전통적인 크랭크‑니콜슨 스키마는 암시적이며, 매 시간 단계마다 복소수 선형 시스템을 풀어야 하는데, 특히 $y$ 방향에 변수가 있는 경우 효율적인 해법이 부족했다. 이를 극복하기 위해 저자들은 포텐셜 $V$ 를 두 부분 $eV$와 $\Delta V=V-eV$ 로 분리하고, 스트랭 분할을 적용한다. 구체적으로는 (1) $\Delta V$에 대한 순수 진동 연산을 절반 시간 동안 수행하고, (2) $eV$와 $\mathcal H_0$를 포함한 전체 연산을 전체 시간 동안 수행한 뒤, (3) 다시 $\Delta V$ 연산을 절반 시간 동안 수행한다. 이 과정은 $eV$가 $x\ge X_0$에서 상수이므로, 첫 번째와 세 번째 단계는 정확히 지수 연산 $E_m=\exp!\big(-i\tau_m\Delta V/(2\hbar\rho)\big)$ 로 구현 가능해 계산 비용을 크게 낮춘다.

시간 이산화는 비균일 시간 격자를 허용하면서도 $\partial_t$를 중앙 차분으로, 평균 연산을 $s_t$ 로 정의한다. 공간 이산화는 $x$와 $y$ 각각에 비균일 격자를 두고, 중앙 차분과 평균 연산을 적절히 조합해 2차 정확도를 유지한다. 특히 $x$ 방향 무한 영역에 대해 $h_j=h$ 로 일정하게 두어, $j>J$ 구간에서는 격자 전이 연산이 일정해 DTBC를 Fourier 변환으로 명시적으로 표현할 수 있다.

DTBC는 연산자 $S_m$ 로 나타내며, 이는 시간에 대한 이산 컨볼루션과 $y$ 방향에 대한 이산 푸리에 급수의 결합 형태이다. 저자들은 $S_m$ 가 기존 연구