物理参数变化对短脉冲激光激励温度场的影响

 为研究多物理参数（耦合系数、电子热导率、电子热容、晶格热容）同时随温度变化对短脉冲激光辐照金属材料产生温度场分布的影响,基于双温耦合理论,建立了短脉冲激光辐照金属材料金的加热过程的有限元求解模型。在同时考虑脉冲激光的空间、时间分布和多参数同时随温度变化的情况下,得到短脉冲激光辐照金属材料金激励产生的温度场二维瞬态分布,并进一步比较了多物理参数同时随温度变化和采用室温物理参数两种情况下温度场分布的区别。数值结果表明：多物理参数同时随温度变化使电子温度和晶格温度的上升变快,最大值变大,而且使得材料中激光穿透直接辐照到的区域温度变高。     

Femtosecond laser ablation of aluminium

We investigate femtosecond laser ablation of aluminium using a hybrid simulation scheme. Two equations are solved simultaneously: one for the electronic system, which accounts for laser energy absorption and heat conduction, the other for the dynamics of the lattice where the ablation process takes place. For the electron-temperature a generalized heat-conduction equation is solved by applying a finite difference scheme. For the lattice properties, e.g. pressure, density or temperature, we use common molecular dynamics. Energy transfer between the subsystems is allowed by introducing an electron-phonon coupling term. This combined treatment of the electronic and atomic systems is an extension of the well known two-temperature model [Anisimov, Kapeliovich, Perel’man, Electron emission from metal surfaces exposed to ultra short laser pulses, JETP Lett. 39 (2)].     

Multi-material two-temperature model for simulation of ultra-short laser ablation

We investigate the interaction of 100 fs laser pulses with metal targets at moderate intensities (10¹² to 5 × 10¹³ W/cm²). To take into account effects of laser energy absorption and relaxation we develop a multi-material two-temperature model based on a combination of different approaches. The backbone of the numerical model is a high-order multi-material Godunov method in a purely Eulerian form. This formulation includes an interface-tracking algorithm and treats spallation at high strain rates and negative pressures. The model consistently describes the hydrodynamic motion of a two-temperature plasma and accounts for laser energy absorption, electron-phonon/ions coupling and electron heat conductivity. In particular, phase transitions are accurately taken into account by means of a wide-range two-temperature multi-phase equation of state in a tabular form. The dynamics of the phase transitions and the evolution of the heat-affected zone are modeled and analyzed. We have found that a careful treatment of the transport coefficients, as well as consideration of phase transitions is of a great importance in obtaining reliable numerical results. Calculation results are furthermore compared for two metals with different electron-phonon coupling parameters (Au and Al). We have found that the main part of ablated material results from fragmentation of melted phase caused by tensile stresses. A homogeneous nucleation mechanism alone does not explain experimentally observed ablation depth.     

Temperature dependences of the electron-phonon coupling, electron heat capacity and thermal conductivity in Ni under femtosecond laser irradiation

The electron temperature dependences of the electron-phonon coupling factor, electron heat capacity and thermal conductivity are investigated for Ni in a range of temperatures typically realized in femtosecond laser material processing applications, from room temperature up to temperatures of the order of 10⁴ K. The analysis is based on the electronic density of states obtained through the electronic structure calculations. Thermal excitation of d band electrons is found to result in a significant decrease in the strength of the electron-phonon coupling, as well as large deviations of the electron heat capacity and the electron thermal conductivity from the commonly used linear temperature dependences on the electron temperature. Results of the simulations performed with the two-temperature model demonstrate that the temperature dependence of the thermophysical parameters accounting for the thermal excitation of d band electrons leads to higher maximum lattice and electron temperatures achieved at the surface of an irradiated Ni target and brings the threshold fluences for surface melting closer to the experimentally measured values as compared to the predictions obtained with commonly used approximations of the thermophysical parameters.     

Direct simulation monte carlo of pulsed laser ablation of metals with clusterization processes in vapor

A mathematical model is proposed to describe laser ablation of metals in vacuum under the action of nanosecond laser pulses of moderate intensity taking into account the processes of cluster formation and decay in the vapor cloud. To describe the laser radiation absorption and metal heating, the thermal model based on the unsteady one-dimensional heat equation with a volume heat source is used, and the method of statistical modelling is employed for modelling vapor expansion and the processes of cluster formation. The efficiency of the proposed complex model is considered by the example of pulsed laser ablation of a niobium target. The work was supported financially by INTAS (grant No. 03-51-5208).     

