New optimized Schwarz algorithms for one dimensional Schrödinger equation with general potential

The aim of this paper is to develop new optimized Schwarz algorithms for the one dimensional Schrodinger equation with linear or nonlinear potential. After presenting the classical algorithm which is an iterative process, we propose a new algorithm for the Schrodinger equation with time-independent linear potential. Thanks to two main ingredients (constructing explicitly the interface problem and using a direct method on the interface problem), the new algorithm turns to be a direct process. Thus, it is free to choose the transmission condition. Concerning the case of time-dependent linear potential or nonlinear potential, we propose to use a pre-processed linear operator as preconditioner which leads to a preconditioned algorithm. Numerically , the convergence is also independent of the transmission condition. In addition, both of these new algorithms implemented in parallel cluster are robust, scalable up to 256 sub domains (MPI process) and take much less computation time than the classical one, especially for the nonlinear case.

[1]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[2]  Christophe Geuzaine,et al.  A quasi-optimal non-overlapping domain decomposition algorithm for the Helmholtz equation , 2012, J. Comput. Phys..

[3]  Christophe Besse,et al.  Schwarz waveform relaxation method for one-dimensional Schrödinger equation with general potential , 2016, Numerical Algorithms.

[4]  Christophe Besse,et al.  Absorbing boundary conditions for the one-dimensional Schrödinger equation with an exterior repulsive potential , 2009, J. Comput. Phys..

[5]  Martin J. Gander,et al.  Optimized Schwarz Methods without Overlap for the Helmholtz Equation , 2002, SIAM J. Sci. Comput..

[6]  Christophe Besse,et al.  Absorbing Boundary Conditions for General Nonlinear Schrödinger Equations , 2011, SIAM J. Sci. Comput..

[7]  A. Durán,et al.  The numerical integration of relative equilibrium solutions. The nonlinear Schrödinger equation , 2000 .

[8]  Frédéric Nataf,et al.  FACTORIZATION OF THE CONVECTION-DIFFUSION OPERATOR AND THE SCHWARZ ALGORITHM , 1995 .

[9]  Martin J. Gander,et al.  Optimized Schwarz Methods , 2006, SIAM J. Numer. Anal..

[10]  Feng Xing,et al.  Domain Decomposition Algorithms for Two Dimensional Linear Schrödinger Equation , 2015, Journal of Scientific Computing.

[11]  Xavier Antoine,et al.  Domain decomposition methods and high-order absorbing boundary conditions for the numerical simulation of the time dependent Schrödinger equation with ionization and recombination by intense electric field , 2014 .

[12]  Laurence Halpern,et al.  Optimized and Quasi-optimal Schwarz Waveform Relaxation for the One Dimensional Schrödinger Equation , 2010 .

[13]  Martin J. Gander,et al.  Schwarz Methods over the Course of Time , 2008 .

[14]  F. Xing,et al.  Domain decomposition algorithms for the two dimensional nonlinear Schr{\"o}dinger equation and simulation of Bose-Einstein condensates , 2016, 1603.05110.