A split-step numerical method for the time-dependent Dirac equation in 3-D axisymmetric geometry

Journal of Computational Physics

Volumn 272, Issue , 2014, Pages 559-587

  Fillion Gourdeau, François ^a   Lorin, Emmanuel ^a,b   Bandrauk, André D ^a,c  

^a Center de Recherches Mathematiques   (Canada)

^b CARLETON UNIVERSITY   (Canada)

^c UNIVERSITÉ DE SHERBROOKE   (Canada)

Author keywords

Alternate direction iteration; Cylindrical coordinates; Dirac equation; Splitting methods

Indexed keywords

INTEGRAL EQUATIONS; ITERATIVE METHODS; LINEAR EQUATIONS; NUMERICAL METHODS; POLYNOMIALS;

ALTERNATE DIRECTION ITERATION; AXISYMMETRIC; AXISYMMETRIC GEOMETRIES; CYLINDRICAL COORDINATES; DIRAC'S EQUATION; SPLIT-STEP; SPLITTING METHOD; SPLITTINGS; TIME DEPENDENT; TIME EVOLUTIONS;

WAVE FUNCTIONS;

EID: 84900544223     PISSN: 00219991     EISSN: 10902716     Source Type: Journal    
DOI: 10.1016/j.jcp.2014.03.068     Document Type: Article

References (67)

