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:2604.14812

Time-Dependent Logarithmic Perturbation Theory for Quantum Dynamics: Formulation and Applications

Juan Carlos del Valle, Paul Bergold, Karolina Kropielnicka

quant-phmath-phmath.APSubmitted Apr 16, 2026• Announced Apr 17, 2026
v1 · Apr 16, 2026v2 · Sep 3, 2026
arXivHTMLPDF

Abstract

We present a time-dependent extension of logarithmic perturbation theory for nonrelativistic quantum dynamics governed by the Schrödinger equation, in which the logarithm of the wave function is expanded in powers of a coupling constant. The resulting hierarchy of equations defining the perturbative corrections is governed by a gauge-rotated Hamiltonian of the unperturbed system and leads to closed-integral expressions for the time-dependent corrections based on Duhamel's formula. This closed-integral structure of corrections is a hallmark of time-independent logarithmic perturbation theory and is preserved in the present extension. This structure, in particular, provides a computable expression for the instantaneous energy shift. Furthermore, dynamic energy shifts arise naturally within this framework in the form of time-averaged expectation values of pseudopotentials and can be related, for example, to AC Stark shifts and electric polarizabilities. As an illustration, we apply the method to the harmonic oscillator and the hydrogen atom, both driven by a time-dependent laser field. The harmonic oscillator provides a proof of principle for which the exact solution is recovered, while the hydrogen atom illustrates the method applied to atomic systems. Supported by numerical simulations, we demonstrate the applicability to obtain relevant physical observables with high accuracy. The present approach offers a promising alternative for analytical studies of time-dependent multi-photon processes in the perturbative regime.

DOI
10.1088/1751-8121/ae8c92
Comments
First version, 31 pages, 4 figures

Related papers

Logarithmic perturbation expansions in nonrelativistic quantum mechanics

Tom D. Imbo, U. Sukhatme

We discuss an alternative method of performing coupling constant perturbation expansions in nonrelativistic quantum mechanics. This method, called logarithmic perturbation theory, yields new expressions for any-order corrections En to an unperturbed bound-state energy which do not involve cumbersome sums over intermediate unperturbed states. In one dimension, these corrections En can be evaluated using a simple explicit form containing a small number of integrals. In more than one dimension the approach is systematic but computations require the solution of well-defined partial differential equations. For order n=2 this equation is identical to that appearing in the method of Dalgarno and Lewis. Numerous illustrative examples are presented.

American Journal of Physics•Feb 1, 1984•44 citations

Aspects of Time-Dependent Perturbation Theory

Peter W. Langhoff, Saul T. Epstein, Martin Karplus

The Dirac variation-of-constants method has long provided a basis for perturbative solution of the time-dependent Schr\"odinger equation. In spite of its widespread utilization, certain aspects of the method have been discussed only superficially and remain somewhat obscure. The present review attempts to clarify some of these points, particularly those related to secular and normalization terms. Secular terms arise from an over-all time-dependent phase in the wave function, while normalization terms preserve the norm of the wave function. A proper treatment of the secular terms is essential in the presence of a physically significant level shift that can produce secular divergences in the time-dependent perturbation functions. The normalization terms are always important, although the formulation of a simple method for including them is of greatest utility in applications requiring higher-order perturbation theory (e.g., nonlinear optical phenomena), where difficulties have arisen in previous treatments. Although the Dirac perturbation technique includes the correct secular and normalization terms when properly executed, it is convenient to reinterpret the perturbation problem so that the secular and normalization terms can be factored from the wave function to all orders. It is shown that an appropriate over-all multiplicative, time-dependent normalization and phase factor can be obtained, and that it is simply the amplitude for finding the system in the unperturbed eigenstate at any time $t$. The regular part of the wave function remaining after this factorization provides a complete description of the physical properties of the system of interest and determines the over-all normalization and phase, as well. Most important, the regular function and its perturbation expansion satisfy equations which are more convenient for computational applications than are the customary Dirac equations, and, in contrast to the latter, they reduce directly to the familiar time-independent perturbation equations in the static limit. To illustrate the general development, the model problem of a linearly perturbed harmonic oscillator and the static, harmonic, and electromagnetic perturbations of arbitrary quantum-mechanical systems are treated explicitly. In the case of an adiabatically applied static perturbation, the familiar adiabatic theorem is recovered with the over-all phase factor giving the perturbed eigenvalue, while in the case of an harmonic perturbation, the overall phase factor obtained includes the system level shift appropriate for a quasiperiodic state. For an electromagnetic perturbation, compact expressions are obtained for various nonlinear optical susceptibilities in forms suitable for computations. Time-dependent Hartree-Fock approximations are treated explicity to demonstrate that difficulties can arise when normalization and secular terms are not extracted prior to application of the perturbation formalism. Connection is also made with other methods which can be employed to eliminate secular and normalization terms from the wave function; these include a projection procedure and multiple-time-scales perturbation theory. The elimination of secular divergences from the perturbation functions is shown to be important for the construction of a valid Fourier transform. Secular and normalization terms also arise in connection with variational principles for the time-dependent Schr\"odinger equation. By employing the Frenkel variational principle and an ansatz for the total wave function that explicitly isolates the secular and normalization terms, a computationally convenient variational functional is obtained. This form of the Frenkel principle provides a bound to the system level shift induced by an oscillatory perturbation and is equivalent to the Ritz variational principle in the static limit. Explicit expressions for the variational functional in the Hartree-Fock approximations are derived in forms suitable for computational applications to the interactions of radiation and matter.

