Trichromatic π-Pulse for Ultrafast Total Inversion of a Four-Level Ladder System

We present a numerical solution for complete population inversion in a four-level ladder system obtained by using a full π-pulse illumination scheme with resonant ultrashort phase-locked Gaussian laser pulses. We find that a set of pulse areas such as √ 3π, 2π, and √ 3π completely inverts the four-level system considering identical effective dipole coupling coefficients. The solution is consistent provided the involved electric fields are not too strong and it is amply accurate also in the case of diverse transition dipole moments. We study the effect of detuning and chirp of the laser pulses on the complete population inversion using the level structure of atomic sodium interacting with ps and fs pulses as an example. Our result opens the door for multiple applications such as efficient ultrashort pulse lasing in the UV or the engineering of quantum states for quantum computing.


Introduction
Control of population transfer in multiphoton transitions in atoms and molecules with laser pulses is a fundamental tool that has been extensively investigated in quantum optics since the advent of ultrashort pulses due to its value in several areas, including nonlinear spectroscopy [1], femtochemistry and biology [2], or in quantum information and in the quantum engineering of light states [3], where quantum systems need to be fully controlled.Several techniques have been reported for control of population transfer between a pair of quantum states involving multiphoton transitions, such as spectral shaping of the fields illuminating the medium [4,5], adiabatic methods [6][7][8][9] and π-pulse polychromatic control [10][11][12][13].
The coherent interaction between a three-level system and two fields has been studied theoretically for different processes [14][15][16][17][18][19][20][21].In a three-level ladder system with resonant fields the solution for π-pulse bichromatic control is analytical [17].The pulse area for complete transfer of the occupation probability from the ground state to the highest state is in this case √ 2π for both fields, with the pulse area defined as ∞ −∞ Ω(t)dt, where Ω(t) = µE(t)/h is the Rabi frequency of the transition, µ is the electric dipole moment and E(t) is the electric field amplitude of the laser pulse.The analysis of coupled differential equations arising from coherent multiphoton interactions however, even in simple multilevel quantum systems, is a complex task that often leads to expressions that lack any physical insight, so that very few have addressed the general problem analytically [22].It is therefore desirable to develop efficient numerical algorithms for the solution of the more complicated systems.
For the case of a four-level ladder system, population transfer schemes combining adiabatic and π-pulse techniques have been proposed [23].Full π-pulse polychromatic control techniques are however preferred in some cases because the intensities required for optimal excitation are on the order of the intensities for resonantly pumped single-photon transitions, and therefore they result energetically more efficient than adiabatic methods.In this paper we present a numerical solution for complete population inversion in a four-level ladder system by using full π-pulse control.

