Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Numerical experiments with optical pulses in V-type media

Open Access Open Access

Abstract

We present the results of a series of numerical experiments on the propagation of pairs of coupled short optical pulses in media with in-homogeneous broadening. The main results apply to propagation in three-level media in the V configuration.

©1999 Optical Society of America

Introduction

We have previously described a method for the analytic solution of the Maxwell-Bloch equations that govern the shape-preserving (self-similar) propagation of coupled pairs of short optical pulses in both lambda and V media, including the influence of inho-mogeneous (e.g., Doppler) broadening [1]. Our self-similar solutions were found to have the typical soliton-like sech shape. A restriction to matched pulses may be considered severe, and the sech shape is very specific, so the generality of the results was not easy to judge.

In this paper we present results under more general conditions. We relax both restrictions and show results for the propagation of non-sech pulses and pulses that are not matched, while restricting attention to V media. We present numerical results in five short movies. The movies demonstrate the anticipated result that unmatched pulses have the ability to match themselves to each other. They also show that the sech solution shape is not mathematically singular but the natural steady-state form for the V-medium pulses. Finally, and perhaps most interesting, we demonstrate the quantitative validity of the two-pulse Area Theorem presented recently [2], and show that it is obeyed even under conditions of pulse breakup.

 figure: Figure 1.

Figure 1. Schematic diagram of a V-type three-level atom and a two-level atom. The pulses Ωa and Ωb (driving transitions 1–2 and 2–3 respectively) are in two-photon resonance with intermediate detuning Δ in the V-type medium. The pulse Ω drives the 1–2 transition in the two-level medium.

Download Full Size | PDF

A two-level medium and a V-type three-level medium are illustrated in Fig. 1. The electric field vector for the two optical pulses in a V-type medium can be written as:

E=x̂εa(z,t)ei(kazωat)+x̂εb(z,t)ei(kbzωbt)+c.c.,

where kac = ωa and kbc = ωb ; and εa and εb , are the amplitudes of the electric fields of the two pulses. The Rabi frequencies that correspond to these fields are given by: Ωa (≡ 2d 12 εa /ħ) and Ωb (≡ 2d 23 εb /ħ), where dij is the dipole moment between levels i and j. In the Rotating Wave Approximation, the evolution equations for the atomic levels with complex amplitudes C 1, C 2, and C 3 can be obtained for V-type media from Schrodinger’s equation as:

iτC1=ΔC112ΩaC2
iτC2=12Ωa*C112Ωb*C3
iτC3=ΔC312ΩbC2

where Δ is the detuning shown in Fig. 1. The evolution equations for the fields can be written from the Maxwell equations as:

ζΩa=C2*C1g(Δ)dΔ
ζΩb=C2*C3g(Δ)dΔ

where g(Δ) is the distribution of detunings resulting from inhomogeneous (e.g. Doppler) broadening. Let T2* stand for the inhomogeneous broadening lifetime. The standard propagation coefficient is assumed to be equal for all transitions and is given by: μ = 4πd 2 /ħc, where N is the density of atoms in the medium. The equations (2–6) are written in local-time coordinates ζ and τ in the frame propagating with velocity c in the medium: ct - z and ζz. The medium is in the ground state so, for all depths ζ of the medium, at the initial time τ = 0 we have,

C1(ζ,0)=0,C2(ζ,0)=1,C3(ζ,0)=0.

The propagation equation of an optical pulse and the atomic response equations in a two level medium can be deduced from the V-type medium equations by ignoring level 3 and pulse ‘b’. Namely, the propagation equation of a pulse with Rabi frequency Ω in a two-level medium can be found from equations (2–6) by putting C 3 = 0 and Ωb = 0 and assuming Ω ≡ Ωa. In the following section we quickly review previous studies in two-level media.

Two level review

The propagation problem in a two-level medium has been studied extensively in the past [3]. McCall and Hahn [4] found that a pulse in an inhomogeneously broadened two-level medium obeys an “area theorem” [5]. If Ω is the Rabi frequency of the optical pulse then the area θ is defined by:

θ(ζ)=Ω(ζ,τ).

The area θ of evolves as function of the propagation distance ζ according to:

ζθ(ζ)=α2sinθ(ζ),

