Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Design of a passively mode-locking whispering-gallery-mode microlaser

Open Access Open Access

Abstract

Ultrahigh-repetition-rate lasers will become vital light sources for many future technologies; however, their realization is challenging because the cavity size must be minimized. Whispering-gallery-mode (WGM) microresonators are attractive for this purpose since they allow the strong light–matter interaction usually needed to enable mode locking. However, the optimum parameter ranges are entirely unknown since no experiments have yet been conducted. Here, we numerically investigate pulsed operation in a toroidal WGM microresonator with gain and saturable absorption (SA) to study the experimental feasibility. We show that dispersion is the key parameter for achieving passive mode locking in this system. Moreover, the design guideline provided in this work can apply to any small resonators with gain and SA and is not limited to a specific cavity system.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. INTRODUCTION

Mode-locked ultrashort-pulse lasers are vital light sources for laser processing [1], optical communication [2,3], signal processing [4,5], LiDAR application, and remote sensing [6,7]. Glass-based fiber lasers are excellent platforms for these applications because of their high beam quality, robustness, and simple configuration [8]. The two key elements needed to build a mode-locked laser are a gain medium and a mode-locker. Especially erbium (Er)-doped fiber and carbon nanotubes (CNTs) are frequently used for these purposes in a fiber laser.

For the applications mentioned above, the use of a pulsed laser with a high repetition rate will open a new avenue to faster processing, higher capacity, and more effective signal acquisition. However, realizing a repetition rate exceeding 1 GHz remains a challenge [9] since the laser cavity length must be minimized so that the light pulses are densely packed in time. Multi-gigahertz operation has been demonstrated by using a 5 mm Fabry–Perot Er-doped fiber cavity, where the end surface is coated with CNTs as a saturable absorber (SA) [10,11], but further scaling of the repetition rate is not possible because the high pump power needed to compensate for the low gain damages the CNTs.

On the other hand, an ultrahigh-$Q$ whispering-gallery-mode (WGM) microresonator has proven to be an excellent platform that allows light to interact strongly with materials [12]. A WGM microresonator made of ${\rm{Si}}{{\rm{O}}_{{2}}}$ on silicon has exhibited one of the highest $Q$ values of larger than ${\sim}{10^8}$. Since it is made of ${\rm{Si}}{{\rm{O}}_{{2}}}$, it is also a good host for Er doping and will eventually provide a good platform on which to demonstrate a WGM microlaser [13,14]. Continuous-wave (CW) lasing has been demonstrated at an ultralow threshold power thanks to the high $Q$.

Moreover, dissipative soliton generation has also been demonstrated with this type of WGM resonator by applying CW pumping to one of the longitudinal modes. Four-wave mixing in a high-$Q$ WGM microresonator generates an optical frequency comb. If a specific condition is met, we can lock all the generated comb modes and form ultrashort pulses. They are attractive because dissipative solitons from a WGM microresonator have an extremely high repetition rate thanks to their small cavity size. However, such solitons sit on top of a CW background, making it difficult to achieve the pulse amplification that is often required if we are to use this light source as the seed for high-power applications. Moreover, sophisticated wavelength sweeping and feedback are needed to generate and stabilize the soliton pulse [15], and this makes the system complex and expensive.

If we can directly obtain mode-locked pulses from a small WGM microlaser, we would achieve background-free operation without using a complex feedback system. In addition to the CW lasing operation in a silica toroid WGM microresonator, SA is also confirmed with the same device by using CNTs. So the combination of Er-doped ${\rm{Si}}{{\rm{O}}_{{2}}}$ toroid WGM microresonator with CNT as the SA is of interest and has great potential in terms of realizing an ultracompact high-repetition-rate mode-locked pulse laser.

However, it is unclear whether we can demonstrate mode-locked operation in such a device since no experiments have yet been conducted with ${\rm{Si}}{{\rm{O}}_{{2}}}$ WGM microresonators for such a purpose, and the optimum parameter ranges are entirely unknown. Therefore, we aim to reveal the optimum range of parameters, such as the cavity size, $Q$, Er ion concentration, dispersion, and nonlinear loss of the CNTs needed to achieve mode locking, and we discuss the experimental feasibility.

Similar studies have already been published where researchers have tried to achieve an integrated mode-locked laser on a silicon substrate, but we would like to emphasize that we are dealing with different regimes. The technological advances made with microresonators allow us to use ultrahigh $Q$ and sophisticated Er doping in silica on silicon. We expect the high $Q$ to compensate for the limited gain of Er and hence allow us to make the cavity size extremely small.

The paper is organized as follows. Section 2 describes the theoretical model and explains how we determined the parameters that we used in our numerical analysis. Some of the essential parameters taken from our experiments are explained in detail. Section 3 is the core part of this paper, where gain, nonlinear loss, and dispersions are studied as mode-locking parameters. In particular, we reveal the trade-off between the gain and larger dispersion when we increase the cavity size. Additional numerical calculations are performed, which set design guidelines to overcome the trade-off. Finally, we propose a design for passively mode locking a WGM microlaser.

2. THEORETICAL MODEL AND EXPERIMENTAL VERIFICATION OF PARAMETERS

A. Model

Conceptually, a WGM microlaser is a miniaturized fiber ring laser. The system that we studied is shown in Fig. 1. The resonator is made of Er-doped ${\rm{Si}}{{\rm{O}}_{{2}}}$ to enable gain and is fabricated by using the solgel method. CNTs are deposited on the surface of the resonator as a saturable absorber. A CW pump is injected into the cavity, and the output light is coupled out by using a tapered fiber. Throughout this paper, we define the intrinsic and loaded $Q$s as follows. The intrinsic $Q$ (${Q_{{\rm{int}}}}$) is the $Q$ determined by the losses, such as the scattering and absorption of the resonator. The excess absorption losses caused by the CNTs and Er ion doping are normally excluded unless specified because they are usually considered separately in equations. The loaded $Q$ is defined as $Q_{{\rm{load}}}^{- 1} = Q_{{\rm{int}}}^{- 1} + Q_{{\rm{ext}}}^{- 1}$, where ${Q_{{\rm{ext}}}} = {\omega _0}{\kappa ^{- 2}}$ is the $Q$ value determined by the amplitude coupling coefficient between the resonator and the waveguide ($\kappa$) and the angular frequency of the cavity resonance (${\omega _0}$).

 figure: Fig. 1.

