OSA's Digital Library

Optics Express

Optics Express

  • Editor: J. H. Eberly
  • Vol. 9, Iss. 10 — Nov. 5, 2001
  • pp: 506–518
« Show journal navigation

Role of group velocity in tracking field energy in linear dielectrics

M. Ware, S. A. Glasgow, and J. Peatross  »View Author Affiliations


Optics Express, Vol. 9, Issue 10, pp. 506-518 (2001)
http://dx.doi.org/10.1364/OE.9.000506


View Full Text Article

Acrobat PDF (517 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

A new context for the group delay function (valid for pulses of arbitrary bandwidth) is presented for electromagnetic pulses propagating in a uniform linear dielectric medium. The traditional formulation of group velocity is recovered by taking a narrowband limit of this generalized context. The arrival time of a light pulse at a point in space is defined using a time expectation integral over the Poynting vector. The delay between pulse arrival times at two distinct points consists of two parts: a spectral superposition of group delays and a delay due to spectral reshaping via absorption or amplification. The use of the new context is illustrated for pulses propagating both superluminally and subluminally. The inevitable transition to subluminal behavior for any initially superluminal pulse is also demonstrated.

© Optical Society of America

1 Introduction

We recently introduced [1

1. J. Peatross, S. A. Glasgow, and M. Ware, “Average Energy Flow of Optical Pulses in Dispersive Media,” Phys. Rev. Lett. 84, 2370–2373 (2000). [CrossRef] [PubMed]

] a context wherein the group delay function always retains meaning for electromagnetic pulses propagating in linear homogeneous media. This includes broadband pulses with arbitrarily complicated temporal structure and with spectra in regions of anomalous dispersion. The context we provide differs from traditional pedagogy [2

2. L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).

, 3

3. M. Born and E. Wolf, Principles of Optics, 7th Ed. (Cambridge, 1999), pp. 19–24.

, 4

4. J. D. Jackson, Classical Electrodynamics, 3rd Ed. (Wiley, New York, 1998), pp. 323, 330–335.

] wherein group velocity has been understood within a strict narrowband limit. (The new context has increasing experimental relevance since ultrashort laser pulses are routinely generated and characterized in the laboratory [5

5. R. Trebino, K. W. DeLong, D. N. Fittinghoff, J. N. Sweetser, M. A. Krumbugel, and B. A. Richman, “Measuring Ultrashort Laser Pulses in the Time-Frequency Domain Using Frequency-Resolved Optical Gating,” Rev. Sci. Instrum. 68, 3277–3295 (1997). [CrossRef]

].) In this article, we provide further details regarding this context, including a full derivation.

Before proceeding, we briefly review the solutions of Maxwell’s equation and how they relate to the traditional understanding of group delay [4

4. J. D. Jackson, Classical Electrodynamics, 3rd Ed. (Wiley, New York, 1998), pp. 323, 330–335.

]. In a linear, non-magnetic, non-conducting medium Maxwell’s equations can be written as

×E(r,t)+B(r,t)t=0,
(1)
×B(r,t)μ00E(r,t)t=P(r,t)t.
(2)

E(r,ω)=12πeiωtE(r,t)dt,
(3)
E(r,t)=12πeiωtE(r,ω)dω.
(4)

Analogous expressions for P and B give the time/frequency connection for the polarization and magnetic field. In our convention, we take all fields to be real in the time domain, so the following symmetry holds in the frequency domain:

E(r,ω)=E*(r,ω)
(5)

(with analogous expressions for B and P).

In the frequency domain, Maxwell’s equations have as a solution

E(r0+Δr,ω)=E(r0,ω)eik·Δr,
(6)

where the initial pulse form E(r 0, ω) is chosen at a point r 0. The solution renders the pulse form E(r 0r, ω) (in terms of the initial pulse form) at any other point shifted by the displacement Δr. (Implicit in (6) is the restricting assumption that each frequency is associated with a single wave vector k(ω), although the wave vectors need not be all in the same direction.) The magnetic field is obtained in terms of the electric field through

B(r,ω)=k(ω)×E(r,ω)ω.
(7)

In accordance with the linear, isotropic, and homogeneous assumptions, the susceptibility χ (ω) defines the (temporally non-local, but frequency local) constitutive relation between the electric field and the polarization:

P(r,ω)=0χ(ω)E(r,ω).
(8)

Maxwell’s equations together with (8) also impose the constraint

k2(ω)=ω2c2[1+χ(ω)].
(9)

Here we take the traditional vantage point of considering the frequency ω to be real, whereas the wave vector k (ω) can be complex when the susceptibility χ (ω) is complex.

In order to obtain the temporal profile of the pulse, the solution (6) must be transformed to the time domain via (4). The traditional concept of group velocity is obtained by expanding the frequency-dependent phase delay as

k(ω)·Δr[kω-·Δr]+[kωω¯·Δr](ωω-)+
(10)

before carrying out the inverse transform. After performing the inverse transform with this approximate phase delay, the real part of the coefficient of the linear frequency term is seen to give (to first order) the time required for the pulse to traverse the displacement Δr. The function Rek/∂ω·Δr is known as the group delay function, and in the traditional context is evaluated at a single ‘carrier’ frequency ω¯. Group velocity is obtained by dividing the displacement Δr by the group delay function evaluated at ω¯.