Reviews of Modern Physics•Jul 1, 1972•619 citations

Quantum mechanics in general quantum systems (I): Exact solution

An Min Wang

Starting from our idea of combining the Feynman path integral spirit and the Dyson series kernel, we find an explicit and general form of time evolution operator that is a $c$-number function and a power series of perturbation including all order approximations in the unperturbed Hamiltonian representation. Based on it, we obtain an exact solution of the Schrödinger equation in general quantum systems independent of time. Comparison of our exact solution with the existed perturbation theory makes some features and significance of our exact solution clear. The conclusions expressly indicate that our exact solution is obviously consistent with the usual time-independent perturbation theory at any order approximation, it explicitly calculates out the expanding coefficients of the unperturbed state in the non-perturbation method, and it fully solves the recurrence equation of the expansion coefficients of final state in the unperturbed Hamiltonian representation from a view of time-dependent perturbation theory. At the same time, the exact solution of the von Neumann equation is also given. Our results can be thought of as theoretical developments of quantum dynamics, and are helpful for understanding the dynamical behavior and related subjects of general quantum systems in both theory and application. Our exact solution, together with its sequence studies on perturbation theory [An Min Wang, quant-ph/0611217] and open system dynamics [An Min Wang, quant-ph/0601051] can be used to establish the foundation of theoretical formulism of quantum mechanics in general quantum systems. Further applications of our exact solution to quantum theory can be expected.

arXiv:quant-ph/0611216•Nov 21, 2006•0 citations

A Time-Dependent Analog of the Ritz Variational Principle in Nonrelativistic Quantum Mechanics

Peter Sigmund, George Basbas

A variational functional has been derived for the expectation value of the energy of a quantum system governed by a time-dependent Hamiltonian. Application to a stationary system perturbed by a time-dependent force shows that a perturbational wave function correct to first order in the perturbation, when inserted into the functional, delivers an expression for the energy that is correct up to third order in the perturbation. The use of the variational functional by means of a trial function is illustrated on the forced harmonic oscillator where well-known results are reproduced easily.

Physica Scripta•Nov 1, 1985•3 citations

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

A topical review on time-independent perturbation theory in one-dimensional quantum systems

Leonardo Lavagna, Sandra Carillo, Massimo Panella

Abstract This review presents time-independent perturbative methods for solving the one-dimensional Schrödinger equation, highlighting representative cases that reveal key aspects of the theory. The focus is on their relevance to quantum computing applications, particularly in systems with finite-dimensional state spaces. The main link between these methods and quantum processing emerges when we consider systems with only a finite number of well-defined energy levels, e.g., two levels corresponding to the ‘zero’ and ‘one’ states of a quantum bit (qubit), for which the Schrödinger equation describes a quantum algorithm corresponding to the unitary evolution of the qubits. Perturbative techniques can then be applied to this evolution to gain deeper insight into the algorithm or to simplify it in a controlled manner, as described in the last sections of this review. We initially discuss representative applications using benchmark cases that exemplify the essential features of perturbation theory while avoiding unnecessary technical complications. The methods presented involve direct approaches to solutions for quantum systems with small potential perturbations, with consideration of indirect and approximation techniques. Key cases examined include the particle in a potential well, the harmonic oscillator, and systems with polynomial or Dirac delta distribution perturbations, each chosen to highlight unique aspects of perturbative analysis. Starting with an overview of the unperturbed Schrödinger equation, we derive canonical solutions and then use perturbative techniques to approximate energy shifts and wavefunction corrections in both non-degenerate and degenerate settings. By combining theoretical discussion with case studies, this review underscores how perturbative methods are applied across a range of elementary quantum mechanical problems especially for applications in the qubit model, which, in many cases of interest, can be reduced to the explicit calculations shown here.

