Skip to main content Accessibility help


  • Access
  • Open access


MathJax is a JavaScript display engine for mathematics. For more information see
      • Send article to Kindle

        To send this article to your Kindle, first ensure is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

        Note you can select to send to either the or variations. ‘’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

        Find out more about the Kindle Personal Document Service.

        High-repetition-rate and high-power picosecond regenerative amplifier based on a single bulk Nd:GdVO4 crystal
        Available formats

        Send article to Dropbox

        To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

        High-repetition-rate and high-power picosecond regenerative amplifier based on a single bulk Nd:GdVO4 crystal
        Available formats

        Send article to Google Drive

        To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

        High-repetition-rate and high-power picosecond regenerative amplifier based on a single bulk Nd:GdVO4 crystal
        Available formats
Export citation


We report on a high-repetition-rate, high-power continuously pumped Nd:GdVO4 regenerative amplifier. Numerical simulations successfully pinpoint the optimum working point free of bifurcation instability with simultaneous efficient energy extraction. At a repetition rate of 100 kHz, a maximum output power of 23 W was obtained with a pulse duration of 27 ps, corresponding to a pulse energy of $230~\unicode[STIX]{x03BC}\text{J}$ . The system displayed an outstanding stability with a root mean square power noise as low as 0.3%. The geometry of the optical resonator and the pumping scheme enhanced output power in the $\text{TEM}_{00}$ mode with a single bulk crystal. Accordingly, nearly diffraction-limited beam quality was produced with $M^{2}\approx 1.2$ at full pump power.

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[15]. 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[69].

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[10]. 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[11]. 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[1214]. 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[15].

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[16]. This bifurcation-termed phenomenon constitutes a major impediment for RAs operated at high repetition rates[3], where the operating parameters of the laser systems have to be carefully optimized to maintain equilibrium[1720]. 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[19, 21], pumping at higher intensity[17], operating at the second operating point outside the unstable region if the effect could not be thoroughly avoided[22], or working in a bistable regime with a negligible lower bifurcation branch[23]. 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[2429]. 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[30]. 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[31]. For high repetition rate reports, an average output power of 25.5 W at 200 kHz was realized in 2008 by Coherent Inc.[32]. 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[33]. 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[21]. 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. [14, 34], 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 [3542]. 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[35, 43]. 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. [3542]).

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[20]. 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[19, 20, 44]. 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. [19, 20, 44], 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[19, 44]. 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. [44], which is based on one-dimensional discrete-time dynamical systems[45]. 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[20]. 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[20, 44].

(1) $$\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}$$
(2) $$\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}$$
(3) $$\begin{eqnarray}\displaystyle G_{0} & = & \displaystyle R_{p}T_{1}\unicode[STIX]{x1D70E}L_{a},\end{eqnarray}$$
(4) $$\begin{eqnarray}\displaystyle g & = & \displaystyle N\unicode[STIX]{x1D70E}L_{a}/G_{0},\end{eqnarray}$$
(5) $$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D700} & = & \displaystyle \unicode[STIX]{x1D711}\unicode[STIX]{x1D70E}/(A_{a}G_{0}),\end{eqnarray}$$
(6) $$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D70F} & = & \displaystyle tG_{0}\unicode[STIX]{x1D6FD}/T_{0},\end{eqnarray}$$
(7) $$\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[19]. 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[46]. 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[47].

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[47]. 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.


This work was supported by the Ministry of Science and Technology of the People’s Republic of China (MOST) (No. 2017YFB0405202).