where α = μπg(0). Equation (9) is the area theorem. It predicts that pulses whose areas are multiples of 2π (i.e., 0, 2π, 4π, ⋯) are stable. Input pulses of other areas will approach one of these values during propagation. Moreover, we know that in an inhomogeneously broadened two-level medium, a 2π pulse gives rise to a sech-type soliton during propagation [4–7]:

Ω(ζ,τ)2τpsech(ττp),

where K is a propagation constant and τp is the pulse width. The area theorem and the formation of sech-type solitons in two-level media are illustrated in the following movie linked from Fig. 2.

The movie in Fig. 2 shows the evolution of a pulse as it propagates inside the medium. The pulse at the entry surface of the medium ζ = 0 has a super-Gaussian (“smooth square”) shape and an area of 2.5π. As it propagates inside, the shape changes and it becomes a sech-shaped pulse. Simultaneously the pulse area becomes 2π.

 figure: Figure 2.

Figure 2. (312KB) Snapshot of a frame of a movie showing the changing shape of an initially nearly square pulse as it propagates in a two-level medium. The pulse shape (≡ ΩT2*) is plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*)-The inset shows the pulse area as a function of space Z.

Download Full Size | PDF

Matched pulses in V-type media

Now consider a V-type three-level system and send two pulses (Ωa and Ωb,) inside the medium. We look at the propagation of a special kind of input pair called ‘matched’ pulses [8]. Two input pulses are called matched if they have proportional time envelopes when they enter the medium. So at the entry surface ζ = 0 the matched pulses are of the form:

Ωb(0,τ)=rΩa(0,τ),

where r is a proportionality constant. By direct substitution, we can show that the following relations are satisfied by Eqns. (2–6) and are consistent with the initial and boundary conditions in Eqns. (7–11):

Ωb(ζ,τ)=rΩa(ζ,τ)
C3(ζ,τ)=rC1(ζ,τ).

Using Eqn. (13) we can reduce Eqns. (2–6) and show that they are equivalent to the Maxwell-Schrödinger equations for two-level media [2]. This extends in a significant way an earlier result [9] showing the correspondence of the two-level and matched three-level Schrödinger equations. Thus a V-type medium with two matched pulses becomes equivalent to a two-level medium with a single pulse. This extended correspondence lets us translate the results for propagation from the two-level medium to the V-type medium. For example, the area theorem in Eqn. (9) can be translated to [2]:

ζθa2+θb2=α2sinθa2+θb2,

where θa and θb are the areas of pulse a and b respectively. Eqn. (14) shows that, for two matched pulses in V-type media, the quantity θa2+θb2 will become a multiple of 2π during propagation. Also by extending Eqn. (10) we find that pulses in a V-type medium will approach a soliton pair of the following type during propagation:

Ωa(ζ,τ)2τp1+r2sech(ττp)
Ωb(ζ,τ)2rτp1+r2sech(ττp).

The movies linked from Fig. 3 and Fig. 4 illustrate the propagation of two matched pulses inside V-type media.

 figure: Figure 3.

Figure 3. (398KB) Snapshot of a frame of a movie showing changing pulse shape during propagation in a V-type medium. The pulse shapes (≡ Ωa T2*, Ωb T2*) are plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*). The inset shows the pulse areas as a function of space Z.

Download Full Size | PDF

The movie in Fig. 3 shows two pulses that have super-Gaussian temporal shapes at the entry surface of the medium ζ = 0. As the pulses travel together, they change their shape and become a pair of sech pulses (checked by curve fitting). Also it is interesting to observe the pulse areas as they propagate. At the entry surface the pulses have θa = 2.24π and θb = 1.12π so that θa2+θb2=2.5π. As the pulses travel inside the medium, we observe that θa2+θb2 approaches 2π.

 figure: Figure 4.

Figure 4. (400KB) Snapshot of a frame of a movie showing changing pulse shape during propagation in a V-type medium. The pulse shapes (≡ Ωa T2*, Ωb T2*) are plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*). The inset shows the pulse areas as a function of space Z.

Download Full Size | PDF

