Open Access
Add to BibSonomyAbstract
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
The asymptotic description of ultrawideband dispersive pulse propagation in a Lorentz model dielectric has its origin in the classic research by Sommerfeld [1] and Brillouin [2, 3] based on Debye’s steepest descent method [4]. This seminal work established the physical phenomena of forerunners, or precursor fields [5], that were originally associated with the evolution of a Heaviside step-function signal. This description is derived from the exact Fourier-Laplace integral representation of the propagated linearly polarized plane wave field [6, 7, 8]
for z ≥ 0 with fixed a larger than the abscissa of absolute convergence [5, 6, 7, 8] for the spectrum of the initial plane wave pulse E(0, t) = f(t) at z = 0. The temporal Fourier spectrum Ẽ (z,ω) of the optical wave field E(z,t) satisfies the Helmholtz equation with complex wave number k̃(ω) = β(ω) + iα(ω) = (ω/c)n(ω) in the temporally dispersive medium with complex index of refraction n(ω) = nr(ω)+ini(ω) whose real nr(ω) ≡ ℜ{n(ω)} and imaginary ni(ω) ≡ ℑ{n(ω)} parts are related through the Kramers-Kronig relations [7]. Here β(ω) ≡ ℜ{k̃(ω)} is the propagation (or phase) factor and α(ω) ≡ ℑ{k̃(ω)} the attenuation factor for plane wave propagation in the dispersive attenuative medium. In addition, with θ ≡ ct/z a nondimensional space-time parameter defined for all z > 0.Based upon this exact integral representation, the modern asymptotic theory [9, 10] provides a uniform asymptotic description of dispersive pulse dynamics in a single-resonance Lorentz model dielectric with causal, complex index of refraction given by [7]
with resonance frequency ω0, damping δ, and plasma frequency ωp. The absorption band, which corresponds to the region of anomalous dispersion, then extends from to , where . Sommerfeld’s relativistic causality theorem then holds [1, 8], stating that if E(0, t) = 0 for all t < 0, then E(z,t) = 0 for all θ < 1 with z > 0.For θ ≥ 1, the asymptotic theory [2, 3, 6, 9] shows that the propagated wave field of an ultrawideband signal in a single-resonance Lorentz model dielectric is due to two independent sets of saddle points and pole contributions from the initial pulse spectrum so that it may then be expressed either as a superposition of component fields as
or as a linear combination of fields of this form [8]. Here ES(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, EB(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 Ec(z, t) is due to the residue of the pole contribution [8, 10] describing the signal contribution (if any). The observed pulse distortion for a Heaviside step function signal is then seen to be due to both the precursor fields and their interference with the signal contribution. For a gaussian envelope pulse, the dynamical pulse evolution is comprised of just the precursor fields [8, 11] because its spectrum ẼG(0, ω) is an entire function of complex ω. Of particular importance here is that the peak amplitude of the Brillouin precursor EB(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, 10]. Part of the purpose of this paper is to show that the precursor field components of an ultrawideband pulse persist when nonlinear effects are included.By comparison, the group velocity description of dispersive pulse propagation was formulated by Havelock [12, 13] based upon Kelvin’s stationary phase method [14]. In this approach, the wave number k(ω), assumed there to be real-valued, is expanded in a Taylor series about some wave number value k0 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, 13] is a significant departure from Kelvin’s stationary phase method [14], where k0 is the stationary phase point of the wavenumber k(ω). This apparently subtle change in the value of k0 results in significant consequences on the accuracy of the resulting group velocity description based on this method [8, 15, 16].
The accuracy of the quadratic Taylor series approximation of the complex wave number
where k̃′(ω) ≡ ∂k̃(ω)/∂ω (the inverse of the group velocity) and k̃″(ω) ≡ ∂2k̃(ω)/∂ω2 (the so-called group velocity dispersion), is illustrated in Fig. 1 with ωc = 1.0 × 1015r/s for a moderately dispersive resonance line at ω0 = 2.4 × 1015r/s with ωp/ω0 ≃ 0.00127 and δ/ω0 = 0.025., the upper plot describing the scaled real part (c/ω)[β(ω) – β(ωc)] and the lower plot the imaginary part α(ω) 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 ωc. Because of this, the frequency dependence of α(ω) is typically neglected in the group velocity approximation where α(ω) ≈ α(ωc). The inclusion of higher-order terms in the Taylor series approximation of k̃(ω) beyond the quadratic dispersion approximation given in Eq. (6) has been shown [8, 15, 16] to yield no improvement in accuracy, and hence, is meaningless.
Fig. 1 Exact (blue curves) and quadratic Taylor series approximation (green curves) of the scaled propagation factor (c/ω)β(ω) (upper plot) and the attenuation coefficient α(ω) of the complex wave number for a single resonance Lorentz model dielectric with ω0 = 2.4 × 1015r/s, δ = 6.0 × 1013r/s, ωp = 3.05 × 1012r/s.
As a consequence of these error sources, which are fundamental to the group method, the accuracy of the group velocity approximation rapidly decreases in the ultrawideband signal, ultrashort pulse limit as the pulse rise- or fall-time decreases below the characteristic relaxation time of the medium resonance [8, 15, 16]. This result is now extended to the nonlinear regime where the relative importance of the precursor fields has so far been neglected.
2. Formulation of the nonlinear problem
The formulation of the nonlinear dispersive pulse propagation problem begins with the inhomogeneous scalar wave equation (see §11.3 of [8] and §1.3 of [17])
for a linearly polarized pulse in a simple dispersive medium characterized by the constitutive relation D(r, t) = E(r, t) + 4πP(r, t). With substitution of the complex phasor representation Eωc(r,t) ≡ A(r, t)ei(φ(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)ei(φ(r,t)–ωct), where P(r, t) = ℜ{Pωc (r,t)} and χe(ω) denotes the electric susceptibility, and the wave equation (7) becomesWith the inclusion of a cubic nonlinearity [17], this wave equation may be generalized as
where γ is the nonlinear-index coefficient.It is essential that the nonlinear wave equation (9) be solved as it stands without any unjustified simplifying assumptions as is typically done in nonlinear optics (see, for example, §11.3 of Ref. [8]). This is accomplished here using the split-step Fourier method introduced by Hardin and Tappert [18] by separating it into linear (dispersive and attenuating) and nonlinear parts as
respectively. The linear part (10) can then be transformed back into Eq. (7) and from there into the Helmholtz equation [cf. Eq. (2)] where the subscript 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].Attention is now restricted to plane wave propagation in the positive 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 = zd/5, where zd ≡ α−1(ωc) = c/(ωcni(ωc)) is the e−1 absorption depth. Notice that zd ≃ 12.18m for the material parameters used in this paper with γ = 1.In each numerical simulation presented here, the pulse was propagated 5 absorption depths in order to fully display the influence of the precursor fields on the nonlinear pulse dynamics. In this numerical simulation, Eq. (14) is first applied to the pulse spectrum at z = zj and the result is then transformed to give E(zj + 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
Consider first a comparison of the dynamical field evolution of a Heaviside unit-step function modulated signal EH (0, t) = uH (t) sin (ωct) with uH (t) = 0 for t < 0 and uH (t) = 1 for t > 0 in the linear and cubic nonlinear cases, illustrated in Fig. 2 at five absorption depths (z = 5zd). The linear result is precisely that described by the modern asymptotic theory [8, 10] and, as the conditions of Sommerfeld’s relativistic causality theorem are satisfied in the linear case, EH (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 EHS(z,t), followed by the slower Brillouin precursor EHB(z,t) whose peak amplitude point travels at the velocity vB = c/θ0 = c/n(0), decaying only as z−1/2, which is then followed by the main signal AHc(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ωc 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 ω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.
Fig. 2 Heaviside step-function signal response of the dispersive material at 5 absorption depths with (green curve) and without (blue curve) the cubic nonlinearity included.
Fig. 3 Propagated spectra of the Heaviside step function signal illustrated in Fig. 2 for the linear (solid blue curve) and nonlinear (solid green curve) dispersion cases.
4. Gaussian envelope pulse evolution
The second type of input signal considered here is the gaussian envelope modulated signal
where φ = 0 when Nosc < 1 (in order to ensure zero pulse area) and φ = 3π/2 when Nosc ≥ 1 (so that the maxima of the envelope and signal coincide when an integer number Nosc of oscillations are contained in the pulse width). Here 2τ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 to 100Tc, where Tc ≡ 2π/ωc is the oscillation period of the carrier wave. For the below resonance carrier frequency ωc = 1 × 1015r/s, this corresponds to initial pulse widths ranging from the ultrashort to the narrowband 100Tc ≃ 628 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 pulse spectrum (violet shaded region) is ultrawideband in comparison to the material dispersion, whereas the ultrashort 3Tc 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.Fig. 4 Comparison of the spectral magnitude |ũ(ω)| of the gaussian pulse envelope for the ultrawideband 2τ0 ≃ 3.14 fs (magenta), the wideband 2τ0 ≃ 18.8 fs (red), and the narrow-band 2τ0 ≃ 628 fs (blue) pulse cases. The linear material phase dispersion (green curve) and ultrawideband spectrum |ω – ωc|–1 (dashed curve) plots are included for reference.
It is expected that the precursor fields will begin to become negligible in the total pulse evolution when the maximum slope of the initial pulse envelope function
is on the order of or greater than the medium relaxation time Tr ∼ δ−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 Nosc = 4.6 for ωc = 1 × 1015r/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 ωc in comparison to the medium resonance frequency ω0. On the other hand, as τ0 is increased above τc, the observed pulse distortion will approach that described by the group velocity approximation, as is now shown.A comparison of the computed linear and nonlinear pulse structure due to an input 3 oscillation gaussian pulse (2τ0 ≃ 18.8 fs) is presented in Fig. 5 at 5 absorption depths (z = 5zd) 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 ωc is sufficiently below the material resonance frequency ω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]; as the initial pulse width 2τ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.
Fig. 5 Comparison of the propagated linear (blue curve) and nonlinear (green curve) gaussian pulses at 5 absorption depths (z = 5zd) in a single resonance Lorentz-model dielectric with below resonance carrier frequency ωc = 0.416̄ω0.
Fig. 6 Comparison of the nonlinear gaussian pulse distortion using the full (blue curve) and quadratic approximation (green curve) of the linear material dispersion at 5 absorption depths for the Nosc = 0.5 gaussian envelope case (2τ0 ≃ 3.14 fs) with τ0/τc ≃ 0.11.
The sequence of figures presented in Figs. 6–12 provides a comparison of the nonlinear gaussian pulse evolution at 5 absorption depths as the initial pulse width 2τ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]. Notice that the group velocity approximation accurately describes just the trailing edge of the gaussian Brillouin precursor component, its peak amplitude being ∼ 10% of the actual peak amplitude value.
Fig. 7 Comparison of the nonlinear gaussian pulse distortion using the full (blue curve) and quadratic approximation (green curve) of the linear material dispersion at 5 absorption depths for the Nosc = 3 gaussian envelope case (2τ0 ≃ 18.8 fs) with τ0/τc ≃ 0.66.
Fig. 8 Comparison of the nonlinear gaussian pulse distortion using the full (blue curve) and quadratic approximation (green curve) of the linear material dispersion at 5 absorption depths for the Nosc = 5 gaussian envelope case (2τ0 ≃ 31.4 fs) with τ0/τc ≃ 1.10.
Fig. 9 Comparison of the nonlinear gaussian pulse distortion using the full (blue curve) and quadratic approximation (green curve) of the linear material dispersion at 5 absorption depths for the Nosc = 10 gaussian envelope case (2τ0 ≃ 62.8 fs) with τ0/τc ≃ 2.20.
Fig. 10 Comparison of the nonlinear gaussian pulse distortion using the full (blue curve) and quadratic approximation (green curve) of the linear material dispersion at 5 absorption depths for the Nosc = 15 gaussian envelope case (2τ0 ≃ 94.3 fs) with τ0/τc ≃ 3.30.
Fig. 11 Comparison of the nonlinear gaussian pulse distortion using the full (blue curve) and quadratic approximation (green curve) of the linear material dispersion at 5 absorption depths for the Nosc = 20 gaussian envelope case (2τ0 ≃ 125.7 fs) with τ0/τc ≃ 4.39.
Fig. 12 Comparison of the nonlinear gaussian pulse distortion using the full (blue curve) and quadratic approximation (green curve) of the linear material dispersion at 5 absorption depths for the Nosc = 100 gaussian envelope case (2τ0 ≃ 628.3 fs) with τ0/τc ≃ 22.0.
As the initial pulse width 2τ0 is increased with ωc 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 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 τc is exceeded, a gradual transition from the wideband, precursor dominated behavior to the narrowband quasimonochromatic behavior is observed, as evidenced in Figs. 9–12.
In Figs. 10–12 as the initial pulse spectrum becomes increasingly quasimonochromatic, the group velocity description provides a reasonable approximation to the actual pulse shape, the major two errors appearing in the underestimated pulse amplitude due to the overestimated material attenuation in the quadratic approximation and the slight phase shift appearing at the leading edge of the pulse. Finally, notice that special care had to be taken in modeling the 100 oscillation gaussian pulse case (Fig. 12) as high frequency numerical noise in the initial pulse spectrum was amplified by the cubic nonlinearity to produce a secondary pulse. This numerical error source was eliminated here by filtering the initial pulse spectrum to eliminate this noise, as illustrated in Fig. 13.
A comparison of the nonlinear space-time evolution of the envelope of the 3 oscillation gaussian pulse with initial pulse width 2τ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 zd, and in Figs. 16 and 17 for the mature pulse evolution [6, 8] as z increases above zd. 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 z/zd ∈ [0, 5] considered in the numerical study presented here. The 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/zd ∼ 0.5 to p ∼ 0.067 at z/zd ∼ 5. Again, this is due to the cubic nonlinearity which also significantly effects the peak amplitude decay in the group velocity approximation as z/zd increases above unity.
Fig. 14 Initial nonlinear space-time evolution of the envelope of a 3 oscillation gaussian pulse with initial pulse width 2τ0 ≃ 18.8 fs using the full material dispersion. The time axis spans the interval Δt = 125 fs.
Fig. 15 Initial nonlinear space-time evolution of the envelope of a 3 oscillation gaussian pulse with initial pulse width 2τ0 ≃ 18.8 fs using the quadratic approximation of the linear material dispersion. The spurious side lobes are due to the manner in which the envelope is numerically constructed from the pulse shape.
Fig. 16 Nonlinear space-time evolution of the envelope of a 3 oscillation gaussian pulse with initial pulse width 2τ0 ≃ 18.8 fs using the full material dispersion. The time axis spans the interval Δt = 125 fs.
Fig. 17 Nonlinear space-time evolution of the envelope of a 3 oscillation gaussian pulse with initial pulse width 2τ0 ≃ 18.8 fs using the quadratic approximation of the linear material dispersion. The time axis spans the interval Δt = 125 fs.
Fig. 18 Peak amplitude decay of a 3 oscillation gaussian pulse with initial pulse width 2τ0 ≃ 18.8 fs using the full material dispersion (blue curve) and using the quadratic approximation of the linear material dispersion (green curve). The black dashed line describes the Beer’s law exponential decay limit e−z/zd for comparison.
5. Other pulse shapes
Similar results are obtained for other canonical pulse shapes. As an example, the propagated field structure due to an input hyperbolic secant envelope pulse
is illustrated in Fig. 19 when τ0 = 3 × 10−15s and ωc = 1 × 1015r/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.Fig. 19 Comparison of the nonlinear hyperbolic secant pulse distortion using the full (blue curve) and quadratic approximation (green curve) of the linear material dispersion at 5 absorption depths for the Nosc = 3 envelope case (2τ0 ≃ 18.8 fs) with τ0/τc ≃ 0.66.
6. Conclusion
The detailed numerical results presented here have served to establish the following conclusions for the case of a causally dispersive absorptive medium with a nondispersive cubic nonlinearity:
- 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τc, where τ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], these precursor fields result in the propagated pulse breaking up into above resonance (Sommerfeld precursor) and below resonance (Brillouin precursor) sub-pulses when the initial pulse width is sufficiently small that 2τ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τc (the accuracy does not improve with the inclusion of higher-order terms);
- in addition, because the group velocity approximation overestimates the linear material attenuation away from the pulse carrier frequency ωc, the odd-order harmonics introduced by the cubic nonlinearity are absent in the group velocity description of the pulse evolution.
It is expected that the precursor fields will strongly influence the predicted pulse dispersion for other types of nonlinearity in the ultrashort/ultrawideband pulse regime. The intent of this paper is to show the central role that the precursor fields play in nonlinear pulse dynamics and thereby instigate a careful reformulation of nonlinear pulse propagation in the ultrashort pulse limit by other researchers in the nonlinear optics community.
Acknowledgments
The research presented here has been supported, in part, by the University of Vermont and by the United States Air Force Office of Scientific Research under AFOSR Grant # FA9550-08-1-0097.
References and links
1. A. Sommerfeld, “Über die fortpflanzung des lichtes in disperdierenden medien,” Ann. Phys. (Leipzig) 44, 4, 177–202 (1914).
2. L. Brillouin, “Über die fortpflanzung des licht in disperdierenden medien,” Ann. Phys. (Leipzig) 44, 4, 177–202 (1914).
3. L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).
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. 67, 535–558 (1909). [CrossRef]
5. J. A. Stratton, Electromagnetic Theory, (McGraw-Hill, New York, 1941), §5.18.
6. K. E. Oughstun and G. C. Sherman, Electromagnetic Pulse Propagation in Causal Dielectrics (Springer-Verlag, Berlin-Heidelberg, 1994).
7. K. E. Oughstun, Electromagnetic and Optical Pulse Propagation1: Spectral Representations in Temporally Dispersive Media (Springer, New York, 2006). [PubMed]
8. K. E. Oughstun, Electromagnetic and Optical Pulse Propagation2: Temporal Pulse Dynamics in Dispersive, Attenuative Media (Springer, New York, 2009).
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]
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]
12. T. H. Havelock, “The propagation of groups of waves in dispersive media,” Proc. R. Soc. A LXXXI, 398 (1908). [CrossRef]
13. T. H. Havelock, The Propagation of Disturbances in Dispersive Media (Cambridge U. Press, Cambridge, 1914).
14. L. Kelvin, “On the waves produced by a single impulse in water of any depth, or in a dispersive medium,” Proc. Roy. Soc . XLII, 80 (1887).
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]
17. G. P. Agrawal, Nonlinear Fiber Optics (Academic Press, San Diego, 1989).
18. R. H. Hardin and F. D. Tappert, SIAM Review Chronicle15, 2, 423 (1973).
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]
