We study the case of two simultaneous three-wave-mixing processes, where one frequency is converted to another through an intermediate frequency. The common assumption is that these processes can occur only when the material is transparent at all participating frequencies. Here we show experimentally that, under appropriate conditions, the intermediate frequency remains dark throughout the interaction. This means that even if the material is opaque at the intermediate frequency, the conversion will remain efficient. New possibilities of frequency conversion are therefore available, e.g. through absorptive bands in the ultraviolet or mid-infrared. Moreover, though it was hitherto assumed that the phase mismatch value is governed only by dispersion, we show here that phase matching also depends on light intensity. These findings promise novel all optical switching techniques.
© 2012 OSA
In the field of nonlinear optics, available lasers are mixed to produce new frequencies, using three-wave-mixing processes in quadratic nonlinear optical crystals. For example, in sum frequency generation (SFG), two frequencies ω1 and ωp1 are summed to produce a higher frequency ω2 = ω1 + ωp1. In order for this process to have significant efficiency, momentum conservation needs to be maintained. For the given example, this would require that k⃗2 = k⃗1 + k⃗p1 where k⃗i is the momentum (wave-vector) of the wave with frequency ωi, a requirement also known as phase-matching. Three-wave-mixing processes can also be cascaded , e.g. the generated sum-frequency can be further combined with another frequency in another SFG process, which produces ω3 = ω2 + ωp2 = ω1 + ωp1 + ωp2.
Recently, a wideband and robust frequency conversion method was developed, by using an analogy between the dynamics of sum or difference frequency generation (DFG) process and a two state quantum system induced by electromagnetic (EM) field [2, 3].
Here, a method to perform frequency conversion via simultaneous three-wave-mixing (STWM) is developed and demonstrated experimentally. Most notably, phase-matching is found to be dependent on light intensity. This is in contrast to the hitherto accepted view that phase-matching depends only on the dispersion properties of the nonlinear medium . Moreover, the intermediate frequency remains dark, hence conversion can be preformed through frequencies at which the nonlinear medium is opaque. The method presented here relies on an analogy drawn between an EM-excited three state quantum system and two STWM processes. This analogy leads to mapping of atomic phenomena such as dark states and Stark shifts to the field of nonlinear optics.
2. Theoretical analysis
2.1. Dynamical equations and atomic analogy
We consider two simultaneous three wave mixing processes (STWM), where the generated frequency is further combined with another frequency, ωp2, in an additional process, generating a new frequency, ω3. For the case of two SFG processes, we obtain ω3 = ω2 + ωp2 = ω1 + ωp1 + ωp2. Generally, two distinct pumps can drive the interaction: the first transfers energy between ω1 and ω2, while the second connects ω2 with ω3. We consider the case where the pumps are much more intense than the other waves and thereby are negligibly affected by the interaction (undepleted pumps approximation). The procedure of adiabatic elimination [5, 6] is used to reduce the three wave system into an effective two coupled waves system. Notably, a work conducted in parallel to ours applied this procedure to the case termed ’direct third harmonic generation’ . Notwithstanding, here we analyze STWM by analogy to a three level quantum system, following the approach of Suchowski et al. [2, 3].
Let us start by writing the coupled equations that govern the evolution of the two STWM processes along a nonlinear crystal:
The same set of equations describes the dynamics of a three level atom interacting with two EM fields [5, 8, 9]. In the context of atomic physics, the complex amplitude of each field Ai (z) corresponds to ai (t), the probability amplitude of population in each state. Each effective coupling coefficient σij represents the strength of the dipole interaction between levels i and j, usually described by the Rabi frequency Ωij = dij · ε(t)/h̄, where dij is the dipole moment between states i and j, ε(t) is the induced EM field, and h̄ is Planck’s constant. Respectively, the terms Δk1 and Δk2 correspond to the detunings, Δ1 and Δ2, with the two atomic transitions.
We note that, for the optical case, the notation used here assumes quasi-monochromatic laser beams, which in the visible and near infrared means pulses longer than about 1ps . However, since the interaction is coherent, femtosecond pulses can also be used. Transform limited femtosecond pulses can be stretched before entering the crystal, and the generated frequency pulse can be compressed to a transform limited duration .
2.2. Adiabatic elimination
Utilizing the adiabatic elimination procedure is done by assuming that each of the SFG or DFG processes exhibits a very large phase mismatch, but their sum is rather small, i.e. |Δk1| L ≫ 1, |Δk2| L ≫ 1 while |Δk1 + Δk2| L ≪ 1, where L is the length of the nonlinear crystal. The sum of the phase-mismatches, which is a central notion here, will be defined as the two-photon phase mismatch, and denoted ΔkTP = Δk1 + Δk2. Under these conditions, it is reasonable to assume that the oscillating terms in Eq. (1) vary much faster than the fields amplitudes at ω1 and ω3, allowing us to write an approximation for the intermediate frequency complex amplitude:
By substituting the outcome of Eq. (2) into Eq. (1) and moving to the rotating frame using A1 = Ã1 (z) exp [−i(Δkeff/2 – σ12σ21/Δk1)z] and A3 = Ã3 (z) exp [i(Δkeff/2 – σ23σ32/Δk2)z], a reduced two coupled wave equations set for ω1 and ω3 is obtained,
The physical interpretation of Eqs. (2) and (3) is as follows: since each of the two STWM processes is greatly phase-mismatched, the amplitude of the intermediate frequency, |A2|, will oscillate very rapidly in comparison to variation in the amplitudes of ω1 and ω3. As a result, |A2| will never build up to a significant value and thus remain low throughout the interaction. Nevertheless, energy can still be transferred between ω1 and ω3 via ω2, effectively imitating a four-wave-mixing process with undepleted pumps, where ω3 = ω1 + ωp1 + ωp2 in the case of two simultaneous SFG processes. We emphasize that the efficiency of conversion from ω1 to ω3 does not change even if the material is opaque at ω2.
The rate of energy transfer and the amount of energy participating in this process are determined by the effective phase-mismatch Δkeff. In this mechanism, two additional contributions to the conventional phase mismatch parameters appear. These contributions, which depend on the coupling coefficients (and thus on the intensities of the pumps, Ip1 and Ip2), are the equivalent of the atomic Stark shift . In the nonlinear optics context, we choose to define it as δkS = σ12σ21/Δk1 + σ23σ32/Δk2.
This term is critical where significant energy transfer between ω1 and ω3 is desired, which occurs when the effective phase-matching parameter is equal to zero, i.e. when Δkeff = Δk1 + Δk2 +δkS = 0. Since δkS ∝ Ip it introduces a pump intensity dependence into the phase-matching condition. This is contrary to the conventional wave-mixing phase mismatch, which is considered to depend solely on the dispersion properties of the nonlinear material. Such dependence exists since the rate of each of the two STWM processes depends on pump intensity. To obtain a net flow of energy from one frequency to another, these rates must be compatible over a significant interaction length. The mechanism here is thus fundamentally different from phase-modulation in four-wave-mixing processes, which stems from intensity dependent refractive index.
The analytical solution of Eq. (3), for the initial condition A3 (0) = 0, is11] and Feynman et al. .
The effective phase-matching condition can be satisfied using the quasi-phase-matching (QPM) technique . For a periodically poled crystal with period Λ, the phase-mismatch of each of the two STWM processes should be modified with the addition of the periodicity pattern, i.e. in the first order Fourier approximation, Δk̃1 = Δk1 + 2π/Λ and Δk̃2 = Δk2 + 2π/Λ. Substituting these expressions into the phase-matching condition produces a cubic equation with one solution that adheres to the conditions of adiabatic elimination.
3. Numerical simulation
A QPM solution as described above was used to numerically simulate the interaction. A single pump at λp1 = λp2 = 1064nm was used to drive both processes. We have used a periodically poled KTiOPO4 (PPKTP) crystal with a poling period of Λ = 8.6μm, which was calculated to phase-match conversion from λ1 = 3010nm to λ3 = 452nm, via the two implicit SFG processes 3010nm + 1064nm → 786nm and 786nm + 1064nm → 452nm, when the crystal temperature is 125°C. The simulation does not assume undepleted pump or adiabatic elimination.
The resulting evolution of the intensities of the interacting waves, with pump intensity of Ip = 10GW/cm2, is plotted in Fig. 1(a). Energy oscillations between the λ1 = 3010nm input and the λ3 = 452nm output are clearly seen. Since each photon converted from λ1 to λ3 is paired with two pump photons, the peak output intensity at λ3 is higher than that at λ1. The inset shows the fast oscillations of the low intensity at the adiabatically eliminated intermediate wavelength λ2 = 786nm, which are over 3 orders of magnitude lower than the peak output intensity.
Using the same parameters, the Stark effect in frequency conversion was analyzed. In Fig. 1(b) the photon conversion efficiency from λ1 to λ3 is calculated for two cases: low pump intensity of 1GW/cm2 and high pump intensity of 12.4GW/cm2. The results are plotted as a function of the two-photon phase-mismatch ΔkTP = Δk1 + Δk2. For low pump intensity the Stark shift is negligible, so the peak efficiency is obtained near the two-photon resonance condition ΔkTP = 0. However, when the pump intensity is high, the Stark shift can no longer be neglected and the efficiency peak shifts by δkS to ΔkTP = −δkS, which is the exact shift amount predicted analytically using Eq. (4).
Following the numerical simulation, STWM with adiabatic elimination had been experimentally realized. The experimental system is depicted in Fig. 2. A 4.4nsec Nd:YAG laser beam with λp = 1064nm, which served as the strong pumps, was passed through a temperature controlled PPKTP crystal with a Λ = 8.6μm period. A laser with λ1 = 3010nm and an estimated 10nm FWHM bandwidth, which was produced from the same pump through an optical parametric oscillator, served as our input. The emerging power at the λ3 = 452nm output and the intermediate wavelength λ2 = 786nm were independently measured, using power detectors or a spectrometer, and peaks of both wavelengths were observed, corresponding to the expected output and intermediate wavelengths. When either the pump or the input were prevented from entering the crystal, both peaks would vanish.
Next, we have chosen to examine the conversion efficiency of the process as a function of the two-photon phase mismatch. This was done by tuning the crystal temperature from 90°C to 175°C, affecting the refractive index , and thus the phase-mismatch parameter. The input and pump intensities were kept constant at 71kW/cm2 and 44.6MW/cm2, respectively. Figure 3(a) shows the experimentally measured 452nm output power vs. the PPKTP crystal temperature, after decrementing the background noise, and accounting for the transmission and reflection of the optical components in the apparatus. The detected noise, defined as the power measured when only the pump beam was incident on the PPKTP crystal, was found to be around 152nW and thus negligible compared to the measured output power which reaches 38μW. As seen, optimum two-photon phase matching was obtained at crystal temperature of 133°C. The small deviation from the expected value of 125°C can be accounted for by the inaccuracy in the Sellemeier equations that were used to calculate the wavelength and temperature dependence of the crystal’s refractive index [14, 15]. The power at λ2 was below the noise level and therefore could not be measured. Its upper bound was found to be 51.8nW, which is 733 times lower than the peak 452nm output power, agreeing with the theoretical prediction of adiabatic elimination.
Finally, the predicted dependence of the 452nm output power on pump power was experimentally checked. A half wave plate was placed in the path of the pump beam before it entered the PPKTP crystal, and gradually rotated through 45 degrees. Due to the fact that KTiOPO4 is a biaxial birefringent crystal, its phase-matching condition is polarization-dependent: only the pump power polarized along the PPKTP c axis contributes to the STWM processes. Figure 3(b) shows the experimentally measured 452nm output power obtained with various values of c-polarized pump power. In this experiment the noise (defined above) was measured at each wave plate angle and decremented from the detected power. The experimental measurements fit well on the analytically calculated output power (using Eq. (5)), showing a quadratic nature, where the nonlinear coefficient χ(2) was taken to be 9% higher than the value reported for other processes and wavelengths . This is consistent with the low conversion efficiency limit of Eq. (5), . This experimental observation is a clear indication that STWM is in action, since a single TWM process would result in a linear rather than quadratic dependence, while direct third harmonic generation would yield a cubic dependence.
We have experimentally demonstrated frequency conversion through an intermediate frequency which remains dark throughout the interaction, enabling conversion through frequencies at which the nonlinear medium is absorptive. Additionally, phase-matching was shown to depend not only on dispersion, but also on pump intensity. This phenomenon enables new methods of all optical switching. We presented an analogy with EM-excited three state quantum systems which lead to these new effects. The same mathematical treatment can also be applied to cases involving three or more STWM processes, analogous to multi-state quantum systems . Other atomic physics schemes, involving temporal variation of the induced EM fields [5, 18], can be transferred to frequency conversion using aperiodically poled structures.
This work was supported by the Israel Science Foundation.
References and links
1. S. M. Saltiel, A. A. Sukhorukov, and Y. S. Kivshar, “Multistep parametric processes in nonlinear optics,” Prog. Opt. 47, 1–73 (2005). [CrossRef]
2. H. Suchowski, D. Oron, A. Arie, and Y. Silberberg, “Geometrical representation of sum frequency generation and adiabatic frequency conversion,” Phys. Rev. A 78, 063821 (2008). [CrossRef]
4. R. W. Boyd, Nonlinear Optics, 3rd ed.(Academic Press, 2008).
5. D. Tannor, Introduction to Quantum Mechanics: A Time-Dependent Perspective (University Science Books, 2007).
6. D. Grischkowsky, M. M. T. Loy, and P. F. Liao, “Adiabatic following model for two-photon transitions: nonlinear mixing and pulse propagation,” Phys. Rev. A 12, 2514–2533 (1975). [CrossRef]
7. J. C. Delagnes and L. Canioni, “Third harmonic generation in periodically poled crystals,” Proc. SPIE 7917, 79171C (2011). [CrossRef]
8. C. Trallero-Herrero and T. C. Weinacht, “Transition from weak- to strong-field coherent control,” Phys. Rev. A 75, 063401 (2007). [CrossRef]
9. F. T. Hioe, “Dynamic symmetries in quantum electronics,” Phys. Rev. A 28, 879–886 (1983). [CrossRef]
10. H. Suchowski, B. D. Bruner, A. Ganany-Padowicz, I. Juwiler, A. Arie, and Y. Silberberg, “Adiabatic frequency conversion of ultrafast pulses,” Appl. Phys. B 105, 697–702 (2011). [CrossRef]
11. F. Bloch, “Nuclear induction,” Phys. Rev. 70, 460–474 (1946). [CrossRef]
12. R. P. Feynman, F. L. Vernon, and R. W. Hellwarth, “Geometrical representation of the Schrödinger equation for solving maser problems,” J. Appl. Phys. 28, 49–52 (1957). [CrossRef]
13. D. S. Hum and M. M. Fejer, “Quasi-phasematching,” C. R. Physique 8, 180–198 (2007). [CrossRef]
15. K. Fradkin, A. Arie, A. Skliar, and G. Rosenman, “Tunable midinfrared source by difference frequency generation in bulk periodically poled KTiOPO4,” Appl. Phys. Lett. 74, 914–916 (1999). [CrossRef]
16. I. Shoji, T. Kondo, A. Kitamoto, M. Shirane, and R. Ito, “Absolute scale of second-order nonlinear-optical coefficients,” J. Opt. Soc. Am. B 14, 2268–2294 (1997). [CrossRef]
17. N. V. Vitanov, B. W. Shore, and K. Bergmann, “Adiabatic population transfer in multistate chains via dressed intermediate states,” Eur. Phys. J. D 4, 15–29 (1998). [CrossRef]
18. J. Keeler, Understanding NMR Spectroscopy (John Wiley & Sons, 2005).