Laser beam quality metrics like M2 can be used to describe the spot sizes and propagation behavior of a wide variety of non-ideal laser beams. However, for beams that have been diffracted by limiting apertures in the near-field, or those with unusual near-field profiles, the conventional metrics can lead to an inconsistent or incomplete description of far-field performance. This paper motivates an alternative laser beam quality definition that can be used with any beam. The approach uses a consideration of the intrinsic ability of a laser beam profile to heat a material. Comparisons are made with conventional beam quality metrics. An analysis on an asymmetric Gaussian beam is used to establish a connection with the invariant beam propagation ratio.
© 2012 OSA
CorrectionsHarold C. Miller, "A laser beam quality definition based on induced temperature rise: erratum," Opt. Express 21, 5635-5635 (2013)
High power laser beam characterization data often include a far-field power-in-the-bucket (PIB) curve [1–3]. The PIB data is compared with an ideal curve calculated assuming a perfect (i.e. uniform phase and intensity) near-field distribution with the same circumscribed profile and power as the beam being measured. A beam quality number is assigned by comparing the curves using some choice of an acceptable bucket width or useful power level. Much more common in the laser community is the use of camera, knife-edge, or other beam-size measurement approaches to acquire data that are fit to Gaussian-beam propagation formulas to arrive at an M2 or related beam quality metric. The primary issue with the current situation is the lack of consistency in defining laser beam quality and the possibility that a reported beam quality number may not adequately describe how the power in a beam is distributed. For example, it is possible for a PIB beam quality number to disregard a large fraction of the power in a beam, or for a second-moment measurement to misrepresent the tangible size of an intensity distribution. Rather than using beam-size measurements or arbitrary comparisons of PIB curves, it might be useful to instead arrive at a beam quality number by looking at how a far-field beam would interact with a material. The following sections motivate an alternative beam quality definition based on a consideration of the maximum-possible laser-induced temperature rise in an ideal absorbing volume.
The PIB method reduces a 2D far-field intensity distribution to a one-dimensional integrated-intensity curve, providing a simple way to visualize the distribution of power in a beam. Formally, the PIB is an area integral over a far-field intensity distribution behind a real or imagined “bucket” (aperture) of a given size. Restricting ourselves to circular buckets,3]. In order to keep the analysis independent of the physical size of the far-field distribution, the bucket radii are usually expressed in terms of an angle coordinate (θ), instead of the radial coordinate, r [cm]. θ × in Eq. (1) is the corresponding integration variable. The scaling of θ depends on the nominal diameter (D) of the near-field beam, the laser wavelength (λ), and the far-field distance (f). For example; an ideal Gaussian beam distribution can be expressed in terms of the angular coordinate using,2]. For Gaussian-like beams, the size of the near-field can be assessed using a second-moment, knife-edge, or other standard measurement approach. However, applying beam-size metrics to arbitrary beam shapes (arrays, annular beams, etc.) can lead to an improper description of the near-field profile. In such cases, it is preferable that the near-field beam profile be defined by a physical aperture, whose size and shape can be optimized for best far-field performance . It is important to note that the ideal far-field intensity distribution and a corresponding beam quality number will always depend on how the ideal near-field beam is defined. For example, a particular application might seek to compare a measured annular beam with a top-hat beam instead of an ideal annular beam.
Iave in Eq. (1) is the azimuthally-averaged intensity at a particular angle,
3. Laser-induced temperature rise
In the pioneering work of Lax , the heat equation was used to calculate the steady-state temperature rise of a semi-infinite absorbing volume exposed to an axially-symmetric laser beam. At the time, the interest was in light-induced semiconductor material degradation. For the purposes of assessing beam quality, our imaginary absorber will be placed at a far-field plane of the laser beam. The geometry is defined such that a beam propagates in air along the z-axis from a near-field plane at z = f<0 before striking an infinite flat absorbing surface on the x-y plane at z = 0. The material occupies all points z≥0.
Lax found that the steady-state induced temperature-rise distribution in the material could be written,5].
ΔTmax is the temperature rise at the origin in the extreme case when all of the power is absorbed by an infinitesimally thin surface layer, meaning that,
It turns out that the underlying formalism is applicable to all beam shapes and not just to those with axial symmetry. This is due to the fact that temperature rise, ΔT(x,y,z), at any point in the material point depends on the superposition of the contributions of each differential laser-induced heat source element, regardless of where it is located. The heat equation, which reduces to a Poisson equation in steady-state, can be written,Eq. (12) for the semi-infinite geometry ,
When all of the power is absorbed on the input surface, the heat source has the form,Eq. (15) into Eq. (13), one finds that the maximum attainable temperature rise is,Eq. (13) and Eq. (15) can be used to locate the temperature-based center of the beam. In cylindrical coordinates, the maximum temperature is given by,Eq. (9), scaling as the intensity-weighted average inverse distance from the center of the beam to each differential heat source element on the surface. It is not necessary for the purposes of this paper to solve Eq. (13) for the temperature profile throughout the entire absorbing volume. Lax and others have performed the calculation for a variety of configurations associated with materials processing [4–8].
4. Laser beam quality definition
Laser beam interactions with materials can be complex, involving details like reflectance, temperature dependencies of the material parameters, melting, cooling, exposure times, etc. However, the intrinsic ability of a far-field distribution to induce the maximum-possible centerline temperature-rise in an ideal absorber is a property that is related by Eq. (9) to an effective spot-size of the beam. Our semi-infinite absorber is “ideal” in the sense that 1) its material properties do not change with temperature, 2) all of the power in the beam is absorbed, 3) there is no loss of heat across any boundary, and 4) the rapid flow of heat produces a steady-state temperature profile, regardless of the energy absorbed. In such an absorber, the scale factor ΔTmax depends explicitly on the shape of the intensity distribution and scales linearly with the laser power. Also, the dependence of ΔTmax on the beam-shape remains the same regardless of the material properties, κ and W. For these reasons, one might consider the alternative beam quality definition,
The beam quality (BQ) definition in Eq. (18) depends only on the measured and ideal beam-shape functions and not on the beam power or any property of the absorbing material. The definition can be used with far-field beams of any shape. For unclipped circular Gaussian beams, it is easy to show using Eq. (17) and Eq. (18) that BQ = M2, consistent with the standard second-moment-based definition. For general beam shapes, BQ≠M2, so a more general kind of beam analysis is required. From the induced temperature rise perspective, the effective far-field performance of any beam scales as BQ−1. The definition captures the contributions due to all of the power (heat sources) in the beam. Therefore, a beam quality number determined using Eq. (18) only depends on how the ideal near-field distribution is defined.
5. Interpretation of power-in-the-bucket data
Since BQ = M2 for axially-symmetric Gaussian beams, standard measurement methodologies and instrumentation are entirely consistent with the induced-temperature-rise beam quality definition. PIB analysis using circular buckets is compatible with all beams since the method samples the far-field distribution in a radial fashion, consistent with the definition of ΔTmax. Using Eq. (5), Eq. (17) and Eq. (18), the beam quality can be expressed,Eq. (16) or Eq. (17). When analyzing a measured curve using Eq. (19), attention must be paid to ensure that scatter in the data does not upset the natural convergence of the integrands to a finite number as θ approaches zero. A simple curve fit using, for example, a Gaussian-based PIB function could be used in this region. The θ −1 dependence in Eq. (19) reduces the contribution due to power at large radial distances from the center of the beam. Therefore, in cases where the entire beam profile is not captured by a PIB measurement , a reasonable extrapolation (i.e. linear, Gaussian-based PIB) of the PIB curves will usually introduce a small and predictable level of uncertainty to a BQ determination.
6. Comparison with conventional beam quality definitions
In this section, a simulated “measured” beam is used to review how conventional PIB analysis is used to arrive at beam quality numbers. We then demonstrate how the various metrics behave as the shape of a far-field beam is changed. The ideal PIB curve in Fig. 1 was generated by assuming that the near-field is a circular top-hat beam of diameter, D, producing an Airy diffraction pattern in the far-field. The measured PIB curve is derived from a “core-pedestal” far-field distribution. The core is approximated by a narrow Gaussian beam with width θc = 0.85 λ/D, which is close to the width of the central lobe of the Airy pattern, while the pedestal is a much broader Gaussian (θp = 12.2 λ/D). The measured PIB curve in Fig. 1 assumes that 25% of the total power is in the core beam and 75% of the power is in the pedestal.
The conventional vertical beam quality (VBQ) metric  defined in Eq. (20) below is often used to compare the amount of power contained within in the first zero (θ0 = 1.22 λ/D) of the ideal-beam Airy pattern. For the curves shown in Fig. 1, the VBQ number would be,3] metric specifies an arbitrary threshold power level (PIBth), then compares the bucket radii (θmeas(PIBth) and θideal(PIBth)) required to capture the specified power. With PIBth set at 0.5, the HBQ of our simulated measured beam is, by inspection of the curves in Fig. 1,Eq. (19) to the same curves gives BQ≈3. The M2 of the beam cannot be calculated since the second-moment width of the Airy distribution is infinite. From this simple example, it is clear how a choice of metric, bucket size or threshold power level can dramatically change the interpretation of a beam quality number.
To further clarify the behavior of the various definitions, calculations were done on an even simpler case where the ideal beam is a Gaussian distribution with width, w0 = 1.0 (arbitrary units) and the measured beam is again fabricated using the sum of a broad pedestal (wp = 10.0) and a narrow core (wc = 1.0) Gaussian beam. The near-field is assumed to be a circular Gaussian sized such that w0 = 1.0. By varying the fractional power in the core beam, one can scan the shape (and beam quality) of the simulated measured far-field intensity distribution. Curves showing the trends in the various beam quality definitions as the fractional power in the core beam are shown in Fig. 2 . PIBth used in the HBQ calculation was again set at 50% of the total power. The bucket radius used for the VBQ determination was arbitrarily set at w = 1.22.
Note the convergence in Fig. 2 of the VBQ2, HBQ and BQ curves above about 60% of the power in the core beam. The VBQ metric clearly over-predicts beam performance from both the average intensity (VBQ2) and temperature-rise (BQ) perspectives. The VBQ2 curve diverges quite dramatically from the BQ curve as the fractional power in the core is reduced. Although not shown Fig. 2, the VBQ2 converges to a limiting value of ~30 as core-beam power goes to zero. This value is consistent with the peak intensity level of the pedestal Gaussian averaged over the small bucket area. At this point, less that 1.5% of the total power in the beam is considered by the metric.
The HBQ curve in Fig. 2 displays a dramatic inflection as the total power in the core region falls below PIBth. Since the beam intensity distribution is changing smoothly, this behavior demonstrates how inconsistencies in an HBQ measurement can arise if PIBth is not well-matched to the shape of the beam being measured. As the core power approaches zero, the HBQ converges toward a more sensible value determined only by the width of the Gaussian pedestal beam.
The M2 curve is generated by using a second-moment calculation on the core-pedestal beam to determine an effective far-field beam waist (wmeas), which is compared with the ideal Gaussian beam waist using M2 = wmeas/w0. The curve agrees with expected beam quality at the endpoints where the far-field beam is a single Gaussian. However, the rest of the curve clearly demonstrates how second-moment-based beam size determinations on unusually-shaped distributions can over-predict beam sizes and lead to incorrect assessments of beam quality.
The BQ curve calculated using Eq. (19) is the only metric that appears to remain consistent over the entire range of beam shapes. On both endpoints, the metric agrees with the M2 of the respective Gaussians, while in between there are no inflections or divergent behavior as the shape of the far-field is scanned. The BQ also displays good agreement with the conventional HBQ and VBQ2 metrics in the core-dominated region where we would expect them to have the most validity.
7. Elliptical Gaussian beam example
This section will show how the temperature-based beam quality is connected with the invariant beam propagation ratio in Gaussian beams with non-circular far-field profiles. The treatment of general astigmatic beams and testing methods is outlined in an International Standards Organization (ISO) document . The approach supports the determination the ten second-order moments to arrive at an effective beam propagation ratio, the intrinsic astigmatism, and the twist parameter. For simple astigmatic beams, the invariant effective beam propagation ratio (usually considered the beam quality) is given by,
We first consider the case of a Gaussian beam propagated from a circular near-field but having an elliptical far-field intensity distribution due to a difference in M2 along the two principal axes of the beam. Any elliptical Gaussian intensity distribution can be written,Eq. (23) into Eq. (17), one arrives at the expression [7,8],Eq. (19) for a circular near-field beam using Eq. (24) is given by,Eq. (25) does not agree with the prediction of Eq. (22) is expected, since the ability of a beam to raise the temperature of a material is not, by definition, a beam invariant. Intuitively, one would expect the centerline heating to be reduced by beam asymmetry. Therefore, it is instructive to consider the ratio,Eq. (25) and Eq. (26), the beam quality can now be written,Equation (27) provides a connection between the BQ the beam invariant for a circular near-field.
In the more general case, where both the measured and the ideal beams are elliptical in the far-field plane, Eq. (24) can again be used to arrive at a formula for the beam quality. Using similar arguments as above, one can also show that impact on heating due to beam asymmetry is now, more generally,Equation (27) provides that if and are measured using standard techniques, the BQ can be calculated using Eq. (28).
An alternative laser beam quality definition (BQ) has been proposed that offers an unambiguous way to compare a measured far-field beam with an ideal beam. The definition is based on the intrinsic ability of a laser beam to induce the maximum-possible steady-state temperature rise in an ideal semi-infinite absorbing volume. The definition weights the contributions to the temperature rise over the entire far-field beam profile, which can be of any shape. The BQ of an axially-symmetric Gaussian beam under the new definition is equal to the invariant beam propagation ratio (M2). For elliptical Gaussian beams, the beam invariant is related to the BQ through the far-field beam asymmetries. For general beam shapes, the definition supports the use of circular far-field power-in-the-bucket analysis to determine laser beam quality. The concepts also provide a straightforward way to extract a BQ number from power-in-the-bucket data.
The author would like to thank Jack Slater and Alan Paxton for their useful comments.
References and links
1. A. E. Siegman, “How to (maybe) measure laser beam quality,” in DPSS (Diode Pumped Solid State) Lasers: Applications and Issues, M. Dowley, ed., 17, OSA Trends in Optics and Photonics (OSA, 1998), paper MQ1.
2. C. X. Yu, S. J. Augst, S. M. Redmond, K. C. Goldizen, D. V. Murphy, A. Sanchez, and T. Y. Fan, “Coherent combining of a 4 kW, eight-element fiber amplifier array,” Opt. Lett. 36(14), 2686–2688 (2011). [CrossRef] [PubMed]
3. J. M. Slater and B. Edwards, “Characterization of high power lasers,” Proc. SPIE 7686, 76860W, 76860W-12 (2010). [CrossRef]
4. M. Lax, “Temperature rise induced by a laser beam,” J. Appl. Phys. 48(9), 3919–3924 (1977). [CrossRef]
5. M. Lax, “Temperature rise induced by a laser beam II. the nonlinear case,” Appl. Phys. Lett. 33(8), 786–788 (1978). [CrossRef]
6. Y. Lu, “Square-shaped temperature distribution induced by a Gaussian-shaped laser beam,” Appl. Surf. Sci. 81(3), 357–364 (1994). [CrossRef]
7. Y. I. Nissim, A. Lietoila, R. B. Gold, and J. F. Gibbons, “Temperature distributions produced in semiconductors by a scanning elliptical or circular cw laser beam,” J. Appl. Phys. 51(1), 274–279 (1980). [CrossRef]
8. J. E. Moody and R. H. Hendel, “Temperature profiles induced by a scanning cw laser beam,” J. Appl. Phys. 53(6), 4364–4371 (1982). [CrossRef]
9. International Standards Organization, “Lasers and laser-related equipment-test methods for laser beam widths, divergence angles and beam propagation ratios. part 2: general astigmatic beams,” Ref. ISO 11146–2:2005(E) (2005).