Laser is a fascinating technology with a wide range of applications. There are many different physical effects associated with lasers, for example properties like high intensity beam power in case of ultra-short pulses. These effects make both studying of laser physics and optimization of lasers a difficult task. Therefore, a laser simulation program is a helpful tool to look inside a laser and analyze the properties of a specific laser configuration. For solid state lasers, we developed laser simulation program ASLD, which meets the requirements for both educational and industrial applications. It also allows us to model and simulate different physical effects, which are important for an accurate laser simulation. These physical effects are (see ):
- thermal lens effect including a detailed thermal and mechanical analysis of a laser crystal
- gain and dynamic behavior of the laser light, including spectral and temperature dependence of all physical parameters
- stability of the resonator cavity and influence of optical elements on beam radius
- birefringence inside a laser crystal and impact on depolarization of the beam and losses
- spectral analysis including stretching and compressing of a short laser pulse
- analysis of the pump light configuration
Several simulation techniques have to be applied in order to simulate these different physical effects. Furthermore, a suitable software concept is required to couple these simulation techniques in a software package. ASLD is designed and implemented using a software concept designed for this specific application. Each software module has its input and output parameters, which can be specified and modified by the user. A change in these parameters causes only affected modules to be recalculated. Figure … depicts the main modules of ASLD and their relations. Let explain simulation techniques of these modules in more detail.
Modelling of pump light
Depending on properties of the pump light, different techniques have to be applied to simulate the pump light. Diode light is usually a multimode light. Then, ray tracing is an optimal simulation technique to calculate the absorbed pump light distribution. However, if the pump light is a laser beam, super Gaussian functions are more appropriate for describing pump light propagation inside a laser crystal. It is important is to take into account the crystal’s doping concentration. By bonding crystals with different doping concentrations, we can obtain a crystal with different doping concentrations at its different sections (see ). Figure 1-a) shows the simulated pump light absorption for such an end- pumped laser rod.
Finite element analysis of temperature and deformation
Finite element is a well-known method for the calculation of temperature and deformation caused by pumping a crystal. The thermal boundary conditions have to be specified accurately. Figure 1-b) shows the temperature distribution calculated for an end-pumped laser rod, where the crystal was cooled only at an interior section. The resulting “von Mises” stress is illustrated in Figure 2-a). Using the stress inside the crystal, one can calculate the stress induced birefringence, which influences the depolarization of beam (see Figure 3).
Gauss mode analysis
Gauss mode analysis applies ABCD matrices for each optical element inside a laser resonator. The ABCD matrix of the laser crystal is obtained by a product of several ABCD matrices corresponding to end faces of the crystal and a sequence of Gaussian ducts inside the crystal. To this end, the deformation of the end faces and the refractive index inside the crystal has to be fitted by parabolic functions. The results of a Gauss mode analysis are stability information of the resonator and the beam radius inside cavity for stable resonators. Beam radius is needed to calculate overlap integrals for the dynamic mode analysis.
Dynamic Mode Analysis DMA
The dynamic mode analysis is a highly efficient simulation technique that calculates output power and beam quality of a laser resonator. It applies a finite volume discretization for the population inversion (see ). Figure 2-b) depicts the calculated 3-dimensional population inversion obtained for the pump configuration in Figure 1 a). Using these 3-dimensional simulation results, one can modify pump configuration and beam radius to obtain an optimal overlap of both. By this approach, one can optimize output power and beam quality of the laser resonator. The numerical algorithm of DMA solves a system of rate equations for population inversion and high and low orders modes. This system is solved by a highly efficient adaptive time stepping method. The dynamic mode analysis can be applied for 3-level and 4-level laser systems. Parameters like stimulated emission cross section and reabsorption cross section can be modeled temperature-, wavelength- and space-dependent, since population inversion is discretized on a finite volume mesh. This leads to more accurate simulation results, than taking a spatial average of these parameters. ASLD applies dynamic mode analysis for passive and active Q-switches, passively Q-switched SESAM, intra-cavity frequency doubling and regenerative amplifiers.
Amplifier simulation using Gaussian functions
There are different simulation techniques to calculate the gain of a laser amplifier. Similar to resonators it is important to consider 3-dimensional behavior of the pumping configuration and population inversion. Then, one can either apply a beam propagation method or restrict himself to the gain of high and low order Gaussian functions. This second approach is less computationally intensive and gives clear information how beam quality changes through the amplification of laser beam. The input beam can be chosen to be the output of a laser resonator. There are two approaches to include polarization effects. One is to calculate the depolarization of single rays through the amplifier and combine this with the amplification of each Gaussian functions. The other approach is the vectorial beam propagation method, which is more accurate but also more computationally intensive.
Vectorial Beam Propagation Method
The vectorial beam propagation method applies a paraxial approximation of the vector Helmholtz equation obtained by Maxwell’s equations (see ). It is discretizes the electrical field vector (Ex,Ey) in xy-plane and propagates it step by step in propagation direction (z-direction). Using this approach, the birefringence effects, like in Figure 3b) can be accurately calculated (see ).
Simulation of ultrashort pulses
ASLD can simulate chirped pulse amplification of short pulses. Due to the high energy of the pulse, it has to be stretched before amplification in order to avoid damage to the laser system. The pulse is re-compressed after amplification. Both stretcher and compressors consist of mirrors and gratings. Simulations take into account that the spectrum of the pulse might be narrowed by limited gain bandwidth (see ).
ASLD is a multi-functional software tool and can be used for educational purpose, research institutes as well as industries. It assists developers and engineers to build lasers with less number of expensive try and error experiments. It is also helpful for students who want to analyze and understand various effects associated with laser beam generation and propagation inside a cavity and amplifier.
Graupeter, T., Hartmann, R., Pflaum, C. “Calculations of Eigenpolarization in Nd:YAG Laser Rods Due to Thermally Induced Birefringence”. IEEE Journal of Quantum Electronics, 50 (12), 1035 - 1043, (2014).Google Scholar
Hartmann R., Pflaum C., Graupeter T., “Analysis of thermal depolarization compensation using full vectorial beam propagation method in laser amplifiers”, Proc. SPIE 9343, (2015).Google Scholar
Huabao Cao, Xingqiang Lu, Dianyuan Fan, “Numerical simulation of gain narrowing control by hybrid amplifiers chain based on Ti:sapphire and Ti:chrysoberyl”, Optics Communications, 284 (6), 1622–1625, (2011).Google Scholar
Koechner, W., [Solid-state laser engineering], Springer, Berlin, (2006).Google Scholar
Rahimi, Z., Pflaum, C., “Simulation of solid-state lasers with composite and ceramic crystals”, Proc. SPIE 8599, (2013).Google Scholar
Wohlmuth, Matthias; Pflaum, Christoph; Hahn, C., “Dynamic multimode analysis of Q-switched solid state laser cavities”, Optics Express, 17, 17303–17316 (2009).Google Scholar