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

Nonlinear polarizations and excitations and their time dependence in discrete multilevel systems

H. B. Puell, C. R. Vidal

Published Dec 1, 1976• 32 citations
Publisher

Abstract

For the interaction of coherent radiation with a nondegenerate, homogeneously broadened multilevel system the nonlinear susceptibilities and excitations and their time dependence are evaluated using a perturbation approach to the density matrix up to fourth order. This is done for the "off-resonance" case, ${(\ensuremath{\Delta}{\ensuremath{\omega}}^{2}+{\ensuremath{\Gamma}}^{2})}^{\frac{1}{2}}\ensuremath{\gg}\frac{1}{\ensuremath{\tau}}$, as well as for the "on-resonance" case, ${(\ensuremath{\Delta}{\ensuremath{\omega}}^{2}+{\ensuremath{\Gamma}}^{2})}^{\frac{1}{2}}\ensuremath{\ll}\frac{1}{\ensuremath{\tau}}$, with $\ensuremath{\tau}$ being the pulse duration, by taking into account the temporal change of the electric field amplitude and the phase and energy relaxation times, ${T}_{2}$ and ${T}_{1}$. It is shown that the transient one-photon excitation known as adiabatic following, is not only important for a pulse duration $\ensuremath{\tau}<{T}_{2}$, but may also be important for ${T}_{2}\ensuremath{\lesssim}\ensuremath{\tau}<{T}_{1}$. In addition, expressions for the transient two-photon excitation in a multilevel system are derived. Finally, the two-level limit is evaluated to all orders of the perturbation including damping effects, which gives results valid also for large incident light intensities. For the "off-resonance" case an expression is obtained which describes not only the transient adiabatic following process, but also the usual one-photon absorption including saturation effects giving rise to power broadening.

DOI
10.1103/physreva.14.2225

Related papers

Time-dependent collisional approximations for two-level systems

Michael E. Crenshaw, Cyrus D. Cantrell

We derive two types of perturbative approximation for the response of collisionally broadened two-level atoms to a time-dependent field. The perturbative approximations for the inversion and polarization are valid to any desired order in the derivatives of the field. Comparison with a numerical solution of the full density-matrix equations of motion demonstrates that the new approximations significantly improve upon the adiabatic-following and collisional steady-state approximations for time-dependent pulses in the appropriate damping regime. Application of the approximations to pulse propagation, sideband instability, and population pulsations is discussed. The collisional steady-state and adiabatic-following approximations find extensive use in many areas of nonlinear optics. The deficiencies of these approximations are well known: The collisional steady-state approxirnation is valid only for monochromatic (cw) radiation and the adiabatic-following approximation is valid only in the absence of damping and only if the product of the detuning and a characteristic time is much larger than one. In spite of their manifest faults, the collisional steady-state and adiabatic-following approximations are often exceptionally good approximations, in the appropriate damping regime, for the local response of a vapor of two-level atoms to a slowly varying field. However, timedependent nonlinear processes modify a propagating field zuch that it varies less slowly and the neglected portions of the atomic response become more important. ' Therefore improved approximations for the damped atomic response to a time-dependent field are desirable. In this paper we present two improved perturbative approximations for the inversion and polarization of collisionally broadened two-level atoms induced by a time-dependent field. The adiabatic-following approximation for two-level systems was derived phenomenologically by Grischkowsky, who used it to explain certain experimental results of near-resonance laser-pulse propagation. Motivated by the work of Grischkowsky, Crisp provided an analytic derivation for the adiabatic-following approximation and obtained a quadrature for the error. The essential feature of Crisp's derivation is the expansion of the kernel of the integral representation of the off-diagonal density-matrix element p2, in a Fourier series with integration proceeding term by term. Both the adiabatic-following approximation and the collisional steady-state approximation were obtained in first order. Both Crisp and Grischkowsky used the adiabatic inversion to derive a correction to the adiabatic-following approximation for the polarization. We have shown how Crisp's method can be extended to obtain corrections to the adiabatic inversion and to improve corrections to the polarization in the absence of damping. In this paper we show how the integrated Fourier series can be used to obtain two types of perturbative approximation for the inversion and polarization of collisionally broadened twolevel systems that are valid to any desired order in the derivatives of the field. The two types of perturbative approximation to the inversion are (i) a formal solution of the rate equation for the inversion which is obtained from the density-matrix equations for the populations; and (ii) an approximation obtained using the approximation (p22 — p») —I — 4p2, P, 2, valid for weak damping. The approximations for the polarization are obtained by substituting either (i) or (ii) for the inversion in the integrated-Fourier-series expansion of p2, . The approximation obtained using the first type of approximate inversion converges most rapidly for strong damping and reduces to the collisional steady-state approximation when the magnitude and phase of the field are constant. The approximation obtained using the second type of approxirnate inversion converges most rapidly for weak damping and in first order reduces to the adiabaticfollowing approximation for zero damping. Comparison with a numerical solution of the density-matrix equations demonstrates that retention of higher-order derivatives by retaining more terms in the integrated Fourier series results in a significant improvement in accuracy for a time-dependent field. However, Rabi oscillations, which can be an important component of the nonlinear response, ' cannot be obtained by finite-order approxi