Fig. 1. Er-doped silica microtoroid resonator coupled with a tapered fiber that we studied numerically. CNTs are attached to the surface of the resonator to enable saturable absorption.

Download Full Size | PDF

We solve the model based on the nonlinear Schrödinger equation as given by [16,17]:

$$\begin{split}{T_{\rm{r}}}\frac{\partial}{{\partial T}}A({t,T} \,) &= \left({- iL\frac{{{\beta _2}}}{2}\frac{{{\partial ^2}}}{{\partial {t^2}}} + iL\gamma {{| A |}^2}} \right)A({t,T}\, )\\&\quad + \left\{{{g_{{{\rm{T}}_{\rm{r}}}}}(T\, ) - \left[{{l_{{{\rm{T}}_{\rm{r}}}}} + {\alpha _{{{\rm{T}}_{\rm{r}}}}}({t,T}\, )} \right]} \right\}A({t,T} \,),\end{split}$$
where $A$, $t$, and $T$ are the slowly varying field envelope in the microcavity, short time, and long time, respectively. $L$, ${T_{\rm{r}}}$, ${\beta _2}$, $\gamma$, $l$, ${g_{{{\rm{T}}_{\rm{r}}}}}$, and ${\alpha _{{{\rm{T}}_{\rm{r}}}}}$ are the cavity length, one round-trip time, second-order dispersion, nonlinear coefficient, linear loss per round trip (determined by the loaded cavity $Q$), net gain of Er per round trip, and nonlinear loss of the CNTs, respectively. Here, the nonlinear coefficient $\gamma$ $(= {n_2}{\omega _0}{(c{A_{{\rm{eff}}}})^{- 1}})$ was calculated using the nonlinear refractive index of silica ${n_2} = 2.2 \times {10^{- 20}}\; {{\rm{m}}^{{2}}}{\rm{/W}}$.
$${g_{{{\rm{T}}_{\rm{r}}}}}(T\, ) = {g_0}\left({\frac{1}{{1 + \frac{{\overline {\left| {A{{(t,T\,)}^2}} \right|}}}{{P_{{\rm{sat}}}^g}}}}} \right)\left({1 + \frac{1}{{{\omega _{\rm{g}}}^2}}\frac{{{\partial ^2}}}{{\partial {t^2}}}} \right),$$
$${\alpha _{{{\rm{T}}_{\rm{r}}}}}({t,T} ) = {\alpha _{{\rm{ns}}}} + {\alpha _0}\frac{1}{{1 + \frac{{\overline {\left| {A{{(t,T)}^2}} \right|}}}{{P_{{\rm{sat}}}^\alpha}}}},$$
where ${g_{{0}}}$, $P_{{\rm{sat}}}^g$, ${\omega _{\rm{g}}}$, ${\alpha _{{\rm{ns}}}}$, ${\alpha _0}$, and $P_{{\rm{sat}}}^\alpha$ are the saturated gain coefficient per round trip, gain saturation power, gain bandwidth, nonsaturable loss per round trip, modulation depth per round trip, and saturation power, respectively. ($\overline {|A{{(t,T)}^2}|}$) is the field intensity averaged over a short time $t$ (i.e., average power in the cavity). We use the standard split-step Fourier method [18] to perform numerical simulations with a step size equal to the round-trip time ${T_{\rm{r}}}$, for $1.5 \times {10^6}$ round trips.

We must verify our numerical parameters with references and experiments before performing our numerical simulation.

B. Parameters

Here we describe the parameters of the system that we are going to use. The values are summarized in Table 1.

Tables Icon

Table 1. Parameter Values Used in Simulation

Among various parameters, dispersion, gain, and modulation depth are critical and need to be carefully verified since they significantly influence the mode-locking behavior.

1. Cavity Parameters: Dispersion and Q

Dispersion is given when the material and geometry of the cavity are determined [20]. We use an ${\rm{Si}}{{\rm{O}}_{{2}}}$ toroid; hence, we can calculate the dispersion by determining the major diameter $D$ and minor diameter $d$ (Fig. 1). The device is usually fabricated by laser reflow, and $D/d$ is typically kept at about 10. Hence, we assumed $d = D/10$ and calculated the dispersion, as shown in Fig. 2, and use these values for further calculations.

 figure: Fig. 2.

Fig. 2. Second-order dispersion and effective mode area at 1550 nm for an ${\rm{Si}}{{\rm{O}}_{{2}}}$ microtoroid with a different major diameter $D$. The minor diameter $d = D/10$. Inset is the dispersion as a function of wavelength.

Download Full Size | PDF

As we describe in Section 2.B.2, the intrinsic $Q$ is usually on the order of ${10^6}$ to ${10^7}$. Furthermore, a cavity with an intrinsic $Q$ of ${10^8}$ is also possible. Hence, we study in the range from ${10^6}$ to ${10^8}$ for the intrinsic $Q$ in the ${10^6}$ to ${10^8}$ range.

2. Er Ion Parameter: Gain

