Abstract:The interactions of femtosecond lasers with gold targets were investigated with a numerical method combining molecular dynamics (MD) and the two-temperature model (TTM). Previous works using MD-TTM method did not consider all the thermodynamic parameters and the interatomic potential dependent of the electron temperature simultaneously. Therefore, we developed a LAMMPS function to achieve this. To accurately capture the physics behind the interactions, we also included the electron blast force from free electron pressure and the modified Fourier law with steep electron temperature gradient in our model. For bulk materials, a stress non-reflecting and heat conducting boundary is added between the atomistic and the continuum parts. The modified boundary force in our study greatly reduces the reflectivity of the atomistic-continuum boundary compared with its original form. Our model is the first to consider all these factors simultaneously and manage to predict four femtosecond laser ablation phenomena observed in the experiments. In this dissertation, the thermodynamic parameters in the two-temperature model were extensively explored. We considered three different approaches to calculate these parameters: namely interpolation, ab initio calculation, and analytical expression. We found that simple interpolation between solid state and plasma state could lead to high level of inaccuracy, especially for electron thermal conductivity. Therefore, ab initio calculation and analytical expression were used for the calculation of the thermodynamic parameters in more advanced studies. The effects of electron thermal conductivity and electron-phonon coupling factor on electron and lattice temperatures were analyzed.Our studies considered electron temperature dependent (ETD) and electron temperature independent (ETI) interatomic potentials. The ETI interatomic potential is easier to implement and therefore it is used in our phase change study to investigate the effects of target thickness on melting. Homogeneous melting occurred for thin films, while melting can be observed through the movement of the solid-liquid interface in thick or bulk materials. However, the ETI potential overestimated the bond strength at high temperatures. Therefore, ablation process was studied with the ETD potential. Three ablation mechanisms were found in our simulations at different laser fluences. Short nonthermal ablation was only observed at the ablation threshold. With increasing laser fluence, spallation was then seen. In high laser fluence regime, phase explosion occurred on the surface and coexisted with spallation.Lastly, we researched on the effects of the delay time between two femtosecond laser pulses. Various delay times did not have much influence on melting depth. In low laser fluence regime, with increasing delay time, the target went through nonthermal ablation, to spallation and to no ablation. In high laser fluence regime, longer delay time encouraged phase explosion while suppressed spallation.

Calculation of Argon-Aluminum Interatomic Potential and Its Application in Molecular Dynamics Simulation of Femtosecond Laser Ablation

Molecular Dynamics Simulation Of Ultra-Fast Thermal Transition In Thin Foil Irradiated By High Power Ultra-Short Laser

Atomistic simulations of ultra-short pulse laser ablation of aluminum: validity of the Lambert-Beer law

Development of the temperature-dependent interatomic potential for molecular dynamics simulation of metal irradiated with an ultrashort pulse laser

Coupling of Material Point Method and Molecular Dynamics for Modeling Ultrafast Laser Interaction with Metals

Numerical simulation study of ultraviolet nanosecond laser ablation of Cu-Al2O3 dispersion strengthened alloys

Mechanisms of femtosecond laser ablation of Ni3Al: Molecular dynamics study

Study on the Interaction Between Nanosecond Laser and 6061 Aluminum Alloy Considering Temperature Dependence

Kinetic simulations of laser-induced plume expansion from a copper target into a vacuum or argon background gas based on ab initio calculation of Cu–Cu, Ar–Ar, and Ar–Cu interactions

Simulation of expansion of aluminum plasmas produced by a small focal spot nanosecond laser irradiation

Hybrid atomistic-macroscale modeling of long-time material behavior in nanosecond laser-material interaction

Kinetic Evolution of Laser Ablating Alloy Materials

Efficient modeling of metal ablation irradiated by femtosecond laser via simplified two-temperature model coupling molecular dynamics

Influence of Laser-Induced Air Breakdown on Femtosecond Laser Ablation of Aluminum.

A fractional dual-phase-lag generalized thermoelastic model of ultrashort pulse laser ablation with variable thermal material properties, vaporization and plasma shielding

Laser Ablation of Al-Ni Alloys and Al-Ni Layer Systems simulated with Molecular Dynamics and the Two-Temperature Model

PHASE CHANGE AND ABLATION STUDY OF METALS BY FEMTOSECOND LASER IRRADIATION USING HYBRID TTM/MD SIMULATIONS

Implementation of kinetics of phase transitions into hydrocode for simulation of laser ablation

Numerical Simulation for Influence of Pulse Duration on Femtosecond Laser Ablation of Alloy

Nanoparticles Formed in Picosecond Laser Argon Crystal Interaction

Accurate Path-Integral Molecular Dynamics Calculation of Aluminum with Improved Empirical Ionic Potentials