1 Introduction
Compact, stable and high-power diode-pumped short-pulse laser amplifier systems with excellent spatial quality are ideal sources for high-power optical parametric chirped pulse amplification (OPCPA) and efficient laser processing[Reference Luther-Davies, Kolev, Lederer, Madsen, Rode, Giesekus, Du and Duering1–Reference Danilevicius, Zaukevicius, Budriunas, Michailovas and Rusteika5]. Among various configurations, regenerative amplifiers (RAs) are routinely used to enhance the output of mode-locked oscillators because of their ability to provide gains of several orders of magnitude and a resonator structure that maintains the spatial quality of the seed[Reference Luhrmann, Theobald, Wallenstein and L’Huillier6–Reference Fattahi, Alismail, Wang, Brons, Pronin, Buberl, Vamos, Arisholm, Azzeer and Krausz9].
In the past two decades, various solutions and architectures have been proposed for power scaling of regenerative amplifiers. The thin-disk geometry was particularly impressive for heat dissipation. Nubbemeyer et al. employed a Kerr-lens mode-locked thin-disk oscillator to seed two stages of regenerative amplifiers, including a thin-disk preamplifier and a main amplifier comprising two thin-disk modules, generating
${>}200~\text{mJ}$
pulses at a 5 kHz repetition rate and
${>}100~\text{mJ}$
at a 10 kHz repetition rate[Reference Nubbemeyer, Kaumanns, Ueffing, Gorjan, Alismail, Fattahi, Brons, Pronin, Barros, Major, Metzger, Sutter and Krausz10]. The cryogenic cooling technique offered another alternative to obtain a high pulse energy at a 1 kHz repetition rate, although measures were necessary to mitigate the reduced gain bandwidth[Reference Morrissey, Fan, Miller and Rand11]. Although not so aggressive in heat removal, laser systems based on bulk materials generally employ simpler structures and work more efficiently. These characteristics make bulk materials favored choices as laser media in high-power diode-pumped solid state lasers[Reference Pouysegur, Delaigue, Zaouter, Honninger, Mottay, Jaffres, Loiseau, Viana, Georges and Druon12–Reference Bai, Fan, Bai, Lian, Kang and Lin14]. An average output power of 34 W with pulses as short as 140 fs at a 500 kHz repetition rate has been achieved by a bulk Yb:CALGO regenerative amplifier with the implementation of the chirped pulse amplification (CPA) technique[Reference Caracciolo, Guandalini, Pirzio, Kemnitzer, Kienle, Agnesi and der Au15].
High-repetition-rate regenerative amplifiers are beneficial for promoting industrial throughput and reducing data accumulation or processing time in scientific applications. However, high repetition rates come with bifurcation and even chaotic pulse train dynamics, in which the pulse energy is unstable[Reference Dorring, Killi, Morgner, Lang, Lederer and Kopf16]. This bifurcation-termed phenomenon constitutes a major impediment for RAs operated at high repetition rates[Reference Fattahi, Schwarz, Geng, Keiber, Kim, Krausz and Karpowicz3], where the operating parameters of the laser systems have to be carefully optimized to maintain equilibrium[Reference Kroetz, Ruehl, Calendron, Chatterjee, Cankaya, Murari, Kartner, Hartl and Miller17–Reference Grishin, Gulbinas and Michailovas20]. A detailed analysis concerning this problem is elaborated in the following section. For more than a decade, several methods have been proposed to eliminate this deleterious effect. The methodologies include raising the seed energy[Reference Grishin, Gulbinas and Michailovas19, Reference Gao, Lin, Li, Guo, Yu, Zhang and Liang21], pumping at higher intensity[Reference Kroetz, Ruehl, Calendron, Chatterjee, Cankaya, Murari, Kartner, Hartl and Miller17], operating at the second operating point outside the unstable region if the effect could not be thoroughly avoided[Reference Kroetz, Ruehl, Chatterjee, Calendron, Murari, Cankaya, Li, Kartner, Hartl and Miller22], or working in a bistable regime with a negligible lower bifurcation branch[Reference von Grafenstein, Bock, Steinmeyer, Griebner and Elsaesser23]. Consequently, efforts to optimize system parameters are necessary to achieve stable and efficient performance.
RAs based on bulk neodymium gain media show specific advantages. A high stimulated-emission cross-section simplifies the system design, and requirements such as the reflectivity for optical components are allayed. RAs based on them are generally free of CPA because of the relatively longer pulse durations. Diode-pumped RAs built with Nd:YAG or Nd:YLF crystals are commonly side-pumped for high pulse energy generation[Reference Mecseki, Bigourd, Patankar, Stuart and Smith24–Reference Noom, Witte, Morgenweg, Altmann and Eikema29]. End-pumping schemes provide simpler structures and are always realized by neodymium-doped vanadate crystals. Lasers based on neodymium-doped vanadate crystals deliver moderate pulse energies as well as average power, but less complex cavities and pumping schemes are required. In 2013, a regenerative amplifier based on a Nd:YVO
$_{4}$
crystal delivering a 28 mJ pulse energy at a 1 kHz repetition rate and a 112.1 ps long pulse was reported[Reference Bai, Bai, Yang, Chen, Chen and Li30]. The whole setup was rather complex because a CPA configuration had to be used and the multi-pulse output limited its applications. Seeded by a cavity-dumped laser, a 10.5 W output power was generated with a Nd:YVO
$_{4}$
RA by the same group in 2015, although the repetition rate was just 10 kHz[Reference Peng, Chen, Yang, Chang and Li31]. For high repetition rate reports, an average output power of 25.5 W at 200 kHz was realized in 2008 by Coherent Inc.[Reference Clubley, Bell and Friel32]. In 2012, a continuous-wave (CW) beam of a distributed feedback (DFB) diode laser modulated by an ultra-fast semiconductor modulator was injected into a regenerative amplifier based on a Nd:YVO
$_{4}$
crystal to produce 47.7 W at 816 kHz and 26.4 W at 99.8 kHz, with pulse durations adjustable from 400 ps to 1 ns[Reference Harth, Ulm, Luhrmann, Knappe, Klehr, Hoffmann, Erbert and L’Huillier33]. However, this scheme suffered from amplitude fluctuations and CW background noise, which are disadvantageous for many applications. Among vanadate crystals, Nd:LuVO
$_{4}$
has the highest stimulated emission cross-section and was applied to achieve 25.1 W at 1.43 MHz and 20.5 W at 100 kHz by our group in 2016[Reference Gao, Lin, Li, Guo, Yu, Zhang and Liang21]. However, beam quality degradation in long-term operation was observed in this Nd:LuVO
$_{4}$
dual-crystal configuration RA, which may be attributed to the immature manufacturing technology of this commercially unavailable gain medium. In Refs. [Reference Bai, Fan, Bai, Lian, Kang and Lin14, Reference Chen, Chang, Yang, Chen and Li34], Nd:YVO
$_{4}$
crystals as long as 20 mm and 30 mm were used, respectively, to produce relatively low extraction efficiencies. The authors experimentally observed the bifurcation phenomena for their specific working conditions by adjusting the number of cavity round trips (NRTs) in the RA resonators; however, no detailed study on the physical process was carried out. As a result, numerous trials had to be implemented in practice until the output pulse became stable.
Here we report a 100 kHz high-power regenerative amplifier based on a single bulk Nd:GdVO4 crystal. The basic properties of three Nd-doped single crystals are compared in Table 1
[Reference Zhang, Li, Chen, Zhang and Wang35–Reference McDonagh, Wallenstein, Knappe and Nebel42]. Among these, Nd:GdVO4 is an attractive candidate for high-power lasers because of its considerably higher thermal conductivity, which endows it with an ideal combination of the merits of Nd:YVO
$_{4}$
and Nd:YAG[Reference Zhang, Li, Chen, Zhang and Wang35, Reference Kleinbauer, Knappe and Wallenstein43]. Numerical simulations were conducted to analyze the system stability and efficiency performance of the continuously pumped high-repetition-rate regenerative amplifier. In virtue of the analysis, optimized operating conditions were proposed and further realized in our experiment with a highly stable and efficient output. An average output power of 23 W at a repetition rate of 100 kHz was delivered with an excellent output beam quality, characterized by
$M^{2}\approx 1.2$
. To the best of our knowledge, this is the highest output power that has been generated at
${\sim}100~\text{kHz}$
repetition rates for bulk crystal end-pumped all-solid-state picosecond regenerative amplifiers with pulse durations below 30 ps.
Table 1. Basic properties of three Nd-doped single crystals (see Refs. [Reference Zhang, Li, Chen, Zhang and Wang35–Reference McDonagh, Wallenstein, Knappe and Nebel42]).
2 Numerical simulations
Looking at its working principle, one RA operation cycle can be divided into two successive stages: the pump stage and the amplification stage. No voltage is applied to the Pockels cell during the pump stage. The amplification stage starts when the seed is injected into the amplification cavity and a high voltage is applied to the Pockels cell. When the pulse energy reaches a certain level, the high voltage is switched off, the amplified pulse is ejected and the next cycle begins[Reference Grishin, Gulbinas and Michailovas20]. For continuously pumped regenerative amplifiers, this operation cycle division is justified when the amplification stage duration is considerably shorter than the inverse of the repetition rate so that the impact of the pump during amplification can be neglected. Nevertheless, for continuously pumped regenerative amplifiers, the operation cycles become interdependent when the pump stage duration becomes comparable to or even shorter than the upper state lifetime of the respective transition. Let us consider a 100 kHz RA based on a Nd:GdVO4 crystal as an example. The pump stage duration (
${<}10~\unicode[STIX]{x03BC}\text{s}$
) is shorter than the upper state lifetime of the
$\,^{4}\text{F}_{3/2}{-}\,^{4}\text{I}_{11/2}$
transition (
$90~\unicode[STIX]{x03BC}\text{s}$
). Then the population inversion accumulated during one operation cycle cannot be fully dissipated by spontaneous emission and nonradiative processes to reset the inversion state before the next pulse arrives. As a result, the output pulse energy depends on the previous energy extraction condition. In this situation, a stable output pulse energy can still be obtained by an appropriate choice of system parameters. A preliminary picture including crucial system parameters will be helpful to achieve a stable and efficient operating performance.
Numerical simulations were performed to fully exploit the system’s capability: to ensure stable operation and extract as much stored energy as possible. Diagrams presented in the parameter space are very helpful to realize this prospect[Reference Grishin, Gulbinas and Michailovas19, Reference Grishin, Gulbinas and Michailovas20, Reference Grishin and Michailovas44]. The parameter space separatrix and curve of
$\text{NRT}^{\text{MAX}}$
are depicted in the same diagram. The separatrix divides the parameter space into two sections: corresponding to the stable and the bifurcation regimes. As described in Refs. [Reference Grishin, Gulbinas and Michailovas19, Reference Grishin, Gulbinas and Michailovas20, Reference Grishin and Michailovas44], the separatrix shrinks to an ellipse for higher seed energy.
$\text{NRT}^{\text{MAX}}$
presents the number of cavity round trips, potentially providing the highest output energy[Reference Grishin, Gulbinas and Michailovas19, Reference Grishin and Michailovas44]. Optimum operation was reached when the curve of
$\text{NRT}^{\text{MAX}}$
was outside the bifurcation region. The aforementioned parameter space includes the repetition rate and NRT. Other fundamental factors influencing the operating regimes involve the seed energy, pumping intensity and intracavity loss. Our simulations are established on the model developed by Grishin et al.
[Reference Grishin and Michailovas44], which is based on one-dimensional discrete-time dynamical systems[Reference Alligood, Sauer and Yorke45]. This model is found to have a fast computation speed and is adequate for our picosecond system with a narrow spectrum bandwidth. The parameters used in our simulations are listed in Table 2. These parameters are consistent with the conditions we planned to utilize in our experiment.
The basic equations describing the pump and amplification stages are listed as Equations (1) and (2), which were re-arranged from rate equations[Reference Grishin, Gulbinas and Michailovas20]. For the pump stage, the initial and final gains were connected by the analytical solution to this set of equations, because no lasing was assumed. The cavity loss was included in the amplification stage, and solutions to Equations (1) and (2) were calculated numerically. The definitions and expressions for nondimensional terms are listed after the equations[Reference Grishin, Gulbinas and Michailovas20, Reference Grishin and Michailovas44].
$$\begin{eqnarray}\displaystyle {\displaystyle \frac{\text{d}\unicode[STIX]{x1D700}(\unicode[STIX]{x1D70F})}{\text{d}\unicode[STIX]{x1D70F}}} & = & \displaystyle \unicode[STIX]{x1D700}(\unicode[STIX]{x1D70F})[g(\unicode[STIX]{x1D70F})-g_{t}],\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle {\displaystyle \frac{\text{d}g(\unicode[STIX]{x1D70F})}{\text{d}\unicode[STIX]{x1D70F}}} & = & \displaystyle -\unicode[STIX]{x1D700}(\unicode[STIX]{x1D70F})g(\unicode[STIX]{x1D70F})+\frac{1-g(\unicode[STIX]{x1D70F})}{\unicode[STIX]{x1D70F}_{1}},\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle G_{0} & = & \displaystyle R_{p}T_{1}\unicode[STIX]{x1D70E}L_{a},\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle g & = & \displaystyle N\unicode[STIX]{x1D70E}L_{a}/G_{0},\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D700} & = & \displaystyle \unicode[STIX]{x1D711}\unicode[STIX]{x1D70E}/(A_{a}G_{0}),\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D70F} & = & \displaystyle tG_{0}\unicode[STIX]{x1D6FD}/T_{0},\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle g_{t} & = & \displaystyle -\ln (1-l)/G_{0}.\end{eqnarray}$$
Here
$R_{p}$
is the pumping rate,
$L_{a}$
is the gain medium length,
$N$
is the population inversion density,
$\unicode[STIX]{x1D711}$
is the photon number,
$A_{a}$
is the mode cross-section in the gain medium,
$\unicode[STIX]{x1D6FD}$
is the number of passes through the gain medium during one round trip (
$\unicode[STIX]{x1D6FD}=2$
for linear cavity;
$\unicode[STIX]{x1D6FD}=1$
for ring cavity),
$T_{0}$
is the cavity round trip time,
$t$
is the time,
$\unicode[STIX]{x1D70F}_{1}$
is the normalized upper state lifetime (
$\unicode[STIX]{x1D70F}_{1}=T_{1}G_{0}\unicode[STIX]{x1D6FD}/T_{0}$
),
$\unicode[STIX]{x1D700}$
is the normalized energy, and
$\unicode[STIX]{x1D70F}$
is the normalized round trip number. The definitions of other parameters can be found in Table 2.
Table 2. Key parameter values for simulation.
Figure 1. Parameter separatrix (colored blue and yellow) and curve of
$\text{NRT}^{\text{MAX}}$
(colored bold black). The red star indicates the optimum working point for 100 kHz repetition rate operation.
In the simulation, the repetition rate of the RA operation range was set between 1 kHz and 150 kHz. For the seed energy, the mode mismatch between the seed and the RA in the spatial and spectral domains should be taken into account[Reference Grishin, Gulbinas and Michailovas19]. The effective seed energy was estimated to be
${\sim}40\%$
of the seed energy available (800 nJ). The simulation results are illustrated in Figure 1. The lower-right elliptical-like separatrix splits the stable and bifurcation regions, the inner of which is the bifurcation region. That is, if a combination of system parameters (repetition rate and round trip number) is chosen outside of the region surrounded by the elliptical-like separatrix then stable operation would be expected. The separatrix has a jagged appearance because in reality the NRT should be integers; this is also the reason for the form of the bold black curve of
$\text{NRT}^{\text{MAX}}$
. Diagrams of this type help us to estimate the optimum NRT when a specific repetition rate is required to exclude ranges inappropriate for operation. It can be seen from the diagram that outside a repetition rate range of 56–78 kHz, optimum operation could be achieved in theory. Consequently, with the seed source available, efficient and stable operation at a 100 kHz repetition rate is attainable at a predicted round trip number of 10, which is indicated with the red star.
3 Experimental setup
The experimental setup of the regenerative amplifier system is shown in Figure 2. The front-end seed source was based on a robust home-built, continuous wave passively mode-locked picosecond Nd:GdVO4 oscillator delivering an output power of 7.6 W at a 7.7 MHz repetition rate. The oscillator avoided the complexity of a preamplifier and delivered excellent beam quality and long-term stability, ensuring reliable and flexible seeding of the amplifier with ample energy. An optical isolator was also included to protect the oscillator against potential feedback from the amplifier. Then a RbTiOPO
$_{4}$
(RTP) pulse picker reduced the repetition rate to 100 kHz to enhance pulse contrast and suppress parasitic consumption of the stored energy. In addition, the lens between M2 and M3 matches the beam divergence and waist location between the modes of the oscillator and RA.
For the amplifier, the gain medium was a 0.5 at.% doped Nd:GdVO4 crystal with dimensions 4 mm (a)
$\times$
4 mm (c)
$\times$
15 mm (a). Each polished facet of the crystal was antireflection coated at 879 nm and 1063 nm. The crystal was wrapped with indium foil and mounted in a water-cooled copper heatsink. The cooling water temperature was maintained at
$15\,^{\circ }\text{C}$
during operation. A linear polarized pump was imaged onto the Nd:GdVO4 crystal of diameter
$1100~\unicode[STIX]{x03BC}\text{m}$
by a relay imaging system, in which a
$\unicode[STIX]{x1D70E}$
polarized pumping configuration was adopted[Reference Lin, Guo, Gao, Yu and Liang46]. The pumping wavelength was centered at 879 nm. Approximately 85% of the pump power was absorbed during a single pass through the crystal.
The RA cavity comprised a dichroic mirror M7, a quarter-wave plate (QWP), a thin-film polarizer (TFP), a Pockels cell (PC) based on BBO crystal with an aperture of 5 mm, and two highly reflective concave mirrors M5 (
$R=600~\text{mm}$
) and M6 (
$R=300~\text{mm}$
) along with two plane mirrors. Concave mirrors M5 and M6 were applied to get an appropriate beam size in the PC. The total cavity length of 1.96 m corresponds to a round trip time of about 13.1 ns, which was adapted to the rise and fall times (about 10 ns) of the PC driver. Although the high-repetition-rate RA benefits productivity in industrial and scientific applications, it also poses challenges for the speed of high-voltage devices. A repetition rate of 100 kHz was chosen for the investigation because it is a suitable compromise.
Figure 2. Schematic of the experimental setup: PP, pulse picker; TFP, thin-film polarizer; FR, Faraday rotator; HWP, half-wave plate; QWP, quarter-wave plate; PC, Pockels cell. The beam inside the RA cavity propagates along the 15 mm long
$a$
-axis of the Nd:GdVO4 crystal.
Figure 3. CW output power versus absorbed pump power.
4 Experimental results and discussions
Continuous wave (CW) operation was first conducted to confirm the optimum cavity performance, in which the combination of TFP2 and QWP served as the output coupler and the PC driver was switched off. Figure 3 shows the continuous-wave output power versus the absorbed pump power. A maximum output power of 28 W was achieved at 66 W absorbed pump power, which corresponds to 42% optical efficiency. The moderate efficiency was due to losses caused by the insertion of the PC, as the output power increased to 33 W once the PC was removed. The roll-over phenomenon in Figure 3 was attributed to thermal lensing of the Nd:GdVO4 crystal, which destabilized the cavity. A similar phenomenon has been observed previously[Reference Sun, Lin, Guo, Wang and Liang47].
In the RA operation regime, the seed with 800 nJ pulse energy was injected into the RA for amplification. Its output power is depicted in Figure 4(a). The onset of the parasitic cavity-dumped background was observed when the absorbed pump power was more than 60 W. Therefore, the operation point was fixed around 60 W pumping power, which generated an average output power of 23 W. A fast photodiode was used to monitor the intracavity signal through the leakage radiation behind M4, the evolution of which is shown in Figure 4(b) (only the last five signals could be observed in this way). The optimum round trip number was experimentally found to be 10, which agreed well with the previous numerical simulations (indicated with the red star in Figure 1). System operation under altered conditions was also surveyed. When the round trip numbers were set as 9 and 11, the output power decreased by 300 mW and 500 mW, respectively; however, pulse trains remained stable. Nevertheless, when the round trip number was set as 13, the system started converting to the period-doubling regime. These performances confirmed the feasibility of our numerical simulations.
Figure 4. (a) RA regime output power versus absorbed pump power; (b) the last five intracavity signals of the RA.
The temporal characteristics of the oscillator and RA output pulse were confirmed by an intensity autocorrelator. As shown in Figure 5, the pulse duration (FWHM) of the RA was determined to be 27 ps, which was almost the same as that of the seed. We also tried to measure the spectral properties of the seed and RA output separately. They were observed to be
${<}0.1~\text{nm}$
in width. Because of the limited resolution of our spectrometer (
${\sim}0.04~\text{nm}$
), accurate spectral widths and spectra of the oscillator and RA output are not provided here. Additionally, the calculated B-integral amounts to 2.1 rad, accounting for the polarizer, Pockels cell, and crystal. This value can be reduced to 1.9 rad if the number of cavity round trips is reduced to 9.
Figure 5. Intensity autocorrelation traces of the oscillator and RA output.
The long-term stability of the RA system at full pump power for 30 min is characterized and presented in Figure 6(a). The fluctuations of the average output power amount to 0.3% and the pulse-to-pulse stability was measured to be 1.7% RMS (root mean square), in which 1000 consecutive pulses were recorded within a time window of 10 ms. The beam profile is demonstrated in Figure 6(b) with
$M^{2}\approx 1.2$
at the maximum output power. The inset in Figure 6(b) displays the far-field beam profile at the focus.
Figure 6. (a) Long-term power stability measurement of the RA output. (b) RA output beam quality.
5 Conclusion
In conclusion, we have demonstrated a 100 kHz high-power continuously pumped Nd:GdVO4 regenerative amplifier. The numerical analysis of the continuously pumped high-repetition-rate RA prior to the experiment facilitated optimization of the parameters in our experiment. This helped eliminate bifurcation instability and achieve efficient energy extraction. With a single bulk crystal, a maximum output pulse energy of
$230~\unicode[STIX]{x03BC}\text{J}$
at a repetition rate of 100 kHz was obtained with a pulse duration of 27 ps and root mean square power noise as low as 0.3%. The
$\unicode[STIX]{x1D70E}$
polarized pumping scheme and careful design of the optical resonator enhanced the output power in the
$\text{TEM}_{00}$
mode with a single bulk crystal. Accordingly, a nearly diffraction-limited beam quality was produced with
$M^{2}\approx 1.2$
at full pump power. This is, to the best of our knowledge, the highest output power that has been attained at
${\sim}100~\text{kHz}$
repetition rates for a bulk crystal end-pumped all-solid-state picosecond regenerative amplifier with pulse durations below 30 ps. Power scaling and beam quality improvement are feasible by applying a double-crystal orthogonal thermal compensation scheme to distribute the thermal load[Reference Sun, Lin, Guo, Wang and Liang47]. Because the beam size needs adapting to avoid damage, the lack of an appropriate optical imaging system for the pump prevented us from implementing further amplification based on the double-crystal geometry. Work based on this scheme is underway.
Acknowledgements
This work was supported by the Ministry of Science and Technology of the People’s Republic of China (MOST) (No. 2017YFB0405202).
Open access
