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

Numerical analysis of synchronously pumped solid-state Raman lasers

Open Access Open Access

Abstract

Considering the spatial distribution of laser beams and phonon waves, the SRS coupling wave equations in the transient regime are derived and normalized for the first time. The synchronously pumped solid-state Raman laser is simulated numerically to investigate the influences of the cavity length detuning, output coupling rate, dispersion, Raman gain and dephasing time of Raman mode on laser performances. It is found that the intensive pulse compression of first Stokes laser in synchronously pumped solid-state Raman laser stems from pulse width gain narrowing and intensity oscillation effects. The cavity length detuning, dispersion, Raman gain and dephasing time considerably affect the pulse width gain narrowing and intensity oscillation processes. The theoretical results can help the design and optimization of synchronously pumped solid-state Raman laser to generate ultrafast Raman laser output efficiently.

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

1. Introduction

Picosecond pulses have been widely used in industrial, medical, military and information fields [14]. Solid-state Raman lasers in different configurations can efficiently generate IR, visible and UV output [510]. For the transient regime where the pump laser pulse duration is similar or shorter than the Raman transition dephasing time T2 (typically T2 is 1∼10 ps for crystalline Raman materials), stimulated Raman scattering (SRS) process has smaller Raman gain and higher threshold compared to the steady state regime [10]. To overcome those withdraws, synchronously pumped Raman lasers have been demonstrated to realize efficient Raman conversion with picosecond pumping laser [1113]. For synchronously pumped Raman lasers, the resonator length of external Raman laser matches the repetition rate of the pump laser, and the pumping pulse and intracavity Stokes pulse propagate in step inside the resonator.

Synchronously pumped solid-state Raman laser is a kind of frequency convertor with good features, such as simple structure, high conversion efficiency, good stability, no need for phase matching, etc. The most attracting feature is that pulse width of Stokes lasers can be greatly compressed resulting in much enhanced peak intensity compared to pumping laser pulses. Recent years, many experimental works on synchronously pumped solid-state Raman laser in picosecond regime have been reported in IR, visible and UV spectrum span [1318]. The research group from Macquarie University used 532nm 28ps mode-locked laser to synchronously pump KGW crystal. It was found that the pulse width of first Stokes laser can be greatly compressed by cavity length detuning (CLD) [13,17], and the shortest pulse width of first Stokes laser obtained was 6.5ps [13]. The same group adopted the 532nm 26 ps pulsed laser as pumping source to simultaneously pump 6.7mm long diamond crystal, and the pulse width was reduced to 9ps with CLD of 200μm. M. Frank et al. of the Czech Republic reported GdVO4, BaWO4 and SrMoO4 as Raman crystals to conduct synchronously pumped Raman laser experiments, and obtained the extremely compressed pulses of Stokes laser by shifting the pumping laser with two different Raman modes in sequence [1922]. In 2019, Frank M. synchronously pumped a-cut SrMoO4 external Raman laser with ring type cavity, 1173nm first Stokes laser generated by the Raman mode of 887cm−1 was successively shifted by the Raman mode of 327cm−1 to generate 1220nm Raman laser. By cavity length detuning, the minimum pulse duration of 1220nm was 1.4ps, which was compressed by 26 times compared to the pumping pulse duration of 36ps [19].

Nonlinear coupling equations of transient SRS were adopted to investigate Raman processes with gas Raman media [11,12,23,24]. In 2010, Granados E. et al. numerically simulated synchronously pumped solid-state Raman lasers with the transient coupled wave equations [14,25], and it was found that the pulse compression effect was produced by a combination of group velocity walk-off and strong pump pulse depletion [14,25]. However, full understanding of synchronously pumped solid-state Raman laser is still lacked, and there is no research to investigate the influences of main parameters on the laser performance, i.e. output coupling rate, dispersion, Raman gain, dephasing time of Raman mode and etc. Moreover, as a third nonlinear process, SRS conversion is sensitive to the beams intensities, and plane wave theories bring about considerable discrepancy [11,12,24,25]. In [26], Boyd G. et al. investigated the threshold condition for steady-state stimulated Raman oscillation by considering the paraxial Gaussian beam distributions, and got results of better precision. Our group deduced the space-dependent rate equations for intracavity Q-switched Raman lasers in the steady-state regime, and also found the new equations have better approximation than the plane-wave ones [9]. In this paper, assuming the laser beam transverse distributions to be of the Gaussian profile, considering laser beam longitudinal spacial distribution, the space-dependent transient SRS coupling wave equations are derived and normalized for the first time. By numerically solving the normalized coupling wave equations, synchronously pumped Raman laser is simulated for the full insight of its operation. The influences of the cavity length detuning, output coupling rate of the first Stokes laser, dispersion, dephasing time of Raman mode and Raman gain on Raman laser performances are investigated systematically. The new understandings of the pulse compression mechanism are revealed. The research results can help the design and optimization of synchronously pumped solid-state Raman laser to generate ultrafast Raman laser output efficiently.

2. Theory

The optical fields are described classically in terms of Maxwell’s equations, and the molecular system is treated quantum mechanically. The SRS coupling equations in the transient regime are deduced quantum-mechanically in the Gaussian unit system [27].

In the Raman medium, the upper state population changes induced by the resonant interaction can be neglected. ωv,L,S and kv,L,S respectively denote the angular frequency and wave number of the material excitation, pumping laser and first Stokes laser, and ω0 and k0 is the frequency and wave number of the material excitation at resonance. With relations of ωv =ωL -ωS and kv = kL - kS, assuming pulse durations and dephasing time T2 of phonon is in the picosecond range or longer, the transient Raman coupling equations at the resonant condition (ωv=ω0) are deduced. For the new coordinates (t, z), which are in a frame moving with the velocity of first Stokes laser vS, the following equations are obtained in [25].

