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

Physical Review A

Multiple-time-scale perturbation theory applied to laser excitation of atoms and molecules

Jack Wong, J. C. Garrison, T. H. Einwohner

Published Feb 1, 1976• 29 citations
Publisher

Abstract

The multiple-linear-time-scale method is used to construct a perturbation theory for $N$-level quantum systems subjected to time-dependent perturbations. The secular and small-denominator terms which plague conventional time-dependent perturbation theory are avoided; consequently, the theory is useful for treating long-term behavior and resonant interactions. The introduction of a self-energy operator allows the shifts in energy levels to be displayed explicitly. When the perturbation is the dipole interaction with an electromagnetic field, the successive approximations yield the well-known rotating-wave approximation, corrections due to counter rotating terms, and the Bloch-Siegert frequency shift. The formalism is applicable to arbitrary pulse shapes, provided that $\ensuremath{\parallel}{H}_{1}\ensuremath{\parallel}\ensuremath{\ll}\ensuremath{\hbar}\overline{\ensuremath{\omega}}$, where $\ensuremath{\parallel}{H}_{1}\ensuremath{\parallel}$ represents the inter-action strength and $\overline{\ensuremath{\omega}}$ is a characteristic frequency of the field. Rabi oscillations induced by multiphoton resonances are automatically included, and this effect is demonstrated in the case of a square pulse with frequency approximately one half the resonance frequency for a transition between two vibrational levels of a molecule. These calculations are compared to an exact numerical solution in order to find the limits of validity of the approximation. Even at high intensities ($I\ensuremath{\sim}{10}^{14}$ W/${\mathrm{cm}}^{2}$) the shifted resonance frequency is quite accurate; however, the approximate and exact solutions are slightly out of phase, so that the approximate solution can only be trusted for a few Rabi cycles. For much lower intensities ($I\ensuremath{\sim}{10}^{10}$ W/${\mathrm{cm}}^{2}$), the approximate solution is valid for many thousands of Rabi cycles.

DOI
10.1103/physreva.13.674

Related papers

Perturbative corrections to the rotating-wave approximation for two-level molecules and the effects of permanent dipoles on single-photon and multiphoton spectra

Mary Ann Kmetic, William J. Meath

Perturbative corrections are derived for the rotating-wave approximation (RWA) for the single-photon and multiphoton resonance profiles due to the interaction between a two-level molecule, with nonzero permanent dipoles, and a sinusoidal time-dependent electric field. The derivation is carried out through the use of a Floquet secular equation that introduces the effects of the permanent dipole moments into the problem through the argument, d\ensuremath{\cdot}e^scrE/\ensuremath{\omega}, of Bessel functions contained in molecule-laser coupling parameters; here d=${\ensuremath{\mu}}_{22}$-${\ensuremath{\mu}}_{11}$, where ${\ensuremath{\mu}}_{\mathrm{ii}}$ is the permanent dipole of molecular state i, and scrE, e^, and \ensuremath{\omega} are the field strength, direction of polarization unit vector, and the circular frequency of the electromagnetic field (EMF). Expansions for the N-photon resonance profiles and the associated resonance frequencies and full widths at half maximum, in powers of the couplings between the transition and permanent dipole moments and the applied EMF, are obtained by expanding the Bessel functions in the perturbative results and are compared with those in the literature (available mostly for d=0 only). The perturbative results, and a series of exactly calculated two-level model absorption spectra, are used to discuss (1) the ability of the explicit perturbative corrections to the RWA to explain the differences between exact and RWA resonance profiles (for example, the positive or negative Bloch-Siegert shifts and dynamic backgrounds absent in the RWA), (2) the usefulness of the perturbative corrections in general, and (3) the effects of permanent dipoles on single-photon and multiphoton absorption spectra.

Physical Review A•Feb 1, 1990•51 citations

Nonperturbative analysis of the two-level atom: Applications to multiphoton excitation

R. E. Duvall, E. Valeo, Carl R. Oberman

Excitation of an atomic system subjected to a slowly varying external electromagnetic field is studied using a two-level model. Time evolution of the system is found using an approach which is nonperturbative in the field strength. There is no constraint to small values of the applied field, that is, the field (in appropriate energy units) need not be small compared to the difference in energies of the two levels. Rather, we use the fact that the situation of interest to us is where the frequency of the exciting field is small compared to the frequency associated with the level difference. Transition probabilities and resonance conditions are found which circumscribe both the large- and small-field limits. In the weak-field limit the previous results of high-order perturbation theory are readily recovered. For a monochromatic field the characteristic features of resonance excitation at high harmonic number of the applied field are (a) extremely narrow resonance widths and (b) shifts in resonance positions which are strong functions of field intensity. Because of this sensitivity, we are able to demonstrate that when slow temporal evolution of the field amplitude is taken into account (e.g., due to finite pulse duration) the appropriate mean excitation rate is that due to the uncorrelated contribution of many resonances. The results of this analysis are used to estimate excitation rates in a specific atomic system, ${\mathrm{Cd}}^{12+}$, which are then compared to multiphoton-ionization rates. Our calculations suggest that the ionization rate exceeds the excitation rate by several orders of magnitude.

