Hostname: page-component-586b7cd67f-g8jcs Total loading time: 0 Render date: 2024-11-29T08:17:47.107Z Has data issue: false hasContentIssue false
  • English
  • Français

A Preconditioned Iterative Solver for the Scattering Solutions of the Schrödinger Equation

Published online by Cambridge University Press:  20 August 2015

Hisham bin Zubair ,
Bram Reps  and
Wim Vanroose
Hisham bin Zubair*
Affiliation:
Department of Mathematical Sciences, Faculty of Computer Science, Institute of Business Administration, University Rd., 75270 Karachi, Pakistan
Bram Reps*
Affiliation:
Department of Mathematics and Computer Science, University of Antwerp, Middelheimlaan 1, B-2020 Antwerpen, Belgium
Wim Vanroose*
Affiliation:
Department of Mathematics and Computer Science, University of Antwerp, Middelheimlaan 1, B-2020 Antwerpen, Belgium
*
Corresponding author.Email:[email protected]
Email address:[email protected]
Email address:[email protected]
Get access

Abstract

The Schrödinger equation defines the dynamics of quantum particles which has been an area of unabated interest in physics. We demonstrate how simple transformations of the Schrödinger equation leads to a coupled linear system, whereby each diagonal block is a high frequency Helmholtz problem. Based on this model, we derive indefinite Helmholtz model problems with strongly varying wavenumbers. We employ the iterative approach for their solution. In particular, we develop a preconditioner that has its spectrum restricted to a quadrant (of the complex plane) thereby making it easily invertible by multigrid methods with standard components. This multigrid preconditioner is used in conjunction with suitable Krylov-subspace methods for solving the indefinite Helmholtz model problems. The aim of this study is to report the feasibility of this preconditioner for the model problems. We compare this idea with the other prevalent preconditioning ideas, and discuss its merits. Results of numerical experiments are presented, which complement the proposed ideas, and show that this preconditioner may be used in an automatic setting.

Keywords