Quantum Electronics and Laser Science Conference•Apr 24, 1989•0 citations

Transient polarization characteristics associated with resonant light scattering by a degenerate two-level system

Masaki Aihara

The time dependence of the polarized and unpolarized parts of the secondary-radiation intensity, ${I}_{\mathrm{p}}$ and ${I}_{\mathrm{up}}$, and the time-dependent degree of polarization $P(t)$ are calculated. A method is based upon the Feynman-diagram technique which is extended to enable the calculation of the reduced-density matrix element. It is found that, for the step excitation, $P(t)$ does not show a simple exponential decay. Especially in the off-resonance case, $P(t)$ shows oscillatory behavior. In the case of pulse excitation, ${I}_{\mathrm{up}}$ takes its maximum value after the excitation pulse has been turned off, and in the long-time limit it decays with the radiative-damping rate. In contrast to ${I}_{\mathrm{up}}$, ${I}_{\mathrm{p}}$ quickly responds to the excitation pulse and shows rapid exponential decay with the transverse-relaxation rate. It is also found that, contrary to the nondegenerate case, the total intensity of the secondary radiation does not exhibit simple exponential decay with the radiative-damping rate.

Physical Review A•Aug 1, 1978•20 citations

Optical coherent transients: Generalized two-level solutions

Axel Schenzle, Richard G. Brewer

Pulsed radio-frequency transient phenomena of resonant two-level spins are often described by the Bloch equations where the phenomenological dipole (${T}_{2}$) and population (${T}_{1}$) decay times appear. A similar description has been carried over recently into the microwave and infrared regions as well. At visible-ultraviolet wavelengths, however, a more general solution is required because the transition levels do not constitute an isolated system and also spontaneous emission can no longer be disregarded. As a consequence, the lower (1) and upper (2) transition levels are characterized by vastly different decay rates (${\ensuremath{\gamma}}_{1}\ensuremath{\ne}{\ensuremath{\gamma}}_{2}$) and the lifetime ${T}_{1}$ is no longer a defined quantity. In view of laser-induced coherent transient phenomena now being observed, we present general solutions of the density-matrix equations of motion for an optically resonant two-level system in contact with a multilevel reservoir where ${\ensuremath{\gamma}}_{1}\ensuremath{\ne}{\ensuremath{\gamma}}_{2}$. Particular and perturbative solutions involving intensity-dependent damping and frequency shifts are included with applications to echo and nutation phenomena.

Physical Review A•Nov 1, 1976•83 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

Adiabatic excitation of multilevel systems

G. L. Peterson, Cyrus D. Cantrell

We derive and discuss the behavior of several multilevel models relevant to the problem of molecular multiphoton excitation under the action of a quasimonochromatic laser field, the amplitude of which rises slowly and falls rapidly on a certain internal time scale. Our method is suitable for the numerical solution of problems involving the adiabatic excitation of large multilevel systems. We find adiabatic inversion to be a general characteristic of the models considered. The fact that adiabatic inversion is accompanied by marked qualitative and quantitative changes in the expectation value of the dipole operator makes our results useful for understanding the propagation of laser pulses. An explicit expression that is especially convenient for numerical evaluation is presented for determining the limiting internal time scale, which is determined by the energy levels and transition moments of the system and the laser frequency. We illustrate the criterion for adiabatic inversion by applying it to a particular multilevel model. From a density-matrix treatment we conclude that collisional damping destroys the coherence required for adiabatic inversion. On the basis of our results we believe that the adiabatic approximation may be appropriate for some experiments of current interest.