We fabricated Er-doped ${\rm{Si}}{{\rm{O}}_{{2}}}$ toroids by the solgel method [21]. We obtained a measured $Q$ of ${\sim}{10^7}$ when no ${\rm{E}}{{\rm{r}}^{3 +}}$ was doped. We made two resonators, both with the same diameter of $D \sim 60\; {\rm{\unicode{x00B5}{\rm m}}}$, but with different ${\rm{E}}{{\rm{r}}^{3 +}}$ concentrations of $1.7 \times {10^{18}} \;{\rm{c}}{{\rm{m}}^{- 3}}$ for device A and $0.85 \times {10^{18}} \;{\rm{c}}{{\rm{m}}^{- 3}}$ for device B. The measured loaded $Q$s (i.e., in the presence of the Er absorption effect) of these two cavities were (device A): ${Q_{{\rm{pump}}}} = 1.4 \times {10^6}$ and ${Q_{{\rm{lase}}}} = 1.1 \times {10^6}$, and (device B): ${Q_{{\rm{pump}}}} = 3.0 \times {10^5}$ and ${Q_{{\rm{lase}}}} = 1.1 \times {10^7}$, where ${Q_{{\rm{pump}}}}$ and ${Q_{{\rm{lase}}}}$ are the $Q$s of the pump and lasing modes. We pumped the cavities at ${\sim}1485 \;{\rm{nm}}$ and observed lasing at ${\sim}1605 \;{\rm{nm}}$, where the lasing power is shown as a function of the pump power in Fig. 3(a).

 figure: Fig. 3.

Fig. 3. (a) Experimentally demonstrated laser output-input power curve for cavities with two different ${\rm{E}}{{\rm{r}}^{+ 3}}$ concentrations ${N_{\rm{T}}}$. Pump power and lasing power are values in tapered fiber. We estimate ${\kappa ^2}(= {\omega _0}/{Q_{{\rm{ext}}}})$ to be ${10^7} \;{{\rm{s}}^{- 1}}$ (under-coupling condition) from the transmittance experiment, with some uncertainty because we were not able to measure the transmittance simultaneously when performing the lasing experiment. Inset is a picture of the ${\rm{E}}{{\rm{r}}^{+ 3}}$-doped silica toroid microresonator ($D = 60\; {\rm{\unicode{x00B5}{\rm m}}}$) fabricated by the solgel method and pumped with 1480 nm laser light. (b) $\eta$ versus ${N_{\rm{T}}}$ for different $\kappa$ values ($\kappa = {\kappa _{\rm{p}}} = {\kappa _{\rm{s}}}$) at a cavity corresponding to ${Q_{{\rm{int}}}} = 2 \times {10^7}$ at 1550 nm. We use previously reported values of $\sigma _{\rm{p}}^{\rm{a}} = 1.5 \times {10^{- 21}} \;{\rm{cm}}^2$, $\sigma _{\rm{s}}^{\rm{a}} = 2.8 \times {10^{- 21}} \;{\rm{cm}}^2$, $\sigma _{\rm{p}}^{\rm{e}} = 0.8 \times {10^{- 21}} \;{\rm{cm}}^2$, $\sigma _{\rm{s}}^{\rm{e}} = 4.8 \times {10^{- 21}} \;{\rm{cm}}^2$, and $\tau = 10 \;{\rm{ms}}$ [22]. (c) Gain $g$ (at saturating pump power) as a function of ${N_{\rm{T}}}$. The vertical axis on the right is the saturated gain ${g_0}$ when $D = 300\; {\rm{\unicode{x00B5}{\rm m}}}$. We use ${A_{{\rm{eff}}}} = 21 \;{{\unicode{x00B5}}}{{\rm{m}}^2}$ (Fig. 2) corresponding to a microresonator with $d = 30 \;{\rm{\unicode{x00B5}{\rm m}}}$.

Download Full Size | PDF

When we take account of the coupled-mode theory and upper-state population of Er ions based on rate equations, the slope $\eta$ of the laser output-input power curve is theoretically given as [14]

$$\begin{split} \eta &= \kappa _{\rm{s}}^2\left({\frac{{{\nu _{\rm{s}}}{n_{\rm{s}}}V_{\rm{m}}^{\rm{s}}}}{{{\nu _{\rm{p}}}{n_{\rm{p}}}V_{\rm{m}}^{\rm{p}}}}} \right)\left({\frac{{{\alpha _{\rm{p}}}({{\alpha _{\rm{s}}} + g_{\rm{s}}^ *} ) - ({{\alpha _{\rm{p}}} + g_{\rm{p}}^ *} )\left({{\alpha _{\rm{s}}} + \alpha _{\rm{s}}^{{\rm{passive}}}} \right)}}{{\alpha _{\rm{s}}^{{\rm{passive}}}({{\alpha _{\rm{s}}} + g_{\rm{s}}^ *} )}}} \right)\\&\cdot \left(\!{\frac{{4n_{\rm{p}}^2\kappa _{\rm{p}}^2{{({{\alpha _{\rm{s}}} + g_{\rm{s}}^ *} )}^2}}}{{{c^2}{{\left[{\left({{\alpha _{\rm{p}}} + \alpha _{\rm{p}}^{{\rm{passive}}}} \right)({{\alpha _{\rm{s}}} + g_{\rm{s}}^ *} ) - ({{\alpha _{\rm{p}}} + g_{\rm{p}}^ *} )\left({{\alpha _{\rm{s}}} + \alpha _{\rm{s}}^{{\rm{passive}}}} \right)} \right]}^2}}}} \!\right)\!,\end{split}$$
where $\nu$, $n$, ${V_{\rm{m}}}$, $\kappa$, and $c$ are the light frequency, refractive index, mode volume, amplitude coupling coefficient between the resonator and the waveguide, and vacuum light velocity, respectively. Scripts p and s denote the pump and signal modes, respectively. $\alpha _{{\rm{s,p}}}^{{\rm{passive}}}$ is the loss (in ${{\rm{m}}^{- 1}}$ units) determined by the loaded $Q$ of a passive cavity without ${\rm{E}}{{\rm{r}}^{3 +}}$ doping for signal and pump modes ($\alpha _{\rm{s}}^{{\rm{passive}}} = 2\pi n{({\lambda _{\rm{s}}}{Q_{{\rm{load}}}})^{- 1}})$. The absorption and gain (both in ${{\rm{m}}^{- 1}}$ units) of ${\rm{E}}{{\rm{r}}^{3 +}}$ at a strong pump (i.e., ${N_{\rm{T}}} = {N_2}$ and ${N_1} = 0$) are theoretically given as
$${\alpha _{{\rm{s,p}}}} = {\Gamma _{{\rm{s,p}}}}\sigma _{{\rm{s,p}}}^{\rm{a}}{N_{\rm{T}}},$$
$$g_{{\rm{s,p}}}^ * = {\Gamma _{{\rm{s,p}}}}\sigma _{{\rm{s,p}}}^{\rm{e}}{N_{\rm{T}}},$$
where ${N_{\rm{T}}}$, ${\sigma ^{\rm{a}}}$, and ${\sigma ^{\rm{e}}}$ are the overlap factor between the optical mode with a normalized ${\rm{E}}{{\rm{r}}^{3 +}}$ distribution, ${\rm{E}}{{\rm{r}}^{3 +}}$ concentration, absorption cross section, and emission cross section, respectively. For simplification we assume ${\Gamma _{{\rm{s,p}}}} = 1$ in our calculations.