$$\frac{{\partial Q}}{{\partial t}} - {v_S}\frac{{\partial Q}}{{\partial x}} + \frac{1}{{{T_2}}}Q = \frac{{{\kappa _Q}}}{{{T_2}}}{E_L}E_S^\ast{+} \hat{F}({x,t} ),$$
$$\frac{-i}{2 k_{S}} \nabla_{\perp}^{2} E_{S}+\frac{1}{v_{S}} \frac{\partial E_{S}}{\partial t}+\frac{\gamma_{S}}{2} E_{S}=g_{S} \kappa_{S} Q^{*} E_{L},$$
$$\frac{-i}{2 k_{L}} \nabla_{\perp}^{2} E_{L}+\left(1-\frac{v_{S}}{v_{L}}\right) \frac{\partial E_{L}}{\partial x}+\frac{1}{v_{L}} \frac{\partial E_{L}}{\partial t}+\frac{\gamma_{L}}{2} E_{L}=g_{S} \kappa_{L} Q E_{S}.$$
with
$${\kappa _Q} = i\frac{{c{\mu _L}{\mu _S}}}{{16\pi {\omega _S}}},\; \; \; {\kappa _S} = i\frac{{{\omega _S}}}{{{\mu _S}}},\; \; \; {\kappa _L} = i\frac{{{\omega _L}}}{{{\mu _L}}}.$$
where transvers Laplace operator ▽2 is included to manipulate the transverse distribution of laser beams and phonon excitation. EL,S is electric field amplitude of pumping and first Stokes lasers, respectively, and Q is the amplitude of phonon excitation. c is the light speed in the vacuum.  μL,S is the refractive index of pumping and first Stokes lasers, and γL,S is the absorption coefficient of pumping and first Stokes lasers, respectively. κQ, κS and κL are coupling coefficients. gS is the steady-state Raman gain, and vL is the velocity of pumping laser.

The quantum statistical Langevin operator $\hat{F}({z,t} )$ is added in (1) corresponding to random fluctuations, which accounts for the spontaneous Raman scattering [28]. The Langevin force is taken to be delta correlated, and the relative deduction obeys

$$\langle{\hat{F}^\dagger }({x,t} )\hat{F}({x^{\prime},t^{\prime}} )\rangle= \frac{{{g_S}\hbar {\mu _L}}}{{2{A_S}{\omega _S}T_2^2}}\; \delta ({x - x^{\prime}} )\delta ({t - t^{\prime}} ).$$
$$\langle\hat{F}({x,t} )\rangle= 0.$$
$$\langle\hat{F}({x^{\prime},t^{\prime}} )\hat{F}({x,t} )\rangle= 0.$$
where AS is the cross section area of first Stokes beam, and $\hbar $ is the Plank constant.

The SRS is third nonlinear process, and sensitive to the field intensity. To precisely model the SRS process in the transient regime, the longitudinal and transverse spatial distributions of laser fields should be taken into consideration. Fields of fundamental and first Stokes beams are assumed to be of Gaussian transverse distribution for the TEM00 mode, and the electric field amplitudes of pumping and first Stokes lasers inside the cavity are expressed as:

$${E_{L,S{\; }}}({r,x,t} )= \frac{{E{0_{L,S{\; }}}({0,t} )}}{{\sqrt {1 + \frac{{{x^2}}}{{ZR_{L,S}^2}}} }}exp \left[ {\frac{{ - {r^2}}}{{w_{L,S}^2(x )}}} \right],$$
where r is radial coordinate, and E0L,S(0, t) is the electric field amplitudes on the axis at the waist positions of pumping and first Stokes lasers, respectively. The waist positions of pumping and first Stokes beams are assumed to be at x=0. wL,S is radii of the pumping and first Stokes beams inside the resonator, respectively. The radii vary along the axis according to
$${w_{L,S}}(x )= w{0_{L,S}}\sqrt {1 + {x^2}/ZR_{L,S}^2} ,$$
$$Z{R_{L,S}} = \pi w0_{L,S}^2/{\lambda _{L,S}},$$
where λL,S, w0L,S and ZRL,S is the wavelength, waist radius and Rayley distance of pumping and first Stokes lasers, respectively.

The electric field amplitude of input pumping laser is expressed as:

$${E_{L\; in}}({r,x,t} )= \frac{{E{0_{L0}}({0,t} )}}{{\sqrt {1 + \frac{{{x^2}}}{{ZR_L^2}}} }}\exp \left[ {\frac{{ - {r^2}}}{{w_L^2(x )}}} \right],$$
$$E{0_{L0}}({0,t} )= E{p_{max}}\exp \left[ {\frac{{ - 2ln2{{({t - {t_0}} )}^2}}}{{t_p^2}}} \right],$$
where E0L0 (0, t) is the electric field amplitudes of input pumping laser on the axis at the waist position, and is assumed to have Gaussian temporal distribution as given in (10). Epmax is the peak field amplitude of the input pumping pulse, tp is the pumping pulse width, and t0 is the pulse peak position of input pumping laser on the temporal axis t.

The phonon excitation is driven by multiplication of ES* and EL according to Eq. (1). When the wavelength of phonon wave is much less than beam radius of pumping and first Stokes lasers, the phonon excitation can be reasonably assumed to be

$$Q({r,x,t} )= Q({0,x,t} )exp \left( {\frac{{ - {r^2}}}{{w_L^2}} + \frac{{ - {r^2}}}{{w_S^2}}} \right),$$
where Q(0, x, t) is the amplitude of phonon excitation of Raman medium on the laser axis.

The transvers Laplace operator in cylindrical coordinate system is expressed as:

$$\nabla _ \bot ^2E = \frac{1}{r}\frac{\partial }{{\partial r}}\left( {r\frac{{\partial E}}{{\partial r}}} \right) + \frac{1}{{{r^2}}}\frac{{{\partial ^2}E}}{{\partial {\theta ^2}}}.$$
Substituting Eqs. (6)-(9) and (11) into Eqs. (1)-(3), integrating the differential equations over the cross section to guarantee the beam spatial distribution during the SRS process, the space-dependent coupling equations can be obtained.

For the generality and simplicity, we introduce the normalized time τ, normalized coordinate ζ, normalized electric field amplitudes EL and ES, normalized phonon amplitude Θ, normalized Raman gain coefficient G, normalized light velocities VL and VS, normalized dephasing time Τ of Raman mode as:

$$\tau = t/{t_p},\; \; \zeta = x/Z{R_L},\; \; {{{\rm E}}_L} = E{0_L}({0,t} )/E{p_{max}},\; \; {{{\rm E}}_S} = E{0_S}({0,t} )/E{p_{max}},$$
$$\mathrm{\Theta } = {Q^{\prime}}({0,x,t} )/Ep_{max}^2,\; \; \; G = {g_S} \cdot Ep_{max}^2 \cdot Z{R_L},$$
$$\; {\; \; }{V_L} = {v_L} \cdot {t_p}/Z{R_L},\; \; \; \; {V_S} = {v_S} \cdot {t_p}/Z{R_L},{{\rm T}} = {T_2}/{t_p},$$
ZRL is assumed to be equal to ZRS, and w0L equals to w0S. With normalized parameters in (13), the normalized equations can be obtained as:
$$\frac{{\partial \mathrm{\Theta }}}{{\partial \tau }} - {V_S}\frac{{\partial \mathrm{\Theta }}}{{\partial \zeta }} - {V_S}\frac{{2\zeta }}{{({1 + {\zeta^2}} )}}\mathrm{\Theta } + \frac{1}{{{\rm T}}}\mathrm{\Theta } = \frac{{{\kappa _Q}}}{{{\rm T}}}\frac{1}{{({1 + {\zeta^2}} )}}{{{\rm E}}_L}{{\rm E}}_S^\ast{+} \hat{\textrm{F}}({\zeta ,\tau } ),$$
$$\frac{{\partial {{{\rm E}}_S}}}{{\partial \tau }} + {\mathrm{\Gamma }_S}{{{\rm E}}_S} = \frac{1}{3}{\kappa _S}G{V_S}{\mathrm{\Theta }^\ast }{{{\rm E}}_L},$$
$$({{V_L} - {V_S}} )\frac{{\partial {{{\rm E}}_L}}}{{\partial \zeta }} + ({{V_L} - {V_S}} )\frac{\zeta }{{({1 + {\zeta^2}} )}}{{{\rm E}}_L} + \frac{{\partial {{{\rm E}}_L}}}{{\partial \tau }} + {\mathrm{\Gamma }_L}{{{\rm E}}_L} = \frac{1}{3}{\kappa _L}G{V_L}\mathrm{\Theta }{{{\rm E}}_S},$$
where the normalized absorption coefficient ΓL, S and normalized Langevin operator is defined as:
$${\mathrm{\Gamma }_L} = \frac{{{\gamma _L}}}{2} \cdot {V_L} \cdot Z{R_L},\; \; {\mathrm{\Gamma }_S} = \frac{{{\gamma _S}}}{2} \cdot {V_S} \cdot Z{R_L},\; \; \hat{\textrm{F}}({\zeta ,\tau } )= \hat{F}({x,t} )\cdot {t_p}/Ep_{max}^2.$$
In the normalized space-dependent coupling equations, there are several composite variables, i.e. G, VL, VS, Τ, ΓS and ΓL. By solving the normalized coupling equations numerically, we can investigate the influences of the composite variables on the performance of synchronously pumped Raman lasers.

In the following simulation, the normalized intensity of first Stokes pulse IS is defined as:

$${I_S} = {{{\rm E}}_S} \cdot {{\rm E}}_S^\ast .$$
ISmax is the normalized peak intensity of first Stokes pulse. The normalized intensity integration of first Stokes pulse ISinteg is defined as:
$${I_{Sinteg}} = \smallint ({{{{\rm E}}_S} \cdot {{\rm E}}_S^\ast } )d\tau ,$$
which is proportional to the pulse energy and conversion efficiency of first Stokes laser. The normalized pulse duration of the first Stokes laser τS=τSr+τSf is obtained from (εS·εS*), here τSr is the normalized rising time from ISmax/2 to ISmax, and τSf is the falling time from ISmax to ISmax/2.

3. Setup for simulation

Simulations are carried out for the synchronously pumped z-fold external ring SRS laser as shown in Fig. 1. The cavity consists of four mirrors, which are pumping mirror (PM), M1, M2 and flat output coupler (OC). PM is coated HT for pumping laser and HR for first Stokes laser, M1 and M2 have coating HR for first Stokes laser, and OC has partial reflectivity for first Stokes laser. All mirrors are coated HT for the second Stokes laser to prohibit its generation through cascaded SRS. The pumping source adopted is usually CW mode-locked lasers, for example, a 532nm laser with an average power of 2W, pulse duration of 10ns and repetition rate of 80 MHz in [25]. The Raman crystal is assumed to be KGW, and the wavelength of the pumping and first Stokes lasers are 532nm and 559nm, respectively. For the cavity length detuning (CLD), Δx=0 corresponds to the perfect synchronization. For positive Δx, the cavity length of Raman laser becomes longer than that of perfect synchronization, and the first Stokes pulse lags behind the pumping pulse. For negative Δx, the cavity length becomes shorter, and the first Stokes pulse goes before the pumping pulse. CLD is realized by adjusting M2.

 figure: Fig. 1.

Fig. 1. Cavity setup of a synchronously pumped solid-state Raman laser

Download Full Size | PDF

The normalized space-dependent coupling Eqs. (14)-(16) are solved numerically to model the Raman interaction between the incident infinite train of pumping pulses and first Stokes intracavity pulse. The first Stokes laser is generated through the spontaneous Raman scattering brought about by quantum fluctuations of the Raman medium at the initial stage, then is amplified through SRS process at the following stage. The interactions of circulating first Stokes pulse with a sequence of pumping pulses are modeled iteratively, and the simulation is terminated when the Stokes pulse has reached its steady state profile after about 300 pumping pulses. The dispersive broadening of the Stokes pulse is approximately modelled by applying a discrete broadening after each pass through the Raman crystal, and is sufficiently accurate in the picosecond regime [25]. The values of normalized variables are estimated with the experiments parameters as given in Table 1.

Tables Icon

Table 1. Estimated values of the normalized parametersa

4. Results and discussion

4.1 Effects of normalized CLD Δζ on output characteristics

The effects of normalized CLD Δζ on output characteristics of synchronously pumped solid-state Raman laser are investigated numerically as shown in Fig. 2. The normalized intensity integration ISinteg, normalized maximum intensity ISmax, and normalized pulse width τS of first Stokes pulse are shown in Fig. 2(a) and (b), with RS=80% or 90%, G=4.0×10−10 or 8.0×10−10, and other parameters given in Table 1.

 figure: Fig. 2.

