Abstract
The performance of regenerative amplifiers at high repetition rates is often limited by the occurrence of bifurcations induced by a destabilization of the pulse-to-pulse dynamics. While bifurcations can be suppressed by increasing the seed energy using dedicated pre-amplifiers, the availability of adjustable filters and control electronics in modern pulse amplifiers allows to exploit feedback strategies to cope with these instabilities. In this paper, we present a theoretical and experimental analysis of active feedback methods to stabilize otherwise unstable operational regimes of regenerative amplifiers. To this end, the dynamics of regenerative amplifiers are investigated starting from a general space-dependent description to obtain a generalization of existing models from the literature. Suitable feedback strategies are then developed utilizing measurements of the output pulse energies or the transmitted pump light, respectively. The effectiveness of the proposed approach is highlighted by experimental results for a Yb:CaF2-based regenerative amplifier.
Published by The Optical Society under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.
1. Introduction
Regenerative amplifiers (RAs) are a widely used tool for the generation of high-intensity laser pulses. Combined with chirped pulse amplification techniques to reduce the B-integral and avoid amplifier optical damage, RAs are used to create ultra-fast high energy pulses with a wide range of potential applications such as strong field physics [1], coherent control [2], and ablation-based material processing [3] including various medical applications. Although high repetition rates are considered beneficial for several applications [4], RAs are known to become unstable quite easily for a range of repetition rates around the inverse life time of the laser medium’s upper lasing energy level [5]. When operated within the unstable region, the sequence of output pulses either converges to a periodic orbit whereby the limit set exhibits so-called period-doubling bifurcations or remains aperiodic in the case of deterministic chaos. While it is always possible to find stable operating points at a desired repetition rate, the operating points with highest power output and highest energy efficiency are typically located in the vicinity of the unstable region.
As a result, this topic has been of significant scientific interest since the first experimental results were documented in [6]. A numerical study using a simple rate equation model in [4] highlights the coupling mechanism of subsequent pulses that leads to period-doubling bifurcations and ultimately deterministic chaos. By using analytic solutions of the rate equation model under quite restrictive assumptions, a detailed parameter study of the underlying problem in [7] and [8] finally shows that the bifurcations can be suppressed by increasing the seed energy. Although this requires a dedicated pre-amplifier, it has become the standard approach today. Recently, there have been attempts to avoid the bifurcation regime by saturating the gain medium and thus operating the amplifier beyond the unstable region [9]. Since this may result in extremely high transient pulse energies, it requires a careful design of the RA to avoid pulses above the damage threshold of the optics in the cavity. Alternatively, operation at a specifically chosen limit cycle and picking the desired pulses from the amplified pulse train afterwards is suggested in [10].
By increasing the seed energy in a saturating amplifier as proposed in [8], one effectively reduces the total gain, which is clearly beneficial for the amplifier’s stability properties as it is also known from electronic amplifiers. However, circuit engineers developed another highly successful and far-reaching strategy to cope with emerging instabilities in high-gain amplifiers: feedback. In fact, [4] already points to the possibility of active stabilization techniques. Since many amplifiers already utilize controllable input filters to compensate for effects such as gain narrowing or nonlinear phase distortions [11–15], it seems natural to investigate strategies to stabilize the amplifier’s dynamics by modulating the energy of the seed pulses. While simple feedback schemes are commonly employed in CW lasers or Q-switched lasers to statically compensate for external disturbances, the authors are not aware of any applications of feedback to actively modify the dynamics of RAs.
In this paper, we present active feedback methods to stabilize otherwise unstable operational regimes of RAs and thus suppress bifurcations of the uncontrolled amplifier. To this end, we will start with a general space-dependent mathematical model of regenerative amplifiers and show that a detailed description of the amplifier’s pulse-to-pulse dynamics for different levels of (continuous) pumping can be systematically derived based on mild assumptions only. The resulting discrete-time dynamic model can be considered a generalization of the model used in [7]. Using this description, two feedback approaches will be developed that utilize measurements of the output pulse energy or the residual pumping beam, respectively. Finally, the validity of the proposed feedback strategy is highlighted by experimental results.
2. Dynamics of regenerative amplifiers
In case the spectral properties of the pulse amplification process can be neglected, instead of using the rather complex Maxwell-Bloch equations, the laser light is usually treated as a transport process of monochromatic photons conveying no additional phase information described by the intensity $I$ known as Frantz-Nodvik model [16]. For simplicity, we will focus on the case of a RA using a ring cavity design as shown in Fig. 1, although the proposed approach can be directly transferred to linear cavities. Thus, identical seed pulses $I_{\textrm{seed}}(t)$ from a source with repetition rate $f_{\textrm{rep}}$ pass some kind of modulation device that allows to adjust the amount of additional spectrally homogeneous damping for each individual seed pulse. Using an acousto-optic modulator (AOM) with an adjustable loss factor $\eta _{\textrm{AOM}}^{n}$, the resulting input pulse $I_{\textrm{in}}^{n}(t)$ is given by
It is assumed that further injection losses are already included in $I_{\textrm{seed}}(t)$ and thus $0 \leq \eta _{\textrm{AOM}}^{n} \leq 1$. The resulting pulse $I_{\textrm{in}}^{n}(t)$ is then injected into the resonator cavity with concentrated losses $\eta _{\textrm{RC}}$, whereby the $k$-th round-trip of the $n$-th input and output pulse is denoted by $I_{\textrm{in}}^{n,k}(t)$ and $I_{\textrm{out}}^{n,k}(t)$, respectively. The final amplified pulse after $N_{\textrm{RT}}$ round-trips is then released by switching the Pockel’s cell (PC) and all losses due to the pulse extraction are again subsumed in $\eta _{\textrm{out}}$, i.e.,To describe the amplification and regeneration process, we assume that the gain medium is operated in a so-called end-pumped configuration, i.e. that the continuous pumping beam is fed through the crystal axially aligned with the laser beam and in the same direction of propagation. Additionally, we assume that the transversal intensity profile of the pulses and the pumping beam is homogeneous over a diameter $d$ and zero outside. Aligning the propagation direction with the $z$-axis and assuming a homogeneously broadened quasi-three-level medium with the population inversion $\Delta N = N_2 - N_1$ of the lasing transition, the interaction of the pulse $I$ and the pumping beam $I_{\textrm{P}}$ during each round trip can be described by [17]
2.1 Amplification during a single pass
For a single pass of the laser pulse through the laser material, the influence of the comparatively weak pumping beam and the relaxation of the excited populations via $\gamma _{21}$ can be neglected due to the extremely short time scales. Applying these assumptions to (3) and transforming the equations onto a time frame $ t \mapsto t - {z \mathord{\left/ {\vphantom {z v}} \right. } v} $ travelling with the pulse, one obtains
2.2 Gain restoration due to pumping
Now that a simple expression for the pulse amplification and the resulting saturation of the available gain has been found, a compatible description of the gain restoration due to the pumping process of duration $ \Delta t = {1 \mathord{\left/ {\vphantom {1 {{f_{\textrm{rep}}}}}} \right. } {{f_{\textrm{rep}}}}} $ shall be derived. Since no laser pulse is present during the regeneration periods (i.e., $I(z,t) = 0$), the full set of Eqs. (3) can be rewritten as
2.3 Discrete-time dynamical model
Using the solution of a single pass through the laser medium according to (7) and considering the losses in the cavity, one can write the effect of a single round-trip according to Fig. 1 as
Apart from short transient operations during start-up or when switching to a different operating point, RAs are mainly run in steady-state operation. For a given constant input $u^s$ and a stationary pump power $p^s$, the resulting steady state $x^s$ is given by the nonlinear equation
with the corresponding steady-state output Because of the one-dimensional state variable, the solution of the fixed-point Eq. (16) as well as the dynamics of the RA can be illustrated graphically as shown in Fig. 2 for different pump powers and round-trip numbers using the parameter values given in Table 1 and a 10 nJ seed pulse. As one would expect, the qualitative behavior of the dynamic map $f$ is similar to the simpler model with neglected cavity losses in [7]. However, the neglected cavity losses influence the observed dynamics quite significantly for high values of $x$ and are particularly relevant for the output map $h$, where the energy of the output pulses starts to decrease for high values of $x$ as the net gain during the last round-trips in a highly saturated gain medium becomes less than one. As we will see in Section 3, this entails limitations on the control strategy. As noted above, the maxima of the steady-state output power $P_{\textrm{out}}^s = \frac {\hbar \omega A_B}{\sigma (f_1 + f_2)} f_{\textrm{rep}} y^s$ and the steady-state efficiency $ {\eta _{\textrm{eff}}} = {{P_{\textrm{out}}^s} \mathord{\left/ {\vphantom {{P_{\textrm{out}}^s} {{P_{\textrm{pump}}}}}} \right. } {{P_{\textrm{pump}}}}} $ are usually inside or at the border of the unstable region, i.e., all $x^s$ with $|{\frac {\partial f}{\partial x}}(x^s,u^s) | \geq 1$, as illustrated in Fig. 3.3. Feedback stabilization of unstable operating regions
Building on the derived model above, we want to stabilize the dynamics of the RA by modulating the seed pulses based on some measured quantity of the amplifier. While using measurements of the output pulses $y^n$ is probably the most obvious quantity, a simpler alternative is to use the transmitted pump light instead.
Since the pump power is usually not changed dynamically, but only adjusted to reach a desired operating point, we will utilize a so-called gain scheduling approach [18] in the following. Thereby, one determines a family of feedback laws for an arbitrary pump power $p^s$ and switches between them during operation according to the current pump power. To simplify the design of a suitable feedback law, we linearize the system dynamics in the vicinity of an arbitrary operating point $(x^s(p^s),u^s)$, which is a solution of (16), to obtain a family of linear systems
with $\Delta u^n = u^n - u^s$, $\Delta x^n = x^n - x^s(p^s)$, and $\Delta y^n = y^n - y^s(p^s)$ as well as3.1 Feedback stabilization using pump light measurements
Like incoming laser pulses are amplified due to the population inversion in the laser medium, so is the light of the pumping beam absorbed. As a result, one can see fluctuations in the transmitted pump light due to the extraction and regeneration of the population inversion. Solving (9a) by separation of variables yields
3.2 Feedback stabilization using output energy measurements
While approaching the stabilization task by measuring the transmitted pump light is a simple and elegant method, such measurements are not always available, as for broadband pumping sources such as flash lamps or arc lamps. Additionally, if the energy of the amplified pulses are significantly below the capabilities of the gain medium, the fluctuation in the transmitted pumping beam may be too small to be measured reliably. In such cases, one can still follow the straightforward way to stabilize the regenerative amplifier using measurements of the pulse energy. Doing so will require to either infer the current state from output measurements by using state observers [19] or by applying some kind of output feedback. In the following, we will proceed along the latter path as it comes up with an easily implementable and computationally inexpensive control law by transforming (14) into a description based on input and output quantities only. In the vicinity of any $x$ with ${\frac {\partial h}{\partial x}}(x,u) \neq 0$ there exists an inverse map $h^{-1}(y,u)$ such that $h(h^{-1}(y,u),u) = y$ for all $u$. Using this inverse map, the state dynamics (14) can be rewritten as
However, one is still able to pass through such singular points by limiting the resulting feedback gains in its vicinity and defining a region around $\tilde {p}$ within which one interpolates between the values at the boundary in some sense. Assigning the nearest boundary values and using an LQR formulation as pointed out above, one obtains a feedback gain ${\textbf{k}}_{\textrm{FB}}^T(p)$ illustrated in Fig. 5. Writing the control law with original variables, the output feedback law reads as
3.3 Simulation scenarios
The performance of the two feedback laws presented above shall be illustrated with simulation scenarios using the full mathematical model (3). To cover most of the interesting operational regimes of a regenerative amplifier, the optical pump power $P_{\textrm{pump}}(t)$ (and thus $p^n$) is chosen to start from a weakly pumped state at 12.5 W and increases linearly up to 87.5 W during the first 80 ms. After having passed the unstable regime completely, the amplifier is held at this highly pumped state for 40 ms, after which the pump power is decreased rapidly to 32.5 W within 15 ms and held there for the remaining time.
Figure 6 shows simulation scenarios for both feedback laws in comparison to the uncontrolled open-loop case. In the uncontrolled case, small perturbations build up shortly after the amplifier enters the unstable regime while both feedback laws are able to successfully stabilize the amplifier and follow the desired steady-state output energies $W^s_{\textrm{out}}$. The comparatively slow increase in the pump power at the beginning elicits only minor corrective action by both feedback laws. However, the rapid decrease in pump power is fast enough to operate the output feedback law (30) close to its limits as indicated by the large input corrections necessary to keep the deviations of the output energy $W_{\textrm{out}}^n$ close to the desired values. The availability of direct state measurements renders the feedback law (22) more robust to errors introduced by the gain scheduling approach. However, notice that the output feedback law (30) manages to keep the amplifier stable while passing through the singular pumping power $\tilde {p}$ where the amplifier produces the highest possible output energies.
4. Experimental results
To demonstrate the effectiveness of the proposed feedback strategy in practice, it is applied to a Yb:CaF2-based regenerative amplifier utilizing a linear cavity design with a double-pass pumping scheme which is seeded by 4 nJ pulses at $f_{\textrm{rep}} = 1$ kHz. By choosing $N_{\textrm{RT}} = 102$ (i.e., each pulse passes the gain medium 204 times due to the linear cavity), one obtains an amplifier that allows sustained unstable operation at low power outputs without damaging the cavity optics by transient pulses. Due to the low energy extraction at such operating points, the fluctuations of the transmitted pump light are too small to be measured reliably. As a result, the feedback stabilization using pulse energy measurements presented in Section 3.2 is the only viable option.
The energy of the individual output pulses $W_{\textrm{out}}^n$ is measured using a photo diode and observing the decay of the resulting voltage $v_{\textrm{ph}}$ due to the parasitic resistances and capacitors. Since the photo current produced by the diode is proportional to the optical intensity, the jump height of the measured voltage signal is directly proportional to the energy of the incident laser pulse. Online estimates of the jump height can be obtained using polynomial or exponential least-squares approximations. In particular, least-squares techniques building on a discrete-time formulation of the underlying RC network compared favorable in terms of reliability and computational effort. As illustrated in Fig. 7, simulation results of the mathematical model (14) calibrated to the given setup (see, e.g., [20] for the relevant parameters of Yb:CaF2 crystals) agree very well with measurements of $W_{\textrm{out}}$. In particular, the model accurately predicts the energies of the $\omega$-limit set $W_{\textrm{out}}^\omega$ for unstable operating points above $P_{\textrm{pump}} = 50\;\textrm{W}$.
Since the gain-scheduled feedback law (30) is only a scalar dynamic system, it can be easily implemented on existing signal processing units. Notice that ${\textbf{k}}_{\textrm{FB}}(p)$ and $y^s(p)$ can be calculated in advance for all $p$ from the given system model and do not need to be computed online. Choosing $u^s = 0.7 {\mkern 1.5mu\overline {\mkern -1.5mu u \mkern -1.5mu}\mkern 1.5mu}$ and applying the feedback law to the experimental setup, one is able to stabilize the regenerative amplifier and steer it into otherwise unstable operational regions. Unlike in the open-loop case, the feedback law suppresses the onset of bifurcations and keeps the output pulses $W_{\textrm{out}}^n$ close to the desired steady-state energy $W_{\textrm{out}}^s$ for a given pump power $P_\textrm{pump}$ (see Fig. 7).
To illustrate the temporal effects of the feedback, we deactivate the feedback law by setting $u^n = u^s$ (i.e., $\eta _\textrm{AOM}^n = 0.7$) while operating the amplifier at $P_\textrm{pump} = 52.36\;\textrm{W}$. Small perturbations start to build up immediately as shown in Fig. 8 (top) and the amplifier approaches the limit cycle $W_\textrm{out}^\omega$ again (cf. to Fig. 7). Details of the corresponding measured photo diode voltage $v_\textrm{ph}(t)$ are shown at the bottom of Fig. 8 for the case with and without feedback, respectively.
5. Conclusions and outlook
This paper demonstrates that unstable operational regimes and the resulting bifurcations of regenerative amplifiers can be suppressed by using feedback methods to actively modify the energy of incoming seed pulses. Starting from a general space-dependent mathematical description of the amplification and regeneration processes, it was shown that the dynamic behavior of RAs can be captured by a scalar discrete-time dynamic system under very mild assumptions only. Based on this dynamic model, two different control strategies were proposed using either output energy or transmitted pump light measurements. Both approaches are able to successfully stabilize RAs for general operating points. As a result, the utilization of feedback methods allows the operation of high-gain optical amplifiers without requiring multiple amplification stages to avoid bifurcations.
Using measurements of the output pulses entails limitations close to the operating point with maximum output energies as shown in Section 3.2. Nevertheless, it is possible to steer the amplifier through this singular point by traversing its vicinity fast enough while truncating the feedback gains. Those limitations do not apply to feedback approaches based on pump light measurements as presented in Section 3.1. Simulation studies furthermore show that such schemes are more robust due to the direct relation of the measurement to the system’s state. However, the fluctuations of the transmitted pumping beam might be too small to obtain reliable measurements if the energy extracted by the amplification process is small compared to the stored energy. The choice between both feedback schemes thus strongly depends on the specific properties of the given amplifier. To demonstrate the practical effectiveness of active feedback approaches, the output feedback scheme was successfully applied to a Yb:CaF2-based RA.
Apart from stabilization purposes, altering the natural dynamics of RAs can be beneficial for the suppression of disturbances and to ensure a fast decay of perturbations to the steady state. Additionally, one could possibly apply similar feedback techniques to instabilities arising from the cavity dynamics in mode-locked cavity-dumped lasers [21]. While both feedback schemes presented in this thesis can be easily implemented on existing signal processing units, they rely on a mathematical model that is calibrated for the specific setup. Since the open-loop operation of unstable amplifiers close to their damage threshold is quite problematic, the development of closed-loop identification methods that can be performed online seems a promising direction for future investigations.
Appendix
Introducing the Jacobian of ${\textbf{f}}_{\textrm{SP}}$ with respect to $x$ and $u$ as
Funding
Technische Universität Wien Bibliothek; Österreichische Forschungsförderungsgesellschaft (859211).
Acknowledgments
The authors acknowledge the TU Wien University Library for financial support through its Open Access Funding Program.
Disclosures
The authors declare no conflicts of interest.
References
1. F. Krausz and M. Ivanov, “Attosecond physics,” Rev. Mod. Phys. 81(1), 163–234 (2009). [CrossRef]
2. D. Goswami, “Optical pulse shaping approaches to coherent control,” Phys. Rep. 374(6), 385–481 (2003). [CrossRef]
3. S. Nolte, C. Momma, H. Jacobs, A. Tünnermann, B. N. Chichkov, B. Wellegehausen, and H. Welling, “Ablation of metals by ultrashort laser pulses,” J. Opt. Soc. Am. B 14(10), 2716–2722 (1997). [CrossRef]
4. J. Dörring, A. Killi, U. Morgner, A. Lang, M. Lederer, and D. Kopf, “Period doubling and deterministic chaos in continuously pumped regenerative amplifiers,” Opt. Express 12(8), 1759–1768 (2004). [CrossRef]
5. J. G. Mance, M. N. Slipchenko, and S. Roy, “Regenerative amplification and bifurcations in a burst-mode Nd:YAG laser,” Opt. Lett. 40(21), 5093–5096 (2015). [CrossRef]
6. D. Mueller, A. Giesen, and H. Huegel, “Picosecond thin disk regenerative amplifier,” in 14th International Symposium on Gas Flow, Chemical Lasers, and High-Power Lasers, (Wroclaw, 2003), pp. 281–287.
7. M. Grishin, V. Gulbinas, and A. Michailovas, “Dynamics of high repetition rate regenerative amplifiers,” Opt. Express 15(15), 9434–9443 (2007). [CrossRef]
8. M. Grishin, V. Gulbinas, and A. Michailovas, “Bifurcation suppression for stability improvement in Nd:YVO4 regenerative amplifier,” Opt. Express 17(18), 15700–15708 (2009). [CrossRef]
9. P. Kroetz, A. Ruehl, G. Chatterjee, A.-L. Calendron, K. Murari, H. Cankaya, P. Li, F. X. Kärtner, I. Hartl, and R. J. D. Miller, “Overcoming bifurcation instability in high-repetition-rate Ho:YLF regenerative amplifiers,” Opt. Lett. 40(23), 5427–5430 (2015). [CrossRef]
10. L. von Grafenstein, M. Bock, G. Steinmeyer, U. Griebner, and T. Elsaesser, “Taming chaos: 16 mJ picosecond Ho:YLF regenerative amplifier with 0.7 kHz repetition rate,” Laser Photonics Rev. 10(1), 123–130 (2016). [CrossRef]
11. A. Efimov, M. D. Moores, N. M. Beach, J. L. Krause, and D. H. Reitze, “Adaptive control of pulse phase in a chirped-pulse amplifier,” Opt. Lett. 23(24), 1915–1917 (1998). [CrossRef]
12. J. Prawiharjo, N. K. Daga, R. Geng, J. H. V. Price, D. C. Hanna, D. J. Richardson, and D. P. Shepherd, “High fidelity femtosecond pulses from an ultrafast fiber laser system via adaptive amplitude and phase pre-shaping,” Opt. Express 16(19), 15074–15089 (2008). [CrossRef]
13. P. Malevich, G. Andriukaitis, T. Flöry, A. J. Verhoef, A. Fernández, S. Ališauskas, A. Pugžlys, A. Baltuška, L. H. Tan, C. F. Chua, and P. B. Phua, “High energy and average power femtosecond laser for driving mid-infrared optical parametric amplifiers,” Opt. Lett. 38(15), 2746–2749 (2013). [CrossRef]
14. Y. Chiba, H. Takada, K. Torizuka, and K. Misawa, “65-fs Yb-doped fiber laser system with gain-narrowing compensation,” Opt. Express 23(5), 6809 (2015). [CrossRef]
15. A. Deutschmann, P. Malevich, A. Baltuška, and A. Kugi, “Modeling and Iterative Pulse-Shape Control of Optical Chirped Pulse Amplifiers,” Automatica 98, 150–158 (2018). [CrossRef]
16. L. M. Frantz and J. S. Nodvik, “Theory of Pulse Propagation in a Laser Amplifier,” J. Appl. Phys. 34(8), 2346–2349 (1963). [CrossRef]
17. A. E. Siegman, Lasers (University Science Books, Mill Valley, 1986).
18. D. J. Leith and W. E. Leithead, “Survey of gain-scheduling analysis and design,” Int. J. Control 73(11), 1001–1025 (2000). [CrossRef]
19. T. Kailath, Linear Systems, Prentice-Hall Information and System Sciences Series (Prentice Hall, Upper Saddle River, NJ, 1980).
20. F. Druon, S. Ricaud, D. N. Papadopoulos, A. Pellegrina, P. Camy, J. L. Doualan, R. Moncorgé, A. Courjaud, E. Mottay, and P. Georges, “On Yb:CaF2 and Yb:SrF2: Review of spectroscopic and thermal properties and their impact on femtosecond and high power laser performance,” Opt. Mater. Express 1(3), 489–502 (2011). [CrossRef]
21. M. S. Pshenichnikov, W. P. de Boeij, and D. A. Wiersma, “Generation of 13-fs, 5-MW pulses from a cavity-dumped Ti:sapphire laser,” Opt. Lett. 19(8), 572–574 (1994). [CrossRef]