By substituting Eqs. (5) and (6) into Eq. (4), we obtain graphs providing the relationship between ${N_{\rm{T}}}$ and $\eta$ for different $\kappa$ values as shown in Fig. 3(b). These graphs enable us to estimate the effective ${N_{\rm{T}}}$ from an experimentally obtained $\eta$.

Once ${N_{\rm{T}}}$ is determined, we can obtain ${g_0}$ from Fig. 3(c). The gain can be estimated by solving the rate equation and propagation equation of an erbium-doped waveguide numerically. Hence, Fig. 3(c) is calculated as [22]

$$\frac{{\partial {P_{{\rm{s,p}}}}({z,t} )}}{{\partial z}} = {\Gamma _{{\rm{s,p}}}}\left[{\sigma _{{\rm{s,p}}}^{\rm{e}}{N_2}(t) - \sigma _{{\rm{s,p}}}^{\rm{a}}{N_1}(t)} \right]{P_{{\rm{s,p}}}}({z,t} ),$$
$$\frac{{\partial {N_2}(t)}}{{\partial t}} = - \frac{{{N_2}(t)}}{\tau} - \frac{1}{{{A_{{\rm{eff}}}}}}\left({\frac{{\partial {P_{\rm{s}}}({z,t} )}}{{\partial z}} + \frac{{\partial {P_{\rm{p}}}({z,t} )}}{{\partial z}}} \right){N_2}(t),$$
$${N_1}(t) + {N_2}(t) = {N_{\rm{T}}},$$
where $P$, ${N_1}$, ${N_2}$, and $\tau$ are the light energy (in units of photons), ground-level carrier density, excited-level carrier density, and carrier lifetime, respectively. The gain saturates as we increase the pump power. Figure 3(c) is obtained by plotting the saturated gain ${g_0}$ as a function of ${N_{\rm{T}}}$, where
$${g_0} = {\Gamma _{\rm{s}}}\left[{\sigma _{\rm{s}}^{\rm{e}}{N_2}(t) - \sigma _{\rm{s}}^{\rm{a}}{N_1}(t)} \right]L.$$

${N_{1,2}}({{t_{\rm{s}}}})$ are the carrier densities in a steady state calculated by using Eqs. (7)–(9), when the pump is sufficiently strong.

Figures 3(b) and 3(c) could be used as follows. If ${N_{\rm{T}}}$ is given, we can immediately obtain a theoretical gain $g$ ($= {g_0}/L$) in units of ${{\rm{m}}^{- 1}}$ from Fig. 3(c). On the other hand, Fig. 3(b) allows us to double-check the ${N_{\rm{T}}}$ value if we perform a lasing experiment, as we demonstrated in Fig. 3(a). Figure 3(a) has a slope efficiency of $\eta = 1\%$ at an experimental ${\rm{E}}{{\rm{r}}^{3 +}}$ doping concentration of ${N_{\rm{T}}} = 1.7 \times {10^{18}} \;{\rm{c}}{{\rm{m}}^{- 3}}$. This corresponds to an effective ${N_{\rm{T}}}$ of $5.0 \times {10^{18}} \;{\rm{c}}{{\rm{m}}^{- 3}}$ according to Fig. 3(b) and an experimentally estimated gain of $g = 1 \times {10^{- 2}} \;{\rm{c}}{{\rm{m}}^{- 1}}$ from Fig. 3(c). On the other hand, the theoretical curve in Fig. 3(b) suggests that we should obtain a gain of $g = 3 \times {10^{- 3}} {\rm{c}}{{\rm{m}}^{- 1}}$ at an ${N_{\rm{T}}} = 1.7 \times {10^{18}} {\rm{c}}{{\rm{m}}^{- 3}}$ concentration. Although we do not know exactly why we obtained a slightly higher effective ${N_{\rm{T}}}$, we think the values are in reasonably good agreement and within the experimental error. Taking some experimental uncertainty into consideration, we decided to investigate the ${g_0}= 10^{- 4}$ to ${10^{- 1}}$ range for a $D = 300 \;{\rm{\unicode{x00B5}{\rm m}}}$ cavity.

Equations (7)–(10) also allow us to calculate ${g_0}$ versus ${P_{\rm{s}}}$ at different ${P_{\rm{p}}}$ values, from which we can obtain $P_{{\rm{sat}}}^{\rm{g}}$, which is given as

$$P_{{\rm{sat}}}^{\rm{g}} = \frac{{1 + \frac{{{P_{\rm{p}}}\left({\sigma _{\rm{p}}^{\rm{a}} + \sigma _{\rm{p}}^{\rm{e}}} \right)\tau}}{{{A_{{\rm{eff}}}}}}}}{{\frac{{\left({\sigma _{\rm{s}}^{\rm{a}} + \sigma _{\rm{s}}^{\rm{e}}} \right)\tau}}{{{A_{{\rm{eff}}}}\hbar {\omega _0}}}}}.$$

