Computing solutions of Schrödinger equations on unbounded domains- On the brink of numerical algorithms

Abstract

We address the open problem of determining which classes of time-dependent linear Schrödinger equations and focusing and defocusing cubic and quintic non-linear Schrödinger equations (NLS) on unbounded domains that can be computed by an algorithm. We demonstrate how such an algorithm in general does not exist, yielding a substantial classification theory of which problems in quantum mechanics that can be computed. Moreover, we establish classifications on which problems that can be computed with a uniform bound on the runtime, as a function of the desired $\epsilon$-accuracy of the approximation. This include linear and nonlinear Schrödinger equations for which we provide positive and negative results and conditions on both the initial state and the potentials such that there exist computational (recursive) a priori bounds that allow reduction of the IVP on an unbounded domain to an IVP on a bounded domain, yielding an algorithm that can produce an $\epsilon$-approximation. In addition, we show how no algorithm can decide, and in fact not verify nor falsify, if the focusing NLS will blow up in finite time or not, yet, for the defocusing NLS, solutions can be computed given mild assumptions on the initial state and the potentials. Finally, we show that solutions to discrete NLS equations (focusing and defocusing) on an unbounded domain can always be computed with uniform bounds on the runtime of the algorithm. The algorithms presented are not just of theoretical interest, but efficient and easy to implement in applications. Our results have implications beyond computational quantum mechanics and are a part of the Solvability Complexity Index (SCI) hierarchy and Smale's program on the foundations of computational mathematics. For example our results provide classifications of which mathematical problems may be solved by computer assisted proofs.

arXiv commentComments welcome

References in corpus (8)

• Numerical Solution of the Gross-Pitaevskii Equation for Bose-Einstein Condensation
• A sharp condition for scattering of the radial 3d cubic nonlinear Schroedinger equation
• Controllability of the discrete-spectrum Schrodinger equation driven by an external field
• Computing quantum dynamics in the semiclassical regime
• Global well - posedness and scattering for the focusing, energy - critical nonlinear Schrödinger problem in dimension $d=4$ for initial data below a ground state threshold
• Computing Spectra -- On the Solvability Complexity Index Hierarchy and Towers of Algorithms
• Convergence rates for the quantum central limit theorem
• Compact schemes for laser-matter interaction in Schrödinger equation