Gas-assisted laser cutting of medium-section metals using spherically aberrated beams

We analyzed the effects of the focal point aberrational offset in optical transport systems for highpower lasers. Transverse and near-axial laser intensity distribution transformations in the presence of both positive and negative spherical aberrations were numerically calculated and experimentally demonstrated for different strengths. We show that spherical aberration yields considerable asymmetry of the focused beam’s caustic. Several optical transport systems with identical optical parameters (excluding the noncorrected axial beam spherical aberration) were designed. We examined the effects of the laser intensity profiles produced by these systems on the quality of oxygen-assisted laser cutting of medium-section mild steel. We show that high-quality cuts can be obtained for different shapes of laser intensity distribution. However, the greater the refocusing magnitude introduced by the spherical aberration correction, the more precisely the focal point position must be maintained during the laser cutting process. © 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. [DOI: 10.1117/1.OE.54.4.044103]


Introduction
Laser cutting of metals and alloys is one of the most commercially valuable laser technology processes.3][4][5] Omitting the energy and resource efficiency issues, 6 the performance of laser cutting can be captured by only two output parameters.These are the cutting speed and kerf quality.Thus, laser cutting optimization reduces to achieving the maximal cutting speed with the best kerf quality.As far as the material separation is concerned, 7 the speed can be dramatically increased.On the other hand, the cutting quality depends on many parameters.The most important and frequently discussed are the kerf geometry (accuracy, width, and taper), the surface quality (cut-edge roughness, striations morphology, and dross and burr inclusions), and the mechanical and metallurgical properties (hardness and strength, heat-affected zone, and oxide layer). 8The number of parameters dramatically affecting the cutting quality ranges from 25 9 to 75. 4,5These are nearly evenly divided between the workpiece properties, assist gas parameters, laser machine, and laser beam characteristics.Input-output parameters relationships are mostly nonlinear and exhibit strong interdependence.In this regard, maximal precision is provided by the studies that focus on only one input-output parameter pair.Yet, this "one parameter at one time" approach does not provide complete information on the process as a whole.
Relatively little is known regarding the aberrational distortion of laser beams due to the specifics of design and operation of transport systems of high-power laser cutters.This knowledge gap becomes critical when considering the dependence of industrial cutting quality on a particular laser intensity distribution (LID) formed by an optical system of the cutter in the laser-matter interaction region.The significance of focal point position (FPP) distance and possible range of its variation during gas-assisted laser cutting remain poorly understood.This motivated us to perform systematic studies of LID transformation in caustics with spherical aberration (SA).The obtained data were used for designing an array of optical transport systems with identical optical parameters, excluding the noncorrected axial beam SA.These systems, yielding considerably different LIDs, were placed in the same laser cutter.The root-mean-square magnitude of an unbalanced spherical wave aberration (WA rms ) was used as the input parameter.Dross height was chosen as the output parameter.It exhibits solidified residual melt which remains attached to the kerf-edge bottom.Together with the kerf roughness, the dross comprises a serious fault, demanding considerable postlaser treatment of the processed tools.
This paper is organized as follows.In Sec. 2, we briefly review the state of the field and the problem.Section 3 presents the study framework.Simulations of near-caustic LID shapes in the presence of SAs of various strengths and signs are reported for several fiber sources typically used in laser cutting machines.The basic findings are discussed.In Sec. 4, we describe the experimental results.First, we describe the experimental setup and methodology.Next, we characterize the beams tailored using SA.Finally, we discuss trial cuts made using these beams.Summary and conclusions are presented in Sec. 5.

