[논문 리뷰] Transparent Boundary Conditions for the Time-Dependent Schrödinger Equation with a Vector Potential

Jason Kaye, Leslie Greengard|arXiv (Cornell University)|2018. 12. 11.

Electromagnetic Simulation and Numerical Methods참고 문헌 30인용 수 7

한 줄 요약

이 논문은 임의의 조각별로 미분 가능한 도메인에서 벡터 포텐셜을 가진 시간에 의존하는 슈뢰딩거 방정식에 대해 잠재이론과 그린의 표현 정리(Green's representation theorem)를 사용하여 정확한 투명 경계 조건(TBCs)을 유도하는 일반적인 방법을 제시한다. 1차원 문제에서 벡터 포텐셜이 0이 아닐 경우, TBCs는 시간에 대해 비컨볼루션 형태가 되며, 저자들은 이러한 비국소 조건을 효율적으로 계산하기 위해 빠른 버터플라이 알고리즘을 개발하여 최소한의 경계 반사로 정확하고 효율적인 수치 시뮬레이션을 가능하게 한다.

ABSTRACT

We consider the problem of constructing transparent boundary conditions for the time-dependent Schrödinger equation with a compactly supported binding potential and, if desired, a spatially uniform, time-dependent electromagnetic vector potential. Such conditions prevent nonphysical boundary effects from corrupting a numerical solution in a bounded computational domain. We use ideas from potential theory to build exact nonlocal conditions for arbitrary piecewise-smooth domains. These generalize the standard Dirichlet-to-Neumann and Neumann-to-Dirichlet maps known for the equation in one dimension without a vector potential. When the vector potential is included, the condition becomes non-convolutional in time. For the one-dimensional problem, we propose a simple discretization scheme and a fast algorithm to accelerate the evaluation of the boundary condition.

연구 동기 및 목표

• 임의의 조각별로 미분 가능한 도메인에서 컴팩트하게 지정된 포텐셜과 시간에 의존하는 벡터 포텐셜을 가진 시간에 의존하는 슈뢰딩거 방정식에 대해 정확한 투명 경계 조건(TBCs)을 개발하는 것.
• 이전에 스칼라 포텐셜 또는 1차원 사례에 국한된 TBC 프레임워크를 고차원에서 전자기 벡터 포텐셜을 포함하도록 확장하는 것.
• 벡터 포텐셜에 의해 유도되는 비컨볼루션 시간 의존성으로 인한 계산적 과제를 해결하기 위해 효율적인 구현을 위한 빠른 알고리즘을 설계하는 것.
• 크랭크-니콜슨 스킴을 사용한 수치 실험을 통해 TBCs와 딜리클레 경계 조건을 모두 적용한 경우의 정확성과 효율성을 입증하는 것.

제안 방법

• 임의의 도메인에서 벡터 포텐셜이 있는 슈뢰딩거 방정식에 대해 정확한 비국소 TBCs를 도출하기 위해 잠재이론과 그린의 표현 정리를 사용한다.
• TBCs는 경계 위의 적분을 포함하는 딜리클레-투-노이만 유형의 맵으로 표현되며, 커널 함수는 문제의 전체 그린 함수로부터 유도된다.
• 벡터 포텐셜이 0이 아닌 1차원 문제에서는 유도된 TBCs가 시간에 대해 비컨볼루션 형태가 되며, 이는 벡터 포텐셜이 없는 경우의 표준 TBCs와 구별된다.
• 비컨볼루션 시간 적분 연산자를 효율적으로 계산하기 위해 새로운 빠른 버터플라이 알고리즘을 제안한다. 이 알고리즘을 통해 시간 단계 수와 거의 선형적인 스케일링을 달성할 수 있다.
• 반디케이션 조건을 가진 반연속적인 크랭크-니콜슨 스킴에 TBCs를 구현하며, 공간 이산화에 고차수 유한차분 방법을 사용한다.
• 행렬 압축 기법과 사전 계산을 활용하여 메모리 및 계산 비용을 감소시키며, 특히 벡터 포텐셜의 시간적 구조가 단순할 경우 효과적이다.

실험 결과

연구 질문

• RQ1 임의의 조각별로 미분 가능한 도메인에서 벡터 포텐셜이 있는 시간에 의존하는 슈뢰딩거 방정식에 대해 정확한 투명 경계 조건을 어떻게 도출할 수 있는가?
• RQ2 벡터 포텐셜이 0이 아닐 경우 TBCs의 수학적 구조는 어떻게 되며, 이는 벡터 포텐셜이 없는 경우의 표준 컨볼루션 형태와 어떻게 다를까?
• RQ3 벡터 포텐셜 사례에서 나타나는 비컨볼루션 시간 적분 연산자를 효율적으로 계산하기 위한 빠른 알고리즘을 개발할 수 있는가?
• RQ4 TBC 기반 스킴의 정확도는 동일한 공간 이산화 조건을 사용할 때 표준 딜리클레 경계 조건과 비교해 어떻게 되는가?
• RQ5 TBCs를 적용하는 데 드는 계산 비용은 얼마이며, 현대의 빠른 알고리즘을 사용해 대규모 시뮬레이션에 대해 확장 가능한가?

주요 결과

• 제안된 TBCs는 임의의 조각별로 미분 가능한 도메인에서 벡터 포텐셜이 있는 슈뢰딩거 방정식에 대해 정확하며, 그린의 표현 정리를 통해 유도되었다.
• 벡터 포텐셜이 0이 아닐 경우 TBCs는 시간에 대해 비컨볼루션 형태가 되며, 이는 표준 컨볼루션 기반 알고리즘을 무효화하고 새로운 계산 전략이 필요하게 한다.
• 1차원 문제에서는 시간 단계 수와 거의 선형적인 스케일링을 달성하는 빠른 버터플라이 알고리즘을 개발하여 TBC 계산 비용을 총 시뮬레이션 시간의 0.02% 미만으로 줄였다.
• 수치 실험 결과, 전체 스킴에서 시간에 대해 2차 수렴을 보였으며, 동일한 공간 이산화 조건을 사용할 경우 TBC 및 딜리클레 경계 조건 스킴 간의 오차는 거의 동일했다.
• 벡터 포텐셜의 부호가 변화하는 경우, 외부 영역에서의 오차 누적가 줄어들기 때문에 TBC 스킴이 딜리클레 경계 조건 스킴보다 더 정확했다.
• 벡터 포텐셜의 시간적 의존성이 단순할 경우, 시스템 행렬 $S_N$ 및 $D_N$의 행렬 압축이 효과적이며, 효율적인 사전 계산과 저장을 가능하게 했다.