Fig. 2. Effects of normalized CLD Δζ on output characteristics of first Stokes laser with RS=80% or 90%, G=4.0×10−10 or 8.0×10−10, and other parameters given in Table 1.

Download Full Size | PDF

CLD has considerable effects on the pulse energy, pulse width and peak power of first Stokes laser outputs for the synchronously pumped Raman laser [13,17]. As shown in Fig. 2(a), normalized CLD Δζ ranges from -0.04∼0.006, and there is a maximum ISinteg for each line, which corresponds the maximum SRS conversion efficiency. From the maximum value, ISinteg declines rapidly as Δζ increases to the positive direction, and decreases gradually as Δζ varies to the negative regime. For G=4.0×10−10, ISinteg of RS=80% is greater than that of RS=90% for negative Δζ, however, for the positive Δζ, ISinteg of RS=80% declines more rapidly, and is smaller than the result of RS=90%. ISinteg of G=8.0×10−10 is much higher than that of G=4.0×10−10. As shown in Fig. 2(b), lines of normalized pulse duration τS have the similar trend. τS declines gradually as Δζ varying from -0.04 to -0.003, and has an abruptly decrease within Δζ =-0.003∼0.0009 to reach the minimum value, where ISmax rises sharply to reach the maximum value. Compared to G=4.0×10−10, the maximum value of ISmax is much higher for G=8.0×10−10. As Δζ further increase from 0.0009 to 0.006, τS increases obviously, and ISmax declines sharply. Lines of ISmax all have an abrupt peak corresponding to the minimum value of τS.

There are two interesting working points adjacent to Δζ = 0 as shown in Fig. 2(a). Under working point 1, which usually has a small negative value of Δζ, the maximum Raman conversion efficiency can be achieved, when the pulse width and peak power of first Stokes laser are close to those of pumping pulses. Under working point 2, which commonly has a small positive value of Δζ, the first Stokes laser pulse is compressed considerably, where it has the maximum peak power, minimum pulse width and moderate conversion efficiency. The relative parameters of working point 1 and 2 are listed in Table 2. The experimental parameters, such as output coupling rate, Raman gain and etc., have little impact on the value of Δζ for working point 2, which ranges from 0.0005 to 0.001. As the absolute value of Δζ increases, the performance of synchronously pumped Raman lasers deteriorates. However, synchronously pumped Raman lasers have higher tolerance for positive Δζ than for negative one.

4.2 Effects of normalized CLD Δζ on pulse temporal distributions

In Fig. 3, normalized intensity temporal distributions of first Stokes laser output and residual pumping laser are given for different Δζ with G=8.0×10−10, RS=80% and other parameters given in Table 1. Δζ is equal to 0.006, 0.0009, 0, -0.008, -0.02 and -0.004 in (a), (b), (c), (d), (e) and (f), respectively. In Fig. 3, pulses profiles are in a frame moving with the velocity of first Stokes laser VS to the left, and left edges of pulses are the leading ones as shown in Fig. 3(a).

 figure: Fig. 3.

Fig. 3. Normalized intensity temporal distributions of first Stokes laser output and residual pumping laser for different Δζ with G=8.0×10−10, RS=80% and other parameters given in Table 1. Δζ is equal to 0.006, 0.0009, 0, -0.008, -0.02 and -0.04 in (a), (b), (c), (d), (e) and (f), respectively.

Download Full Size | PDF

It is found that pulse compression of first Stokes laser can be realized for working point 2 as depicted in Fig. 2(a) with appropriate CLD. In the transient SRS, the gain at a point on the first Stokes envelope depends on the integral of the pumping power up to that point, which leads to the delayed appearance of the transient phonon wave and first Stokes laser pulse with respect to the pumping pulse. SRS process is a third nonlinear process, and has threshold effect. The Stokes pulse usually has a leading edge which rises sharply to a maximum with some delay with respect to the maximum of the pumping pulse. For positive CLDs as given in Fig. 3(a) and (b), the first Stokes pulse temporally lags behind the pumping laser pulse by an additional detuning delay, and the first Stokes pulse will gradually catch up with the pumping pulse due to dispersion. Thus, the leading edge of the first Stokes pulse sweeps through the pumping pulse temporally in the Raman medium, and always meets the undepleted part of pumping pulse to obtain much higher Raman gain than the falling edge. The unbalanced gain leads to pulse width gain narrowing effect, which results in intensive pulse compression, and forms giant Stokes pulse of high peak intensity and small pulse width. As shown in Fig. 3(b), the first Stokes pulse has the smallest pulse width and maximum peak intensity with Δζ = 0.0009, which is unexpectedly smaller because the Stokes pulse usually lag behind the pumping pulse as mentioned above. When CLD further increases in the positive direction, the peak intensity and conversion efficiency of first Stokes laser drop dramatically, because the first Stokes pulse misses the pumping pulse peak, and large portion of pumping energy cannot be converted as shown in Fig. 3(a). Thus, the pulse compression of first Stokes laser comes from the pulse width gain narrowing, which utilizes the group velocity difference of pumping and first Stokes laser due to linear dispersion under appropriate CLD.

As shown in Figs. 3(a), (b) and (c), the pulse compression is usually accompanied by intensity modulation. As intense ultrafast laser pulses go through Raman-active media, the first Stokes and phonon fields grow so large that the pumping field is driven toward zero. Since the phonon responds slowly, the transfer of power from the pumping laser to first Stokes reverses, resulting in a flow of energy back into the pumping [23]. Thus, the depleted pumping pulse and first Stokes pulse show a series of decreasing oscillations, and more oscillations follow if the depletion level is larger as shown in Fig. 3(b). The intensity oscillation can further narrow the pulse width of first Stokes laser, however, clamp its peak intensity.

As CLD goes into negative regime with relatively small value, the pumping and first Stokes laser pulses have a good temporal overlap in the Raman crystal, resulting in high conversion efficiency as shown in Figs. 3(d) and (e). It is noted that the best temporal overlap does not occur at the Δζ = 0 due to the velocity dispersion and the behind-lagging of first Stokes pulses, and the maximum ISinteg is obtained with Δζ =-0.008 as shown in Fig. 3(d). As the absolute value of CLD further increases in the negative regime, the temporal shape of the first Stokes pulse is reshaped by the high Raman gain for its trailing edge. With negative CLD, the first Stokes pulse has a large pulse width with a long leading edge and relatively short trailing edge as shown in Fig. 3(f), and the peak intensity decreases sharply.