Effect of phonon spectrum heating on the formation of laser-induced electron plasma in wide-gap insulators

The effect of lattice heating by laser pulses on the dynamics of electron plasma generation in transparent solids has been theoretically studied. Several ways of taking into account the contribution of the phonon spectrum heating to the electron avalanche dynamics, depending on the type of the effective (with respect to the field energy transfer to electrons) phonons and laser pulse duration, have been proposed. A comparative analysis of the results of Monte Carlo computation of electron gas heating in the laser pulse field, which were obtained for cold and heated lattices, has been performed. It is shown that the consideration of the effect of lattice heating on the probabilities of electron-phonon and electron-phonon-photon scattering leads to an increase in the avalanche rate, which is more pronounced at longer wavelengths of the incident radiation and under longer laser pulses. Some qualitative features of the redistribution of the energy, absorbed during a pulse, between the electron plasma and lattice are revealed, which suggest initiation of irreversible microscopic changes in the insulator. In particular, the ratio R of the energy accumulated in the electron subsystem to the excess (with respect to the initial equilibrium state) energy in the phonon subsystem has been calculated for different initial lattice temperatures. It is shown that this ratio increases with a decrease in the laser wavelength in the computation scheme with lattice heating disregarded and decreases at all pulse durations when the lattice heating is taken into account.     

Thermal Analysis of Gold Nanorods Heated with Femtosecond Laser Pulses

We present an axisymmetric computational model to study the heating processes of gold nanoparticles, specifically nanorods, in aqueous medium by femtosecond laser pulses. We use a two-temperature model for the particle, a heat diffusion equation for the surrounding water to describe the heat transfer processes occurring in the system, and a thermal interface conductance to describe the coupling efficiency at the particle/water interface. We investigate the characteristic time scales of various fundamental processes, including lattice heating and thermal equilibration at the particle/surroundings interface, the effects of multiple laser pulses, and the influence of nanorod orientation relative to the beam polarization on energy absorption. Our results indicate that the thermal equilibration at the particle/water interface takes approximately 500 ps, while the electron-lattice coupling is achieved at approximately 50 ps when a 48×14 nm gold nanorod is heated to a maximum temperature of 1270 K with the application of a laser pulse having 4.70 J/m(2) average fluence. Irradiation by multiple pulses arriving at 12.5 ns time intervals (80 MHz repetition rate) causes a temperature increase of no more than 3 degrees during the first few pulses with no substantial changes during the subsequent pulses. We also analyze the degree of the nanorods' heating as a function of their orientation with respect to the polarization of the incident light. Lastly, it is shown that the temperature change of a nanorod can be modeled using its volume equivalent sphere for femtosecond laser heating within 5-15% accuracy.     

Melting and resolidification of gold film irradiated by nano- to femtosecond lasers

Rapid melting and resolidification of a free-standing gold film subject to nano- to femtosecond laser pulses are investigated using the two-temperature model in conjunction with an interfacial tracking method. The interfacial velocity, as well as elevated melting temperature and depressed solidification temperature, in the ultra-fast phase-change process are obtained by considering the interfacial energy balance and nucleation dynamics. A nonlinear electron heat capacity and a temperature-dependent electron–lattice coupling factor for the rapid phase change are taken into account. Effects of laser pulse width and fluence on melting and resolidification are also studied. PACS 42.62.Eh; 63.20.Kr; 64.70.Dv     

Effects of particle size and laser wavelength on heating of silver nanoparticles under laser irradiation in liquid

Laser energy absorption results in significant heating of metallic nanoparticles and controlling the heating of nanoparticles is one of the essential stages of selective cell targeting. It is necessary to note that the laser action should be done by laser pulses with a wavelength that is strongly absorbed by the particles and it is important to select wavelengths that are not absorbed by the medium. Laser pulse duration must be chosen sufficiently short to minimize heat flow emitted from absorbing particles. Numerical calculations based on Mie theory were used to obtain the effect of laser wavelength and particle size on absorption factor for colloidal silver nanoparticles with radii between 5 and 50 nm. Calculations for acquiring temperatures under irradiations of pulsed KrF laser and pulsed Nd:YAG laser were performed. We showed that for low wavelengths of the laser, smaller nanoparticles have larger absorption efficiency compared to larger nanoparticles and in high wavelengths, temperature of all particles increased in the same way.     