$P_{{\rm{sat}}}^{\rm{g}}$ is ${\sim}0.145 \;{\rm{W}}$ when we assume a 500 mW pump and $D = 300 \;{\rm{\unicode{x00B5}{\rm m}}}$. We use this value for the numerical calculation.

3. SA Parameters: Saturation Power and Modulation Depth

Saturable absorption is achieved by using CNTs. Here we explain how we verified the SA operation and obtained the SA parameters that we needed for our calculation. The CNTs are dispersed in polydimethylsiloxane (PDMS) at a concentration of 0.24 mg/ml. A PDMS/CNT droplet is transferred from a thin fiber to the ${\rm{Si}}{{\rm{O}}_{{2}}}$ toroid as shown in the inset of Fig. 4. Once the PDMS has been cured by applying heat, the absorption coefficient is measured by a transmittance measurement. The measurement method is described in further detail elsewhere [23]. The result is shown in Fig. 4. The fitting curve gives us a saturation power of $15 \;{\rm{MW}}/{\rm{cm}}^2$ and a modulation depth ${\alpha _0}$ of $3.2 \times {10^{- 3}}$. These values are also consistent with experiments reported by other groups [24]. Therefore, we used $P_{{\rm{sat}}}^\alpha /{A_{{\rm{eff}}}} = 15\; {\rm{MW}}/{\rm{cm}}^2$ and studied a modulation depth range of ${10^{- 4}}$ to ${10^{- 1}}$.

 figure: Fig. 4.

Fig. 4. Absorption coefficient measurement. CNTs are dispersed in PDMS, and a droplet is transferred to the fabricated ${\rm{Si}}{{\rm{O}}_{{2}}}$ toroid as shown in the inset figure.

Download Full Size | PDF

3. NUMERICAL CALCULATION RESULTS AND DISCUSSIONS

This numerical study aims to reveal the optimum parameters needed to achieve mode-locking operation in a tiny Er-doped WGM microresonator. When we monitored the waveform under different conditions, we found four different states as shown in Figs. 5(a)–5(d). These four states are chaotic pulses (CP), multiple pulses (MP), a stable mode-locking regime (ML), and a continuous wave (CW). The regime that we are interested in is the ML regime, where a single pulse circulates in the WGM resonator.

 figure: Fig. 5.

Fig. 5. Operation regimes were identified based on our three criteria: (a) chaotic pulses (CP), (b) multiple pulses (MP), (c) stable mode locking (ML), and (d) continuous wave (CW).

Download Full Size | PDF

We investigated the way in which the gain, loss, and dispersion of the cavity play a role in achieving stable mode locking.

A. Q-Factor Dependence

In this section, the relationship between saturated gain per round trip ${g_0}$ and nonlinear loss (modulation depth ${\alpha _0}$) is investigated for an ${\rm{Si}}{{\rm{O}}_{{2}}}$ toroid microresonator where $D = 300 \;{\rm{\unicode{x00B5}{\rm m}}}$. Figure 6 shows the simulation results for intrinsic $Q= 10^7$ and ${10^8}$. The red line indicates the estimation of ${g_0}$ when ${N_{\rm{T}}} = 5.0 \times {10^{18}}\; {\rm{c}}{{\rm{m}}^{- 3}}$, which we obtained from our lasing experiment (Fig. 3).

 figure: Fig. 6.

Fig. 6. Mode-locking investigation of a $D = 300 \;{\rm{\unicode{x00B5}{\rm m}}}$ toroidal microresonator with gain and nonlinear loss as parameters. (a) Results for $Q= 10^7$ and (b) for $Q= 10^8$. The red line shows the position of the theoretical gain at ${N_{\rm{T}}} = 5.0 \times {10^{18}} \;{\rm{c}}{{\rm{m}}^{- 3}}$.

Download Full Size | PDF

The regimes for the four different states (CP, MP, ML, and CW) are indicated in the figure. When the gain is low, the loss is larger than the gain, and the device will not lase. When the gain is greater than the loss, the device starts to lase. At a higher gain, along with an adequate modulation depth, the system exhibits ML operation. When we increase the gain even further, the pulse splits into multiple components and exhibits an MP state. Finally, the system is in an unstable regime at CP state when the gain is too high.

Figure 6(a) shows that the cavity exhibits CW lasing when $Q$ is ${10^7}$ at ${g_0}= 10^{- 3}$, which explains our CW lasing demonstration in Fig. 3(a). However, it also shows that it is very challenging to achieve the ML state due to insufficient gain if we use this cavity. This situation will change if we use a cavity with a higher $Q$ of ${10^8}$, with which we should be able to reach the ML regime, as shown in Fig. 6(b). Although a $Q$ of ${10^8}$ is experimentally feasible, it is not easy to realize.

Therefore, we will look for another way of achieving ML, such as increasing the WGM microcavity laser diameter or designing the cavity dispersion. By increasing the cavity diameter, we expect to obtain increased gain per round trip. On the other hand, we expect that the mode-locking threshold will decrease if we can make a cavity that has a smaller dispersion. Both approaches will relax the cavity loss conditions for mode locking.

B. Size and Dispersion Dependence

As discussed previously, a sufficiently large ${g_0}$ is crucial for achieving ML. An alternative way of increasing ${g_0}$ is to increase the size $D$ of the microresonator. This approach appears to be straightforward, but the presence of dispersion makes the optimization more complex [20]. Figure 2 shows that the dispersion is at its minimum at $D \sim 150 \;{\rm{\unicode{x00B5}{\rm m}}}$ but becomes anomalous when we increase the cavity size. A smaller dispersion is usually advantageous for easy mode locking because of the stronger self-phase modulation that occurs at a lower intracavity power. The presence of self-phase modulation is normally needed to achieve mode locking. On the other hand, we prefer a larger cavity size because of the larger gain. Therefore, it appears that there is a trade-off between gain and dispersion when we change the cavity size, and there may be an optimum point.