1. Luther-Davies, B. Kolev, V. Z. Lederer, M. J. Madsen, N. R. Rode, A. V. Giesekus, J. Du, K. M. and Duering, M. Appl. Phys. A 79, 1051 (2004).
2. Mero, M. Heiner, Z. Petrov, V. Rottke, H. Branchi, F. Thomas, G. M. and Vrakking, M. J. J. Opt. Lett. 43, 5246 (2018).
3. Fattahi, H. Schwarz, A. Geng, X. T. Keiber, S. Kim, D. E. Krausz, F. and Karpowicz, N. Opt. Express 22, 31440 (2014).
4. Ahrens, J. Prochnow, O. Binhammer, T. Lang, T. Schulz, B. Frede, M. and Morgner, U. Opt. Express 24, 8074 (2016).
5. Danilevicius, R. Zaukevicius, A. Budriunas, R. Michailovas, A. and Rusteika, N. Opt. Express 24, 17532 (2016).
6. Luhrmann, M. Theobald, C. Wallenstein, R. and L’Huillier, J. A. Opt. Express 17, 22761 (2009).
7. von Grafenstein, L. Bock, M. Ueberschaer, D. Griebner, U. and Elsaesser, T. Opt. Lett. 41, 4668 (2016).
8. Ueffing, M. Lange, R. Pleyer, T. Pervak, V. Metzger, T. Sutter, D. Major, Z. Nubbemeyer, T. and Krausz, F. Opt. Lett. 41, 3840 (2016).
9. Fattahi, H. Alismail, A. Wang, H. Brons, J. Pronin, O. Buberl, T. Vamos, L. Arisholm, G. Azzeer, A. M. and Krausz, F. Opt. Lett. 41, 1126 (2016).
10. Nubbemeyer, T. Kaumanns, M. Ueffing, M. Gorjan, M. Alismail, A. Fattahi, H. Brons, J. Pronin, O. Barros, H. G. Major, Z. Metzger, T. Sutter, D. and Krausz, F. Opt. Lett. 42, 1381 (2017).
11. Morrissey, F. X. Fan, T. Y. Miller, D. E. and Rand, D. Opt. Lett. 42, 707 (2017).
12. Pouysegur, J. Delaigue, M. Zaouter, Y. Honninger, C. Mottay, E. Jaffres, A. Loiseau, P. Viana, B. Georges, P. and Druon, F. Opt. Lett. 38, 5180 (2013).
13. Caracciolo, E. Kemnitzer, M. Guandalini, A. Pirzio, F. Agnesi, A. and der Au, J. A. Opt. Express 22, 19912 (2014).
14. Bai, Z. Fan, Z. Bai, Z. Lian, F. Kang, Z. and Lin, W. Appl. Sci. 5, 359 (2015).
15. Caracciolo, E. Guandalini, A. Pirzio, F. Kemnitzer, M. Kienle, F. Agnesi, A. and der Au, J. A. Proc. SPIE 10082, 100821E (2017).
16. Dorring, J. Killi, A. Morgner, U. Lang, A. Lederer, M. and Kopf, D. Opt. Express 12, 1759 (2004).
17. Kroetz, P. Ruehl, A. Calendron, A. L. Chatterjee, G. Cankaya, H. Murari, K. Kartner, F. X. Hartl, I. and Miller, R. J. D. Appl. Phys. B 123, 126 (2017).
18. von Grafenstein, L. Bock, M. and Griebner, U. IEEE J. Sel. Top. Quantum Electron. 24, 3000213 (2018).
19. Grishin, M. Gulbinas, V. and Michailovas, A. Opt. Express 17, 15700 (2009).
20. Grishin, M. Gulbinas, V. and Michailovas, A. Opt. Express 15, 9434 (2007).
21. Gao, P. Lin, H. Li, J. Guo, J. Yu, H. Zhang, H. and Liang, X. Opt. Express 24, 13963 (2016).
22. Kroetz, P. Ruehl, A. Chatterjee, G. Calendron, A. L. Murari, K. Cankaya, H. Li, P. Kartner, F. X. Hartl, I. and Miller, R. J. D. Opt. Lett. 40, 5427 (2015).
23. von Grafenstein, L. Bock, M. Steinmeyer, G. Griebner, U. and Elsaesser, T. Laser Photon. Rev. 10, 123 (2016).
24. Mecseki, K. Bigourd, D. Patankar, S. Stuart, N. H. and Smith, R. A. Appl. Opt. 53, 2229 (2014).
25. Okishev, A. V. Dorrer, C. Smirnov, V. I. Glebov, L. B. and Zuegel, J. D. Opt. Express 15, 8197 (2007).
26. Lin, H. Li, J. He, J. and Liang, X. Chin. Opt. Lett. 9, 081404 (2011).
27. Kornev, A. F. Balmashnov, R. V. Kuchma, I. G. Davtian, A. S. and Oborotov, D. O. Opt. Lett. 43, 4394 (2018).
28. Long, M. Chen, L. Chen, M. and Li, G. Appl. Phys. B 122, 142 (2016).
29. Noom, D. W. E. Witte, S. Morgenweg, J. Altmann, R. K. and Eikema, K. S. E. Opt. Lett. 38, 3021 (2013).
30. Bai, Z. Bai, Z. Yang, C. Chen, L. Chen, M. and Li, G. Opt. Laser Technol. 46, 25 (2013).
31. Peng, Z. Chen, M. Yang, C. Chang, L. and Li, G. Jpn. J. Appl. Phys. 54, 028001 (2015).
32. Clubley, D. A. Bell, A. S. and Friel, G. Proc. SPIE 6871, 68711D (2008).
33. Harth, F. Ulm, T. Luhrmann, M. Knappe, R. Klehr, A. Hoffmann, T. Erbert, G. and L’Huillier, J. A. Opt. Express 20, 7002 (2012).
34. Chen, M. Chang, L. Yang, C. Chen, L. and Li, G. Chin. J. Lasers 40, 0602010 (2013).
35. Zhang, B. Li, G. Chen, M. Zhang, Z. and Wang, Y. Opt. Lett. 28, 1829 (2003).
36. Vodchits, A. I. Orlovich, V. A. and Apanasevich, P. A. J. Appl. Spectrosc. 78, 918 (2012).
37. Koechner, W. Solid State Laser Engineering, 6th edition (Springer, New York, 2006).
38. Yu, H. Liu, J. Zhang, H. Kaminskii, A. A. Wang, Z. and Wang, J. Laser Photon. Rev. 8, 847 (2014).
39. Lavi, R. and Jackel, S. Appl. Opt. 39, 3093 (2000).
40. Lavi, R. Jackel, S. Tzuk, Y. Winik, M. Lebiush, E. Katz, M. and Paiss, I. Appl. Opt. 38, 7382 (1999).
41. Jensen, T. Ostroumov, V. G. Meyn, J. P. Huber, G. Zagumennyi, A. I. and Shcherbakov, I. A. Appl. Phys. B 58, 373 (1994).
42. McDonagh, L. Wallenstein, R. Knappe, R. and Nebel, A. Opt. Lett. 31, 3297 (2006).
43. Kleinbauer, J. Knappe, R. and Wallenstein, R. Appl. Phys. B 81, 163 (2005).
44. Grishin, M. and Michailovas, A. Advances in Solid State Lasers: Development and Applications (Intech, Vukovar, 2010).
45. Alligood, K. T. Sauer, T. D. and Yorke, J. A. Chaos: An Introduction to Dynamical Systems (Springer, New York, 1996).
46. Lin, H. Guo, J. Gao, P. Yu, H. and Liang, X. Opt. Express 24, 13957 (2016).
47. Sun, D. Lin, H. Guo, J. Wang, W. and Liang, X. Opt. Lett. 43, 4346 (2018).