Theoretical investigations of the processes of selective laser interaction with melanin granules in pigmented tissues for laser applications in medicine

Theoretical investigations and the results of computer modeling of optical, thermophysical, thermochemical, and hydrodynamical processes during selective laser interaction with melanoprotein granules (melanosomes) in heterogeneous pigmented tissues (retinal pigment epithelium) are reviewed in this paper. Physico-mathematical models and system of equations are formulated which describe interaction processes for “short” laser pulses of duration t _p < 10^?6 s and for “long” pulses of duration t _p > 10^?6 s. The results of numerical simulation of the processes give the space-time distributions of temperature and degrees of thermodenaturation of the protein molecules inside and around melanosomes and in the volume of irradiated tissue. Energy absorption, heat transfer, and thermochemical processes occurring during the interaction of laser pulses with pigmented spherical and spheroidal granules in heterogeneous tissues are theoretically investigated. The possibility for selective interaction of short laser pulses with pigmented granules, which results in the formation of denaturation microregions inside and near the pigmented granules (granular thermodenaturation) without origination of a continuous macroscopic thermodenaturation lesion in tissue, is discussed. An analytical model of heating of a single spherical and spheroidal granule by a laser pulse is presented. Simple equations for the time dependences of particle temperature are obtained. Vapor generation under the action of a laser pulse on pigmented spherical granules in a water-containing tissue and the formation and dynamics of a vapor blanket are theoretically investigated. The values of pulse energy which give rise to granular and ophthalmoscopically visible thermodenaturation lesions on the retina and to vapor generation are discussed, as well as laser-induced breakdown on granules in pigmented tissues, on the basis of experimental results and numerical and analytical calculations. The comparison and agreement of the numerical results with the experimental data validate the models and techniques developed. The presented results are of essential interest for laser applications in ophthalmology and can be used to investigate laser interaction with heterogeneous tissues in dermatology and various fields of laser medicine.     

Thermal description of laser annealing

Heating and cooling rates as well as melting and solidification velocity of surface layers of irradiated samples by laser pulses are semi-quantitatively described in terms of heat flow concepts based on the assumption that laser light is directly converted in lattice heating. The range of validity of this approach compared with a more complete scheme taking into account the free carrier plasma evolution is sketched and the importance of Auger effect in the plasma-lattice coupling mechanism is detailed. The most important consequences of the quencing rates achievable by short laser pulse irradiation on the structure modification of semiconductor surface layers are reviewed with more details on the liquid to amorphous silicon transition. This is in fact the more new and less understood fast solidification process induced by pulsed laser irradiation.     

Nonlinear absorption,scattering, and extinction of laser radiation by two-layered spherical system-gold nanoparticle and vapor shell in water

Nonlinear absorption, scattering and extinction of laser radiation with wavelengths 532, 633 nm by spherical gold nanoparticles (NPs) with radii in the range of 5–100 nm placed in water and heated by laser radiation with formation and expansion of vapor nanoshells is theoretically investigated. Decrease of absorption, decrease and subsequent increase of scattering and extinction with increasing of shell radius beginning from the initial period of shell expansion is established. Optical indicatrixes and nonlinear behavior of scattered radiation are investigated including the examination of these characteristics during the adiabatic expansion of vapor shell. Formation of vapor nanoshells (bubbles) as a result of the action of short laser pulses on NPs placed in tissue was proposed for cutting of tissue.     

Short-pulse laser ablation of solids: from phase explosion to fragmentation

The mechanisms of laser ablation in silicon are investigated close to the threshold energy for pulse durations of 500 fs and 50 ps. This is achieved using a unique model coupling carrier and atom dynamics within a unified Monte Carlo and molecular-dynamics scheme. Under femtosecond laser irradiation, isochoric heating and rapid adiabatic expansion of the material provide a natural pathway to phase explosion. This is not observed under slower, nonadiabatic cooling with picosecond pulses where fragmentation of the hot metallic fluid is the only relevant ablation mechanism.     

Transient effects in pulsed laser irradiation

