Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

Double-autoionization decay of resonantly excited single-electron states

Perturbation theory in the lowest non-vanishing order in interelectron interaction has been applied to the theoretical investigation of double-ionization decays of resonantly excited single-electron states. The formulae for the transition probabilities were derived in the LS coupling scheme, and the orbital angular momentum and spin selection rules were obtained. In addition to the formulae, which are exact in this order, three approximate expressions, which correspond to illustrative model mechanisms of the transition, were derived as limiting cases of the exact ones. Numerical results were obtained for the decay of the resonantly excited Kr 1 3d^{-1}5p[^1P] state which demonstrated quite clearly the important role of the interelectron interaction in double-ionization processes. On the other hand, the results obtained show that low-energy electrons can appear in the photoelectron spectrum below the ionization threshold of the 3d shell. As a function of the photon frequency, the yield of these low-energy electrons is strongly amplified by the resonant transition of the 3d electron to 5p (or to other discrete levels), acting as an intermediate state, when the photon frequency approaches that of the transition.

Two-photon decay of the 1s2s ^1S_0 state in _36Kr^34+ produced by resonant transfer and excitation

The doubly excited 2s2p ^1P_1 level of Kr^{34+} populated via resonant transfer and excitation (RTE) feeds selectively the metastable ls2s ^1 S_0 state which can only decay via simultaneous emission of two photons to the ground state 1s^2 ^1 S_0. X-ray/X-ray coincidence measurements in heavy ionatom collisions enable the direct measurement of the spectral distribution of the two-photon decay in He-like ions. In addition, we observe strong photon cascades indueed by radiative electron capture.