Many-body theory of laser-induced ultrafast demagnetization and angular momentum transfer in ferromagnetic transition metals

An electronic theory is developed, which describes the ultrafast demagnetization in itinerant ferromagnets following the absorption of a femtosecond laser pulse. The present work intends to elucidate the microscopic physics of this ultrafast phenomenon by identifying its fundamental mechanisms. In particular, it aims to reveal the nature of the involved spin excitations and angular-momentum transfer between spin and lattice, which are still subjects of intensive debate. In the first preliminary part of the thesis the initial stage of the laser-induced demagnetization process is considered. In this stage the electronic system is highly excited by spin-conserving elementary excitations involved in the laser-pulse absorption, while the spin or magnon degrees of freedom remain very weakly excited. The role of electron-hole excitations on the stability of the magnetic order of one- and two-dimensional 3d transition metals (TMs) is investigated by using ab initio density-functional theory. The results show that the local magnetic moments are remarkably stable even at very high levels of local energy density and, therefore, indicate that these moments preserve their identity throughout the entire demagnetization process. In the second main part of the thesis a many-body theory is proposed, which takes into account these local magnetic moments and the local character of the involved spin excitations such as spin fluctuations from the very beginning. In this approach the relevant valence 3d and 4p electrons are described in terms of a multiband model Hamiltonian which includes Coulomb interactions, interatomic hybridizations, spin-orbit interactions, as well as the coupling to the time-dependent laser field on the same footing. An exact numerical time evolution is performed for small ferromagnetic TM clusters. The dynamical simulations show that after ultra-short laser pulse absorption the magnetization of these clusters decreases on a time scale of hundred femtoseconds. In particular, the results reproduce the experimentally observed laser-induced demagnetization in ferromagnets and demonstrate that this effect can be explained in terms of the following purely electronic non-adiabatic mechanism: First, on a time scale of 10–100 fs after laser excitation the spin-orbit coupling yields local angular-momentum transfer between the spins and the electron orbits, while subsequently the orbital angular momentum is very rapidly quenched in the lattice on the time scale of one femtosecond due to interatomic electron hoppings. In combination, these two processes result in a demagnetization within hundred or a few hundred femtoseconds after laser-pulse absorption.

Electroproduction of hyperons at low momentum transfer

A high resolution study of the H(e,e'K+)Λ,Σ 0 reaction was performed at Hall A, TJNAF as part of the hypernuclear experiment E94-107. One important ingredient to the measurement of the hypernuclear cross section is the elementary cross section for production of hyperons, Λ and Σ0. This reaction was studied using a hydrogen (i.e. a proton) target. Data were taken at very low Q2 (∼0.07 (GeV/c) 2) and W∼2.2 GeV. Kaons were detected along the direction of q, the momentum transferred by the incident electron (&thetas;CM∼6°). In addition, there are few data available regarding electroproduction of hyperons at low Q2 and &thetas;CM and the available theoretical models differ significantly in this kinematical region of W. The measurement of the elementary cross section was performed by scaling the Monte Carlo cross section (MCEEP) with the experimental-to-simulated yield ratio. The Monte Carlo cross section includes an experimental fit and extrapolation from the existing data for electroproduction of hyperons. Moreover, the estimated transverse component of the electroproduction cross section of H(e,e'K+)Λ was compared to the different predictions of the theoretical models and exisiting data curves for photoproductions of hyperons. None of the models fully describe the cross-section results over the entire angular range. Furthermore, measurements of the Σ 0/Λ production ratio were performed at &thetas; CM∼6°, where data are not available. Finally, data for the measurements of the differential cross sections and the Σ 0/Λ production were binned in Q2, W and &thetas;CM to understand the dependence on these variables. These results are not only a fundamental contribution to the hypernuclear spectroscopy studies but also an important experimental measurement to constrain existing theoretical models for the elementary reaction.

Electroproduction of hyperons at low momentum transfer