The failings of (10) are well documented [6

6. K. E. Oughstun and H. Xiao, “Failure of the Quasimonochromatic Approximation for Ultrashort Pulse Propagation in a Dispersive, Attenuative Medium,” Phys. Rev. Lett. 78, 642–645 (1997). [CrossRef]

]. In particular, if the bandwidth of an electromagnetic pulse encompasses a substantial portion of a resonance structure, it becomes necessary to retain a large number of terms in the expansion (10) to accurately describe the phase delay. Moreover, as was pointed out by Garret and McCumber [7

7. C. G. B. Garrett and D. E. McCumber, “Propagation of a Gaussian Light Pulse through an Anomalous Dispersion Medium,” Phys. Rev. A 1, 305–313 (1970). [CrossRef]

], if the bandwidth exceeds the modulus of the difference between the carrier frequency and the nearest pole in the refractive index, the series fails to converge altogether. Because of these difficulties, textbooks on the subject state that the traditional concept of group velocity “has no longer any appreciable physical significance” [3

3. M. Born and E. Wolf, Principles of Optics, 7th Ed. (Cambridge, 1999), pp. 19–24.

] or “may even lack precise meaning” [4

4. J. D. Jackson, Classical Electrodynamics, 3rd Ed. (Wiley, New York, 1998), pp. 323, 330–335.

] in these situations. In this article (see also [1

1. J. Peatross, S. A. Glasgow, and M. Ware, “Average Energy Flow of Optical Pulses in Dispersive Media,” Phys. Rev. Lett. 84, 2370–2373 (2000). [CrossRef] [PubMed]

]), we endeavor to supply this “precise meaning” that has been absent in the broadband context. In particular, we demonstrate that the time required for a pulse to traverse a displacement Δr is described by a linear superposition of group delays, weighted by the surviving spectral content of the pulse. Mathematically, this delay is given by

ΔtG(r)=ρ(r,ω)(Rekω·Δr)dω,
(11)

where ρ (r, ω) is a normalized spectral distribution of Poynting flux (at the final point) defined by

ρ(r,ω)η̂·S(r,ω)η̂·S(r,ω)dω.
(12)

(η^ is a unit vector normal to the detector surface.) Another delay term due to spectral reshaping also plays a role in determining the total delay time. This additional delay term, which is evaluated without propagation, becomes important when an initially chirped pulse undergoes an alteration in its spectral profile. For an unchirped pulse, the reshaping effect is generally negligable. No approximations or expansions of the type (10) are used in the new formulation, making the results general.

We note that the fact that group velocity can become superluminal and/or negative in the vicinity of a resonance structure does not indicate a failure of the expansion (10) [7

7. C. G. B. Garrett and D. E. McCumber, “Propagation of a Gaussian Light Pulse through an Anomalous Dispersion Medium,” Phys. Rev. A 1, 305–313 (1970). [CrossRef]

, 8

8. R. Y. Chiao, “Superluminal (but Causal) Propagation of Wave Packets in Transparent Media with Inverted Atomic Populations,” Phys. Rev. A. 48, R34–R37 (1993). [CrossRef] [PubMed]

, 9

9. E. L. Bolda, J. C. Garrison, and R. Y. Chiao, “Optical Pulse Propagation at Negative Group Velocities due to a Nearby Gain Line,” Phys. Rev. A 49, 2938–2947 (1994). [CrossRef] [PubMed]

, 10

10. R. Y. Chiao and A. M. Steinberg, “Tunneling Times and Superluminality,” Progress in Optics 37, pp. 347–406 (Emil Wolf ed., Elsevier, Amsterdam, 1997). [CrossRef]

]. This type of exotic behavior (well verified experimentally [11

11. S. Chu and S. Wong, “Linear Pulse Propagation in an Absorbing Medium,” Phys. Rev. Lett. 48, 738–741 (1982). [CrossRef]

, 12

12. L. J. Wang, A. Kuzmmich, and A. Dogariu, “Gain-Assisted Superluminal Light Propagation,” Nature 406, 277–279 (2000). [CrossRef] [PubMed]

]) occurs in a narrowband limit. The context we present describes the gradual (and inevitable) transition from superluminal to subluminal behavior as the initial temporal profile of a pulse becomes more sharply defined, and the expansion (10) cannot be used. Since the new formalism is exact, the extreme broadband limit can also be taken which leads to the well-known Sommerfeld-Brillouin [2

2. L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).

] result that a sharply defined signal travels at exactly c, the speed of light in vacuum.

2 Average arrival time of poynting flux to a point

In order to concretely discuss the time required for a pulse to traverse a given displacement, we first define the arrival time of a pulse to a point (say where a detector is located). This definition need only involve the Poynting vector, since it alone transports energy. To deal with arbitrary broadband pulses, the arrival time should avoid presupposing a specific pulse shape, since the pulse may evolve in complicated ways during propagation. For example, the peak or the midpoint on the rising edge of a pulse is a poor indicator of arrival time if the pulse contains multiple peaks or a long and non-uniform rise time.

For the reasons given, we use a time expectation integral (proposed by Smith in 1970 [13

13. R. L. Smith, “The Velocities of Light,” Am. J. Phys. 38, 978–983 (1970). [CrossRef]

]) to specify the arrival of the pulse:

trtρ(r,t)dt,
(13)

