The splitting in potential Crank-Nicolson scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip
Ducomet, Bernard ; Zlotnik, Alexander ; Zlotnik, Ilya
ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 48 (2014), p. 1681-1699 / Harvested from Numdam

We consider an initial-boundary value problem for a generalized 2D time-dependent Schrödinger equation (with variable coefficients) on a semi-infinite strip. For the Crank-Nicolson-type finite-difference scheme with approximate or discrete transparent boundary conditions (TBCs), the Strang-type splitting with respect to the potential is applied. For the resulting method, the unconditional uniform in time L2-stability is proved. Due to the splitting, an effective direct algorithm using FFT is developed now to implement the method with the discrete TBC for general potential. Numerical results on the tunnel effect for rectangular barriers are included together with the detailed practical error analysis confirming nice properties of the method.

Publié le : 2014-01-01
DOI : https://doi.org/10.1051/m2an/2014004
Classification:  65M06,  65M12,  35Q40
@article{M2AN_2014__48_6_1681_0,
     author = {Ducomet, Bernard and Zlotnik, Alexander and Zlotnik, Ilya},
     title = {The splitting in potential Crank-Nicolson scheme with discrete transparent boundary conditions for the Schr\"odinger equation on a semi-infinite strip},
     journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
     volume = {48},
     year = {2014},
     pages = {1681-1699},
     doi = {10.1051/m2an/2014004},
     mrnumber = {3264369},
     language = {en},
     url = {http://dml.mathdoc.fr/item/M2AN_2014__48_6_1681_0}
}
Ducomet, Bernard; Zlotnik, Alexander; Zlotnik, Ilya. The splitting in potential Crank-Nicolson scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 48 (2014) pp. 1681-1699. doi : 10.1051/m2an/2014004. http://gdmltest.u-ga.fr/item/M2AN_2014__48_6_1681_0/

[1] X. Antoine, A. Arnold, C. Besse, M. Ehrhardt and A. Schädle, A review of transparent and artificial boundary conditions techniques for linear and nonlinear Schrödinger equations. Commun. Comp. Phys. 4 (2008) 729-796. | MR 2463189

[2] X. Antoine, C. Besse and V. Mouysset, Numerical schemes for the simulation of the two-dimensional Schrödinger equation using non-reflecting boundary conditions. Math. Comp. 73 (2004) 1779-1799. | MR 2059736 | Zbl 1053.65072

[3] A. Arnold, M. Ehrhardt and I. Sofronov, Discrete transparent boundary conditions for the Schrödinger equation: fast calculations, approximation and stability. Commun. Math. Sci. 1 (2003) 501-556. | MR 2069942 | Zbl 1085.65513

[4] J.F. Berger, M. Girod and D. Gogny, Microscopic analysis of collective dynamics in low energy fission. Nuclear Physics A 428 (1984) 23-36.

[5] J.-F. Berger, M. Girod and D. Gogny, Time-dependent quantum collective dynamics applied to nuclear fission. Comp. Phys. Commun. 63 (1991) 365-374. | Zbl 0850.65332

[6] S. Blanes and P.C. Moan, Splitting methods for the time-dependent Schrödinger equation. Phys. Lett. A 265 (2000) 35-42. | MR 1739100 | Zbl 0941.81025

[7] C.R. Chinn, J.F. Berger, D. Gogny and M.S. Weiss, Limits on the lifetime of the shape isomer of 238U. Phys. Rev. C 45 (1984) 1700-1708.

[8] L. Di Menza, Transparent and absorbing boundary conditions for the Schrödinger equation in a bounded domain. Numer. Funct. Anal. Optimiz. 18 (1997) 759-775. | MR 1472153 | Zbl 0895.65041

[9] B. Ducomet and A. Zlotnik, On stability of the Crank-Nicolson scheme with approximate transparent boundary conditions for the Schrödinger equation. Part I. Commun. Math. Sci. 4 (2006) 741-766. | MR 2264818 | Zbl 1119.65085

[10] B. Ducomet and A. Zlotnik, On stability of the Crank-Nicolson scheme with approximate transparent boundary conditions for the Schrödinger equation. Part II. Commun. Math. Sci. 5 (2007) 267-298. | MR 2334843 | Zbl 1141.65069

[11] B. Ducomet, A. Zlotnik and A. Romanova, On a splitting higher order scheme with discrete transparent boundary conditions for the Schrödinger equation in a semi-infinite parallelepiped. Appl. Math. Comp. To appear (2014). | MR 3316597

[12] B. Ducomet, A. Zlotnik and I. Zlotnik, On a family of finite-difference schemes with discrete transparent boundary conditions for a generalized 1D Schrödinger equation. Kinetic Relat. Models 2 (2009), 151-179. | MR 2472153 | Zbl 1195.65106

[13] M. Ehrhardt and A. Arnold, Discrete transparent boundary conditions for the Schrödinger equation. Riv. Mat. Univ. Parma 6 (2001) 57-108. | MR 1881105 | Zbl 0993.65097

[14] Z. Gao and S. Xie, Fourth-order alternating direction implicit compact finite difference schemes for two-dimensional Schrödinger equations. Appl. Numer. Math. 61 (2011) 593-614. | MR 2754579 | Zbl 1221.65220

[15] L. Gauckler, Convergence of a split-step Hermite method for Gross-Pitaevskii equation. IMA J. Numer. Anal. 31 (2011) 396-415. | MR 2813177 | Zbl 1223.65079

[16] H. Goutte, J.-F. Berger, P. Casoly and D. Gogny, Microscopic approach of fission dynamics applied to fragment kinetic energy and mass distribution in 238U. Phys. Rev. C 71 (2005) 4316.

[17] H. Hofmann, Quantummechanical treatment of the penetration through a two-dimensional fission barrier. Nuclear Physics A 224 (1974) 116-139.

[18] C. Lubich, On splitting methods for Schrödinger-Poisson and cubic nonlinear Schrödinger equations. Math. Comput. 77 (2008) 2141-2153. | MR 2429878 | Zbl 1198.65186

[19] C. Lubich, From quantum to classical molecular dynamics. Reduced models and numerical analysis. Zürich Lect. Adv. Math. EMS, Zürich (2008). | MR 2474331 | Zbl 1160.81001

[20] C. Neuhauser and M. Thalhammer, On the convergence of splitting methods for linear evolutionary Schrödinger equations involving an unbounded potential. BIT Numer. Math. 49 (2009) 199-215. | MR 2486135 | Zbl 1162.65385

[21] P. Ring, H. Hassman and J.O. Rasmussen, On the treatment of a two-dimensional fission model with complex trajectories. Nuclear Physics A 296 (1978) 50-76.

[22] P. Ring and P. Schuck, The nuclear many-body problem. Theoret. Math. Phys. Springer-Verlag, New York, Heidelberg, Berlin (1980). | MR 611683

[23] S.G. Rohozinski, J. Dobaczewski, B. Nerlo-Pomorska, K. Pomorski and J. Srebny, Microscopic dynamic calculations of collective states in Xenon and Barium isotopes. Nuclear Physics A 292 (1978) 66-87.

[24] A. Schädle, Non-reflecting boundary conditions for the two-dimensional Schrödinger equation. Wave Motion 35 (2002) 181-188. | MR 1869342 | Zbl 1163.74435

[25] M. Schulte and A. Arnold, Discrete transparent boundary conditions for the Schrödinger equation, a compact higher order scheme. Kinetic Relat. Models 1 (2008) 101-125. | MR 2383718 | Zbl 1145.65058

[26] G. Strang, On the construction and comparison of difference scheme. SIAM J. Numer. Anal. 5 (1968) 506-517. | MR 235754 | Zbl 0184.38503

[27] J. Szeftel, Design of absorbing boundary conditions for Schrödinger equations in Rd. SIAM J. Numer. Anal. 42 (2004) 1527-1551. | MR 2114289 | Zbl 1094.35037

[28] N.N. Yanenko, The method of fractional steps: solution of problems of mathematical physics in several variables. Springer, New York (1971). | MR 307493 | Zbl 0209.47103

[29] S.B. Zaitseva and A.A. Zlotnik, Error analysis in L2(Q) for symmetrized locally one-dimensional methods for the heat equation. Russ. J. Numer. Anal. Math. Model. 13 (1998) 69-91. | MR 1611638 | Zbl 0943.65104

[30] S.B. Zaitseva and A.A. Zlotnik, Sharp error analysis of vector splitting methods for the heat equation. Comput. Math. Phys. 39 (1999) 448-467. | MR 1683579 | Zbl 0970.65099

[31] A.A. Zlotnik, Some finite-element and finite-difference methods for solving mathematical physics problems with non-smooth data in n-dimensional cube. Sov. J. Numer. Anal. Math. Modell. 6 (1991) 421-451. | MR 1137008 | Zbl 0816.65063

[32] A. Zlotnik and S. Ilyicheva, Sharp error bounds for a symmetrized locally 1D method for solving the 2D heat equation. Comput. Meth. Appl. Math. 6 (2006) 94-114. | MR 2220487 | Zbl 1094.65094

[33] A. Zlotnik and A. Romanova, On a Numerov-Crank-Nicolson-Strang scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip. Appl. Numer. Math. To appear (2014). | MR 3323462

[34] A.A. Zlotnik and I.A. Zlotnik, Family of finite-difference schemes with transparent boundary conditions for the nonstationary Schrödinger equation in a semi-infinite strip. Dokl. Math. 83 (2011) 12-18. | MR 2810155 | Zbl 1236.65116

[35] I.A. Zlotnik, Computer simulation of the tunnel effect. Moscow Power Engin. Inst. Bulletin 17 (2010) 10-28 (in Russian).

[36] I.A. Zlotnik, Family of finite-difference schemes with approximate transparent boundary conditions for the generalized nonstationary Schrödinger equation in a semi-infinite strip. Comput. Math. Phys. 51 (2011) 355-376. | MR 2839572 | Zbl 1249.35313