Raman scattering has found renewed interest owing to the development of gas-filled hollow-core fibers, which constitute a unique platform for exploration of novel ultrafast nonlinear phenomena beyond conventional solid-core-fiber and free-space systems. Much progress has been made through models for particular interaction regimes, which are delineated by the relation of the excitation pulse duration to the time scales of the Raman response. However, current experimental settings are not limited to one regime, prompting the need for tools spanning multiple regimes. Here, we present a theoretical framework that accomplishes this goal. The theory allows us to review recent progress with a fresh perspective, makes new connections between distinct temporal regimes of Raman scattering, and reveals new degrees of freedom for controlling Raman physics. Specific topics that are addressed include transient Raman gain, the interplay of electronic and Raman nonlinearities in short-pulse propagation, and interactions of short pulses mediated by phonon waves. The theoretical model also accommodates vector effects, which have been largely neglected in prior works on Raman scattering in gases. The polarization dependence of transient Raman gain and vector effects on pulse interactions via phonon waves is investigated with the model. Throughout this Perspective, theoretical results are compared to the results of realistic numerical simulations. The numerical code that implements the new theory is freely available. We hope that the unified theoretical framework and numerical tool described here will accelerate the exploration of new Raman-scattering phenomena and enable new applications.
I. INTRODUCTION
Raman scattering is a type of inelastic scattering in which photons undergo frequency up- or down-conversion through interactions with material vibrational or rotational excitations (phonons). It was predicted by Smekal in 19231 and observed experimentally in liquids by Raman and Krishnan in 19282–4 and independently in quartz crystals by Mandelstam and Landsherg during the same period.5–7 Raman was awarded the 1930 Nobel prize, although the full and undivided recognition of his contribution has been contentious.8–11 There was even a suggested designation of the phenomenon as “the effect of Raman, Mandelstam, and Landsherg.”12 Within two years, its discovery had led to numerous subsequent observations in 60 different liquids and gases, which were recognized as providing support for the correctness of the then-new quantum theory.13 Moreover, it has given rise to a diverse range of applications. For example, Raman spectroscopy enables identification and characterization of various materials,14–17 and Raman amplifiers and Raman lasers extend wavelengths beyond inherent limitations of natural lasing media.18–25
With the advent of hollow-core fibers (HCFs), studies of optical interactions with atomic or molecular gases have attracted much attention.26–30 In particular, Raman scattering in molecular gases31–35 has experienced a resurgence in popularity since the pioneering work by Benabid et al.36 HCFs offer a substantial reduction of the Raman threshold through the combination of long interaction length and high intensity. Raman frequency down-conversion in H2,37–45 D2,46–52 CH4,53–59 CO2,60,61 and N262 has been demonstrated to generate a wide range of colors. Similarly, frequency up-conversion has been achieved with control of wave-vector matching63–65 and successfully applied in the generation of quantum-state-preserving photons.66 Further tunability of the converted wavelengths has been realized through the Raman-induced soliton self-frequency shift (SSFS)67–69 or its combination with photoionization-induced blue-shifting of the spectrum.70,71 Raman-enhanced supercontinuum generation has been demonstrated to cover from vacuum ultraviolet to near-infrared wavelengths,72–76 which makes possible pulse compression down to few-cycle durations.77–82
The analysis of stimulated Raman scattering (SRS) is commonly classified into multiple regimes based on the relationship of the excitation pulse duration to the material response time (Fig. 1). The most extensively examined regime is the “steady-state” regime, in which the pulse duration (△tp) is significantly longer than the phonon dephasing time T2.25,83–85 In this regime, the Raman response of a medium depends on each temporal segment of a pulse, so the Raman interaction is determined by the pulse’s instantaneous intensity. Moreover, the independence of the temporal segments leads to an incoherent Raman pulse that originates from noise. If the pulse duration becomes shorter than the dephasing time, the interaction enters the “transient” regime. This regime allows long-lasting phonon interactions to drive the SRS process within a pulse. When the Raman gain is high and saturated, “Raman memory” enables coherent generation of Raman pulses that tend to exhibit the same temporal phase profile as the pump pulse.53 Ultrashort femtosecond Raman pulses can also be produced without Raman spectral narrowing by scattering the pulse with existing phonons.86–90 If the pulse is even shorter than the phonon oscillation period TR {equivalently, the inverse of the Raman transition frequency }, phonons are “impulsively” excited91 and persist in the medium after the pulse has left. The excited phonons allow for nonlocal interactions between the first and subsequent excitation pulses, which leads to controllable nonlinear dynamics, such as blue-shifting, red-shifting, and pulse compression, with varying delays between the pulses.92–95 By probing the nuclear motions (phonons), the ultrafast dynamics of molecules can be monitored, which underlies time-resolved Raman spectroscopy.96–100 In addition, so-called time crystals formed by periodic phonon waves establish an analogy with condense-matter physics, where Bloch oscillations and Zener tunneling are observed.101 Generalization of the “phonon” concept from solid-state physics to refer to the “coherence wave” (the off-diagonal term of the density matrix of molecular motion) in molecular gases was introduced by Russell and co-workers.42,66,102,103 The use of “phonon” to refer to the coherence wave also minimizes potential confusion in discussions about the “coherent” generation of a Raman process, which pertains to interactions involving temporal phases between pump and Raman pulses. Details of the generalization and justification are presented in Sec. 2 of the supplementary material.
Various Raman regimes depending on the temporal relation between the pulse (blue lines) and the Raman response (orange lines).
Various Raman regimes depending on the temporal relation between the pulse (blue lines) and the Raman response (orange lines).
To establish a complete physical picture of SRS, it will be valuable to develop a unified theory that applies across temporal regimes. Despite the long history of Raman scattering and its widespread use in both the scientific and industrial communities, such a theory is still lacking. As an example of where it will be useful, we highlight techniques for ultrashort Stokes-pulse generation based on the interactions of multiple pulses with phonons.53,86–90 In this process, it is essential for the pump pulse duration to be long enough to prevent distortions and avoid competing nonlinear processes. Only with the unified theory can the interaction of impulsively excited phonons and a transient Raman process be fully understood. As another example, most prior theoretical works on Raman scattering neglect other effects that can play critical roles.102,104–110 A delicate balance between phonon annihilation (anti-Stokes processes) and creation (Stokes processes) can occur through wave-vector matching of the interaction of Stokes, pump, and anti-Stokes waves. In this situation, no new phonons are generated, resulting in no new photons as well. This phenomenon is referred to as “Raman gain suppression,” where both wave-vector matching and Kerr-induced four-wave mixing (FWM) play a non-negligible role. While its impact has been thoroughly explored in the steady-state regime,85,111–114 it has only been examined in the vicinity of perfect wave-vector matching in the transient regime.115 Even recent investigations into transient Raman gain suppression in HCFs rely on the steady-state gain model,90,102,103,110 which necessitates an explanation for its effectiveness. Unlike the case of steady-state Raman gain with FWM, a straightforward analytical expression for transient Raman gain with FWM is currently lacking.
In addition to the absence of a theory that can address multiple temporal regimes, there is a deficiency in our understanding of vector (i.e., polarization) effects in SRS. This knowledge gap can lead to uncontrolled physical phenomena, as exemplified by the depolarization observed during the SSFS process in a H2-filled anti-resonant fiber.69 Almost all prior studies employ scalar models. Vector aspects have been addressed within the steady-state formalism,84,116–121 but a model that can elucidate the temporal dynamics of vector Raman interactions does not exist. Typical Raman-scattering mechanisms involve transitions between vibrational and rotational states. Vibrational SRS in solids and liquids has been the focus of many prior studies, and its vector aspects in crystalline materials have been addressed extensively.122–126 However, vibrational SRS is isotropic in gaseous environments. On the other hand, rotational SRS exhibits pronounced anisotropy owing to its inherent connection with exchange of angular momentum, and, as a result, exhibits distinct responses to light with varying polarizations. Only recently has a tensor formulation been introduced to simulate rotational SRS in air to mitigate laser-plasma instabilities through nonlinear spectral broadening with elliptically polarized light.127 General aspects of vector SRS interactions are still not fully understood.
This Perspective aims to address these gaps in our current understanding of SRS, with a focus on processes relevant to Raman generation with picosecond- or femtosecond-duration pulses in gas-filled hollow waveguides. Although steady-state Raman gain has been investigated extensively, transient Raman gain and its associated dynamics have not been thoroughly examined. Previous studies have focused on Raman effects influenced by electronic-induced nonlinear phase modulations.128–130 However, the potentially significant contribution of the Raman response to nonlinear phases has been overlooked. Although nonlocal interactions through phonon waves have been extensively investigated,91–95,101,131–143 a complete picture only emerges through analysis that can handle different temporal regimes. Finally, the vector Raman response of gases has largely been neglected.
Here, we present a theoretical model of pulse propagation in gas-filled HCFs that covers all temporal regimes (steady-state, transient, and impulsive) of SRS as well as vector effects. Application of the theory to previous experiments will serve as a partial review of recent activities in the area and set the context for future directions. The theory aids conceptual and intuitive understanding of previously observed phenomena, in some cases from new points of view, and makes predictions about processes that are currently under investigation, or will be in the future. The results also illustrate new approaches to controlling Raman scattering for wavelength-conversion applications. Analytic expressions provide qualitative and quantitative descriptions of example phenomena, and detailed numerical simulations with representative experimental parameters back up and complement the analytical results. Although the presented model is capable of treating spontaneous Raman scattering (and all numerical results presented in this Perspective were calculated with the inclusion of spontaneous Raman scattering), the quantum nature of the initiation process does not play a key role in the processes examined here, which are in the nonlinear regime of Raman scattering (i.e., when the gain is high and saturated). Thus, we do not discuss the initiation process. For treatments of spontaneous processes and the initiation of Raman scattering, readers are referred to excellent prior works.13,107,144–149
The rest of this Perspective is organized as follows. Sec. II presents the nonlinear pulse propagation equation, which is the foundation of the unified theory, numerical simulations, and subsequent analysis presented in this Perspective. Section III applies the scalar model to phenomena in all three temporal regimes, with emphases on the behavior of the Raman gain and the interplay of Raman and electronic nonlinearities that govern a variety of ultrafast phenomena. In Sec. IV, the vector Raman response is derived and applied to polarization effects in vibrational and rotational SRS in diatomic molecules. Section V offers some perspectives on possible future directions for the field of SRS physics in gas-filled HCF.
The computational process is described in Secs. 11–13 of the supplementary material, along with some technical issues that are encountered in simulations of SRS in gases. The numerical code used to produce all of the simulation results presented in this Perspective is implemented in the MATLAB computing platform and is freely available.
II. PULSE PROPAGATION EQUATION
Raman scattering is generally analyzed with one of the two models. One model is based on the assumption of a harmonic oscillator for the excited vibrational mode,85,139–141 while the other treats the density matrix for the vibrational or rotational states.68,73,90,101,103,105,107–110,115,142,150–153 Although they are not discussed in this Perspective, those models underlie our theoretical approach. Readers are referred to Sec. 3 of the supplementary material for detailed descriptions and comparisons of the two models.
Shot noise of the input field is modeled semi-classically by including one noise photon per spectral discretization bin.162,169–171 Section 4 of the supplementary material includes more details about its numerical implementation. It is worth mentioning that there are other input-pulse noise models, obtained by either carefully modeling the buildup of lasing of the source or adding a Lorentzian noise spectrum derived from a phase-diffusion model.170 However, they aim to realistically model the input-pulse noise spectrum. Gases usually have large Raman transition frequencies, and such relatively narrowband noise spectra have insignificant impact on nonlinear processes. Thus, more-realistic noise models are unnecessary for the purposes of this work.
Two types of hollow-core fibers are considered below. In HCF, the dispersion has components of anomalous waveguide dispersion and normal gas dispersion.159 An anti-resonant fiber with a small 30-μm core diameter (and 300-nm tube-wall thickness) is assumed for analysis of processes with anomalous dispersion, while a Ag-coated capillary with a large 300-μm core diameter is assumed for those with normal dispersion.
Numerical simulations assume that the HCFs are filled with H2, N2, or O2. These gases are currently the most prevalent Raman-active gases in use. Numerical values of their parameters, as well as initial validation of our model, are provided in Secs. 11 and 12 of the supplementary material.
III. SCALAR RAMAN RESPONSE WITH FWM
In this section, we aim to present Raman effects across different regimes using a unified model based on a scalar UPPE. Several known phenomena will be reexamined with this unified framework. It covers not only Raman responses to Stokes generation but also coupling with anti-Stokes waves arising from both electronic- and Raman-induced FWM. FWM underlies Raman gain suppression, Raman-induced phase modulations, and other related effects, which will be discussed later.
Raman gain is a fundamental quantity in SRS physics, so Raman gain suppression is a critical concept. Since the groundbreaking research conducted by Bloembergen and Shen on Raman gain suppression,111 there has been a significant focus on understanding parametric effects on Raman gain.112–114,116–118,120,121,172–177 The small dispersion of gases translates into small wave-vector mismatch among Stokes, pump, and anti-Stokes waves, so the effect of Raman gain suppression (which is the strongest at perfect wave-vector matching) in gas-filled HCF cannot be ignored. Argon has been added to hydrogen to modify the overall dispersion profile, creating wave-vector mismatch and mitigating gain suppression,110 for example. Higher-order transverse modes have been introduced to excite intermodal coherence waves to disrupt the Stokes/anti-Stokes balance and allow the Stokes wave to grow.42 In this paper, we focus mostly on the Raman gain of Stokes waves; for anti-Stokes evolution, please see Sec. 8 of the supplementary material for details.
In gases, the Raman response is narrowband, and frequencies detuned from ωR do not play an important role. We assume △ω = ωR to simplify the subsequent analysis. The integral is responsible for the Raman-induced index change (which underlies phase modulations), and is usually responsible for the Raman gain, where i*, j*, and k* represent the complex conjugates of the corresponding fields. Equations (7) and (8) are the basis of the analysis of distinct temporal regimes below.
A. Steady-state regime (△tp ≫ T2)
We start with the steady-state regime. While this regime has been extensively explored, it serves as a reference point for subsequent discussions on the less-explored transient and impulsive regimes. Additionally, most prior studies ignore the electronic response.85,102,103,105,110 The few studies that include the electronic nonlinearity treat silica fiber, where SRS is weak.172,173,178 A thorough investigation of strong Raman response plus electronic nonlinearity in gases is still lacking.
Steady-state Raman gain g [Eq. (11)] with respect to wave-vector mismatch △β. The arrow indicates the trend of increasing (blue → red → black → green → pink) with fixed . The thick black line is the line with equal magnitudes of real and imaginary parts of .
Steady-state Raman gain g [Eq. (11)] with respect to wave-vector mismatch △β. The arrow indicates the trend of increasing (blue → red → black → green → pink) with fixed . The thick black line is the line with equal magnitudes of real and imaginary parts of .
This section analytically calculates the steady-state Raman integrals [Eq. (9)] and Raman gain [Eq. (10)], incorporating FWM arising from electronic and Raman nonlinearities. The variation of Raman gain with the strength of the electronic and Raman responses is illustrated, ranging from negligible Raman response to the pronounced response in Raman-active gases (Fig. 2).
B. Transient regime (△tp ≪ T2, △tp ≫ TR)
Because gases typically possess dephasing times greater than ∼100 ps at moderate gas pressures, experiments performed with picosecond or femtosecond pulses are commonly in the transient Raman regime. Additionally, gases usually exhibit large rotational and vibrational Raman transition frequencies. Consequently, the condition of having a pulse duration longer than the phonon oscillation period is frequently met.
The biggest challenge in analyzing transient SRS is the calculation of the Raman integrals, and [Eq. (8)]. During pulse propagation, AP acquires a z-dependent phase through self-phase modulations (SPM), while Ai (i = S or AS) acquires it through cross-phase modulations (XPM) from the pump. These phase modulations are intensity- and, thus, time-dependent, which leave the Raman integrals analytically unsolvable. To simplify the problem, we assume that all three pulses have an identical flat-top temporal structure, leading to time-independent phase modulations within the pulse. In addition, it is well-known that, when the gain is high and saturated, both the Stokes and anti-Stokes signals exhibit the same temporal phase structure as the pump.105 This results from amplifying only those frequencies that can maximize the Raman integrals, which always include either or [Eq. (7)]. As a result, Ai = CiAP (i = S or AS, and |Ci| ≪ 1 with Ci being real-valued) is reasonably assumed. Furthermore, the flat-top assumption facilitates the analysis of temporal characteristics of the Raman processes. While we make the “flat-top” assumption here, the validity of the following discussions extends more generally to pulses with negligible nonlinear phase modulations in comparison to Raman responses.
The similarity of the transient equations [Eqs. (13) and (14)] to the steady-state equations [Eqs. (9) and (10)] allows several phenomena in the steady-state regime to be transferred to the transient regime, such as the Raman gain shape and the resulting Raman gain suppression at △β = 0. This explains why the expression for the steady-state Raman gain has proven useful in predicting the qualitative features of transient Raman gain suppression.42,90,102,103,110
Despite the apparent similarity of the steady-state and transient regimes, they exhibit some important differences. First, the connection of between the steady-state value and the transient value shows that the transient Raman gain is smaller than the steady-state gain by a factor of △tp/T2. This leads to the potential existence of a Raman gain peak for an ultrashort pulse due to the reduced or increased . In addition, monotonically increasing with time results in not only stronger gain at the pulse trailing edge but also a varying gain shape throughout the pulse [Fig. 3(a)]. There is also a significant difference regarding where the Raman growth is the strongest. The Raman gain in the steady-state regime depends on pulse instantaneous intensity, whereas it relies on integrated pulse energy in the transient regime. This leads to the strongest Raman generation at the pulse center in the steady-state regime [Fig. 3(b)] but at the trailing edge in the transient regime [Fig. 3(c)].
(a) Varying Raman gain (extracted from Fig. 2) at different positions of a pulse in the transient Raman regime. Left represents the leading edge. Temporal profiles of the pump (P) and the Stokes (S) waves from the rotational Raman-Stokes generation in (b) the steady-state and (c) the transient regimes are shown. (b) A 20-ns and 70-μJ and (c) a 10-ps and 1-μJ Gaussian pulses at 1030 nm are launched into a 1-m-long anti-resonant fiber, filled with 30-bar H2. At this pressure, the dephasing time is 100 ps. Despite the linearly polarized fields, the dominant process is S(1) rotational SRS181 due to the stronger vibrational Raman gain suppression. All pulses in (b) and (c) are normalized to a peak power equal to one.
(a) Varying Raman gain (extracted from Fig. 2) at different positions of a pulse in the transient Raman regime. Left represents the leading edge. Temporal profiles of the pump (P) and the Stokes (S) waves from the rotational Raman-Stokes generation in (b) the steady-state and (c) the transient regimes are shown. (b) A 20-ns and 70-μJ and (c) a 10-ps and 1-μJ Gaussian pulses at 1030 nm are launched into a 1-m-long anti-resonant fiber, filled with 30-bar H2. At this pressure, the dephasing time is 100 ps. Despite the linearly polarized fields, the dominant process is S(1) rotational SRS181 due to the stronger vibrational Raman gain suppression. All pulses in (b) and (c) are normalized to a peak power equal to one.
In transient Raman scattering with high and saturated gain, the Raman pulses tend to inherit the temporal phase of the pump pulse. This results in periodic temporal modulations of participating pulses,105,106,109 self-similarity in pulse evolutions,182–184 Kerr-induced Raman suppression, and Raman pulse compression. In this paragraph, we will specifically focus on the last two effects. Konyashchenko et al. have discussed the suppression of Raman generation,128 which arises from the nonlinear wave-vector mismatch of pump and Raman pulses that do not have flat-top intensity profiles. This Kerr-induced suppression of Raman generation from SPM and XPM occurs because different rates of spectral broadening detune the frequency difference between the pump and the Stokes pulses from the Raman transition frequency (△ω = ωP − ωS ≠ ωR), where Raman generation is most efficient. Mathematically, a differential nonlinear phase between the pump and Stokes pulses results in a reduction of the Raman temporal integral , which contains either (i = S or AS) or its complex conjugate. Since the pump pulse acquires a phase of through SPM and SRS (contributed by ) and the Stokes/anti-Stokes pulse acquires through XPM and SRS (contributed by and ), the difference is , assuming that κ and κe are frequency-independent. This differential nonlinear phase increment is eliminated when the expression within the square brackets is equal to zero, a condition that can only be met in Stokes generation. To avoid suppression of Raman generation by SPM or XPM for non-flat-top pulses, there must be no differential nonlinear phase accumulation between the pump and Stokes pulses. For vibrational SRS in hydrogen with ∼125-THz transition frequency, this condition can be met with a linearly polarized pump wavelength around 970–1000 nm, which almost aligns with that of established and readily available Yb-based laser systems. Neglect of the Raman-induced nonlinear phases that arise from and leads to the conclusion that Raman generation is most effective at 1205 nm.128 As will be discussed below, SRS in gases induces significant nonlinear phases, even larger than the electronic-induced nonlinear phases (Fig. 7) and is thus non-negligible. Here, we treat Raman suppression using a simplified model for the differential nonlinear phase. The most-accurate treatment of Raman gain suppression will come from direct incorporation of nonlinear-phase effects in the Raman gain equation [Eq. (14)]; that will require generalization of the theory beyond the assumption of “flat-top” pulses, which remains for future work.
The coherent nature of transient SRS can be illustrated by the simultaneous compression and Raman shifting of femtosecond-duration pulses. In numerical simulations, positively chirped pulses with a 200-fs transform-limited duration were launched into a H2-filled capillary. The pulse energy is fixed at 500 μJ for constant Raman gain (at huge |△β|); the varying chirped-pulse duration affects only the Kerr-induced nonlinear phase. Initially, the Stokes generation efficiency increases monotonically with increasing chirped duration due to weakening Kerr-induced Raman suppression [orange line in Fig. 4(a)]. When the gain is high and saturated, Raman pulses attempt to acquire the same phase as the pump to maximize the Raman integrals.105 This includes the acquisition of the nonlinear phase of the pump before the Kerr-induced suppression becomes significant. Spectral broadening through nonlinear phase accumulations can be used to compress pulses in a dispersive delay line. In this case, the SPM-induced nonlinear phase accumulated by the pump is transferred to the Stokes pulses, which can be compressed to durations about 3–15 times shorter than the input pump pulse [Fig. 4(b)]. In the transient Raman regime, Raman spectral narrowing occurs in the single-pulse approach, where only the trailing edge of the chirped pump pulse is transformed into Raman pulses and results in their narrower bandwidths than the pump’s. By chirping the pump appropriately, Raman spectral narrowing can be overcome by this pulse-compression effect. It has also been experimentally demonstrated by Konyashchenko et al.129,130 Figure 2 plots Raman-gain curves with a fixed real part and varying imaginary part of in the gain equation, which, in the transient regime, becomes . This is useful in visualizing the temporally varying Raman gain over a pulse [Fig. 3(a)]. On the other hand, it can be visualized differently with the varying real part and fixed imaginary part, which is the case of varying the chirp of a pulse here [Fig. 4(c)]. Because the real part of depends on peak power, whereas its imaginary part on integrated pulse energy [Eqs. (13) and (14)], reducing the chirped-pulse duration but fixing the pulse energy creates the FWM peak in the Raman gain [Fig. 4(d)]. Therefore, with a proper wave-vector mismatch, Raman gain becomes stronger with reducing duration and the resulting higher peak power. This peak-power-dependent transient Raman-gain phenomenon is observed in the highly chirped regime in Fig. 4(a), where the efficiency starts to drop with increasing chirped-pulse duration. Fluctuations in generation efficiency arise from input-pulse shot noise and spontaneous Raman scattering, with increased significance observed in regimes of weaker Raman generation. However, in our example, heightened fluctuations are only evident in the long-pulse regime. In the short-pulse regime, Stokes generation is efficient due to the prominent Raman gain peak (with proper wave-vector mismatch) before Kerr-induced spectral detuning of pump and Stokes waves becomes significant. This initially efficient Raman generation quickly amplifies the pulse beyond the linear Raman regime where spontaneous emission noise dominates. After the initial generation, the Stokes pulse stops growing due to increased spectral detuning. On the other hand, in the long-pulse regime, slow Stokes generation, despite achieving a rather high efficiency after a long propagation, renders a Stokes pulse susceptible to spontaneous emission noise. Figure 4(e) summarizes the effects on Raman generation of varying nonlinear phase modulations: in the short-pulse regime, Kerr-induced Raman suppression dominates; increasing the pulse duration reduces the Kerr-induced suppression effect, enabling Raman pulse compression, and potentially introduces a rising Raman gain (if wave-vector mismatch meets where the gain peak is); further increasing the duration sees constant or dropping Raman gain based on the amount of wave-vector mismatch. More details about peak-power dependence, pulse compression, and quantum fluctuations are provided in Sec. 9 of the supplementary material.
Nonlinear-phase effects on transient vibrational Stokes generation by chirped pulses with varying duration. Pulses at 1030 nm and with 500-μJ energy and 200-fs transform-limited pulse duration are launched into a 50-cm-long capillary, filled with H2 to 20-bar pressure. (a) Generation quantum efficiency (QE) and peak power of the dechirped output Stokes pulses, produced from the Raman process. Mean values (center lines) and 1σ standard deviations (shaded areas) are calculated from ten simulations for each duration. The x axis is plotted on logarithmic scale to expose the variation at small durations. (b) Temporal profile of the dechirped Stokes pulse with the highest peak power in (a). (TL: transform-limited pulse; D: dechirped pulse). (c) Transient Raman gain gtr of a flat-top pulse, which ignores Kerr-induced Raman suppression from intensity-dependent SPM and XPM, with respect to wave-vector mismatch △β [Eq. (14)]. The arrow indicates the trend of increasing (blue → red → black → green → pink) with fixed . (d) Varying Raman gain at the trailing edge of a pulse with varying chirped duration and a fixed pulse energy. (e) Summary of the influence of nonlinear phase modulations on the Raman process with different pulse durations in the transient regime.
Nonlinear-phase effects on transient vibrational Stokes generation by chirped pulses with varying duration. Pulses at 1030 nm and with 500-μJ energy and 200-fs transform-limited pulse duration are launched into a 50-cm-long capillary, filled with H2 to 20-bar pressure. (a) Generation quantum efficiency (QE) and peak power of the dechirped output Stokes pulses, produced from the Raman process. Mean values (center lines) and 1σ standard deviations (shaded areas) are calculated from ten simulations for each duration. The x axis is plotted on logarithmic scale to expose the variation at small durations. (b) Temporal profile of the dechirped Stokes pulse with the highest peak power in (a). (TL: transform-limited pulse; D: dechirped pulse). (c) Transient Raman gain gtr of a flat-top pulse, which ignores Kerr-induced Raman suppression from intensity-dependent SPM and XPM, with respect to wave-vector mismatch △β [Eq. (14)]. The arrow indicates the trend of increasing (blue → red → black → green → pink) with fixed . (d) Varying Raman gain at the trailing edge of a pulse with varying chirped duration and a fixed pulse energy. (e) Summary of the influence of nonlinear phase modulations on the Raman process with different pulse durations in the transient regime.
In the steady-state regime, these coherent phenomena play out differently. Both Raman-enhanced SPM and nonlinear-phase-induced Raman suppression are observable in both the steady-state and the transient regimes [Eqs. (9a) and (13a)]. However, in the steady-state regime, spectral broadening can only result from the Raman process attempting to maintain a constant frequency difference ωP − ωS = ωR, where the highest Raman gain is, rather than from the transfer of the pump’s nonlinear phase as in the transient regime. If a Raman pulse originates from noise and lacks coherence, such incoherent spectral broadening does not lead to pulse compression. Since the steady-state Raman gain depends on , determined only by materials, steady-state Raman gain in Raman-active gases has no FWM-induced gain peak and shows no varying gain shape with different chirped-pulse durations as in the transient regime.
C. Impulsive regime (△tp ≪ T2, △tp ≪ TR)
In both the steady-state and transient regimes, the nuclear motions of a medium exhibit a response analogous to the electronic nonlinearity due to either a short phonon decay time or fast molecular oscillations. As a result, they both exhibit Kerr-like characteristics [Eqs. (9a) and (13a)] in addition to Raman gains [Eqs. (9b) and (13b)]. However, the response of the medium is no longer stationary if the pulse is shorter than the molecular oscillation period TR. The medium response is delayed relative to the pulse (Fig. 5). In addition to the delayed response, the medium is impulsively excited, which creates long-lived phonon waves after the pulse and allows for nonlocal interactions between pulses (Fig. 6).
Normalized index change (orange) induced by rotational SRS (represented by △ɛ) of hydrogen in response to a pulse of varying pulse duration (blue). The oscillation period of the S(1) rotational transition is 56.8 fs. The index change is computed from , with being the analytic signal of the electric field.185
Normalized index change (orange) induced by rotational SRS (represented by △ɛ) of hydrogen in response to a pulse of varying pulse duration (blue). The oscillation period of the S(1) rotational transition is 56.8 fs. The index change is computed from , with being the analytic signal of the electric field.185
Impulsively excited index change in 1 bar of (a) H2 or (b) N2 due to rotational SRS. (c) Values of Rcoeff for 1 bar of (blue) H2, (black) N2, and (red) O2 under the linearly polarized scalar assumption. The “revivals” in N2 and O2 at large time delays arise from beating of phonon waves from different rotational transitions,186 whereas in H2, the S(1) transition at 17.6 THz dominates, which leads to a constantly oscillating sinusoidal-like index wave.
Impulsively excited index change in 1 bar of (a) H2 or (b) N2 due to rotational SRS. (c) Values of Rcoeff for 1 bar of (blue) H2, (black) N2, and (red) O2 under the linearly polarized scalar assumption. The “revivals” in N2 and O2 at large time delays arise from beating of phonon waves from different rotational transitions,186 whereas in H2, the S(1) transition at 17.6 THz dominates, which leads to a constantly oscillating sinusoidal-like index wave.
Due to the delayed response of the medium to the pulse, corresponds to the buildup of the Raman-induced index change, which leads to a time-dependent phase increment. This causes the pulse spectrum to red-shift .139 Traditionally, red-shifting is treated as a consequence of impulsive phonon-wave generation after the pulse, whereas here it is explained by the nonlinear dynamics resulting from the index change. Greater red-shifting at the trailing edge of the pulse induces the negative chirp.140 If the dispersion is anomalous, the pulse will stretch temporally. On the other hand, in the normal-dispersion regime, this nonlinear phase will compete with the positive chirp from dispersion and SPM, potentially leading to pulse compression. However, in initial simulations, we find that SPM-induced chirp dominates over Raman-induced nonlinear phases. More investigations into this effect are thus required.
The Raman-enhanced SPM effect in the transient regime, resulting from and , has been exploited to generate ultra-broadband continua and few-cycle pulses.73,75 However, with shorter pulses, the Raman-induced index changes undergo a transition from SPM and XPM in the transient regime [Eqs. (9a) and (13a)] to red-shifting in the impulsive regime [Eq. (18)], which diminishes the SPM enhancement. Figure 5 illustrates the transition from pulse-following (∝|A|2) to a rising index with reducing pulse duration. In H2, the transition occurs at around 50 fs, where the dominant S(1) rotational transition starts to become impulsive, while in N2, it occurs at 500 fs [Figs. 6(c) and 7]. Because N2 has more available Raman transitions, which are also at smaller frequencies [Fig. 6(c)], the Raman-induced SPM enhancement in N2 is about times as strong as in H2, leading to the final times stronger Raman-enhanced SPM. This effect has been studied with 30 and 280-fs pulses in H273 and N2,75 respectively, and greater enhancements should be achievable by using pulses longer than the transition durations. On the other hand, simultaneous spectral broadening with red-shifting can produce few-cycle pulses at long wavelengths.78 The optimal pulse duration for this process lies between the transient and impulsive regimes, where the distinctive characteristics of Kerr and red-shifting effects both come into play. It is worth noting that the Raman-induced SPM is 4.5 times larger than the electronic-induced SPM in N2, while they have nearly the same magnitude in H2, values that are consistent with measurements of the total and electronic-induced nonlinear refractive indices.186
Accumulation of a nonlinear phase from numerical simulations with varying pulse duration in (a) H2 and (b) N2. A linearly polarized Gaussian pulse at 800 nm and with a fixed 6.67-MW peak power is launched into a 2-m-long anti-resonant fiber pressurized to 1 bar. The peak power is chosen so that the nonlinearity is sufficiently weak to avoid further amplitude-related Raman-gain effects, and the dispersion is artificially set to zero to isolate the nonlinear effects. The nonlinear phase ϕmax is calculated from the spectral broadening of the pulse.187,188
Accumulation of a nonlinear phase from numerical simulations with varying pulse duration in (a) H2 and (b) N2. A linearly polarized Gaussian pulse at 800 nm and with a fixed 6.67-MW peak power is launched into a 2-m-long anti-resonant fiber pressurized to 1 bar. The peak power is chosen so that the nonlinearity is sufficiently weak to avoid further amplitude-related Raman-gain effects, and the dispersion is artificially set to zero to isolate the nonlinear effects. The nonlinear phase ϕmax is calculated from the spectral broadening of the pulse.187,188
It is common to model the Raman response of a medium effectively as a single damped harmonic oscillator, with amplitude determined by the “Raman fraction” fr of the total nonlinear response, which can be found from Fig. 7: .77–79,189,190 We are now in position to examine the validity of the model while shedding some insights on Raman processes. This model is especially useful if the Raman response can be approximated as a harmonic oscillator with strong damping, which mitigates the potential numerical aliasing resulting from long-lived phonon waves. This is the case in N2 and O2. Due to their reviving indices, sufficiently short pulses essentially experience only the first spike of index change [Fig. 6(b)], which can be approximated by a highly damped harmonic oscillator. When the pulse becomes long enough to be affected by further reviving indices, this simplified model fails. In other words, the model operates by smearing the clustered sharp spectral Raman lines around 2 THz into one broadband Raman spectrum [Fig. 6(c)] and will fail when the smearing effect becomes invalid.
Weaker Raman-induced nonlinear phase accumulation can be advantageous for some processes. An example is SSFS, which has been widely studied in solid glass fibers.191–197 Since the soliton number is proportional to γeff△tp,198 where γeff is the effective or total nonlinear coefficient, a pulse in the impulsive regime has a smaller soliton number than that in the transient regime due to the reduced pulse duration. A smaller soliton number mitigates energy loss during soliton fission and thus enhances the efficiency of generating the reddest Raman soliton during SSFS.199 Furthermore, the decrease in Raman enhancement of SPM leads to a reduction in γeff, which also decreases the soliton number. The decrease can be as much as a factor of 4.5 times in N2, for example [Fig. 7(b)]. Finally, the reduced nonlinear phase is transformed to more-pronounced red-shifting through the rising Raman-induced index change [Eq. (18)]. Together, these factors make the SSFS a compelling option for wavelength-shifting in gases when it can be driven in the impulsive regime.
The temporal behavior [Figs. 6(a) and 6(b)] of Raman-induced index change can play an important role in Raman processes. We continue with the example of the SSFS. In gases with ultra-narrowband Raman responses, only pulses in the impulsive regime have the required bandwidth for SSFS to occur. The Raman soliton continues to red-shift due to the rising Raman-induced index change. Moreover, the soliton duration increases and its bandwidth decreases due to the soliton area theorem.198 The SSFS significantly slows down when the Raman soliton becomes long enough that the Raman-induced index change starts to catch up with the pulse, and the nonlinear-phase-induced red-shifting ceases. As a result, the Raman-soliton duration is approximately bounded above by 1–2 times the duration of the first spike of the Raman-induced index change. In this situation, the Raman soliton can also become too narrowband for intrapulse SSFS to occur. In H2 where the dominant S(1) rotational transition has TR = 57.8 fs [Fig. 6(a)], the soliton duration quickly evolves to 63 fs when the index starts to follow the pulse [Fig. 8(a)]. Similarly in N2, where the first spike is 118-fs-long [Fig. 6(b)], the soliton duration quickly evolves to 183 fs [Fig. 8(b)]. A slow increase of the duration follows due to a weakly delayed Raman-induced index change. In H2, the increase of the duration slows down also at its transition duration between two Raman regimes, where SSFS basically stops [Fig. 7(a)]. However, in N2, it slows down at a duration smaller than its transition duration of 500 fs [Fig. 7(b)]. Since the duration of temporal spikes of Raman response is inversely correlated with the collective spectral width of clustered Raman responses (as in N2 and O2), the ability of pulse-following is determined by not only the spectral positions of clustered Raman responses but also their collective spectral width. Higher frequencies and larger collective spectral width of clustered Raman responses correspond to a sharper temporal spike of the overall Raman response, leading to a shorter soliton duration after SSFS. On the contrary, since only one S(1) rotational transition dominates in H2, there is no such thing as collective spectral width and only the Raman frequency determines the temporal behavior of the Raman response that affects the Raman-soliton duration.
Evolutions of the (full-width-at-half-maximum) pulse duration during the SSFS process in a 100-m-long anti-resonant hollow-core fiber, filled with 10-bar (a) H2 or (b) N2. Black dashed lines are where the increase of the duration drops to 1% for every 1-m propagation. Insets are temporal profiles of the Raman-induced index changes (△ɛ; orange) and Raman solitons (blue) when the increase of the soliton duration significantly slows down (after black dashed lines). The injected linearly polarized fundamental solitons at 1030 nm have 10-fs duration, with 350-nJ energy in H2 and 270-nJ energy in N2 due to different dispersion and electronic nonlinearities.
Evolutions of the (full-width-at-half-maximum) pulse duration during the SSFS process in a 100-m-long anti-resonant hollow-core fiber, filled with 10-bar (a) H2 or (b) N2. Black dashed lines are where the increase of the duration drops to 1% for every 1-m propagation. Insets are temporal profiles of the Raman-induced index changes (△ɛ; orange) and Raman solitons (blue) when the increase of the soliton duration significantly slows down (after black dashed lines). The injected linearly polarized fundamental solitons at 1030 nm have 10-fs duration, with 350-nJ energy in H2 and 270-nJ energy in N2 due to different dispersion and electronic nonlinearities.
Double-pulse approach to Raman generation that avoids Raman spectral narrowing. The first pulse excites the phonon waves, which scatter the second pulse starting with its leading edge so that the entire pulse undergoes Raman generation, rather than only the trailing edge that results in Raman spectral narrowing for a chirped pulse.
Double-pulse approach to Raman generation that avoids Raman spectral narrowing. The first pulse excites the phonon waves, which scatter the second pulse starting with its leading edge so that the entire pulse undergoes Raman generation, rather than only the trailing edge that results in Raman spectral narrowing for a chirped pulse.
Effects of phonon dynamics in the transient regime. Two 50-fs pulses of different wavelengths are chirped to 10-ps duration and are launched into a 1-m-long capillary filled with 20-bar H2. Their wavelengths are chosen to yield the process of interest, with the varying first-pulse wavelength and fixed second-pulse 2-μm wavelength. Because of the unique S-shaped dispersion curve of gas-filled HCF, pumps of different wavelengths create different phonon wave vector βph.64,90 Evolutions of photon numbers of Stokes (red), pump (black), and anti-Stokes (blue) waves under different wave-vector-matching conditions are displayed for the cases of weak (10 μJ) and strong (5 mJ) second-pulse energy; the first pulse is fixed at 2 mJ. (a) and (c) satisfy the phonon amplification condition, where the first pulse is at 0.7 μm, while (b) and (d) satisfy the phonon absorption condition, where the first pulse is at 1.09 μm. The second pulse is weak in (a) and (b), while it is sufficiently strong to induce its own SRS in (c) and (d). Mean values (center lines) and 1σ standard deviations (shaded areas) are calculated from ten repeated simulations. (e) and (f) The temporal profiles of the dechirped anti-Stokes pulses from (b); in (f), the pulse quality is degraded by back-conversion. (f) The result from one of the ten simulations, which fluctuates significantly due to the fluctuating phonon-beating pattern. (g) and (h) The temporal profiles of the dechirped Stokes pulses from (c) and (d), respectively. The Stokes pulse in (g) has more than twice the peak power as that in (h). Despite being not shown here, the pulse quality of the dechirped Stokes pulse in (a) is the same as that in (g) from (c) (TL: transform-limited pulse, D: dechirped pulse).
Effects of phonon dynamics in the transient regime. Two 50-fs pulses of different wavelengths are chirped to 10-ps duration and are launched into a 1-m-long capillary filled with 20-bar H2. Their wavelengths are chosen to yield the process of interest, with the varying first-pulse wavelength and fixed second-pulse 2-μm wavelength. Because of the unique S-shaped dispersion curve of gas-filled HCF, pumps of different wavelengths create different phonon wave vector βph.64,90 Evolutions of photon numbers of Stokes (red), pump (black), and anti-Stokes (blue) waves under different wave-vector-matching conditions are displayed for the cases of weak (10 μJ) and strong (5 mJ) second-pulse energy; the first pulse is fixed at 2 mJ. (a) and (c) satisfy the phonon amplification condition, where the first pulse is at 0.7 μm, while (b) and (d) satisfy the phonon absorption condition, where the first pulse is at 1.09 μm. The second pulse is weak in (a) and (b), while it is sufficiently strong to induce its own SRS in (c) and (d). Mean values (center lines) and 1σ standard deviations (shaded areas) are calculated from ten repeated simulations. (e) and (f) The temporal profiles of the dechirped anti-Stokes pulses from (b); in (f), the pulse quality is degraded by back-conversion. (f) The result from one of the ten simulations, which fluctuates significantly due to the fluctuating phonon-beating pattern. (g) and (h) The temporal profiles of the dechirped Stokes pulses from (c) and (d), respectively. The Stokes pulse in (g) has more than twice the peak power as that in (h). Despite being not shown here, the pulse quality of the dechirped Stokes pulse in (a) is the same as that in (g) from (c) (TL: transform-limited pulse, D: dechirped pulse).
Control of the Raman processes through the wave-vector matching relies on a fixed phonon wave vector βph. This assumption holds only when phonons are excited through Stokes generation with a long pulse in the transient regime, resulting in phonons with a fixed wave vector equal to the difference between the pump and the Stokes fields, [Fig. 11(a)]. If the excitation pulse is short enough to enter the impulsive regime, the phonon wave vector is no longer fixed, but is instead determined by the nonlinear red-shifting process [Fig. 11(b)]. Since the phonon wave vector is not constant, the wave-vector-matching condition cannot be consistently met throughout the evolution, which results in an uncontrolled scenario of Stokes and anti-Stokes generation from the second pulse. This knowledge not only directs us toward controllable Raman generation through nonlocal phonon interactions but also offers more information of various Raman processes previously investigated, such as the phonon amplification conducted by Bustard et al.202 If phonons are impulsively excited by an ultrashort pulse, the upcoming pulse for amplifying the phonons with transient SRS must be sufficiently strong to guarantee the Stokes generation.
Evolutions of phonon wave vector βph during the SRS process in two different Raman regimes. Black, red, and blue lines represent total wave vectors of pump, Stokes, and phonons, respectively. Since the total momentum should be conserved, phonon wave vector is the difference between the pump and the Stokes waves in the transient regime: βtotal = NPβP + NSβS + Nphβph, where NP + NS + Nph = Ntotal represents the number of photons/phonons. In the impulsive regime, no discrete Stokes wave is generated and the pump pulse shifts to the red with varying propagation constant: , where NP + Nph = Ntotal. In (b), βP ∼ cos(Ntotal − NP) is assumed to artificially create a varying βP to simulate the red-shifting of the pump. Because of the red-shifting process, phonon wave vector is not a constant, in contrast to the transient regime.
Evolutions of phonon wave vector βph during the SRS process in two different Raman regimes. Black, red, and blue lines represent total wave vectors of pump, Stokes, and phonons, respectively. Since the total momentum should be conserved, phonon wave vector is the difference between the pump and the Stokes waves in the transient regime: βtotal = NPβP + NSβS + Nphβph, where NP + NS + Nph = Ntotal represents the number of photons/phonons. In the impulsive regime, no discrete Stokes wave is generated and the pump pulse shifts to the red with varying propagation constant: , where NP + Nph = Ntotal. In (b), βP ∼ cos(Ntotal − NP) is assumed to artificially create a varying βP to simulate the red-shifting of the pump. Because of the red-shifting process, phonon wave vector is not a constant, in contrast to the transient regime.
The story is different when the second pulse is in the impulsive Raman regime. In this situation, there is no discrete generation of Stokes or anti-Stokes fields by the second pulse. However, control of the delay between the two pulses permits for the control of the index change experienced by the second pulse. This can lead to controlled red-shifting through the rising index, as observed with in the discussed impulsive regime [Eq. (18)], or blue-shifting for the falling index.92–95
IV. VECTOR PROPERTIES OF RAMAN RESPONSES
Vector (i.e., polarization) effects in SRS from molecules have been largely neglected. Vibrational SRS is almost isotropic and thus exhibits good consistency between experiments and scalar models.203,204 The isotropy results from in the perturbative regime, where and are shown in Fig. 12. (Details are in Sec. 5 of the supplementary material.) Rotational SRS in gases requires a vector model due to its anisotropic nature arising from the exchange of angular momentum. Although there are a few prior works deriving the scalar model for rotational SRS in gases, they neglect its anisotropic nature.76,186,190,205–208 Here, we aim to address this gap by introducing a vector UPPE for gases.
Illustration of spherical coordinates, where the diatomic molecular axis is aligned along and the field propagates along . The field is polarized in the xy-plane.
Illustration of spherical coordinates, where the diatomic molecular axis is aligned along and the field propagates along . The field is polarized in the xy-plane.
There is a long history of investigation of rotational SRS in gases. Since Yoshikawa and Imasaka proposed the idea of phase-locking an SRS-induced multi-frequency spectrum to produce femtosecond pulses,209 simultaneous generation of multiple orders of Stokes/anti-Stokes waves has received substantial attention.72,76,141,151–153,210–213 One of the promising approaches is through molecular phase modulation, which requires impulsive excitation of gas molecules through SRS.63,92,141,151–153,210,212–215 While vibrational SRS can potentially produce sub-femtosecond pulses, rotational SRS is usually preferred in experiments because the longer period makes it easier to reach the impulsive regime. In addition, it can fill the spectral gap between vibrational Stokes/anti-Stokes wavelengths to create a smooth supercontinuum76 for potential generation of a single ultrashort pulse. It has long been known that rotational SRS is sensitive to the polarization of the field. A circularly polarized field preferentially drives rotational transitions, while linearly polarized light preferentially drives vibrational transitions.84,203,211 How these two scattering processes compete under different polarization conditions, and how to quantify the outcomes, have remained unclear. Rotational SRS, particularly in N2, is commonly observed with a linearly polarized field,37–40,46,47,68,72,75–78,90,95,138,143,151,153,201,211,214 but how its anisotropic nature affects the polarization and the subsequent nonlinear dynamics is unknown.
A. Derivation of the vector UPPE
Derivation of the vector UPPE begins with the (real-valued) polarization , averaged over microscopic molecular orientations and summed over many molecules (Ng is the number density of gas molecules in 1/m3). Here, we consider only diatomic molecules to simplify the derivation (Fig. 12), but an analytic expression should also be possible for more-complicated molecules by following the same process.
The unperturbed part of Eq. (23) can be solved with a simple integral over solid angle. By observing that , the unperturbed polarization corresponds to the linear susceptibility, which is usually replaced by a Sellmeier equation.218–221
B. Analysis of vector Raman gain
where . The gain shape and its evolution with varying imaginary part of are basically the same as those found in the scalar scenario (Fig. 2), but translated to negative △β (Fig. 13). Perfect Raman gain suppression occurs at , which depends on pump wavelength and peak power.
Raman gain shape gcross-linear [Eq. (38)] vs wave-vector mismatch △β. The arrow indicates the trend of increasing (blue → red → black → green → pink) with fixed . The thick black line is the Raman gain with equal magnitudes of real and imaginary parts of the quantity .
Raman gain shape gcross-linear [Eq. (38)] vs wave-vector mismatch △β. The arrow indicates the trend of increasing (blue → red → black → green → pink) with fixed . The thick black line is the Raman gain with equal magnitudes of real and imaginary parts of the quantity .
The Raman gain for the cross-circularly polarized case behaves distinctly from all other cases: it is unaffected by FWM [Eq. (36a)] because only phase-modulation terms, , appear in its evolution equation. In contrast to the other Raman gains, cross-circularly polarized Raman generation cannot be FWM-suppressed. This effect has been experimentally observed117,118,121 and modeled with a different theoretical formalism that considers only the rotational Raman response.116,120 In gases, wave-vector mismatch is usually small and the Raman gain is generally reduced substantially. Without a means to suppress it, the cross-circularly polarized Raman gain can be 100 times stronger than the Raman gain for other combinations of polarizations (Fig. 14).
Transient S(1) rotational Raman gain of H2 at the pulse trailing edge for co-circularly polarized [Eq. (14)] and cross-circularly polarized [Eq. (40b)] pump and Stokes waves. A 1.5-mJ and 1-ps pulse is launched into a capillary under the same conditions as in Fig. 4. To calculate the Raman gain, the Raman-induced index change must include contributions from all Raman responses in the sum , which appears in and , while only one is required in due to the narrowband feature of Raman response functions. Suppression of the co-polarized gain occurs at 2.04 μm, where △β = 0. The ratio of the gain magnitudes is 90 at 1030 nm.
Transient S(1) rotational Raman gain of H2 at the pulse trailing edge for co-circularly polarized [Eq. (14)] and cross-circularly polarized [Eq. (40b)] pump and Stokes waves. A 1.5-mJ and 1-ps pulse is launched into a capillary under the same conditions as in Fig. 4. To calculate the Raman gain, the Raman-induced index change must include contributions from all Raman responses in the sum , which appears in and , while only one is required in due to the narrowband feature of Raman response functions. Suppression of the co-polarized gain occurs at 2.04 μm, where △β = 0. The ratio of the gain magnitudes is 90 at 1030 nm.
Overall, the gain behaves similarly to the steady-state gain but with the differences pointed out in the scalar case, such as the varying gain shape throughout the pulse and the coherent properties.
Cross-polarized SRS can play an important role in experiments, with behavior that deviates from the predictions of a scalar model. As an example, for a 300-μm-core capillary, filled with H2 to 20-bar pressure and exposed to a 1.5-mJ and 1-ps pump pulse at 1030 nm, the cross-linearly polarized rotational Raman gain is ten times stronger than the vibrational Raman gain [Fig. 15(a)]. This finding challenges the scalar prediction that vibrational SRS dominates. To further examine this phenomenon, we investigated the Raman gain for varying pulse duration and energy. Vibrational SRS exceeds the co-linearly polarized rotational SRS, as predicted by scalar calculations.90 However, the scalar calculation does not foresee the dominance of cross-polarized rotational SRS for wide ranges of pulse duration and energy [Fig. 15(b)]. The distinct Raman-suppression conditions for cross-linearly polarized (which depends on the pump peak power and wavelength) and co-polarized SRS (which occurs at △β = 0) underlie the possible control for dominant vibrational SRS [Fig. 15(c)] or cross-linearly polarized rotational SRS [Fig. 15(d)].
Vector aspects of Raman gain. (a) Transient Raman gain at the pulse trailing edge for S(1) rotational and Q(1) vibrational SRS co-polarized with a linearly polarized pump [Eq. (14)], as well as for S(1) rotational SRS cross-polarized with it [Eq. (40a)]. An x-polarized, 1.5-mJ and 1-ps pulse is launched into a capillary filled with H2 to a pressure of 20 bar. The inset shows the variation of co-polarized Raman gains. (b) Ratio of vibrational transient Raman gain to cross-polarized rotational transient Raman gain at 1030 nm. The vibrational and rotational gains are equal along the black dashed line. Stars indicate conditions where (c) the vibrational SRS (0.3 mJ, 8 ps) or (d) the cross-polarized rotational SRS (1.5 mJ, 1 ps) dominate. PSD: power spectral density (spectral intensity) of the pulse. The lengths of the capillary in (c) and (d) are 1 m and 10 cm, respectively, chosen to produce enough Raman generation for visualization.
Vector aspects of Raman gain. (a) Transient Raman gain at the pulse trailing edge for S(1) rotational and Q(1) vibrational SRS co-polarized with a linearly polarized pump [Eq. (14)], as well as for S(1) rotational SRS cross-polarized with it [Eq. (40a)]. An x-polarized, 1.5-mJ and 1-ps pulse is launched into a capillary filled with H2 to a pressure of 20 bar. The inset shows the variation of co-polarized Raman gains. (b) Ratio of vibrational transient Raman gain to cross-polarized rotational transient Raman gain at 1030 nm. The vibrational and rotational gains are equal along the black dashed line. Stars indicate conditions where (c) the vibrational SRS (0.3 mJ, 8 ps) or (d) the cross-polarized rotational SRS (1.5 mJ, 1 ps) dominate. PSD: power spectral density (spectral intensity) of the pulse. The lengths of the capillary in (c) and (d) are 1 m and 10 cm, respectively, chosen to produce enough Raman generation for visualization.
The generation of a cross-polarized Stokes or anti-Stokes wave can naturally produce depolarization of a pulse in propagation. For a linearly polarized pulse, the growth of the cross-linearly polarized Raman pulses is slow due to gain suppression. As a consequence, the linearly polarized pulse can reasonably maintain its polarization [Fig. 16(a)]. In contrast, the growth of cross-circularly polarized waves is unimpeded and is unaffected by FWM-induced gain suppression. In N2, with its numerous low-frequency transitions [Fig. 6(c)], SRS becomes an intrapulse effect and ultimately produces significant depolarization [Fig. 16(b)].
Depolarization effects in SRS. A 100-μJ and 300-fs pulse at 1030 nm is launched into a 2-m-long capillary, filled with N2 to a pressure of 3 bar. Output spectra of (a) linearly and (b) circularly polarized input pulses. PSD: power spectral density (spectral intensity) of the pulse. The inset shows the initial growth of the circularly polarized Stokes pulse in cross-polarized with the pump.
Depolarization effects in SRS. A 100-μJ and 300-fs pulse at 1030 nm is launched into a 2-m-long capillary, filled with N2 to a pressure of 3 bar. Output spectra of (a) linearly and (b) circularly polarized input pulses. PSD: power spectral density (spectral intensity) of the pulse. The inset shows the initial growth of the circularly polarized Stokes pulse in cross-polarized with the pump.
These equations illustrate that when excited by a circularly polarized pulse, the phonon waves exhibit an isotropic behavior in relation to the second pulse [Eq. (42b)]. In contrast, if the excitation pulse is linearly polarized, the phonon waves are anisotropic for the second pulse. In this situation, the second pulse experiences a twofold increase in the strength of the rotational phonon waves, along with a π phase delay, when it is co-polarized as opposed to when it is cross-polarized [Eq. (42a)]. In addition, the impact of rotational phonon waves is the weakest when they are excited by a circularly polarized field. The polarization dependence that results from these effects can be exploited to drive vibrational or rotational transitions. As an explicit example, an intense ultrashort pulse was launched into a H2-filled capillary to impulsively excite rotational phonon waves, followed by a weak delayed pulse. The results exhibit the twofold enhancement in rotational SRS for the second pulse if the excitation pulse is linearly polarized [Fig. 17(a)], with no significant difference if it is circularly polarized [Fig. 17(b)].
Vector effects on propagating rotational phonon waves. Two pulses are launched into a capillary filled with H2. The first pulse drives the S(1) rotational transition impulsively, and the second pulse scatters from the phonon waves. The energy of the first pulse is 50 μJ in (a) and 200 μJ in (b) to compensate for the four times weaker rotational SRS in circular polarization. The output spectrum of the second pulse is shown for the cases where the first pulse is (a) linearly polarized or (b) circularly polarized . The second pulse is either x-polarized or y-polarized. PSD: power spectral density (spectral intensity) of the pulse. In (b), the two polarizations do not have exactly equal magnitudes because both the Stokes and anti-Stokes waves grow from noise, and the wave-vector matching is not specifically controlled [Eq. (20)].
Vector effects on propagating rotational phonon waves. Two pulses are launched into a capillary filled with H2. The first pulse drives the S(1) rotational transition impulsively, and the second pulse scatters from the phonon waves. The energy of the first pulse is 50 μJ in (a) and 200 μJ in (b) to compensate for the four times weaker rotational SRS in circular polarization. The output spectrum of the second pulse is shown for the cases where the first pulse is (a) linearly polarized or (b) circularly polarized . The second pulse is either x-polarized or y-polarized. PSD: power spectral density (spectral intensity) of the pulse. In (b), the two polarizations do not have exactly equal magnitudes because both the Stokes and anti-Stokes waves grow from noise, and the wave-vector matching is not specifically controlled [Eq. (20)].
V. PERSPECTIVES
Femtosecond wavelength conversion is investigated to extend various ultrafast science applications beyond natural lasing wavelengths. As just one example, nonlinear microscopy requires high peak power at 1300 or 1700 nm for three-photon imaging.225,226 Raman generation in gas-filled HCF is a promising candidate for the generation of wavelengths from the ultraviolet to the mid-infrared regimes. In addition, it enables power scaling due to low nonlinearity of gases, high damage threshold, and low propagation loss.26,227 However, current research is mostly restricted to pulse durations from hundreds of picoseconds to a few nanoseconds, where SRS can dominate over competing nonlinearities. A relatively small number of works have targeted ultrafast (femtosecond) Raman generation.42,48,57,68,78,86–90,128–130,228 Detrimental competition from the Kerr nonlinearity can be avoided by appropriately stretching an ultrashort pulse in time (Fig. 4) and removing the frequency chirp from the generated Raman pulse, albeit with Raman spectral narrowing. Although femtosecond Raman generation without Raman spectral narrowing has been achieved in a two-pulse approach, it might also be possible with a simpler single-pulse approach through interference with a continuous-wave field.53 Additionally, Raman spectral narrowing in a single-pulse approach can be overcome by operating in the temporal regime where Raman-induced pulse compression occurs (Fig. 4). Gas-based optical parametric amplification (OPA) is another candidate for femtosecond wavelength conversion.229,230 Raman-enhanced SPM can potentially boost the OPA process. It also provides a new degree of freedom for controlling the process through the pulse duration, which, for example, stabilizes OPA against variations of pulse duration in a regime where increasing the duration leads to competition from the Kerr nonlinearity; the effects of lower peak power can be offset by stronger SPM enhancement.
With a deeper understanding of SRS in N2 and O2, gas-based nonlinear optics may advance toward air photonics, first envisioned in the terahertz regime.231,232 To date, most complex physical phenomena, such as phonon dynamics, have been studied in H2 due to its simple Raman response.63–65,87,90,92–94,101,153,210,214,233,234 On the other hand, investigations of N2 and O2 allow for progress toward photonics based on ambient-air-filled HCFs that can be handled the same as solid-core fibers. Prior studies in N2 or air include spectral broadening75–77,127,208 and the resulting ultrashort pulse generation,78,79 phonon-driven spectral control143,200 and transient Raman generation,228 soliton compression235,236 and self-frequency shift.189 Only recently has vibrational Stokes generation in N2 at around 1.3–1.4 μm been realized in the nanosecond regime.62 Optimization and extension of these studies into various temporal regimes with varying pulse energies, as well as under ambient pressure, will be interesting. As an example, pulse compression in ambient-air-filled hollow-core fiber covering a broad range of pulse duration and energy can potentially offer the advantages of a waveguide (e.g., alignment-free operation and compactness) over multipass cells237–242 and multiplate compressors.243–246 Most current gas-filled pulse compressors use inert gases to avoid SRS.239–242,247–252 However, with specific combinations of pulse energy and duration, it will be possible to exploit SPM enhancement from SRS or to generate ultrashort spatially stabilized pulses through the formation of Raman-induced multidimensional solitary states.79,80,82 It will be straightforward to extend the theoretical framework of this Perspective to these other gas-based nonlinear optical platforms.
The gas pressure is a degree of freedom with HCF. This enables novel physical phenomena based on longitudinally varying pressure, which is not realizable in free-space or solid-core-fiber platforms. Self-focusing can be prevented with a negative pressure gradient, and ionization-induced defocusing can be prevented with a positive pressure gradient.253,254 Moreover, pressure gradients can extend the physics of tapered solid-core fiber192,255 to gas-filled HCF. As an example, SSFS with a positive pressure gradient exhibits several advantages over constant pressure. A flattened anomalous dispersion profile during SSFS suppresses pulse temporal broadening through the soliton area theorem Ns = γ△tp/|β2|, if the pulse experiences constant anomalous dispersion (|β2|) despite red-shifting, due to higher gas pressure. Increasing nonlinearity can compensate the loss of the Raman process; furthermore, it can reduce the pulse duration by maintaining Ns = 1. These effects combined might help push SSFS performance beyond what can be achieved with constant pressure.68,69 Moreover, the excitation of higher-order modes, which has been observed in recent SSFS experiments,69 at the fiber input can be suppressed. It is worth noting that although peak-power-induced ionization can be suppressed by gradient pressures, scaling of average power is ultimately limited by thermal effects in Raman processes,256 which can be resolved by reducing the repetition rate or working with pulse bursts.
Regarding technological developments, gas-filled HCF may have impact beyond research laboratories if they are pressurized and sealed, so gas-handling apparatus is not required.257 Such fiberized cells may enhance SRS by splicing both ends to fiber Bragg gratings to form resonant cavities.258,259 With continuous-wave light, 99.99% of the output light at the Stokes wavelength has been obtained,259 and Stokes waves at 1.7 μm have been generated in H2260 or D2.50 Beyond the previous single-pass structure, an all-fiber gas Raman oscillator has been recently implemented to produce 1.8 W of continuous Stokes waves at 1693 nm through rotational SRS in H2.261 With further development of the fusion splicing technology,262–268 a simple and monolithic gas-based all-fiber system that supports the aforementioned SRS phenomena will be possible.
Recent years have witnessed a surge of effort in the field of quantum information science.269–273 Quantum frequency conversion (QFC) serves the role of connecting quantum nodes274,275 that operate at diverse wavelengths. It has been achieved through three-wave mixing in quadratic (χ(2)) nonlinear crystals276–280 and through Bragg-scattering FWM in solid-core photonic crystal fibers.281–283 Although solid-core fibers overcome several challenges in crystals, such as stringent phase-matching requirements and undesirable optical background, further improvements have recently been reached with a H2-filled hollow-core fiber.66 In hydrogen, long-lasting phonon waves are established through SRS, which can frequency up-convert a signal pulse with proper wave-vector matching. The theoretical framework presented above will facilitate the optimization of phonon-driven nonlocal interactions and competing Raman gains. Temporal overlap of the pump and signal pulses may be unnecessary, both frequency down- and up-conversions are possible by changing gas pressure,284 and controllable use of either rotational or vibrational SRS for different amounts of frequency shift is possible. To date, QFC in gases is limited to hydrogen, which offers both the huge transition frequency and simple Raman response [Fig. 6(c)];66,284 implementations with other gases may be desirable for various reasons. Ultimately, we foresee a promising future where ultra-tunable QFC with hollow-core fiber, filled with a wide range of gases, is achieved.
In optical quantum communication, it is crucial to create photonic quantum states with a controlled degree of entanglement and preserved coherence among the modes over long-distance transmission. Among the four degrees of freedom available for encoding a photon (polarization, two transverse spatial dimensions, and time/frequency), time-frequency encoding not only spans an unbounded high-dimensional Hilbert space but is also compatible with existing single-mode waveguide platforms. During studies of spontaneous Raman scattering, fluctuations of the generated Stokes pulse were explained as resulting from the number of excited statistically independent “coherent temporal modes,”146,147 or so-called “time-frequency Schmidt modes.”285 Recently, temporal modes have garnered increasing interest in quantum information science because they constitute an orthogonal broadband wave-packet/pulsed basis and enable the use of the time/frequency degree of freedom.286–301 Photon pairs are generated in various platforms, such as bulk nonlinear crystals,288 solid-core photonic crystal fibers,302–306 and tapered,307 birefringent,308 dispersion-shifted,309–311 or highly nonlinear312 fibers and waveguides.313,314 However, they are limited by spontaneous Raman scattering, which creates uncorrelated Raman photons.309–311,314 Liquid-filled315,316 and Xe-filled317,318 HCFs have been used for Raman-free photon-pair generation. As demonstrated for quantum frequency conversion by Tyumenev et al.,66 Raman interactions in gases are controllable and can exhibit (quantum-state-preserving) coherence in the transient regime. Raman interactions within gas-filled HCFs present an opportunity to avoid the imposition of adverse noisy effects on quantum communication and conversion. For example, phonon-driven processes might be considered for controllable coherent generation and conversion of temporal modes. Aforementioned ambient-air-filled or sealed pressurized HCF is a promising platform for the simplest operation. All-fiber operations based on both solid-core and hollow-core fibers might be possible with a fiberized quantum pulse gate and pulse shaper.293,300,301 As HCF proves its utility in achieving high-fidelity and low-latency single-photon transmission319 and photon-pair generation,317,318 the exploration of SRS in gases holds the promise of unveiling quantum applications based on temporal modes in the future.
In this Perspective, we have only considered forward SRS. However, a Stokes pulse can be generated in the opposite direction from the pump pulse. The counter-propagating Stokes pulse can not only extract energy throughout the pump pulse but also experiences strong temporal compression.320–322 This has been observed in several liquids,320,321,323–325 methane,326–328 and D2329 in free-space geometries. Investigations of backward SRS have also been conducted in solid-core silica fiber330–333 and H2-filled HCF.334–338 The transient self-similar nature of the evolution, which exhibits in forward SRS as well,182–184 and the corresponding periodically modulated Stokes pulse have been studied.335,336 Prior works have relied on a seeded process, but only recently has noise-initiated backward SRS been demonstrated to dominate the Raman process.337 This occurs when the pump spectral linewidth is much smaller than the Raman linewidth so that the forward and backward Raman gains are comparable. Equivalently, this corresponds to the steady-state Raman regime if the pulse is transform-limited (Fig. 1). The temporal phase profile of the backward-SRS Stokes pulse can be retained in silica fiber, as manifested by successful dechirping of the Stokes pulse from ∼23 ps to 500 fs.333 These observations suggest that backward SRS will have similar coherence properties as forward SRS, which motivates extension of the model presented in this Perspective to the backward situation. To date, most of the investigations of backward SRS have been confined to the steady-state Raman regime. However, the transient regime, as previously shown, can unveil diverse, unexplored, coherent physical phenomena. As demonstrated by Abdolvand et al. in a H2-filled HCF, transient amplification through backward SRS is possible.334 This observation raises questions about the phase relation between the pump and Stokes waves within the transient regime of backward SRS. For example, if either wave is produced by chirping a femtosecond pulse, can a femtosecond Stokes pulse be generated through backward SRS (after dechirped)? An affirmative answer could lead to more-efficient Raman generation without Raman spectral narrowing than the forward situation. Due to the dominance of forward Raman gain, it will be a challenge to demonstrate backward SRS in the transient regime, which might be overcome with gas-filled HCF.
The analysis of the vector properties of SRS in gases sheds insights on the evolution of the polarizations of the interacting fields. On the other hand, preservation of the polarization is essential to various applications, such as interferometric sensors,339–342 frequency metrology,343–345 and quantum communications.346–350 Hollow-core fibers with a symmetric structure can exhibit exceptional polarization purity, with up to 70-dB polarization extinction ratio,351 but only for linear pulse propagation. Nonlinear interactions, prevalent in ultrashort pulse propagation, can introduce significant polarization coupling through XPM, FWM, or SRS as illustrated in this Perspective. Highly birefringent hollow-core fiber preserves the polarization state by introducing a short beat length between polarization modes, effectively suppressing any polarization coupling through either linear or nonlinear interactions. Remarkably, it has achieved birefringence comparable to commercial solid-core polarization-maintaining fibers, accompanied by ultra-low loss.352–356 Further development of such fibers will be desirable for applications that must avoid polarization effects in Raman-active gases.
The vector physics with two polarization modes described here can be considered the initial phase of a broader exploration of multimode gas-based nonlinear optics. Multimode nonlinear effects have spurred numerous studies and potential applications. The trend began in solid-core fibers176,357–365 and has recently been extended to hollow-core fibers. Interested readers are referred to a recent review for multimode research.366 To date, multimode nonlinear effects in hollow-core fibers have been primarily studied in capillaries due to their large core size, where extreme red-shifted spectra,78 multidimensional solitary states,79–82 and few-cycle visible pulse generation367 have been observed. There are also recent developments in multimode hollow-core photonic crystal fibers,355,368–373 as well as the exploration of multimode nonlinear physics, such as circumvention of Raman gain suppression through multimode propagation.42,103 The unified theory covering all temporal regimes can be a starting point for understanding the potential complexities that arise from multimode interactions. As an example, we expect that the Raman-gain equations for co- and cross-polarized multimode fields will resemble those discussed here, but may exhibit unique gain-suppression relations. Intermodal SRS may play out differently in the different temporal regimes. As an example, the wave-vector-matching relations for intermodal nonlocal phonon interactions will be different from the single-mode versions. These differences may open new possibilities for controlling Raman scattering and inspire new applications.
VI. SUMMARY
The analytic theory described in this Perspective relies solely on two fundamental Raman integrals [Eq. (8)]. The theory helps elucidate a wide range of Raman phenomena, including Raman gain suppression in diverse regimes, nonlinear-phase-induced Raman suppression and Raman-pulse compression, and the interplay of Raman-enhanced SPM and its transitional behavior. In addition, interactions of pulses with Raman-generated phonons in various Raman regimes can be investigated. The vector model of Raman interactions presented here facilitates the exploration of nonlinear dynamics influenced by the anisotropic nature of SRS, such as cross-polarized Raman generation, cross-polarized Raman gain, and vector nonlocal interactions. The development of the theoretical framework presented here was motivated by current directions in the interaction of ultrashort pulses with Raman-active gases, and we believe that it will be a valuable tool for future investigations in this area. We hope that this Perspective will be instructive for individuals who are just beginning the investigation of Raman phenomena and will also foster a deeper understanding of Raman physics in experienced researchers, all with the goal of advancing the basic science and applications of Raman scattering.
SUPPLEMENTARY MATERIAL
See the supplementary material for supporting content.
ACKNOWLEDGMENTS
This work was supported by the Office of Naval Research (Grant No. N00014-19-1-2592) and the National Institutes of Health (Grant Nos. R01EB033179 and U01NS128660). The authors acknowledge motivating discussions with Yishai Eisenberg.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Yi-Hao Chen: Conceptualization (lead); Data curation (lead); Formal analysis (lead); Investigation (lead); Software (lead); Validation (lead); Visualization (lead); Writing – original draft (lead). Frank Wise: Funding acquisition (lead); Supervision (lead); Writing – review & editing (lead).
DATA AVAILABILITY
The code used in this work has been made publicly available at https://github.com/AaHaHaa/gas_UPPE. It allows for modeling not only scalar and vector but also single-mode and multimode (transverse modes + polarization modes) situations, as well as in inert (Raman-inactive) and Raman-active gases with either a constant or gradient pressure. In addition, photoionization of a single-mode scalar field is supported. Scripts to generate simulation results in this paper are also included.