Speed-Power-Pressure Concept
A literature review suggests that the main factors affecting a laser cutting quality are the cutting speed, the laser power, and the assist gas pressure.Radovanovic and Madic 10 reported that different combinations of these parameters have been considered in up to 67% of the cited papers.Studies addressing the joint impact of five input parameters comprised only 6%.Some studies addressed optimization by considering different combinations of input parameters.Hashemzadeh et al. 11 introduced specific point energy, SPE ¼ Power × Beam Size∕Speed and ðSpeed∕PowerÞ 2 .Orishich et al. 12 used Peclet number, Pe ¼ Speed × Kerf Width/Heat Diffusivity and density of the absorbed laser energy averaged over the kerf volume.The working hypothesis of the final kerf relief formation is solidification of the unstable melt flow. 13,14Detailed investigations by Hirano and Fabbro 15,16 show that this instability is a sum of thermodynamic instability (related to the localized melting) and hydrodynamic instability (caused by the imbalance of assist gas flux and surface tension forces).The instability pattern varies along the melt flow, causing the complicated final relief of the kerf surface.Considerable physical complexity justifies the existence of the literature dedicated to empirical investigations of quality parameters.Techniques such as analysis of variances, [17][18][19][20] response surface method, [20][21][22][23] multiple regression analysis, 9,24,25 artificial neural networks, 9,[26][27][28] and genetic algorithms 22 have been used in these studies.Results are captured as polynomials that relate a single output parameter to a set of input parameters.The studies using the cutting speed, laser power, and assist gas pressure as input parameters comprising, respectively, 100%, 77%, and 62%.Kerf surface roughness was considered as the output parameter in 70% of the studies.Kerf taper and width were considered in 46% of the studies.Yet, the obtained polynomials had different functional forms.Thus, Zaied et al. 19 reported results using a single quadratic term, but others obtained polynomials that included from 5-10 9,21,24,28 to 15 20,22,23 terms.Discordant polynomial forms are particularly appealing in light of the dramatic effect of unaccounted factors such as the environment, 24 workpiece parameters, 29 and laser beam characteristics.

