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

Reviews of Modern Physics

Aspects of Time-Dependent Perturbation Theory

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

Published Jul 1, 1972• 619 citations
Publisher

Abstract

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.

DOI
10.1103/revmodphys.44.602

Related papers

A modification of the Dirac time-dependent perturbation theory

Kaushik Bhattacharya

Abstract Within the framework of the method of variation of constants, time-dependent perturbation theory is presented in a form that naturally eliminates the undesirable ‘secular ’ terms for a time-independent perturbation, permits an order-by-order calculation of the exact wavefunction for a ‘resonant’ harmonic perturbation and offers a simple route to establish the adiabatic hypothesis in quantum mechanics. The formulation rests on the introduction of a flexibility in the choice of the phase factors associated with the varying amplitudes, together with development of a perturbation procedure that closely resembles the Brillouin-Wigner formalism in the static case, and which exploits the flexibility of the phase factors at each order of the theory.

Proceedings of the Royal Society of London A Mathematical and Physical Sciences•Aug 8, 1984•12 citations

New Formulation of Time-Dependent Perturbation Theory

Kwong T. Chung

The previous time-dependent variation-perturbation theory, which is derived from the conventional time-dependent perturbation equations, contains two inconvenient features. When the perturbation becomes static, the wave function so obtained does not satisfy the time-dependent Schr\"odinger equation, and the perturbation equations cannot be reduced to the ordinary time-independent perturbation equations. To remove these difficulties, we introduce a perturbed time-dependent phase factor in deriving a new set of perturbation equations. We then construct the variation-perturbation functionals as in the previous theory; the perturbed wave functions are solved by minimizing these functionals. The same method is also applied to many-electron systems, using the coupled and uncoupled Hartree-Fock scheme.

Physical Review•Nov 25, 1967•29 citations

Perturbative quantum dynamics: variants of the Dirac method

Krishnendu Bhattacharyya, Debashis Mukherjee

A detailed comparative study of the structure and workability of various rearranged versions of the Dirac method of variation of constants for perturbative treatment of quantum dynamics is presented. The variants include a modified strategy proposed recently and also a new rearrangement of this form; both of these contain lambda dependent phase factors associated with the amplitudes which are determined order-by-order by following a Brillouin-Wigner type perturbation procedure, lambda being the perturbation parameter. The following classes of perturbation problems are considered: V, (V exp(i omega t)+V Dagger exp(-i omega t)), f(t)V and f(t)(V exp(i omega t)+V Dagger exp(-i omega t)), where f(t) defines some suitable adiabatic 'switching' function. Special attention is paid to (i) the avoidance of divergent parts in the amplitude-correction terms from the scheme, (ii) applicability to cases involving a degenerate initial state and a 'resonant' harmonic perturbation, (iii) use of both exponential and non-exponential forms for f(t) and (iv) possible sources of convergence difficulties in such variants. The merits and demerits of the schemes concerned for the different problems under investigation are also listed, for convenience, in a tabular form.

Journal of Physics A Mathematical and General•Jan 11, 1986•7 citations

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

Juan Carlos del Valle, Paul Bergold, Karolina Kropielnicka

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.

quant-pharXiv:2604.14812•Apr 17, 2026

Adiabatic and nonadiabatic contributions to the energy of a system subject to a time-dependent perturbation: Complete separation and physical interpretation

Anirban Mandal, Katharine L. C. Hunt

When a time-dependent perturbation acts on a quantum system that is initially in the nondegenerate ground state ∣0> of an unperturbed Hamiltonian H(0), the wave function acquires excited-state components ∣k> with coefficients c(k)(t) exp(-iE(k)t/ℏ), where E(k) denotes the energy of the unperturbed state ∣k>. It is well known that each coefficient c(k)(t) separates into an adiabatic term a(k)(t) that reflects the adjustment of the ground state to the perturbation--without actual transitions--and a nonadiabatic term b(k)(t) that yields the probability amplitude for a transition to the excited state. In this work, we prove that the energy at any time t also separates completely into adiabatic and nonadiabatic components, after accounting for the secular and normalization terms that appear in the solution of the time-dependent Schrödinger equation via Dirac's method of variation of constants. This result is derived explicitly through third order in the perturbation. We prove that the cross-terms between the adiabatic and nonadiabatic parts of c(k)(t) vanish, when the energy at time t is determined as an expectation value. The adiabatic term in the energy is identical to the total energy obtained from static perturbation theory, for a system exposed to the instantaneous perturbation λH'(t). The nonadiabatic term is a sum over excited states ∣k> of the transition probability multiplied by the transition energy. By evaluating the probabilities of transition to the excited eigenstates ∣k'(t)> of the instantaneous Hamiltonian H(t), we provide a physically transparent explanation of the result for E(t). To lowest order in the perturbation parameter λ, the probability of finding the system in state ∣k'(t)> is given by λ(2) ∣b(k)(t)∣(2). At third order, the transition probability depends on a second-order transition coefficient, derived in this work. We indicate expected differences between the results for transition probabilities obtained from this work and from Fermi's golden rule.

The Journal of Chemical Physics•Oct 26, 2012•12 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
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.
Abstract Within the framework of the method of variation of constants, time-dependent perturbation theory is presented in a form that naturally eliminates the undesirable ‘secular ’ terms for a time-independent perturbation, permits an order-by-order calculation of the exact wavefunction for a ‘resonant’ harmonic perturbation and offers a simple route to establish the adiabatic hypothesis in quantum mechanics. The formulation rests on the introduction of a flexibility in the choice of the phase factors associated with the varying amplitudes, together with development of a perturbation procedure that closely resembles the Brillouin-Wigner formalism in the static case, and which exploits the flexibility of the phase factors at each order of the theory.
The previous time-dependent variation-perturbation theory, which is derived from the conventional time-dependent perturbation equations, contains two inconvenient features. When the perturbation becomes static, the wave function so obtained does not satisfy the time-dependent Schr\"odinger equation, and the perturbation equations cannot be reduced to the ordinary time-independent perturbation equations. To remove these difficulties, we introduce a perturbed time-dependent phase factor in deriving a new set of perturbation equations. We then construct the variation-perturbation functionals as in the previous theory; the perturbed wave functions are solved by minimizing these functionals. The same method is also applied to many-electron systems, using the coupled and uncoupled Hartree-Fock scheme.
A detailed comparative study of the structure and workability of various rearranged versions of the Dirac method of variation of constants for perturbative treatment of quantum dynamics is presented. The variants include a modified strategy proposed recently and also a new rearrangement of this form; both of these contain lambda dependent phase factors associated with the amplitudes which are determined order-by-order by following a Brillouin-Wigner type perturbation procedure, lambda being the perturbation parameter. The following classes of perturbation problems are considered: V, (V exp(i omega t)+V Dagger exp(-i omega t)), f(t)V and f(t)(V exp(i omega t)+V Dagger exp(-i omega t)), where f(t) defines some suitable adiabatic 'switching' function. Special attention is paid to (i) the avoidance of divergent parts in the amplitude-correction terms from the scheme, (ii) applicability to cases involving a degenerate initial state and a 'resonant' harmonic perturbation, (iii) use of both exponential and non-exponential forms for f(t) and (iv) possible sources of convergence difficulties in such variants. The merits and demerits of the schemes concerned for the different problems under investigation are also listed, for convenience, in a tabular form.
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.
When a time-dependent perturbation acts on a quantum system that is initially in the nondegenerate ground state ∣0> of an unperturbed Hamiltonian H(0), the wave function acquires excited-state components ∣k> with coefficients c(k)(t) exp(-iE(k)t/ℏ), where E(k) denotes the energy of the unperturbed state ∣k>. It is well known that each coefficient c(k)(t) separates into an adiabatic term a(k)(t) that reflects the adjustment of the ground state to the perturbation--without actual transitions--and a nonadiabatic term b(k)(t) that yields the probability amplitude for a transition to the excited state. In this work, we prove that the energy at any time t also separates completely into adiabatic and nonadiabatic components, after accounting for the secular and normalization terms that appear in the solution of the time-dependent Schrödinger equation via Dirac's method of variation of constants. This result is derived explicitly through third order in the perturbation. We prove that the cross-terms between the adiabatic and nonadiabatic parts of c(k)(t) vanish, when the energy at time t is determined as an expectation value. The adiabatic term in the energy is identical to the total energy obtained from static perturbation theory, for a system exposed to the instantaneous perturbation λH'(t). The nonadiabatic term is a sum over excited states ∣k> of the transition probability multiplied by the transition energy. By evaluating the probabilities of transition to the excited eigenstates ∣k'(t)> of the instantaneous Hamiltonian H(t), we provide a physically transparent explanation of the result for E(t). To lowest order in the perturbation parameter λ, the probability of finding the system in state ∣k'(t)> is given by λ(2) ∣b(k)(t)∣(2). At third order, the transition probability depends on a second-order transition coefficient, derived in this work. We indicate expected differences between the results for transition probabilities obtained from this work and from Fermi's golden rule.
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.