Impact of cavity spectrum on span in microresonator frequency combs

We experimentally study the factors that limit the span in frequency combs derived from the crystalline whispering gallery mode resonators. We observe that cavity dispersion is the key property that governs the parameters of the combs resulting from cascaded four wave mixing process. Two different regimes of comb generation are observed depending on the precise cavity dispersion behavior at the pump wavelength. In addition, the comb generation efficiency is found to be affected by the crossing of modes of different families. The influence of Raman lasing is discussed.

Femtosecond frequency combs have revolutionized precision measurement, optical clocks, communications and spectroscopy. While there have been notable breakthroughs, several years after the demonstrations of the first micro-comb based on miniature whispering gallery mode (WGM) resonators [1][2][3] many questions remain open, preventing realization of its full potential [4]. For example, understanding of the mode locking mechanisms, stability of the RF beat note signals, and relationship between the comb span and cavity parameters is only beginning to emerge [5][6][7][8][9][10][11][12][13][14]. While various aspects of the comb generation have been studied, the comb span limitations have not been explored in detail. A nearly octave comb span and mode locking will be required for optical clocks and optical metrology. Comb spacing also needs to be small enough for electronic detection. It was observed that the comb span generally grows with pumping power, however the coherence of the beatnote is always lost in a broad comb due to generation of the sub-combs [6,15,16]. Moreover, observed comb spans do not always increase with higher power and the causes of this bandwidth limitation are not well understood. In this paper we report results of our experimental investigation of how the pump wavelength and dispersion, mode crossing and Raman lasing influence the comb span and generation efficiency. Our findings suggest a path towards the realization of an octave spanning, coherent microcomb.
To investigate factors limiting the span of microcombs we have fabricated a number of MgF 2 WGM resonators of various size and shape. Each resonator has its axis aligned with the crystalline optical axis (z-cut). We use analytical approximations and finite element method (FEM) [17] to analyse spectrum and dispersion of our resonators. The resonators are used to generate combs for the experimental part of the study.
To study the influence of dispersion and mode crossings on comb span we fabricated a resonator supporting * Corresponding author: grudinin@jpl.nasa.gov only a few families of modes [18]. We used FEM [19] to calculate the free spectral range (FSR) of the first three families of modes as shown in Fig. 1. The resonator profile image was taken through a microscope and digitally processed to extract its shape. The experimentally measured FSR of the combs produced by this resonator is F =46.027±0.003 GHz. During the FEM analysis the cavity radius was changed to fit the fundamental FSR to F . Our FEM solutions used adaptive meshing with over 30k elements, providing absolute frequency precision of a few MHz. It was found that the FSR of the l-m=1 and l-m=2 mode families are 0.02 % and 0.06% larger than the fundamental mode FSR. This means that mode crossings in these basic families of WGMs are possible in the comb pump wavelength region. Optical image processing in Fig. 1 is not accurate enough to determine the frequencies of the non-fundamental modes precisely. The error stems from determining the exact resonator boundary shape. In addition, these families may be shifted in a complex way during the comb generation due to partial overlap with the fundamental modes. These are the reasons why we did not attempt to identify exactly the crossing modes in this case. To study the dispersion the resonator is modelled as an oblate ellipsoid, for which we can compute the spectrum using known approximations [20], iteratively tak-ing the wavelength dependence of refractive index into account. The minor radius of curvature of our resonator r=54 µm is chosen as shown in Fig. 1. Here the major ellipsoid axis is a=750 µm and the minor axis can be approximated as b= √ ar=200 µm. The computed effective group velocity dispersion is shown in Fig. 2. To estimate the precision of this approach we note that changing b from 200 to 150 µm leads to zero dispersion wavelength (ZDW) shift from 1.459 to 1.458 µm. To investigate the role of cavity dispersion in comb dynamics we pumped the fundamental modes of this resonator at 1530 nm, 1560 nm and 1590 nm wavelengths. A fiber laser with 20 kHz linewidth was used at 1560 nm while others are distributed feedback (DFB) semiconductor lasers having larger linewidths. Thermal locking was used to pump the resonator and the stable spectra were recorded. The DFB lasers used for 1530 nm and 1590 nm pumping have linewidth of around 1 MHz. Since laser linewidth exceeded cavity linewidth in these cases, the thermal lock didn't work well and we obtained the comb spectra by repeatedly scanning the lasers around the cavity mode. The comb spectra obtained this way for the case of 1560 nm pump reproduce the envelope of the comb obtained with the thermal lock, although some comb lines are missing.
The comb pumped at 1590 nm started with N=12 FSR primary comb, which evolved into a comb shown in Fig. 3 as the laser was tuned into the resonance. The comb pumped at 1560 nm started at N=14 FSR, and evolved similarly to the 1590 nm comb. The important observation is that a comb changes dynamics notably as the pump wavelength approached the ZDW. The 1530 nm comb did not start along the primary-secondary path. Close to threshold pump power it showed a number of sidebands spaced by 1 FSR with random amplitudes. As the laser pump was increased, a mixture of 1 FSR and 39 FSR peaks was present. Eventually at 9 mW pump power the two strong lines at 79 FSR (80 FSR at 100mW pump) developed as shown in Fig. 3. We note that the increase in N with decreasing cavity dispersion is consistent with the theoretical model for the comb excited in negative GVD regime [6,[11][12][13]. However, the dynamics of a comb pumped near zero cavity GVD have not been extensively analyzed and is expected to be a transition towards the normal GVD combs [2,3,21].
As shown in Fig. 4, there are local variations in FWM efficiency at around 1505 nm and 1585 nm. These spectral features may be explained by the coupling of modes from two different WGM families. The mode crossing induces local change of dispersion and thus phase matching for the FWM process, leading to increase or decrease in comb amplitude. While we have not directly confirmed that the mode crossing is present in our case, several findings point to this explanation of the spectral features. First, the FEM provides the evidence that the modes indeed cross at some wavelength in our resonator. Second, the combs in Fig. 4 A and B are centered around their corresponding pump wavelength, while the spectral features remain at a constant wavelength, indicating the intrinsic property of the cavity under investigation. Third, the spectral features in the comb excited in the fundamental family of WGMs and the features in the comb of the non-fundamental family of WGMs have different wavelength (compare Fig. 4 A-B and C).
We have also carried out experiments aimed at extending the comb span by increasing the pump power. However, we found that the comb span only weakly changes with the pump power. For example, the span of a 45 GHz repetition rate comb in a MgF 2 cavity increased 50% with the 4-fold increase in the pump power from 100 mW to 400 mW at 1560 nm wavelength.
Mode crossing is known to shift the frequencies of the modes locally, thus affecting dispersion [22]. The apparent improvement of comb generation efficiency near the mode crossing feature in Fig.4-B and -A means that the resulting dispersion change strongly influences the FWM process. Along these lines we can attribute the gradual decrease of comb generation efficiency away from the pump wavelength to the change in dispersion as shown in Fig. 2. The mode crossing, while not affecting the comb span directly, provides an important hint: the cavity dispersion (spectrum) needs to be engineered to generate a comb with a broader span. This idea is supported by the significant increase of span when the dispersion is flattened in other nonlinear FWM-based supercontinuum and comb generation examples [15,23].
It is also known that parametric gain is slightly higher in the perfectly phase matched regime than the Raman gain, which does not require phase matching [24]. We found that in larger MgF 2 resonators only parametric generation of comb is observed, while in smaller resonators the Raman lasing notably competes with the generation of combs. To study the concurrent FWM and Raman lasing experimentally, we fabricated a MgF 2 microcavity with FSR comparable to Raman gain linewidth, then changed the cavity temperature to tune the relative offset of the Raman line and mode positions. This cavity is nearly single mode and its dispersion can be modelled with an ellipsoid having the major axis a=190 µm and minor axis b=105 µm. The GVD of such an ellipsoid is shown in Fig. 5.
The FSR of around 180 GHz is derived from the comb line spacing. The Raman gain for the strongest phonon mode in MgF 2 at temperatures 296 and 307.7 K is a lorentzian having the linewidth of 241 and 253 GHz and offset from the pump of 12.148 and 12.139 THz correspondingly [25]. The shift of the cavity modes upon heating was measured to be -50±8 GHz and the change in FSR was negligible. Thus between these two temper- atures the mode closest to the Raman gain peak shifts by around 10 GHz relative to the peak and thus experiences different Raman gain at the two temperatures as shown in Fig. 6. The exact gain is not known due to the error in FSR measurements resulting in the uncertainty (17 GHz) of the frequency of the mode near the Raman peak. The change in gain was enough, however, to observe both Raman lasing and comb generation at room temperature and only comb generation at the elevated temperature. Thus, we find that by changing the cavity temperature we can adjust the relative strengths of these two competing nonlinear processes, benefiting comb generation.
We also note that the larger resonators have anomalous GVD at pump, comb and Raman stokes wavelengths, while the smaller resonators have normal GVD for the pump and the comb wavelengths and near zero GVD for the Raman wavelength (compare Fig. 2 and 5). Due to this behavior of dispersion in small resonators, we were able to observe a mode-locked Raman laser as shown in Fig. 5. This comb is observed in the regime when only the Raman lasing is initially present and is similar to the results of Ref. [26].
Thus, parametric gain always prevails in larger resonators where mode spectrum is dense enough. In smaller resonators the GVD becomes normal due to geometric contribution and the usual FWM leading to combs can not be phase matched. We attribute the comb observed in smaller resonators [27] and the combs observed in near zero GVD pump regime (Fig. 3-top) to a different FWM generation mechanism similar to that observed in CaF 2 resonators also in the normal GVD regime [2,3].
In conclusion, we have investigated several mechanisms which play an important role in the span of microcombs generated in crystalline cavities.
Our findings suggest that total cavity dispersion -as well as its wavelength dependence -determine the comb span. We observe that while mode crossings may not directly limit comb width, the resulting local changes in dispersion may increase or decrease FWM efficiency in that region. We show that by tuning cavity temperature, we may suppress the competing process of stimulated Raman scattering. Our observations suggest a route towards the broad coherent comb in crystalline WGM resonators: the spectrum must be engineered to have a flat anomalous GVD in the desired spectral region and a few if any mode crossings. The flatness of the GVD is equivalent to low third order dispersion.
This work was carried out at the Jet Propulsion Laboratory, California Institute of Technology under a contract with the National Aeronautics and Space Administration.