Physical Review A•Feb 1, 1985•30 citations

Phase-resolved time-domain nonlinear optical signals

Sarah M. Gallagher Faeder, David M. Jonas

A systematic theoretical and computational investigation of the microscopic factors which determine the phase of the signal field in time-resolved quasidegenerate three-pulse scattering experiments is presented. The third-order phase-matched response is obtained by density-matrix perturbation theory using a Green-function formalism for a system composed of two well-separated sets of closely spaced energy levels. Equations for calculating the electric field of four-wave-mixing signals generated by path-length delayed pulses are given. It is found that the phase of the signal field is determined by the excitation pulse phases, the dynamics of the nonlinear polarization decay, the product of four transition dipole matrix elements, and by a pulse-delay-dependent phase modulation at the frequency of the first dipole oscillation in the four-wave-mixing process. Analytic results for a two-level Bloch model show the phase shift from rapid nonlinear polarization decay. The product of dipole matrix elements is real and positive for three-level processes (bleached ground-state absorption and excited-state emission), but can be real and negative for some four-level Raman processes. The pulse-delay-dependent phase modulation treated here is closely related to the interferometric pulse-delay-dependent amplitude modulation observed in some collinear experiments, and plays a role in producing photon echos in inhomogeneously broadened samples. Numerical calculations of phase-resolved electric fields for finite duration pulses using a Brownian oscillator model appropriate for condensed-phase dynamics are presented. The ability of pulse-delay-dependent phase modulation to encode the frequency of the initially excited dipole onto the phase of the signal field can be exploited to examine energy-level connectivity, reveal correlations hidden under the inhomogeneous lineshape, and probe relaxation pathways in multilevel systems.