Laser Beam Shape, Position, and Distortion:
Causes, Effects, and Definitions The most explored area of laser beam characteristics impact relates to the energy balance studies utilizing "cutting speedcut front inclination angle-absorption of polarized lightthermal dynamics versus quantity of absorbed energy-quality parameters" models.Summarizing the simulations, 15,30 experiments, 13,16 and articles combining experiments and theory, 31,32 the best speed-quality relationship is achieved when the melt front inclination angle provides absorbed energy supporting maximally uniform melt flow in all melt pool regions, with minimal melt accumulation.The relationship between the averaged melt front inclination angle, kerf geometry, and laser beam size governs the averaged fraction of the beam power involved in the cutting process. 12,13 facto, the melt pool surface is a complex three-dimensional wavy humped structure with strongly varying local characteristics.Existing studies [7][8][9][10][11][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29] do not consider the spatial characteristics of laser beams. Neverteless, these are accounted for by the FPP, which is one of the key factors of the cutting process.The FPP denotes the beam waist location relative to the workpiece's front surface.Accordingly, previous studies were performed with FPPs that yielded the best laser cutting performance.The details of this approach are described by Orishich et al. 12 Only in one study was the FPP used as the input optimization parameter.28 It is obvious that the beam waist should always be located on the thin-section workpiece's surface (FPP ¼ 0) during cold ablation cutting at extremely high intensities produced by sharp beam focusing.33 In gas-assisted laser cutting of medium-and thick-section metals, the beam waist rarely coincides with any of the material surfaces of the cutting sheet or melt pool.In these cases, determining the optical components yielding the minimal waist size of the focused beam 34 is tangential.Wandera 35,36 showed that in the case of high beam quality fiber lasers, the FPP scales with the sheet thickness, consistent with gas and fluid dynamical aspects of melt removal from the wider kerf (see also Ref. 37).At the same time, for any fixed FPP value, an excellent cut, good cut, dross, and no cut occur with nearly equal frequency.In the case of disk lasers 35,36 characterized by considerably worse beam quality, no "cutting quality-thickness-FPP" correlation exists.Inconsistent results have been obtained for cutting aluminum alloys.38 FPP variation affects the Gaussian beam size in the interaction region without affecting the transverse LID.39 In real applications, the LID shape strongly varies along the propagation axis because the laser beam and the optical system are not ideal.Hence, the best FPP search amounts to finding the beam characteristics that provide an optimal (not maximal) load of the laser power into the laser-matter interaction region.40 The effect of the LID shape on the performance of laser technology processes has been repeatedly confirmed in experiments.Some papers directly attributed the cutting quality to fiber [2][3][4][5]18,35 and conventional 41,42 laser beam quality.Conversely, many studies reported considerable performance improvement by using non-Gaussian beams.Thus, considerably increased cutting speed is predicted when using coma-like beams. 43A 20% increase in maximal dross-free cutting speed was experimentally demonstrated using bifocal optics.][48][49] Different detectable keyhole shapes were observed for Gaussian and tophat beams.50,51 Applications of doughnut and annular beams yielded a 20% increase in stainless steel cutting speed and higher quality holes drilled in titanium foils.52 Common to all these studies is the absence of the beam characteristics choice concept as well as the lack of information on its actual shape for defocusing range used in every actual process.
It is well known 53,54 that aberrations in optical systems powerfully affect caustic shape distortion.In actual laser cutters, SAs of considerable strength are likely present when other primary aberrations are negligibly small.So-called design and manufacturing SAs are usually retained in laser cutter transport systems during optical design and manufacturing.Furthermore, thermally induced distortions in high-average power laser optical systems 55,56 in typical cases of circular beams and optics are mostly characterized by spherical components. 57Scarce research on aberrational distortions in laser beams 58,59 is likely owing to the apparent triviality of the problem.In addition, the beam quality factor 60 is widely used as a universal parameter.Thus, some studies address M 2 degradation 61,62 instead of the beam shape transformation.The M 2 value in question of the beam propagation reasonably describes near-Gaussian beams. 51,63,64n this case, the beam size is adequately determined using second-order intensity moments. 65Determination of 4σsize of noticeably non-Gaussian beams becomes too time consuming.Express techniques 66,67 are not regularly applied.Introducing M 2 into the Gaussian beam equations suggests the caustic symmetry 37,62,64 that is never asserted for an aberrated beam.Furthermore, different beams may be characterized by the same M 2 . 68These facts generated a certain scepticism as to the versatility of the M 2 approach, expressed even by the author of this concept in his subsequent papers. 69ualitatively, it is evident that any M 2 degradation 60 as well as the dramatically rarely used Strehl ratio, 70,71 the overlapping integral, 71 and the V-parameter 66 witness some LID distortion.However, knowing the actual LID shape in an interaction region is desired.None of the above parameters, including the aberration strength in the units of Seidel and Zernike 57,72 factors, provide this knowledge.
It is clear that an LID shape distorted by aberration critically determines the laser cutting performance.An initially high beam quality and diffraction-limited optics do not satisfactorily determine the cutting quality.At the same time, the concept of speed-power-pressure implies neither no requirements on the transport system nor knowing the actual beam shape.

Spherical Aberration and Laser Intensity
Distribution Near the Beam Waist