Theoretical analysis of the influence of the temporal profile (rectangular, triangular, Gaussian) of the laser pulse on heating/cooling and phase transition velocities and quantity of ablated material was performed on the basis of a multifront Stephan problem. Modeling showed that material removal under stationary conditions (that correspond to long pulses) is entirely controlled by specific heat and material density, while in the case of transient regimes (short pulses) thermal conductivity and heat capacity play a predominant role. Interaction of the melting and evaporation fronts characterized by an evaporation front velocity far exceeding the melting front one is one of the examples of the transient nature of the phenomena influenced by the laser pulse parameters.     

Nanotexturing of gold films by femtosecond laser-induced melt dynamics

Using tightly focused femtosecond laser pulses one can produce sub-micrometer holes and periodic structures in almost all solid materials. But with structure sizes below one micron, femtosecond laser structuring forfeits its universal character. In this paper we concentrate on the discussion of the influence of electron-phonon coupling strength in metal ablation morphology. In metals with weak electron-phonon coupling, results of femtosecond laser ablation can be dominated by solidified flow dynamics in molten material. One can make use of this phenomenon for ablation-free nanotexturing of gold films. PACS  42.62.-b; 42.65.Re; 52.38.hf     

大口径谐振腔式固体激光中的动力学不稳定性研究

大口径谐振腔式固体激光定标到足够高功率后,由于非稳定腔内激光不对称振荡导致介质上激光强度横向不均匀,激光提取与介质产热耦合会引起激光输出功率和光束质量在时域上的不稳定性。针对三种激光器构型:Nd:YAG薄片、Yb:YAG薄片和浸入式液冷叠片激光器,通过理论分析与数值模拟揭示了其各自不同的光热耦合机制及其影响因素,给出了激光输出的瞬态演化规律。结果表明,前两种构型中基于激光动力学的光热耦合具有饱和效应,其动力学不稳定性只出现在特定参数区间,可以通过恰当的设计加以避免;第三种构型中动力学不稳定性表现出明显的阈值特征,在强光状态下只能通过减少冷却液吸收系数来抑制。     

Effects of the ZnSe concentration on the structural and optical properties of ZnSe/PVA nanocomposite thin film

This study investigated the effects of ZnSe nanoparticles (NPs) on the structural and (linear and nonlinear) optical properties of polyvinyl alcohol (PVA) thin film. Three samples of ZnSe NP-doped PVA thin films with different concentrations of ZnSe were produced on a glass substrate. The ZnSe NPs were synthesized by pulsed laser ablation of the ZnSe bulk target immersed in distilled water using a 1064 nm wavelength and a high frequency pulsed Nd:YAG laser. The optical bandgap energies of the films were extracted from their UV-Vis-NIR absorption spectra. The corresponding energy bandgaps of the nanocomposite films declined as the ZnSe NPs doping concentration increased. X-ray diffraction analysis was used to characterize the crystalline phases of the ZnSe/PVA nanocomposite films. The concentration-dependent nonlinear optical absorption and nonlinear refraction behaviors of the films after exposure to 532-nm nanosecond laser pulses were investigated using the Z-scan technique. The nonlinear absorption response of the films was positive when measured using an open aperture scheme, which was attributed to the two-photon absorption mechanism. In addition, the nonlinear refraction indices had a negative value and they increased as the concentration of ZnSe NPs in the films increased.     

Mathematical models of microwave heating of a coal mass with release of absorbed energy by the heat radiation law

Innovation development of coal power engineering, besides qualitative improvement of traditional technologies, should be based also on invoking novel methods and means both in preparation and efficient combustion of coal and for reduction of toxic emissions. The most activating action on coal fuel is provided by microwave radiation. A microwave field intensifies energy-consuming processes such as drying of coal, its dispersion, removal of impurities, qualitative heating, combustion with improved ecological characteristics, etc. Base information for development of microwave technologies in coal power engineering is a temperature field. In this paper we study approximately analytically a mathematical model of flat coal mass heating when the heat release source, due to absorption of microwave energy, obeys the Bouguer law. A specific feature of this model is the presence of nonlinearity caused by dissipation of energy from the mass surface due to heat radiation. The obtained dependences on the temperature field are the base for calculation of ignition and burnout of coal substance, search for the optimal microwave regimes, estimation of thermal destruction parameters, etc. The present paper is a continuation of authors’ publication in this journal [1].