Spin transport and spin dynamics after femtosecond laser pulse irradiation of iron (Fe) are studied using a kinetic Monte Carlo model. This model simulates spin dependent dynamics by taking into account two interaction processes during nonequilibrium: Elastic electron -lattice scattering, where only the direction of the excited electrons changes neglecting the energy loss, and inelastic electron -electron interaction, where secondary electrons are generated. An analysis of the particle kinetics inside the material shows that a smaller elastic scattering time affects the spin dynamics by leading to a larger spatial spread of electrons in the material, whereas generation of secondary electrons affects the spin transport with a larger time of propagation of homogeneous spin polarization.

This paper is concerned with a theoretical analysis of the behavior of optically excited spin currents in bilayer and multilayer systems of ferromagnetic and normal metals. As the propagation, control and manipulation of the spin currents created in ferromagnets by femtosecond optical pulses is of particular interest, we examine the influence of different thicknesses of the constituent layers for the case of electrons excited several electronvolts above the Fermi level. Using a Monte-Carlo simulation framework for such highly excited electrons, we first examine the spatio-temporal characteristics of the spin current density driven in a Fe layer, where the absorption profile of the light pulses plays an important role. Further, we examine how the combination of light absorption profiles, spin-dependent transmission probabilities, and iron layer thicknesses affect spin current density in a Fe/Au bilayer system. For high-energy electrons studied here, the interface and secondary electron generation have a small influence on spin transport in the bilayer system. However, we find that spin injection from one layer to another is most effective within a certain range of iron layer thicknesses.