## Optical precursor fields in nonlinear pulse dynamics |

Optics Express, Vol. 18, Issue 22, pp. 23104-23120 (2010)

http://dx.doi.org/10.1364/OE.18.023104

Acrobat PDF (2991 KB)

### Abstract

Under certain conditions, ultrashort pulse dynamics in a linear dispersive medium with absorption result in the appearance of optical precursors that dominate the pulse evolution for large propagation distances as the peak amplitude in the initial pulse spectrum decays exponentially. The effects of a nonlinear medium response on this precursor formation is considered using the split-step Fourier method. Comparison of the nonlinear pulse evolution when the full dispersion is used to that when a quadratic Taylor series approximation of the wave number is used shows that the group velocity approximation misses the precursor fields entirely.

© 2010 Optical Society of America

## 1. Introduction

9. K. E. Oughstun and G. C. Sherman, “Propagation of electromagnetic pulses in a linear dispersive medium with absorption (the Lorentz medium),” J. Opt. Soc. Am. B **5**, 4, 817–849 (1988). [CrossRef]

10. N. A. Cartwright and K. E. Oughstun, “Uniform asymptotics applied to ultrawideband pulse propagation,” SIAM Review **49**, 4, 628–648 (2007). [CrossRef]

7. K. E. Oughstun, *Electromagnetic and Optical Pulse Propagation1: Spectral Representations in Temporally Dispersive Media* (Springer, New York, 2006). [PubMed]

*ω*

_{0}, damping

*δ*, and plasma frequency

*ω*. The absorption band, which corresponds to the region of anomalous dispersion, then extends from

_{p}*E*(0,

*t*) = 0 for all

*t*< 0, then

*E*(

*z*,

*t*) = 0 for all

*θ*< 1 with

*z*> 0.

