Generalized Pseudospectral Method for Solving the Time-Dependent Schrödinger Equation Involving the Coulomb Potential
ZENG Si-Liang1, ZOU Shi-Yang1, YAN Jun1,2
1The Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, PO Box 8009, Beijing 1000882Center for Applied Physics and Technology, Peking University, Beijing 100871
Generalized Pseudospectral Method for Solving the Time-Dependent Schrödinger Equation Involving the Coulomb Potential
ZENG Si-Liang1, ZOU Shi-Yang1, YAN Jun1,2
1The Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, PO Box 8009, Beijing 1000882Center for Applied Physics and Technology, Peking University, Beijing 100871
摘要We present an accurate and efficient generalized pseudospectral method for solving the time-dependent Schrödinger equation for atomic systems interacting with intense laser fields. In this method, the time propagation of the wave function is calculated using the well-known second-order split-operator method implemented by the numerically exact, fast transform between the grid and spectral representations. In the grid representation, the radial coordinate is discretized using the Coulomb wave discrete variable representation (CWDVR), and the angular dependence of the wave function is expanded in the Gauss-Legendre-Fourier grid. In the spectral representation, the wave function is expanded in terms of the eigenfunctions of the field-free zero-order Hamiltonian. Calculations on the high order harmonic generation and ionization dynamics of hydrogen atom in strong laser pulses are presented to demonstrate the accuracy and efficiency of the present method. This new algorithm will be found more computationally attractive than the close-coupled wave packet method using CWDVR and/or methods based on evenly spaced grids.
Abstract：We present an accurate and efficient generalized pseudospectral method for solving the time-dependent Schrödinger equation for atomic systems interacting with intense laser fields. In this method, the time propagation of the wave function is calculated using the well-known second-order split-operator method implemented by the numerically exact, fast transform between the grid and spectral representations. In the grid representation, the radial coordinate is discretized using the Coulomb wave discrete variable representation (CWDVR), and the angular dependence of the wave function is expanded in the Gauss-Legendre-Fourier grid. In the spectral representation, the wave function is expanded in terms of the eigenfunctions of the field-free zero-order Hamiltonian. Calculations on the high order harmonic generation and ionization dynamics of hydrogen atom in strong laser pulses are presented to demonstrate the accuracy and efficiency of the present method. This new algorithm will be found more computationally attractive than the close-coupled wave packet method using CWDVR and/or methods based on evenly spaced grids.
[1] Tong X M and Chu S I 2000 Phys. Rev. A 61031401 [2] Colgan J, Pindzola M S and Robicheaux F 2007 Phys.Rev. Lett. 98 153001 [3] Balint-Kurti G G 2003 Adv. Chem. Phys. 128 249 [4] Althorpe S C 2001 J. Chem. Phys. 114 1601 [5] Brabec T and Krausz F 2000 Rev. Mod. Phys. 72545 [6] Balint-Kurti G G, Zou S and Brown A 2008 Adv. Chem.Phys. 138 43 [7] Zou S, Balint-Kurti G G and Manby F R 2007 J. Chem.Phys. 127 044107 [8] Zou S, Ren Q, Balint-Kurti G G and Manby F R 2006 Phys. Rev. Lett. 96 243003 [9] Zou S, Sanz C and Balint-Kurti G G 2008 J. Chem.Phys. 129 124307 [10] Abu-samha M, Dimitrovski D and Madsen L B 2008 J.Phys. B 41 245601 [11] Steeg G L V, Bartschat K and Bray I 2003 J. Phys. B 36 3325 [12] Corey G C and Lemoine D 1992 J. Chem. Phys. 97 4115 [13] Corey G C and Tromp J W 1995 J. Chem. Phys. 103 1812 [14] Dunseath K M et al 2002 J. Phys. B 35 3539 [15] Peng L Y et al 2006 J. Chem. Phys. 125 154311 [16] Feit M D and Fleck J A 1983 J. Chem. Phys. 78301 [17] Feit M D, Fleck J A and Steiger A 1982 J. Comput.Phys. 47 412 [18] Gonzalez-Lezana T et al 2004 J. Chem. Phys. 120 2247 [19] Manolopoulos D E 2002 J. Chem. Phys. 117 9552 [20] Chu S I and Cooper J 1985 Phys. Rev. A 322769 [21] Pont M, Proulx D and Shakeshaft R 1991 Phys. Rev. A 44 4486 [22] Kulander K C 1987 Phys. Rev. A 35 445
2009, Vol. 26 