Physical Review A•Jun 1, 1988•24 citations

Multiple-time-scale perturbation theory: Radiative decay of coupled atomic states

G. L. Brooks, Leonard M. Scarfone

The spontaneous radiative decay to the ground state of an atomic system initially in the higher of two excited states, coupled in a radiationless fashion by an external field, is investigated by the method of multiple-time-scale perturbation theory. The coupled differential equations of motion for the probability amplitudes are solved to second order, while the usual secular behavior of conventional time-dependent perturbation theory is eliminated in terms of two natural time scales: a fast time scale corresponding to the reciprocal of the level transition frequency and a slower time scale associated with the inverse of the radiation linewidth. Expressions are obtained to this order for the energy shifts and decay characteristics of the excited states and the energy distribution of the final-state photons, quantities previously determined by Fourier analysis and contour integration and by the phenomenological approach of Weisskopf and Wigner. A detailed comparison of these various methods of solution includes conditions under which the present results agree to second order with the previous calculations.

Physical Review A•Dec 1, 1982•4 citations

Adiabatic following in two-photon transition

Munir Hasan Nayfeh, Ali Hasan Nayfeh

The coherent interaction of two smoothly varying, near-resonant, two-photon pulses with a three-level system can be described by "two-photon damped Bloch equations" which are analogous to those for a one-photon transition in a two-level system except for the presence of a two-photon coupling and a frequency shift. These equations are solved for the cases ${\ensuremath{\gamma}}_{1}$, ${\ensuremath{\gamma}}_{2}\ensuremath{\ll}\ensuremath{\Omega}$, ${\ensuremath{\gamma}}_{1}={\ensuremath{\gamma}}_{2}$, and $\frac{{\ensuremath{\gamma}}_{2}{k}^{2}{\ensuremath{\epsilon}}^{4}}{{\ensuremath{\Omega}}^{2}}$, ${\ensuremath{\gamma}}_{1}\ensuremath{\ll}\ensuremath{\Omega}$, where ${\ensuremath{\gamma}}_{1}$ and ${\ensuremath{\gamma}}_{2}$ are the atomic energy and phase relaxation widths, respectively, and $\ensuremath{\Omega}$ is the Rabi frequency. The leading contribution to the refractive index is intensity dependent, caused by the level shifts inherent in multiphoton processes; it includes a relaxation dependent part which is important at times shorter than ${\ensuremath{\gamma}}_{1}^{\ensuremath{-}1}$. The second-order contributions depend on the square of the intensity and the time-integrated square of the intensity. The latter contribution, which is relaxation dependent, causes line asymmetry at the long-wavelength wing; it consists of a term proportional to ${\ensuremath{\gamma}}_{2}\ensuremath{-}{\ensuremath{\gamma}}_{1}$ and only important at early times and a term proportional to $2{\ensuremath{\gamma}}_{2}\ensuremath{-}{\ensuremath{\gamma}}_{1}$.

Physical Review A•Mar 1, 1977•10 citations

Nonlinear resonance effects in multilevel quantum systems weakly coupled to a thermal reservoir. II. Multiple-quantum transitions

Clayton W. Bates

This is the second in a series of three papers on nonlinear resonance effects in multilevel quantum systems in which quantum-mechanical systems are considered that are under the influence of any number of classically described driving fields (well defined in amplitude and phase) and are in weak contact with a thermal reservoir, such that relaxation may be described by means of phenomenological decay constants. In the first paper, referred to as I, a projection-operator formalism was developed for separating the equation of motion for the density operator into a part which may be treated by conventional time-dependent perturbation techniques and another part which, because it would give rise to secular terms using a time-dependent perturbation theory, must be treated by more exact means. In this paper three applications of the formalism developed in I are considered. The first application is to a system in which a multiple-quantum transition is taking place between two levels of the system which contains many levels between the two levels under consideration. In particular, the transition probability is calculated for these two levels, which are not directly connected by matrix elements of the interaction but are connected indirectly through intermediate states. The next application considered is the Raman effect in a four-level system. As a final application of the formalism developed in I an expression is derived for the magnetization ${M}_{+}(t)={M}_{x}(t)+i{M}_{y}(t)$ of a system of $N$ noninteracting spins under the influence of a static field in the $z$ direction and an arbitrary number of rf fields in the $x$ direction. As in I, a phenomenological form is assumed for the relaxation in all of these applications. It is shown in this problem that fields which are not involved in directly causing transitions in this system give rise to frequency shifts in the resonance conditions, which are accounted for in this treatment in a very natural and straightforward manner. For these fields, the rotating as well as the counter-rotating fields give shifts in the resonance condition.