The movie in Fig. 4 shows another example of matched pulses in a V-type medium. In this case, the pulses start as super-Gaussians with input areas of θa = 3.12π and θb = 1.56π so that θa2+θb2=3.5π. As the pulses travel together, they change their shape and break up into two pairs of sech pulses (checked by curve fitting). We also observe that θa2+θb2 approaches 4π during propagation, which is in accordance with Eqn. (14).

In both of the simulations of Fig. 3 and Fig. 4, if we measure the amplitude of each output sech pulse and denote the amplitudes of pulse a and b by A and B respectively then we observe that they satisfy an amplitude relation given by:

A2+B2=4τp2

where τp is the pulse width. This is consistent with the predictions in Eqn. (15–16), and reproduces a result first found for simulton propagation [10].

Arbitrary pulses in V-type media

So far we have discussed the special case of matched pulses in V-type media and showed the correspondence with two-level media. Now we will show the propagation of arbitrary pulses in V-type media. Fig. 5 is linked with a movie that shows the propagation of two unmatched pulses. The input pulses have Gaussian temporal shape. They start at different times and overlap for some time but they have their peaks at different times. As they propagate inside the medium, the pulses break into three secft-shaped pulse pairs. The inset of the movie shows that the pulse areas θa , θb do not follow the three-level area theorem as in Fig. 3 or Fig. 4 because the quantity θa2+θb2 does not approach a multiple of 2π. Nevertheless, if we measure the amplitudes of each sech pulse pair we find that Eqn. (17) is satisfied by all the pairs separately.

 figure: Figure 5.

Figure 5. (522KB) Snapshot of a frame of a movie showing changing pulse shape during propagation in a V-type medium. The pulse shapes (≡ Ωa T2*, Ωb T2*) are plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*). The inset shows the pulse areas as a function of space Z.

Download Full Size | PDF

A last example of two-pulse propagation in a V-type medium is shown in the movie linked from Fig. 6. We start with two sech-shaped pulses that do not overlap in time. These pulses have 2π area but different widths. As the pulses propagate in the medium, the faster pulse overtakes the slower pulse and the slower pulse gets a 180° phase shift.

 figure: Figure 6.

Figure 6. (462KB) Snapshot of a frame of a movie showing changing pulse shape during a pulse collision in a V-type medium. The pulse shapes (≡ Ωa T2*, Ωb T2*) are plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*). The inset shows the pulse areas as a function of space Z.

Download Full Size | PDF

Summary

We have shown the results of numerical experiments on propagation of different types of input pulse inside inhomogeneously broadened V-type media. When the input pulses are matched, the propagation dynamics of V-type three-level media is equivalent to that of two-level media. We have also shown that the input pulses become sech shaped pulse pairs during propagation, and obey a new area theorem for pulses in three-level media. When the pulses are not matched, the area theorem is not obeyed globally, but is still obeyed by the individual pairs of output pulses, which are still sech shaped and still obey the amplitude relation A 2 + B 2 = 4/τp2. Finally, we have shown a collision of two sech shaped pulses in a V-type medium.

This research was partially supported by NSF through grants PHY94-15583 and PHY97-22024.

References and links

1. A. Rahman and J. H. Eberly, “Theory of shape-preserving short pulses in inhomogeneously broadened three-level media,” Phys. Rev. A 58, R805 (1998). [CrossRef]  

2. A. Rahman, “Optical pulse propagation in V-type media,” submitted for publication.

3. L. Allen and J. H. Eberly, Optical Resonance and Two-Level Atoms (Dover Publications, Inc., New York,1987).

4. S. L. McCall and E. L. Hahn, “Self-induced transparency,” Phys. Rev. 183, 457 (1969). [CrossRef]  

5. For a recent new derivation of the Area Theorem in two-level media, see J. H. Eberly, “Area Theorem rederived,” Opt. Express 2, 173 (1998);http://epubs.osa.org/oearchive/source/4295.htm [CrossRef]   [PubMed]  

6. G. L. Lamb Jr., “Analytical descriptions of ultrashort optical pulse propagation in a resonant medium,” Rev. Mod. Phys. 43, 99 (1971). [CrossRef]  

7. G. L. Lamb Jr., “Coherent-optical-pulse propagation as an inverse problem,” Phys. Rev. A 9, 422 (1974). [CrossRef]  

