[논문 리뷰] Some remarks on discrete and semi-discrete transparent boundary conditions for solving the time-dependent Schrödinger equation on the half-axis

Alexander Zlotnik, Ilya Zlotnik|arXiv (Cornell University)|2014. 06. 19.

Electromagnetic Simulation and Numerical Methods참고 문헌 1인용 수 5

한 줄 요약

이 논문은 유한차분 스킴을 사용하여 半축 위에서 시간에 의존하는 슈뢰딩거 방정식을 풀기 위해 이산적이고 반연속적인 투명경계조건(TBCs)에 대한 새로운 분석을 제시한다. 이는 이산적 TBC와 반연속적 TBC의 합성 커널 간의 차이에 대한 첫 번째 경계를 증명하며, 이산적 TBC를 반연속적 TBC로 대체할 경우 절대 오차가 매우 작아짐을 보여주며, 특히 향상된 반연속적 TBC(ISDTBC)를 사용할 경우, $L^2$-노름에서 해의 차이가 $4.42\cdot10^{-5}$ 이하로 줄어든다.

ABSTRACT

We consider the generalized time-dependent Schrödinger equation on the half-axis and a broad family of finite-difference schemes with the discrete transparent boundary conditions (TBCs) to solve it. We first rewrite the discrete TBCs in a simplified form explicit in space step $h$. Next, for a selected scheme of the family, we discover that the discrete convolution in time in the discrete TBC does not depend on $h$ and, moreover, it coincides with the corresponding convolution in the semi-discrete TBC rewritten similarly. This allows us to prove the bound for the difference between the kernels of the discrete convolutions in the discrete and semi-discrete TBCs (for the first time). Numerical experiments on replacing the discrete TBC convolutions by the semi-discrete one exhibit truly small absolute errors though not relative ones in general. The suitable discretization in space of the semi-discrete TBC for the higher-order Numerov scheme is also discussed.

연구 동기 및 목표

• 반축 위에서 시간에 의존하는 슈뢰딩거 방정식을 풀이할 때 이산적 TBC와 반연속적 TBC 간의 관계를 분석하는 것.
• 공간 스텝 $h$에 따라 표현되는 형태로 이산적 TBC를 단순화하는 것.
• 일부 스킴에 대해 TBC 내 이산적 합성의 $h$에 독립성과 반연속적 TBC 합성과의 일치 여부를 조사하는 것.
• 절대 오차 및 상대 오차 측면에서 이산적 TBC를 반연속적 TBC로 대체했을 때의 수치 성능 평가.
• 고차수 Numerov 스킴에 적합한 반연속적 TBC의 적절한 공간 이산화를 제안하고 분석하는 것.

제안 방법

• 공간 스텝 $h$에 의존하는 간단하고 명시적인 형태로 이산적 TBC를 재작성하는 것.
• 이 가족 내 특정 스킴($\theta = \frac{1}{12}$)에서 시간에 대한 이산적 합성이 $h$에 독립적이며 반연속적 TBC 합성과 정확히 일치함을 규명하는 것.
• 이산적 TBC와 반연속적 TBC 합성 커널 간의 차이에 대한 이론적 경계를 증명하는 것 — 이는 처음으로 확립된 결과이다.
• 수치 실험을 통해 이산적 TBC, 반연속적 TBC, 그리고 향상된 반연속적 TBC(ISDTBC) 간의 $L^2$ 및 $C$-노름 오차를 비교하는 것.
• 공간 평균 매개변수 $\theta$를 가진 이중단계 대칭 유한차분 스킴(크랭크-니콜슨 유형)의 가족을 사용하는 것.
• 선형화 방법을 적용하여 반연속적 TBC를 유도하고, 이를 공간적으로 이산화하여 고차수 스킴용 ISDTBC를 구성하는 것.

실험 결과

연구 질문

• RQ1 선택된 스킴에 대해 이산적 TBC의 이산적 합성은 공간 스텝 $h$에 의존하지 않는가?
• RQ2 이산적 TBC 합성 커널은 반연속적 TBC 합성 커널에 대해 경계를 가질 수 있는가?
• RQ3 이산적 TBC를 사용한 수치적 해와 반연속적 TBC를 사용한 해 사이의 절대 오차 및 상대 오차는 어떻게 비교되는가?
• RQ4 향상된 반연속적 TBC(ISDTBC)를 사용할 경우 정확도와 계산 비용에 어떤 영향을 미치는가?
• RQ5 고차수 Numerov 스킴은 적절히 이산화된 반연속적 TBC와 효과적으로 결합될 수 있는가?

주요 결과

• 스킴에서 $\theta = \frac{1}{12}$일 경우, 이산적 TBC 합성은 $h$에 독립적이며 반연속적 TBC 합성과 정확히 일치하며, 이는 두 커널 간의 차이에 대한 첫 번째 이론적 경계를 가능하게 한다.
• 수치 실험 결과, 이산적 TBC를 반연속적 TBC로 대체하면 절대 오차가 매우 작아지며, 해 간의 최대 $L^2$-노름 차이는 $4.42 \cdot 10^{-5}$ 이하, $C$-노름 차이는 $8.87 \cdot 10^{-5}$ 이하로 나타난다.
• 향상된 반연속적 TBC(ISDTBC)는 상대 오차를 크게 감소시키며, 특히 $C$-노름에서 두드러진다. 이는 $\theta = \frac{1}{12}$ 스킴이 모든 테스트된 구성 중에서 가장 정확한 것으로 나타났다.
• $\theta = \frac{1}{12}$ 이며 $J=3200$일 경우, 이산적 TBC 및 ISDTBC 모두 절대 오차 감소 순서가 약 2에 가까운 반면, 표준 반연속적 TBC의 경우 중간 수준의 $M$ 이후 상대 오차가 정체되어 수렴 개선이 제한됨을 나타낸다.
• 알고리즘의 런타임은 $M$에 거의 비례하며, TBC 커널 계산 비용은 선형 시스템을 푸는 데 비해 무시할 수 있을 정도로 낮아, ISDTBC 접근법은 계산적으로 효율적이다.
• ISDTBC를 사용한 수치적 해는 정확한 해와 비교해 매우 유사하며, 이는 적절히 구현된 반연속적 접근법이 높은 정밀도를 유지함을 확인한다.