Dynamic and inner-dressing control of four-wave mixing in periodically-driven atomic system
Li Yuan-Yuan, Li Li, Zhang Yun-Zhe, Zhang Lei
Institute of Physics, Xi’an University, Xi’an 710065, China

 

† Corresponding author. E-mail: liyynxcn@aliyun.com

Abstract

Four-wave-mixing (FWM) process is examined by using density matrix formalism in a periodically-driven atomic medium. Numerical result shows that FWM signals can be controlled by selecting different dynamic parameters of the probe field and strengths of the inner-dressing fields. It is also shown that the controllable FWM process is dominantly influenced by the evolution of atomic population difference and two-photon coherence. This dynamic and inner-dressing control of FWM is probably used for optimizing the optical nonlinear process and information processing.

1. Introduction

The suppression and enhancement of four-wave-mixing (FWM) process in electromagnetically induced transparency (EIT) windows[18] can be implemented in many multi-level systems, such as ladder-type,[2,9,10] Y-type,[8] and double lambda-type[46,1113] level configurations. In a ladder-type system with hyperfine ground states, the magnitude of FWM signal is found to be dependent on the transition route, and dominantly related to the residual two-photon coherence according to the degree of optical pumping to the other ground state.[9] While these features in the cycling transition are dramatically influenced by ladder-type and V-type two-photon coherence components, respectively.[9] Highly efficient FWM and six-wave mixing (SWM) signals coexisting in EIT windows were reported in an open four-level Y-type atomic system, in which the interplay between these two nonlinear processes is due to three- and five-photon destructive interferences.[8] Recently, in a backward double-lambda system, the conversion efficient in cold rubidium atoms was observed to be 63%, and the conversion efficient was predicted to be 96% by using a medium with a large optical depth.[13]

The control of FWM process can be achieved by changing several parameters of the fields applied to the system, e.g., the detuning and the strengths,[8,14] the relative phases,[15] and the pulse dynamics.[16] This control of FWM process can be applied to the quantum entanglement,[17] the fabrication of scalable multimode quantum resources,[18] and the enhancement of bright-seeded SU(1,1) interferometer.[19,20] In the present work, we consider a reverse Y-type atomic system driven by one probe and two coupling fields. The probe field is periodically modulated by a phase function. The two coupling fields are involved in both FWM and dressing process, thus we term the coupling fields inner-dressing fields hereafter. We will show that the FWM process can be successfully controlled by selecting different dynamic parameters of the probe field, and strengths of the inner-dressing fields. We will also show that the controllable FWM process is dominantly influenced by the evolution of atomic population difference and two-photon coherence. This type of dynamic and inner-dressing control of FWM can be probably used in optical information processing.

2. Theoretical model

A reverse Y-type four-level atomic system is shown in Fig. 1(a), where and are selected as hyperfine levels with a splitting difference Δ; and are the intermediate and the excited states of the atom, respectively. Beams diagram is shown in Fig. 1(b), where beam 1 represents a weak probe field , while beams 2 and propagating with a small angle difference (∼0.5°) represent strong inner-dressing fields and , respectively. The transition frequency of , , and are , , and respectively. The frequency, the detuning, the Rabi frequency and the wave vector of the field are denoted as , , gi , and respectively. The is assumed to be periodically modulated by a phase function , where A and ν are the modulation amplitude and frequency, respectively.[21] If the phase-matching condition is satisfied, the FWM signal propagating in the direction almost opposite to the beam can be generated.

Fig. 1. Diagram of (a) energy levels and (b) beam configuration in reverse Y-type four-level system.

The FWM signal is governed by the third-order density matrix elements and . Considering that exhibits almost the same feature as , we calculate only the density matrix element by using the following coupling equations:

where the complex detunings , , , , , with ( ); the decay rate from level to ( ) is defined as with being the decay rates of level ( ).

To obtain the above coupling equations, we assume that and . It is obvious that the modification of is mainly influenced by the evolution of population difference , two-photon coherences and . To solve these coupling equations, Fourier expansion of the field is usually used to find the resonant components involved in the coupling of relative transitions.[19] However, this method is cumbersome and incomplete for solving the dynamic and inner-dressing problem. In the following discussion, we numerically solve the above coupling equations by using the algorithm developed in Refs. [21] and [22]

