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

Tutorial: Perturbation Approximations to Asymptotic Quantum Transition Rates

Steven Kenneth Kauffmann

Published Aug 1, 2019• 0 citations
OpenAlex

Abstract

The orders of the perturbation approximation to the Schroedinger-picture time evolution operator in powers of the presumed-small perturbation part of the Hamiltonian operator are developed by the iteration of an identity; the possibly confusing switch to the interaction picture isn't needed. Those time-evolution operator approximations are sandwiched between two orthogonal, normalized eigenstates of the unperturbed part of the Hamiltonian operator to produce the transition amplitude approximations, whose calculation is reduced to quadrature when every occurrence of the perturbation part of the Hamiltonian operator in them is expanded in the unperturbed basis. That expansion also reveals the time-dependent parts of those approximations to be multiply nested integrals which in the long-time limit approach simple products of non-singular inverses (principle value plus delta function) of differences of unperturbed energy eigenvalues. Closely related to the long-time limits of transition amplitudes are the long-time averages of transition rates. When the unperturbed basis is that of free-particle states, sums over relevant final states of transition rates from an initial state, divided by the initial state's particle flux, produce cross sections. Here the perturbation approximations to generic quantum transition rates are parlayed to the corresponding approximations to differential cross sections for nonrelativistic-particle potential scattering.

Related papers

Light-Matter Interactions

Ying Fu, Min Qiu

We start the discussion by a formal theory of perturbation, which will be further elucidated by straightforward calculations. 2.1.1 The Golden Rule Assume that H0 is a time-independent Hamiltonian and all solutions of the eigenvalue problem H0|k〉 = Ek|k〉 are known; it is straightforward to construct the solution of the equation of motion i ∂|ψ(t)〉 ∂t = H0|ψ(t)〉 (2.1) which coincides at t = 0 with the given initial state |ψ(0)〉. Under these circumstances, the state vector |ψ(t)〉 = e−iH0t/|ψ(0) (2.2) represents the general solution of the equation of motion, Eq. 2.1. Since∑ |k〉〈k| = 1 we obtain |ψ(t)〉 = ∑ k e−iH0t/|k〉〈k|ψ(0)〉 = ∑ k e−iEkt/|k〉〈k|ψ(0)〉 (2.3) We now calculate transition amplitudes between the relevant unperturbed eigenstates due to the appearance of an external perturbation V at t = 0. We transform the equation of motion in the Schro¨dinger picture i ∂|ψ(t)〉 ∂t = (H0 + V )|ψ(t)〉 (2.4) to the interaction picture by a time-dependent unitary operator |ψˆ(t)〉 = exp (iH0t/)|ψ(t)〉 (2.5) so that in the interaction picture, the equation of motion for the state is i ∂|ψˆ(t)〉 ∂t = Vˆ (t)|ψˆ(t)〉 (2.6) where the new interaction operator is given by Vˆ (t) = exp (iH0t/)V exp (−iH0t/) (2.7) Next we introduce a time development operator Tˆ (t) such that |ψˆ(t)〉 = Tˆ (t)|ψˆ(0)〉 (2.8) The time development operator satisfies the integral equation Tˆ (t) = 1− i Vˆ (t)Tˆ (t′)dt′ (2.9) The relationship between the operator expressions in the two pictures is Tˆ (t) = exp (iH0t/)T (t) exp (−iH0t/) (2.10) The transition matrix elements of the time development operator between eigen- states of the unperturbed Hamiltonian are 〈q|Tˆ (t)|k〉 = δ(q − k)− i 〈q|Vˆ (t)Tˆ (t′)|k〉dt′ = δ(q − k)− i 〈q|Vˆ (t′)|k′〉〈k′|Tˆ (t′)|k〉dt′ (2.11) By Eq. 2.10, and note H0|k〉 = Ek|k〉 and H0|q〉 = Eq|q〉, 〈q|Tˆ (t)|k〉 = δ(q − k)− i ei(Eq−Ek′ )t ′/〈q|V (t′)|k′〉〈k′|Tˆ (t′)|k〉dt′ (2.12) The strategy of time-dependent perturbation theory is to assume that V is small and proceed with iterating Eq. 2.9 as a power series in terms of Vˆ Tˆ (t) = 1− i Vˆ ′dt′ − 1 2 Vˆ ′′dt′′ + · · · (2.13) and the transition amplitude from initial state |k〉 to final state |q〉 in the perturbation expansion becomes 〈q|Tˆ (t)|k〉 = δ(q − k)− i − 1 2 × ∫ t′ 0 ei(Ek′−Ek)t ′′/〈k′|V (t′′)|k〉dt′dt′′ + · · · (2.14) It is equivalent to obtain the above equation by iteratively replacing the term 〈k′|Tˆ (t′)|k〉 on the right side of Eq. 2.12 by the whole right-side expression of Eq. 2.12, i.e., ∑ k′ = ∫ t′ 0 + ∑ k′ =k ei(Eq−Ek′ )t ′′/〈q|V (t′′)|k′〉〈k′|Tˆ (t′′)|k〉dt′′ (2.15) By Eq. 2.12 〈k|Tˆ (t′)|k〉 = 1− i ei(Ek−Ek′′ )t ′′/〈k|V (t′′)|k′′〉〈k′′|Tˆ (t′′)|k〉dt′′ (2.16) The first term on the right side of the above equation results in the second term on the right side of Eq. 2.14. For k′ = k and by Eq. 2.12, the second term on the right side of Eq. 2.15 becomes 〈k′|Tˆ (t′)|k〉 = − i ei(Ek′−Ek′′ )t ′′/〈k′|V (t′′)|k′′〉〈k′′|Tˆ (t′′)|k〉dt′′ (2.17) We use Eq. 2.15 to separate the term k′′ = k from k′′ = k in the above equation, which becomes the third term on the right side of Eq. 2.14.

