The authors present a concurrent Monte Carlo (MC)-molecular dynamics (MD) approach to modeling matter response to excitation of its electronic system at nanometric scales. The two methods are combined on-the-fly at each time step in one code, TREKIS-4. The MC model describes the arrival of irradiation (a photon, an electron, or a fast ion). It traces induced cascades of secondary electrons and holes, and their energy exchange with atoms due to scattering. The excited atomic system is simulated with an MD model. An efficient way is proposed to account for nonthermal effects in the electron-atom energy transfer in covalent materials via the conversion of the potential energy of the electronic ensemble into the kinetic energy of atoms. Such a combined MC-MD approach enables a time-resolved tracing of the excitation kinetics of both, the electronic and atomic systems, and their simultaneous response to a deposited dose. As a proof-of-principle, it is shown that the proposed method describes atomic dynamics after X-ray irradiation in good agreement with tight-binding MD. The model also allows gaining insights into the atomic system behavior during the energy deposition from a nonequilibrium electronic system excited by an ion impact.

Ultrafast laser irradiation of metals can often be described theoretically with the two-temperature model. The energy exchange between the excited electronic system and the atomic one is governed by the electron–phonon coupling parameter. The electron–phonon coupling depends on both, the electronic and the atomic temperature. We analyze the effect of the dependence of the electron–phonon coupling parameter on the atomic temperature in ruthenium, gold, and palladium. It is shown that the dependence on the atomic temperature induces nonlinear behavior, in which a higher initial electronic temperature leads to faster electron–phonon equilibration. Analysis of the experimental measurements of the transient thermoreflectance of the laser-irradiated ruthenium thin film allows us to draw some, albeit indirect, conclusions about the limits of the applicability of the different coupling parametrizations.

We performed the experimental and theoretical study of the heating and damaging of ruthenium thin films induced by femtosecond laser irradiation. Results of an optical pump-probe thermoreflectance experiment with rotating sample allowing to significantly reduce heat accumulation in irradiated spot are presented. We show the evolution of surface morphology from growth of a heat-induced oxide layer at low and intermediate laser fluences to cracking and grooving at high fluences. Theoretical analysis of pumpprobe signal allows us to relate behavior of hot electrons in ruthenium to the Fermi smearing mechanism. The analysis of heating is performed with the two-temperature modeling and molecular dynamics simulation, results of which demonstrate that the calculated melting threshold is higher than experimental damage threshold. We attribute it to heat-induced surface stresses leading to cracking which accumulates to more severe damage morphology. Our results provide an upper limit for operational conditions for ruthenium optics and also direct to further studies of the Fermi smearing mechanism in other transition metals.

No abstract

Since a few breakthroughs in the fundamental understanding of the effects of swift heavy ions (SHIs) decelerating in the electronic stopping regime in the matter have been achieved in the last decade, it motivated us to review the state-of-the-art approaches in the modeling of SHI effects. The SHI track kinetics occurs via several well-separated stages and spans many orders of magnitude in time: from attoseconds in ion-impact ionization depositing an extreme amount of energy in a target to femtoseconds of electron transport and hole cascades, to picoseconds of lattice excitation and response, to nanoseconds of atomic relaxation, and even longer times of the final macroscopic reaction. Each stage requires its own approaches for quantitative description. We discuss that understanding the links between the stages makes it possible to describe the entire track kinetics within a hybrid multiscale model without fitting procedures. The review focuses on the underlying physical mechanisms of each process, the dominant effects they produce, and the limitations of the existing approaches, as well as various numerical techniques implementing these models. It provides an overview of the ab initio-based modeling of the evolution of the electronic properties, Monte Carlo simulations of nonequilibrium electronic transport, molecular dynamics modeling of atomic reaction including phase transformations and damage on the surface and in the bulk, kinetic Mote Carlo of atomic defect kinetics, and finite-difference methods of track interaction with chemical solvents describing etching kinetics. We outline the modern methods that couple these approaches into multiscale and combined multidisciplinary models and point to their bottlenecks, strengths, and weaknesses. The analysis is accompanied by examples of important results, improving the understanding of track formation in various materials. Summarizing the most recent advances in the field of the track formation process, the review delivers a comprehensive picture and detailed understanding of the phenomenon. Important future directions of research and model development are also outlined.