To simplify the calculation, the atoms are assumed to be initially populated at the ground state , i.e., and ( ), one can obtain immediately . The coupling equations (3)–(6), and (8) are reduced to Eqs. (11)–(14), and (15) respectively:

The FWM signal is proportional to the third-order density matrix element , which can be obtained by numerically solving Eqs. (1), (2), and (11)–(15). In the following, we analyze these numerical results of FWM signals modified by dynamic parameters of the probe field, and the strength of the inner-dressing fields.

3. Numerical results and discussion

Parameters of 85Rb atom corresponding to possible experiments are chosen in our numerical calculation, where 5S1/2 (F = 2), 5S1/2 (F = 3), 5P1/2, and 5D1/2 levels are corresponding to states , , , and , respectively. These resonances at the structure are feasible to accomplish a chip-scale setup, which have a great practical importance. At the same time the harmonic phase modulation (frequency deviation) of miniature semiconductor laser (so-called VICSEL) is ordinarily used in practice.[23] The decay rates are , , , , and , respectively. It is assumed that the probe field is always very weak, with taking , the strengths of the inner-dressing fields and are much greater than that of the probe field. We also assume that the dressing fields are detuned exactly to transition, i.e., . Considering that satisfies (Fig. 1(a)), we are concerned about only the detuning in the following. A typical modulation frequency is usually , and the FWM lines versus around are examined in all cases.

