[논문 리뷰] An Efficient and Accurate Grid Method for Solving the Time-Dependent Schroedinger Equation: Application of Coulomb Wave Function DVR to Atomic Systems in Strong Laser Fields

Liang-You Peng, Anthony F. Starace|CERN Bulletin|2006. 04. 21.

Laser-Matter Interactions and Applications인용 수 30

한 줄 요약

이 논문은 강한 레이저장에 노출된 원자 시스템에서 시간에 의존하는 슈뢰딩거 방정식을 해결하기 위한 새로운 격자 방법인 쿠론트 웨이브 함수 이산 변수 표현(CWDVR)을 제시한다. 유한 차분 대신 쿠론트 웨이브 함수를 사용하여 근방 좌표를 이산화함으로써, 쿠론트 특이점과 연속 상태를 정확히 다루면서도 격자 점 수를 3~10배 줄일 수 있다. 이로 인해 수소 기반 이온(H 및 H–)의 이온화율은 기존 이론 결과와 뛰어난 일치를 보인다.

ABSTRACT

We present an efficient and accurate grid method for solving the time-dependent Schrödinger equation of atomic systems interacting with intense laser pulses. As usual, the angular part of the wave function is expanded in terms of spherical harmonics. Instead of the usual finite difference (FD) scheme, the radial coordinate is discretized using the discrete variable representation which is constructed from the Coulomb wave function. For an accurate description of the ionization dynamics of atomic systems, the Coulomb wave function discrete variable representation (CWDVR) method needs 3-10 times less grid points than the FD method. The resultant grid points of CWDVR distribute unevenly so that one has finer grid near the origin and coarser one at larger distances. The other important advantage of the CWDVR method is that it treats the Coulomb singularity accurately and gives a good representation of continuum wave functions. The time propagation of the wave function is implemented using the well-known Arnoldi method. As examples, the present method is applied to the multiphoton ionization of both H and H$^-$ in intense laser fields. Short-time excitation and ionization dynamics of H by static electric fields is also investigated. For a wide range of photon energies and laser intensities, ionization rates calculated using this method are in excellent agreement with those from other theoretical calculations.

연구 동기 및 목표

• 강한 레이저 펄스에 노출된 원자 시스템에서 시간에 의존하는 슈뢰딩거 방정식(TDSE)을 더 효율적이고 정확하게 해결하기 위한 방법을 개발하기 위해.
• 유한 차분(FD) 기반 방법이 쿠론트 특이점과 연속 상태 파동함수를 표현하는 데 가지는 한계를 극복하기 위해.
• 수렴에 필요한 균일한 반경 좌표 격자 점 수를 최소화하여 계산 비용을 줄이기 위해.
• 강한 레이저장과 정적 전기장 하에서 수소 기반 시스템의 다광자 이온화 및 여기 동역학을 정확하게 시뮬레이션할 수 있도록 하기 위해.
• 실제 쿠론트 포텐셜을 고려한 다전자 원자 시스템에 대한 아비니타 시뮬레이션을 확장할 수 있는 기반을 마련하기 위해.

제안 방법

• 양성 에너지 쿠론트 웨이브 함수로부터 구성된 이산 변수 표현(DVR)을 사용하여 균일한 반경 좌표를 이산화함으로써 CWDVR 방법을 구성한다.
• CWDVR 격자는 본질적으로 균일하지 않으며, 핵 부근에서는 밀도가 높고, 더 큰 r에서는 간격이 넓어져 쿠론트 포텐셜의 중요도 변화를 반영한다.
• 이 방법은 자연스럽고 정확하게 쿠론트 특이점을 다루며, 이온화 동역학에 있어 핵심적인 연속 상태를 잘 표현한다.
• 파동함수의 시간 진화는 Krylov 부분공간 전파를 위한 아르놀디 방법을 사용하여 구현한다.
• 이 접근법은 단일 활성전자(SAE) 근사를 사용하며, 적절한 모델 포텐셜을 통해 중성 원자(H)와 음이온(H–) 모두에 적용 가능하다.
• H–의 경우, 핵의 극성 효과를 반영하기 위해 수정된 형태인 D = [1 - α_d / r^3 * W_3(r/r_c)] r을 사용한 디폴드 연산자를 적용한다.

실험 결과

연구 질문

• RQ1 쿠론트 웨이브 함수 DVR 기반 격자 방법이 강한 레이저장에서 TDSE를 해결하는 데 있어 표준 유한 차분 기반 방법보다 정확도와 효율성 면에서 뛰어나게 작용할 수 있는가?
• RQ2 기존 방법과 비교해 CWDVR 방법은 쿠론트 특이점과 연속 상태 파동함수를 어떻게 다루는가?
• RQ3 CWDVR 방법은 H 및 H–에서 수렴된 이온화율을 확보하기 위해 필요한 균일한 반경 격자 점 수를 얼마나 줄일 수 있는가?
• RQ4 다양한 레이저 강도와 광자 에너지 범위에서 H 및 H–의 다광자 이온화에 대해 CWDVR 방법이 알려진 이온화율을 얼마나 잘 재현하는가?
• RQ5 CWDVR 방법은 정적 전기장 하에서 H의 짧은 시간 동안의 여기 및 이온화 동역학을 정확하게 기술할 수 있는가?

주요 결과

• CWDVR 방법은 H 및 H–에서 이온화 동역학의 수렴 결과를 확보하기 위해 유한 차분(FD) 방법보다 3~10배 적은 균일한 반경 격자 점 수가 필요하다.
• CWDVR를 사용해 계산한 H 및 H–의 이온화율은 복소 회전 방법 및 비섭동 다전자 이론과 같은 다른 이론 계산 결과와 뛰어난 일치를 보인다.
• 정적 전기장 하에서 H의 경우, CWDVR로 계산한 기본 상태 생존 확률과 이온화율은 복소 회전 방법 결과와 높은 정밀도로 일치한다.
• CWDVR 격자에서 H–의 기본 상태 에너지는 -0.027730 a.u.로 측정되었으며, Telnov와 Chu의 계산 결과인 -0.027733 a.u.와 실험값인 -0.027716 a.u.와 매우 유사한 결과를 보였다.
• 레이저 강도 1×10^11 W/cm², 파장 1064 nm 조건에서 H–의 탈리프율은 4.43×10^(-4) a.u.로 이전 이론 결과와 양호한 일치를 보였다.
• 디폴드 연산자에 핵의 극성 효과를 포함시킨 결과는 표준 r 연산자를 사용한 결과보다 비섭동 MEMPT 계산 결과에 더 가까운 값을 보였으며, 물리적 정확도 향상이 확인되었다.