65F0865F1065F1565M55ScatteringSchrödinger equationexterior complex scalingpreconditioningmulti-gridcomplex-shifted Laplacian (CSL)complex-scaled grid (CSG)quadrant-definite(QD)
Type
Research Article
Information
Communications in Computational Physics , Volume 11 , Issue 2 , February 2012 , pp. 415 - 434
Copyright
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Vanroose, W., Martin, F., Rescigno, T.N., McCurdy, C. W., Complete photo-induced breakup of the H2 molecule as a probe of molecular electron correlation, Science 310 (2005) 1787–1789.CrossRefGoogle ScholarPubMed
[2]Taylor, K., Computational challenges in atomic, molecular and optical physics, Philosophical Transactions: Mathematical, Physical and Engineering Sciences (2002) 1135–1147.Google Scholar
[3]Simon, B., The definition of molecular resonance curves by the method of exterior complex scaling, Physics Letters A 71 (1979) 211.CrossRefGoogle Scholar
[4]Moiseyev, N., Quantum theory of resonances: Calculating energies, widths and cross-sections by complex scaling, Physics Reports 302 (5) (1998) 211.CrossRefGoogle Scholar
[5]Bérenger, J.-P., A perfectly matched layer for the absorption of electromagnetic waves, Journal of Computational Physics 114 (1994) 185–200.Google Scholar
[6]Chew, W. C., Weedon, W. H., A 3d perfectlymatched medium from modified Maxwell’s equations with stretched coordinates, Microwave and Optical Technology Letters 7 (13) (1994) 599–604.CrossRefGoogle Scholar
[7]Reps, B., Vanroose, W., H. bin Zubair, , On the indefinite Helmholtz equation: Exterior complex scaled absorbing boundary layers, iterative analysis, and preconditioning, Journal of Computational Physics 229 (2010) 8384–8405.Google Scholar
[8]Antoine, X., Arnold, A., Besse, C., Ehrhardt, M., Schaedle, A., A review of transparent and artificial boundary conditions techniques for linear and nonlinear schrödinger equations, Communications in Computational Physics 4 (4) (2008) 729–796.Google Scholar
[9]Bayliss, A., Goldstein, C. I., Turkel, E., An iterative method for the Helmholtz equation, Journal of Computational Physics 49 (1983) 443–457.CrossRefGoogle Scholar
[10]Bayliss, A., Goldstein, C. I., Turkel, E., On accuracy conditions for the numerical computation of waves, Journal of Computational Physics 59 (1985) 396–404.Google Scholar
[11]Erlangga, Y., Vuik, C., Oosterlee, C. W., On a class of preconditioners for solving the Helmholtz equation, Applied Numerical Mathematics 50 (2004) 409–425.CrossRefGoogle Scholar
[12]Erlangga, Y., Oosterlee, C. W., Vuik, C., A novel multigrid based preconditioner for heterogeneous Helmholtz problems, SIAM Journal on Scientific Computing 27 (2006) 1471–1492.Google Scholar
[13]Brandt, A., Multi-level adaptive solutions to boundary-value problems, Mathematics of Computation 31 (1977) 333–390.Google Scholar
[14]Stüben, K., Trottenberg, U., Multigrid Methods: Fundamental algorithms, model problem analysis and applications, Springer Berlin, 1982.Google Scholar
[15]Trottenberg, U., Oosterlee, C. W., Schöller, A., Multigrid, Academic Press, 2001.Google Scholar
[16]Elman, O. C. E. H. C., O’Leary, D.P., A multigrid method enhanced by Krylov subspace iteration for discrete Helmholtz equations, SIAM Journal on Scientific Computing 23 (2001) 1291–1315.CrossRefGoogle Scholar
[17]Meerbergen, K., Coyette, J., Connection and comparison between frequency shift time integration and a spectral transformation preconditioner, Numerical Linear Algebra with Applications 16 (1).Google Scholar
[18]Rescigno, T., Baertschy, M., Isaacs, W., McCurdy, C., Collisional breakup in a quantum system of three charged particles, Science 286 (5449) (1999) 2474.Google Scholar
[19]Baertschy, M., Rescigno, T., Isaacs, W., Li, X., McCurdy, C., Electron-impact ionization of atomic hydrogen, Physical Review A 63 (2) (2001) 22712.Google Scholar
[20]Vanroose, W., Horner, D., Martín, F., Rescigno, T., McCurdy, C., Double photoionization of aligned molecular hydrogen, Physical Review A 74 (5) (2006) 52702.Google Scholar
[21]Arfken, G. B., Weber., H. J., Mathematical Methods for Physicists, Academic Press, 1995.Google Scholar
[22]Baertschy, M., Li, X., Solution of a three-body problem in quantum mechanics using sparse linear algebra on parallel computers, in: Proceedings of the 2001 ACM/IEEE conference on Supercomputing (CDROM), ACM, 2001, p. 47.Google Scholar
[23]McCurdy, C., Horner, D., Rescigno, T., Martin, F., Theoretical treatment of double photoioniza-tion of helium using a B-spline implementation of exterior complex scaling, Physical Review A 69 (3) (2004) 32707.Google Scholar
[24]Horner, D., Morales, F., Rescigno, T., Martín, F., McCurdy, C., Two-photon double ionization of helium above and below the threshold for sequential ionization, Physical Review A 76 (3) (2007) 30701.Google Scholar
[25]Rappaport, C. M., Perfectly matched absorbing boundary conditions based on anisotropic lossy mapping of space, IEEE Microwave and Guided Wave Letters 5 (3) (1995) 90–92.Google Scholar
[26]Teixeira, F. L., Chew, W. C., General closed-form PML constitutive tensors to match arbitrary bianisotropic and dispersive linear media, IEEE Microwave and Guided Wave Letters 8 (6) (1998) 223–225.CrossRefGoogle Scholar
[27]Nicolaides, C. A., Beck, D. R., The variational calculation of energies and widths of resonances, Physics Letters A 65 (1978) 11.Google Scholar
[28]McCurdy, C. W., Baertschy, M., Rescigno, T. N., Solving the three-body Coulomb breakup problem using exterior complex scaling, Journal of Physics B: Atomic, Molecular and Optical Physics 37 (17) (2004) 137–187.CrossRefGoogle Scholar