Physica Scripta•Sep 23, 2025•0 citations
We present a time-dependent extension of logarithmic perturbation theory for nonrelativistic quantum dynamics governed by the Schrödinger equation, in which the logarithm of the wave function is expanded in powers of a coupling constant. The resulting hierarchy of equations defining the perturbative corrections is governed by a gauge-rotated Hamiltonian of the unperturbed system and leads to closed-integral expressions for the time-dependent corrections based on Duhamel's formula. This closed-integral structure of corrections is a hallmark of time-independent logarithmic perturbation theory and is preserved in the present extension. This structure, in particular, provides a computable expression for the instantaneous energy shift. Furthermore, dynamic energy shifts arise naturally within this framework in the form of time-averaged expectation values of pseudopotentials and can be related, for example, to AC Stark shifts and electric polarizabilities. As an illustration, we apply the method to the harmonic oscillator and the hydrogen atom, both driven by a time-dependent laser field. The harmonic oscillator provides a proof of principle for which the exact solution is recovered, while the hydrogen atom illustrates the method applied to atomic systems. Supported by numerical simulations, we demonstrate the applicability to obtain relevant physical observables with high accuracy. The present approach offers a promising alternative for analytical studies of time-dependent multi-photon processes in the perturbative regime.
We discuss an alternative method of performing coupling constant perturbation expansions in nonrelativistic quantum mechanics. This method, called logarithmic perturbation theory, yields new expressions for any-order corrections En to an unperturbed bound-state energy which do not involve cumbersome sums over intermediate unperturbed states. In one dimension, these corrections En can be evaluated using a simple explicit form containing a small number of integrals. In more than one dimension the approach is systematic but computations require the solution of well-defined partial differential equations. For order n=2 this equation is identical to that appearing in the method of Dalgarno and Lewis. Numerous illustrative examples are presented.
The Dirac variation-of-constants method has long provided a basis for perturbative solution of the time-dependent Schr\"odinger equation. In spite of its widespread utilization, certain aspects of the method have been discussed only superficially and remain somewhat obscure. The present review attempts to clarify some of these points, particularly those related to secular and normalization terms. Secular terms arise from an over-all time-dependent phase in the wave function, while normalization terms preserve the norm of the wave function. A proper treatment of the secular terms is essential in the presence of a physically significant level shift that can produce secular divergences in the time-dependent perturbation functions. The normalization terms are always important, although the formulation of a simple method for including them is of greatest utility in applications requiring higher-order perturbation theory (e.g., nonlinear optical phenomena), where difficulties have arisen in previous treatments. Although the Dirac perturbation technique includes the correct secular and normalization terms when properly executed, it is convenient to reinterpret the perturbation problem so that the secular and normalization terms can be factored from the wave function to all orders. It is shown that an appropriate over-all multiplicative, time-dependent normalization and phase factor can be obtained, and that it is simply the amplitude for finding the system in the unperturbed eigenstate at any time ttt. The regular part of the wave function remaining after this factorization provides a complete description of the physical properties of the system of interest and determines the over-all normalization and phase, as well. Most important, the regular function and its perturbation expansion satisfy equations which are more convenient for computational applications than are the customary Dirac equations, and, in contrast to the latter, they reduce directly to the familiar time-independent perturbation equations in the static limit. To illustrate the general development, the model problem of a linearly perturbed harmonic oscillator and the static, harmonic, and electromagnetic perturbations of arbitrary quantum-mechanical systems are treated explicitly. In the case of an adiabatically applied static perturbation, the familiar adiabatic theorem is recovered with the over-all phase factor giving the perturbed eigenvalue, while in the case of an harmonic perturbation, the overall phase factor obtained includes the system level shift appropriate for a quasiperiodic state. For an electromagnetic perturbation, compact expressions are obtained for various nonlinear optical susceptibilities in forms suitable for computations. Time-dependent Hartree-Fock approximations are treated explicity to demonstrate that difficulties can arise when normalization and secular terms are not extracted prior to application of the perturbation formalism. Connection is also made with other methods which can be employed to eliminate secular and normalization terms from the wave function; these include a projection procedure and multiple-time-scales perturbation theory. The elimination of secular divergences from the perturbation functions is shown to be important for the construction of a valid Fourier transform. Secular and normalization terms also arise in connection with variational principles for the time-dependent Schr\"odinger equation. By employing the Frenkel variational principle and an ansatz for the total wave function that explicitly isolates the secular and normalization terms, a computationally convenient variational functional is obtained. This form of the Frenkel principle provides a bound to the system level shift induced by an oscillatory perturbation and is equivalent to the Ritz variational principle in the static limit. Explicit expressions for the variational functional in the Hartree-Fock approximations are derived in forms suitable for computational applications to the interactions of radiation and matter.
Starting from our idea of combining the Feynman path integral spirit and the Dyson series kernel, we find an explicit and general form of time evolution operator that is a ccc-number function and a power series of perturbation including all order approximations in the unperturbed Hamiltonian representation. Based on it, we obtain an exact solution of the Schrödinger equation in general quantum systems independent of time. Comparison of our exact solution with the existed perturbation theory makes some features and significance of our exact solution clear. The conclusions expressly indicate that our exact solution is obviously consistent with the usual time-independent perturbation theory at any order approximation, it explicitly calculates out the expanding coefficients of the unperturbed state in the non-perturbation method, and it fully solves the recurrence equation of the expansion coefficients of final state in the unperturbed Hamiltonian representation from a view of time-dependent perturbation theory. At the same time, the exact solution of the von Neumann equation is also given. Our results can be thought of as theoretical developments of quantum dynamics, and are helpful for understanding the dynamical behavior and related subjects of general quantum systems in both theory and application. Our exact solution, together with its sequence studies on perturbation theory [An Min Wang, quant-ph/0611217] and open system dynamics [An Min Wang, quant-ph/0601051] can be used to establish the foundation of theoretical formulism of quantum mechanics in general quantum systems. Further applications of our exact solution to quantum theory can be expected.
A variational functional has been derived for the expectation value of the energy of a quantum system governed by a time-dependent Hamiltonian. Application to a stationary system perturbed by a time-dependent force shows that a perturbational wave function correct to first order in the perturbation, when inserted into the functional, delivers an expression for the energy that is correct up to third order in the perturbation. The use of the variational functional by means of a trial function is illustrated on the forced harmonic oscillator where well-known results are reproduced easily.
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.
Abstract This review presents time-independent perturbative methods for solving the one-dimensional Schrödinger equation, highlighting representative cases that reveal key aspects of the theory. The focus is on their relevance to quantum computing applications, particularly in systems with finite-dimensional state spaces. The main link between these methods and quantum processing emerges when we consider systems with only a finite number of well-defined energy levels, e.g., two levels corresponding to the ‘zero’ and ‘one’ states of a quantum bit (qubit), for which the Schrödinger equation describes a quantum algorithm corresponding to the unitary evolution of the qubits. Perturbative techniques can then be applied to this evolution to gain deeper insight into the algorithm or to simplify it in a controlled manner, as described in the last sections of this review. We initially discuss representative applications using benchmark cases that exemplify the essential features of perturbation theory while avoiding unnecessary technical complications. The methods presented involve direct approaches to solutions for quantum systems with small potential perturbations, with consideration of indirect and approximation techniques. Key cases examined include the particle in a potential well, the harmonic oscillator, and systems with polynomial or Dirac delta distribution perturbations, each chosen to highlight unique aspects of perturbative analysis. Starting with an overview of the unperturbed Schrödinger equation, we derive canonical solutions and then use perturbative techniques to approximate energy shifts and wavefunction corrections in both non-degenerate and degenerate settings. By combining theoretical discussion with case studies, this review underscores how perturbative methods are applied across a range of elementary quantum mechanical problems especially for applications in the qubit model, which, in many cases of interest, can be reduced to the explicit calculations shown here.