[논문 리뷰] Second order gauge invariant discretizations to the Schrödinger and Pauli equations

Snorre H. Christiansen, Tore Gunnar Halvorsen|arXiv (Cornell University)|2015. 05. 29.

Electromagnetic Simulation and Numerical Methods참고 문헌 15인용 수 3

한 줄 요약

이 논문은 유한요소와 격자 게이지 이론을 사용하여 슈뢰딩거 방정식과 파울리 방정식을 해결하기 위한 두 번째 차수 정확도를 갖는 게이지 불변 유한요소 방법을 제시한다. 게이지 포텐셜에 의존하는 유한요소 공간을 구성하고, 이를 게이지에 의존하는 노름을 통해 분석함으로써, 이 방법은 이산 게이지 불변성을 유지하면서도 두 번째 차수 수렴을 달성한다. 이는 이전의 첫 번째 차수 방법에 비해 핵심적인 개선이다.

ABSTRACT

We introduce a numerical method, based on finite elements and lattice gauge theory, to compute approximate solutions to Schrödinger and Pauli equations. The crucial geometric property of the method is discrete gauge invariance. The main new achievement is second order convergence. This is proved by interpreting the method as defined on gauge potential dependent finite element spaces and providing an analysis of such spaces in terms of gauge potential dependent norms on simplices of all dimensions.

연구 동기 및 목표

• 지역 U(1) 변환에 대해 이산 게이지 불변성을 유지하는 슈뢰딩거 및 파울리 방정식을 위한 수치적 방법을 개발하기.
• 표준 갈레르킨 유한요소 방법이 1차 수렴을 보이는 것과는 달리, 이산 H^1 및 L^2 노름에서 두 번째 차수 수렴을 달성하기.
• 위트니 형식과 노드 보간자를 사용하여 게이지 포텐셜에 의존하는 함수 공간을 다룰 수 있도록 유한요소 방법을 확장하기.
• 에너지 준위와 확률 밀도와 같은 관측 가능량이 이산 설정에서 게이지 변환에 대해 불변을 유지하도록 보장하기.
• 모든 차원의 단체에 정의된 노름을 통해 게이지 포텐셜 A에 의존하는 유한요소 공간의 엄밀한 분석을 제공하기.

제안 방법

• 이 방법은 단체 메esh 위에서 연속적인 조각별 아핀 복소값 유한요소를 사용하며, 양자 기울기(코바리언트 기울기)는 $\nabla_A u = \nabla u + i A u$로 정의된다.
• 이산 게이지 불변성은 파동함수의 게이지 변환을 $u' = \Pi_h e^{i\alpha} u$로 수정함으로써 달성되며, 여기서 $\Pi_h$는 노드 보간자로, $u'$가 유한요소 공간에 남아 있도록 보장한다.
• 변분 설정은 이차형식 $a[u,v] = \int (\nabla_A u)^\dagger \cdot (\nabla_A v)$ 와 $\langle u,v \rangle = \int u^\dagger v$를 사용하며, 이는 수정된 변환 규칙을 통해 이산 게이지 불변성이 유지된다.
• 모든 차원의 단체에 대해 게이지 포텐셜에 의존하는 노름을 사용하여 분석을 수행함으로써, 두 번째 차수 수렴의 엄밀한 증명이 가능해진다.
• 파울리 방정식의 경우, 스핀어 웨이브함수와 $\vec{\sigma} \cdot \operatorname{curl} A$ 항을 요소 위에서 조각별 상수 근사로 이산화함으로써 일致성과 게이지 불변성을 확보한다.
• 이 방법은 게이지 포텐셜 $A$를 위트니 1형식으로 근사하는 데 의존하며, 이는 일정한 자석장에서 정확한 표현을 가능하게 하여 최적의 수렴을 가능하게 한다.

실험 결과

연구 질문

• RQ1 슈뢰딩거 방정식을 위한 유한요소 방법이 이산 게이지 불변성을 유지하면서도 두 번째 차수 수렴을 달성할 수 있는가?
• RQ2 표준 유한요소 공간이 국소 U(1) 변환에 대해 불변이 아니므로, 이산 설정에서 게이지 불변성이 어떻게 유지될 수 있는가?
• RQ3 게이지 포텐셜에 의존하는 유한요소 공간과 노름이 더 높은 수렴 차수를 가능하게 하는 데 어떤 역할을 하는가?
• RQ4 이 방법은 스핀을 포함한 파울리 방정식으로 확장될 수 있으며, 게이지 불변성과 두 번째 차수 수렴을 유지하는가?
• RQ5 게이지 포텐셜 $A$를 위트니 형식으로 근사할 때, 어떤 조건에서 수렴 속도가 떨어지지 않는가?

주요 결과

• 제안된 방법은 슈뢰딩거 방정식에 대해 H^1 노름에서 두 번째 차수 수렴과 L^2 노름에서 두 번째 차수 수렴을 달성한다. 이는 표준 갈레르킨 방법이 h의 차수로 수렴하는 것과 비교해 큰 향상이다.
• 파동함수의 게이지 변환을 노드 보간을 사용하도록 수정함으로써 이산 게이지 불변성이 유지되며, 이는 변환된 해가 유한요소 공간에 남아 있도록 보장한다.
• 모든 차원의 단체에 대해 정의된 노름을 사용하여 게이지 포텐셜에 의존하는 유한요소 공간을 분석함으로써, 두 번째 차수 수렴을 증명하는 이론적 기초를 확립한다.
• 일정한 자석장에서는 $A$의 위트니 1형식에 의한 근사가 정확하므로, 포텐셜 근사에서 발생하는 추가 오차 없이 전체 두 번째 차수 수렴이 달성된다.
• 이 방법은 스핀을 포함한 파울리 방정식으로 확장 가능하며, $\vec{\sigma} \cdot \operatorname{curl} A$ 항과 스칼라 포텐셜을 조각별 상수 및 노드 근사를 통해 정확히 이산화함으로써 두 번째 차수 수렴과 게이지 불변성을 유지한다.
• 노드 보간이 사용되더라도, 적절한 함수 공간 프레임워크 내에서 분석을 수행함으로써, 이산 이차형식의 일致성이 게이지 변환 하에 유지되어 변분 범죄를 피한다.