We consider first the case that the strengths of the inner-dressing fields are far less than , with taking . In this case, both coherent population trapping (CPT)[24] in lambda-type sub-system ( ) and EIT in ladder-type sub-system ( ) contribute to the formation of a dark state, and FWM signals can transmit in the medium almost without any absorption. The FWM lines are shown in Fig. 2 for the different probe detuning values: , 0.48, 0.50, 0.52, and 0.54 (in units of , respectively, where the modulation index A = 1. The peaks of these lines are at the resonant position of , i.e., if is detuned away from transition with a value of , a probe field with a modulation frequency can introduce an enhanced FWM process at . When the modulation frequency ν is tuned to a fixed value , the peak values of these signals exhibit a damped oscillation behavior (Fig. 3(a)), which indicates that the FWM process can be optimized by appropriately choosing modulation frequency and indices. To verify this, FWM lines are plotted in Fig. 3(b) for A = 0.5, 0.8, 1, 1.95, 2.5, and 3, respectively. It is shown apparently that the FWM signals can be controlled by choosing different modulation parameters of the probe field.

Fig. 2. Plots of FWM lines versus modulation frequency ν near for different probe detuning values: , 0.48, 0.50, 0.52, and 0.54, respectively. Other parameters are A = 1, .
Fig. 3. (a) Modulation index-dependent peak amplitude of FWM lines at (b) FWM lines varying with ν for A = 0.5 (square), 0.8 (solid circle), 1 (upward triangle), 1.95 (downward triangle), 2.5 (rhombus), and 3 (leftward triangle), respectively. Other parameters are , .

Next, we examine the modulation dependent FWM lines for relative larger inner-dressing strength (Figs. 4(a)4(d)). In this case, the EIT in ladder-type sub-system ( ) dominantly contributes to the formation of a dark state, and FWM signals transmit in the EIT windows. Figure 4(a) and 4(b) show FWM lines for the inner-dressing strengths , 1, 1.5, 2, and 3 (in units), respectively, where the modulation index A = 1, and the probe detuning for Figs. 4(a) and 4(b) are and 0.48 (in units of ), respectively. It can be found that an increasing dressing intensity can suppress FWM signals at the resonant frequency of . As the dressing strength becomes greater than , a largely suppression of FWM signal happens at , and several enhanced peaks can be observed at (p = 1, 2, 3, …). As shown in Fig. 4(c), this type of suppression and enhancement of FWM signals can persist for different probe detuning values: , 0.48, 0.50, 0.52, and 0.54 (in units of ). One can find that the red detuning of the probe field can introduce a larger suppression of FWM signals when the resonant condition is satisfied. The modulation index A can also dramatically modify FWM signals in a strong coupling regime (Fig. 4(d)). The greater the peak amplitude, the larger the suppression at the resonant position will be. We believe that this suppression of FWM signals results dominantly from destructive interference between two channels induced by inner-dressing fields,[14] while the enhancement of FWM signals is induced by both dynamic parameters of the probe field and the strong inner-dressing effect.

Fig. 4. Modulation-dependent FWM lines for relative larger inner-dressing strength: (a) and (b) for different dressing strengths (square), 1 (solid circle), 1.5 (upward triangle), 2 (downward triangle), and 3 (rhombus), respectively; the probe detunings for panels (a) and (b) are and 0.48, respectively; (c) for different modulation frequencies (square), 0.48 (solid circle), 0.5 (upward triangle), 0.52 (downward triangle), and 0.54 (rhombus), respectively; and (d) for different modulation amplitudes A = 1 (square), 1.8 (solid circle), and 3 (upward triangle), respectively. Other parameters are taken as A = 1 in panels (a), (b) and (c), in panel (c), and in panel (d), respectively.

The controllable FWM signals can be explained by analysing the evolution of atomic level population difference and two-photon coherences. As stated, this evolution of FWM signals is mainly influenced by population difference , two-photon coherences and . However, is of third order and less than in our calculation, thus the modification of FWM process is influenced dominantly by and . An increasing coupling intensity will suppress the FWM signals at resonant position, while an enhanced contribution induced by the population change can counteract part of the suppression, e.g., the peak amplitude for is larger than that for .

As the dressing strength is greater than , is dominantly determined by two-photon coherence , since is far less than . In this case, the signal features can be explained by analysing the evolution of . The real (imaginary) part of is shown in Figs. 5(a) 5((b)) and 5(c) 5((d)) for different coupling strengths and modulation indices, respectively. Where the parameters used in Figs. 5(a) (5(b)) and 5(c) (5(d)) are corresponding to Figs. 4(a) and 4(d), respectively. At , the imaginary part of is zero, and the FWM process is governed mainly by the real part of , which gives rise to a largely suppression of FWM signals for an increasing dressing strength. As the modulation frequency is detuned away from the resonant position, the real part of this coherence decreases fast to zero, and then changes to a gain region with a damped oscillation behaviour, both the real part and imaginary part contribute to the splitting of peaks and the oscillation fringes. A similar evolution process can also be observed in the change of modulation index A as shown in Figs. 5(c) and 5(d). Thus, the control of FWM process can be achieved by tuning the dynamic parameters of the probe field and strengths of the inner-dressing fields.

Fig. 5. Plots of real and imaginary part of versus ν for different values of ((a), (b)) g2 and ((c), (d)) A. Parameters in panels (a)–(d) are corresponding to those in Figs. 4(a)4(d) respectively.

Here we should point out that this dynamic and inner-dressing controlling method suffers the competition between EIT and CPT induced by probe and coupling fields, and by the modulation of the probe field as well, respectively. It is different from the process in the out-dressing system investigated in Refs. [25]–[28] where the control process is achieved mainly by the modulation of electromagnetically induced gratings (EIG) formed by the dressing fields via switching between bright and dark states.

4. Conclusions

We theoretically investigate the FWM process in a periodically-driven atomic system. It is shown that an increasing dressing intensity can suppress FWM signal at the resonant position of . As the dressing strength becomes greater than , the almost totally suppression of FWM signals can happen at the resonant position, accompanied by several enhanced peaks at . This dynamic and inner-dressing control of FWM process is influenced dominantly by the evolution of atomic population difference and two-photon coherence, and can be probably used in optical information processing. By comparing with the advantages of the closed lambda-type system, N-type system and ladder-type system, the advantage of a reverse Y-type system in our discussion can be used to accomplish the coexistence of FWM and SWM processes since the former can be suppressed and the latter can be enhanced.[8,14]

Reference
[1] Fleischhauer M Imamoglu A Marangos J P 2005 Rev. Mod. Phys. 77 633 https://doi.org/10.1103/RevModPhys.77.633
[2] Harris S E Field J E Imamoglu A 1990 Phys. Rev. Lett. 64 1107 https://doi.org/10.1103/PhysRevLett.64.1107
[3] Petch J C Keitel C H Knight P L Marangos J P 1996 Phys. Rev. A 53 543 https://doi.org/10.1103/PhysRevA.53.543
[4] Li Y Xiao M 1996 Opt. Lett. 21 1064 https://doi.org/10.1364/OL.21.001064
[5] Zibrov A S Lukin M D Scully M O 1999 Phys. Rev. Lett. 83 4049 https://doi.org/10.1103/PhysRevLett.83.4049
[6] Wu Y Yang X X 2004 Phys. Rev. A 70 053818 https://doi.org/10.1103/PhysRevA.70.053818
[7] Li H J Huang G X 2007 Phys. Rev. A 76 043809 https://doi.org/10.1103/PhysRevA.76.043809
[8] Zhang Y Brown A W Xiao M 2007 Phys. Rev. Lett. 99 123603 https://doi.org/10.1103/PhysRevLett.99.123603
[9] Lee Y S Moon H S 2016 Opt. Express 24 10723 https://doi.org/10.1364/OE.24.010723
[10] Khadka U Zheng H Xiao M 2012 Opt. Express 20 6204 https://doi.org/10.1364/OE.20.006204
[11] Zhang Z Wen F Che J Zhang D Li C Zhang Y Xiao M 2015 Sci. Rep. 5 15058 https://doi.org/10.1038/srep15058
[12] Lee C Y Wu B H Wang G Chen Y F Chen Y C Yu I A 2016 Opt. Express 24 1008 https://doi.org/10.1364/OE.24.001008
[13] Liu Z Y Xiao J T Lin J K Wu J J Juo J Y Cheng C Y Chen Y F 2017 Sci. Rep. 7 15796 https://doi.org/10.1038/s41598-017-16062-5
[14] Li Y Li L Zhang Y 2010 J. Mod. Opt. 57 885 https://doi.org/10.1080/09500340.2010.495457
[15] Zhang H J Li X M Sun D Li H J Sun H 2018 Appl. Opt. 57 567 https://doi.org/10.1364/AO.57.000567
[16] Xue Y He Q Y LaRocca G C Artoni M Xu J H Gao J Y 2006 Phys. Rev. A 73 013816 https://doi.org/10.1103/PhysRevA.73.013816
[17] Zheng Z Wang H Cheng B Jing J 2017 Opt. Lett. 42 2754 https://doi.org/10.1364/OL.42.002754
[18] Wang H Fabre C Jing J 2017 Phys. Rev. A 95 051802 (R) https://doi.org/10.1103/PhysRevA.95.051802
[19] Xin J Wang H Jing J 2016 Appl. Phys. Lett. 109 051107 https://doi.org/10.1063/1.4960585
[20] Liu S Lou Y Xin J Jing J 2018 Phys. Rev. Appl. 10 064046 https://doi.org/10.1103/PhysRevApplied.10.064046
[21] Yudin V I Taichenachev A V Basalaev M Y 2016 Phys. Rev. A 93 013820 https://doi.org/10.1103/PhysRevA.93.013820
[22] Yudin V I Taichenachev A V Basalaev M Y Kovalenko D V 2017 Opt. Express 25 2742 https://doi.org/10.1364/OE.25.002742
[23] Post A B Jau Y Y Kuzma N N Happer W 2005 Phys. Rev. A 72 033417 https://doi.org/10.1103/PhysRevA.72.033417
[24] Hemmer P R Katz D P Donoghue J Cronin-Golomb M Shahriar M S Kumar P 1995 Opt. Lett. 20 982 https://doi.org/10.1364/OL.20.000982
[25] Du Y Zhang Y Zuo C Li C Nie Z Zheng H Shi M Wang R Song J Lu K Xiao M 2009 Phys. Rev. A 79 063839 https://doi.org/10.1103/PhysRevA.79.063839
[26] Zhang Y Wang Z Nie Z Li C Chen H Lu K Xiao M 2011 Phys. Rev. Lett. 106 093904 https://doi.org/10.1103/PhysRevLett.106.093904
[27] Wu Z Zhang Y Yuan C Wen F Zheng H Zhang Y 2013 Phys. Rev. A 88 063828 https://doi.org/10.1103/PhysRevA.88.063828
[28] Li C Jiang Z Zhang Y Zhang Z Wen F Chen H Zhang Y Xiao M 2017 Phys. Rev. Appl. 7 014023 https://doi.org/10.1103/PhysRevApplied.7.014023