Paraxial Focus, the Best Focal Spot, and Focal Point Position in Aberrated Systems
The analytical description of an aberrated beam spatial shape distortion was provided by Born and Wolf in Sec.9.4 of their monograph. 54The aberration is most evident as a deviation of the best focus spot (BFS) from the paraxial focus plane (PFP).Assuming this deviation is constant, it does not affect laser cutting.Critically, the axial and transverse LID in the aberrated caustic dramatically differ from those obtained in the aberration-free case.The aberrations' effect can be demonstrated simply by alluding to the ray model of geometrical optics. 53The rays normal to the surface of the distorted wavefront traveling from the exit pupil of an optical system intersect in the image space.The intersection density (i.e., intensity) depends on the wavefront distortion strength and is a function of both the axial and transverse coordinates.Consequently, cross sections at different points along the beam propagation direction are characterized by different transverse LID shapes.We define the BSF as the region with the maximal intensity and refer to "the spot" instead of "the point" because a region of indistinguishably equal intensity is always of finite size.Alternative methods of a BFS definition involve the averaging of the beam energetic characteristics.The size and position of a BFS defined in different ways may not be consistent but the difference is small compared to the actually used FPP ranges.An optical configuration of the refractive laser cutter heads is typically telecentric 73 and the image space is limited by the nozzle size.In principle, field aberrations do not relate to the cutting head optical performance parameters.Refocusing is easily performed to add defocus for achieving the required correction of the SA for maximizing the system's point-spread function.Coma and astigmatism are close to zero in the vicinity of axial ray propagation, even in cases of certain lens misalignment. 53,54,74,75The number of refractive surfaces of a typical transport system is sufficient for producing a diffraction-limited beam for apertures up to 0.20, assuming that the optical path is computed as a whole. 74,76The basic problem here is that the collimating unit and the processing head are usually purchased separately, often from different vendors.Assuming that the collimating unit produces a diffraction-limited beam, a WA rms of ∼0.05λ is introduced using an optimized focusing singlet with an output NA of 0.06.For an optimized focusing doublet, the same WA rms is achieved at an output NA of approximately 0.20. 74ccumulation of optical components in a transport system is prevented by nonthreshold thermal phenomena related to the laser radiation absorption in the bulk and on the surfaces.The simplest and most investigated phenomenon is the laser-induced focal shift, 64,77,78 which is directly proportional to the laser fluence.Among the optical component parameters, the most significant are the absorption factor and temperature-induced index variation, ∂n∕∂T. 56,77,7957]72 Chow et al. 72 estimated the wavefront distortion as a few hundredths of the wavelength per kilowatt per 1 cm length for state-of-the-art fused silica plano-plano windows.The same measurements for a commercial fourlens transport system 57 yield a wavefront error up to several wavelengths per kilowatt per total thickness of anti-reflection coated glass.Thermal phenomena significantly constrain the choice of glass, its sum length, and the number of refractive surfaces. 77,80This tradeoff between the optical and thermal performance of a system often does not allow for a complete correction of aberration, even during preliminary optical design.Contamination of optical surfaces crucially contributes to laser-induced beam deterioration. 80All contaminations, including the invisible ones, are a powerful source of laser heat dissipation.In this case, the speed-power-pressure-FPP correlation is lost and process failure occurs.Contaminations of reflective surfaces yield considerably smaller beam deterioration 72 because ∂n∕∂T is irrelevant to thermal lensing.Different scenarios of BFS and PFP relative locations and FPP range are schematically shown in Fig. 1.For an ideal beam propagating through an ideal optical system, the BFS 0 and PFP are colocalized at the axial origin, the transverse LID shape does not depend on the axial coordinate, and the beam is congruent over the entire propagation length. 39he waist size 2ω 0 defines the Rayleigh range z R .With spherical aberration, the BFS p , BFS N move in the direction opposite to the sign of wave aberration.This is also the direction in which the beam congruence is lost.We denote this region as the aberration location range (ALR).Loss of congruence causes the formation of an LID with complicated axial and transverse patterns.Outside the ALR, the beam remains congruent.Nonetheless, both inside and outside the ALR, the Gaussian transverse shape as well as the hyperbolic envelope are frustrated proportionally to the aberration strength.The transverse LID shape associated with the SA has rotational symmetry.Inside the ALR, an initially diffraction-limited beam transforms into a sharp-edged modulated pattern.Outside the ALR, it assumes a smooth wide-winged shape with a sharp maximum.The Rayleigh range of an aberrated beam z RP , z RN has no strict physical definition and is introduced here only as a measure of beam size variation.6][37][38] Depending on the waist size it can exceed z R by an order of magnitude.The FPP ¼ 0 case always implies that the BFS is placed on the upper surface of the sheet, as shown by the dashed line in Fig. 1.The FPP sign is determined by the coordinate choice.In our interpretation, the beam propagation direction was taken along the z axis with the zero FPP coordinate placed on the beam waist.In this way, the FPP was taken as positive when the beam waist was placed above the upper surface of the sheet.This means of framing the axes agrees with the signs' convention used in optics and consolidates the signs of all z, z o , and FPP axes for all of the figures in the present paper.[37][38]

