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.
Mots clés : the time-dependent Schrödinger equation, the Crank-Nicolson finite-difference scheme, the strang splitting, approximate and discrete transparent boundary conditions, stability, tunnel effect
@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 }, pages = {1681--1699}, publisher = {EDP-Sciences}, volume = {48}, number = {6}, year = {2014}, doi = {10.1051/m2an/2014004}, mrnumber = {3264369}, language = {en}, url = {http://www.numdam.org/articles/10.1051/m2an/2014004/} }
TY - JOUR AU - Ducomet, Bernard AU - Zlotnik, Alexander AU - Zlotnik, Ilya TI - The splitting in potential Crank-Nicolson scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2014 SP - 1681 EP - 1699 VL - 48 IS - 6 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an/2014004/ DO - 10.1051/m2an/2014004 LA - en ID - M2AN_2014__48_6_1681_0 ER -
%0 Journal Article %A Ducomet, Bernard %A Zlotnik, Alexander %A Zlotnik, Ilya %T The splitting in potential Crank-Nicolson scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2014 %P 1681-1699 %V 48 %N 6 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an/2014004/ %R 10.1051/m2an/2014004 %G en %F 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 , Tome 48 (2014) no. 6, pp. 1681-1699. doi : 10.1051/m2an/2014004. http://www.numdam.org/articles/10.1051/m2an/2014004/
[1] A review of transparent and artificial boundary conditions techniques for linear and nonlinear Schrödinger equations. Commun. Comp. Phys. 4 (2008) 729-796. | MR
, , , and ,[2] Numerical schemes for the simulation of the two-dimensional Schrödinger equation using non-reflecting boundary conditions. Math. Comp. 73 (2004) 1779-1799. | MR | Zbl
, and ,[3] Discrete transparent boundary conditions for the Schrödinger equation: fast calculations, approximation and stability. Commun. Math. Sci. 1 (2003) 501-556. | MR | Zbl
, and ,[4] Microscopic analysis of collective dynamics in low energy fission. Nuclear Physics A 428 (1984) 23-36.
, and ,[5] Time-dependent quantum collective dynamics applied to nuclear fission. Comp. Phys. Commun. 63 (1991) 365-374. | Zbl
, and ,[6] Splitting methods for the time-dependent Schrödinger equation. Phys. Lett. A 265 (2000) 35-42. | MR | Zbl
and ,[7] Limits on the lifetime of the shape isomer of 238U. Phys. Rev. C 45 (1984) 1700-1708.
, , and ,[8] Transparent and absorbing boundary conditions for the Schrödinger equation in a bounded domain. Numer. Funct. Anal. Optimiz. 18 (1997) 759-775. | MR | Zbl
,[9] 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 | Zbl
and ,[10] 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 | Zbl
and ,[11] 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
, and ,[12] 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 | Zbl
, and ,[13] Discrete transparent boundary conditions for the Schrödinger equation. Riv. Mat. Univ. Parma 6 (2001) 57-108. | MR | Zbl
and ,[14] Fourth-order alternating direction implicit compact finite difference schemes for two-dimensional Schrödinger equations. Appl. Numer. Math. 61 (2011) 593-614. | MR | Zbl
and ,[15] Convergence of a split-step Hermite method for Gross-Pitaevskii equation. IMA J. Numer. Anal. 31 (2011) 396-415. | MR | Zbl
,[16] Microscopic approach of fission dynamics applied to fragment kinetic energy and mass distribution in 238U. Phys. Rev. C 71 (2005) 4316.
, , and ,[17] Quantummechanical treatment of the penetration through a two-dimensional fission barrier. Nuclear Physics A 224 (1974) 116-139.
,[18] On splitting methods for Schrödinger-Poisson and cubic nonlinear Schrödinger equations. Math. Comput. 77 (2008) 2141-2153. | MR | Zbl
,[19] From quantum to classical molecular dynamics. Reduced models and numerical analysis. Zürich Lect. Adv. Math. EMS, Zürich (2008). | MR | Zbl
,[20] On the convergence of splitting methods for linear evolutionary Schrödinger equations involving an unbounded potential. BIT Numer. Math. 49 (2009) 199-215. | MR | Zbl
and ,[21] On the treatment of a two-dimensional fission model with complex trajectories. Nuclear Physics A 296 (1978) 50-76.
, and ,[22] The nuclear many-body problem. Theoret. Math. Phys. Springer-Verlag, New York, Heidelberg, Berlin (1980). | MR
and ,[23] Microscopic dynamic calculations of collective states in Xenon and Barium isotopes. Nuclear Physics A 292 (1978) 66-87.
, , , and ,[24] Non-reflecting boundary conditions for the two-dimensional Schrödinger equation. Wave Motion 35 (2002) 181-188. | MR | Zbl
,[25] Discrete transparent boundary conditions for the Schrödinger equation, a compact higher order scheme. Kinetic Relat. Models 1 (2008) 101-125. | MR | Zbl
and ,[26] On the construction and comparison of difference scheme. SIAM J. Numer. Anal. 5 (1968) 506-517. | MR | Zbl
,[27] Design of absorbing boundary conditions for Schrödinger equations in Rd. SIAM J. Numer. Anal. 42 (2004) 1527-1551. | MR | Zbl
,[28] The method of fractional steps: solution of problems of mathematical physics in several variables. Springer, New York (1971). | MR | Zbl
,[29] 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 | Zbl
and ,[30] Sharp error analysis of vector splitting methods for the heat equation. Comput. Math. Phys. 39 (1999) 448-467. | MR | Zbl
and ,[31] 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 | Zbl
,[32] Sharp error bounds for a symmetrized locally 1D method for solving the 2D heat equation. Comput. Meth. Appl. Math. 6 (2006) 94-114. | MR | Zbl
and ,[33] 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
and ,[34] 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 | Zbl
and ,[35] Computer simulation of the tunnel effect. Moscow Power Engin. Inst. Bulletin 17 (2010) 10-28 (in Russian).
,[36] 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 | Zbl
,Cité par Sources :