A high resolution study of the H(e,e'K+)Λ,Σ0 reaction was performed at Hall A, TJNAF as part of the hypernuclear experiment E94-107. One important ingredient to the measurement of the hypernuclear cross section is the elementary cross section for production of hyperons, Λ and Σ0. This reaction was studied using a hydrogen (i.e. a proton) target. Data were taken at very low Q2 (∼0.07 (GeV/c)2) and W∼2.2 GeV. Kaons were detected along the direction of q, the momentum transferred by the incident electron (θCM~6°). In addition, there are few data available regarding electroproduction of hyperons at low Q2 and θCM, and the available theoretical models differ significantly in this kinematical region of W. The measurement of the elementary cross section was performed by scaling the Monte Carlo cross section (MCEEP) with the experimental-to-simulated yield ratio. The Monte Carlo cross section includes an experimental fit and extrapolation from the existing data for electroproduction of hyperons. Moreover, the estimated transverse component of the electroproduction cross section of H(e,e'K+)Λ was compared to the different predictions of the theoretical models and exisiting data curves for photoproductions of hyperons. None of the models fully describe the cross-section results over the entire angular range. Furthermore, measurements of the Σ0/Λ production ratio were performed at θCM, where data are not available. Finally, data for the measurements of the differential cross sections and the Σ0/Λ production were binned in Q2, W and θCM to understand the dependence on these variables. These results are not only a fundamental contribution to the hypernuclear spectroscopy studies but also an important experimental measurement to constrain existing theoretical models for the elementary reaction.

Topological charge and angular momentum of light beams carrying optical vortices

We analyze the properties of light beams carrying phase singularities, or optical vortices. The transformations of topological charge during free-space propagation of a light wave, which is a combination of a Gaussian beam and a multiple charged optical vortex within a Gaussian envelope, are studied both in theory and experiment. We revise the existing knowledge about topological charge conservation, and demonstrate possible scenarios where additional vortices appear or annihilate during free propagation of such a combined beam. Coaxial interference of optical vortices is also analyzed, and the general rule for angular-momentum density distribution in a combined beam is established. We show that, in spite of any variation in the number of vortices in a combined beam, the total angular momentum is constant during the propagation. [S1050-2947(97)09910-1].

Entanglement in the steady state of a collective-angular-momentum (Dicke) model

We model the behavior of an ion trap with all ions driven simultaneously and coupled collectively to a heat bath. The equations for this system are similar to the irreversible dynamics of a collective angular momentum system known as the Dicke model. We show how the steady state of the ion trap as a dissipative many-body system driven far from equilibrium can exhibit quantum entanglement. We calculate the entanglement of this steady state for two ions in the trap and in the case of more than two ions we calculate the entanglement between two ions by tracing over all the other ions. The entanglement in the steady state is a maximum for the parameter values corresponding roughly to a bifurcation of a fixed point in the corresponding semiclassical dynamics. We conjecture that this is a general mechanism for entanglement creation in driven dissipative quantum systems.

Calculation of bound and resonance states of HO2 for nonzero total angular momentum

Bound and resonance states of HO2 have been calculated quantum mechanically by the Lanczos homogeneous filter diagonalization method [Zhang and Smith, Phys. Chem. Chem. Phys. 3, 2282 (2001); J. Chem. Phys. 115, 5751 (2001)] for nonzero total angular momentum J = 1,2,3. For lower bound states, agreement between the results in this paper and previous work is quite satisfactory; while for high lying bound states and resonances these are the first reported results. A helicity quantum number V assignment (within the helicity conserving approximation) is performed and the results indicate that for lower bound states it is possible to assign the V quantum numbers unambiguously, but for resonances it is impossible to assign the V helicity quantum numbers due to strong mixing. In fact, for the high-lying bound states, the mixing has already appeared. These results indicate that the helicity conserving approximation is not good for the resonance state calculations and exact quantum calculations are needed to accurately describe the reaction dynamics for HO2 system. Analysis of the resonance widths shows that most of the resonances are overlapping and the interferences between them lead to large fluctuations from one resonance to another. In accord with the conclusions from earlier J = 0 calculations, this indicates that the dissociation of HO2 is essentially irregular. (C) 2003 American Institute of Physics.