•Aug 8, 2011•0 citations

Quantum Dynamics

Uri Peskin

Explicitly time-dependent Hamiltonians are ubiquitous in applications of quantum theory. It is therefore necessary to solve the time-dependent Schrödinger equation directly. The system’s dynamics is associated with a unitary time-evolution operator (a propagator), formally given as an infinite Dyson series. Time-dependent observables are invariant under unitary time-dependent transformations, where it is sometimes useful to transform the time-evolution from the states into the corresponding operators. This is carried out in part (in full) by transforming to the interaction (Heisenberg) picture. The corresponding equations of motion for the time-dependent operators are introduced. For quadratic potential energy functions, the time evolution of quantum expectation values coincides with the corresponding classical dynamics. This is demonstrated and analyzed in detail for Gaussian wave packets and a coherent state. Finally, we derive exact and approximate expressions for time-dependent transition probabilities and transition rates between quantum states. The validity of time-dependent perturbation theory is analyzed by comparison to exact dynamics.

Cambridge University Press eBooks•May 11, 2023•0 citations

Quantum transition

Tao Xiang

A small parameter is used to identify the orders of perturbation. The unperturbed Hamiltonian is assumed time-independent and can be diagonalized. The interaction of an electron with an applied electromagnetic field is governed by the gauge-invariant Hamiltonian. If both the initial and final energy levels are discretized, the delta function indicates that the transition happens only when the energy difference between the two levels matches the frequency of the perturbation. The delta function suggests that the absorption or stimulated emission happens only when the energy difference between the final and initial states matches precisely the frequency of the applied electromagnetic wave. The quantum transition induced by a monochromatic perturbation is determined by the matrix element between the initial and final states.

CRC Press eBooks•Mar 31, 2022•0 citations

Time Dependent Perturbation Theory

Asim Gangopadhyaya, Constantin Rasinariu

In time-independent perturbation theory, we were interested in finding the eigenvalues and eigenstates of a system whose Hamiltonian could be written as the sum of a known Hamiltonian and another “small” term describing the perturbation. In time-dependent perturbation theory, the goal is different. We determine the time evolution operator of the perturbed system and use it to estimate the transition probabilities between two states of the unperturbed system.