where ρ (r, t) is a normalized temporal distribution of the Poynting flux at r:

ρ(r,t)η̂·S(r,t)η̂·S(r,t)dt.
(14)

The Poynting vector is defined as as usual:

S(r,t)E(r,t)×B(r,t)μ0.
(15)

As Smith explained, this definition “is not restricted to a quasi monochromatic pulse or to a situation which one has always only a single pulse. … It yields a unique value regardless of the dispersion present.” The unit vector η^ refers to the direction in which the energy flow is detected (normal to a detector surface). It has importance for angularly dispersive systems such as grating pairs used for pulse compression where the result of the integral in the numerator is not necessarily parallel to that in the denominator. We have addressed the angularly dispersive case elsewhere [14

14. M. Ware, W. E. Dibble, S. A. Glasgow, and J. Peatross, “Energy Flow in Angularly Dispersive Optical Systems,” J. Opt. Soc. Am. B 18839–845 (2001). [CrossRef]

].

Motivated by a desire to make connection with the traditional group velocity, we rewrite the arrival time (13) of the pulse in terms of the frequency representations of the fields involved. The denominator of (14) may be transformed immediately via Parseval’s theorem:

S(r,t)dt=S(r,ω)dω,
(16)

where the frequency representation of the Poynting vector is given by:

S(r,ω)E(r,ω)×B*(r,ω)μ0,
(17)

Employing the Fourier relations (3) and (4) in the numerator of (13) (with (14) inserted), we find