Physical Review A•Dec 1, 1974•0 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
The multiple-linear-time-scale method is used to construct a perturbation theory for NNN-level quantum systems subjected to time-dependent perturbations. The secular and small-denominator terms which plague conventional time-dependent perturbation theory are avoided; consequently, the theory is useful for treating long-term behavior and resonant interactions. The introduction of a self-energy operator allows the shifts in energy levels to be displayed explicitly. When the perturbation is the dipole interaction with an electromagnetic field, the successive approximations yield the well-known rotating-wave approximation, corrections due to counter rotating terms, and the Bloch-Siegert frequency shift. The formalism is applicable to arbitrary pulse shapes, provided that \ensuremath∥H1\ensuremath∥\ensuremath≪\ensuremathℏ\ensuremathω‾\ensuremath{\parallel}{H}_{1}\ensuremath{\parallel}\ensuremath{\ll}\ensuremath{\hbar}\overline{\ensuremath{\omega}}\ensuremath∥H1​\ensuremath∥\ensuremath≪\ensuremathℏ\ensuremathω​, where \ensuremath∥H1\ensuremath∥\ensuremath{\parallel}{H}_{1}\ensuremath{\parallel}\ensuremath∥H1​\ensuremath∥ represents the inter-action strength and \ensuremathω‾\overline{\ensuremath{\omega}}\ensuremathω​ is a characteristic frequency of the field. Rabi oscillations induced by multiphoton resonances are automatically included, and this effect is demonstrated in the case of a square pulse with frequency approximately one half the resonance frequency for a transition between two vibrational levels of a molecule. These calculations are compared to an exact numerical solution in order to find the limits of validity of the approximation. Even at high intensities (I\ensuremath∼1014I\ensuremath{\sim}{10}^{14}I\ensuremath∼1014 W/cm2{\mathrm{cm}}^{2}cm2) the shifted resonance frequency is quite accurate; however, the approximate and exact solutions are slightly out of phase, so that the approximate solution can only be trusted for a few Rabi cycles. For much lower intensities (I\ensuremath∼1010I\ensuremath{\sim}{10}^{10}I\ensuremath∼1010 W/cm2{\mathrm{cm}}^{2}cm2), the approximate solution is valid for many thousands of Rabi cycles.
Perturbative corrections are derived for the rotating-wave approximation (RWA) for the single-photon and multiphoton resonance profiles due to the interaction between a two-level molecule, with nonzero permanent dipoles, and a sinusoidal time-dependent electric field. The derivation is carried out through the use of a Floquet secular equation that introduces the effects of the permanent dipole moments into the problem through the argument, d\ensuremath{\cdot}e^scrE/\ensuremath{\omega}, of Bessel functions contained in molecule-laser coupling parameters; here d=${\ensuremath{\mu}}_{22}$-${\ensuremath{\mu}}_{11}$, where ${\ensuremath{\mu}}_{\mathrm{ii}}$ is the permanent dipole of molecular state i, and scrE, e^, and \ensuremath{\omega} are the field strength, direction of polarization unit vector, and the circular frequency of the electromagnetic field (EMF). Expansions for the N-photon resonance profiles and the associated resonance frequencies and full widths at half maximum, in powers of the couplings between the transition and permanent dipole moments and the applied EMF, are obtained by expanding the Bessel functions in the perturbative results and are compared with those in the literature (available mostly for d=0 only). The perturbative results, and a series of exactly calculated two-level model absorption spectra, are used to discuss (1) the ability of the explicit perturbative corrections to the RWA to explain the differences between exact and RWA resonance profiles (for example, the positive or negative Bloch-Siegert shifts and dynamic backgrounds absent in the RWA), (2) the usefulness of the perturbative corrections in general, and (3) the effects of permanent dipoles on single-photon and multiphoton absorption spectra.
Excitation of an atomic system subjected to a slowly varying external electromagnetic field is studied using a two-level model. Time evolution of the system is found using an approach which is nonperturbative in the field strength. There is no constraint to small values of the applied field, that is, the field (in appropriate energy units) need not be small compared to the difference in energies of the two levels. Rather, we use the fact that the situation of interest to us is where the frequency of the exciting field is small compared to the frequency associated with the level difference. Transition probabilities and resonance conditions are found which circumscribe both the large- and small-field limits. In the weak-field limit the previous results of high-order perturbation theory are readily recovered. For a monochromatic field the characteristic features of resonance excitation at high harmonic number of the applied field are (a) extremely narrow resonance widths and (b) shifts in resonance positions which are strong functions of field intensity. Because of this sensitivity, we are able to demonstrate that when slow temporal evolution of the field amplitude is taken into account (e.g., due to finite pulse duration) the appropriate mean excitation rate is that due to the uncorrelated contribution of many resonances. The results of this analysis are used to estimate excitation rates in a specific atomic system, Cd12+{\mathrm{Cd}}^{12+}Cd12+, which are then compared to multiphoton-ionization rates. Our calculations suggest that the ionization rate exceeds the excitation rate by several orders of magnitude.
The spontaneous radiative decay to the ground state of an atomic system initially in the higher of two excited states, coupled in a radiationless fashion by an external field, is investigated by the method of multiple-time-scale perturbation theory. The coupled differential equations of motion for the probability amplitudes are solved to second order, while the usual secular behavior of conventional time-dependent perturbation theory is eliminated in terms of two natural time scales: a fast time scale corresponding to the reciprocal of the level transition frequency and a slower time scale associated with the inverse of the radiation linewidth. Expressions are obtained to this order for the energy shifts and decay characteristics of the excited states and the energy distribution of the final-state photons, quantities previously determined by Fourier analysis and contour integration and by the phenomenological approach of Weisskopf and Wigner. A detailed comparison of these various methods of solution includes conditions under which the present results agree to second order with the previous calculations.
The coherent interaction of two smoothly varying, near-resonant, two-photon pulses with a three-level system can be described by "two-photon damped Bloch equations" which are analogous to those for a one-photon transition in a two-level system except for the presence of a two-photon coupling and a frequency shift. These equations are solved for the cases ${\ensuremath{\gamma}}_{1}$, ${\ensuremath{\gamma}}_{2}\ensuremath{\ll}\ensuremath{\Omega}$, ${\ensuremath{\gamma}}_{1}={\ensuremath{\gamma}}_{2}$, and $\frac{{\ensuremath{\gamma}}_{2}{k}^{2}{\ensuremath{\epsilon}}^{4}}{{\ensuremath{\Omega}}^{2}}$, ${\ensuremath{\gamma}}_{1}\ensuremath{\ll}\ensuremath{\Omega}$, where ${\ensuremath{\gamma}}_{1}$ and ${\ensuremath{\gamma}}_{2}$ are the atomic energy and phase relaxation widths, respectively, and $\ensuremath{\Omega}$ is the Rabi frequency. The leading contribution to the refractive index is intensity dependent, caused by the level shifts inherent in multiphoton processes; it includes a relaxation dependent part which is important at times shorter than ${\ensuremath{\gamma}}_{1}^{\ensuremath{-}1}$. The second-order contributions depend on the square of the intensity and the time-integrated square of the intensity. The latter contribution, which is relaxation dependent, causes line asymmetry at the long-wavelength wing; it consists of a term proportional to ${\ensuremath{\gamma}}_{2}\ensuremath{-}{\ensuremath{\gamma}}_{1}$ and only important at early times and a term proportional to $2{\ensuremath{\gamma}}_{2}\ensuremath{-}{\ensuremath{\gamma}}_{1}$.
This is the second in a series of three papers on nonlinear resonance effects in multilevel quantum systems in which quantum-mechanical systems are considered that are under the influence of any number of classically described driving fields (well defined in amplitude and phase) and are in weak contact with a thermal reservoir, such that relaxation may be described by means of phenomenological decay constants. In the first paper, referred to as I, a projection-operator formalism was developed for separating the equation of motion for the density operator into a part which may be treated by conventional time-dependent perturbation techniques and another part which, because it would give rise to secular terms using a time-dependent perturbation theory, must be treated by more exact means. In this paper three applications of the formalism developed in I are considered. The first application is to a system in which a multiple-quantum transition is taking place between two levels of the system which contains many levels between the two levels under consideration. In particular, the transition probability is calculated for these two levels, which are not directly connected by matrix elements of the interaction but are connected indirectly through intermediate states. The next application considered is the Raman effect in a four-level system. As a final application of the formalism developed in I an expression is derived for the magnetization M+(t)=Mx(t)+iMy(t){M}_{+}(t)={M}_{x}(t)+i{M}_{y}(t)M+​(t)=Mx​(t)+iMy​(t) of a system of NNN noninteracting spins under the influence of a static field in the zzz direction and an arbitrary number of rf fields in the xxx direction. As in I, a phenomenological form is assumed for the relaxation in all of these applications. It is shown in this problem that fields which are not involved in directly causing transitions in this system give rise to frequency shifts in the resonance conditions, which are accounted for in this treatment in a very natural and straightforward manner. For these fields, the rotating as well as the counter-rotating fields give shifts in the resonance condition.
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.