Number, time of appearance and thickness of spalled fragments.
The basic mathematical models, computational algorithms, and results of mathematical modeling of various modes of laser action on metals are considered. It is shown that for mathematical description and analysis of the processes of laser heating, melting, and evaporation of condensed media, various theoretical approaches are used: continuum, kinetic, atomistic, etc. Each of them has its own field of applicability, its advantages, and disadvantages. Mathematical description of ns-laser ablation is usually carried out within the framework of continuum approach in the form of hydrodynamic models that take into account reaction of irradiated material to varying density, pressure, and energy both in the target and in the vapor-gas medium. Within the framework of continuum approach, a multiphase, multifront hydrodynamic model and computational algorithm were constructed that were designed for modeling ns-PLA of metal targets embedded in gaseous media. It is shown that proposed model and computational algorithm allow to carry out the simulation of interrelated mechanisms of heterogeneous and homogeneous evaporation of metals manifested as a series of explosive boiling. Modeling has shown that explosive boiling in metals occurs due to the presence of a near-surface temperature maximum. It has been established that in ns-PLA, exposure regimes can be realized in which a phase explosion is the main mechanism of material removal. The verification of reliability of obtained results was carried out by comparing experimental data and calculations with atomistic models.
- nanosecond pulse
- laser action
- hydrodynamic model
- mathematical modeling
- explosive boiling
- phase explosion
- subsurface temperature maximum
The pulsed laser ablation (PLA) of condensed media has been intensively studied over the past few decades [1, 2]. The increased interest in PLA is determined by the increasing possibilities of its use in a variety of applications, beginning with the already traditional ones: microprocessing [3, 4], pulsed laser deposition (PLD) [5, 6], laser-induced breakdown spectroscopy (LIBS) [7, 8], and new rapidly developing areas: production of nanomaterials [9, 10], surface nanostructuring , chemical, and physical synthesis . The total effect of nanotechnology and photonics led to the emergence of a new direction—laser synthesis of colloids , which draws general attention to its extensive applicability, primarily in biomedicine [14, 15].
Numerous applications make PLA an attractive direction for fundamental investigations. Despite extensive studies of fundamental properties of laser ablation performed earlier, a number of important physical phenomena still remain insufficiently well studied and understood. Previous studies have established that the nature of interaction of laser radiation depends both on the modes of action: wavelength [16, 17], duration [18, 19], and laser pulse intensity [20, 21] and on thermophysical and optical properties of the target [22, 23], presence of the surrounding gas  and its pressure [25, 26]. The greatest differences in the physical mechanisms of laser ablation of metals are observed between short (ns) and ultrashort (ps, fs) pulsed modes [18, 19].
In ultrashort range (fs, pc) of influence, laser radiation freely reaches surface of the target. Absorption of laser radiation by a degenerate electron gas followed by a slowed-down energy exchange between electron and phonon components leads to strong deviation from locally thermodynamic equilibrium of the system as a whole. As a result, laser ablation and its accompanying processes develop after the end of the pulse.
Laser ablation in the nanosecond range is a more complex phenomenon involving many interrelated processes both during action and after the end of laser pulse. Such processes include laser target heating, heterogeneous, and homogeneous phase transformations, taking place in the evaporated matter, formation and expansion of plasma plume, heat transfer, laser and intrinsic plasma radiation transfer, generation and propagation of shock waves, and contact boundaries in the vaporized matter and the surrounding gas environment. In contrast to ultrashort regime, in nanosecond range, with a certain choice of parameters of the action, two experimentally observed and explored phenomena arise—volumetric boiling (phase explosion) of liquid phase of the target and formation of laser plasma in vaporized matter and surrounding gas.
Putting an irradiated target into an external gas environment, which is typical for most PLA applications, significantly complicates the situation. In these cases, the long-lived processes of laser-plasma plume and associated generation and propagation of interacting fronts of shock waves and contact boundaries in the vapor-gas medium play an important role in overall picture of laser ablation . The presence of such a large number of interconnected physical processes creates additional difficulties in determination and investigation of the basic mechanisms of ablation. The information obtained by the methods of instrumental diagnostics is not sufficient due to the lack of data on fundamental phenomena associated with rapidly changing energy (thermal and laser radiation), hydrodynamic fields, and the kinetics of heterogeneous and homogeneous phase transformations in the solid and liquid phases, ionization of vaporized matter, and gas environment. At the same time, the understanding of fundamental physics of internal structure of plasma plumes and their spatiotemporal evolution in the process of plasma expansion in the background gas at atmospheric pressure is of decisive importance for many engineering applications. For this reason, nanosecond laser ablation continues to be an area of active research in which mathematical modeling plays an increasing role [28–31].
Any simulation begins with a choice of mathematical model, construction and development of which is given a paramount importance in the computing experiments. For theoretical description and analysis of PLA process of condensed media, various theoretical approaches are used: continuum, kinetic, atomistic (molecular dynamics, etc.). Each of them has its own field of applicability, its advantages, and disadvantages.
Atomistic models allow us to conduct research at the atomic level and obtain fundamental knowledge about structure, thermodynamic, and mechanical properties of crystalline materials [32, 33], about physical mechanisms of various processes [34, 35], including the kinetics of heterogeneous and homogeneous phase transitions [36, 37]. The basic methods of atomistic modeling—molecular dynamics (MD) and Monte Carlo (MC), which use as a rule, semiempirical interaction potentials, operate with tens and hundreds of millions of atoms, and allow calculations in the time range of nanosecond duration.
However, even with the use of parallel computer platforms, computational costs are enormous, and the space-time scales inherent for PLA processes are beyond the limits of accessibility for atomistic modeling methods. Therefore, in spite of constant progress in the field of designing interatomic potentials and increasing the power of computing systems, the final overcoming of computational constraints is hardly achievable, and continuum models will always remain relevant.
Continuous models based on the equations of continuous medium mechanics are realized, as a rule, in the form of hydrodynamic models [18, 24, 28, 30, 31, 38–42] and use the minimum of information and operate with average values of physical characteristics calculated on infinitesimal volume. The methods for solving them are more compact, they have higher accuracy and a relatively small amount of computation. The main shortcomings of continuum approach are manifested in the absence of the possibility of direct investigation of elementary processes in materials and limited possibilities of mathematical description of homogeneous mechanisms of phase transitions of the first kind and calculation of thermophysical, thermodynamic, optical, and other characteristics of matter in a wide range of parameters. These problems are much easier and more fully solved within the framework of atomistic modeling, the results of which can be used as input parameters in meso- and macrolevel models.
In this paper, the application of continuum approach to modeling in preplasma regime of processes dynamics and of main mechanisms of ns-PLA of metal target (Al) in air is considered. The mechanisms of heterogeneous and homogeneous phase transitions interacting with each other are analyzed in detail. Primary attention is paid to model the dynamics of phase explosion of liquid phase of aluminum and the expansion of its fragments in the air, since explosive boiling is considered to be one of the most effective thermal mechanisms of ns laser ablation of materials. Various aspects of this problem have been studied in a number of theoretical and experimental studies [43–56], but there is still no consensus on the mechanism of the phase explosion in metals. In order to obtain detailed information on interaction of heterogeneous and homogeneous mechanisms and data on laser plume morphology, simulation of laser heating, melting, surface evaporation, and evolution of plume in the vapor-gas medium is performed within the framework of new hydrodynamic model with temperature dependences of material properties of the target and explicit tracking of interphase boundary fronts, contact boundary, and shock wave. The release of liquid phase fragments into the atmosphere as a result of phase explosion is modeled by the procedure for introducing a quasinucleus of a new phase (vapor) in the region of near-surface maximum of temperature reaching the value
2. Hydrodynamic model
Laser ablation of Al target placed in the air atmosphere is considered. In preplasma regimes of the action in ablation, the processes in three phase states (solid, liquid, and vapor) and two states (perturbed, unperturbed) of external gaseous medium are taken into account. Laser radiation propagates from right to left. The air for the selected action mode is completely transparent to the laser flow, which is partially reflected from the metal surface and partially absorbed in the near-surface layer of the target.
The mathematical description of ns laser ablation within the continuum approach is realized in the form of hydrodynamic models that allow one to take into account the reaction of condensed (target) and vapor-gas (vaporized matter, gas) media to varying density, pressure, and energy. The processes in each medium are described by a system of nonstationary equations of gas-hydrodynamics supplemented by equations of energy with thermal conductivity, equation of laser radiation transfer, and corresponding equations of state. The 1-D approximation is used for spatial variables. For laser action regimes under consideration, conditions of locally thermodynamic equilibrium (LTE) are assumed for all processes both in the irradiated target and in the vapor-gas medium. Accordingly, hydrodynamic model is formulated in one-temperature approximation.
Phase states and vapor-gas medium are separated among themselves by moving interphase boundaries solid-liquid
In order to obtain complete information on the kinetics of phase transformations, the morphology and dynamics of laser plume, all moving fronts and boundaries in the course of the solution are subjected to explicit tracking using appropriate relationships, which are simultaneously the boundary conditions for hydrodynamic and energy equations.
Complete system of equations is represented as:
2.1. Boundary conditions
In the hydrodynamic models that easily combine with the kinetic ones, the heterogeneous mechanisms of phase transitions: melting/crystallization and evaporation/condensation are described naturally. A distinctive feature of heterogeneous phase transitions is the presence of sharp interfaces at which the main thermophysical and optical characteristics: enthalpy
In the case of rapid phase transformations that typical for pulsed nanosecond action, a complete set of equations for mass, momentum, and energy flows is used to describe heterogeneous phase transitions. However, considering that the absence of phase equilibrium at interphase boundary is a necessary condition for phase transition to occur, the expressions for conservation laws must be supplemented by appropriate kinetic relations characterizing the degree of nonequilibrium of phase transition.
2.1.1. Left fixed boundary,
x = Гs
The condition that flow of mass and heat be equal to zero is used as boundary conditions on the left fixed boundary:
2.1.2. Model of heterogeneous (surface) melting: crystallization, õ =
For fast phase transitions, heterogeneous melting model consists of a system of equations expressing three conservation laws: mass, momentum, and energy supplemented by kinetic condition for the velocity of melting front
where is nonequilibrium heat of melting,
2.1.3. Model of heterogeneous evaporation
x = Гlv
Investigation of heterogeneous evaporation process began already in the century before last from experimental work: Hertz  and theoretical work: Knudsen  and continues at the present time, which is determined by its practical importance and not fully clarified features of nonequilibrium behavior of the substance during its evaporation. The mechanism of heterogeneous laser evaporation is realized in subcritical region of surface temperature
Heterogeneous evaporation is characterized by a sharp phase boundary. Thin (several mean path lengths) nonequilibrium Knudsen layer (KL) is adjacent to the boundary surface. The nonequilibrium of the KL is determined by the flow of matter through the phase boundary. The total flow consists of flow of particles (atoms, molecules) leaving the condensed phase and flow of particles (atoms of molecules) returning from the evaporated matter as a result of collisions. These flows have different distribution functions,
In the phase equilibrium state, when the pressure of saturated vapor
where . With these expressions taken into account, energy conservation law can be represented as:
where is the nonequilibrium heat of evaporation,
To determine parameters
the value m is determined from the equation
2.1.4. Model of surface condensation
When the inequality
When passing through the point
Thus, the description of the kinetics of heterogeneous evaporation is carried out by two one-parameter dependencies used as boundary conditions (for a known parameter
2.1.5. Model of volumetric boiling of liquid phase heated by a laser pulse
The greatest difficulty in the continuum approach is the description of homogeneous mechanisms of phase transformations: melting crystallization and evaporation. Homogeneous mechanisms of phase transformations are characterized by nucleation of a new phase in a certain volume of superheated/supercooled matter. Representing them in continuum hydrodynamic models requires considerable additional efforts [57, 76] associated with formation of a cavity filled with vapor within a condensed medium.
The simplest scheme for simulating volumetric boiling in one-dimensional approximation can be represented by introducing, when certain criteria are satisfied into superheated liquid phase of artificial quasi-nuclei with thickness
The values of limit of superheat temperature depend on the rate of energy input and can be determined in advance from the molecular dynamics simulation. For Al, depending on the rate of heating, limit of superheat temperature of liquid phase
It should be noted that the use of gas dynamic equations inside the cavities which dimensions are small in the initial moments of time in comparison with the mean free path is a simplifying approximation that makes it easy to take into account the influence of expansion velocity of the cavity.
2.1.6. Moving contact boundary,
x = Гvg( t)
On contact boundary and front of the shock wave, well-known standard relations are used . At vapor-air interface, the boundary conditions were set in the form of equal values of velocity, pressure, and temperature:
2.1.7. Moving shock wave,
x = Г sh, g( t)
A shock wave in air
The indices 0 and 1 denote the values of the quantities on side of the background and the shock wave, respectively.
2.1.8. Right moving boundary,
x = Гg( t)
The right boundary on the side of unperturbed gas is declared moving in order to improve economic effectiveness of the computational algorithm . The speed of its motion is found from the differential equation of momentum.
2.2. Computational algorithm
The differential model (1)–(8) was approximated by a family of conservative finite-difference schemes [80, 81] written on computational grids with dynamic adaptation [79, 82, 83]. The method is based on the idea of transition to an arbitrary nonstationary coordinate system that allows calculations with an arbitrary number of discontinuous solutions, such as shock waves, propagating phase and temperature fronts, contact boundaries, and spalled fragments.
3. Temperature dependences of aluminum properties
In the mathematical modeling of the process of laser action on metals with an absorption coefficient of
With the advent of lasers capable of evaporating metals and producing plasma on their surfaces, it has become possible to observe laser-induced phase transitions metal-dielectric with the formation of transparency waves in the massive targets  and thin metallic films  in the subcritical temperature range. In other experimental and theoretical studies of the interaction of ns-laser pulses with an intensity of 107–108 W/cm2 with metallic targets, the results of the appearance of metal-dielectric transitions accompanied by the formation of transparency waves are reported [53, 54, 85, 88–91].
Nevertheless, the physical mechanisms of interaction of laser radiation with metals taking into account metal-insulator transition have not yet been fully studied. This causes great difficulties in determining thermophysical and optical characteristics of metals in the vicinity of the critical region.
Figure 2 shows temperature dependences of thermophysical
4. Results discussion
One of the purposes of this work is a detailed study of mechanism of explosive boiling in metals since explosive boiling is considered to be the most efficient thermal mechanism for the laser ablation of materials. At the same time, the difficulties associated with understanding of the mechanism of homogeneous phase transitions in metals occurring under the action of ns-laser pulses are known. In , based on analytical solution of thermal model, it was established that laser removal of material from a solid target, with a certain choice of irradiation parameters and material, is determined by the presence of near-surface temperature maximum. The validity of this position was confirmed by the results of numerical solution of thermal model for low-absorbing liquids irradiated by laser pulses  and nonmetallic solid materials . For strongly absorbing media mainly metallic, calculations based on thermal model [44, 45, 47, 56] have shown that the magnitude of near-surface temperature maximum is several degrees. On this basis, overheating was excluded from consideration, up to the statement  that in metals, maximum temperature is always on the surface of the target, and sub-surface superheating is impossible. Explosive boiling was interpreted as a surface spitting of liquid phase when a critical temperature is reached on the surface . This interpretation is not convincing, since it does not allow to determine even approximately the parameters of explosive boiling.
4.1. Mode of exposure
Let us consider ns-PLA process of aluminum target under the mode of exposure close to experimental conditions . Al target with the thickness of 10−4 m is placed in the air under the normal conditions with pressure 1 bar and room temperature (
4.2. Evolution of processes
At the initial stage, time evolution of the processes at the surface of the target, in the target and in the gas medium near the target such as appearance of phase fronts (melting, evaporation), contact boundary, and shock wave and associated with formation of new phase media (liquid, vapor) occurs at the leading edge of laser pulse. Heterogeneous melting begins at the moment
By the moment
4.3. Formation of near-surface temperature maximum
A consistent study of laser ablation is complicated by the fact that heterogeneous and homogeneous mechanisms of both melting and evaporation prove to be interrelated, and this interaction must be taken into account explicitly. Pulsed laser action on materials (including metals) has volumetric nature of energy release. Thus, when target is heated strongly in the vicinity of the critical point (
Figure 4a shows a fragment of the spatial temperature distribution in near-surface layer of the target with overheating
4.4. Explosive boiling
The process of explosive boiling begins at the backside front of the laser pulse. First boiling occurs at the moment
|Fragment number, n||1||2||3||4||5|
|Moment of occurrence, ns||1.5||1.7||1.9||2.1||2.4|
|Fragment thickness d, nm||30||37||26||21||20|
Figure 6a – d shows the state of all liquid fragments of explosive boiling at the time
The simulation results are in good qualitative agreement with the experimental data . Quantitatively, the total depth due to release of superheated liquid phase (
Within the framework of continuum approach, a multiphase, multifront hydrodynamic model and computational algorithm were constructed that were designed for modeling ns-PLA of metal targets embedded in gaseous media. The model contains temperature dependences of optical, thermal, and thermodynamic properties of metal target (Al). The temperature dependences of heat capacity
The computational algorithm is based on the method of dynamic adaptation, which ideally fit to the study of problems with heterogeneous phase transitions and allows the production of numerical solutions with an explicit allocation of unlimited number of interface boundaries in condensed medium, contact boundaries, and shock waves in a gaseous medium.
The calculations have shown that under the influence of temperature dependences of optical and thermophysical properties of Al the nature of laser heating of the target in high-temperature region changes substantially. The energy release has a volume nature, and irradiated surface is markedly cooled by the process of surface evaporation. An inhomogeneity arises in the spatial temperature profile—a near-surface maximum of temperature. When the limiting superheating temperature is reached, the conditions for explosive boiling were realized with formation of a cavity in the region of the maximum, where the maximum penetration into the region of superheated metastable liquid near the critical point of liquid-vapor transition is achieved. For Al in action modes under the consideration, the typical depths of temperature maximum are
An approach was proposed to formulation in continuum models of homogeneous nucleation in a superheated liquid phase based on the introduction into hydrodynamic model of quasi-nuclei with thickness of 1–5 nm. Generation of quasi-nuclei is carried out by the criterion of maximum permissible overheating of initial region. The criterion of limit overheating is determined from molecular dynamics simulation.
Mathematical modeling using the developed technique allowed to obtain a sequence of five explosive boilings for the threshold fluence
In metals, as in nonmetals, the transition from surface evaporation to volume removal of mass occurs in the region of temperature maximum near the critical temperature, where the influence of overheated metastable liquid phase determines the competition between surface evaporation and explosive boiling.
It should be noted that in explosive boiling regimes under investigation, cavity formation takes place in the region of positive values of pressure generated by homogeneous nucleation in superheated liquid phase, in contrast to the case of a series of spalles  arising in the region of negative pressures in unloading wave upon action to metal targets of ultrashort pico-femtosecond laser pulses.
The research was funded by the Russian Scientific Foundation (grant No 15-11-00032).