Raman-Kerr Comb Generation Based on Parametric Wave Mixing in Strongly Driven Raman Molecular Gas Medium
Aurélien Benoît, Anton Husakou, Benoît Beaudou, Benoît Debord, Frédéric Gérôme, Fetah Benabid
RRaman-Kerr Comb Generation Based on Parametric Wave Mixing in Strongly Driven Raman Molecular Gas Medium
Aurélien Benoît, Anton Husakou, Benoît Beaudou, Benoît Debord,
Frédéric Gérôme, and Fetah Benabid GPPMM Group, Xlim Research Institute, CNRS UMR 7252, University of Limoges, Limoges, France College of Max Born Institute, Forschungsverbund Berlin e.V.RudowerChaussee 1712489 Berlin, Germany GLOphotonics SAS, 123 avenue Albert Thomas 87060 Limoges Cedex, France * Corresponding author: [email protected]
We report on experimental and theoretical demonstration of an optical comb spectrum based on a combination of cascaded stimulated Raman scattering and four-wave mixing mediated by Raman-induced nonresonant Kerr-type nonlinearity. This combination enabled to transform a conventional quasi-periodic Raman comb into a comb with a single and smaller frequency spacing. This new phenomenon is realized using a hollow-core photonic crystal fiber filled with 40 bars of deuterium, and pumped with a high-power picosecond laser. The resultant comb shows more than 100 spectral lines spanning over 220 THz from 800 nm to 1710 nm, with a total output power of 7.1 W. In contrast to a pure Raman comb, a 120 THz wide portion of the spectrum exhibits denser and equally-spaced spectral lines with a frequency spacing of around 1.75 THz, which is much smaller than the lowest frequency of the three excited deuterium Raman resonances. A numerical solution of the generalized nonlinear Schrödinger equation in the slowly varying envelope approximation provides very good agreement with the experimental data. The additional sidebands are explained by cascaded four-wave mixing between pre-existing spectral lines, mediated by the large Raman induced optical nonlinearity. The results open a new route to the generation of optical frequency combs that combine large bandwidth and high power controllable frequency spacing. I. I NTRODUCTION
The generation of optical combs is turning into an extremely important topic in several scientific and technological fields, such as waveform synthesis [1], the optical clock [2], communication systems [3] or biomolecular imaging [4]. Historically, generation of optical combs was associated with mode-locked lasers [5] or with high harmonic generation using high-peak-power and ultra-short pulse laser sources [6-7]. The former has matured to be the cornerstone in optical frequency metrology thanks to the outstanding frequency stability of the comb spectral components [8], and the latter is a powerhouse in synthesizing attosecond pulse in VUV and Soft X-ray thanks to its multi-octave bandwidth [9]. During the last two decades, photonics’ solutions to generate optical combs have emerged with potential to address current limitations such as compactness or multi-octave comb generation in more convenient spectral ranges. Among these, we list Kerr optical nonlinearity based combs, which are generated by micro-resonators, and pumped with a continuous-wave laser [10]. Here, the generated spectrum can be up to one octave wide and the frequency spacing of its components, set by the micro-resonator diameter, ranges from a few hundred GHz up to 1 THz. Another means of generating discrete spectrum is based on stimulated Raman scattering (SRS) in gas-filled hollow-core photonics-crystal-fiber (HC-PCF) [11]. Compared to mode-locked laser or micro-resonator mechanisms, the Raman-gas-filled HC-PCF based combs stand out with their large bandwidth [11], which can span more than five octaves from the ultraviolet to mid-infrared range, and with their high optical power and energy handling [12]. The generated spectra using Raman-gas-filled HC-PCF with nanosecond or picosecond pump pulses show intra-pulse phase coherence [13] and are typically composed of lines with several frequency spacings exclusively set by the rotational and/or vibrational Raman resonances [11-12, 14-15]. Consequently, the typical spectral structure of a HC-PCF based Raman comb consists of a central frequency ν set by the pump laser, and a cascade of sidebands with frequencies in the form of ν l,m,n = ν + lν R1 + mν R2 + nν R3 for the case of a medium with three Raman transitions. Here, l , m , n are integers corresponding to the negative (Stokes) and positive (anti-Stokes) orders of the Raman transitions with frequencies ν R1 , ν R2 and ν R3 respectively. Note that the spectral structure of such SRS comb is set by resonant third-order susceptibility in the vicinity of the Raman resonance [16], with no measurable influence from the non-resonant non-linear susceptibility ( i.e. Kerr effect) [17]. As such, the HC-PCF based Raman comb suffers from two major drawbacks that hinder their use in a number of applications. It exhibits multiple frequency spacing (quasi-periodic comb), and its Raman frequency is too large (> 10 THz for most of the representative molecular gases). In this paper, we show a new approach to overcome this limit of Raman combs and to generate a small and single frequency-spacing comb in Raman gas media of the form ν q = ν + q∆ . The approach relies on using a Raman gas containing at least two nearly commensurable resonance frequencies, so that, for example, ν R2 and ν R3 can be written as ν R2 = N∆ + δ and ν R3 = M∆ + δ' with N and M being positive integers, δ,δ' ≪ ∆. The proof of principle of the approach is undertaken by experimentally and theoretically demonstrating the generation of a wide comb in a Raman gas-filled HC-PCF that exhibits spectral lines with a frequency-spacing as low as 1.75 THz. The observed spectrum, which is generated from a high-pressure deuterium molecular gas, consists of up to 100 spectral lines overing a width of 220 THz. In addition to lines associated with the three excited Raman resonances of deuterium at frequencies of 89.6 THz, 12.4 THz and 5.3 THz respectively, a good part of spectrum includes a denser comb of 52 lines covering a range of 120 THz with a frequency spacing of 1.75 THz. We attribute these denser spectral lines to an action of Raman-induced parametric four wave-mixing (FWM) ( i.e. the out-of-resonance contribution of the Raman susceptibility), and to a FWM process between close Raman sidebands to generate frequency components outside the Raman resonance frequencies. II. E XPERIMENTAL SETUP AND R ESULTS
The experimental set-up consists of a high-power picosecond pulsed Yb-fiber laser coupled into a photonic microcell (PMC) comprising a length of HC-PCF whose ends are attached to gas cells [18]. The laser emits at 1030 nm, has an average power of 16 W, a repetition rate of 1 MHz, pulse duration of 27 ps and a self-phase modulation broadened spectrum with a bandwidth of 200 GHz [19]. The PMC is filled with deuterium molecules (D ) at a pressure of 40 bar. The choice of D is motivated by the frequency values of its two ro-vibrational Raman transitions. At room temperature, D gas is a 2:1 molecular mixture of ortho- deuterium (o-D ) and para-deuterium (p-D ). It exhibits three Raman resonances: a vibrational one around ν R1 = 89.6 THz with a steady-state Raman gain of 0.37 cm/GW, and two ro-vibrational ones at ν R2 = 12.44 THz from o-D [20], and at ν R3 = 5.33 THz from p-D [21] with gains of ~0.12 cm/GW and ~0.05 cm/GW respectively. As a result, the deuterium Raman shifts for para- and ortho- molecules are close to be commensurable, with their ratio being equal to 3:7 within less than 0.1% of relative error. As such, we have 5.33 THz ~ 3Δ, and 12.44 THz ~ 7Δ, where Δ = ~ 1.78 THz coincides with the | ν R2 - 2 ν R3 | value, which is also the frequency spacing between the Raman lines ν l,m+1,-1 and ν l,m,1 , and is three times smaller than the smallest Raman frequency. The impact of this peculiarity will be apparent later. Furthermore, the gains of the o-D and p-D ro-vibrational resonances have closer values (2:1 ortho- to para- ratio of gains) than for the hydrogen gas (3:1 ratio) [22], favoring thus the excitation of p-D Raman resonance when compared to H . Finally, operating at a pressure as high as 40 bar ensures higher Raman gain but broader gain linewidth, typically ~ 10 GHz. The latter is notably useful to fulfilling the above-mentioned commensurability. The PMC fiber consists of a 3 m long inhibited-coupling single-ring tubular-lattice HC-PCF [23] with an inner diameter of 40 µm [Fig. 1(a)]. Fig. 1. (a) Scanning electron microscope of the inhibited-coupling single-ring tubular-lattice HC-PCF used for this experiment. (b) Evolution of losses (black curve) and effective refractive index in the vacuum (blue curve) and hydrogen (red curve) as a function of wavelength. (c) Evolution of group index as a function of the wavelength in the vacuum (blue curve) and in hydrogen (red curve). (d) Evolution of dispersion as a function of wavelength in the vacuum (blue curve) and in hydrogen (red curve). (e) Output spectrum after a 16 W power pump propagation in 3 m long HC-PCF filled with 40 bars of deuterium. The pump laser line is identified by black color, the vibrational Stokes and anti-Stokes by purple, the Stokes and anti-Stokes lines due to ortho-ro-vibrational Raman by blue, the Stokes and anti-Stokes lines due to para-ro-vibrational Raman by green, and four-wave-mixing generated sidebands by red. Inset: near field intensity distribution of the output beam. (f) Zoomed-in spectrum of the black-dashed rectangle. The fiber transmission loss in the visible-infrared range is characterized by an ultrawide (over one octave) band from 600 to 1200 nm with loss around 11 dB/km [Fig. 1(b)]. Figure. 1(b) presents also the evolution of the refractive index as a function of the wavelength. Figs. 1(c-d) show the used fiber spectra of the group index and the dispersion respectively. A detailed account of all the optical properties of the HC-PCF is given in [23]. The laser beam is coupled into the PMC in a manner to match its beam size with the fiber fundamental core mode. Inspection of the fiber output beam profile over the operating spectral range of 800-1600 nm shows single mode propagation [inset of Fig. 1(e)]. Figure 1(e) presents the typical PMC output spectrum hen excited with 16 W average power pump laser. The power of the generated and transmitted spectrum was measured to be 7.1 W. The spectrum spans a wide range from 170 to 390 THz and exhibits a large number of lines in a comb-like structure with more than 100 lines. The spectral components reveal that two purple-colored lines correspond in their frequency shifts (89.63 THz) to the vibrational Stokes ν -1,0,0 and anti-Stokes ν lines of the pump. The vibrational comb lines are accompanied by the comb of blue lines, which appear around both pump and the 1 st order vibrational Stokes and anti-Stokes lines. The blue lines are readily associated with ortho-ro-vibrational Raman resonances, as corroborated with their frequency offset from their associated pump (average experimental shift of 12.45 ± 0.06 THz as compared to the theoretical value of 12.44 THz). Similarly, the green lines correspond to the Stokes and anti-Stokes sidebands of either the purple or the blue lines, and are generated by the para-ro-vibrational Raman resonance (average experimental shift of 5.30 ± 0.20 THz as compared to the tabulated value of 5.33 THz). Finally, the spectrum includes red lines, which, even though they can be described in terms of the Raman frequency expression ν l,m,n , are generated by more subtle mechanism than the simple cascaded Raman picture effect, as discussed below. Limiting the further discussion to only the pump laser sidebands of orders (0, m , n ), we see that this part of the comb extends over 220-340 THz, and its Raman frequencies take the form of ν with no overlapping lines from the vibrational Stokes and anti-Stokes sidebands. Remarkably, the spectrum stands out with equally spaced spectral components, with a frequency-spacing Δ = 1.75 ± 0.05 THz, which is also equal to ν R2 - 2 ν R3 , the spacing of two Raman lines located at ν and ν . Thus, this part of the spectrum is transformed from a quasi-periodic Raman comb into a single-frequency comb whose component takes the form ν ̅ q = ν + q∆ , with q being an integer. Figure. 1(f) shows a close-up of a portion of this spectrum, which is also representative of any range between ν and ν . It shows that in addition to Raman blue lines ( ν and ν ), and green lines ( ν and ν ), we have four red lines (labeled ν FWM1 , ν FWM2 , ν FWM3 , ν FWM4 ). We start by recalling that the generation of the comb components is chronologically ordered by the Raman gain strength. Consequently, the vibrational lines (purple lines) appear first, then the ortho-ro-vibrational (blue lines) [denoted by Stage i) in Fig. 2], and finally the para-ro-vibrational Raman lines ν (green) [Stage ii) in Fig. 2]. The latter are sidebands of the black line ν (pump), and blue lines ν (ortho-ro-vibrational Raman). This sequence was corroborated theoretically (see below). In order to determine the dominating mechanism between cascaded SRS and Kerr effect in the red lines generation, we qualitatively assessed the Raman gain in the transient regime [24-25] and Kerr gain [26]; or alternatively their associated characteristic lengths (see Appendix A: Competition between SRS and FWM). The results show that the Kerr characteristic length is 110 shorter than the Raman characteristic length, thus indicating that the generation of lines at ν FWM1 , ν FWM2 , ν FWM3 , ν FWM4 is dominated by Kerr nonlinearity. Consequently, the red lines result from a cascaded parametric FWM process, which is further facilitated by the combination of three favorable circumstances: (i) The near-commensurability of ν R2 and ν R3 , (ii) The small frequency spacing between two adjacent green lines, and (iii) The PMC weak dispersion (see below and figure 3). Below, we discuss the involved comb lines behind this wave-mixing, and the origin of such a Kerr-type nonlinearity. These red-colored lines can result from wave-mixing involving several combinations of ν . For example, looking at the ν FWM2 spectral line 269.8 THz in Fig. 1(f), we identify several possible FWM processes to generate this line, like ν FWM2 = 2 ν - ν or ν FWM2 = ν FWM1 + ν - ν or ν FWM2 = ν + ν - ν FWM1 , and so on. Our aim now is to determine which of the possible FWM processes and spectral combinations lead to the generation of the red lines. Due to complicated dynamics involving high-order nonlinearities, the cascaded nature of the considered process, and the influence of phase-matching, we found it useful to start by recalling the expression for the intensity of the new frequency ν by FWM process from frequencies ν , ν and ν [27]: 𝐼 = 9 𝜋 𝜒 (3) 2 (𝜈 , 𝜈 , 𝜈 , 𝜈 )𝐼 𝐼 𝐼 𝜈 𝜀 𝑠𝑖𝑛 (𝛥𝑘𝑧) (1) Where χ (3) ( ν , ν , ν , ν ) is the frequency-dependent nonlinear susceptibility, I to I are the intensities at the participating frequencies, z is the propagation length, and Δk is the associated wavenumber mismatch. Expression (1) shows three possible dominating effects in the generation of red spectral lines: i) the strength of frequency-dependent nonlinear susceptibility (see Appendix B: Susceptibility in Raman-active medium); ii) the powers of the participating lines; and iii) the mismatch Δk between the participating frequencies. Before trying to identify the leading effects in the observed FWM, we can easily disregard the nonresonant Kerr susceptibility, as the estimated characteristic length L NL = c ⁄ (2 π . ν . I . n ) necessary for the generation of FWM sidebands by nonresonant Kerr process is ~ 150 m for n = 0.3x10 -8 cm /TW, which is much larger than the fiber length. On the other hand, we can calculate a rough estimate of Raman-induced effective value of n by comparing the characteristic length observed experimentally for the red lines with that of the Kerr effect. We obtain an effective value of n = 3x10 - 8 cm /TW, which corresponds to the value of 5.5x10 -8 measured for D [28]. In the Appendix B, we have derived a rigorous expression for the nonlinear polarization in a Raman medium using the standard adiabatic elimination of the upper Raman states. The resultant expression is cumbersome to investigate, however, for the aim of current analysis it suffices to extract the terms that are not in resonance with the Raman transition, and keep the dominant from the remaining terms. The expression of a representative term is: 𝑃 𝑅−𝐹𝑊𝑀 (𝑡) = 𝑁𝜇 𝑗𝑎2 (8𝜋ℏ) |𝜇 𝑗𝑎 | ∑ ∑ 𝐸 𝑞∗ 𝐸 𝑞′ 𝐸 𝑞" 𝑒 𝑖2𝜋(𝜈𝑞−𝜈𝑞′−𝜈𝑞′′)𝑡 (𝜈 𝑞 −𝜈 𝑞′ )(𝜈 𝑗 −𝜈 𝑞′ )(𝜈 𝑗 −𝜈 𝑞" )𝑞"𝑞,𝑞′≠𝑞 (2) This Raman induced non-resonant polarization clearly shows that the polarization responsible for the FWM process increases for fields with close frequencies due to the ν q - ν q' term appearing in the denominator (see Appendix B for notations). Based on the above, we conclude that the FWM process contributing to the red line generation stems from Raman induced enhanced nonlinearity. Second, by observing that the ratio of the values of the susceptibility χ (3) ( ν , ν , ν , ν ) for two different possible FWM processes is determined by the ratio of the smallest difference between the participating frequencies, we can deduce that among all the possible spectral lines fulfilling the required nergy conservation, only those with small frequency difference dominate. For example, we found this ratio to have a maximum value of 2 for the relevant FWM processes from neighboring lines, while mechanisms involving non-neighboring lines will be suppressed by the factor in the denominator. Further shortlisting of the dominant FWM frequency combinations can be done using the phase-matching condition. Figure. 3 shows the coherence length, L coh ≡ π . ∆k -1 , for three representative FWM frequency combinations. The first two correspond to ν FWM2 = 2 ν - ν and ν FWM1 = ν + ν - ν , respectively, and include frequencies from only neighboring lines. The third is of type ν FWM1 = ν + ν - ν , ν FWM1 = ν + ν - ν containing larger spaced frequencies. Fig. 2 . Schematic sequence of the processes during the comb generation. The pump laser line is identified by black color, the vibrational Stokes and anti-Stokes by purple (not visible), the Stokes and anti-Stokes lines due to ortho-ro-vibrational Raman by blue, the Stokes and anti-Stokes lines due to para-ro-vibrational Raman by green, and four-wave-mixing generated sidebands by red. The results show that for the majority of relevant processes from neighboring lines L coh is longer than the fiber length, so that phase mismatch is fulfilled, and does not play a significant role in selecting the FWM mechanism except for suppressing contributions from non-neighboring lines. On the other hand, the difference between the intensities of the different line series is typically one order of magnitude. Therefore, we can conclude that the power of the participating frequencies will play the dominant role in selecting the FWM process which provides the maximum contribution to the generation of new lines. Based on that, we expect that the strongest process which contributes to the generation of the ν FWM1 line is ν FWM1 = ν + ν - ν process, while the subsequent ν FWM2 = ν + ν - ν FWM1 and ν FWM2 = - ν processes are mainly responsible for the generation of ν FWM2 lines. These processes are indicated as Stage iiia) and Stage iiib) in Fig. 2. Other processes certainly provide additional and coherent contribution to the abovementioned ones, since the offsets of all frequencies are roughly multiples of Δ . For deuterium gas we get a clear comb-like spectrum with a step of approximately Δ = 1.75 THz. The generation lines at ν FWM3 and ν FWM4 follow similar process as for ν FWM1 and ν FWM2 . Finally, the spectrum shows very weak orange lines in Fig. 1(e-f), which are equally spaced from their adjacent lines by ~ Δ /2. Because of the same arguments above, these lines result from further mixing of the lines. Due to their weakness, we do not find it worthwhile to identify their originating lines. Fig. 3 . The coherence length for the ν FWM2 = 2 ν - ν type of FWM process (red curve), for the ν FWM1 = ν + ν - ν type of FWM process (green curve), and for the ν FWM1 = ν + ν - ν type of FWM process (blue curve). The dashed black line indicates the 3 m fiber length. III. N UMERICAL AND THEORETICAL APPROACH
To corroborate the above conclusions and to achieve a deeper understanding of the processes which govern the generation of FWM sidebands in Raman-gas medium, we have developed a numerical model. The simulation of the Raman comb generation was performed by numerically solving the generalized nonlinear Schrödinger equation in the slowly varying envelope approximation which is written as [23, 29]: 𝜕𝐴(𝜔)𝜕𝑧 = 𝑖 [𝑘(𝜔) − 𝑘(𝜔 ) − (𝜔 − 𝜔 ) 𝜕𝑘(𝜔)𝜕𝜔 | 𝜔 ] 𝐴(𝜔) + 𝑖𝜔2𝑐𝜀 ∑ 𝑃 𝑅,𝑗 (𝜔) 𝑗=2,3 (3) where
A(ω) is the field envelope in the spectral presentation, k(ω) is the wavenumber which includes both gas dispersion (calculated using the Sellmeyer-type expressions), the waveguide dispersion (determined from [23]), the phase mismatch, and the fiber loss spectrum. The quantity ω is the pump central angular frequency, and P R,j is the Raman polarization induced by the Raman transition Rj . In our case, we consider only the two ro-vibrational resonances of ortho- and para-deuterium, and ignore the vibrational transition. This choice is motivated by the large value of the vibrational Raman frequency, the higher fiber loss at ν ±1,0,0 , and correspondingly the lower number of sidebands generated by the vibrational Stokes and anti-Stokes. In this simulation, we have deliberately neglected the nonresonant Kerr response since it is weak in comparison to nonresonant Raman response, and to determine the role of Raman in the non-resonant polarization. The Raman olarization is given by: 𝑃 𝑅,𝑗 (𝑡) = 𝐴(𝑡)𝜀 𝑝̃ 𝑗 (𝛼 𝜌,𝑗 𝜌 + 𝛼 𝑤,𝑗 𝑤 𝑗 ) (4) where p ̃ j are pressures, ρ and w j are the off-diagonal element of the density matrix and the inversion, correspondingly, and α p,j and α w,j are coefficients which are related to the dipole moment of the Deuterium molecules. The temporal evolution of the quantities ρ and w j is described by standard two-photon two-level Bloch equations: 𝜕𝜌 𝜕𝑡 = 𝜌 [− 𝑗 + 𝑖ω 𝑅𝑗 + 𝛼̃ |𝐴(𝑡)| ] +𝛼 𝑤 𝑗 |𝐴(𝑡)| (5) 𝜕𝑤 𝑗 𝜕𝑡 = − 𝑤 𝑗 +1𝑇 𝑗 + 𝛼 (𝜌 − 𝜌 )|𝐴(𝑡)| (6) Here t j and T j are the polarization and inversion decay times, α , α , and α are coefficients which describe the coupling between electric field and the state of a molecule. All of the above coefficients are determined phenomenologically, as described in e.g. [29]. Fig. 4 . (a-c) Evolution of the numerical output spectra for the case of a coupled pump power of 12.8 W (a), 14.4 W (b) and 16 W (c). The fiber length was 3 m long, the D gas pressure was 40 bar. The pump laser line is shown with black color, the ortho-ro-vibrational shifts with blue color, the para- ro-vibrational with green color, and all the FWM-generated sidebands with red color. A single carrier frequency along with envelope varying on a few-fs time scale were used to treat all the comb components simultaneously and self-consistently. Figure 4 shows the calculated evolution of the output spectra for the case of coupled pump power of 12.8 W [Fig. 4(a)], 14.4 W [Fig. 4(b)] and 16 W [Fig. 4(c)] respectively. These numerical output spectra present the same pattern as the experimental one and in a good agreement with the generation sequence of Fig. 2. The color scheme is the same as in Fig. 1, with the pump laser line ν shown by black color, the ortho-D ν lines by blue, the para-D ν lines by green, and all the FWM sidebands ν FWM shown by red color. Besides the very good agreement between the numerical and experimental data, the main conclusion of the numerical approach is the mechanism and sequence of processes, which yield the comb in a molecular Raman gas medium. In the simulation, as the propagation distance increases, ortho-D ro-vibrational comb (blue) appears first during Stage i). The lines of this comb are sufficiently strong, that during the Stage ii) para-D ro-vibrational comb (green) appears on both sides of each line in ortho-D ro-vibrational comb. Figure. 4(a) shows perfectly the generated and classic Raman comb only based on the pure Raman transitions. After that, as a third step, FWM processes set in. They include the ν FWM1 = ν + ν - ν type of FWM [Stage iiia) and Fig. 4(b)] as explained above, followed by ν FWM2 = ν + ν - ν FWM1 and ν FWM2 = ν - ν [Stage iiib) and Fig. 4(c)]. This sequence is a cascaded one, since the ν FWM1 is required for the ν FWM2 to be generated. Figure 4 also confirms such a sequence, where one can see that the ν FWM1 lines are stronger than ν FWM2 . Finally, very weak orange lines are not visible in the numerical simulations. This sequence of the elementary processes, in particular the mechanism of the generation of the red FWM lines, agrees with the analysis presented in the previous section, which shows that the small and single frequency spacing is related to the two-step Raman process, followed by cascaded parametric FWM, which is based on Kerr-type nonlinearity provided by the medium Raman response.
IV. C ONCLUSION
In conclusion, we reported on the generation of a high power, broadband and single-frequency spaced comb based on a sequence of SRS and parametric four-wave in a multiple Raman resonance molecular gas medium. The comb frequency spacing Δ results from the fact that two Raman frequencies of the gas are multiple of Δ . This new optical comb generation process unlocks several limitations of current comb generation techniques. It can combine the multi-octave bandwidth of HC-PCF based Raman combs with single and smaller frequency spacing than in current spectra used for waveform synthesis [1] and attosecond pulse generation [9]. Also, unlike with micro-resonator or SRS comb, the frequency spacing can be actively controlled. This can be achieved for example by exciting a Raman medium, having at least two nearly commensurable Raman resonance frequencies, with two pump lasers whose frequency difference can be tuned to a value equal to a fraction of one of the two commensurable resonance frequencies. Finally, the results represent an important milestone towards exploring this Raman-Kerr comb into a frequency ruler [8]. However, unlike the presented pumping scheme, this requires specific pumping schemes such as the use of pump laser ith stabilized carrier-envelop phase, and with sufficiently high repetition rate so the required small number of Raman excited molecules survive the duration between two adjacent pump pulses [30]. A CKNOWLEDGMENTS
Funding: Agence Nationale de la Recherche (ANR) (PhotoSyn th, Labex SigmaLim, UV factor); Région Limousin; Air Force Office of Scientific Research (AFOSR) (FA9+550-14-1-0024); National Science Foundation (NSF) (PHY-1068865). The authors thank the PLATINOM platform in their help in the fiber fabrication. A PPENDIX A: C OMPETITION BETWEEN
SRS
AND
FWM
In order to have a qualitative assessment on the dominant process between SRS and FWM in the generation of the red line (Fig. 1), we proceed by estimating their effective gain or alternatively their generation characteristic lengths. For both SRS and Kerr effect, we can write the following for Raman net gain: 𝐺 𝑅 = 𝑔 𝑅 (𝑃 𝑅 /𝐴)𝐿 = 𝐿𝐿 𝑅−1 (A1) with 𝑔 𝑅 (𝑃 𝑅 𝐴⁄ ) = 𝐿
𝑅−1 (A2) and for the Kerr process, the characteristic length required to reach the power P si can be derived as [27]: 𝐿 𝐾 = 𝜆2𝜋𝑛 √ 𝐴 𝑃 𝑠𝑖 𝑃 𝐾3 (A3) Although the generation of a new frequency from three existing components is not a gain-type process, the following gain coefficient can be formally associated with the relevant four-wave mechanism: 𝑔 𝐾 = 𝐿 𝐾−1 (A4)
Here, 𝑃 𝑅 and 𝑃 𝐾 are the optical powers of the exciting fields involved in SRS and FWM respectively. In the expression of the Kerr gain, we ignore the phase matching condition between the fields in the FWM. 𝐿 𝑅 and 𝐿 𝐾 are the generation characteristic lengths for the process of SRS and Kerr respectively. The stronger the generation process the shorter the characteristic length. Consequently, finding out whether FWM is dominant or not is determined by the ratio ( L K / L R ). The quantity λ is the wavelength of the pumping field. n is the nonlinear index of the Raman driven D , and g R is the Raman gain coefficient. Recalling that in the transient regime, the Raman gain coefficient g R is equal to g R,t which is related to the steady-state Raman gain coefficient g R,ss 𝑔 𝑅,𝑠𝑠 by [24]: 𝑔 𝑅,𝑡 = 𝑅 × 𝑔
𝑅,𝑠𝑠 (A5) with
𝑅 = 2√ Γ 𝜏 𝑔 𝑅,𝑠𝑠 (𝑃 𝑅 /𝐴)𝐿 (A6) And, the pulse duration of the Stokes, τ S is given by [25]: 𝜏 𝑆 = √ 𝜏16𝑔 𝑅,𝑠𝑠
Γ(𝑃 𝑅 /A)L (A7) Here, Γ is the Raman dephasing rate, τ is the pulse duration of the Raman pumping field. Let us take the case of cascade SRS generation of one of the spectral lines of interest, the line at frequency ν (labelled FWM1). The latter is the second order Stokes of the blue line at ν or the first order Stokes of the green line at ν . Using the above expressions of the gain coefficient and the pulse duration (equations (A6) and (A7)), we readily find the Raman gain coefficient for the generation of the red line ν . Setting the average power of the blue line ν to 1.6 W, and that of the green line ν to 0.16 W, which are consistent with the measured output spectrum when 16 W of the laser power is sent to the fiber, we find g R,t ≈ 10 -6 cm.GW -1 and L R ≈ 495 cm for the generation of the red line ν . Similarly, we estimate the Kerr gain and characteristic length by setting the pumping average power to 1.6 W, and taking n = 3x10 -8 cm /TW, the value estimated as explained in the text of the paper [28]. We found g K ≈ 3.8x10 -6 cm.GW -1 and L K ≈ 4.5 cm. These approximate calculations show that L R is 110 times larger than L K , thus indicating that Kerr dominates the red line ν . Similar approach shows that this is valid for all red lines. A PPENDIX B: S USCEPTIBILITY IN R AMAN - ACTIVE MEDIUM
The full expression for the susceptibility in the Raman-active medium is discussed. Our derivation of the expression for the nonlinear susceptibility of the Raman medium follows the general line and denotations of Ref. [29]. We begin the derivation of this expression by considering a system which can be described by a three-level system with two lower levels |𝑎⟩ and |𝑏⟩ with energies ħω a and ħω b , ω a < ω b , and higher level |𝑗⟩ with energy ħω j . Without the loss of generality we ω a = 0 and ω a = ω R , where ω R is the Raman shift. The higher level j can be understood as a combination of multiple higher levels. The external field has three frequency components: 𝐸(𝑡) = ∑ 𝐸 𝑗 𝑒 −𝑖𝜔 𝑗 𝑡𝑗=1,3 + 𝑐. 𝑐. (B1) which couple the energy levels of the system through the Hamiltonian 𝐻 = ( 0 0 𝑉 𝑎𝑗 𝜔 𝑅 𝑉 𝑏𝑗 𝑉 𝑗𝑎 𝑉 𝑗𝑏 ħ 𝜔 𝑗 ) (B2) Here 𝑉 𝑎𝑗 = 𝑉 𝑗𝑎∗ = −𝜇 𝑎𝑗 𝐸 (B3) 𝑉 𝑏𝑗 = 𝑉 𝑗𝑏∗ = −𝜇 𝑏𝑗 𝐸 (B4) As one can see from these expressions, the levels |𝑎⟩ and |𝑗⟩ , as well as |𝑎⟩ and |𝑗⟩ are coupled by electric dipole transition, while the transition between |𝑎⟩ and |𝑏⟩ is forbidden. n this case, the amplitudes C n of the wavefunction of the system |𝜓⟩ = ∑ 𝐶 𝑛 𝑒 −𝑖𝜔 𝑛 𝑡𝑛=𝑎,𝑏,𝑗 |𝑛⟩ (B5) satisfy the following equations which follow from the Schrödinger equation: 𝜕𝐶 𝑛 𝜕𝑡 = − 𝑖ħ ∑ 𝑒 −(𝜔 𝑛 −𝜔 𝑘 )𝑡 𝑉 𝑛𝑘 𝐶 𝑘𝑘 (B6) This set of linear equations is solved by usual methods, utilizing the adiabatic elimination of C j [29]. The details of the derivation are standard and too cumbersome to be presented here. After the expression for the coefficients C n are obtained, the polarization of the medium can be expressed through 𝑃 = 𝑁 ∑ 𝜇 𝑗𝑎 𝐶 𝑎 𝐶 𝑗∗ 𝑒 −𝑖(𝜔 𝑎 −𝜔 𝑗 ) + 𝑐. 𝑐. 𝑛=𝑎,𝑏 (B7) N being the particle density. The final expression for the polarization reads 𝑃(𝑡) = 𝑁𝜇 𝑗𝑎2 (|𝜇 𝑗𝑎 | 𝐴 + |𝜇 𝑗𝑏 | 𝐵) (B8) Where
𝐴 = {∑ [− 𝐸 𝑞 𝐸 𝑞′ 𝑒 −𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑞 +𝜔 𝑞′ )(𝜔 𝑗 −𝜔 𝑞′ ) + 𝐸 𝑞∗ 𝐸 𝑞′∗ 𝑒 −𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑞 −𝜔 𝑞′ )(𝜔 𝑗 −𝜔 𝑞′ ) ] 𝑞,𝑞 ′ +∑ [− 𝐸 𝑞 𝐸 𝑞′∗ 𝑒 −𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑞 −𝜔 𝑞′ )(𝜔 𝑗 +𝜔 𝑞′ ) + 𝐸 𝑞∗ 𝐸 𝑞′ 𝑒 −𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑞 −𝜔 𝑞′ )(𝜔 𝑗 −𝜔 𝑞′ ) ] 𝑞,𝑞 ′ ≠𝑞 +𝑐. 𝑐. } × {∑ [ 𝐸 𝑞′′ 𝑒 −𝑖𝜔𝑞′′𝑡 𝜔 𝑗 −𝜔 𝑞′′ + 𝐸 𝑞′′∗ 𝑒 𝑖𝜔𝑞′′𝑡 𝜔 𝑗 +𝜔 𝑞′′ ] 𝑞 ′′ + 𝑐. 𝑐. } (B9) And
𝐵 = {∑ [− 𝐸 𝑞 𝐸 𝑞′ 𝑒 −𝑖(𝜔𝑞+𝜔𝑞′)𝑡 (𝜔 𝑅 −𝜔 𝑞 −𝜔 𝑞′ )(𝜔 𝑗 −𝜔 𝑞′ ) + 𝑞,𝑞 ′ 𝐸 𝑞∗ 𝐸 𝑞′∗ 𝑒 𝑖(𝜔𝑞+𝜔𝑞′)𝑡 (𝜔 𝑅 +𝜔 𝑞 +𝜔 𝑞′ )(𝜔 𝑗 +𝜔 𝑞′ ) ] +∑ [− 𝐸 𝑞 𝐸 𝑞′∗ 𝑒 −𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑅 −𝜔 𝑞 +𝜔 𝑞′ )(𝜔 𝑗 +𝜔 𝑞′ ) + 𝑞,𝑞 ′ ≠𝑞𝐸 𝑞∗ 𝐸 𝑞′ 𝑒 𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑅 +𝜔 𝑞 −𝜔 𝑞′ )(𝜔 𝑗 −𝜔 𝑞′ ) ] + 𝑐. 𝑐. } × {∑ [ 𝐸 𝑞′′ 𝑒 −𝑖𝜔𝑞′′𝑡 𝜔 𝑗 −𝜔 𝑅 −𝜔 𝑞′′ + 𝑞 ′′ 𝐸 𝑞′′∗ 𝑒 𝑖𝜔𝑞′′𝑡 𝜔 𝑗 −𝜔 𝑅 +𝜔 𝑞′′ + 𝐸 𝑞′′∗ 𝑒 𝑖𝜔𝑞′′𝑡 𝜔 𝑗 −𝜔 𝑞′′ + 𝐸 𝑞′′ 𝑒 −𝑖𝜔𝑞′′𝑡 𝜔 𝑗 +𝜔 𝑞′′ ] + 𝑐. 𝑐. } (B10) In this expression, indices q , q' , q'' denote pump fields and can take values from 1 to 3. In deriving this expression, we have neglected the decay rates which are for our case much lower than the few-THz difference between frequencies. This expression includes all contributions provided by the Raman nonlinearity, including Raman contributions to self- and cross-phase modulation and in particular four-wave mixing. Equation (B10) contains in total 1440 terms, therefore in the main text we present a susceptibility corresponding to only one of these contributions. Since the values of dipole momenta µ ja and µ ja are similar, the dominant terms are determined by the smallest value of the denominator. For typical conditions, ω j = ω R and ω j = ω q , therefore the factors in the denominator which contain ω j will have similar value for all terms, with terms containing ω j – ω q in the denominator being larger. On the other hand, the terms in A which contain in the denominator the combination ω q – ω q' , and the terms in B which contain in the denominator combinations ω R ± (ω q – ω q' ) , are larger than other terms and will provide dominant contribution to the respective nonlinear processes. Considering the above arguments, we identify the dominant contributions A dom and B dom to A and B as follows: 𝐴 𝑑𝑜𝑚 = {∑ 𝐸 𝑞∗ 𝐸 𝑞′ 𝑒 𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑞 −𝜔 𝑞′ )(𝜔 𝑗 −𝜔 𝑞′ )𝑞,𝑞 ′ ≠𝑞 +𝑐. 𝑐. } {∑ 𝐸 𝑞′′ 𝑒 −𝑖𝜔𝑞′′𝑡 𝜔 𝑗 −𝜔 𝑞′′ 𝑞 ′′ + 𝑐. 𝑐. } (B11) And 𝐵 𝑑𝑜𝑚 = {∑ 𝐸 𝑞 𝐸 𝑞′ 𝑒 𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑅 −𝜔 𝑞 −𝜔 𝑞′ )(𝜔 𝑗 −𝜔 𝑞′ )𝑞,𝑞 ′ +∑ 𝐸 𝑞∗ 𝐸 𝑞′ 𝑒 𝑖(𝜔𝑞−𝜔𝑞′)𝑡 (𝜔 𝑅 +𝜔 𝑞 −𝜔 𝑞′ )(𝜔 𝑗 −𝜔 𝑞′ )𝑞,𝑞 ′ + 𝑐. 𝑐. } × {∑ [ 𝐸 𝑞′′ 𝑒 −𝑖𝜔𝑞′′𝑡 𝜔 𝑗 −𝜔 𝑅 −𝜔 𝑞′′ + 𝑞 ′′ 𝐸 𝑞′′∗ 𝑒 𝑖𝜔𝑞′′𝑡 𝜔 𝑗 −𝜔 𝑞′′ ] + 𝑐. 𝑐. } (B12) REFERENCES
1. H.-S. Chan, Z.-M. Hsieh, W.-H. Liang, A. H. Kung, C.-K. Lee, C.-J. Lai, R.-P. Pan, and L.-H. Peng, “
Synthesis and measurement of ultrafast waveforms from five discrete optical harmonics ” Sci. , 1165 (2011)
2 S. B. Papp, K. Beha, P. Del’Haye, F. Quinlan, H. Lee, K. J. Vahala, and S. A. Diddams, “
Microresonator frequency comb optical clock,”
Optica (1) 10-14 (2014) 3. J. Pfeifle, V. Brasch, M. Lauermann, Y. Yu, D. Wegner, T. Herr, K. Hartinger, P. Schindler, J. Li, D. Hillerkuss, R. Schmogrow, C. Weimann, R. Holzwarth, W. Freude, J. Leuthold, T. J. Kippenberg, and C. Koosl., “ Coherent terabit communications with microresonator Kerr frequency combs, ” Nat. Phot. , 375–380 (2014) 4. T. Ideguchi, S. Holzner, B. Bernhardt, G. Guelachvili, N. Picque, and T. W. Hansch, “ Coherent Raman spectro-imaging with laser frequency combs, ” Nature , 355-358 (2013) 5. H. W. Mocker, and R. J. Collins, “
Mode competition and self-locking effects in a Q-switched ruby laser,”
Appl. Phys. Lett. (270) (1965) 6. S. E. Harris and A. V. Sokolov, “ Broadband spectral generation with refractive index control,”
Phys. Rev. A (6) 4019-4022 (1997) 7. H-S. Chan, Z.-M. Hsieh, W.-H Liang, A. H. Kung, C.-K. Lee, C.-J. Lai, R.-P. Pan, and L.-H. Peng, “ Synthesis and measurement of ultrafast waveforms from five discrete optical harmonics, ” Science (6021), 1165–1168 (2011) T. Udem, R. Holzwarth, and T. W. Hänsch, “
Optical frequency metrology ”, Nature , 233 (2002) 9 F. Krausz and M. Ivanov, “
Attosecond physics ”, Rev. Mod. Phys. , 163 (2009) 10 P. Del’Haye, A. Schliesser, O. Arcizet, T. Wilken, R. Holzwarth, and T. J. Kippenberg, “ Optical frequency comb generation from a monolithic microresonator, ” Nature , 1214-1217 (2007) 11. F. Couny, F. Benabid, P. J. Roberts, P. S. Light and M. G. Raymer, “
Generation and photonic guidance of multi-octave optical-frequency combs, ” Science , 1118-1121 (2007) 12. A. Benoît, B. Beaudou, M. Alharbi, B. Debord, F. Gerome, F. Salin, and F. Benabid, “
Over-five octaves wide Raman combs in high-power picosecond-laser pumped H2-filled inhibited coupling Kagome fiber,”
Opt. Exp. (11), 14002-14009 (2015) 13. Y. Y. Wang, C. Wu, F. Couny, M. G. Raymer, and F. Benabid, “ Quantum-Fluctuation-Initiated Coherence in Multioctave Raman Optical Frequency Combs,”
Phys. Rev. Lett. , (2010) 14. P. Hosseini, A. Abdolvand, and P. St.J. Russell, “
Generation of spectral clusters in a mixture of noble and Raman-active gases,”
Opt. Lett. (23), 5543-5546 (2016) 15. F. Tani, F. Belli, A. Abdolvand, J. C. Travers and P. S. Russell, “ Generation of three-octave-spanning transient Raman comb in hydrogen-filled hollow-core PCF ,” Opt. Lett. (6) (2015) 16. Y. R. Shen and N. Bloembergen, “ Theory of Stimulated Brillouin and Raman scattering, ” Phys. Rev. , A1787 (1965) 17. R. W. Boyd, M. G. Raymer, P. Narum and D. J. Harter, “
Four-wave parametric interactions in a strongly driven two-level systems ,” Phys. Rev. A. (1), 411 (1980) 18. F. Benabid, J. C. Knight, G. Antonopoulos, and P. S. J. Russell, “ Stimulated Raman scattering in Hydrogen-filled Hollow-core photonic crystal fiber, ” Science , 399 (2002) 19 A. Benoît, E. Ilinova, B. Beaudou, B. Debord, F. Gérôme, and F. Benabid, “
Spectral-temporal dynamics of high power Raman picosecond pulse using H -filled Kagome HC-PCF ”, Opt. Lett. (19), 3896-3899 (2017) 20. R. Z. Martínez, D. Bermejo, P. Wcisło, and F. Thibault, “ Accurate wavenumber measurements for the S (0), S (1), and S (2) pure rotational Raman lines of D , ”. Journal of Raman Spectroscopy, (1), 127-129 (2019) 21. D. K. Veirs and G. M. Rosenblatt, “ Raman line positions in molecular hydrogen: H , HD, HT, D , DT, and T , ” J. Mol. Spectrosc. , 401 (1987) 22. A. Maczek, and A. Meijer, “ Statistical thermodynamics, ” Oxford Chemistry Primers Series, Oxford University Press, pp.63-64 (2017) 23. B. Debord, A. Amsanpally, M. Chafer, A. Baz, M. Maurel, J.-M. Blondy, E. Hugonnot, F. Scol, L. Vincetti, F. Gerome, and F. Benabid, “
Ultralow transmission loss in inhibited-coupling guiding hollow fibers ,” Optica (2), 290-217 (2017) 24 R. J. Heeman and H. P. Godfried, “ Gain reduction measurements in transient Stimulated Raman Scattering ” IEEE J. Quantum Electron. , 358 (1995) 25 M. G. Raymer, I. A. Walmsley, J. Mostowski, and B. Sobolewska, “ Quantum theory of spatial and temporal coherence properties of stimulated Raman scattering ” Phys. Rev. A , 332 (1985) 26 Y. Okawachi, M. Yu, V. Venkataraman, P. M. Latawiec, A. G. Griffith, M. Lipson, M. Lončar, and A. L. Gaeta, “ Competition between Raman and Kerr effects in microresonator comb generation ” Opt. Lett. , 2786 (2017) 27. R. W. Boyd, “
Nonlinear Optics. 3 rd edition”, Academic press, (2003) 28 S. Xu, B. Lavorel, O. Faucher, and R. Chaux, “
Characterization of self-phase modulated ultrashort optical pulses by spectral phase interferometry ”, Jour. Opt. Soc. Of A. B (1), 165-168 (2002) 29 V. P. Kalosha and J. Herrmann, “ Phase Relations, Quasicontinuous Spectra and Subfemtosecond Pulses in High-Order Stimulated Raman Scattering with Short-Pulse Excitation ,” Phys. Rev. Lett., (6), 1226 (2000) 30 D. T. Smithey, M. Belsley, K. Wedding, and M. G. Raymer, “ Near quantum-limited phase memory in a Raman amplifier ”. Phys. Rev. Lett.67