Arxiv:2011.01656V1 [Physics.Optics] 3 Nov 2020 K (3) 3 (3) of Them Combined

From Raman frequency combs to supercontinuum generation in nitrogen-ﬁlled hollow-core anti-resonant ﬁber

Shou-Fei Gao,1, 2 Ying-Ying Wang,1, 2, ∗ Federico Belli,3 Christian Brahms,3 Pu Wang,2 and John C. Travers3, † 1Guangdong Provincial Key Laboratory of Optical Fiber Sensing and Communication, Institute of Photonics Technology, Jinan University, Guangzhou, 511443, China 2National Center of Laser Technology, Institute of Laser Engineering, Beijing University of Technology, Beijing 100124, China 3School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, United Kingdom (Dated: November 4, 2020) We demonstrate a route to supercontinuum generation in gas-filled hollow-core anti-resonant fibers through the creation of a broad vibrational Raman frequency comb followed by continuous broadening and merging of the comb lines through either rotational Raman scattering or the optical Kerr effect. Our demonstration experiments, utilizing a single pump pulse with 20 ps duration at 532 nm in a nitrogen-filled fiber, produce a supercontinuum spanning from 440 nm to 1200 nm, with an additional deep ultraviolet continuum from 250 nm to 360 nm. Numerical results suggest that this approach can produce even broader supercontinuum spectra extending from the ultraviolet to mid-infrared.

