13 April 2018 Modeling of the laser beam shape for high-power applications
Author Affiliations +
Aperture losses and thermo-optic effects (TOE) inside optics as well as the effective beam width in far field should be taken into account in the analysis of the most appropriate laser beam profile for high-power applications. We have theoretically analyzed such a problem for a group of super-Gaussian beams taking first only diffraction limitations. Furthermore, we have investigated TOE on far-field parameters of such beams to determine the influence of absorption in optical elements on beam quality degradation. The best compromise gives the super-Gaussian profile of index p  =  5, for which beam quality does not decrease noticeably and the thermo-optic higher order aberrations are compensated. The simplified formulas were derived for beam quality metrics (parameter M2 and Strehl ratio), which enable estimation of the influence of heat deposited in optics on degradation of beam quality. The method of dynamic compensation of such effect was proposed.
Jabczyński, Kaskow, Gorajek, Kopczyński, and Zendzian: Modeling of the laser beam shape for high-power applications



One of the specific tasks of laser engineering is the transformation of high-power laser beam satisfying special requirements on beam waist size in far field and additionally the restrictions on beam forming optics (sizes, aperture losses, thermal effects, etc.). From the practical point of view, the main limitations in near-field (i.e., inside optics) are posed by thermal effects especially for multi kWatt, complex, multielement optical systems. The “rule of thumb” in engineering practice is that optics should have aperture of about 2 times larger than the laser beam diameter defined for Gaussian profile at 1/e2 level, which is not a good choice for mitigation of thermal effects by the way. Thus, the question on the optimal laser beam profile which mitigates thermal effects in near-field optics not degrading brightness and beam width in far field is still open. Let us note that the laser beam quality has to be determined according to D86.5 definitions (see Ref. 1), Strehl ratio and second moment definitions are not appropriate measures for such a case. Therefore, beam parameter product (BPP)=πw86.5θ86.5/λ defined on the basis of D86.5 will be used here as a metrics.

To determine the best laser beam profile, the minimization of aperture losses and thermo-optic effects (TOEs) in near-field as well as the minimal effective width of laser beam in far field should be taken into account. Let us notice that the best from the thermo-optic limitations “top-hat” profile results in far field in multilobe, “sombrero”-like shape, which has wider D86.5 diameter than equivalent Gaussian one. On the other hand, for Gaussian profile of given 1/e2 diameter in the near-field, the clear aperture should be wider (2 times), which significantly increases TOEs. Thus, it is an evident trade-off between the above requirements in near and far fields. Moreover, we have to take into account the manufacturability of beam shaper destined to transform Gaussian beam to the most appropriate shape.

To start the analysis, we have to choose the appropriate basis of beam profiles (see Refs. 23.4.5). The most convenient for our purposes seems to be the super-Gaussian beams (SGB)2,3 and flattened Gaussian beams (FGB).4,5 In both cases, the analytical method for M2 parameter calculations was found,3,4 whereas for BPP calculations the numerical approach is required. Both approaches (SGB and FGB) should lead to the same results in principle.5 We have taken for our purposes the SGBs because of the simpler mathematical description important for numerical simulations and wide diversity of shapes especially for close to Gaussians profiles. Moreover, the technical realization of diffraction limited beam shapers transforming Gaussian beam into SGB is well known (see e.g., Refs. 6 and 7).

In Sec. 2, we have analyzed such a problem for SGBs taking only diffraction effects into account. In Sec. 3, we have analyzed TOEs for such cases. The typical laser system consists of several elements, including mirrors and refractive elements, and each experiences the TOEs as a result of residual absorption on surfaces and in the volume. We have modeled the influence of TOEs on the far-field parameters of SGB beams applying COMSOL software8 for two cases: volume heat density source and surface heat source for idealized heat contacts. The results of numerical analysis undertaken for given heat power can be rescaled to other cases applying the simple relations given at the end of paper.


Diffraction Analysis for Super-Gaussian Beams

To analyze the parameters of the laser beam of a different profile in far field, we have chosen a family of SGB defined in near-field as follows:


where SGp is the amplitude function of SGB of index p, Np(wp)=πwp2Γ(1/p)21/pp1 norm of SGp accomplishing the unit power of each beam, Γ(x) is the gamma function, and wp is the radius of beam.

Furthermore, we calculated the amplitude functions SGff,p in far field applying Hankel transform valid for cylindrical symmetry (see details in Ref. 6) as follows:


where cp is the constant and f is the focal length.