tS(r,t)dt=t[12πE(r,ω)eiωtdω]×[12πB(r,ω')eiω'tdω']dt.
(18)

By introducing a frequency derivative to account for the factor t, the expression (18) can be recast as

tS(r,t)dt=idωE(r,ω)×ω[dωB(r,ω)μ012πei(ω+ω)tdt].
(19)

The final integral is recognized as a delta function, which permits the integration over ω′ to be performed:

tS(r,t)dt=iE(r,ω)×ωB(r,ω)μ0dω=iE(r,ω)ω×B*(r,ω)μ0dω.
(20)

Here we have invoked the symmetry relation (5). We have also integrated by parts in order to place (by preference) the derivative on the electric field rather than on the magnetic field. In doing this, we have assumed a pulse of a finite duration. Equations (16) and (20) lead to the frequency domain representation of the arrival time:

tr=T[E(r,ω)],
(21)

where

T[E(r,ω)]iη̂·E(r,ω)ω×B*(r,ω)μ0dωη̂·S(r,ω)dω.
(22)

In (22) we have indicated that the function T depends only on the electric field since the magnetic field can be written in terms of the electric field through (7).

3 Delay time between two points

The arrival time of the pulse as expressed in (22) is not very useful in itself. Its usefulness is revealed when applied to the difference in arrival times of a pulse to two different locations, say r 0 and r=r 0r. The difference in the arrival times of the two points is defined as

Δttrtr0=T[E(r0+Δr,ω)]T[E(r0,ω)].
(23)

Our goal in this section is to recast (23) into a form that makes clear the connection between Δt and the group delay function Rek/∂ω·Δr.

In order to arrive at the desired result, we first obtain an explicit form for T [E(r 0r, ω)]. With the solution (6) we can evaluate the integrand of the numerator of (22) at r 0r as follows:

E(r0,ω)eik·Δrω×B*(r0,ω)eik*·Δrμ0=iRek·ΔrωE(r0,ω)×B*(r0,ω)e2Imk·Δr+E(r0,ω)eImk·Δrω×B*(r0,ω)eImk·Δrμ0.
(24)

The integrand of the denominator in (22) is also evaluated at r 0r, giving

S(r0+Δr,ω)=E(r0,ω)eik·Δr×B*(r0,ω)eik*·Δrμ0
=E(r0,ω)eImk·Δr×B*(r0,ω)eImk*·Δrμ0.
(25)

Upon using (24) and (25) with the arrival time (22) we obtain

T[E(r0+Δr,ω)]=η̂·S(r,ω)(Rekω·Δr)dωη̂·S(r,ω)dω
iη̂·E(r0,ω)eImk·Δrω×B*(r0,ω)eImk*·Δrμ0dωη̂·e2Imk·ΔrS(r0,ω)dω.
(26)

Using the expressions (22) and (26), we can now write the delay time (23) as the sum of two intuitive terms:

Δt=ΔtG(r)+ΔtR(r0).
(27)

The first term in (27), the net group delay, is a spectral average of the group delay function over the spectral intensity that is experienced at the final point r (as already mentioned in connection with):

ΔtG(r)=ρ(r,ω)(Rekω·Δr)dω,
(28)

where we have defined a normalized spectral distribution of field energy:

ρ(r,ω)η̂·S(r,ω)η̂·S(r,ω)dω.
(29)

Note the close resemblance between (13) and (28). Both are expectation integrals, the former being executed as a ‘center of mass’ integral in time, the later being executed in the frequency domain on Rek/∂ω·Δr. Clearly, the group delay function evaluated at every frequency present in the pulse influences the result. The net group delay depends only on the spectral content of the pulse, independent of its temporal organization (i.e. the phases of E(r, ω) and B(r, ω) do not contribute). Only the real part of k plays a direct role in (28).

The second term in (27) represents a delay which arises solely from a reshaping of the spectrum through absorption (or amplification) and is given by

ΔtR(r0)T[eImk·ΔrE(r0,ω)]T[E(r0,ω)].
(30)

This reshaping delay is evaluated at r 0, before propagation takes place. The reshaping delay is the difference between the pulse arrival time at the initial point r 0 evaluated without and with the spectral amplitude that will be lost (or gained) during propagation. Both terms in (30) utilize the phase of the fields at r 0. Thus, in contrast to the net group delay, the reshaping delay is sensitive to how the pulse is organized. The reshaping delay is zero if the spectral amplitude of the pulse is unaltered during propagation (i.e. if the imaginary part of k can be neglected). The reshaping delay is also relatively unimportant if the pulse is unchirped before propagation and in the narrowband limit even if pulses experience strong absorption (or amplification).

Before proceeding, it is interesting to note that in (27) the ordering for the evaluation of the net group delay and reshaping delay can be reversed. In other words, the arguments r and r 0 in (27) can be interchanged which gives

ΔtG(r)+ΔtR(r0)=ΔtG(r0)+ΔtR(r).
(31)

This commutative property can dramatically alter the group and reshaping delays taken individually. Nevertheless, their sum Δt is unaffected by the ordering. The group delay computed over the initial spectrum is

ΔtG(r0)=ρ(r0,ω)(Rekω·Δr)dω.
(32)

When the group delay is computed using the initial spectrum before the propagation, the reshaping delay is evaluated at the end of propagation without and with the lost spectral amplitude:

ΔtR(r)=T[E(r0+Δr,ω)]T[eImk·ΔrE(r0+Δr,ω)].
(33)

In this expression, the initial spectral amplitude is restored to the second term. The proof of (31) is straightforward. Directly by substitution into (26) we have

T[eImk·rE(r0+Δr,ω)]=η̂·S(r0,ω)(Rekω·Δr)dωη̂·S(r0,ω)dω
iη̂·E(r0,ω)ω×B*(r0,ω)μ0dωη̂·S(r0,ω)dω.
(34)

When this and (26) are used in (33), the commutative relationship (31) is immediately verified.

Of course in the narrowband limit we expect the delay time to give the traditional result presented in section 1. In this limit the reshaping delay tends to zero so that the total delay is dominated by the net group delay. The normalized spectral distribution approaches a delta function (ρ(ω)→δ(ω-ω¯)) so that the total delay reduces to Δt→ΔtG =Rek/∂ω|ω¯·Δr. This is in agreement with the well verified observation [7

7. C. G. B. Garrett and D. E. McCumber, “Propagation of a Gaussian Light Pulse through an Anomalous Dispersion Medium,” Phys. Rev. A 1, 305–313 (1970). [CrossRef]

, 8

8. R. Y. Chiao, “Superluminal (but Causal) Propagation of Wave Packets in Transparent Media with Inverted Atomic Populations,” Phys. Rev. A. 48, R34–R37 (1993). [CrossRef] [PubMed]

, 9

9. E. L. Bolda, J. C. Garrison, and R. Y. Chiao, “Optical Pulse Propagation at Negative Group Velocities due to a Nearby Gain Line,” Phys. Rev. A 49, 2938–2947 (1994). [CrossRef] [PubMed]

, 10

10. R. Y. Chiao and A. M. Steinberg, “Tunneling Times and Superluminality,” Progress in Optics 37, pp. 347–406 (Emil Wolf ed., Elsevier, Amsterdam, 1997). [CrossRef]

, 11

11. S. Chu and S. Wong, “Linear Pulse Propagation in an Absorbing Medium,” Phys. Rev. Lett. 48, 738–741 (1982). [CrossRef]

] that a pulse travels at the group velocity evaluated at ω¯ in the narrowband limit.

4 The effect of pulse bandwidth on delay time

[Ren(ω)+iImn(ω)]2=1+χ(ω)
(35)
Fig. 1. (a) Real and imaginary parts of the refractive index. (b) Group delay function for the given index

For simplicity, we consider propagation in only one dimension, so that k=kẑ. We will consider the delay time of a pulse as it propagates through a displacement Δr=0.1 (c/γ) . Figure 1(b) depicts the group delay function ∂k/∂ω·Δr for these parameters.

Fig. 2. (a) Normalized spectral distribution ρ for the narrowband pulse centered on resonance before (dotted) and after (solid) propagation. (b) Total delay as a function of ω¯ for the narrowband pulse. (c) Overall pulse transmission as a function of ω¯

The electric field at r 0 is chosen to be Gaussian, E(r 0, t)= E 0exp(-t 2/τ 2)cos (ω¯ t), with initial duration τ=10/γ. In all of the following examples the pulses are initially unchirped, so that the reshaping delay ΔtR is negligable. Figure 2(a) depicts the normalized spectral distribution function (Eq. (29)) of the narrowband pulse before (dotted line) and after (solid line) propagation in the medium, with the initial spectrum centered on resonance (i.e. ω¯=ω 0). Because the initial pulse spectrum is narrow compared to the resonance, ρ changes little during propagation. Since the reshaping delay is negligable, the total delay time is given by the expectation of the group delay function (shown in Fig. 1(b)) weighted by the spectral distribution (shown in Fig. 2(a)) as in Eq. (28). Figure 2(b) plots the delay time as the pulse’s central frequency ω¯ is varied relative to the resonance frequency ω 0. Since ρ(r, ω) is narrow compared to the width of the resonance, the total delay closely resembles the group delay function in Fig. 2 (b). Figure 2 (c) plots the ratio of the pulse energy after propagation to that before propagation as the pulse’s central frequency is varied.

Fig. 3 repeats Fig. 2 for the case of the broadband pulse. Note that in the broadband case, the negative delay times seen in Fig. 2(b) disappear. This illustrates the general principle that superluminal effects are a narrowband phenomenon which disappears as bandwidth is increased.

The difference between the narrowband and broadband examples just discussed can be clearly seen when one constructs the spatial profile of the pulses’ field energy using the methods summarized in the introduction. Figure 4 is an animation of the narrowband (top) and broadband (bottom) pulses from Figs. 2 and 3 as they traverse the amplifying medium (represented by the vertical bar). For this illustration we have centered the bandwidth of the pulses at ω¯-ω 0=γ where the pulse in Fig. 2 exhibits the most dramatic superluminal behavior. The dotted line represents the pulse as it would have travelled in vacuum and the solid line indicates the actual pulse profile. In the animation it is clear that the centroid of the narrowband pulse is shifted forward relative to the vacuum pulse upon traversing the medium. In contrast, the broadband pulse is delayed relative to the vacuum pulse.

Fig. 3. (a) Normalized spectral distribution ρ for the broadband pulse centered on resonance before (dotted) and after(solid) propagation. (b) Total delay as a function of ω¯ for the broadband pulse. (c) Overall pulse transmission as a function of ω¯
Fig. 4. Animation of narroband (top) and broadband (bottom) gaussian pulses traversing an amplifying medium (1.4 MB)
Fig. 5. Total delay as the width of the pulse is changed from broadband to narrow band.

We now illustrate the transition from superluminal to subluminal delay times as the pulse bandwidth increases. For this example we use the medium described in relation to Fig. 1 and center the bandwidth of the pulse at ω¯-ω 0=γ. Figure 5 plots the delay time as the pulse duration τ (and hence the bandwidth) is varied. As expected, the pulses with narrow bandwidths (large τ) exhibit superluminal effects. The behavior becomes subluminal as the bandwidth increases (τ decreases). This is easily understood in terms of (28). As the bandwidth increases, ρ(r, ω) includes more and more of the ‘slow’ spectral components near ω 0. When the expectation (28) is taken, these ‘slow’ spectral components outweigh the effect of the frequencies with a negative group delay function resulting in a subluminal delay time.

As the pulse duration goes to zero in Fig. 5, the total delay time decreases and approaches an exactly luminal delay time. This behavior can be clearly explained using the context developed in sections 2–3. When the duration of the pulse goes to zero, becoming a delta function in time, the spectral content of the pulse extends uniformly over all frequencies (ρ(ω)→constant). By comparison, any resonance structure encompasses a finite bandwidth, and the reshaping delay tends to zero. If, for simplicity, we consider all wave vectors to be parallel (i.e. k(ω)=k(ω)) we have from (28) or (32) that Δt→ΔtG ≈Rek/ω|·Δr=·Δr/c, assuming a refractive index which approaches unity at high frequency. This demonstrates (without lengthy analyses or numerical simulations) the well-known fact that the velocity for a sharply defined signal is exactly c. Thus, the Sommerfeld-Brillouin result [2

2. L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).

] of luminality for pulses of definite support is consistent with and even implied by the present context.

Fig. 6. Animation of narroband (top) and broadband (bottom) gaussian pulses traversing an absorbing medium (1.7 MB)

As has been known for several decades [7

7. C. G. B. Garrett and D. E. McCumber, “Propagation of a Gaussian Light Pulse through an Anomalous Dispersion Medium,” Phys. Rev. A 1, 305–313 (1970). [CrossRef]

, 11

11. S. Chu and S. Wong, “Linear Pulse Propagation in an Absorbing Medium,” Phys. Rev. Lett. 48, 738–741 (1982). [CrossRef]

, 4

4. J. D. Jackson, Classical Electrodynamics, 3rd Ed. (Wiley, New York, 1998), pp. 323, 330–335.

], superluminal effects are also observed in absorbing dielectric media. In an absorbing medium, the group delay function is similar to that plotted in Fig. 2(b), except it is inverted (i.e. the group delay function is negative near resonance and becomes positive away from resonance). Thus, in an absorbing medium pulses whose spectra are centered on resonance exhibit the most dramatic superluminal delay times. Here again, if the pulse’s spectrum is too broad, the effects of the slow off-resonance spectral components in ρ(r, ω) will outweigh the on resonance components resulting in luminal delay times. Figure 6 is an animation illustrating the difference between narrowband and broadband pulses in an absorbing medium. The pulse parameters are the same as in fig. 4, and the medium parameters are the same except f=+1 in order to make an absorbing medium. (These parameters are identical to those used in [1

1. J. Peatross, S. A. Glasgow, and M. Ware, “Average Energy Flow of Optical Pulses in Dispersive Media,” Phys. Rev. Lett. 84, 2370–2373 (2000). [CrossRef] [PubMed]

].) Again it is the narrowband pulse that is seen to be advanced relative to the vacuum pulse.

5 Propagation distance and delay time

Although we have demonstrated that superluminal behavior is most pronounced for narrowband pulses, when Δr is small enough it is possible to observe superluminal effects for broadband pulses. This was demonstrated recently by Talukder and coworkers [16

16. Md. Aminul Islam Talukder, Yoshimitsu Amagishi, and Makoto Tomita, “Superluminal to Subluminal Transition in the Pulse Propagation in a Resonantly Absorbing Medium,” Phys. Rev. Lett. 86, 3546–3549 (2001). [CrossRef]

] when they performed a superluminal propagation experiment in an absorbing medium under relatively broadband conditions (where the traditional narrowband expansion of group velocity fails). Although their femtosecond laser pulses initially propagated superluminally through a dye cell, as the effective length of the propagation medium was increased, they observed a transition to subluminal behavior. The formalism developed in the present paper may be used to explain why this transition inevitably occurs as an initially superluminal pulse is allowed to propagate longer distances in an absorbing medium.

If the pulse is initially unchirped, the reshaping delay ΔtR (r 0) is negligable, so that the total delay is given by ΔtG (r). When the pulse spectrum is centered on-resonance, ρ(r, ω) will primarily include the ‘fast’ spectral components near ω 0 (remember that in an absorbing medium, frequencies near ω 0 have the shortest group delay function). However, as it travels, the fast spectral components near ω 0 are absorbed from the pulse spectrum. Eventually only spectral components away from the resonance (i.e. frequencies where the group delay function gives a luminal delay) are left in the spectrum. The delay time, given by the net group delay (28), then becomes subluminal.

The situation in an amplifying medium is similar. In this case, it is the off resonance frequency components that have the superluminal group delay. As the pulse propagates, the amplifying medium adds the slow near-resonance frequencies to the pulse spectrum. Thus, an initially superluminal pulse whose spectrum is centered off-resonance will eventually have a luminal delay time as ρ (r, ω) changes to include more near-resonance frequencies.

Fig. 7. (a) Total delay, scaled by cr, as a function of displacement. (b) The group delay function divided by displacement.

To illustrate this transition to luminal delay times as propagation distance is increased, Figure 7(a) plots delay time as a function of displacement in an amplifying medium (medium parameters are the same as in Fig. 1 except for the variable displacement). The initial pulse is chosen to be a narrowband gaussian as in fig. 2, with τ=10/γ and ω¯-ω 0=1γ. For ease of comparison, the delay time is scaled by cr so that propagation at the speed of light has a scaled delay time of 1. Figure 7 (b) shows the group delay function (also scaled by cr). For small displacements, the delay time matches the group delay function evaluated at ω¯. As the pulse propagates and acquires more on-resonance spectral intensity, the delay time shifts to the on-resonance value of the group delay function.

The preceding example illustrates the principle that superluminal effects in amplifying media will be observed over the longest distances when the pulse’s spectrum has little spectral content near the resonance. In absorbing media, as mentioned above, superluminal effects will persist longest for pulses with the least spectral content away from resonance. This is why narrowband pulses are generally used in both amplifying and absorbing experiments designed to observe these effects. Another technique was proposed by Bolda et. al. [9

9. E. L. Bolda, J. C. Garrison, and R. Y. Chiao, “Optical Pulse Propagation at Negative Group Velocities due to a Nearby Gain Line,” Phys. Rev. A 49, 2938–2947 (1994). [CrossRef] [PubMed]

] which allows for longer superluminal propagation distances in amplifying media. They proposed to initially send the pulse through an absorbing medium to remove the near-resonance components from the pulse spectrum. This modified pulse is then sent through an amplifier with the same resonance frequency as the absorber. Since the slow near-resonance spectral components were attenuated in the absorber, a relatively large propagation distance is required before they enter into the expectation (28). (Of course there is no free lunch here: the extra long delay time in the absorber counterbalances the short delay time in the amplifier.)

Fig. 8. Animation of a gaussian pulse traversing an absorbing medium followed by an amplifying medium as proposed by Bolda [9] (1.7 MB)

Figure 8 reproduces a time sequence of the spatial distribution of field energy which appeared in the article describing the experiment proposed by Bolda. Again, a dashed pulse is shown which propagates at exactly c (i.e. as if the absorber and amplifier were not there). The Lorentz model is used for the two media with ω 0=2×107γ, ωp =100γ, and Δr=0.5c/γ. The oscillator strength is f=1 for the absorbing medium and f=-1 for the amplifier. The pulse is chosen to be gaussian as above with τ=0.264 √2γ and ω¯-ω 0=ωp /3.

Fig. 9. (a) Spectral distribution before (dashed) and after (solid) traversing the absorber. The distribution after the amplifier is the same as the initial distribution. (b) Group delay function for the amplifying medium.

Figure 9(a) plots the spectral distribution of the pulse at different points in the medium. Figure 9(b) shows the group delay function for the amplifying medium. As usual, the delay time is predicted by the expectation of the group delay function weighted by the spectral distribution. (In the amplifying medium the reshaping delay ΔtR also plays a role since the pulse was chirped in the absorber. The reshaping delay is still small, however, amounting to about a 5% correction to ΔtG .) Because the near-resonance spectral components were removed, the large values of the group delay function near resonance do not play a role in the expectation. Thus, the pulse passing through the amplifier can exhibit superluminal behavior over a relatively long distance. If the full gaussian is sent through the amplifier without this preparation, superluminal behavior persists over a displacement of only Δr=7×10-3 c/γ before becoming subluminal (and experiencing extreme amplification). This transition occurs when the amplitude of the on-resonance peak in ρ(r, ω) has roughly the same amplitude as the peak at the carrier frequency ω¯ .

In addition to providing an example of how superluminal behavior can be extended over longer displacements, the Bolda example discussed above demonstrates the elegance and utility of the result (27). In Ref. [9

9. E. L. Bolda, J. C. Garrison, and R. Y. Chiao, “Optical Pulse Propagation at Negative Group Velocities due to a Nearby Gain Line,” Phys. Rev. A 49, 2938–2947 (1994). [CrossRef] [PubMed]

] the prediction of a negative delay time is made using the traditional series expansion approach. The analytic motivation for introducing the absorber is simply that “frequencies very close to resonance are amplified by the gain medium to such an extent that the power-series expansion breaks down.” The new context gives clear insight into the physical reasons why the power series approach cannot be used to describe the pulse in the amplifier without first sending it through the absorber. If the unabsorbed pulse were sent into the amplifier, as the near-resonance frequency components were added to the pulse its velocity would change rapidly from superluminal to highly subluminal (as in fig. 7). Since the expansion terms evaluated at a single carrier frequency are useful to describe roughly constant velocity, they fail to describe the situation.

6 Conclusion

To summarize, we have found a context wherein the group delay function always retains meaning in a uniform linear medium. The results apply to pulses of arbitrary bandwidth and temporal complexity, even in the presence of anomalous dispersion where the pulse shape can evolve dramatically during propagation. The formalism distinguishes the roles of group velocity and spectral reshaping in determining total delay time. The net group delay is averaged over the final spectrum after the propagation takes place. The reshaping delay or displacement is the effect of eliminating from (or adding to) the original pulse the spectral amplitude that is lost (or gained) during propagation. This formalism does not attempt to identify a specific velocity as being representative of the overall pulse propagation as can be done for narrowband pulses. Rather it shows that a linear superposition of the group delays is necessary to determine total delay time.

We have illustrated the use of this formalism in situations where electromagnetic pulses propagate with superluminal delay times. We note that the formalism tracks the temporal ‘center of mass’ of the arrival of a pulse’s electromagnetic field energy only. It is important to realize that a superluminal delay time for the field energy center of mass does not require superluminal energy transport. For example, energy stored in the medium downstream from a pulse’s centroid may convert to field energy. The appearance of downstream field energy can make the pulse’s center of mass move superluminally while no actual energy transport occurs faster than c. In the absorbing case, the medium can preferentially absorb the trailing portion of a pulse, resulting in an advance of the pulse’s center of mass. In both of these cases the medium treats the leading and trailing portions of the pulse differently. In a companion article [17

17. M. Ware, S. A. Glasgow, and J. Peatross “Energy Transport in Linear Dielectrics,” Opt. Express (Submitted).

], we discuss why dielectric media behave in this manner as well as addressing the issue of energy transport velocity.

References and links

1.

J. Peatross, S. A. Glasgow, and M. Ware, “Average Energy Flow of Optical Pulses in Dispersive Media,” Phys. Rev. Lett. 84, 2370–2373 (2000). [CrossRef] [PubMed]

2.

L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).

