Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2003
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/math/0303158 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866914099491766272 |
|---|---|
| author | Bao, Weizhu Jaksch, Dieter |
| author_facet | Bao, Weizhu Jaksch, Dieter |
| contents | This paper introduces an extension of the time-splitting sine-spectral (TSSP) method for solving damped focusing nonlinear Schrödinger equations (NLS). The method is explicit, unconditionally stable and time transversal invariant. Moreover, it preserves the exact decay rate for the normalization of the wave function if linear damping terms are added to the NLS. Extensive numerical tests are presented for cubic focusing nonlinear Schrödinger equations in 2d with a linear, cubic or a quintic damping term. Our numerical results show that quintic or cubic damping always arrests blowup, while linear damping can arrest blowup only when the damping parameter $\dt$ is larger than a threshold value $\dt_{\rm th}$. We note that our method can also be applied to solve the 3d Gross-Pitaevskii equation with a quintic damping term to model the dynamics of a collapsing and exploding Bose-Einstein condensate (BEC). |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_math_0303158 |
| institution | arXiv |
| publishDate | 2003 |
| record_format | arxiv |
| spellingShingle | An explicit unconditionally stable numerical method for solving damped nonlinear Schrödinger equations with a focusing nonlinearity Bao, Weizhu Jaksch, Dieter Numerical Analysis This paper introduces an extension of the time-splitting sine-spectral (TSSP) method for solving damped focusing nonlinear Schrödinger equations (NLS). The method is explicit, unconditionally stable and time transversal invariant. Moreover, it preserves the exact decay rate for the normalization of the wave function if linear damping terms are added to the NLS. Extensive numerical tests are presented for cubic focusing nonlinear Schrödinger equations in 2d with a linear, cubic or a quintic damping term. Our numerical results show that quintic or cubic damping always arrests blowup, while linear damping can arrest blowup only when the damping parameter $\dt$ is larger than a threshold value $\dt_{\rm th}$. We note that our method can also be applied to solve the 3d Gross-Pitaevskii equation with a quintic damping term to model the dynamics of a collapsing and exploding Bose-Einstein condensate (BEC). |
| title | An explicit unconditionally stable numerical method for solving damped nonlinear Schrödinger equations with a focusing nonlinearity |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/math/0303158 |