An efficient method is proposed for numerical solutions of nonlinear Schrödinger equations on an unbounded domain. Through approximating the kinetic energy term by a one-way equation and uniting it with the potential energy equation, absorbing boundary conditions are designed to truncate the unbounded domain, which are in nonlinear form and can perfectly absorb waves outgoing from the boundaries of the truncated computational domain. The stability of the induced initial boundary value problem defined on the computational domain is examined by a normal mode analysis. Numerical examples are given to illustrate the stable and tractable advantages of the method.

• Received 14 May 2008

DOI:https://doi.org/10.1103/PhysRevE.78.026709