3.

M. Born and E. Wolf, Principles of Optics, 7th Ed. (Cambridge, 1999), pp. 19–24.

4.

J. D. Jackson, Classical Electrodynamics, 3rd Ed. (Wiley, New York, 1998), pp. 323, 330–335.

5.

R. Trebino, K. W. DeLong, D. N. Fittinghoff, J. N. Sweetser, M. A. Krumbugel, and B. A. Richman, “Measuring Ultrashort Laser Pulses in the Time-Frequency Domain Using Frequency-Resolved Optical Gating,” Rev. Sci. Instrum. 68, 3277–3295 (1997). [CrossRef]

6.

K. E. Oughstun and H. Xiao, “Failure of the Quasimonochromatic Approximation for Ultrashort Pulse Propagation in a Dispersive, Attenuative Medium,” Phys. Rev. Lett. 78, 642–645 (1997). [CrossRef]

7.

C. G. B. Garrett and D. E. McCumber, “Propagation of a Gaussian Light Pulse through an Anomalous Dispersion Medium,” Phys. Rev. A 1, 305–313 (1970). [CrossRef]

8.

R. Y. Chiao, “Superluminal (but Causal) Propagation of Wave Packets in Transparent Media with Inverted Atomic Populations,” Phys. Rev. A. 48, R34–R37 (1993). [CrossRef] [PubMed]