1. Salamin Y.I., Hu S., Hatsagortsyan K.Z., Keitel C.H. Relativistic high-power laser-matter interactions. Phys. Rep. 2006, 427(2-3):41-155. 10.1016/j.physrep.2006.01.002.
2. Katsnelson M.I., Novoselov K.S., Geim A.K. Chiral tunnelling and the Klein paradox in graphene. Nat. Phys. 2006, 2:620-625. 10.1038/nphys384.
3. Kutzelnigg W. Basis set expansion of the Dirac operator without variational collapse. Int. J. Quant. Chem. 1984, 25(1):107-129. 10.1002/qua.560250112.
4. Greiner W. Relativistic Quantum Mechanics, Wave Equations 1987, Springer.
5. Bagrov V., Gitman D. Exact Solutions of Relativistic Wave Equations. Mathematics and Its Applications: Soviet series 1990, Kluwer Academic Publishers.
6. Milosevic N., Krainov V.P., Brabec T. Semiclassical Dirac theory of tunnel ionization. Phys. Rev. Lett. 2002, 89:193001. 10.1103/PhysRevLett.89.193001.
7. Maquet A., Grobe R. Atoms in strong laser fields: challenges in relativistic quantum mechanics. J. Mod. Opt. 2002, 49(12):2001-2018. 10.1080/09500340210140740.
8. Stacey R. Eliminating lattice fermion doubling. Phys. Rev. D 1982, 26(2):468-472. 10.1103/PhysRevD.26.468.
9. Müller C., Grün N., Scheid W. Finite element formulation of the Dirac equation and the problem of fermion doubling. Phys. Lett. A 1998, 242(4-5):245-250. 10.1016/S0375-9601(98)00218-7.
10. Kogut J.B., Susskind L. Hamiltonian formulation of Wilson's lattice gauge theories. Phys. Rev. D 1975, 11:395. 10.1103/PhysRevD.11.395.
11. Braun J.W., Su Q., Grobe R. Numerical approach to solve the time-dependent Dirac equation. Phys. Rev. A 1999, 59(1):604-612. 10.1103/PhysRevA.59.604.
12. Mocken G.R., Keitel C.H. FFT-split-operator code for solving the Dirac equation in 2 + 1 dimensions. Comput. Phys. Commun. 2008, 178(11):868-882. 10.1016/j.cpc.2008.01.042.
13. Mocken G.R., Keitel C.H. Quantum dynamics of relativistic electrons. J. Comput. Phys. 2004, 199(2):558-588. 10.1016/j.jcp.2004.02.020.
14. Momberger K., Belkacem A., Sørensen A.H. Numerical treatment of the time-dependent Dirac equation in momentum space for atomic processes in relativistic heavy-ion collisions. Phys. Rev. A 1996, 53(3):1605-1622. 10.1103/PhysRevA.53.1605.
15. Bauke H., Keitel C.H. Accelerating the Fourier split operator method via graphics processing units. Comput. Phys. Commun. 2011, 182(12):2454-2463. 10.1016/j.cpc.2011.07.003.
16. Huang Z., Jin S., Markowich P.A., Sparber C., Zheng C. A time-splitting spectral scheme for the Maxwell-Dirac system. J. Comput. Phys. 2005, 208(2):761-789. 10.1016/j.jcp.2005.02.026.
17. Bao W., Li X.-G. An efficient and stable numerical method for the Maxwell-Dirac system. J. Comput. Phys. 2004, 199(2):663-687. 10.1016/j.jcp.2004.03.003.
18. Xu J., Shao S., Tang H. Numerical methods for nonlinear Dirac equation. J. Comput. Phys. 2013, 245(0):131-149. 10.1016/j.jcp.2013.03.031.
19. Ackad E., Horbatsch M. Calculation of electron-positron production in supercritical uranium-uranium collisions near the Coulomb barrier. Phys. Rev. A 2008, 78:062711. 10.1103/PhysRevA.78.062711.
20. McConnell S.R., Artemyev A.N., Mai M., Surzhykov A. Solution of the two-center time-dependent Dirac equation in spherical coordinates: application of the multipole expansion of the electron-nuclei interaction. Phys. Rev. A 2012, 86:052705. 10.1103/PhysRevA.86.052705.
21. Bottcher C., Strayer M.R. Numerical solution of the time-dependent Dirac equation with application to positron production in heavy-ion collisions. Phys. Rev. Lett. 1985, 54(7):669-672. 10.1103/PhysRevLett.54.669.
22. U91 +. J. Phys. B, At. Mol. Opt. Phys. 1983, 16(11):1967.
23. Salomonson S., Öster P. Relativistic all-order pair functions from a discretized single-particle Dirac Hamiltonian. Phys. Rev. A 1989, 40(10):5548-5558. 10.1103/PhysRevA.40.5548.
24. Selstø S., Lindroth E., Bengtsson J. Solution of the Dirac equation for hydrogenlike systems exposed to intense electromagnetic pulses. Phys. Rev. A 2009, 79(4):043418. 10.1103/PhysRevA.79.043418.
25. Gelis F., Kajantie K., Lappi T. Quark-antiquark production from classical fields in heavy-ion collisions: 1 + 1 dimensions. Phys. Rev. C 2005, 71(2):024904. 10.1103/PhysRevC.71.024904.
26. Hammer R., Pötz W. Staggered grid leap-frog scheme for the (2 + 1)D Dirac equation. Comput. Phys. Commun. 2014, 185(1):40-52. 10.1016/j.cpc.2013.08.013.
27. Kogut J.B. The lattice gauge theory approach to quantum chromodynamics. Rev. Mod. Phys. 1983, 55(3):775-836. 10.1103/RevModPhys.55.775.
28. Trefethen L.N. Group velocity in finite difference schemes. SIAM Rev. 1982, 24(2):113-136. 10.1137/1024038.
29. Lorin E., Bandrauk A. A simple and accurate mixed P0-Q1 solver for the Maxwell-Dirac equations. Nonlinear Anal., Real World Appl. 2011, 12(1):190-202. 10.1016/j.nonrwa.2010.06.007.
30. Fillion-Gourdeau F., Lorin E., Bandrauk A.D. Numerical solution of the time-dependent Dirac equation in coordinate space without fermion-doubling. Comput. Phys. Commun. 2012, 183(7):1403-1415. 10.1016/j.cpc.2012.02.012.
31. Succi S., Benzi R. Lattice Boltzmann equation for quantum mechanics. Physica D: Nonlinear Phenomena 1993, 69(3-4):327-332. 10.1016/0167-2789(93)90096-J.
32. Huang W., Sloan D.M. Pole condition for singular problems: the pseudospectral approximation. J. Comput. Phys. 1993, 107(2):254-261. 10.1006/jcph.1993.1141.
33. Mohseni K., Colonius T. Numerical treatment of polar coordinate singularities. J. Comput. Phys. 2000, 157(2):787-795. 10.1006/jcph.1999.6382.
34. Constantinescu G., Lele S. A highly accurate technique for the treatment of flow equations at the polar axis in cylindrical coordinates using series expansions. J. Comput. Phys. 2002, 183(1):165-186. 10.1006/jcph.2002.7187.
35. Bandrauk A., Lu H.-Z. A singularity-free method for the time-dependent Schrödinger equation for nonlinear molecules. Proceedings of the 17th Annual International Symposium on High Performance Computing Systems and Applications and the OSCAR Symposium 2003, 91-98. NRC Research Press. D. Sénéchal (Ed.).
36. Mickens R. Applications of Nonstandard Finite Difference Schemes 2000, World Scientific.
37. Feit M., Fleck J., Steiger A. Solution of the Schrödinger equation by a spectral method. J. Comput. Phys. 1982, 47(3):412-433. 10.1016/0021-9991(82)90091-2.
38. Itzykson C., Zuber J.B. Quantum Field Theory 1980, Mcgraw-Hill.
39. Schluter P., Wietschorke K.-H., Greiner W. The Dirac equation in orthogonal coordinate systems. I. The local representation. J. Phys. A, Math. Gen. 1983, 16(9):1999.
40. Kullie O., Kolb D., Rutkowski A. Two-spinor fully relativistic finite-element (fem) solution of the two-center Coulomb problem. Chem. Phys. Lett. 2004, 383(3-4):215-221. 10.1016/j.cplett.2003.11.010.
41. Lewis H.R., Bellan P.M. Physical constraints on the coefficients of Fourier expansions in cylindrical coordinates. J. Math. Phys. 1990, 31(11):2592-2596. 10.1063/1.529009.
42. Polyanin A. Handbook of Linear Partial Differential Equations for Engineers and Scientists 2001, Chapman & Hall/CRC.
43. Evans L. Partial Differential Equations. Grad. Stud. Math. 1997, American Mathematical Society.
44. Bandrauk A.D., Dehghanian E., Lu H. Complex integration steps in decomposition of quantum exponential evolution operators. Chem. Phys. Lett. 2006, 419(4-6):346-350. 10.1016/j.cplett.2005.12.006.
45. Bandrauk A.D., Shen H. Exponential split operator methods for solving coupled time-dependent Schrödinger equations. J. Chem. Phys. 1993, 99(2):1185-1193. 10.1063/1.465362.
46. Suzuki M. General theory of fractal path integrals with applications to many-body theories and statistical physics. J. Math. Phys. 1991, 32(2):400-407. 10.1063/1.529425.
47. Sornborger A.T., Stewart E.D. Higher-order methods for simulations on quantum computers. Phys. Rev. A 1999, 60:1956-1965. 10.1103/PhysRevA.60.1956.
48. Bauke H., Hetzheim H.G., Mocken G.R., Ruf M., Keitel C.H. Relativistic ionization characteristics of laser-driven hydrogenlike ions. Phys. Rev. A 2011, 83(6):063414. 10.1103/PhysRevA.83.063414.
49. Su Q., Smetanko B., Grobe R. Relativistic suppression of wave packet spreading. Opt. Express 1998, 2(7):277-281. 10.1364/OE.2.000277.
50. Krekora P., Su Q., Grobe R. Effects of relativity on the time-resolved tunneling of electron wave packets. Phys. Rev. A 2001, 63(3):032107. 10.1103/PhysRevA.63.032107.
51. Thaller B. Visualizing the kinematics of relativistic wave packets arxiv:quant-ph/0409079.
52. Demikhovskii V.Y., Maksimova G.M., Perov A.A., Frolova E.V. Space-time evolution of Dirac wave packets. Phys. Rev. A 2010, 82(5):052115. 10.1103/PhysRevA.82.052115.
53. Park S.T. Propagation of a relativistic electron wave packet in the Dirac equation. Phys. Rev. A 2012, 86:062105. 10.1103/PhysRevA.86.062105.
54. Ruf M., Mocken G.R., Müller C., Hatsagortsyan K.Z., Keitel C.H. Pair production in laser fields oscillating in space and time. Phys. Rev. Lett. 2009, 102:080402. 10.1103/PhysRevLett.102.080402.
55. Esteban M. A short review on computational issues arising in relativistic atomic and molecular physics. CRM Proceedings and Lecture Notes 2007, vol. 41:105-116. American Mathematical Society. A. Bandrauk, M. Delfour, C.L. Bris (Eds.).
56. Dolbeault J., Esteban M.J., Séré E. A variational method for relativistic computations in atomic and molecular physics. Int. J. Quant. Chem. 2003, 93:149-155. 10.1002/qua.10549.
57. Desclaux J., Dolbeault J., Esteban M., Indelicato P., Séré E. Computational approaches of relativistic models in quantum chemistry. Handbook of Numerical Analysis 2003, vol. 10:453-483. Elsevier. 10.1016/S1570-8659(03)10006-3. C.L. Bris (Ed.).
58. Fillion-Gourdeau F., Lorin E., Bandrauk A.D. Numerical solution of the time-independent Dirac equation for diatomic molecules: b splines without spurious states. Phys. Rev. A 2012, 85:022506. 10.1103/PhysRevA.85.022506.
59. Shabaev V.M., Tupitsyn I.I., Yerokhin V.A., Plunien G., Soff G. Dual kinetic balance approach to basis-set expansions for the Dirac equation. Phys. Rev. Lett. 2004, 93:130405. 10.1103/PhysRevLett.93.130405.
60. Grant I.P., Quiney H.M. Rayleigh-Ritz approximation of the Dirac operator in atomic and molecular physics. Phys. Rev. A 2000, 62:022508. 10.1103/PhysRevA.62.022508.
61. Almanasreh H., Salomonson S., Svanstedt N. Stabilized finite element method for the radial Dirac equation. J. Comput. Phys. 2013, 236(0):426-442. 10.1016/j.jcp.2012.11.020.
62. Fischer C.F., Zatsarinny O. A b-spline Galerkin method for the Dirac equation. Comput. Phys. Commun. 2009, 180(6):879-886. 10.1016/j.cpc.2008.12.010.
63. Marsman A., Horbatsch M. Calculation of supercritical Dirac resonance parameters for heavy-ion systems from a coupled-differential-equation approach. Phys. Rev. A 2011, 84:032517. 10.1103/PhysRevA.84.032517.
64. Domínguez-Adame F., Rodríguez A. A one-dimensional relativistic screened Coulomb potential. Phys. Lett. A 1995, 198(4):275-278. 10.1016/0375-9601(95)00064-A.
65. Cole J.B. High-accuracy Yee algorithm based on nonstandard finite differences: new developments and verifications. IEEE Trans. Antennas Propag. 2002, 50(9):1185-1191.
66. Suzuki M. General decomposition theory of ordered exponentials. Proc. Jpn. Acad. Ser. B, Phys. Biol. Sci. 1993, 69(7):161-166.
67. Watson G. A Treatise on the Theory of Bessel Functions 1962, University Press, Cambridge.