The intensity profiles in near-field for p=1, 2, 5, and 32 are shown in Fig. 1 and corresponding intensity distributions in far field are shown in Fig. 2. For each SGp, the different beam radii wp were taken in such a way that it contains 99.95% of SGp power in the same aperture Daper=2Waper (which corresponds to Waper=2w1 for Gaussian beam p=1).

Fig. 1

Intensity versus radius of SG1, SG2, SG5, and SG32 beams in near-field having the same aperture losses of 0.05% at radius=2.


Fig. 2

Intensity distributions of SG1, SG2, SG5, and SG32 beams in far field.


To define beam diameter in far field, power in bucket (PIB) distributions (see Ref. 1) were calculated as follows:


where Iff,p(rff)=|SGff,p(rff)|2 is the intensity profile of SGp in far field and rff,max is the radius in far field, beyond which the intensity is negligible for all p (see Figs. 2 and 3). Knowing PIB distributions, we have determined the beam radius Rff,86.5,p according to criterion of 86.5% of PIB.

Fig. 3

PIB functions for intensity distributions of SG1, SG2, SG5, and SG32 beams in far field.


Let us notice (see Fig. 1) that for p=1, we have Gaussian beam, whereas for p=32 we have nearly top-hat profile. The criterion of 86.5% of PIB, frequently used in laser engineering (see e.g., Refs. 1, 9, and 10), corresponds to classical 1/e2 diameter definition of Gaussian beam. Strictly according to that criterion SG32 is better, but it has multilobe shape with wide pedestals in far field (see Figs. 2 and 3); moreover, technical realization of cost-effective beam shapers for such profile is problematic. Applying higher level of criterion, e.g., 95%, the lowest diameter has the SG5 beam and similar smooth profile in far field as Gaussian with <5% in pedestal (see Fig. 3).

We can now determine the optical measures being the merit functions in analysis: relative brightness Bp, effective beam parameter product BPPp, and Strehl ratio SRp as follows:






where Anf,p=πRnf,86.5,p2 and Aff,p=πRff,86.5,p2 are the areas of SGp beam in near and far field, respectively, Rnf,86.5,p and Rff,86.5,p are the beam radii of SGp beam in near and far field determined according to 86.5% criterion, respectively.

To assess the flatness of SGp profile in near-field, we have introduced additional merit function MFnf,p as a ratio of maximal intensity of SG1 and SGp as follows:


It is easy to show that for perfect top-hat beam (p) MFnf,=8. In Table 1, the results of calculations for a few SGp beams were collected.

Table 1

Results of calculations of merit functions for p=1, 2, 5, and 32.


Let us notice that taking into account only the maximum of brightness (or equivalent minimum of BPP), the best is Gaussian beam (p=1) evidently. For the almost top-hat (p=32) beam of equivalent aperture, 43% drop in brightness and 33% increase in BPP were found. However, if we additionally have to consider the profile flatness in near-field, which affects temperature profile and thermally induced distortions, the answer becomes more complicated. Taking both near- and far-field requirements, the best compromise seems to be SG5 (p=5) beam (blue curves in Figs. 13), which has smooth profile in far field, slightly increased BPP (BPP5=1.259), but much more flattened profile in near-field (MFnf,5=4.77). Moreover, as we will show in Sec. 3, such profile enables nearly compensation of higher order TOEs.


Modeling of Thermal Optics Effects for SG Beams

To investigate the influence of TOEs on beam quality, we have taken the typical optical element of 50-mm diameter, 10-mm thickness made of fused silica. The “idealized” heat contacts (constant temperature on side, and negligible heat transfer to rear and front facet) were assumed. The eight different heat sources (for each SGp, pair of volume or surface source) were considered.

To compare the effects of such eight different heat sources, we have assumed that the same heat power of 1 W (corresponding to 10-ppm heat conversion for incident 100 kW of laser power) was deposited in optical element. Let us notice that in the best quality fused silica glass, the absorption coefficient is much lower than 105  1/cm.1112.13 The surface absorption is determined by quality of surface itself and the absorption in dielectric coatings. It is possible to achieve the same level of absorption losses in the highest quality mirrors.1112.13 However, as a result of technical imperfections and superposition of different factors, the realistic value of absorption is estimated of a few dozens of ppm.

For each case, the problem was solved applying COMSOL software,8 and three-dimensional (3-D) maps of temperature increase were calculated [see examples for SG1 and SG32 beams in Figs. 4(a)4(d)].

Fig. 4