First, we investigated only the dispersion effect. Namely, we assumed that we can fix the diameter of the toroid microresonator at 300 µm but can change the dispersion. Figure 7(a) shows the numerical results we obtained when we fixed the modulation depth ${\alpha _0}$ and $Q$ at $5 \times {10^{- 4}}$ and ${10^7}$, respectively. The red line again shows an experimental estimation of the ${g_0}$. As expected, ML is possible even with a small ${g_0}$ when the dispersion is close to zero. Unfortunately, in this case, however, the dispersion of a 300 µm diameter toroid is about ${-}20 \;{\rm{p}}{{\rm{s}}^2}/{\rm{km}}$, which is outside the ML regime.

 figure: Fig. 7.

Fig. 7. (a) Mode-locking investigation of a 300 µm diameter toroidal microresonator for $Q= 10^7$, with gain per round trip (${g_{{{\rm{T}}_{\rm{r}}}}}$) and dispersion (${\beta _2}$) as parameters. (b) Mode-locking investigation for toroid microcavities ($Q= 10^7$) with different diameters $D$ while $D/d = 10$ is maintained. Note that the vertical axis is now gain per second (${g_{\rm{t}}}$), where ${g_{\rm{t}}} = {g_0}/{T_{\rm{r}}}$. The value is at ${g_{\rm{t}}} = 2.2 \times {10^8}\; {{\rm{s}}^{- 1}}$ when ${g_0}= 10^{- 3}$ and $D = 300 \;{\rm{\unicode{x00B5}{\rm m}}}$ (${T_{\rm{r}}} = 4.5\; {\rm{ps}}$).

Download Full Size | PDF

Second, we performed further numerical investigations at different cavity diameters $D$ for a cavity where $Q$ was fixed at ${10^7}$. The dispersions for cavities with different $D{\rm{s}}$ are taken from Fig. 2. Increasing the diameter will increase the gain per round trip but result in a larger dispersion. On the other hand, the decreasing diameter will reduce the gain per round trip but achieve a weaker anomalous dispersion. Since these two approaches appear to have a trade-off relationship, we confirmed which design approach was more suitable for passively mode locking a WGM microlaser. The result is shown in Fig. 7(b).

Figure 7(b) shows interesting behavior. As expected, the CW lasing threshold decreases as we increase the cavity diameter since the gain increases. However, the ML threshold is not sensitive to the diameter, and it is almost constant when $D \gt 300 \;{\rm{\unicode{x00B5}{\rm m}}}$. This is encouraging because it is telling us that ML is possible even with a small cavity, which is usually more challenging. Taking into consideration that a larger-diameter WGM toroid resonator is often more difficult to fabricate [25], we concluded that the target diameter of the cavity is $D = 300 \;{\rm{\unicode{x00B5}{\rm m}}}$, but with a slightly higher intrinsic $Q$ and higher doping concentration.

Moreover, our results also suggested another approach, namely, dispersion control. If we can shift the null dispersion point at a large $D$, we should be able to reduce the ML threshold, since we can use a larger cavity. Throughout this study, we kept $D/d = 10$, but we might change this ratio (a larger $D$ with a smaller $d$) by adjusting our laser reflow condition.

Finally, we would like to emphasize that this design guideline is not especially for this specific system but could be applied to any microcavity system with gain and SA.

4. CONCLUSION

In this paper, we numerically investigated the passive mode locking of a toroidal WGM microlaser. The small cavity size means that ultrahigh repetition rates can be achieved; on the other hand, the gain per round trip is minimized, and therefore ultrahigh-$Q$ microresonators are necessary. Moreover, dispersion plays an essential role in mode locking: a weak anomalous dispersion promotes ML operation with a limited gain.

In addition, the trade-off relationship between gain and dispersion was highlighted due to their dependence on microresonator diameter. As a result, a 300 µm diameter WGM microlaser with a $Q$ slightly higher than ${10^7}$ proves to be a promising platform for a microlaser with an ultrahigh repetition rate exceeding 100 GHz.

Finally, this design guideline can be applied to any microcavity system that has gain and SA.

Funding

Japan Society for the Promotion of Science (JP18K19036, JP19H00873, JP18H05238); Amada Foundation; Ministry of Education, Culture, Sports, Science and Technology (Q-LEAP).

Disclosures

The authors declare no conflicts of interest.

Data Availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

REFERENCES

1. C. Kerse, H. Kalaycıoğlu, P. Elahi, B. Çetin, D. K. Kesim, Ö. Akçaalan, S. Yavaş, M. D. Aşık, B. Öktem, H. Hoogland, R. Holzwarth, and F. Ö. Ilday, “Ablation-cooled material removal with ultrafast bursts of pulses,” Nature 537, 84–88 (2016). [CrossRef]  

2. L. F. Mollenauer, P. V. Mamyshev, J. Gripp, M. J. Neubelt, N. Mamysheva, L. Grüner-Nielsen, and T. Veng, “Demonstration of massive wavelength-division multiplexing over transoceanic distances by use of dispersion-managed solitons,” Opt. Lett. 25, 704–706 (2000). [CrossRef]  

3. I. Ogura, H. Kurita, T. Sasaki, and H. Yokoyama, “Precise operation-frequency control of monolithic mode-locked laser diodes for high-speed optical communication and all-optical signal processing,” Opt. Quantum Electron. 33, 709–725 (2001). [CrossRef]  

4. D. Cotter, R. J. Manning, K. J. Blow, A. D. Ellis, A. E. Kelly, D. Nesset, I. D. Phillips, A. J. Poustie, and D. C. Rogers, “Nonlinear optics for high-speed digital information processing,” Science 286, 1523–1528 (1999). [CrossRef]  