*θ*≥ 1, the asymptotic theory [2, 3, 6, 9

9. K. E. Oughstun and G. C. Sherman, “Propagation of electromagnetic pulses in a linear dispersive medium with absorption (the Lorentz medium),” J. Opt. Soc. Am. B **5**, 4, 817–849 (1988). [CrossRef]

*E*(

_{S}*z*,

*t*) is the first forerunner [1, 2, 3] or Sommerfeld precursor due to a pair of distant saddle points describing the high-frequency (|

*ω*| ≥

*ω*

_{1}) response of the dispersive medium,

*E*(

_{B}*z*,

*t*) is the second forerunner [2, 3] or Brillouin precursor due to a pair of near saddle points describing the low-frequency (|

*ω*| ≤

*ω*

_{0}) response, and

*E*(

_{c}*z*,

*t*) is due to the residue of the pole contribution [8, 10

10. N. A. Cartwright and K. E. Oughstun, “Uniform asymptotics applied to ultrawideband pulse propagation,” SIAM Review **49**, 4, 628–648 (2007). [CrossRef]

11. K. E. Oughstun and C. M. Balictsis, “Gaussian Pulse Propagation in a Dispersive, Absorbing Dielectric,” Phys. Rev. Lett . **77**, 11, 2210–2213 (1996). [CrossRef]

*Ẽ*(0,

_{G}*ω*) is an entire function of complex

*ω*. Of particular importance here is that the peak amplitude of the Brillouin precursor

*E*(

_{B}*z*,

*t*) experiences zero exponential decay with propagation distance

*z*> 0, decreasing algebraically as

*z*

^{−1/2}in the dispersive, absorptive medium for

*δ*> 0 bounded away from zero [6, 8, 9

9. K. E. Oughstun and G. C. Sherman, “Propagation of electromagnetic pulses in a linear dispersive medium with absorption (the Lorentz medium),” J. Opt. Soc. Am. B **5**, 4, 817–849 (1988). [CrossRef]

10. N. A. Cartwright and K. E. Oughstun, “Uniform asymptotics applied to ultrawideband pulse propagation,” SIAM Review **49**, 4, 628–648 (2007). [CrossRef]

12. T. H. Havelock, “The propagation of groups of waves in dispersive media,” Proc. R. Soc. A **LXXXI**, 398 (1908). [CrossRef]

*k*(

*ω*), assumed there to be real-valued, is expanded in a Taylor series about some wave number value

*k*

_{0}that the spectrum of the wave group is clustered about, where [13] “the range of integration is supposed to be small and the amplitude, phase and velocity of the members of the group are assumed to be continuous, slowly varying, functions” of the wave number. Notice that Havelock’s group method [12

12. T. H. Havelock, “The propagation of groups of waves in dispersive media,” Proc. R. Soc. A **LXXXI**, 398 (1908). [CrossRef]

*k*

_{0}is the stationary phase point of the wavenumber

*k*(

*ω*). This apparently subtle change in the value of

*k*

_{0}results in significant consequences on the accuracy of the resulting group velocity description based on this method [8, 15

15. K. E. Oughstun and H. Xiao, “Failure of the quasimonochromatic approximation for ultrashort pulse propagation in a dispersive, attenuative medium,” Phys. Rev. Lett. **78**, 4, 642–645 (1997). [CrossRef]

16. H. Xiao and K. E. Oughstun, “Failure of the group velocity description for ultrawideband pulse propagation in a double resonance Lorentz model dielectric,” J. Opt. Soc. Am. B **16**, 10, 1773–1785 (1999). [CrossRef]

*k̃*′(

*ω*) ≡

*∂k̃*(

*ω*)/

*∂ω*(the inverse of the group velocity) and

*k̃*″(

*ω*) ≡

*∂*

^{2}

*k̃*(

*ω*)/

*∂*

*ω*

^{2}(the so-called group velocity dispersion), is illustrated in Fig. 1 with

*ω*= 1.0 × 10

_{c}^{15}

*r*/

*s*for a moderately dispersive resonance line at

*ω*

_{0}= 2.4 × 10

^{15}

*r*/

*s*with

*ω*/

_{p}*ω*

_{0}≃ 0.00127 and

*δ*/

*ω*

_{0}= 0.025., the upper plot describing the scaled real part (

*c*/

*ω*)[

*β*(

*ω*) –

*β*(

*ω*)] and the lower plot the imaginary part

_{c}*α*(

*ω*) of

*k̃*(

*ω*). Although

*β*(

*ω*) is accurately modeled over a finite frequency interval below resonance (because

*ω*<

_{c}*ω*

_{0}) in this case, diverging from the actual behavior both when

*ω*increases sufficiently above or below the medium resonance frequency

*ω*

_{0},

*α*(

*ω*) is overestimated for all frequencies above and below the carrier frequency

*ω*. Because of this, the frequency dependence of

_{c}*α*(

*ω*) is typically neglected in the group velocity approximation where

*α*(

*ω*) ≈

*α*(

*ω*). The inclusion of higher-order terms in the Taylor series approximation of

_{c}*k̃*(

*ω*) beyond the quadratic dispersion approximation given in Eq. (6) has been shown [8, 15

15. K. E. Oughstun and H. Xiao, “Failure of the quasimonochromatic approximation for ultrashort pulse propagation in a dispersive, attenuative medium,” Phys. Rev. Lett. **78**, 4, 642–645 (1997). [CrossRef]

16. H. Xiao and K. E. Oughstun, “Failure of the group velocity description for ultrawideband pulse propagation in a double resonance Lorentz model dielectric,” J. Opt. Soc. Am. B **16**, 10, 1773–1785 (1999). [CrossRef]

15. K. E. Oughstun and H. Xiao, “Failure of the quasimonochromatic approximation for ultrashort pulse propagation in a dispersive, attenuative medium,” Phys. Rev. Lett. **78**, 4, 642–645 (1997). [CrossRef]

16. H. Xiao and K. E. Oughstun, “Failure of the group velocity description for ultrawideband pulse propagation in a double resonance Lorentz model dielectric,” J. Opt. Soc. Am. B **16**, 10, 1773–1785 (1999). [CrossRef]

## 2. Formulation of the nonlinear problem

*D*(

**r**,

*t*) =

*E*(

**r**,

*t*) + 4

*π*

*P*(

**r**,

*t*). With substitution of the complex phasor representation

*E*(

_{ωc}**r**,

*t*) ≡

*A*(

**r**,

*t*)

*e*

^{i(φ(r,t)–ωct)}, where

*E*(

**r**,

*t*) = ℜ{

*E*(

_{ωc}**r**,

*t*)}, the linear material polarization response is found as [8]

*P*(

_{ωc}**r**,

*t*) =

*χ*(

_{e}*ω*)

_{c}*A*(

**r**,

*t*)

*e*

^{i(φ(r,t)–ωct)}, where

*P*(

**r**,

*t*) = ℜ{

*P*

_{ωc}(

**r**,

*t*)} and

*χ*(

_{e}*ω*) denotes the electric susceptibility, and the wave equation (7) becomes

*γ*is the nonlinear-index coefficient.

*D*indicates that this describes the linear dispersive part of the wave field. The pair of wave equations (11) and (12) then describe nonlinear dispersive pulse propagation, including diffraction effects, without the usual approximations associated with the group velocity method. An analogous formulation has been presented by Laine and Friberg [19

19. T. A. Laine and A. T. Friberg, “Self-guided waves and exact solutions of the nonlinear Helmholtz equation,” J. Opt. Soc. Am. B **17**, 5, 751–757 (2000). [CrossRef]

*z*-direction through the nonlinear dispersive material. Equations (11) and (12) then admit the respective solutions where

*h*is the numerical step size taken in the split-step method. The step size is taken here as

*h*=

*z*/5, where

_{d}*z*≡

_{d}*α*

^{−1}(

*ω*) =

_{c}*c*/(

*ω*

_{c}*n*(

_{i}*ω*)) is the

_{c}*e*

^{−1}absorption depth. Notice that

*z*≃ 12.18

_{d}*m*for the material parameters used in this paper with

*γ*= 1.

*z*=

*z*and the result is then transformed to give

_{j}*E*(

*z*+

_{j}*h*,

*t*). This linear propagation step is then followed by a nonlinear propagation step over the distance

*h*/2 through Eq. (13), as suggested by Agarwal [17] in order to improve accuracy of the split-step Fourier method. Notice that after the application of the nonlinear operator in Eq. (14), the real part of the resultant field is taken in order to assure physical results.

## 3. Heaviside step-function signal evolution

*E*(0,

_{H}*t*) =

*u*(

_{H}*t*) sin (

*ω*

_{c}*t*) with

*u*(

_{H}*t*) = 0 for

*t*< 0 and

*u*(

_{H}*t*) = 1 for

*t*> 0 in the linear and cubic nonlinear cases, illustrated in Fig. 2 at five absorption depths (

*z*= 5

*z*). The linear result is precisely that described by the modern asymptotic theory [8, 10

_{d}**49**, 4, 628–648 (2007). [CrossRef]

*E*(

_{H}*z*,

*t*) = 0 for all superluminal space-time points

*θ*< 1 with

*z*> 0. The propagated wave field then arrives at the luminal space-time point

*θ*= 1 with the onset of the Sommerfeld precursor

*E*(

_{HS}*z*,

*t*), followed by the slower Brillouin precursor

*E*(

_{HB}*z*,

*t*) whose peak amplitude point travels at the velocity

*v*=

_{B}*c*/

*θ*

_{0}=

*c*/

*n*(0), decaying only as

*z*

^{−1/2}, which is then followed by the main signal

*A*(

_{Hc}*z*,

*t*). The nonlinear signal evolution is remarkably similar, vanishing for superluminal space-time points

*θ*< 1, the propagated wave field arriving at

*θ*= 1 with the onset of a Sommerfeld precursor that is essentially identical with that in the linear case. This is then followed by a Brillouin precursor whose peak amplitude is ∼ 82% of that in the linear case with a nearly identical oscillation frequency. Notice that the cubic nonlinearity generates a small amplitude frequency component at 3

*ω*into the propagated wave field as the signal arrival at

_{c}*ω*=

*ω*is approached. The propagated signal spectrum presented in Fig. 3 shows that the cubic nonlinearity generates frequency components at the odd harmonics 3

_{c}*ω*, 5

_{c}*ω*, 7

_{c}*ω*,... of the carrier frequency as well as filling in the linear spectral loss in the absorption band about

_{c}*ω*

_{0}. This odd-harmonic frequency structure, which is missing when the quadratic wave number approximation in Eq. (6) is employed, may account for the reduced energy evident in the precursor field structure in the nonlinear case.

## 4. Gaussian envelope pulse evolution

*φ*= 0 when

*N*< 1 (in order to ensure zero pulse area) and

_{osc}*φ*= 3

*π*/2 when

*N*≥ 1 (so that the maxima of the envelope and signal coincide when an integer number

_{osc}*N*of oscillations are contained in the pulse width). Here 2

_{osc}*τ*

_{0}denotes the temporal width of the gaussian envelope at the

*e*

^{−1}amplitude points. In the numerical examples presented here, the initial pulse width 2

*τ*

_{0}was varied from

*T*, where

_{c}*T*≡ 2

_{c}*π*/

*ω*is the oscillation period of the carrier wave. For the below resonance carrier frequency

_{c}*ω*= 1 × 10

_{c}^{15}

*r*/

*s*, this corresponds to initial pulse widths ranging from the ultrashort

*T*≃ 628

_{c}*fs*. In all cases considered, the computed initial pulse area is found to be less than the machine epsilon (∼ 2.2 × 10

^{−16}). The spectral magnitude |

*ũ*(

*ω*)| for each extreme case is illustrated in Fig. 4 in reference to the linear material phase dispersion. The dashed curve in the figure describes the magnitude of the scaled Heaviside step function signal spectrum |

*ω*–

*ω*|

_{c}^{–1}. Notice that the ultrashort

*T*pulse spectrum (red shaded region) is wideband below resonance and the narrowband pulse spectrum (blue shaded region) is quasimonochromatic. The Sommerfeld and Brillouin precursor fields that are characteristic of the full material dispersion will clearly dominate the field evolution in the ultrawideband but not in the narrow-band case. Because the intermediate wideband case fills the spectral region below the material resonance, only the low-frequency Brillouin precursor will be present in that case.

_{c}*T*∼

_{r}*δ*

^{−1}[8]. The first and second time derivatives of this function are given respectively by The inflection point of

*u*′(

*t*) is given by the appropriate zero of

*u*″(

*t*), which has zeroes given by

*t*=±

*τ*

_{0}/(2)

^{1/2}. Substitution of the negative root, which corresponds to the maximum, into the expression for

*u*′(

*t*) and equating the result to

*δ*then yields the critical value which corresponds to a minimum

*N*= 4.6 for

_{osc}*ω*= 1 × 10

_{c}^{15}

*r*/

*s*. Below this critical value, either one or both of the precursor fields will be fully realized in the propagated wave field, depending on the value of the pulse carrier frequency

*ω*in comparison to the medium resonance frequency

_{c}*ω*

_{0}. On the other hand, as

*τ*

_{0}is increased above

*τ*, the observed pulse distortion will approach that described by the group velocity approximation, as is now shown.

_{c}*τ*

_{0}≃ 18.8

*fs*) is presented in Fig. 5 at 5 absorption depths (

*z*= 5

*z*) into the dispersive medium whose linear frequency dispersion is illustrated in Fig. 1, the full dispersion response being employed in both cases. Because the pulse carrier frequency

_{d}*ω*is sufficiently below the material resonance frequency

_{c}*ω*

_{0}so that there is negligible spectral energy above

*ω*

_{1}, as seen in Fig. 4, the high-frequency Sommerfeld precursor is essentially absent from the propagated pulse, the pulse structure then being dominated by the Brillouin precursor [11

11. K. E. Oughstun and C. M. Balictsis, “Gaussian Pulse Propagation in a Dispersive, Absorbing Dielectric,” Phys. Rev. Lett . **77**, 11, 2210–2213 (1996). [CrossRef]

*τ*

_{0}is decreased below a single oscillation, however, and the initial pulse spectrum becomes increasingly ultrawideband, the Sommerfeld precursor becomes increasingly dominant in the propagated field structure (see Fig. 6). Comparison of the propagated linear and nonlinear pulse shapes presented in Fig. 5 shows that the nonlinearity primarily decreases the amplitude of the gaussian Brillouin precursor evolution, the peak amplitude being decreased to ∼ 82% of its linear value, the same amount obtained for the step function signal. This decrease may then be attributed to the generation of odd-harmonics by the nonlinearity.

*τ*

_{0}is increased from the ultrawideband to the quasimonochromatic spectral extremes when the full and approximate [see Eq. (6)] dispersion relations are employed in the numerical propagation model. For the ultrawideband 3.14

*fs*pulse case illustrated in Fig. 6, the propagated gaussian pulse has separated into an above resonance gaussian Sommerfeld precursor and a below resonance gaussian Brillouin precursor component, as described in [11

11. K. E. Oughstun and C. M. Balictsis, “Gaussian Pulse Propagation in a Dispersive, Absorbing Dielectric,” Phys. Rev. Lett . **77**, 11, 2210–2213 (1996). [CrossRef]

*τ*

_{0}is increased with

*ω*fixed below resonance, the Sommerfeld precursor component becomes negligible in comparison to the Brillouin precursor component (the opposite will be found when the pulse carrier frequency is situated sufficiently far above resonance), as seen in Fig. 7 for the 18.8

_{c}*fs*pulse case and Fig. 8 for the 31.4

*fs*pulse case. In both cases, the quadratic group velocity approximation only describes the trailing edge of the Brillouin precursor component with any accuracy. At this point when the critical pulse width value

*τ*is exceeded, a gradual transition from the wideband, precursor dominated behavior to the narrowband quasimonochromatic behavior is observed, as evidenced in Figs. 9–12.

_{c}*τ*

_{0}≃ 18.8

*fs*using the full material dispersion and the quadratic approximation of this dispersion is presented in Figs. 14 and 15, respectively, for the initial pulse evolution as

*z*increases from 0 to

*z*, and in Figs. 16 and 17 for the mature pulse evolution [6, 8] as

_{d}*z*increases above

*z*. The evolution of the wideband pulse into a gaussian Brillouin precursor with minimal attenuation is clearly evident in Figs. 14 and 16 (full dispersion model) while it is noticeably absent in Figs. 15 and 17 (quadratic approximation of the dielectric dispersion). In both cases there is a sharp drop in the pulse amplitude as the propagation distance increases to a single absorption depth. This is then followed by a transition to the precursor dominated field evolution in the full dispersion case illustrated in Figs. 14 and 16, where the peak amplitude decay switches from a supra-exponential to a sub-exponential decay, as illustrated by the blue curve in Fig. 18. Notice that the peak amplitude decay for the group velocity approximation, indicated by the green curve in Fig. 18, remains below the Beer’s law limit of pure exponential decay (indicated by the dashed black curve in the figure) over the propagation distance domain

_{d}*z*/

*z*∈ [0, 5] considered in the numerical study presented here. The

_{d}*z*

^{−1/2}algebraic decay that is a characteristic of the Brillouin precursor in the linear case [6, 8] is greatly exceeded in the nonlinear case, the numerically determined slope

*p*of the full dispersion curve in Fig. 18, which describes the

*z*

^{−p}peak amplitude decay, decreasing from

*p*∼ 0.56 at

*z*/

*z*∼ 0.5 to

_{d}*p*∼ 0.067 at

*z*/

*z*∼ 5. Again, this is due to the cubic nonlinearity which also significantly effects the peak amplitude decay in the group velocity approximation as

_{d}*z*/

*z*increases above unity.

_{d}## 5. Other pulse shapes

*τ*

_{0}= 3 × 10

^{−15}

*s*and

*ω*= 1 × 10

_{c}^{15}

*r*/

*s*, which corresponds to the 3 oscillation gaussian pulse case illustrated in Fig. 7. As can be seen, the two results are quite similar in pulse structure, this being due to the dominance of the Brillouin precursor.

## 6. Conclusion

- the precursor fields that are a characteristic of the linear material dispersion persist in the nonlinear case;
- these precursor fields dominate gaussian pulse evolution when the initial pulse width 2
*τ*_{0}is on the order of or less than the critical pulse width 2*τ*, where_{c}*τ*∼_{c}*δ*^{−1}is the relaxation time of the medium response [see Eq. (18)]; - as found in the linear case for the gaussian envelope pulse [11
**77**, 11, 2210–2213 (1996). [CrossRef]*τ*_{0}≪ 2*τ*;_{c} - the characteristic
*z*^{−1/2}peak amplitude decay of the gaussian Brillouin precursor in the linear case is enhanced by the cubic nonlinearity; - the group velocity description using the quadratic dispersion approximation fails to accurately describe gaussian pulse evolution when the initial pulse width 2
*τ*_{0}is on the order of or less than the critical pulse width 2*τ*(the accuracy does not improve with the inclusion of higher-order terms);_{c} - in addition, because the group velocity approximation overestimates the linear material attenuation away from the pulse carrier frequency
*ω*, the odd-order harmonics introduced by the cubic nonlinearity are absent in the group velocity description of the pulse evolution._{c}

## Acknowledgments

## References and links

1. | A. Sommerfeld, “Über die fortpflanzung des lichtes in disperdierenden medien,” Ann. Phys. (Leipzig) |

2. | L. Brillouin, “Über die fortpflanzung des licht in disperdierenden medien,” Ann. Phys. (Leipzig) |

3. | L. Brillouin, |

4. | P. Debye, “Näherungsformeln für die zylinderfunktionen für grosse werte des arguments und unbeschränkt verander liche werte des index,” Math. Ann. |

5. | J. A. Stratton, |

6. | K. E. Oughstun and G. C. Sherman, |

7. | K. E. Oughstun, |

8. | K. E. Oughstun, |

9. | K. E. Oughstun and G. C. Sherman, “Propagation of electromagnetic pulses in a linear dispersive medium with absorption (the Lorentz medium),” J. Opt. Soc. Am. B |

10. | N. A. Cartwright and K. E. Oughstun, “Uniform asymptotics applied to ultrawideband pulse propagation,” SIAM Review |

11. | K. E. Oughstun and C. M. Balictsis, “Gaussian Pulse Propagation in a Dispersive, Absorbing Dielectric,” Phys. Rev. Lett . |

12. | T. H. Havelock, “The propagation of groups of waves in dispersive media,” Proc. R. Soc. A |

13. | T. H. Havelock, |

14. | L. Kelvin, “On the waves produced by a single impulse in water of any depth, or in a dispersive medium,” Proc. Roy. Soc . |

15. | K. E. Oughstun and H. Xiao, “Failure of the quasimonochromatic approximation for ultrashort pulse propagation in a dispersive, attenuative medium,” Phys. Rev. Lett. |

16. | H. Xiao and K. E. Oughstun, “Failure of the group velocity description for ultrawideband pulse propagation in a double resonance Lorentz model dielectric,” J. Opt. Soc. Am. B |

17. | G. P. Agrawal, |

18. | R. H. Hardin and F. D. Tappert, |

19. | T. A. Laine and A. T. Friberg, “Self-guided waves and exact solutions of the nonlinear Helmholtz equation,” J. Opt. Soc. Am. B |

**OCIS Codes**

(260.2030) Physical optics : Dispersion

(320.5550) Ultrafast optics : Pulses

**ToC Category:**

Ultrafast Optics

**History**

Original Manuscript: September 9, 2010

Revised Manuscript: October 11, 2010

Manuscript Accepted: October 11, 2010

Published: October 18, 2010

**Citation**

Chris L. Palombini and Kurt E. Oughstun, "Optical precursor fields in nonlinear pulse dynamics," Opt. Express **18**, 23104-23120 (2010)

http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-18-22-23104

Sort: Year | Journal | Reset

### References

- A. Sommerfeld, “Über die fortpflanzung des lichtes in disperdierenden medien,” Ann. Phys. (Leipzig) 44, 4, 177–202 (1914).
- L. Brillouin, “Über die fortpflanzung des licht in disperdierenden medien,” Ann. Phys. (Leipzig) 44, 4, 177–202 (1914).
- L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).
- P. Debye, “Näherungsformeln für die zylinderfunktionen für grosse werte des arguments und unbeschränkt verander liche werte des index,” Math. Ann. 67, 535–558 (1909). [CrossRef]
- J. A. Stratton, Electromagnetic Theory, (McGraw-Hill, New York, 1941), §5.18.
- K. E. Oughstun, and G. C. Sherman, Electromagnetic Pulse Propagation in Causal Dielectrics (Springer-Verlag, Berlin-Heidelberg, 1994).
- K. E. Oughstun, Electromagnetic and Optical Pulse Propagation1: Spectral Representations in Temporally Dispersive Media (Springer, New York, 2006). [PubMed]
- K. E. Oughstun, Electromagnetic and Optical Pulse Propagation2: Temporal Pulse Dynamics in Dispersive, Attenuative Media (Springer, New York, 2009).
- K. E. Oughstun, and G. C. Sherman, “Propagation of electromagnetic pulses in a linear dispersive medium with absorption (the Lorentz medium),” J. Opt. Soc. Am. B 5(4), 817–849 (1988). [CrossRef]
- N. A. Cartwright, and K. E. Oughstun, “Uniform asymptotics applied to ultrawideband pulse propagation,” SIAM Rev. 49(4), 628–648 (2007). [CrossRef]
- K. E. Oughstun, and C. M. Balictsis, “Gaussian Pulse Propagation in a Dispersive, Absorbing Dielectric,” Phys. Rev. Lett. 77(11), 2210–2213 (1996). [CrossRef]
- T. H. Havelock, “The propagation of groups of waves in dispersive media,” Proc. R. Soc. London, Ser. A LXXXI, 398 (1908). [CrossRef]
- T. H. Havelock, The Propagation of Disturbances in Dispersive Media (Cambridge U. Press, Cambridge, 1914).
- L. Kelvin, “On the waves produced by a single impulse in water of any depth, or in a dispersive medium,” Proc. Roy. Soc., London XLII, 80 (1887).
- K. E. Oughstun, and H. Xiao, “Failure of the quasimonochromatic approximation for ultrashort pulse propagation in a dispersive, attenuative medium,” Phys. Rev. Lett. 78(4), 642–645 (1997). [CrossRef]
- H. Xiao, and K. E. Oughstun, “Failure of the group velocity description for ultrawideband pulse propagation in a double resonance Lorentz model dielectric,” J. Opt. Soc. Am. B 16(10), 1773–1785 (1999). [CrossRef]
- G. P. Agrawal, Nonlinear Fiber Optics (Academic Press, San Diego, 1989).
- R. H. Hardin, and F. D. Tappert, “Applications of the split-step Fourier method to the solution of nonlinear and variable coefficient wave equations,” SIAM Rev. 15(2), 423 (1973).
- T. A. Laine, and A. T. Friberg, “Self-guided waves and exact solutions of the nonlinear Helmholtz equation,” J. Opt. Soc. Am. B 17(5), 751–757 (2000). [CrossRef]

## 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.

### Figures

Fig. 1 |
Fig. 2 |
Fig. 3 |

Fig. 4 |
Fig. 5 |
Fig. 6 |

Fig. 7 |
Fig. 8 |
Fig. 9 |

Fig. 10 |
Fig. 11 |
Fig. 12 |

Fig. 13 |
Fig. 14 |
Fig. 15 |

Fig. 16 |
Fig. 17 |
Fig. 18 |

Fig. 19 |
||

« Previous Article | Next Article »

OSA is a member of CrossRef.