(a) 3-D map of temperature increase, SG1_vol is the volume absorption for Gaussian beam SG1 with 2w1=25  mm, heat power 1 W, ΔT=4.5  K. (b) 3-D map of temperature increase, SG1_surf is the surface absorption for Gaussian beam SG1 with 2w1=25  mm, heat power 1 W, ΔT=14  K. (c) 3-D map of temperature increase, SG32_surf is the surface absorption for “near top-hat” beam SG32 2w32=47  mm, heat power 1 W, ΔT=0.2  K. (d) 3-D map of temperature increase, SG32_surf is the surface absorption for near top-hat beam SG32 with 2w1=47  mm, heat power 1 W, ΔT=3.5  K


Next, we calculated the profile of averaged temperature ΔTavg,i versus radius ri integrating temperature profile over z-depth. Then, multiplying by effective thermal dispersion coefficient χT and b—thickness of sample, the thermally induced optical path difference (OPD) was determined


where the effective thermal dispersion for transmission is defined as



The OPD can be divided into paraxial component OPDpar and residual nonparaxial OPDnp according to the following equation:


where MT is the paraxial thermal optical power.

To determine OPD dependence on radius, the mean square approximation of data array (ri,ΔTavg,i) calculated in COMSOL was applied as follows:



After a few simple transformations


where a1 is the quadratic coefficient of power series of mean square approximation of data array (ri,ΔTavg,i).

Residual, nonparaxial part of OPDnp corresponds to higher order thermally induced distortions [see Figs. 5(a) and 5(b)] resulting in beam quality degradation.

Fig. 5

(a) Nonparaxial OPDnp versus radius for four cases of volume heat source, heat power 1 W, SG1 - dotted, SG2 - dashed, SG5 -dash dot, SG32 - continuous. (b) Nonparaxial OPDnp versus radius for four cases of surface heat source, heat power 1 W, SG1 - dotted, SG2 - dashed, SG5 -dash dot, SG32 - continuous.


The surface heat sources [Fig. 5(b)] result in larger OPDs comparing to volume absorption, which agrees well with intuition and engineering practice. The magnitude of near top-hat OPDnp [SG32 black curves in Figs. 5(a) and 5(b)] is much smaller and has the opposite sign regarding to OPDnp calculated for Gaussian beam. Therefore, we can suppose that the proper choice of beam profile can result in near compensation of nonparaxial OPDs at least. It was shown in Figs. 5(a) and 5(b) (green curves) that for SG5 beam such effect exists, which gives another argument for its beneficial properties in a case of high-power applications. Moreover, having beam shaper transforming SG1SG5 and playing with reflecting and refractive elements, we can achieve dynamic compensation of nonparaxial OPDs of all system for variable incident laser power.

To determine the impact of TOEs on laser beam metrics, we have to calculate in first step variance σ=rms(OPDnp) taking into account weighting functions corresponding to given beam profile. Furthermore, we calculated Strehl ratio SRapr and Mapr2 parameter applying the following approximated equation:1


where σ=rms(OPDnp) of nonparaxial OPDnp calculated with appropriate weighting function.

The results of calculations for 1-W heat power (10 ppm for 100-kW incident power) were collected in Table 2. The calculations made for lower power of 0.5 W showed the similar dependencies, only magnitudes have changed proportionally to absorbed power. Thus, we can conclude that the OPDs shown in Figs. 5(a) and 5(b) are typical for given type of heat source in the framework of linear heat diffusion approximation.

Table 2

Results of calculations for eight types of heat sources.

Type of sourceMT (1/km)rms OPDnpMapr2SRapr

Note: MT, optical power; rms OPDnp, root mean square of nonparaxial OPDnp; SRapr, Strehl ratio; Mapr2, beam quality parameter; Pheat=1  W, cylinder sample made of fused silica, Daper=50  mm, and b=10  mm

According to theoretical predictions and numerical examples showed here, the values of resulting rms(OPD) and MT are proportional to effective absorption or more generally to heat power deposited in the element. Thus, knowing parameter MT,0 calculated for given heat power Pheat,0, we can determine the optical power MT,1 for the different heat power Pheat,1 as follows:


where κ1,0=Pheat,1/Pheat,0 is the ratio of heat powers for cases 0 and 1.

The remaining metrics of optical quality (Strehl ratio SRapr and parameter Mapr2) are highly nonlinear with respect to σ. However, applying the same approach, we can determine the similar relations for SRapr,1, Mapr,12 knowing the parameters SRapr,0, Mapr,02 as follows:



Let us notice that the simplest way to mitigate TOEs is increasing the laser beam size in aperture, e.g., taking collimator of longer focal length and wider aperture. As a rule, lowering power density, the temperature increase basically diminishes.10,14 As was shown in Refs. 15 and 16, the average temperature increase is proportional to absorbed power and has nonlinear dependence on laser beam diameter.

An alternative method of TOEs diminishing is application of more flattened than Gaussian laser beam profiles. As shown above, the profile SG5 should give much lower rms(OPDnp) comparing to Gaussian one. On the other hand, the typical beam shaper used to transform SG1SG5 consists of at least two optical elements (see Refs. 6 and 7) with the first one exposed to Gaussian beam. For complex, multielement optical trains, or multipass systems, e.g., beam cleaners,10,11 dynamic mitigation of TOEs of all system can be achieved by application of appropriate beam shaper and playing with mirrors and refractive elements. It is also important for transient, nonstationary regime of operation typical, e.g., laser weapon engagement.

In summary, we have prepared the numerical method to estimate the TOEs for wide class of optical elements valid for linear heat equation and small absorption approximation. In first step, we have to solve numerically the problem for given type of heat source, absorption efficiency, and sample geometry. The results can be rescaled for different absorptions or heat powers applying Eqs. (14) and (15). The paraxial thermal lensing (determined by MT) can be compensated by defocusing. The thermo-optic distortions can be transient, time-dependent functions, and for short operation time and relatively low duty, factor can be low as well. The beam quality deteriorates during operation up to the worst case of stationary value typical, e.g., for industry applications. Equations (14) and (15) can be applied also for the estimation of those unstationary, transient effects, knowing additionally time constants of elements. Let us insist on the limits of above approximations. It is valid only for a case of linear, heat equation of constant coefficients, small stresses, and low absorption. Moreover, for higher (comparable with wavelength) variances of OPD, the approximated Eq. (13) is not valid.



To achieve the high beam quality and flattened profile in near-field, several beam profiles were analyzed. The SGB SGp of index p=5 was found as the best compromise.

Furthermore, we have developed the simplified numerical–analytical model of TOEs to estimate acceptable level of heat power dissipated on surfaces and in volumes of optics. Such analysis was performed for several SGp profiles. In first step, we have to solve problem numerically for given type of heat source, absorption, and sample geometry. The results can be rescaled for different absorption or heat powers applying approximated formulas derived in the paper. The model can be applied also for the estimation of unstationary, transient effects. Moreover, the method of dynamic compensation of TOEs due to application of the appropriate beam profile and combinations of lenses and mirrors was proposed.


The work was supported by the National Centre for Research and Development of Poland in the framework of strategic program DOB-1-6/1/PS/2014 and project PBS3/B3/27/2015. We thank Dr. Zbigniew Zawadzki from the Institute of Optoelectronics MUT for critical discussion of the work.


1. T. S. Ross, Laser Beam Quality Metrics, SPIE Press, Bellingham, Washington (2013). Google Scholar

2. S. De Silvestri et al., “Solid-state laser unstable resonators with tapered reflectivity mirrors: the super-Gaussian approach,” IEEE J. Quantum Electron. 24, 1172–1177 (1988). https://doi.org/10.1109/3.240 Google Scholar

3. A. Parent, M. Morin and P. Lavigne, “Propagation of super-Gaussian field distributions,” Opt. Quantum Electron. 24, S1071 (1992).OQELDI0306-8919 https://doi.org/10.1007/BF01588606 Google Scholar

4. V. Bagini et al., “Propagation of axially symmetric flattened Gaussian beams,” J. Opt. Soc. Am. A 13, 1385–1393 (1996).JOAOD60740-3232 https://doi.org/10.1364/JOSAA.13.001385 Google Scholar

5. M. Santarsiero and R. Borghi, “Correspondence between super-Gaussian and flattened Gaussian beams,” J. Opt. Soc. Am. A 16, 188–190 (1999).JOAOD60740-3232 https://doi.org/10.1364/JOSAA.16.000188 Google Scholar

6. M. F. Dickey, Laser Beam Shaping, Theory and Techniques, CRS Press Taylor & Francis Group, London (2014). Google Scholar

7. J. A. Hoffnagle and C. M. Jefferson, “Design and performance of a refractive optical system that converts a Gaussian to a flattop beam,” Appl. Opt. 39, 5488–5499 (2000).APOPAI0003-6935 https://doi.org/10.1364/AO.39.005488 Google Scholar

8. “COMSOL multiphysics v. 5.3,”  https://www.comsol.com/products Google Scholar