Off-axis vortices in trapped Bose-condensed gases: Angular momentum and frequency splitting

We consider noncentered vortices and their arrays in a cylindrically trapped Bose-Einstein condensate at zero temperature. We study the kinetic energy and the angular momentum per particle in the Thomas-Fermi regime and their dependence on the distance of the vortices from the center of the trap. Using a perturbative approach with respect to the velocity field of the vortices, we calculate, to first order, the frequency shift of the collective low-lying excitations due to the presence of an off-center vortex or a vortex array, and compare these results with predictions that would be obtained by the application of a simple sum-rule approach, previously found to be very successful for centered vortices. It turns out that the simple sum-rule approach fails for off-centered vortices.

Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores

Time correlation functions between the velocity of a tagged particle and velocities of particles within specified ranges of initial separations have been obtained by molecular dynamics simulation. These correlation functions have allowed us to analyze the momentum transfer between particles in different coordination shells. Two simple liquids at very different densities and two purely repulsive potentials with very different softnesses have been considered. The longitudinal correlations, which are the velocity cross-correlations along the initial direction defined by the centers of two given particles, have been calculated separately. It has been proven that these correlations should be attributed to particles both in front of and behind the central one. As with propagating longitudinal modes, they are strongly dependent on the softness of the potential core. Some characteristic features of the velocity correlation functions after the initial rise should be related to nonlongitudinal correlations. It has been shown that velocity cross-correlations between distinct particles cannot only be attributed to the direct interactions among particles, but also to the motions induced by the movement of a tagged particle on their neighbors.

Angular momentum conservation and gravity wave drag parameterization: implications for climate models

The robustness of the parameterized gravity wave response to an imposed radiative perturbation in the middle atmosphere is examined. When momentum is conserved and for reasonable gravity wave drag parameters, the response to a polar cooling induces polar downwelling above the region of the imposed cooling, with consequent adiabatic warming. This response is robust to changes in the gravity wave source spectrum, background flow, gravity wave breaking criterion, and model lid height. When momentum is not conserved, either in the formulation or in the implementation of the gravity wave drag parameterization, the response becomes sensitive to the above-mentioned factors—in particular to the model lid height. The spurious response resulting from nonconservation is found to be nonnegligible in terms of the total gravity wave drag–induced downwelling.

The angular momentum constraint on climate sensitivity and downward influence in the middle atmosphere

It is shown that under reasonable assumptions, conservation of angular momentum provides a strong constraint on gravity wave drag feedbacks to radiative perturbations in the middle atmosphere. In the time mean, radiatively induced temperature perturbations above a given altitude z cannot induce changes in zonal mean wind and temperature below z through feedbacks in gravity wave drag alone (assuming an unchanged gravity wave source spectrum). Thus, despite the many uncertainties in the parameterization of gravity wave drag, the role of gravity wave drag in middle-atmosphere climate perturbations may be much more limited than its role in climate itself. This constraint limits the possibilities for downward influence from the mesosphere. In order for a gravity wave drag parameterization to respect the momentum constraint and avoid spurious downward influence, any nonzero parameterized momentum flux at a model lid must be deposited within the model domain, and there must be no zonal mean sponge layer. Examples are provided of how violation of these conditions leads to spurious downward influence. For planetary waves, the momentum constraint does not prohibit downward influence, but it limits the mechanisms by which it can occur: in the time mean, downward influence from a radiative perturbation can only arise through changes in reflection and meridional propagation properties of planetary waves.

Simulations of angular momentum evolution in wind-fed CV precursors

It is thought that the secondary stars in cataclysmic variables (CVs) may undergo a period of mass loss in the form of a wind during the evolution of the system (Mullan et al. 1992). This wind is thought to magnetically brake the secondary star with a time-scale ~ 10^8 yr (e.g. van Paradijs 1986). When the secondary’s spin has been brought close to synchronism with the orbit it is possible for tidal torques to lock the secondary in synchronous rotation.