4.3 Effects of dispersion on pulse compression

The influence of dispersion on the transient nonlinear process of synchronously pumped solid-state Raman laser was investigated in [25], and need further study. Based on the numerical simulation, it is found that the linear dispersion has important effects on pulse compression. Figure 4 shows the effects of dispersion on pulse compression of the first Stokes laser with RS = 90%, G=1.0×10−9, Τ=0.5 and other parameters given in Table 1. Solid lines are the results with dispersion, and dashed lines are those without dispersion. With dispersion, the pulse width of first Stokes laser is compressed to the minimum value of 0.0052 for Δζ = 0.0003 reaching the maximum peak intensity of 5.51 as depicted in Fig. 4(a). Without dispersion, there is no pulse compression effect. The line of ISmax is flat, and τS decrease gradually as Δζ increasing in Fig. 4(b).

 figure: Fig. 4.

Fig. 4. Effects of dispersion on pulse compression of the first Stokes laser with RS = 90%, G=1.0×10−9, Τ=0.5 and other parameters given in Table 1. Solid lines are the results with dispersion, and dashed lines are those without dispersion.

Download Full Size | PDF

To illustrate the effect of dispersion on pulse compression more clearly, Fig. 5 depicts the effects of dispersion (VS-VL) on pulse compression with RS = 90%, Δζ = 0.0005 and other parameters given in Table 1. VS-VL is equal to 0, 0.0013, 0.0026 and 0.0052 in (a), (b), (c) and (d), respectively. It is found that SRS amplification is favorably affected by a velocity mismatch between the pumping and first Stokes pulses due to linear dispersion. As shown in Fig. 5, as the dispersion is increasing, the greater displacement for the pumping laser pulse occurs relative to the current reference frame, and the first Stokes pulse goes through the pumping pulse profile to extract more energy from the pumping pulse. Consequently, as the dispersion increases, the intensity of the residing pumping pulse decreases, and the conversion efficiency of first Stokes laser is improved.

 figure: Fig. 5.

Fig. 5. Effects of dispersion (VS-VL) on pulse compression of the first Stokes laser with RS=90%, Δζ = 0.0005 and other parameters given in Table 1. VS-VL is equal to 0, 0.0013, 0.0026 and 0.0052 in (a), (b), (c) and (d), respectively.

Download Full Size | PDF

The group velocity dispersion is the critical factor for the pulse compression [23,25]. Having velocity dispersion, the Stokes pulse can propagate forward through the pump pulse, and the leading edge of the Stokes field sweep through the pumping pulse profile which have not been depleted. The Stokes pulse extracts more energy from the pumping pulse, and develops into a giant pulse with a greatly compressed pulse width and a high peak power. Without dispersion, there is no obvious pulse compression of the first Stokes laser pulse (τS=0.808) as shown in Fig. 5(a). With dispersion, the pulse width compression of the first Stokes pulse occurs as shown in Fig. 5(b), (c) and (d). As dispersion (VS-VL) increases from 0.0013 to 0.0026, τS decreases from 0.296 to 0.068 as depicted in Figs. 5(b) and (c). On the other hand, the dispersion brings about broadening of pulse width, and limits the narrowest pulse width of first Stokes pulses. As dispersion (VS-VL) increases to 0.0052, τS increases to 0.096 as shown in Fig. 5(d). In a nondispersive medium, the Stokes pulse profile has no intensity oscillation as shown in Fig. 5(a). As the dispersion increases, intensity oscillation of the pumping and first Stokes laser pulses get more intensive as shown in Figs. 5(b), (c) and (d).

With dispersion, the pumping pulse moves to right relative to the current reference system to generate the time delay, and the first Stokes laser pulse actually moves to right along with the pumping pulse through the SRS amplification. The higher the dispersion is, the larger the time delay is. Passing through the Raman crystal with parameters given in Table 1, the time delay of the pumping pulse is calculated to be 0.23 relative to the current reference system for VS-VL = 0.0013 (half dispersion). For VS-VL = 0.0052 (double dispersion), the time delay is estimated to be 0.92. The first Stokes pulse with double dispersion in Fig. 5(d) lags behind that of half dispersion in Fig. 5(a) obviously. Therefore, the first Stokes pulse and the pumping pulse can propagate basically synchronously under different dispersion conditions. The experimental conditions, such as output coupling rate, Raman gain, dephasing time and dispersion, have little effect on Δζ of the working point 1 and 2 as given in Table 2. It is also noted that there are different definitions of CLD. In this paper, Δζ = 0 corresponds to the perfect synchronization, for which the round trip time of first Stokes laser is equal to the time interval between two adjacent pumping pulses. For the experimental reports, CLD was defined so that the resonator had the lowest laser threshold for Δx=0 [1417].

Pulse width gain narrowing and intensity oscillation lead to the pulse compression of synchronously pumped solid-state Raman lasers. The dispersion can bring about pulse width gain narrowing and intensity modulation, and higher dispersion is beneficial for the pulse compression. The evidence can be found in the experiment reports. The Australia group used the 7.5 W 532nm 26 ps coherent pulsed laser [14] and 4.8 W 1064 nm 15 ps CW mode-locked Nd:YVO4 laser [16] pumped synchronously diamond Raman lasers, and minimum pulse duration of first Stokes obtained experimentally was 9 ps (pulse compression factor=2.9) and 9 ps (pulse compression factor=1.7) for appropriate CLD, respectively. Relative to infrared laser, frequency-doubled laser has higher dispersion, and can achieve higher pulse compression factor of first Stokes laser output [14].

4.4. Effects of normalized Raman gain G on output characteristics

Figure 6 shows effects of normalized Raman gain G on output characteristics of the first Stokes laser with RS=80% or 90%, Δζ = 0.0005, and other parameters given in Table 1.

 figure: Fig. 6.

Fig. 6. Effects of normalized Raman gain G on output characteristics of the first Stokes laser with RS=80% or 90%, Δζ = 0.0005, and other parameters given in Table 1.

Download Full Size | PDF