Numerical Simulations
We have solved numerically the Liouville equations for a four-level ladder system coherently interacting with a tricolor pulse without the assumptions of the rotating-wave approximation (RWA).Avoiding the RWA in the density-matrix formalism allows to account for the interaction of all the Fourier components of the pulses with all the permitted transitions in the system, and therefore in our study we can evaluate the effect of large detuning and chirp of ultrashort pulses with a broad range of intensities.The density matrix equations that we numerically solve are as follows [24].where ρ ii are the populations of the levels i, the off-diagonal elements ρ ij represent the coherences between levels i, j, µ ij are the dipole coupling coefficients of the corresponding electronic transitions, and ω ij = ( j − i )/h are the angular frequencies of the transitions i − j, with i being the energies of the corresponding quantum states.The decay rates of the transitions i−j are given by γ ij , and the population relaxation rates of the levels satisfy γ 22 = γ 21 , γ 33 = γ 31 + γ 32 , γ 44 = γ 41 + γ 42 + γ 43 .The overall rate γ ii , which is inversely proportional to the population decay time T 1 of level i (T 1 = γ −1 ii ), is thus given by the sum of each deactivation pathway in the four-level system.The decay rates of the coherences are defined by Γ12 = γ 22 /2 + Γ 12 , Γ13 = γ 33 /2 + Γ 13 , Γ14 = γ 44 /2 + Γ 14 , Γ23 = (γ 22 + γ 33 )/2 + Γ 23 , Γ24 = (γ 22 + γ 44 )/2 + Γ 24 and Γ34 = (γ 33 + γ 44 )/2 + Γ 34 , which include the relaxation of the coherences and the pure coherence dephasing rates Γ ij , which are defined as the inverse of the corresponding pure coherence dephasing times (Γ ij = T −1 2 ).The laser pulses used in our simulations are of Gaussian shape, with the electric field E(t) given by where E 01 , E 02 and E 03 are the peak amplitudes of the three pulse components E 1 (t), E 2 (t) and E 3 (t), and ω 1 , ω 2 and ω 3 are the respective optical angular frequencies.The duration of the pulses is given by τ = t p /( √ 2 ln 2), where t p is the full width at half maximum of the pulse intensity profile.In the present study we consider ultrashort pulses of ps and fs durations, so that considering an effective dipole coupling coefficient such as µ ∼ 10 −29 Cm for all transitions, the peak intensities for 2π pulses of 10 ps and 100 fs durations result in ∼ 10 5 W/cm 2 and ∼ 10 10 W/cm 2 , respectively.
The level structure that we have chosen in our simulations corresponds to the energy levels of sodium 23 11 Na, in particular our ground state represents the sodium 3s state, the intermediate states are the 3p and 4s states, and the upper state is the 7p state.We numerically integrate Equation ( 1) with an algorithm based on a 4-5 order Runge-Kutta solver which scans the areas of the three pulse components E 1 (t), E 2 (t) and E 3 (t) between 0 to 4π in steps of 0.01π.Starting with the system in the ground 3s state the goal is to obtain complete population transfer to the 7p state.We first consider a trichromatic pulse with a duration of 10 ps composed of three resonant components with central wavelengths of 589 nm, 1141 nm and 781 nm (see Figure 1).We obtain that the only solution for total population transfer from the ground state 3s to the highest state 7p is the set of pulse areas √ 3π, 2π, and √ 3π, as it is illustrated in Figure 1.As it is shown in Figure 1, the intermediate levels 3p and 4s are temporally occupied during the interaction process, but at the end of the interaction all the population is efficiently transferred to the 7p state.By varying the area of pulses E 1 (t) and E 3 (t) in a 4% we obtain that the population transfer to the 7p state is reduced to roughly 98%, which gives an estimation of the accuracy of our result.In the inset of Figure 1 we show how doubling the area of each pulse the 2π-pulse effect is retrieved, with all the population being back to the ground state at the end of the interaction.
. Time evolution of the populations of the four levels, as indicated, for a "π-pulse" transition produced by a tricolor resonant 10 ps Gaussian pulse of central angular frequencies ω 1 , ω 2 and ω 3 .The inset on the left shows a scheme of the four-level ladder system with the pulse areas that completely invert the population.Atomic sodium has been considered in our study as an example, with transition energies such as ω 12 = ω 1 = 16, 978 cm −1 (589 nm for 3s-3p), ω 23 = ω 2 = 8762 cm −1 (1141 nm for 3p-4s), and ω 34 = ω 3 = 12, 801 cm −1 (781 nm for 4s-7p).The inset on the right shows a complete oscillation of a "2π-pulse", which is produced with the corresponding double pulse areas 2 √ 3π, 4π, and 2 √ 3π.
The solution that we report can be scaled to other pulse durations, which is also described in the case of a three-level ladder system [12].Indeed, we find that complete population inversion is reproduced for pulses longer than 10 ps and also for shorter pulses.The limitation of the result is determined, on the one hand, on the relaxation times and line-broadening factors of the system, which prevent the usability of long pulses, and on the other hand, on the value of the transition dipole moments combined with the duration of the pulse, since weak coupling dipole moments with ultrashort pulses result in large peak intensities involving carrier field dynamics [26] and ac-Stark shifts that preclude the population inversion effect.As an example, for values of the transition dipole moments in sodium as µ 3s−3p = 2.11 × 10 −29 Cm, µ 3p−4s = 2.09 × 10 −29 Cm, and µ 4s−7p = 1.32 × 10 −30 Cm, and considering 10 ps laser pulses, the peak intensity of the corresponding laser pulses for √ 3π, 2π, and √ 3π illumination is I 1 = 4.3 × 10 5 W/cm 2 , I 2 = 5.8 × 10 5 W/cm 2 , I 3 = 1.1 × 10 8 W/cm 2 , which results in a population inversion to the 7p state of 99.3%.This population inversion is reduced to 98% for 5 ps pulses and to 84% for 1 ps pulses, which require pulse peak intensities of two more orders of magnitude.Dephasing times of the considered transitions in sodium are in general of the order of ns and therefore they do not influence the results obtained in the present study, which considers ps and fs pulses, but as commented above the complete density matrix equations including dephasing rates need to be computed for longer pulses.
Figure 2. Study of the effect of field detuning on the population inversion from level 3s to level 7p in atomic sodium (see Figure 1) for 100 fs pulses.Note that in (a-c) the most yellow regions correspond to 100% population transfer; while in (d-f), they correspond to 80%, 40% and 12%, respectively.
We have investigated the robustness of the full π-pulse trichromatic scheme against detuning and chirp of the laser pulses by considering different pulse durations.In order to study the effect of field detuning we define the following detuning parameters: with ω ij corresponding to the energy transition of the i, j levels (1 ≡ 3s, 2 ≡ 3p, 3 ≡ 4s, 4 ≡ 7p, see Figure 1).∆ 3 = 0 corresponds to all field detunings that keep the three-photon resonance between levels 7p and 3s, ∆ 2 = 0 means that the two-photon resonance between 3s and 4s is maintained, and ∆ 1 = 0 means that the E 1 field is in resonance with the 3s-3p transition.Figure 2 shows how complete population inversion from level 3s to level 7p is affected by detuning in the case of using 100 fs pulses.Figure 2a shows the results obtained when the three-photon resonance between 3s and 7p is kept (∆ 3 = 0).In the case shown in Figure 2b the two-photon resonance between 3s and 4s is maintained while the frequency of the three pulses is varied (∆ 2 = 0).In Figure 2c the results corresponding to having resonance in the lowest transition 3s-3p are shown (∆ 1 = 0).In these three first cases we observe a roughly symmetrical geometry around the center of the figure where the detuning is zero (∆ 1 = ∆ 2 = ∆ 3 = 0), which shows the sensitivity of the complete inversion effect.We observe that for inversion to stay at values larger than 95% the system supports detuning of the order of about 100 cm −1 , which roughly corresponds to the width of the Gaussian 100 fs pulses used (147 cm −1 ).In Figure 2d-f the three detuning parameters have been varied simultaneously.Figure 2d shows that for ∆ 3 = 75 cm −1 a comparable detuning ∆ 2 produces a population inversion as high as 80%.The combined effect of ∆ 2 and ∆ 3 is further shown in Figure 2e for a higher value of ∆ 3 = 150 cm −1 , where a corresponding higher value of ∆ 2 reduces the inversion to about 40%.For larger values of ∆ 3 population inversion is basically suppressed, as it is shown in Figure 2f, where the yellow color in the plot represents only a 12% of inversion.For longer pulses the supported detuning is scaled accordingly with the duration of the pulses, as it was also reported in the case of a three-level system interacting with two phase-locked pulses [12].This is, if a factor m is multiplied to the duration of the pulse, the same factor m divides the scale of the supported detuning, which concurs with the fact the detuning should be comparable to the spectral width of the pulses.The optimal pulse duration for a particular application and the corresponding effect of detuning will therefore be dependent basically on the line-broadening factors of the system under consideration and the value of the dipole moments of the transitions, which for ultrafast population inversion have to be high enough not to involve too strong laser fields, as commented above.
We next have considered pulse chirp on the three pulses interacting with the four-level system and have computed its effect for different transform limited (TL) pulse durations.In Figure 3 the effect of pulse chirp on the final populations is shown.The populations of the four levels after the pulse interaction are plotted as a function of the spectral chirp parameter ϕ , which has been implemented in the Fourier domain as ϕ(ω) = 0.5ϕ (ω − ω 0 ) 2 , with ϕ(ω) being the group delay dispersion (GDD) added to the spectral phase of the pulse.We note that, as only the spectral phase of the laser pulse is varied in our study, the frequency content remains as the TL pulse for any value of the spectral chirp ϕ .The rephasing of the spectral components produces a time broadening of the pulse given by t pc = t 4 p + 16(ln 2) 2 ϕ 2 /t p , where t pc is the full width at half maximum of the chirped pulse intensity profile [22], and therefore we implement the regular effect produced by for instance propagating the laser pulse in a linear dispersive medium.For GDD measurements of atmospheric and other gases with ultrashort pulses see, e.g., [27].:&-!-"#-=%-0;<"36###1 4%-!-"#-&% 4%-!-"-&%-0;<4#1 4%-!-"##-=%-0;<"5##1 4%-!-"#-=%-0;<"36###1 Effect of pulse chirp ϕ on the populations of the four levels after the interaction for different durations of the corresponding TL pulses.Indicated are the final populations of levels 7p (red lines) and 4s (blue lines); levels 3p and 3s become basically unpopulated in all cases (green and black lines).Note that the x-axis has been scaled for the different cases as indicated in the legend of the figure.Clearly, the population inversion produced by the longer pulses is less affected by chirp.In our study we have considered the same amount of spectral chirp for the three frequency components of the interacting pulse and with central frequencies resonant to their closest transition (ω 1 = ω 12 , ω 2 = ω 23 , ω 3 = ω 34 ).The red lines in Figure 3 correspond to the final population in level 7p, the blue lines show the final population in level 4s, and the green and black lines (close to zero) are the final populations in levels 3p and 3s.It is clear from Figure 3 that the main effect of pulse chirp in the case that we have considered is to bring part of the final population to level 4s, while levels 3p and 3s remain basically unpopulated.Note that the x-axis in Figure 3 have been scaled for the different pulse durations, as it is indicated in the legend of the figure.We observe that the trichromatic π-pulse effect is robust for spectral chirps of the order of ps 2 for TL pulses of 10 ps.In Table 1 we show the GDD (ϕ ) values supported for the pulses with TL durations of t p = 1 ps, 10 ps and 100 fs for a reduction in the population inversion to values of ≈ 92%, together with the associated pulse broadening (t pc ).Clearly, the shorter pulses support a lower value of chirp, with a pulse broadening that corresponds roughly to three times the TL duration for the three cases.We finally observe that the large electric fields associated to shorter pulses of the considered pulse areas prevent complete population inversion, as shown for the case of 10 fs pulses by the red-dotted line in Figure 3.

Conclusions
In conclusion, we have described the optimal conditions for ultrafast complete inversion of the electronic population of an atomic or molecular four-level ladder system using a full π-pulse control scheme.We have shown how much trichromatic π-pulse control is affected by detuning and pulse chirp.Other parameters that can affect the complete population transfer, such as fluctuations in the energy and duration of the pulses or the ionization rate associated to the shortest pulses, needs to be studied for every specific experimental geometry.Our study hence provides a general numerical solution to an intriguing fundamental theoretical problem which has many potential experimental applications to a variety of fields such as (i) spectroscopy and microscopy, (ii) quantum control of chemical and biochemical reactions, (iii) effective generation of UV ultrashort pulses and pumping of UV lasers, and (iv) quantum computing.