9.

E. L. Bolda, J. C. Garrison, and R. Y. Chiao, “Optical Pulse Propagation at Negative Group Velocities due to a Nearby Gain Line,” Phys. Rev. A 49, 2938–2947 (1994). [CrossRef] [PubMed]

10.

R. Y. Chiao and A. M. Steinberg, “Tunneling Times and Superluminality,” Progress in Optics 37, pp. 347–406 (Emil Wolf ed., Elsevier, Amsterdam, 1997). [CrossRef]

11.

S. Chu and S. Wong, “Linear Pulse Propagation in an Absorbing Medium,” Phys. Rev. Lett. 48, 738–741 (1982). [CrossRef]

12.

L. J. Wang, A. Kuzmmich, and A. Dogariu, “Gain-Assisted Superluminal Light Propagation,” Nature 406, 277–279 (2000). [CrossRef] [PubMed]

13.

R. L. Smith, “The Velocities of Light,” Am. J. Phys. 38, 978–983 (1970). [CrossRef]

14.

M. Ware, W. E. Dibble, S. A. Glasgow, and J. Peatross, “Energy Flow in Angularly Dispersive Optical Systems,” J. Opt. Soc. Am. B 18839–845 (2001). [CrossRef]

15.

R. Loudon, “The Propagation of Electromagnetic Energy through an Absorbing Dielectric,” J. Phys. A 3, 233–245 (1970). [CrossRef]