Normalized Raman gain G varies from 1×10−10 to10×10−10. At the beginning, ISinteg increases rapidly with the increase of G, and then slows down tending to a stable value as shown in Fig. 6(a). ISmax also goes up rapidly with the increase of G at the initial stage. However, further increasingly G brings about more intensive intensity modulation, and impedes the Raman conversion. Thus, ISmax decreases with the further increase of G. τS decreases rapidly with the increase of G initially, then decreases slowly tending to a stable value as shown in Fig. 6(b). It can be seen that the increased Raman gain enhances pulse width compression, and brings about more intensive intensity oscillation of the first Stokes laser, leading to the decrease of pulse width τS.

4.5 Effects of normalized dephasing time T on pulse compression

In the numerical calculation, the normalized dephasing time defined as Τ=T2/tp is varied according to the pulse duration of pumping laser tp with T2 of Raman mode to be constant. To investigate the influence of Τ on the pulse compression effect without being interfered by transient SRS gain reduction, the simulation is carried out for the same pumping pulse energy, i.e. G/Τ=2×10−9. Figure 7 shows the output characteristics of first Stokes laser for Τ=0.2 (G=4.0×10−10) and 0.5 (G=1.0×10−10), respectively. The normalized intensity integration ISinteg, normalized maximum intensity ISmax, and normalized pulse width τS of first Stokes pulses are shown in (a) and (b). ISinteg of Τ=0.5 is slightly higher than that of Τ=0.2. Around Δζ = 0.0005, the smallest pulse width of first Stokes is obtained, and is similar for two dephasing time conditions. However, for Τ=0.5, the maximum value of ISmax is 5.51, and is much larger than that of 3.09 for Τ=0.2. The parameters of two working point1 and 2 are also listed in Table 2 for Τ=0.5 (G=1.0×10−10, RS=90%).

 figure: Fig. 7.

Fig. 7. Effects of normalized dephasing time Τ on pulse compression of the first Stokes laser with RS=90%, Τ=0.2 (G=4.0×10−10) or 0.5 (G=1.0×10−9), and other parameters given in Table 1.

Download Full Size | PDF

Tables Icon

Table 2. Output parameters of first Stokes laser under working point 1 and 2

Figure 8 shows normalized intensity temporal distributions of first Stokes laser output and residual pumping laser for different Τ with RS=90% and Δζ = 0.0005 for the same pumping pulse energy, i.e. G/Τ=2×10−9. Τ is equal to 0.1, 0.2, 0.5 and 1 in (a), (b), (c) and (d), respectively. Intensities of first Stokes laser pulse are modulated, and the first peak of the oscillations usually has the highest intensity. As Τ increases from 0.1 to 1, the transient effect gets more serious, and the first peak possess more energy and higher intensity. The relation between pulse width of first Stokes laser and dephasing time is somewhat complicated. As Τ increases from 0.1 to 0.2, the pulse width decreases. For Τ=0.5 and 1, however, the pulse width increases due to serious dispersion pulse broadening in the highly transient regime.

 figure: Fig. 8.

Fig. 8. For the same pumping pulse energy G/Τ=2×10−9, normalized intensity temporal distributions of first Stokes laser output and residual pumping laser for different Τ with RS=90% and Δζ = 0.0005. Τ is equal to 0.1, 0.2, 0.5 and 1 in (a), (b), (c) and (d), respectively.

Download Full Size | PDF

With the same pumping pulse energy, smaller pulse width of pumping laser favor the generation of first Stokes laser pulse of high intensity. On the other hand, in the highly transient regime, pulse broadening due to group velocity dispersion will get serious, and limits the smallest pulse width of the first Stokes laser. Dispersion compensation is needed to realize shorter pulse width.

5. Conclusions

In this paper, the synchronously pumped solid-state Raman laser is analyzed systematically by numerically solving the normalized nonlinear coupling wave equations in the transient regime. There are two working points interesting to the actual applications. Under working point 1 with relative small negative CLD, the maximum Raman conversion efficiency can be achieved. Under working point 2 with relative small positive CLD, Stokes laser pulse is compressed to have the maximum peak power and minimum pulse. Actually, the first Stokes and the pumping pulses can propagate basically synchronously under different dispersion conditions. The experimental conditions, such as output coupling rate, Raman gain, dephasing time and dispersion, have little effect on the CLD of the working point 1 and 2.

The pulse compression of first Stokes laser stems from the pulse width gain narrowing and intensity oscillation effects. Thus, the early result in [25] is confirmed, which is pulse compression comes from a combination of group velocity walk-off and strong pump pulse depletion. Pulse width gain narrowing utilizes the group velocity difference of pumping and first Stokes laser due to linear dispersion under appropriate CLD. The intensity oscillation further narrows the pulse width of first Stokes laser, and clamp its peak intensity.

As the dispersion increases, the SRS conversion efficiency is improved. The higher dispersion or Raman gain leads to more intensive pulse width gain narrowing and intensity modulation, and favors the pulse compression. However, the dispersion brings about pulse width broadening, and limits the narrowest pulse width of first Stokes laser. The increased Raman gain enhances pulse width gain narrowing and intensity oscillation of the first Stokes laser leading to the decrease of pulse width τS, but too high Raman gain leads to the decrease of maximum peak pulse intensity. As Τ increases, the transient effect gets more serious, and the first peak of modulated first Stokes pulse possess more energy and higher intensity. Smaller pulse width of pumping laser favor the generation of first Stokes laser pulse of high intensity and small pulse width. On the other hand, in the highly transient regime, pulse broadening due to group velocity dispersion will get serious, and dispersion compensation is needed to realize shorter pulse width.

Funding

Natural Science Foundation of Shandong Province (ZR2018LF014).

Disclosures

The authors declare that there are no conflicts of interest.

References

1. J. Berzins, S. Indrisiunas, K. van Erve, A. Nagarajan, S. Fasold, M. Steinert, G. Gerini, P. Gecys, T. Pertsch, S. M. B. Baumer, and F. Setzpfandt, “Direct and high-throughput fabrication of Mie-resonant metasurfaces via single-pulse laser interference,” ACS Nano 14(5), 6138–6149 (2020). [CrossRef]  