Numerical Simulation of the Beam Caustic in the Presence of Spherical Aberration
Simulations of beam caustics with SAs were performed using standard optical design software.An integrated finite optical system was used as the model.It was allowed to inversely alter the sign and magnitude of the aberration by changing the position of one component, resulting in relatively small variations in the system parameters.This allowed analyzing the LID shapes in caustics with similar parameters in different configurations providing wave SAs with WA rms þ0.004, þ0.08, AE0.The blue thin curves refer to nearly zero SA, the black solid curves-to positive wave SA with 0.08λ WA rms , the purple solid and blank triangles-to negative and positive wave SAs with 0.55λ WA rms , the red solid and blank circles-to negative and positive wave SAs with 0.8λ WA rms .BFS positions for the case of the two lowest SAs are placed onto z 0 axis.For larger SAs, they can be found on z axis.The curves in each plot are normalized by the absolute maximal intensity obtained using the given source.
by the absolute maximal intensity obtained using the given source.Moving the component to positions providing maximal positive and negative aberrations caused a AE1.50% focal length variation.The back focal distance variation was as large as AE2.5%.Meanwhile, the maximal BFS variations were much larger owing to considerable defocus, and mostly depended on the source size.In this regard, and for clarity, we folded the curves onto each other and introduced the second longitudinal axis z 0 for the coordinate inside a given caustic.The axis z represents the absolute coordinate of the aberrated caustic.The zeros of both coordinates are placed in the waist of the caustic produced by an aberration-free reference sample.The BFS positions for the cases of the two lowest aberrations are placed on the z 0 axis.For larger aberrations, they can be found on the z axis.As shown in Fig. 2(a), every aberration of the optical system reduces the maximal intensity of the fundamental mode caustic, consistent with the Strehl ratio.However, Figs. 2(b)-2(d) show that the maximal intensity in the aberrated system exceeds those produced in diffraction-limited systems involving multimode-fiber sources for both positive and negative SAs.All aberration plots indicate that the maximal intensity curve is asymmetric with respect to the zero of z 0 coordinate.The ratio of the maximal intensity variation rate inside the ALR to that outside it depends on the source size.In our simulations, this ratio increased from 2 to 6 for the source size varying from 7 to 400 μm for the largest aberrations.This behavior of the calculated curves becomes clear upon considering the transverse LID shape dynamics, as shown in Fig. 3.For an aberration-free system, the transverse LID shapes may be considered symmetrical with respect to z 0 ¼ 0 [Figs.3(a) and 3(c)].In the presence of an aberration, beam truncation occurs for the sequence indicated by the arrows 1 to 4. This yields formation of M and П shapes inside the ALR.Outside ALR, the maximal intensity decreases much slower owing to the energy pumping from an approximately 1∕e 2 level into the peak (arrows 5).At the same time, the beam wings crawl off the axis, i.e., the full width at half maximum and 1∕e 2 beam sizes grow slower than the full size.It is important to note that in the presence of aberrations, the BFS cross section is always characterized by a Λ-shape LID, which is by no means the source's image [Figs.3(b) and 3(d)].These processes are common to both positive and negative aberrations, and the extent of their appearance depends solely on the relationship between the aberration strength and source size.
Our simulation results are consistent with classical analytic estimates 53,54 and agree with the results of other studies. 51,58,59,63Yet, the defocusing ranges are considerably smaller than the ones reported by Pu and Zhang. 58The claim that Gaussian beam focusing can be improved by introducing positive wave aberration 58 appears doubtful and requires reformulation.Simulation results by Park and Chung confirm this doubt. 59Yet these authors have not reported submaximal axial intensity inside the ALR.Experimentally, a focused beam LID transformation has been reported by Kaplan,51,63 but the author had not considered this finding.It is important to emphasize that an asymmetric caustic shape may cause considerable mistakes in Gaussian beam scaling for energetic estimates.