16.

Md. Aminul Islam Talukder, Yoshimitsu Amagishi, and Makoto Tomita, “Superluminal to Subluminal Transition in the Pulse Propagation in a Resonantly Absorbing Medium,” Phys. Rev. Lett. 86, 3546–3549 (2001). [CrossRef]

17.

M. Ware, S. A. Glasgow, and J. Peatross “Energy Transport in Linear Dielectrics,” Opt. Express (Submitted).

OCIS Codes
(260.2030) Physical optics : Dispersion
(260.2110) Physical optics : Electromagnetic optics
(350.5500) Other areas of optics : Propagation

ToC Category:
Research Papers

History
Original Manuscript: September 19, 2001
Published: November 5, 2001

Citation
Michael Ware, S. Glasgow, and Justin Peatross, "Role of group velocity in tracking field energy in linear dielectrics," Opt. Express 9, 506-518 (2001)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-9-10-506


Sort:  Journal  |  Reset  

References

  1. J. Peatross, S. A. Glasgow, and M. Ware, "Average Energy Flow of Optical Pulses in Dispersive Media," Phys. Rev. Lett. 84, 2370-2373 (2000). [CrossRef] [PubMed]
  2. L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).
  3. M. Born and E. Wolf, Principles of Optics, 7th Ed. (Cambridge, 1999), pp. 19-24.
  4. J. D. Jackson, Classical Electrodynamics, 3rd Ed. (Wiley, New York, 1998), pp. 323, 330-335.
  5. R. Trebino, K. W. DeLong, D. N. Fittinghoff, J. N. Sweetser, M. A. Krumbugel, and B. A. Richman, "Measuring Ultrashort Laser Pulses in the Time-Frequency Domain Using Frequency-Resolved Optical Gating," Rev. Sci. Instrum. 68, 3277-3295 (1997). [CrossRef]
  6. K. E. Oughstun and H. Xiao, "Failure of the Quasimonochromatic Approximation for Ultrashort Pulse Propagation in a Dispersive, Attenuative Medium," Phys. Rev. Lett. 78, 642-645 (1997). [CrossRef]
  7. C. G. B. Garrett and D. E. McCumber, "Propagation of a Gaussian Light Pulse through an Anomalous Dispersion Medium," Phys. Rev. A 1, 305-313 (1970). [CrossRef]
  8. R. Y. Chiao, "Superluminal (but Causal) Propagation of Wave Packets in Transparent Media with Inverted Atomic Populations," Phys. Rev. A. 48, R34-R37 (1993). [CrossRef] [PubMed]
  9. E. L. Bolda, J. C. Garrison, and R. Y. Chiao, "Optical Pulse Propagation at Negative Group Velocities due to a Nearby Gain Line," Phys. Rev. A 49, 2938-2947 (1994). [CrossRef] [PubMed]
  10. Y. Chiao and A. M. Steinberg, "Tunneling Times and Superluminality," Progress in Optics 37, pp. 347-406 (Emil Wolf ed., Elsevier, Amsterdam, 1997). [CrossRef]
  11. S. Chu and S. Wong, "Linear Pulse Propagation in an Absorbing Medium," Phys. Rev. Lett. 48, 738-741 (1982). [CrossRef]
  12. L. J. Wang, A. Kuzmmich, and A. Dogariu, "Gain-Assisted Superluminal Light Propagation," Nature 406, 277-279 (2000). [CrossRef] [PubMed]
  13. R. L. Smith, "The Velocities of Light," Am. J. Phys. 38, 978-983 (1970). [CrossRef]
  14. M. Ware, W. E. Dibble, S. A. Glasgow, and J. Peatross, "Energy Flow in Angularly Dispersive Optical Systems," J. Opt. Soc. Am. B 18 839-845 (2001) . [CrossRef]
  15. R. Loudon, "The Propagation of Electromagnetic Energy through an Absorbing Dielectric," J. Phys. A 3, 233-245 (1970). [CrossRef]
  16. Md. Aminul Islam Talukder, Yoshimitsu Amagishi, and Makoto Tomita, "Superluminal to Subluminal Transition in the Pulse Propagation in a Resonantly Absorbing Medium," Phys. Rev. Lett. 86, 3546-3549 (2001). [CrossRef]
  17. M. Ware, S. A. Glasgow, and J. Peatross "Energy Transport in Linear Dielectrics," Opt. Express 9, 519-532 (2001), http://www.opticsexpress.org/oearchive/source/35289.htm

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Supplementary Material


» Media 1: MOV (1424 KB)     
» Media 2: MOV (1746 KB)     
» Media 3: MOV (1742 KB)     

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited