摘要
In this paper we study numerical issues related to the Schr ¨odinger equationwith sinusoidal potentials at infinity. An exact absorbing boundary condition in a formof Dirichlet-to-Neumann mapping is derived. This boundary condition is based on ananalytical expression of the logarithmic derivative of the Floquet solution toMathieu’sequation, which is completely new to the author’s knowledge. The implementationof this exact boundary condition is discussed, and a fast evaluation method is used toreduce the computation burden arising from the involved half-order derivative operator.Some numerical tests are given to showthe performance of the proposed absorbingboundary conditions.
基金
the National Natural Science Foundation of China underGrant No. 10401020.