CRC Press eBooks•Jul 28, 2025•0 citations

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 orders of the perturbation approximation to the Schroedinger-picture time evolution operator in powers of the presumed-small perturbation part of the Hamiltonian operator are developed by the iteration of an identity; the possibly confusing switch to the interaction picture isn't needed. Those time-evolution operator approximations are sandwiched between two orthogonal, normalized eigenstates of the unperturbed part of the Hamiltonian operator to produce the transition amplitude approximations, whose calculation is reduced to quadrature when every occurrence of the perturbation part of the Hamiltonian operator in them is expanded in the unperturbed basis. That expansion also reveals the time-dependent parts of those approximations to be multiply nested integrals which in the long-time limit approach simple products of non-singular inverses (principle value plus delta function) of differences of unperturbed energy eigenvalues. Closely related to the long-time limits of transition amplitudes are the long-time averages of transition rates. When the unperturbed basis is that of free-particle states, sums over relevant final states of transition rates from an initial state, divided by the initial state's particle flux, produce cross sections. Here the perturbation approximations to generic quantum transition rates are parlayed to the corresponding approximations to differential cross sections for nonrelativistic-particle potential scattering.
We start the discussion by a formal theory of perturbation, which will be further elucidated by straightforward calculations. 2.1.1 The Golden Rule Assume that H0 is a time-independent Hamiltonian and all solutions of the eigenvalue problem H0|k〉 = Ek|k〉 are known; it is straightforward to construct the solution of the equation of motion i ∂|ψ(t)〉 ∂t = H0|ψ(t)〉 (2.1) which coincides at t = 0 with the given initial state |ψ(0)〉. Under these circumstances, the state vector |ψ(t)〉 = e−iH0t/|ψ(0) (2.2) represents the general solution of the equation of motion, Eq. 2.1. Since∑ |k〉〈k| = 1 we obtain |ψ(t)〉 = ∑ k e−iH0t/|k〉〈k|ψ(0)〉 = ∑ k e−iEkt/|k〉〈k|ψ(0)〉 (2.3) We now calculate transition amplitudes between the relevant unperturbed eigenstates due to the appearance of an external perturbation V at t = 0. We transform the equation of motion in the Schro¨dinger picture i ∂|ψ(t)〉 ∂t = (H0 + V )|ψ(t)〉 (2.4) to the interaction picture by a time-dependent unitary operator |ψˆ(t)〉 = exp (iH0t/)|ψ(t)〉 (2.5) so that in the interaction picture, the equation of motion for the state is i ∂|ψˆ(t)〉 ∂t = Vˆ (t)|ψˆ(t)〉 (2.6) where the new interaction operator is given by Vˆ (t) = exp (iH0t/)V exp (−iH0t/) (2.7) Next we introduce a time development operator Tˆ (t) such that |ψˆ(t)〉 = Tˆ (t)|ψˆ(0)〉 (2.8) The time development operator satisfies the integral equation Tˆ (t) = 1− i Vˆ (t)Tˆ (t′)dt′ (2.9) The relationship between the operator expressions in the two pictures is Tˆ (t) = exp (iH0t/)T (t) exp (−iH0t/) (2.10) The transition matrix elements of the time development operator between eigen- states of the unperturbed Hamiltonian are 〈q|Tˆ (t)|k〉 = δ(q − k)− i 〈q|Vˆ (t)Tˆ (t′)|k〉dt′ = δ(q − k)− i 〈q|Vˆ (t′)|k′〉〈k′|Tˆ (t′)|k〉dt′ (2.11) By Eq. 2.10, and note H0|k〉 = Ek|k〉 and H0|q〉 = Eq|q〉, 〈q|Tˆ (t)|k〉 = δ(q − k)− i ei(Eq−Ek′ )t ′/〈q|V (t′)|k′〉〈k′|Tˆ (t′)|k〉dt′ (2.12) The strategy of time-dependent perturbation theory is to assume that V is small and proceed with iterating Eq. 2.9 as a power series in terms of Vˆ Tˆ (t) = 1− i Vˆ ′dt′ − 1 2 Vˆ ′′dt′′ + · · · (2.13) and the transition amplitude from initial state |k〉 to final state |q〉 in the perturbation expansion becomes 〈q|Tˆ (t)|k〉 = δ(q − k)− i − 1 2 × ∫ t′ 0 ei(Ek′−Ek)t ′′/〈k′|V (t′′)|k〉dt′dt′′ + · · · (2.14) It is equivalent to obtain the above equation by iteratively replacing the term 〈k′|Tˆ (t′)|k〉 on the right side of Eq. 2.12 by the whole right-side expression of Eq. 2.12, i.e., ∑ k′ = ∫ t′ 0 + ∑ k′ =k ei(Eq−Ek′ )t ′′/〈q|V (t′′)|k′〉〈k′|Tˆ (t′′)|k〉dt′′ (2.15) By Eq. 2.12 〈k|Tˆ (t′)|k〉 = 1− i ei(Ek−Ek′′ )t ′′/〈k|V (t′′)|k′′〉〈k′′|Tˆ (t′′)|k〉dt′′ (2.16) The first term on the right side of the above equation results in the second term on the right side of Eq. 2.14. For k′ = k and by Eq. 2.12, the second term on the right side of Eq. 2.15 becomes 〈k′|Tˆ (t′)|k〉 = − i ei(Ek′−Ek′′ )t ′′/〈k′|V (t′′)|k′′〉〈k′′|Tˆ (t′′)|k〉dt′′ (2.17) We use Eq. 2.15 to separate the term k′′ = k from k′′ = k in the above equation, which becomes the third term on the right side of Eq. 2.14.
Explicitly time-dependent Hamiltonians are ubiquitous in applications of quantum theory. It is therefore necessary to solve the time-dependent Schrödinger equation directly. The system’s dynamics is associated with a unitary time-evolution operator (a propagator), formally given as an infinite Dyson series. Time-dependent observables are invariant under unitary time-dependent transformations, where it is sometimes useful to transform the time-evolution from the states into the corresponding operators. This is carried out in part (in full) by transforming to the interaction (Heisenberg) picture. The corresponding equations of motion for the time-dependent operators are introduced. For quadratic potential energy functions, the time evolution of quantum expectation values coincides with the corresponding classical dynamics. This is demonstrated and analyzed in detail for Gaussian wave packets and a coherent state. Finally, we derive exact and approximate expressions for time-dependent transition probabilities and transition rates between quantum states. The validity of time-dependent perturbation theory is analyzed by comparison to exact dynamics.
A small parameter is used to identify the orders of perturbation. The unperturbed Hamiltonian is assumed time-independent and can be diagonalized. The interaction of an electron with an applied electromagnetic field is governed by the gauge-invariant Hamiltonian. If both the initial and final energy levels are discretized, the delta function indicates that the transition happens only when the energy difference between the two levels matches the frequency of the perturbation. The delta function suggests that the absorption or stimulated emission happens only when the energy difference between the final and initial states matches precisely the frequency of the applied electromagnetic wave. The quantum transition induced by a monochromatic perturbation is determined by the matrix element between the initial and final states.
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.