9. N. Hodgson and H. Weber, Optical Resonators, Fundamentals, Advanced Concepts and Applications, Springer, New York (1997). Google Scholar

10. W. Koechner, Laser Engineering, Springer, New York (1999). Google Scholar

11. K. L. Dooley et al., “Thermal effects in the input optics of the enhanced laser interferometer gravitational-wave observatory interferometers,” Rev. Sci. Instrum. 83, 033109 (2012).RSINAK0034-6748 https://doi.org/10.1063/1.3695405 Google Scholar

12. C. L. Mueller et al., “The advanced LIGO input optics,” Rev. Sci. Instrum. 87, 014502 (2016).RSINAK0034-6748 https://doi.org/10.1063/1.4936974 Google Scholar

13. L. Pinard et al., “Mirrors used in the LIGO interferometers for first detection of gravitational waves,” Appl. Opt. 56, C11–C15 (2017).APOPAI0003-6935 https://doi.org/10.1364/AO.56.000C11 Google Scholar

14. A. K. Cousins, “Temperature and thermal stress scaling in finite-length end-pumped laser rods,” IEEE J. Quantum Electron. 28, 1057–1069 (1992).IEJQA70018-9197 https://doi.org/10.1109/3.135228 Google Scholar

15. J. Penano et al., “Optical quality of high-power laser beams in lenses,” J. Opt. Soc. Am. B 26(3), 503–510 (2009).JOBPDE0740-3224 https://doi.org/10.1364/JOSAB.26.000503 Google Scholar

16. B. Fakizi et al., “Laser heating of uncoated optics in a convective medium,” Appl. Opt. 51(14), 2573–2580 (2012).APOPAI0003-6935 https://doi.org/10.1364/AO.51.002573 Google Scholar


Jan K. Jabczyński is a full professor at Military University of Technology (MUT), Warsaw, Poland. He received his MSc Eng, PhD, and DSc, 1981, 1989, and 1997, respectively. Since 1982, he has been a worker at the Institute of Optoelectronics, MUT. He has been a tutor of applied optics, laser optics, and propagation at MUT. He is the author for about 190 journal papers and conference reports and one handbook. His current research interests include laser optics, coherent optics, high-power laser systems, and diode-pumped lasers.

Mateusz Kaskow received his MSc Eng in electronics from Wroclaw University of Technology in 2009 and his PhD in electronics from MUT in 2015. He was with the Institute of Optoelectronics from 2011 to 2017. He is the author for about 40 journal papers and conference reports. His current research interests include diode-pumped lasers, nonlinear optics, and computer graphics.

Lukasz Gorajek is an assistant professor at MUT. He received his MScEng and PhD degrees in electronics from MUT in 2009 and 2014, respectively. He is with the Institute of Optoelectronics since 2010. He is the author for about 60 journal papers and conference reports. His current research interests include diode-pumped lasers, nonlinear optics, and high-power laser systems.

Krzysztof Kopczyński is an assistant professor at MUT, Warsaw, Poland. He received his MSc Eng and PhD in 1985 and 2000, respectively. Since 1994, he has been with the Institute of Optoelectronics, and for the last eight years, he has worked as a director. He is the author of about 150 journal papers and conference reports and 10 patents. His current research interests include high-power laser systems, laser spectroscopy, lidar systems, diode-pumped lasers, and optoelectronic systems for safety and defense.

Waldemar Zendzian is a full professor at MUT, Warsaw, Poland. He received his MSc Eng degree in technical physics in 1984 and PhD and DSc degrees in electronics in 1993 and 2006, respectively. He is with the Institute of Optoelectronics MUT since 1985. He is the author of about 150 journal papers and conference reports and two handbooks. His current research interests include diode-pumped lasers, nonlinear optics, and quantum electronics.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Jan K. Jabczyński, Mateusz Kaskow, Łukasz F. Gorajek, Krzysztof Kopczyński, Waldemar Zendzian, "Modeling of the laser beam shape for high-power applications," Optical Engineering 57(4), 046107 (13 April 2018). https://doi.org/10.1117/1.OE.57.4.046107 Submission: Received 6 December 2017; Accepted 21 March 2018
Submission: Received 6 December 2017; Accepted 21 March 2018


Incoherent and coherent coupling for multi-slab solid laser
Proceedings of SPIE (January 18 2006)
New type of unstable resonator
Proceedings of SPIE (June 01 1990)
Beam combining with using beam shaping
Proceedings of SPIE (May 07 2012)
Measuring a laser focal spot on a large intensity range...
Proceedings of SPIE (August 24 2009)

Back to Top