8. R. J. Cook and B. W. Shore, “Coherent dynamics of N-level atoms and molecules. III. An analytically soluble periodic case,” Phys. Rev. A 20, 539 (1979). [CrossRef]  

9. C. E. Carroll and F. T. Hioe, “Analytic solutions for three-state systems with overlapping pulses,” Phys. Rev. A 42, 1522 (1990). [CrossRef]   [PubMed]  

10. M. J. Konopnicki and J. H. Eberly, “Simultaneous propagation of short different-wavelength optical pulses,” Phys. Rev. A 24, 2567 (1981). [CrossRef]  

Supplementary Material (5)

Media 1: MOV (312 KB)     
Media 2: MOV (398 KB)     
Media 3: MOV (400 KB)     
Media 4: MOV (522 KB)     
Media 5: MOV (461 KB)     

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Figure 1.
Figure 1. Schematic diagram of a V-type three-level atom and a two-level atom. The pulses Ω a and Ω b (driving transitions 1–2 and 2–3 respectively) are in two-photon resonance with intermediate detuning Δ in the V-type medium. The pulse Ω drives the 1–2 transition in the two-level medium.
Figure 2.
Figure 2. (312KB) Snapshot of a frame of a movie showing the changing shape of an initially nearly square pulse as it propagates in a two-level medium. The pulse shape (≡ ΩT2*) is plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*)-The inset shows the pulse area as a function of space Z.
Figure 3.
Figure 3. (398KB) Snapshot of a frame of a movie showing changing pulse shape during propagation in a V-type medium. The pulse shapes (≡ Ω a T2*, Ω b T2*) are plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*). The inset shows the pulse areas as a function of space Z.
Figure 4.
Figure 4. (400KB) Snapshot of a frame of a movie showing changing pulse shape during propagation in a V-type medium. The pulse shapes (≡ Ω a T2*, Ω b T2*) are plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*). The inset shows the pulse areas as a function of space Z.
Figure 5.
Figure 5. (522KB) Snapshot of a frame of a movie showing changing pulse shape during propagation in a V-type medium. The pulse shapes (≡ Ω a T2*, Ω b T2*) are plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*). The inset shows the pulse areas as a function of space Z.
Figure 6.
Figure 6. (462KB) Snapshot of a frame of a movie showing changing pulse shape during a pulse collision in a V-type medium. The pulse shapes (≡ Ω a T2*, Ω b T2*) are plotted as a function of time T (≡ τ/T2*) and space Z (≡ μζT2*). The inset shows the pulse areas as a function of space Z.

Equations (17)

Equations on this page are rendered with MathJax. Learn more.

E = x ̂ ε a ( z , t ) e i ( k a z ω a t ) + x ̂ ε b ( z , t ) e i ( k b z ω b t ) + c . c . ,
i τ C 1 = Δ C 1 1 2 Ω a C 2
i τ C 2 = 1 2 Ω a * C 1 1 2 Ω b * C 3
i τ C 3 = Δ C 3 1 2 Ω b C 2
ζ Ω a = C 2 * C 1 g ( Δ ) d Δ
ζ Ω b = C 2 * C 3 g ( Δ ) d Δ
C 1 ( ζ , 0 ) = 0 , C 2 ( ζ , 0 ) = 1 , C 3 ( ζ , 0 ) = 0 .
θ ( ζ ) = Ω ( ζ , τ ) .
ζ θ ( ζ ) = α 2 sin θ ( ζ ) ,
Ω ( ζ , τ ) 2 τ p sech ( τ τ p ) ,
Ω b ( 0 , τ ) = r Ω a ( 0 , τ ) ,
Ω b ( ζ , τ ) = r Ω a ( ζ , τ )
C 3 ( ζ , τ ) = r C 1 ( ζ , τ ) .
ζ θ a 2 + θ b 2 = α 2 sin θ a 2 + θ b 2 ,
Ω a ( ζ , τ ) 2 τ p 1 + r 2 sech ( τ τ p )
Ω b ( ζ , τ ) 2 r τ p 1 + r 2 sech ( τ τ p ) .
A 2 + B 2 = 4 τ p 2
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.