Experimental Setup and Methodology
Our experimental technique was based on conducting one parameter (WA rms ) at a time investigations.To this effect, two arrays (10, 11, 12, and 15, 16, 17) of identical transport systems with the same axial glass thickness, focal length, and NA were built.The numbering indicates our experimental schedule.The only difference between the systems within a given array was the value of the uncorrected wave SAs.This value was varied from 0.06 to 0.40 waves rms in BFS by altering the orientation and airspacing of the optical components.In this study, WA rms values were computed for each optical configuration by using standard optical design software.For the systems with smaller WA rms magnitudes, an additional experimental verification of the computed data was performed.For this, complete diffraction-limited reference analogs were assembled for the systems 10 and 15 using a different set of optical components.The Strehl ratio (S) was calculated as the ratio of the maximal waist intensity produced by the model system to the maximal waist intensity produced by its diffraction-limited correlate.
The resulting S value was used to determine the WA rms magnitude according to the famous approximated estimate WA rms ≈ ð2πÞ −1 ffiffiffiffiffiffiffiffiffiffi ffi 1 − S p . 54The computed and measured WA rms magnitudes deviated by 6% and 14% from those of transport systems with WA rms 0.06 and 0.08λ (experiments 10 and 15).For larger aberrations, the above estimate yields incorrect results. 54,70All components of the transport systems were made of Heraeus Suprasil 3002 and Schott N-BK7HT with bulk absorption of 0.25 ppm × cm −1 and 2 × 10 −4 cm −1 , respectively.We used commercial singlemode fiber lasers LC-0.4 and YLR-150/1500-QCW (International Photonics Group) with maximal continue wave powers of 400 and 250 W and output NA of about 0.05 (2ω s ¼ 14 μm, M 2 ¼ 1.17 and 1.05, respectively).In the systems 10-12 with linear magnification −1.70, aberrations were inserted at the collimation stage.In the systems 15-17 with linear magnification −1.05, aberrations were inserted at the focusing stage.The details of the transport systems can be found in our previous paper. 81Wave SAs produced by these systems were positive.
A schematic of the experimental setup is shown in Fig. 4. The monitoring unit was rigidly attached to the cutting head and employed the same optics that were used for the cutting head assembly.The collimated laser beam was split to target the monitoring unit using a system of high-optical-quality wedges.The first of the wedges was periodically inserted into the beam path to allow for quasi-online monitoring.Our estimate of a thermally induced aberration (TIA) effect by combining the heating and probe beams have yielded TIA coefficients as small as 0.5 × 10 −6 λ∕W and 1 × 10 −5 λ∕W per 1 cm thickness for Suprasil 3002 and N-BK7HT, respectively.The TIA on the wedge-reflective surface was smaller than the detection threshold of our apparatus because the expansion coefficient was more than 20 times smaller than ∂n∕∂T. 55,72In these conditions, the maximal TIA magnitudes were smaller than 10 −3 waves for the systems 11 and 12 made of Suprasil 3002 with a 21 mm total axial thickness at 400 W laser power.The systems 10, 15, 16, and 17 that yielded the combination of 14 mm N-BK7HT and 11 mm Suprasil 3002 demonstrated about 10 −2 waves total TIA strength.This allowed considering our trials free of thermally induced phenomena.The difference in optical paths in the working and monitoring arms was about 300 mm and did not affect any caustic shape alterations.The principal experimental complexity was associated with arranging the nonprotected charge-coupled device surface coplanar with the upper surface of the cut sheet.This was achieved by registering the shadow pattern from the sliding beam of the He-Ne laser.The maintenance error was about AE0.20 mm and was mostly caused by the sheet nonflatness.This error was considerably smaller than the length of an LID shape variation even for strongly aberrated beams.

