Dynamics of short pulses and phase matched second harmonic generation in negative index materials

Abstract: We study pulsed second harmonic generation in metamaterials under conditions of significant absorption. Tuning the pump in the negative index range, a second harmonic signal is generated in the positive index region, such that the respective indices of refraction have the same magnitudes but opposite signs. This insures that a forward-propagating pump is exactly phase matched to the backward-propagating second harmonic signal. Using peak intensities of ~500 MW/cm, assuming χ~80pm/V, we predict conversion efficiencies of 12% and 0.2% for attenuation lengths of 50 and 5μm, respectively.

that occurs between counter-propagating beams.Also recently, Gorkunov at al have proposed the use of a doubly resonant system to induce perfect phase matching conditions, and order of magnitude enhancement of SHG [11].
SHG in NIMs has also been studied in cavity environments.It has been pointed out that SHG can occur for standard, centro-symmetric materials near the zero-n gap [12] of a Bragg grating made of alternating NIM/PIM layers [13].It turns out that a combination of field localization and coherent oscillations of the nonlinear dipoles located at the structure's interfaces leads to conversion efficiencies at least an order of magnitude larger than those achievable in the same length of a nonlinear, phase-matched, bulk material.In a related development, coupled cavities composed of nonlinear PIMs and linear NIMs [14] also display advantages over more conventional cavities because the NIM materials make better mirrors, leading to improved field localization effects and enhanced emission.
The overarching theme that characterizes all the works cited above is that SHG is efficient only if absorption inside the NIM is negligible.In fact, all linear and nonlinear theoretical approaches used to model NIMs have relied heavily on the assumption that absorption should be avoided.Although this constraint may not be crucial in cavity surroundings if NIM layers are relatively thin, the issue becomes critical in bulk environments.Naturally, this is a severe restriction that can raise serious questions in the minds of arbiters and ordinary readers alike, and may also lead one to question whether any proposed NIM-related phenomenon is observable in a bulk setting.Any predictions must in fact be ultimately reconciled with the fact that experimental observations [15] strongly suggest that incident waves are significantly attenuated in NIMs.For this reason, in this paper we focus on the proposed phase-matched scheme for backward second harmonic generation, similar to that discussed in references [8,9], and show that the process remains efficient even when absorption plays a significant role.Specifically, using the Drude model, the integration of Maxwell's equations in the time domain reveals that a medium having the following, generic characteristics: (i) a χ (2) ~ 80pm/V; (ii) pumped with pulses at least 2.5ps in duration, with peak intensities ~500MW/cm 2 ; is still capable of generating a backward-propagating (reflected) SH signal at a remarkable ~12% conversion efficiency rate even though the pump attenuation depth is roughly 50μm, and an equally impressive ~0.2% conversion rate for attenuation depths of 5μm.
We begin by writing nonlinear, second order polarization and magnetization as follows: (2) ( ; (2) P χ and (2) M χ are the electric and magnetic nonlinear coefficients, respectively.Assuming linearly polarized fields of the type:

ˆ( , ) . ( , ) . i kz t i kz t i kz t i kz t z t e c c z t e c c z t e c c z t e c c
the nonlinear polarization and magnetization may be written as: where the temporal derivatives arise from the linear dispersion of the medium, , and so on, and both ε and μ are complex functions of frequency (allowing for causal material dispersion) and of the spatial coordinate (allowing for spatial, material discontinuities).We have chosen r λ =1μm as the reference wavelength, and have adopted the following scaling: is the scaled frequency.Although they contain no approximations other than the assumption that the medium is isotropic, under most circumstances Eqs.( 4) may be simplified.For instance, in ordinary dielectric materials the dispersion length (associated with the second order temporal derivative on the left hand sides of Eqs.( 4)) varies from a few millimeters for few-cycle pulses, to a few meters for pulses several tens or a few hundred wave cycles in duration.NIMs, however, are more dispersive than ordinary materials, as described by either a Drude or a Lorentz model near a resonance condition.Nevertheless, typical dispersion lengths in regions of interest range from several hundred to several thousand wave cycles, depending on incident pulse width.In any case, the typical absorption length that we consider, which varies from 5μm to approximately 50μm, may be several orders of magnitude smaller than higher order dispersion lengths.This fact by itself obviates the need for the inclusion of higher order terms beyond the first order temporal derivative, which affects the group velocity of the pulse and should be preserved.As a result, Eqs.( 4) may be simplified to read: An inspection of both Eqs.( 4) and ( 5) suggests that they are symmetric with respect to the introduction of either a nonlinear polarization or magnetization.Therefore, one could just as easily treat a magnetic nonlinearity [19], and obtain qualitatively similar results.Moreover, while we do not dwell on the precise origin of the nonlinearity, one can reasonably assume that the nonlinearity resides in the nonlinear dipoles, and/or that the dipoles (both electric and magnetic) are embedded within a background medium that responds nonlinearly [19].With these considerations in mind, we assume .This choice results in an attenuation length of approximately 50μm, so that the intensity of the peak of the pulse drops to 1/e of its input value within that distance, and it is down to approximately 1% at 200 microns from the surface.One may compare this to the degree of absorption (i.e.imaginary part of the index) exhibited by GaAs at λ~890nm.
In the same figure we also report the general trend of the conversion efficiency as a function of normalized frequency (we do not scale it for simplicity), and we also identify pump and SH tuning ranges.To obtain the conversion efficiency curve shown we used incident pulses approximately 600fs in duration, or 200 wave cycles, and peak field values of ~10 4 Volts/m.We then calculated the total SH energy emitted in both directions, defined as ( ) ( , ) , where the integrand represents the local, instantaneous, energy density, normalized by the initial pump energy.We estimate that the second order dispersion length, defined as  .The phase matching condition between fundamental and SH pulses, which requires ( ) , is nearly exactly satisfied for the backward generated SH pulse [8,9].Neglecting the imaginary part, we have: In Fig. 2 we show the typical scenario that unfolds as a 200 wave-cycle, 500MW/cm 2 incident pulse (corresponding to a field ~3x10 7 V/m) crosses into a nonlinear NIM from vacuum.The pump is quickly attenuated, as a SH pulse is generated in the backward direction and rapidly exits the medium.The Drude dispersion reveals that the SH absorption length is approximately ten times larger compared to the pump absorption length, as can easily be deduced from the relative magnitudes of the part of the index.The study of pulsed dynamics reveals several noteworthy characteristics.For instance, in Fig. 3 we depict the conversion efficiency as a function of incident pulse duration, for two values of the damping coefficient.Given , the attenuation depth is drastically reduced to ~5μm.Quasi-monochromatic pulses, i.e. pulse duration greater than 2.5ps, yield conversion efficiencies of ~12% and ~0.2%, respectively.
In Fig. 4 we show the Fourier decomposition of the incident and scattered pump fields.We have three major components: (i) the input wave packet; (ii) a transmitted component, having negative wave vector; and (iii) a low intensity, reflected component.Reflection is small (less then 1 part in 10 7 ), because we have assumed that ε=μ, a restriction that can easily be removed.
In Fig. 5 we show the Fourier components of the SH signal.Here, too, we find three major components: (i) the one generated backward, that has already exited into free space, at twice the incident pump wave-vector; (ii) the phase matched component, traveling backward inside the medium toward the entrance interface, with wave vector almost precisely at twice the pump wave vector; (iii) a forward-moving, SH component, generated far from the phase matching condition, which is being absorbed away.The conversion efficiencies that we report are large by any measure, considering the amount of absorption we have introduced.Even though the magnitude of the nonlinear coefficient may be somewhat large, all things being equal, i.e. χ (2) and incident peak power, it would be possible to obtain similar conversion efficiencies if, in the absence of absorption, two conditions could be satisfied: (i) exact phase and (ii) group velocity matching.For comparison purposes, if a 200fs, 500MW/cm 2 pulse were incident on a semi-infinite, ideally matched medium having n~1.42, one obtains a conversion efficiency of η~10% after the peak of the pulse propagates ~45μm.In contrast, the introduction of an index mismatch due to normal material dispersion (we consider a typical index variation of ~6% between ω (n~1.42)and2ω (n~1.52)),along with an appropriate group velocity mismatch ), reduces the conversion efficiency dramatically down to η~0.01%, with further decreases if absorption is added.We should note that a group velocity mismatch has a much less pronounced impact on conversion efficiency compared to an index mismatch, as expected, but it can make its impact felt, particularly for shorter pulses.However, it is not feasible to simultaneously have phase and group velocity matching in a bulk medium, as natural material dispersion generally does not allow it.Indeed it is possible to have effective phase matching and achieve approximately similar group velocities for the fundamental and SH waves in photonic band gap structures, where geometrical dispersion can be used to compensate the naturally occurring material dispersion [20,23].Nevertheless, absorption usually has a detrimental effect that significantly reduces conversion efficiency under any circumstances where the beams co-propagate, including resonant conditions.In contrast, the unusual phase matching conditions that naturally occur in a NIM for counter-propagating waves also suffer the effects of absorption, but not to the same detrimental degree because the rapid decay of the pump along the positive direction is accompanied by an equally swift exponential increase of the SH signal toward the nearby entry surface.In addition to nearly exact phase matching conditions, the fields overlap for the entire duration of the pump pulse, and the impact of absorption is minimized.In conclusion, we have shown that second harmonic generation can occur in a NIM and be relatively efficient in the presence of significant absorption.Assuming a Drude dispersive model, we have shown that a phase matching condition established between counterpropagating fundamental and second harmonic waves helps boost conversion efficiencies up to 12% for the reflected SH component when the attenuation depth of the NIM is approximately 50μm, χ (2) ~80pm/V, and peak pump intensity is ~500MW/cm 2 .The interaction is consumed entirely near the surface interface, where both fields are more intense.In the absence of absorption, similar conversion efficiencies are possible in bulk PIMs, only by imposing simultaneous phase and group velocity matching.Finally, our calculations show that conversion efficiencies of order 0.2% are also possible, even when the attenuation depth is approximately 5μm, and suggest that similarly efficient higher harmonic generation and parametric amplification may be also achievable in a multi-wave mixing environment under conditions of significant absorption.
(the plasma frequency corresponds to the reference frequency); the damping coefficient (also scaled by the reference frequency

2 ( 2 Fig. 1 .
Fig. 1.Real (blue curve-left axis) and imaginary (red curve-right axis) parts of the index of refraction for the Drude model, with

Figure 1
Figure1features the larger point we wish to make: at

Fig. 2 .
Fig. 2. A 600fs (200 wave cycles) pump (blue curve) is incident on a NIM (shaded region) from vacuum, and tuned to

Fig. 4 .
Fig. 4. Fourier decomposition of the incident ( k ω ) and scattered pump pulses.The wave number is dimensionless.The transmitted wave packet ( ( ))k nω ω is a superposition of negative wave vectors, even though the pulse propagates in the forward direction.There is a small reflected component propagating backward in free space ( k ω − ).

Fig. 5 .
Fig. 5. Fourier decomposition of the generated SH field.The wave number is dimensionless, as in Fig. 4. The positive component identifies a mismatched signal propagating forward inside the medium ( 2 (2 ) k n ω