Midday auroral breakup events and related energy and momentum transfer from the magnetosheath

Combined observations by meridian-scanning photometers, all-sky auroral TV camera and the EISCAT radar permitted a detailed analysis of the temporal and spatial development of the midday auroral breakup phenomenon and the related ionospheric ion flow pattern within the 71°–75° invariant latitude radar field of view. The radar data revealed dominating northward and westward ion drifts, of magnitudes close to the corresponding velocities of the discrete, transient auroral forms, during the two different events reported here, characterized by IMF |BY/BZ| < 1 and > 2, respectively (IMF BZ between −8 and −3 nT and BY > 0). The spatial scales of the discrete optical events were ∼50 km in latitude by ∼500 km in longitude, and their lifetimes were less than 10 min. Electric potential enhancements with peak values in the 30–50 kV range are inferred along the discrete arc in the IMF |BY/BZ| < 1 case from the optical data and across the latitudinal extent of the radar field of view in the |BY/BZ| > 2 case. Joule heat dissipation rates in the maximum phase of the discrete structures of ∼ 100 ergs cm−2 s−1 (0.1 W m−2) are estimated from the photometer intensities and the ion drift data. These observations combined with the additional characteristics of the events, documented here and in several recent studies (i.e., their quasi-periodic nature, their motion pattern relative to the persistent cusp or cleft auroral arc, the strong relationship with the interplanetary magnetic field and the associated ion drift/E field events and ground magnetic signatures), are considered to be strong evidence in favour of a transient, intermittent reconnection process at the dayside magnetopause and associated energy and momentum transfer to the ionosphere in the polar cusp and cleft regions. The filamentary spatial structure and the spectral characteristics of the optical signature indicate associated localized ˜1-kV potential drops between the magnetopause and the ionosphere during the most intense auroral events. The duration of the events compares well with the predicted characteristic times of momentum transfer to the ionosphere associated with the flux transfer event-related current tubes. It is suggested that, after this 2–10 min interval, the sheath particles can no longer reach the ionosphere down the open flux tube, due to the subsequent super-Alfvénic flow along the magnetopause, conductivities are lower and much less momentum is extracted from the solar wind by the ionosphere. The recurrence time (3–15 min) and the local time distribution (∼0900–1500 MLT) of the dayside auroral breakup events, combined with the above information, indicate the important roles of transient magnetopause reconnection and the polar cusp and cleft regions in the transfer of momentum and energy between the solar wind and the magnetosphere.

Observations at the magnetopause and in the auroral ionosphere of momentum transfer from the solar wind

Recent radar studies of field-perpendicular flows in the auroral ionosphere, in conjunction with observations of the interplanetary medium immediately upstream of the Earth's bow shock, have revealed direct control of dayside convection by the Bz component of the interplanetary magnetic field (IMF). The ionospheric flows begin to respond to both northward and southward turnings of the IMF impinging upon the magnetopause after a delay of only a few minutes in the early afternoon sector, rising to about 15 minutes nearer dawn and dusk. In both the polar cap and the auroral oval, the subsequent rise and decay times are of order 5–10 minutes. We conclude there is very little convection “flywheel” effect in the dayside polar ionosphere and that only newly-opened flux tubes impart significant momentum to the ionosphere, in a relatively narrow region immediately poleward of the cusp. These findings concerning the effects of quasi-steady reconnection have important implications for any ionospheric signatures of transient reconnection which should be considerably shorter-lived than thought hitherto. In order to demonstrate the difficulty of uniquely identifying a Flux Transfer Event (FTE) in ground-based magnetometer data, we present observations of an impulsive signature, identical with that expected for an FTE if data from only one station is studied, following an observed magnetopause compression when the IMF was purely northward. We also report new radar observations of a viscous-like interaction, consistent with an origin on the flanks of the magnetotail and contributing an estimated 15–30kV to the total cross-cap potential during quiet periods.