2. J. Cheng, C. Liu, S. Shang, D. Liu, W. Perrie, G. Dearden, and K. Watkins, “A review of ultrafast laser materials micromachining,” Opt. Laser Technol. 46, 88–102 (2013). [CrossRef]  

3. S. Hu, C. Yang, S. Chang, Y. Huang, Y. Lin, and M. Lee, “Efficacy and safety of the picosecond 755-nm alexandrite laser for treatment of dermal pigmentation in Asians-a retrospective study,” Lasers Med. Sci. 35(6), 1377–1383 (2020). [CrossRef]  

4. C. Yan, D. Xu, Y. Wang, Z. Wang, Z. Wei, L. Tang, Y. He, J. Li, K. Zhong, W. Shi, and J. Yao, “Enhanced terahertz wave generation via Stokes wave recycling in non-synchronously picosecond pulse pumped gerahertz source,” IEEE Photonics J. 11(5), 1–8 (2019). [CrossRef]  

5. Y. C. Liu, C. M. Chen, J. Q. Hsiao, Y. Y. Pan, C. H. Tsou, H. C. Liang, and Y. F. Chen, “Compact efficient high-power triple-color Nd:YVO4 yellow-lime-green self-Raman lasers,” Opt. Lett. 45(5), 1144–1147 (2020). [CrossRef]  

6. Y. F. Chen, K. Y. Chen, Y. C. Liu, C. M. Chen, C. H. Tsou, and H. C. Liang, “Criterion for optimizing high-power acousto-optically Q-switched self-Raman yellow lasers with repetition rates up to 500 kHz,” Opt. Lett. 45(7), 1922–1925 (2020). [CrossRef]  

7. Z. Ding, S. Ding, H. Jia, J. Liu, L. Yang, J. Huangfu, and S. Wang, “Laser diode pumped, actively Q-switched, and mode-locked Nd:YAG/PbWO4 Raman laser,” Appl. Opt. 54(17), 5375–5381 (2015). [CrossRef]  

8. S. Ding, M. Wang, S. Wang, and W. Zhang, “Investigation on LD end-pumped passively Q-switched c-cut Nd: YVO4 self-Raman laser,” Opt. Express 21(11), 13052–13061 (2013). [CrossRef]  

9. S. Ding, W. Zhang, S. Wang, X. Wang, J. Zhang, and M. Wang, “Theoretical and experimental study on passively Q-switched intracavity frequency-doubled solid-state yellow Raman lasers,” Appl. Opt. 52(15), 3583–3590 (2013). [CrossRef]  

10. H. M. Pask, “The design and operation of solid-state Raman lasers,” Prog. Quantum Electron. 27(1), 3–56 (2003). [CrossRef]  

11. M. J. Colles, “Ultrashort pulse formation in a short-pulse-stimulated Raman oscillator,” Appl. Phys. Lett. 19(2), 23–25 (1971). [CrossRef]  

12. R. L. Carman, F. Shimizu, C. S. Wang, and N. Bloembergen, “Theory of Stokes pulse shapes in transient stimulated Raman scattering,” Phys. Rev. A 2(1), 60–72 (1970). [CrossRef]  

13. E. Granados, H. M. Pask, E. Esposito, G. McConnell, and D. J. Spence, “Multi-wavelength, all-solid-state, continuous wave mode locked picosecond Raman laser,” Opt. Express 18(5), 5289–5294 (2010). [CrossRef]  

14. D. J. Spence, E. Granados, and R. P. Mildren, “Mode-locked picosecond diamond Raman laser,” Opt. Lett. 35(4), 556–558 (2010). [CrossRef]  

15. E. Granados, D. J. Spence, and R. P. Mildren, “Deep ultraviolet diamond Raman laser,” Opt. Express 19(11), 10857–10863 (2011). [CrossRef]  

16. A. M. Warrier, J. P. Lin, H. M. Pask, R. P. Mildren, D. W. Coutts, and D. J. Spence, “Highly efficient picosecond diamond Raman laser at 1240 and 1485 nm,” Opt. Express 22(3), 3325–3333 (2014). [CrossRef]  

17. E. Granados, H. M. Pask, and D. J. Spence, “Synchronously pumped continuous-wave mode locked yellow Raman laser at 559 nm,” Opt. Express 17(2), 569–574 (2009). [CrossRef]  

18. M. Murtagh, J. Lin, R. P. Mildren, G. McConnell, and D. J. Spence, “Efficient diamond Raman laser generating 65 fs pulses,” Opt. Express 23(12), 15504–15513 (2015). [CrossRef]  

19. M. Franka, S. N. Smetaninb, M. Jelíneka, D. Vyhlídala, A. A. Kopalkinb, V. E. Shukshinb, L. I. Ivlevab, P. G. Zverevb, and V. Kubeček, “Synchronously-pumped all-solid-state SrMoO4 Raman laser generating at combined vibrational Raman modes with 26-fold pulse shortening down to 1.4 ps at 1220 nm,” Opt. Laser Technol. 111, 129–133 (2019). [CrossRef]  

20. M. Frank, M. Jelínek Jr, D. Vyhlídal, V. Kubeček, L. I. Ivleva, P. G. Zverev, and S. N. Smetanin, “Multi-wavelength picosecond BaWO4 Raman laser with long and short Raman shifts and 12-fold pulse shortening down to 3 ps at 1227 nm,” Laser Phys. 28(2), 025403 (2018). [CrossRef]  

21. M. Franka, S. N. Smetaninb, M. Jelíneka, D. Vyhlídala, L. I. Ivlevab, P. G. Zverevb, and V. Kubeček, “Highly efficient picosecond all-solid-state Raman laser at 1179 and 1227 nm on single and combined Raman lines in a BaWO4 crystal,” Opt. Lett. 43(11), 2527–2530 (2018). [CrossRef]  

22. M. Franka, S. N. Smetaninb, M. Jelíneka, D. Vyhlídala, V. E. Shukshinb, L. I. Ivlevab, P. G. Zverevb, and V. Kubečeka, “Efficient synchronously-pumped all-solid-state Raman laser at 1178 and 1227 nm on stretching and bending anionic group vibrations in a SrWO4 crystal with pulse shortening down to 1.4 ps,” Opt. Laser Technol. 119, 105660 (2019). [CrossRef]  