Cascaded inelastic scattering of laser light from beyond 2000 nm). Experimentally we generate a flat su- molecules, known as stimulated Raman scattering (SRS), percontinuum spanning 440 nm to 1200 nm and an addi- can create broad sets of discrete spectral lines from nar- tional deep ultraviolet continuum from 250 nm to 360 nm. rowband driving lasers [1–3]. Such Raman frequency These results are achieved by pumping a nitrogen-filled combs have been widely studied in both free-space [4, 5] anti-resonant fiber with 20 ps pump pulses at 532 nm. and waveguide geometries [6, 7] and found application Theory and numerical simulations—We model the in the synthesis of trains of light-field transients [8–10] propagation of the electric field through the HC-PCF us- and single ultrashort pulses [11, 12] as well as for effi- ing a unidirectional pulse propagation equation [22, 23], cient wavelength conversion [13, 14]. By tightly confining ω both the gas molecules and laser light inside hollow-core ∂ E˜(z, ω) = i β(z, ω) − E˜(z, ω) z v(z) photonic-crystal fibers (HC-PCF), such as anti-resonant (1) iω2 guiding fibers, Raman frequency combs can be efficiently + P˜(z, ω), 2 generated directly from noise and do not require multi- 2c 0β(z, ω) color pumping [6, 7, 15–18]. However, the Raman combs where E˜(z, ω) = F[E(z, t)] is the Fourier transform produced to date, while spanning multiple octaves, have of the linearly polarized, mode-averaged electric field always consisted of discrete spectral lines with at most E(z, t), z is axial propagation distance through the fiber, moderate spectral broadening or shifting [17, 19]. ω is angular frequency, β(z, ω) is the frequency- and In this Letter we demonstrate, both experimentally pressure-dependent propagation constant of the funda- and numerically, that SRS can create a broadband and mental fiber mode, which we model using the hollow- smooth supercontinuum from a single narrowband pump capillary model [23, 24] with the pressure-dependent dis- laser. The key mechanism is a two-step process in which persion of the gas [25], v(z) is the group velocity of the the large vibrational frequency shift leads to a widely reference frame, c is the vacuum light speed, 0 is the per- spaced but extremely broadband frequency comb, which mittivity of free space, and P˜(z, ω) is the Fourier trans- is then broadened into a smooth supercontinuum by the form of the nonlinear polarization influence of either rotational SRS, the optical Kerr ef- P (z, t) = P k(z, t) + P i(z, t) + P v(z, t) + P r(z, t). (2) fect (instantaneous nonlinear refractive index), or both arXiv:2011.01656v1 [physics.optics] 3 Nov 2020 k (3) 3 (3) of them combined. Interestingly, by combining the vi- The Kerr polarization is P (t) = 0Nγ E where γ brational and rotational response, it is possible for the is the third-order hyperpolarizability and N is the num- continuum to form by the Raman effect alone. This ber density of molecules. The polarization due to ioniza- is an alternative regime of supercontinuum formation in tion and the plasma response, P i(t), is calculated using optical fibers—conventional routes, such as soliton self- the method of Geissler [26], using the Peremolov-Popov- compression and fission dynamics [20], modulational in- Terent’ev ionization rate [27] with an ionization potential stability [21], or self-phase modulation (SPM) and self- of 15.58 eV. The vibrational polarization P v and the ro- steepening, are all dominated by the contribution of the tational response P r will be discussed below. The numer- optical Kerr effect. Our numerical results suggest that ical parameters we use are given in Table I. We neglect our approach can produce supercontinuum spectra span- the loss of the fiber and the dispersion spikes arising due ning from the ultraviolet to the mid-infrared (270 nm to to the fiber resonances [28]. 2

(a) P k (b) P v (c) P r 0 −40

SPD (dB) 10 8 6 SPD 4 (dB) 2

Position (m) 0 0 0.2 0.6 1.0 0.2 0.6 1.0 0.2 0.6 1.0 −30

(d) P k + P v (e) P v + P r (f) P k + P v + P r −60 0 −40

SPD (dB) 10 8 6 4 2 Position (m) 0 0.2 0.6 1.0 0.2 0.6 1.0 0.2 0.6 1.0 Frequency (PHz) Frequency (PHz) Frequency (PHz)

FIG. 1. Numerical simulations of the spectral evolution. We model 20 ps duration, 532 nm wavelength, 80 µJ pump pulses in a 10 m long, 26 µm diameter hollow fiber filled with nitrogen with a pressure gradient from 1 to 40 bar. Each subplot includes the indicated contributions to the nonlinear polarization. SPD = relative spectral power density. These data are an ensemble average of 100 simulations to account for shot-to-shot fluctuations.

Figure 1(a) shows numerical simulation results when length, Aeff is the effective mode area of the fiber core, including only the Kerr effect (only P k in the model) for N¯ is the average number density, and the factor 2/3 ac- a 26 µm diameter hollow fiber and pulses with 20 ps counts for the pressure gradient. For our parameters this duration, 80 µJ energy and a wavelength of 532 nm. gives ∆ν = ∆ω/2π ≈ 2.6 THz, and as a consequence the We considered a nitrogen pressure gradient from 1 to spectrum barely changes. 40 bar (corresponding to our experimental parameters, The generation of a Raman frequency comb from a which are detailed later). Pure Kerr propagation with- single pump pulse can be initiated spontaneously from out the Raman contribution results in negligible broad- noise [35] and subsequently amplified through SRS [7]. ening, because the pump pulse and fiber parameters do The (nonlinear or inelastic) coupling of the light field not lead to either significant self-phase modulation or with the stretching and alignment of homonuclear di- modulational instability. The spectral width induced by atomic molecules is described through first-order pertur- SPM is approximately given by ∆ω ≈ φmax/τ0, with bations of a non-rigid rotor model [36–40]. As the vi- ¯ (3) φmax = (2/3)n2ωP0L/cAeff , where n2 = 3Nγ /40c is brational frequency νv is typically much larger than the J the nonlinear refractive index of the gas, P0 and τ0 are rotational frequencies νr , and large enough that only the the peak power and duration of the pulse, L is the fiber vibrational ground state is occupied at room temperature, the laser field interacts with a single vibrational transition. The time-dependent nonlinear polarization due to the vibrational coupling results as [36] TABLE I. Parameters for numerical simulations. Z ∞ v 2 v 2 Parameter Unit Value Ref. P = E(t)N(4π0) κv h(νv,T2 , τ)E(t − τ) dτ,(3) γ(3) 10−64 Cm4V−3 95.2 [29] 0 2 νv THz 69.85 [30] ∂α 1 κ = − , (4) ∆νv THz ∼0.04 [31] v ∂Q 8πµνv ∂α/∂Q 10−20 m2 1.75 [32] −26 µ 10 kg 1.16 where the Raman response function is h(ν, T2, t) = −1 B m 199 [30] sin(2πνt) exp(−t/T ), T = (π∆ν)−1 is the dephasing −4 −1 2 2 D 10 m 5.76 [30] time with ∆ν the full-width at half-maximum (FWHM) qJ odd 1 linewidth of the Raman transition, ∂α/∂Q is the isotropic qJ even 2 ∆α 10−31 m3 6.7 [33] averaged polarizability derivative with respect to the −1 ensemble-averaged molecular stretch Q, and µ is the re- br GHz amagat ∼3.3 [34] duced molecular mass. 3

(a) In Fig. 1(b), we include only the vibrational Raman 20 response P v. Our parameters correspond to pumping in the steady-state gain regime [35], and because the ra- (MW) 10 tio ∆νv/νv is relatively small (Table I) the Raman response gives rise to a vibrational frequency comb with Power 0 Power (b) well-separated peaks at both lower (Stokes) and higher 1.2 (MW) 20 (anti-Stokes) frequencies separated by νv. This is similar to previous results in light homonuclear gases such as hydrogen [6, 7, 15–18]. Note that each line results from scattering from the previous line using the same (funda- (m) 0.9 10 mental) vibrational transition, and hence the Stokes and anti-Stokes lines are equally spaced in frequency. Position In contrast to the vibrational levels, the spacing of the rotational levels is small enough that a manifold of levels 0 are populated at room temperature and therefore con- 0.6 −30 −20 −10 0 10 20 30 tribute to the nonlinear response. The energy eJ and Delay (fs) thermal population ρJ of each level are given by 2 2 eJ = hc[BJ(J + 1) − DJ (J + 1) ], (5) FIG. 2. Initial temporal dynamics. We focus only on the cen- qJ (2J + 1) exp(−eJ /kBT ) tral 60 fs of the 20 ps pulse, for the same parameters as Fig. 1 ρJ = P , (6) with all contributions included (nearly identical results are qK (2K + 1) exp(−eK /kBT ) K obtained for pure vibrational evolution for the propagation where J is the rotational quantum number, h is Planck’s length shown here). (a) Slices of the instantaneous power at constant, B and D are the rotational and centrifugal the input (blue) and after 0.91 m (orange), where the peak constants for the specific molecular species, kB is Boltz- power is highest (20 MW) and the pulse duration is 2.1 fs. (b) The instantaneous power evolution from 0.6 m to 1.2 m. mann’s constant, T is the temperature, and qJ is a sta- tistical weighting factor that varies with odd and even J to account for nuclear spin (see Table I). Selection rules limit the rotational transitions to pairs of levels with When combining the vibrational Raman and Kerr ef- ∆J = 2. In nitrogen we can neglect the centrifugal con- fects, the vibrational comb lines broaden into each other tribution (B/D > 105), and the rotational line spacing is (Fig. 1(d)), forming a broadband continuum. As noted J previously, the pump pulse does not drive significant νr = (eJ+2−eJ )/h = 4cB(J+3/2) with 4cB = 0.24 THz. Over 20 rotational levels play a role and the linewidth of SPM. However, in the time domain the vibrational fre- each level is pressure-broadened. The Lorentzian gain quency comb corresponds to an ultrashort pulse train. peaks associated with each rotational state have been Fig. 2 shows that after 0.91 m propagation, the vibra- −1 tional comb consists of a train of 2.1 fs pulses with a peak found to have a linewidth ∆νr = (πT2) = ρbr where ρ power of 20 MW. The combination of shorter duration is the gas density in Amagat and br is an experimentally determined constant [34]. For nitrogen at high pressure, and higher peak power enhances self-phase modulation the rotational Raman gain spectrum is continuous: at by a factor of ∼50000 compared to the 20 ps pump pulses with a peak power of 3.8 MW, leading to spectral broad- 40 bar, ∆νr = 0.15 THz, which is similar to the line J ening of the comb. Furthermore, the individual comb spacing νr [41, 42]. The rotational Raman polarization is approximated by lines broaden into each other to form a continuum. Pre- a sum over all allowed rotational transitions [36], vious Raman frequency comb experiments have either worked with much lower peak power and longer pump r 2 P = E(t)N(4π0) pulses [43], producing a pure frequency comb in which X Z ∞ SPM never occurs and no continuum is formed, or with × κJ h(νJ ,T r, τ)E(t − τ)2dτ, (7) r r 2 much shorter pulses, leading to immediate SPM broad- J 0 ening of the pump pulse itself [44]. With our parameters, 4π∆α2 (J + 1)(J + 2) ρ ρ J J+2 J in contrast, SPM occurs only after the formation of the κr = − (8) 15h 2J + 3 2J + 5 2J + 1 frequency comb. where ∆α is the molecular polarizability anisotropy. While broadening of the vibrational comb due to the The spectral evolution obtained for pure rotational Kerr effect leads to a continuum, it is not the only mecha- SRS is shown in Fig. 1(c). It shows different behaviour nism that can achieve this. Remarkably, combining both to the case of pure vibration: the spectrum appears vibration and rotation, but not the Kerr effect, as shown very smooth and continuous, and mostly down-shifted to in Fig. 1(e), leads to a pure Raman supercontinuum span- longer wavelengths; the overall spectral coverage is also ning multiple octaves achieved via spectral broadening much reduced. and merging of several vibrational lines. This works 4 well in nitrogen, because the spacing and linewidth of (a) 0 the rotational levels are similar, which creates a near- anti-Stokes 1 64.4 μJ 5 −20 6 4 continuous rotational Raman gain spectrum. However, 87 −40 we expect that this is a universal feature of molecular SPD (dB) −60 gases. Our numerical results indicate that in hydrogen 400 600 800 1000 1200 (not shown), where the rotational gain spectrum consists Wavelength (nm) of well-separated narrow lines, a smooth supercontinuum (b) 2.1 μJ 18.4 μJ can still be formed. a1 Figure 1(f) shows the full supercontinuum simulations, including all contributions to the nonlinear polarization. p s1 1.3 μJ 8.1 μJ s2 s3 Importantly it is the vibrational comb that establishes s4 the full continuum extent, covering 3 octaves, and then 1.1 μJ 5.6 μJ the Kerr and rotational Raman responses that broaden the vibrational lines to form a true continuum. Ionization SPD (20 dB / div) 0.5 μJ 3.8 μJ plays a negligible role in the dynamics—the low peak intensity (a maximum of 3 × 1013 W/cm2 is reached) leads to a negligible ionisation fraction of 0.002%. This is in 400 600 800 1000 1200 400 600 800 1000 1200 contrast to supercontinuum generation in gas-filled hol- Wavelength (nm) Wavelength (nm) low fibers pumped with shorter pulses based on modulational instability [45], where ionization contributes sig- FIG. 3. Experimental measurements. (a) The measured spec- nificantly. trum at the highest pulse energy of 64.4 µJ. The orange part of Experiment—We verified these dynamics experimen- the spectrum was measured with a dedicated DUV spectrom- eter. The inset shows a photograph of the dispersed spectrum. tally using a single-ring nodeless hollow-core anti- The vibrational anti-Stokes lines are labelled at 426 nm (1), resonant fiber consisting of 6 untouching thin tubes with 356 nm (4), 331 nm (5), 305 nm (6), 287 nm (7) and 269 nm diameter of 12 µm and membrane thickness of 210 nm (8) (b) The spectral evolution as the input pulse energy is forming a negative-curvature core shape with inscribed gradually increased from 0.5 µJ to 18.4 µJ as indicated. The diameter of 26 µm. The fiber exhibits a minimum at- pump line is labelled (p) as are the discrete vibrational Stokes tenuation of 80 dB/km at the pump laser wavelength lines at 607 nm (s1), 707 nm (s2), 847 nm (s3), and 1056 nm of 532 nm and its first-order transmission band spans (s4). from 440 nm to > 1200 nm with single-mode opera- tion [46]. The ultraviolet transmission band was not measured, but results from finite element modelling indicate input pulse energy is gradually increased from 0.5 µJ that it spans 250 nm to 360 nm. We used a 10 m length to 18.4 µJ. At low pulse energies, discrete vibrational of fiber with the input end open to atmosphere and the Stokes lines at 607 nm (s1), 707 nm (s2), 847 nm (s3), output end sealed inside a gas cell filled with 40 bar of and 1056 nm (s4) are observed. A further increase of nitrogen. Our 532 nm pump laser emitted 20 ps duration the input pulse energy shows spectral broadening of the pulses at 1 kHz, and was coupled to the fiber through a vibrational lines. Increasing the pulse energy also results telescope and plano-convex lens. The output spectrum in the generation of an anti-Stokes vibrational line at was collected by an integrating sphere and two optical 473 nm (a1). That these lines require higher energies can spectrum analyzers. be attributed to the fiber loss. All of the vibrational lines Figure 3(a) shows the measured spectrum at the high- generated experience spectral broadening and merge to est pulse energy of 65 µJ coupled into the fibre. A form the supercontinuum at around 8 µJ of pump energy. wide supercontinuum spectrum is generated, in agree- The DUV part of the continuum shown in Figure 3(a) ment with our numerical simulations, spanning more corresponds to vibrational anti-Stokes lines at 356 nm than 775 nm from 440 nm to 1200 nm (the long- (a4), 331 nm (a5), 305 nm (a6), 287 nm (a7), 269 nm wavelength limit of our measurements). Furthermore, we (a8) and their subsequent spectral broadening. These observe deep ultraviolet (DUV) generation from 260 nm dynamics are in excellent agreement with our numerical to 360 nm. The presence of the DUV component is fur- modelling. ther confirmed by the photograph of the dispersed spec- Discussion—Experiments in molecular gases in HC- trum shown as an inset. The spectral gap from 375 nm PCF to date have mostly concentrated on pumping to 440 nm is due to the high resonance loss in the hollow- with pulses of much longer duration than the Raman core fiber [46], resulting in the absence of the 2nd and 3rd oscillation period. Apart from comb generation, this vibrational anti-Stokes (AS) lines (426 nm and 387 nm). system has also been used for extremely efficient Ra- This could be minimized through further optimization of man conversion and generation [13, 47], continuous-wave the uniformity of the anti-resonant tubes of the fiber. pumped gas-Raman lasers [43], solitary pulse genera- Figure 3(b) shows the output spectral evolution as the tion [48], self-similar evolution [49] and supercontinuum 5 up-conversion [50]. With much shorter pump pulses, Ra- man effects combine with a strong Kerr contribution, allowing for enhanced supercontinuum generation [51]. ∗ [email protected] In the impulsive regime, with a pump pulse duration † [email protected] close to the Raman oscillation period, Raman dynam- [1] T. Imasaka, S. Kawasaki, and N. Ishibashi, Generation of ics combined with Kerr-driven soliton self-compression more than 40 laser emission lines from the ultraviolet to have enabled supercontinuum generation from the vac- the visible regions by two-color stimulated raman effect, uum ultraviolet to the near infrared [44], ultrafast coher- Applied Physics B 49, 389 (1989). ent pulse shaping of a dispersive wave in the ultravio- [2] S. E. Harris and A. V. Sokolov, Broadband spectral gen- let [52], and asymmetrical spectral broadening for pulse eration with refractive index control, Phys. Rev. A 55, R4019 (1997). post-compression [53]. [3] S. Baker, I. A. Walmsley, J. W. G. Tisch, and J. P. The use of intermediate pulse durations of a few pi- Marangos, Femtosecond to attosecond light pulses from coseconds leads to completely different nonlinear dynam- a molecular modulator, Nature Photonics 5, 664 (2011). ics characterized by efficient noise-seeded generation of [4] A. V. Sokolov, D. R. Walker, D. D. Yavuz, G. Y. Yin, vibrational lines and a frequency comb followed by spec- and S. E. Harris, Raman generation by phased and tral broadening and supercontinuum generation. Re- antiphased molecular states, Phys. Rev. Lett. 85, 562 (2000). cently a single frequency-shifted vibrational line was gen- [5] H. Kawano, Y. Hirakawa, and T. Imasaka, Generation erated by pumping an atmospheric air-filled HC-PCF of high-order rotational lines in hydrogen by four-wave with 6 ps pulses [54], and this line broadened, along with Raman mixing in the femtosecond regime, IEEE Journal the pump spectrum, into a high-power sub-octave contin- of Quantum Electronics 34, 260 (1998). uum. Here we have shown that even a broad frequency [6] F. Benabid, J. C. Knight, G. Antonopoulos, and P. S. J. comb, consisting of a multitude of vibrational lines (we Russell, Stimulated raman scattering in hydrogen-filled experimentally demonstrate 10 vibrational lines, numer- hollow-core photonic crystal fiber, Science 298, 399 (2002). ically more than 12) can be broadened into a flat con- [7] F. Couny, F. Benabid, P. J. Roberts, P. S. Light, and tinuum. This dramatically increases the width of the M. G. Raymer, Generation and photonic guidance of continuum to cover multiple octaves. For the fiber pa- multi-octave optical-frequency combs, Science 318, 1118 rameters, gas species and pressure that we use, the pump (2007). pulses cannot drive Kerr-based broadening mechanisms, [8] H. Kawano, Y. Hirakawa, and T. Imasaka, Generation such as self-phase modulation, soliton dynamics or mod- of more than 40 rotational raman lines by picosecond ulational instability, by themselves. Instead, it is only and femtosecond ti:sapphire laser for fourier synthesis, Applied Physics B 65, 1 (1997). through the sequential effects of vibrational comb forma- [9] H.-S. Chan, Z.-M. Hsieh, W.-H. Liang, A. H. Kung, C.-K. tion and consequent breakup of the long pump pulse into Lee, C.-J. Lai, R.-P. Pan, and L.-H. Peng, Synthesis and few-femtosecond pulses, which then proceed to drive both Measurement of Ultrafast Waveforms from Five Discrete rotational Raman scattering and self-phase modulation, Optical Harmonics, Science 331, 1165 (2011). that the continuum is created. [10] T. Suzuki, M. Hirai, and M. Katsuragawa, Octave- We expect our observations to be a general result for spanning raman comb with carrier envelope offset con- other molecular gases too, albeit with different charac- trol, Phys. Rev. Lett. 101, 243602 (2008). [11] S. E. Harris, D. R. Walker, and D. D. Yavuz, Raman teristics. This scheme could be expanded for VUV and technique for single-cycle pulses, Phys. Rev. A 65, 021801 mid-infrared supercontinuum generation by altering the (2002). pump wavelength and using a fiber with suitably tuned [12] N. Zhavoronkov and G. Korn, Generation of single in- resonances. Furthermore, this technique can be scaled tense short optical pulses by ultrafast molecular phase to high average power by increasing the laser repetition modulation, Phys. Rev. Lett. 88, 203901 (2002). rate. [13] F. Benabid, G. Bouwmans, J. C. Knight, P. S. J. Russell, and F. Couny, Ultrahigh efficiency laser wavelength con- Acknowledgements—This work is supported by version in a gas-filled hollow core photonic crystal fiber the National Natural Science Foundation of China by pure stimulated rotational raman scattering in molec- (No. 61675011, 61827820), Beijing Nova Program ular hydrogen, Phys. Rev. Lett. 93, 123903 (2004). (No. Z181100006218097), the Research program [14] J. Zheng and M. Katsuragawa, Freely designable opti- of Beijing Municipal Education Commission (No. cal frequency conversion in Raman-resonant four-wave- KZ201810005003), and the State Key Laboratory mixing process, Scientific Reports 5, 8874 (2015). of Advanced Optical Communication Systems and [15] Y. Y. Wang, F. Couny, P. S. Light, B. J. Mangan, and F. Benabid, Compact and portable multiline uv and vis- Networks, China. JCT, FB and CB are funded by ible raman lasers in hydrogen-filled HC-PCF, Opt. Lett. the European Research Council (ERC) under the 35, 1127 (2010). European Union’s Horizon 2020 research and innova- [16] A. Abdolvand, A. M. Walser, M. Ziemienczuk, tion program: Starting Grant agreement HISOL, No. T. Nguyen, and P. S. J. Russell, Generation of a phase- 679649. This work used EPCC’s Cirrus HPC Service locked raman frequency comb in gas-filled hollow-core (https://www.epcc.ed.ac.uk/cirrus). photonic crystal fiber, Opt. Lett. 37, 4362 (2012). 6

[17] F. Tani, F. Belli, A. Abdolvand, J. C. Travers, and P. S. ature and density dependence of the linewidths and line Russell, Generation of three-octave-spanning transient shifts of the rotational raman lines in N2 and H2, Phys. raman comb in hydrogen-filled hollow-core PCF, Opt. Rev. A 34, 1944 (1986). Lett. 40, 1026 (2015). [35] M. G. Raymer and J. Mostowski, Stimulated raman scat- [18] A. Benoˆıt, B. Beaudou, M. Alharbi, B. Debord, tering: Unified treatment of spontaneous initiation and F. Gérôme, F. Salin, and F. Benabid, Over-five oc- spatial propagation, Phys. Rev. A 24, 1980 (1981). taves wide raman combs in high-power picosecond-laser [36] J. K. Wahlstrand, S. Zahedpour, Y.-H. Cheng, J. P. pumped H2-filled inhibited coupling kagome fiber, Opt. Palastro, and H. M. Milchberg, Absolute measurement Express 23, 14002 (2015). of the ultrafast nonlinear electronic and rovibrational re- [19] Z. Cui, M. Chaturvedi, B. Tian, J. Ackert, F. Turner, and sponse in H2 and D2, Phys. Rev. A 92, 063828 (2015). D. Strickland, Spectral red-shifting of multi-frequency ra- [37] J. P. Palastro, T. M. Antonsen, and H. M. Milchberg, man orders, Optics Communications 288, 118 (2013). Compression, spectral broadening, and collimation in [20] A. V. Husakou and J. Herrmann, Supercontinuum gener- multiple, femtosecond pulse filamentation in atmosphere, ation of higher-order solitons by fission in photonic crys- Phys. Rev. A 86, 033834 (2012). tal fibers, Phys. Rev. Lett. 87, 203901 (2001). [38] J. K. Wahlstrand, J. H. Odhner, E. T. McCole, Y.-H. [21] K. Tai, A. Hasegawa, and A. Tomita, Observation of Cheng, J. P. Palastro, R. J. Levis, and H. M. Milchberg, modulational instability in optical fibers, Phys. Rev. Effect of two-beam coupling in strong-field optical pump- Lett. 56, 135 (1986). probe experiments, Phys. Rev. A 87, 053801 (2013). [22] M. Kolesik and J. V. Moloney, Nonlinear optical pulse [39] Y.-H. Chen, S. Varma, A. York, and H. M. Milchberg, propagation simulation: From maxwell’s to unidirec- Single-shot, space- and time-resolved measurement of ro- tional equations, Phys. Rev. E 70, 036604 (2004). tational wavepacket revivals in H2, D2, N2, O2, and N2O, [23] J. C. Travers, T. F. Grigorova, C. Brahms, and F. Belli, Opt. Express 15, 11341 (2007). High-energy pulse self-compression and ultraviolet gener- [40] E. T. J. Nibbering, G. Grillon, M. A. Franco, B. S. Prade, ation through soliton dynamics in hollow capillary fibres, and A. Mysyrowicz, Determination of the inertial contri- Nature Photonics 13, 547 (2019). bution to the nonlinear refractive index of air, N2, and [24] E. Marcatili and R. Schmeltzer, Hollow metallic and di- O2 by use of unfocused high-intensity femtosecond laser electric waveguides for long distance optical transmis- pulses, J. Opt. Soc. Am. B 14, 650 (1997). sion and lasers, Bell System Technical Journal 43, 1783 [41] G. Mikhatlov, The influence of temperature and pressure (1964). on the raman spectrum of nitrogen, Soviet Physics JETP [25] A. Börzsönyi, Z. Heiner, M. P. Kalashnikov, A. P. 36 (1959). Kovács,and K. Osvay, Dispersion measurement of in- [42] K. S. Jammu, G. E. S. John, and H. L. Welsh, Pressure ert gases and gas mixtures at 800 nm, Appl. Opt. 47, broadening of the rotational raman lines of some simple 4856 (2008). gases, Canadian Journal of Physics 44, 797 (1966). [26] M. Geissler, G. Tempea, A. Scrinzi, M. Schnürer, [43] F. Couny, F. Benabid, and P. S. Light, Subwatt thresh- F. Krausz, and T. Brabec, Light propagation in field- old cw raman fiber-gas laser based on H2-filled hollow- ionizing media: Extreme nonlinear optics, Phys. Rev. core photonic crystal fiber, Phys. Rev. Lett. 99, 143903 Lett. 83, 2930 (1999). (2007). [27] A. M. Perelomov, V. S. Popov, and M. V. Terent’ev, [44] F. Belli, A. Abdolvand, W. Chang, J. C. Travers, and Ionization of atoms in an alternating electric field, Soviet P. S. Russell, Vacuum-ultraviolet to infrared supercon- Journal of Experimental and Theoretical Physics 23, 924 tinuum in hydrogen-filled photonic crystal fiber, Optica (1966). 2, 292 (2015). [28] F. Tani, F. Köttig,D. Novoa, R. Keding, and P. S. Rus- [45] F. Tani, J. C. Travers, and P. S. J. Russell, Phz-wide su- sell, Effect of anti-crossings with cladding resonances on percontinua of nondispersing subcycle pulses generated ultrafast nonlinear dynamics in gas-filled photonic crys- by extreme modulational instability, Phys. Rev. Lett. tal fibers, Photon. Res. 6, 84 (2018). 111, 033902 (2013). [29] H. Lehmeier, W. Leupacher, and A. Penzkofer, Nonreso- [46] S. fei Gao, Y. ying Wang, X. lu Liu, C. Hong, S. Gu, nant third order hyperpolarizability of rare gases and N2 and P. Wang, Nodeless hollow-core fiber for the visible determined by third harmonic generation, Optics Com- spectral range, Opt. Lett. 42, 61 (2017). munications 56, 67 (1985). [47] F. Benabid, G. Bouwmans, J. C. Knight, P. S. J. Russell, [30] J. Bendtsen, The rotational and rotation-vibrational ra- and F. Couny, Ultrahigh efficiency laser wavelength con- man spectra of 14N2, 14N15N and 15N2, Journal of Ra- version in a gas-filled hollow core photonic crystal fiber man Spectroscopy 2, 133 (1974). by pure stimulated rotational raman scattering in molec- [31] M. Koszykowski, L. Rahn, R. Palmer, and M. Coltrin, ular hydrogen, Phys. Rev. Lett. 93, 123903 (2004). Theoretical and experimental studies of high-resolution [48] A. Abdolvand, A. Nazarkin, A. V. Chugreev, C. F. inverse raman spectra of molecular nitrogen at 1-10 atm, Kaminski, and P. S. Russell, Solitary pulse generation by Journal of Physical Chemistry 91, 41 (1987). backward raman scattering in H2-filled photonic crystal [32] W. F. Murphy, W. Holzer, and H. J. Bernstein, Gas phase fibers, Phys. Rev. Lett. 103, 183902 (2009). raman intensities: A review of “pre-laser” data, Applied [49] A. Nazarkin, A. Abdolvand, A. V. Chugreev, and Spectroscopy 23, 211 (1969). P. St.J. Russell, Direct observation of self-similarity in [33] J. K. Wahlstrand, Y.-H. Cheng, and H. M. Milchberg, evolution of transient stimulated raman scattering in Absolute measurement of the transient optical nonlin- gas-filled photonic crystal fibers, Phys. Rev. Lett. 105, earity in N2,O2,N2O, and Ar, Phys. Rev. A 85, 043820 173902 (2010). (2012). [50] S. T. Bauerschmidt, D. Novoa, B. M. Trabold, [34] G. C. Herring, M. J. Dyer, and W. K. Bischel, Temper- A. Abdolvand, and P. S. Russell, Supercontinuum up- 7

conversion via molecular modulation in gas-filled hollow- Nagar, Y. Liu, J. Nesper, O. Suarez, F. Rivas, Y. Wu, core pcf, Opt. Express 22, 20566 (2014). B. Shim, and M. Chini, Multioctave supercontinuum [51] B. Debord, M. Maurel, F. Gerome, L. Vincetti, generation and frequency conversion based on rotational A. Husakou, and F. Benabid, Strong nonlinear optical nonlinearity, Science Advances 6, (2020). effects in micro-confined atmospheric air, Photon. Res. [54] S. A. Mousavi, H. C. H. Mulvad, N. V. Wheeler, P. Ho- 7, 1134 (2019). rak, J. Hayes, Y. Chen, T. D. Bradley, S. ul Alam, S. R. [52] F. Belli, A. Abdolvand, J. C. Travers, and P. S. J. Russell, Sandoghchi, E. N. Fokoua, D. J. Richardson, and F. Po- Control of ultrafast pulses in a hydrogen-filled hollow- letti, Nonlinear dynamic of picosecond pulse propagation core photonic-crystal fiber by raman coherence, Phys. in atmospheric air-filled hollow core fibers, Opt. Express Rev. A 97, 013814 (2018). 26, 8866 (2018). [53] J. E. Beetar, M. Nrisimhamurty, T.-C. Truong, G. C.