Measured Profiles of a Cutting Beam
Figure 5 shows the LID shape trains registered in the monitoring channel for all transport systems and various axial coordinates.Apart from the difference caused by the varying NA values, the LID patterns for low-aberration systems 10 and 15 were qualitatively similar; thus, the LID train for the latter was omitted.The beam-waist size for the listed systems varied from 16 to 25 μm for systems 15 and 10 to 41-43 μm for systems 12 and 17.Regarding the differing aberration degradation of an initially Gaussian input beam, a Rayleigh range in its formal sense did not exceed 0.5 mm for each of the mentioned waist sizes.In the figure, this is marked in yellow and is included solely to indicate its correlation with the FPP range.The areas marked in blue and pink represent the FPP range, for which an excellent or acceptable cut quality was observed, respectively.

Experiments with Laser Cutting using Aberrated Beams
Low-pressure oxygen-assisted laser cutting experiments with the beams produced by the aforementioned optical systems were exactly conducted under the same conditions.The cut material was mild steel (C 0.14-0.22,Mn 0.3-0.6,Si < 0.05, P < 0.04, S < 0.05, Cr, Ni, Cu <0.3 mass % each).The assist gas was oxygen of 99.5% purity with 2.0-Bar pressure.We used the standard 1.0-mm Trumpf and Lasag nozzles with a 1.0-mm standoff distance.Material thickness was varied from 1 to 4 mm in 0.5-mm steps.The only parameter violating the one parameter at a time condition was a gradual decrease of the cutting speed with increasing sheet thickness.
The cut speed was 5 to 30 mm∕s and was kept constant for each thickness of the sheet.The cutting quality was estimated according to the minimal dross height.The cuts having dross height under 10 μm were qualified as excellent.
Those with a dross height ranging from 10 to 200 μm were qualified as acceptable.Cuts having a dross height above 200 μm were classified as laser cutting failures.The minimal cut width was registered in FPP ranges termed "best cut FPP" and was as small as 0.3 mm for the maximal cut-sheet thickness.In the ranges termed "acceptable cut FPP," the cut widths were increased by up to 0.5 to 0.7 mm for the thickest sheets.The cut geometry was a 50-mm-long line, and the piercing point was outside the analysis range.Our experimental results are shown in Fig. 6, classified according to the above criteria.The FPP ranges yielding excellent cuts are highlighted in blue, and those yielding acceptable cuts are highlighted in pink.The As in every laser cutting experiment, our data contain an unmistakable degree of noise.Nevertheless, a thorough inspection of the experimental results reveals the following: for each studied LID, it is possible to find a pair of positive and negative ranges of the best cut FPP that produces an excellent laser-cut quality.In a similar manner, the cut quality qualifying as acceptable is allowed within a pair of considerably wider ranges of the FPP.Both positive ranges are considerably wider than the negative ones.The LID profile itself does not noticeably influence the cutting quality if the FPP is adjusted in the best cut and acceptable cut ranges.The width of the ranges is inversely proportional to the extent of SA.These facts may be consistently explained by employing the following model.Because no method has been developed for reliably predicting the focal point position, the best cut FPP is intuitively chosen from the highlighted region.This empirically provides the best thermal and hydrodynamic balance in the melt pool.In the space between a pair of acceptable cut FPPs, a too high laser intensity yields strong thermo-and hydrodynamic instabilities.The cutting quality becomes poor and only high-speed material separation is possible.Outside this space, the cutting process gradually moves from low-speed separation to the lack of throughout cut due to low laser intensities.
Under the constant performance of the laser source and a constant NA, the melt flow destabilization may be caused only by beam intensity variations.As seen from Fig. 5 in the region of positive FPP, the LID shape variation is insignificant so that the intensity variation is mostly caused by changing the divergent beam size.In contrast, for a negative FPP, the melt pool interacts with the part of the convergent beam with a rate of maximum intensity variation several times greater.Furthermore, a high maximum intensity variation is combined with a deep transverse modulation of the LID shape.The sum of these two properties makes the system much more sensitive to focal point position fluctuations and narrows the FPP range.We could not find any reliable difference in the cutting quality provided by the transport systems 10 and 15.This suggests that the crucial magnitude of the wave SA that worsens the cutting quality was about 0.1λ rms in our trials (here we have taken into account the 0.01λ contribution of TIA in the transport system 15).Taking into account the scattering of experimental data cited in Sec. 2, this implies that the given estimate may be characterized by other values under different experimental conditions.Most importantly, it can be caused by using small V-factor fiber lasers with complex and temporally unstable near-field LID shapes.