5. P. T. Callahan, M. L. Dennis, and T. R. Clark Jr., “Photonic analog-to-digital conversion,” Johns Hopkins APL Tech. Digest 30, 280–286 (2012).

6. S. Henderson, P. Suni, C. Hale, S. Hannon, J. Magee, D. Bruns, and E. Yuen, “Coherent laser radar at 2 m using solid-state lasers,” IEEE Trans. Geosci. Remote Sens. 31, 4–15 (1993). [CrossRef]  

7. R. J. D. Young and N. P. Barnes, “Profiling atmospheric water vapor using a fiber laser lidar system,” Appl. Opt. 49, 562–567 (2010). [CrossRef]  

8. A. Martinez and Z. Sun, “Nanotube and graphene saturable absorbers for fibre lasers,” Nat. Photonics 7, 842–845 (2013). [CrossRef]  

9. J. W. Nicholson and D. J. DiGiovanni, “High-repetition-frequency low-noise fiber ring lasers mode-locked with carbon nanotubes,” IEEE Photon. Technol. Lett. 20, 2123–2125 (2008). [CrossRef]  

10. S. Yamashita, Y. Inoue, K. Hsu, T. Kotake, H. Yaguchi, D. Tanaka, M. Jablonski, and S. Set, “5-GHz pulsed fiber Fabry-Pérot laser mode-locked using carbon nanotubes,” IEEE Photon. Technol. Lett. 17, 750–752 (2005). [CrossRef]  

11. A. Martinez and S. Yamashita, “Multi-gigahertz repetition rate passively modelocked fiber lasers using carbon nanotubes,” Opt. Express 19, 6155–6163 (2011). [CrossRef]  

12. T. Tanabe, S. Fujii, and R. Suzuki, “Review on microresonator frequency combs,” Jpn. J. Appl. Phys. 58, SJ0801 (2019). [CrossRef]  

13. L. Yang, T. Carmon, B. Min, S. M. Spillane, and K. J. Vahala, “Erbium-doped and Raman microlasers on a silicon chip fabricated by the sol–gel process,” Appl. Phys. Lett. 86, 091114 (2005). [CrossRef]  

14. B. Min, T. J. Kippenberg, L. Yang, K. J. Vahala, J. Kalkman, and A. Polman, “Erbium-implanted high-Q silica toroidal microcavity laser on a silicon chip,” Phys. Rev. A 70, 033803 (2004). [CrossRef]  

15. T. Herr, V. Brasch, J. D. Jost, C. Y. Wang, N. M. Kondratiev, M. L. Gorodetsky, and T. J. Kippenberg, “Temporal solitons in optical microresonators,” Nat. Photonics 8, 145–152 (2014). [CrossRef]  

16. F. Kurtner, J. der Au, and U. Keller, “Mode-locking with slow and fast saturable absorbers—What’s the difference?” IEEE J. Sel. Top. Quantum Electron. 4, 159–168 (1998). [CrossRef]  

17. H. A. Haus, J. G. Fujimoto, and E. P. Ippen, “Structures for additive pulse mode locking,” J. Opt. Soc. Am. B 8, 2068–2076 (1991). [CrossRef]  

18. G. P. Agrawal, Nonlinear Fiber Optics (Academic, 2007).

19. I. A. Yarutkina, O. V. Shtyrina, A. Skidin, and M. P. Fedoruk, “Theoretical study of energy evolution in ring cavity fiber lasers,” Opt. Commun. 342, 26–29 (2015). [CrossRef]  

20. S. Fujii and T. Tanabe, “Dispersion engineering and measurement of whispering gallery mode microresonator for Kerr frequency comb generation,” Nanophotonics 9, 1087–1104 (2020). [CrossRef]  

21. R. Imamura, S. Fujii, T. S. L. P. Suzuki, R. Suzuki, R. Ishida, M. Ito, H. Maki, L. Yang, and T. Tanabe, “Saturable absorption with CNT coupled WGM and fabrication of Er-doped microresonator for on-chip mode-locked laser,” in Conference on Lasers and Electro-Optics Europe and European Quantum Electronics Conference (Optical Society of America, 2019), paper ck_3_5.

22. P. M. Becker, A. A. Olsson, and J. R. Simpson, Erbium-Doped Fiber Amplifiers: Fundamentals and Technology (Elsevier, 1999).

23. R. Ishida, T. S. L. P. Suzuki, S. Fujii, and T. Tanabe, “The influence of reverse saturable absorption of SWCNTs on mode locking in a fiber-ring resonator,” in Frontiers in Optics + Laser Science APS/DLS (Optical Society of America, 2019), paper JW4A.36.

24. A. Martinez, M. Al Araimi, A. Dmitriev, P. Lutsyk, S. Li, C. Mou, A. Rozhin, M. Sumetsky, and S. Turitsyn, “Low-loss saturable absorbers based on tapered fibers embedded in carbon nanotube/polymer composites,” APL Photon. 2, 126103 (2017). [CrossRef]  

25. J. Ma, X. Jiang, and M. Xiao, “Kerr frequency combs in large-size, ultra-high-Q toroid microcavities with low repetition rates [invited],” Photon. Res. 5, B54–B58 (2017). [CrossRef]  