23. W. H. Lowdermilk and G. I. Kachen, “Coherent transient Raman pulse propagation,” Opt. Commun. 18(1), 68–69 (1976). [CrossRef]  

24. A. I. Vodchits, A. G. Shvedko, and V. A. Orlovich, “Stimulated Raman amplification of ultrashort seed pulses in compressed methane,” J. Opt. Soc. Am. B 22(2), 453–458 (2005). [CrossRef]  

25. E. Granados and D. J. Spence, “Pulse compression in synchronously pumped mode locked Raman lasers,” Opt. Express 18(19), 20422–20427 (2010). [CrossRef]  

26. G. Boyd, W. Johnston, and I. Kaminow, “Optimization of the stimulated Raman scattering threshold,” IEEE J. Quantum Electron. 5(4), 203–206 (1969). [CrossRef]  

27. A. Penzkofer, A. Laubereau, and W. Kaiser, “High intensity Raman interactions,” Prog. Quantum Electron. 6(2), 55–140 (1979). [CrossRef]  

28. M. G. Raymer and J. Mostowski, “Stimulated Raman scattering: Unified treatment of spontaneous initiation and spatial propagation,” Phys. Rev. A 24(4), 1980–1993 (1981). [CrossRef]  

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 (8)

Fig. 1.
Fig. 1. Cavity setup of a synchronously pumped solid-state Raman laser
Fig. 2.
Fig. 2. Effects of normalized CLD Δζ on output characteristics of first Stokes laser with RS=80% or 90%, G=4.0×10−10 or 8.0×10−10, and other parameters given in Table 1.
Fig. 3.
Fig. 3. Normalized intensity temporal distributions of first Stokes laser output and residual pumping laser for different Δζ with G=8.0×10−10, RS=80% and other parameters given in Table 1. Δζ is equal to 0.006, 0.0009, 0, -0.008, -0.02 and -0.04 in (a), (b), (c), (d), (e) and (f), respectively.
Fig. 4.
Fig. 4. Effects of dispersion on pulse compression of the first Stokes laser with RS = 90%, G=1.0×10−9, Τ=0.5 and other parameters given in Table 1. Solid lines are the results with dispersion, and dashed lines are those without dispersion.
Fig. 5.
Fig. 5. Effects of dispersion (VS-VL) on pulse compression of the first Stokes laser with RS=90%, Δζ = 0.0005 and other parameters given in Table 1. VS-VL is equal to 0, 0.0013, 0.0026 and 0.0052 in (a), (b), (c) and (d), respectively.
Fig. 6.
Fig. 6. Effects of normalized Raman gain G on output characteristics of the first Stokes laser with RS=80% or 90%, Δζ = 0.0005, and other parameters given in Table 1.
Fig. 7.
Fig. 7. Effects of normalized dephasing time Τ on pulse compression of the first Stokes laser with RS=90%, Τ=0.2 (G=4.0×10−10) or 0.5 (G=1.0×10−9), and other parameters given in Table 1.
Fig. 8.
Fig. 8. For the same pumping pulse energy G/Τ=2×10−9, normalized intensity temporal distributions of first Stokes laser output and residual pumping laser for different Τ with RS=90% and Δζ = 0.0005. Τ is equal to 0.1, 0.2, 0.5 and 1 in (a), (b), (c) and (d), respectively.

Tables (2)

Tables Icon

Table 1. Estimated values of the normalized parametersa

Tables Icon

Table 2. Output parameters of first Stokes laser under working point 1 and 2

Equations (23)

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

Q t v S Q x + 1 T 2 Q = κ Q T 2 E L E S + F ^ ( x , t ) ,
i 2 k S 2 E S + 1 v S E S t + γ S 2 E S = g S κ S Q E L ,
i 2 k L 2 E L + ( 1 v S v L ) E L x + 1 v L E L t + γ L 2 E L = g S κ L Q E S .
κ Q = i c μ L μ S 16 π ω S , κ S = i ω S μ S , κ L = i ω L μ L .
F ^ ( x , t ) F ^ ( x , t ) = g S μ L 2 A S ω S T 2 2 δ ( x x ) δ ( t t ) .
F ^ ( x , t ) = 0.
F ^ ( x , t ) F ^ ( x , t ) = 0.
E L , S ( r , x , t ) = E 0 L , S ( 0 , t ) 1 + x 2 Z R L , S 2 e x p [ r 2 w L , S 2 ( x ) ] ,
w L , S ( x ) = w 0 L , S 1 + x 2 / Z R L , S 2 ,
Z R L , S = π w 0 L , S 2 / λ L , S ,
E L i n ( r , x , t ) = E 0 L 0 ( 0 , t ) 1 + x 2 Z R L 2 exp [ r 2 w L 2 ( x ) ] ,
E 0 L 0 ( 0 , t ) = E p m a x exp [ 2 l n 2 ( t t 0 ) 2 t p 2 ] ,
Q ( r , x , t ) = Q ( 0 , x , t ) e x p ( r 2 w L 2 + r 2 w S 2 ) ,
2 E = 1 r r ( r E r ) + 1 r 2 2 E θ 2 .
τ = t / t p , ζ = x / Z R L , E L = E 0 L ( 0 , t ) / E p m a x , E S = E 0 S ( 0 , t ) / E p m a x ,
Θ = Q ( 0 , x , t ) / E p m a x 2 , G = g S E p m a x 2 Z R L ,
V L = v L t p / Z R L , V S = v S t p / Z R L , T = T 2 / t p ,
Θ τ V S Θ ζ V S 2 ζ ( 1 + ζ 2 ) Θ + 1 T Θ = κ Q T 1 ( 1 + ζ 2 ) E L E S + F ^ ( ζ , τ ) ,
E S τ + Γ S E S = 1 3 κ S G V S Θ E L ,
( V L V S ) E L ζ + ( V L V S ) ζ ( 1 + ζ 2 ) E L + E L τ + Γ L E L = 1 3 κ L G V L Θ E S ,
Γ L = γ L 2 V L Z R L , Γ S = γ S 2 V S Z R L , F ^ ( ζ , τ ) = F ^ ( x , t ) t p / E p m a x 2 .
I S = E S E S .
I S i n t e g = ( E S E S ) d τ ,
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.