Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2002.04400 · European Journal of Physics

A numerical damped oscillator approach to constrained Schrödinger equations

Magnus Ögren, M Gulliksson

Published Jul 17, 2020• 5 citations
arXivHTMLPDF

Abstract

Abstract This article explains and illustrates the use of a set of coupled dynamical equations, second order in a fictitious time, which converges to solutions of stationary Schrödinger equations with additional constraints. In fact, the method is general and can solve constrained minimization problems in many fields. We present the method for introductory applications in quantum mechanics including three qualitative different numerical examples: the radial Schrödinger equation for the hydrogen atom; the 2D harmonic oscillator with degenerate excited states; and a nonlinear Schrödinger equation for rotating states. The presented method is intuitive, with analogies in classical mechanics for damped oscillators, and easy to implement, either with coding or with software for dynamical systems. Hence, we find it suitable to introduce it in a continuation course in quantum mechanics or generally in applied mathematics courses which contain computational parts. The undergraduate student can, for example, use our derived results and the code (supplemental material ( https://stacks.iop.org/EJP/41/065406/mmedia )) to study the Schrödinger equation in 1D for any potential. The graduate student and the general physicist can work from our three examples to derive their own results for other models including other global constraints.

DOI
10.1088/1361-6404/aba70b

Related papers

A numerical damped oscillator approach to constrained Schrödinger equations

Magnus Ögren, Gulliksson, Mårten

This article explains and illustrates the use of a set of coupled dynamical equations, second order in a fictitious time, which converges to solutions of stationary Schrödinger equations with additional constraints. In fact, the method is general and can solve constrained minimization problems in many fields. We present the method for introductory applications in quantum mechanics including three qualitative different numerical examples: the radial Schrödinger equation for the hydrogen atom; the 2D harmonic oscillator with degenerate excited states; and a nonlinear Schrödinger equation for rotating states. The presented method is intuitive, with analogies in classical mechanics for damped oscillators, and easy to implement, either with coding or with software for dynamical systems. Hence, we find it suitable to introduce it in a continuation course in quantum mechanics or generally in applied mathematics courses which contain computational parts. The undergraduate student can, for example, use our derived results and the code (supplemental material (https://stacks.iop.org/EJP/41/065406/mmedia)) to study the Schrödinger equation in 1D for any potential. The graduate student and the general physicist can work from our three examples to derive their own results for other models including other global constraints.

•Jan 1, 2020•2 citations

Numerical solution to the Schrödinger equation by constrained optimization

James Garner

The nonlinear optimization techniques of dynamic simulated annealing and steepest descents make it feasible to solve numerically the Schrödinger equation for systems that have many degrees of freedom. These methods involve finding the solution to a set of quasi-Newtonian equations of motion derived from a fictitious, time-dependent Schrödinger equation. The formalism is reviewed and the methods are illustrated by solving the Schrödinger equation for three simple models: a two-state system, an infinite square well, and the Mathieu equation.

Computers in Physics•Jul 1, 1990•0 citations

Unitary approach to the quantum forced harmonic oscillator

D. Velasco-Martínez, Victor Guadalupe Ibarra‐Sierra, Juan Carlos Sandoval-Santana +2 more

In this paper we introduce an alternative approach to studying the evolution of a quantum harmonic oscillator subject to an arbitrary time dependent force. With the purpose of finding the evolution operator, certain unitary transformations are applied successively to Schrödinger's equation reducing it to its simplest form. Therefore, instead of solving the original Schrödinger's partial differential equation in time and space the problem is replaced by a system of ordinary differential equations. From the obtained evolution operator we workout the propagator. Even though we illustrate the use of unitary transformations on the solution of a forced harmonic oscillator, the method presented here might be used to solve more complex systems. The present work addresses many aspects regarding unitary transformations and the dynamics of a forced quantum harmonic oscillator that should be useful for students and tutors of the quantum mechanics courses at the senior undergraduate and graduate level.

arXiv:1409.0236•Aug 31, 2014•1 citation

New Developments in Time-Independent Quantum-Mechanical Perturbation Theory

Gerald I. Kerley

This report discusses two new ideas for using perturbation methods to solve the time-independent Schrödinger equation. The first concept begins with rewriting the perturbation equations in a form that is closely related to matrix diagonalization methods. That allows a simple, compact derivation of the standard (Rayleigh-Schrödinger) equations. But it also leads to a new iterative solution method that is not based on the usual power series expansion. The iterative method can also be used when the unperturbed system is two-fold degenerate. The second concept is quite different from the first but compatible with it. It is based on the fact that one can replace the true Hamiltonian with a synthetic Hamiltonian having the same eigenvalues. This approach allows one to cancel part of the perturbation and to reduce the size of the off-diagonal matrix elements, giving better convergence of the series. These methods are illustrated by application to three perturbed harmonic oscillator problems--a 1-D oscillator with a linear perturbation, a 1-D oscillator with a quartic perturbation, and a perturbed 2-D oscillator, which involves degenerate states. The iterative method and the synthetic Hamiltonian method are shown to give significantly better results than the standard method.

arXiv:1306.6912•Jun 24, 2013•0 citations

Derivation and application of quantum Hamilton equations of motion

Jeanette Köppe, Wilfried Grecksch, Wolfgang Paul

Non‐relativistic quantum systems are analyzed theoretically or by numerical approaches using the Schrödinger equation. Compared to the options available to treat classical mechanical systems this is limited, both in methods and in scope. However, based on Nelson's stochastic mechanics, the mathematical structure of quantum mechanics has in some aspects been developed into a form analogous to classical analytical mechanics. We show here that finding the Nash equilibrium for a stochastic optimal control problem, which is the quantum equivalent to Hamilton's principle of least action, allows to derive two things: i) the Schrödinger equation as the Hamilton‐Jacobi‐Bellman equation of this optimal control problem and ii) a set of quantum dynamical equations which are the generalization of Hamilton's equations of motion to the quantum world. We derive their general form for the non‐stationary and the stationary case. For the harmonic oscillator, the stationary equations lead to the coherent states, and we establish a numerical procedure to solve for the ground state properties without using the Schrödinger equation. image

Annalen der Physik•Dec 13, 2016•14 citations

Quantum Harmonic Oscillators in One and Two Dimensions

Guanghan Xia

The Schrödinger equation is a significant achievement on the development of quantum mechanics. By solving the Schrödinger equation, the fundamental behaviors and properties of a microscopic particle can be found in one- to three-dimensions. The article focuses on the derivations of quantum harmonic oscillators in one dimension by solving second order-differential equation. The wave functions, probability densities under different energy levels are presented. The results can be used to estimate different forms of continuous potential experienced by an oscillator. By calculating the uncertainty relation of the oscillator under one specific excited state, the general relation can be confirmed. Based on the output, by splitting the variables, the two-dimensional harmonic oscillator can also be derived. The degenerate energy levels are presented as . Numerical simulations are made to visualize those results and suggests the relation of energy levels and number of maxima of the probability densities.

Highlights in Science Engineering and Technology•Aug 21, 2023•0 citations
Abstract This article explains and illustrates the use of a set of coupled dynamical equations, second order in a fictitious time, which converges to solutions of stationary Schrödinger equations with additional constraints. In fact, the method is general and can solve constrained minimization problems in many fields. We present the method for introductory applications in quantum mechanics including three qualitative different numerical examples: the radial Schrödinger equation for the hydrogen atom; the 2D harmonic oscillator with degenerate excited states; and a nonlinear Schrödinger equation for rotating states. The presented method is intuitive, with analogies in classical mechanics for damped oscillators, and easy to implement, either with coding or with software for dynamical systems. Hence, we find it suitable to introduce it in a continuation course in quantum mechanics or generally in applied mathematics courses which contain computational parts. The undergraduate student can, for example, use our derived results and the code (supplemental material ( https://stacks.iop.org/EJP/41/065406/mmedia )) to study the Schrödinger equation in 1D for any potential. The graduate student and the general physicist can work from our three examples to derive their own results for other models including other global constraints.
This article explains and illustrates the use of a set of coupled dynamical equations, second order in a fictitious time, which converges to solutions of stationary Schrödinger equations with additional constraints. In fact, the method is general and can solve constrained minimization problems in many fields. We present the method for introductory applications in quantum mechanics including three qualitative different numerical examples: the radial Schrödinger equation for the hydrogen atom; the 2D harmonic oscillator with degenerate excited states; and a nonlinear Schrödinger equation for rotating states. The presented method is intuitive, with analogies in classical mechanics for damped oscillators, and easy to implement, either with coding or with software for dynamical systems. Hence, we find it suitable to introduce it in a continuation course in quantum mechanics or generally in applied mathematics courses which contain computational parts. The undergraduate student can, for example, use our derived results and the code (supplemental material (https://stacks.iop.org/EJP/41/065406/mmedia)) to study the Schrödinger equation in 1D for any potential. The graduate student and the general physicist can work from our three examples to derive their own results for other models including other global constraints.
The nonlinear optimization techniques of dynamic simulated annealing and steepest descents make it feasible to solve numerically the Schrödinger equation for systems that have many degrees of freedom. These methods involve finding the solution to a set of quasi-Newtonian equations of motion derived from a fictitious, time-dependent Schrödinger equation. The formalism is reviewed and the methods are illustrated by solving the Schrödinger equation for three simple models: a two-state system, an infinite square well, and the Mathieu equation.
In this paper we introduce an alternative approach to studying the evolution of a quantum harmonic oscillator subject to an arbitrary time dependent force. With the purpose of finding the evolution operator, certain unitary transformations are applied successively to Schrödinger's equation reducing it to its simplest form. Therefore, instead of solving the original Schrödinger's partial differential equation in time and space the problem is replaced by a system of ordinary differential equations. From the obtained evolution operator we workout the propagator. Even though we illustrate the use of unitary transformations on the solution of a forced harmonic oscillator, the method presented here might be used to solve more complex systems. The present work addresses many aspects regarding unitary transformations and the dynamics of a forced quantum harmonic oscillator that should be useful for students and tutors of the quantum mechanics courses at the senior undergraduate and graduate level.
This report discusses two new ideas for using perturbation methods to solve the time-independent Schrödinger equation. The first concept begins with rewriting the perturbation equations in a form that is closely related to matrix diagonalization methods. That allows a simple, compact derivation of the standard (Rayleigh-Schrödinger) equations. But it also leads to a new iterative solution method that is not based on the usual power series expansion. The iterative method can also be used when the unperturbed system is two-fold degenerate. The second concept is quite different from the first but compatible with it. It is based on the fact that one can replace the true Hamiltonian with a synthetic Hamiltonian having the same eigenvalues. This approach allows one to cancel part of the perturbation and to reduce the size of the off-diagonal matrix elements, giving better convergence of the series. These methods are illustrated by application to three perturbed harmonic oscillator problems--a 1-D oscillator with a linear perturbation, a 1-D oscillator with a quartic perturbation, and a perturbed 2-D oscillator, which involves degenerate states. The iterative method and the synthetic Hamiltonian method are shown to give significantly better results than the standard method.
Non‐relativistic quantum systems are analyzed theoretically or by numerical approaches using the Schrödinger equation. Compared to the options available to treat classical mechanical systems this is limited, both in methods and in scope. However, based on Nelson's stochastic mechanics, the mathematical structure of quantum mechanics has in some aspects been developed into a form analogous to classical analytical mechanics. We show here that finding the Nash equilibrium for a stochastic optimal control problem, which is the quantum equivalent to Hamilton's principle of least action, allows to derive two things: i) the Schrödinger equation as the Hamilton‐Jacobi‐Bellman equation of this optimal control problem and ii) a set of quantum dynamical equations which are the generalization of Hamilton's equations of motion to the quantum world. We derive their general form for the non‐stationary and the stationary case. For the harmonic oscillator, the stationary equations lead to the coherent states, and we establish a numerical procedure to solve for the ground state properties without using the Schrödinger equation. image
The Schrödinger equation is a significant achievement on the development of quantum mechanics. By solving the Schrödinger equation, the fundamental behaviors and properties of a microscopic particle can be found in one- to three-dimensions. The article focuses on the derivations of quantum harmonic oscillators in one dimension by solving second order-differential equation. The wave functions, probability densities under different energy levels are presented. The results can be used to estimate different forms of continuous potential experienced by an oscillator. By calculating the uncertainty relation of the oscillator under one specific excited state, the general relation can be confirmed. Based on the output, by splitting the variables, the two-dimensional harmonic oscillator can also be derived. The degenerate energy levels are presented as . Numerical simulations are made to visualize those results and suggests the relation of energy levels and number of maxima of the probability densities.