Data Availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. Er-doped silica microtoroid resonator coupled with a tapered fiber that we studied numerically. CNTs are attached to the surface of the resonator to enable saturable absorption.
Fig. 2.
Fig. 2. Second-order dispersion and effective mode area at 1550 nm for an ${\rm{Si}}{{\rm{O}}_{{2}}}$ microtoroid with a different major diameter $D$ . The minor diameter $d = D/10$ . Inset is the dispersion as a function of wavelength.
Fig. 3.
Fig. 3. (a) Experimentally demonstrated laser output-input power curve for cavities with two different ${\rm{E}}{{\rm{r}}^{+ 3}}$ concentrations ${N_{\rm{T}}}$ . Pump power and lasing power are values in tapered fiber. We estimate ${\kappa ^2}(= {\omega _0}/{Q_{{\rm{ext}}}})$ to be ${10^7} \;{{\rm{s}}^{- 1}}$ (under-coupling condition) from the transmittance experiment, with some uncertainty because we were not able to measure the transmittance simultaneously when performing the lasing experiment. Inset is a picture of the ${\rm{E}}{{\rm{r}}^{+ 3}}$ -doped silica toroid microresonator ( $D = 60\; {\rm{\unicode{x00B5}{\rm m}}}$ ) fabricated by the solgel method and pumped with 1480 nm laser light. (b)  $\eta$ versus ${N_{\rm{T}}}$ for different $\kappa$ values ( $\kappa = {\kappa _{\rm{p}}} = {\kappa _{\rm{s}}}$ ) at a cavity corresponding to ${Q_{{\rm{int}}}} = 2 \times {10^7}$ at 1550 nm. We use previously reported values of $\sigma _{\rm{p}}^{\rm{a}} = 1.5 \times {10^{- 21}} \;{\rm{cm}}^2$ , $\sigma _{\rm{s}}^{\rm{a}} = 2.8 \times {10^{- 21}} \;{\rm{cm}}^2$ , $\sigma _{\rm{p}}^{\rm{e}} = 0.8 \times {10^{- 21}} \;{\rm{cm}}^2$ , $\sigma _{\rm{s}}^{\rm{e}} = 4.8 \times {10^{- 21}} \;{\rm{cm}}^2$ , and $\tau = 10 \;{\rm{ms}}$ [22]. (c) Gain $g$ (at saturating pump power) as a function of ${N_{\rm{T}}}$ . The vertical axis on the right is the saturated gain ${g_0}$ when $D = 300\; {\rm{\unicode{x00B5}{\rm m}}}$ . We use ${A_{{\rm{eff}}}} = 21 \;{{\unicode{x00B5}}}{{\rm{m}}^2}$ (Fig. 2) corresponding to a microresonator with $d = 30 \;{\rm{\unicode{x00B5}{\rm m}}}$ .
Fig. 4.
Fig. 4. Absorption coefficient measurement. CNTs are dispersed in PDMS, and a droplet is transferred to the fabricated ${\rm{Si}}{{\rm{O}}_{{2}}}$ toroid as shown in the inset figure.
Fig. 5.
Fig. 5. Operation regimes were identified based on our three criteria: (a) chaotic pulses (CP), (b) multiple pulses (MP), (c) stable mode locking (ML), and (d) continuous wave (CW).
Fig. 6.
Fig. 6. Mode-locking investigation of a $D = 300 \;{\rm{\unicode{x00B5}{\rm m}}}$ toroidal microresonator with gain and nonlinear loss as parameters. (a) Results for $Q= 10^7$ and (b) for $Q= 10^8$ . The red line shows the position of the theoretical gain at ${N_{\rm{T}}} = 5.0 \times {10^{18}} \;{\rm{c}}{{\rm{m}}^{- 3}}$ .
Fig. 7.
Fig. 7. (a) Mode-locking investigation of a 300 µm diameter toroidal microresonator for $Q= 10^7$ , with gain per round trip ( ${g_{{{\rm{T}}_{\rm{r}}}}}$ ) and dispersion ( ${\beta _2}$ ) as parameters. (b) Mode-locking investigation for toroid microcavities ( $Q= 10^7$ ) with different diameters $D$ while $D/d = 10$ is maintained. Note that the vertical axis is now gain per second ( ${g_{\rm{t}}}$ ), where ${g_{\rm{t}}} = {g_0}/{T_{\rm{r}}}$ . The value is at ${g_{\rm{t}}} = 2.2 \times {10^8}\; {{\rm{s}}^{- 1}}$ when ${g_0}= 10^{- 3}$ and $D = 300 \;{\rm{\unicode{x00B5}{\rm m}}}$ ( ${T_{\rm{r}}} = 4.5\; {\rm{ps}}$ ).

Tables (1)

Tables Icon

Table 1. Parameter Values Used in Simulation

Equations (11)

Equations on this page are rendered with MathJax. Learn more.

T r T A ( t , T ) = ( i L β 2 2 2 t 2 + i L γ | A | 2 ) A ( t , T ) + { g T r ( T ) [ l T r + α T r ( t , T ) ] } A ( t , T ) ,
g T r ( T ) = g 0 ( 1 1 + | A ( t , T ) 2 | ¯ P s a t g ) ( 1 + 1 ω g 2 2 t 2 ) ,
α T r ( t , T ) = α n s + α 0 1 1 + | A ( t , T ) 2 | ¯ P s a t α ,
η = κ s 2 ( ν s n s V m s ν p n p V m p ) ( α p ( α s + g s ) ( α p + g p ) ( α s + α s p a s s i v e ) α s p a s s i v e ( α s + g s ) ) ( 4 n p 2 κ p 2 ( α s + g s ) 2 c 2 [ ( α p + α p p a s s i v e ) ( α s + g s ) ( α p + g p ) ( α s + α s p a s s i v e ) ] 2 ) ,
α s , p = Γ s , p σ s , p a N T ,
g s , p = Γ s , p σ s , p e N T ,
P s , p ( z , t ) z = Γ s , p [ σ s , p e N 2 ( t ) σ s , p a N 1 ( t ) ] P s , p ( z , t ) ,
N 2 ( t ) t = N 2 ( t ) τ 1 A e f f ( P s ( z , t ) z + P p ( z , t ) z ) N 2 ( t ) ,
N 1 ( t ) + N 2 ( t ) = N T ,
g 0 = Γ s [ σ s e N 2 ( t ) σ s a N 1 ( t ) ] L .
P s a t g = 1 + P p ( σ p a + σ p e ) τ A e f f ( σ s a + σ s e ) τ A e f f ω 0 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.