Conclusions
The performance of a laser cutting system is determined not only by the speed-power-pressure concept, but also by the focused beam spatial characteristics.In this paper, we demonstrate that a wave SA of few tenths of wavelength is sufficient for considerably altering the beam characteristics, and elucidate how an aberrational beam distortion effects on laser cutting quality.All aberrated beams produce a strongly asymmetric caustic shape.This offers no possibility of using Gaussian beam propagation for energy estimates of cutting performance and links the laser-melt pool interaction mode to an actual laser intensity distribution on the kerf front.A caustic asymmetry significantly contributes to asymmetrization of both the layout of FPPs pair relative to the BFS and to the possible range of FPP variation.The particularity of the optical design of transport systems for powerful laser cutters suggests using dominantly positive lenses and optical materials with positive ∂n∕∂T.It leads to the situation in which the design and manufacturing SA as well as TIA amplify each other and cause an appreciable beam distortion, especially in the cases of low-M 2 lasers.This first results in the BFS moving toward the transport system.An ALR is always characterized by a truncated beam shape with strong transverse modulation; likewise, by strong axial variation of maximal intensity, as schematically shown in Fig. 7(a).This narrows the range of negative FPPs considerably, making the system much more sensitive to deviations from the best position.In contrast, a smooth appearance of a beam variation outside ALR significantly extends the range of the positive FPPs.
Negative wave SAs are unlikely in typical transport systems of gas-assisted cutters.Yet, they can appear in the case of a marked number of negative lenses.Examples of such systems are telephoto and reversed telephoto objectives focusing the radiation of solid-state and gas lasers and integrated finite optical systems for fiber lasers as well as F-Theta lenses for remote cutting machines.Our results show that in this case the character of the process is inversely changed, as shown in Fig. 7(b).The only difference is that the BFS position and caustic shape are caused by the competition of SA and TIA of opposite signs.Rigorous experimental verification of preliminary results for negative wave SAs' effect on cutting quality will be the direction for future work.

2 .Fig. 1 Fig. 2
Fig.1Best focal spot (BFS) and paraxial focus plane (PFP) locations for: (a) the ideal beam, (b) a beam with positive, and (c) negative wave spherical aberrations (SAs).RC is the ray congruence region, aberration location range (ALR) is the aberration location range, z R is the Rayleigh range.Blue highlighted area denotes the possible range of focal point position (FPP) used in gas-assisted laser cutting of medium-section metals.

Fig. 3
Fig. 3 Computed rotationally symmetric laser intensity distribution (LID) trains for: (a) Gaussian source 2ω S ¼ 7 μm, negligible SA, (b) the same source in the presence of 0.08λ WA rms positive wave SA, (c) uniform source 2ω S ¼ 200 μm, negligible SA, (d) the same source in the presence of 0.8λ WA rms negative wave SA.Object space NA is 0.10 and linear magnification is −1.175 for all figures.

Fig. 4
Fig. 4 Schematic diagram of the experimental setup.

Fig. 5 Fig. 6
Fig.5Cross sections of rotationally symmetric transverse LIDs registered for the array of experimental optical systems.Marked regions show the FPP ranges where excellent and acceptable cuts have been obtained (see also Fig.6).Each LID shape is normalized by its maximal amplitude.

Fig. 7
Fig. 7 Schematic of gas-assisted laser cutting with aberrated beams.(a) Thin red line shows the beam size measured at the 1∕e 2 level, dotted line shows the maximal intensity.(b) The scales for all plots are free.