Physical Review A•Aug 18, 2000•41 citations
For the interaction of coherent radiation with a nondegenerate, homogeneously broadened multilevel system the nonlinear susceptibilities and excitations and their time dependence are evaluated using a perturbation approach to the density matrix up to fourth order. This is done for the "off-resonance" case, (\ensuremathΔ\ensuremathω2+\ensuremathΓ2)12\ensuremath≫1\ensuremathτ{(\ensuremath{\Delta}{\ensuremath{\omega}}^{2}+{\ensuremath{\Gamma}}^{2})}^{\frac{1}{2}}\ensuremath{\gg}\frac{1}{\ensuremath{\tau}}(\ensuremathΔ\ensuremathω2+\ensuremathΓ2)21​\ensuremath≫\ensuremathτ1​, as well as for the "on-resonance" case, (\ensuremathΔ\ensuremathω2+\ensuremathΓ2)12\ensuremath≪1\ensuremathτ{(\ensuremath{\Delta}{\ensuremath{\omega}}^{2}+{\ensuremath{\Gamma}}^{2})}^{\frac{1}{2}}\ensuremath{\ll}\frac{1}{\ensuremath{\tau}}(\ensuremathΔ\ensuremathω2+\ensuremathΓ2)21​\ensuremath≪\ensuremathτ1​, with \ensuremathτ\ensuremath{\tau}\ensuremathτ being the pulse duration, by taking into account the temporal change of the electric field amplitude and the phase and energy relaxation times, T2{T}_{2}T2​ and T1{T}_{1}T1​. It is shown that the transient one-photon excitation known as adiabatic following, is not only important for a pulse duration \ensuremathτ<T2\ensuremath{\tau}<{T}_{2}\ensuremathτ<T2​, but may also be important for T2\ensuremath≲\ensuremathτ<T1{T}_{2}\ensuremath{\lesssim}\ensuremath{\tau}<{T}_{1}T2​\ensuremath≲\ensuremathτ<T1​. In addition, expressions for the transient two-photon excitation in a multilevel system are derived. Finally, the two-level limit is evaluated to all orders of the perturbation including damping effects, which gives results valid also for large incident light intensities. For the "off-resonance" case an expression is obtained which describes not only the transient adiabatic following process, but also the usual one-photon absorption including saturation effects giving rise to power broadening.
We derive two types of perturbative approximation for the response of collisionally broadened two-level atoms to a time-dependent field. The perturbative approximations for the inversion and polarization are valid to any desired order in the derivatives of the field. Comparison with a numerical solution of the full density-matrix equations of motion demonstrates that the new approximations significantly improve upon the adiabatic-following and collisional steady-state approximations for time-dependent pulses in the appropriate damping regime. Application of the approximations to pulse propagation, sideband instability, and population pulsations is discussed. The collisional steady-state and adiabatic-following approximations find extensive use in many areas of nonlinear optics. The deficiencies of these approximations are well known: The collisional steady-state approxirnation is valid only for monochromatic (cw) radiation and the adiabatic-following approximation is valid only in the absence of damping and only if the product of the detuning and a characteristic time is much larger than one. In spite of their manifest faults, the collisional steady-state and adiabatic-following approximations are often exceptionally good approximations, in the appropriate damping regime, for the local response of a vapor of two-level atoms to a slowly varying field. However, timedependent nonlinear processes modify a propagating field zuch that it varies less slowly and the neglected portions of the atomic response become more important. ' Therefore improved approximations for the damped atomic response to a time-dependent field are desirable. In this paper we present two improved perturbative approximations for the inversion and polarization of collisionally broadened two-level atoms induced by a time-dependent field. The adiabatic-following approximation for two-level systems was derived phenomenologically by Grischkowsky, who used it to explain certain experimental results of near-resonance laser-pulse propagation. Motivated by the work of Grischkowsky, Crisp provided an analytic derivation for the adiabatic-following approximation and obtained a quadrature for the error. The essential feature of Crisp's derivation is the expansion of the kernel of the integral representation of the off-diagonal density-matrix element p2, in a Fourier series with integration proceeding term by term. Both the adiabatic-following approximation and the collisional steady-state approximation were obtained in first order. Both Crisp and Grischkowsky used the adiabatic inversion to derive a correction to the adiabatic-following approximation for the polarization. We have shown how Crisp's method can be extended to obtain corrections to the adiabatic inversion and to improve corrections to the polarization in the absence of damping. In this paper we show how the integrated Fourier series can be used to obtain two types of perturbative approximation for the inversion and polarization of collisionally broadened twolevel systems that are valid to any desired order in the derivatives of the field. The two types of perturbative approximation to the inversion are (i) a formal solution of the rate equation for the inversion which is obtained from the density-matrix equations for the populations; and (ii) an approximation obtained using the approximation (p22 — p») —I — 4p2, P, 2, valid for weak damping. The approximations for the polarization are obtained by substituting either (i) or (ii) for the inversion in the integrated-Fourier-series expansion of p2, . The approximation obtained using the first type of approximate inversion converges most rapidly for strong damping and reduces to the collisional steady-state approximation when the magnitude and phase of the field are constant. The approximation obtained using the second type of approxirnate inversion converges most rapidly for weak damping and in first order reduces to the adiabaticfollowing approximation for zero damping. Comparison with a numerical solution of the density-matrix equations demonstrates that retention of higher-order derivatives by retaining more terms in the integrated Fourier series results in a significant improvement in accuracy for a time-dependent field. However, Rabi oscillations, which can be an important component of the nonlinear response, ' cannot be obtained by finite-order approxi
The time dependence of the polarized and unpolarized parts of the secondary-radiation intensity, Ip{I}_{\mathrm{p}}Ip​ and Iup{I}_{\mathrm{up}}Iup​, and the time-dependent degree of polarization P(t)P(t)P(t) are calculated. A method is based upon the Feynman-diagram technique which is extended to enable the calculation of the reduced-density matrix element. It is found that, for the step excitation, P(t)P(t)P(t) does not show a simple exponential decay. Especially in the off-resonance case, P(t)P(t)P(t) shows oscillatory behavior. In the case of pulse excitation, Iup{I}_{\mathrm{up}}Iup​ takes its maximum value after the excitation pulse has been turned off, and in the long-time limit it decays with the radiative-damping rate. In contrast to Iup{I}_{\mathrm{up}}Iup​, Ip{I}_{\mathrm{p}}Ip​ quickly responds to the excitation pulse and shows rapid exponential decay with the transverse-relaxation rate. It is also found that, contrary to the nondegenerate case, the total intensity of the secondary radiation does not exhibit simple exponential decay with the radiative-damping rate.
Pulsed radio-frequency transient phenomena of resonant two-level spins are often described by the Bloch equations where the phenomenological dipole (T2{T}_{2}T2​) and population (T1{T}_{1}T1​) decay times appear. A similar description has been carried over recently into the microwave and infrared regions as well. At visible-ultraviolet wavelengths, however, a more general solution is required because the transition levels do not constitute an isolated system and also spontaneous emission can no longer be disregarded. As a consequence, the lower (1) and upper (2) transition levels are characterized by vastly different decay rates (${\ensuremath{\gamma}}_{1}\ensuremath{\ne}{\ensuremath{\gamma}}_{2}$) and the lifetime T1{T}_{1}T1​ is no longer a defined quantity. In view of laser-induced coherent transient phenomena now being observed, we present general solutions of the density-matrix equations of motion for an optically resonant two-level system in contact with a multilevel reservoir where ${\ensuremath{\gamma}}_{1}\ensuremath{\ne}{\ensuremath{\gamma}}_{2}$. Particular and perturbative solutions involving intensity-dependent damping and frequency shifts are included with applications to echo and nutation phenomena.
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}$.
We derive and discuss the behavior of several multilevel models relevant to the problem of molecular multiphoton excitation under the action of a quasimonochromatic laser field, the amplitude of which rises slowly and falls rapidly on a certain internal time scale. Our method is suitable for the numerical solution of problems involving the adiabatic excitation of large multilevel systems. We find adiabatic inversion to be a general characteristic of the models considered. The fact that adiabatic inversion is accompanied by marked qualitative and quantitative changes in the expectation value of the dipole operator makes our results useful for understanding the propagation of laser pulses. An explicit expression that is especially convenient for numerical evaluation is presented for determining the limiting internal time scale, which is determined by the energy levels and transition moments of the system and the laser frequency. We illustrate the criterion for adiabatic inversion by applying it to a particular multilevel model. From a density-matrix treatment we conclude that collisional damping destroys the coherence required for adiabatic inversion. On the basis of our results we believe that the adiabatic approximation may be appropriate for some experiments of current interest.
A systematic theoretical and computational investigation of the microscopic factors which determine the phase of the signal field in time-resolved quasidegenerate three-pulse scattering experiments is presented. The third-order phase-matched response is obtained by density-matrix perturbation theory using a Green-function formalism for a system composed of two well-separated sets of closely spaced energy levels. Equations for calculating the electric field of four-wave-mixing signals generated by path-length delayed pulses are given. It is found that the phase of the signal field is determined by the excitation pulse phases, the dynamics of the nonlinear polarization decay, the product of four transition dipole matrix elements, and by a pulse-delay-dependent phase modulation at the frequency of the first dipole oscillation in the four-wave-mixing process. Analytic results for a two-level Bloch model show the phase shift from rapid nonlinear polarization decay. The product of dipole matrix elements is real and positive for three-level processes (bleached ground-state absorption and excited-state emission), but can be real and negative for some four-level Raman processes. The pulse-delay-dependent phase modulation treated here is closely related to the interferometric pulse-delay-dependent amplitude modulation observed in some collinear experiments, and plays a role in producing photon echos in inhomogeneously broadened samples. Numerical calculations of phase-resolved electric fields for finite duration pulses using a Brownian oscillator model appropriate for condensed-phase dynamics are presented. The ability of pulse-delay-dependent phase modulation to encode the frequency of the initially excited dipole onto the phase of the signal field can be exploited to examine energy-level connectivity, reveal correlations hidden under the inhomogeneous lineshape, and probe relaxation pathways in multilevel systems.