In absence of external fields, vibrational modes of periodic systems are usually considered as linearly polarized and, as such, they do not carry angular momentum. Our work proves that non-adiabatic effects due to the electron-phonon coupling are time-reversal symmetry breaking interactions for the vibrational field in systems with non-collinear magnetism and large spin-orbit coupling. Since in these systems the deformation potential matrix elements are necessarily complex, a nonzero synthetic gauge field (Berry curvature) arises in the dynamic equations of the ionic motion. As a result, phonon modes are elliptically polarized in the non-adiabatic framework and intrinsic vibrational angular momenta occur even for non-degenerate modes and without external probes. These results are validated by performing fully relativistic ab-initio calculations on two insulating platinum clusters and a metallic manganese compound, with non-collinear magnetism. In both cases, non-adiabatic vibrational modes carry sizeable angular momenta comparable to the orbital electronic ones in itinerant ferromagnets.

<p>The present thesis concerns method development and applications in the field of first principles electronic structure calculations.</p><p>Augmented planewaves combine the simple planewaves with exact solutions of the Schrödinger equation for a spherical potential. This combination yields a very good set of basis functions for describing the electronic structure everywhere in a crystal potential. In the present work, developments of the original augmented planewave (APW) method are presented. It is shown that the exact APW eigenvalues can be found using information from the eigenvalues of the APW secular matrix. This provides a more efficient scheme to solve the APW eigenvalue problem, than the traditional evaluation of the secular determinant. Further, a new way of linearizing the APW method is presented and compared to the traditional linearized APW method (LAPW). Using a combination of the original APW basis functions and the so called local orbitals (lo), the APW+lo linearization is found to reproduce the results of the LAPW method, but already at a smaller basis set size. Another advantage of the new linearization is a faster convergence of forces, with respect to the basis set size, as compared to the LAPW method.</p><p>The applications include studies of the non-collinear magnetic configuration in the fcc-based high-temperature phase of iron, γ-Fe. The system is found to be extremely sensitive to volume changes, as well as to a tetragonal distortion of the cubic unit cell. A continuum of degenerate spin spiral configurations, including the global energy minimum, are found for the undistorted crystal. The in-plane anisotropy of the ideal interface between a ferromagnetic layer of bcc Fe and the semiconducting ZnSe crystal is also investigated. In contrast to the four-fold symmetric arrangement of the atoms at the interface, the in-plane magnetic anisotropy displays a large uniaxiality. The calculated easy axes are in agreement with experiments for both Se and Zn terminated interfaces. In addition, calculations of the hyperfine parameters were performed for Li intercalated battery materials.</p>

In this thesis, non-collinear magnetism has been studied by using density functional theory and the augmented plane wave method with local orbitals (APW+lo). Two conditions for non-collinear instabilities have been identified in this thesis. First, the Fermi energy should cut through both spin up and down states. Secondly, strong nesting between the spin up and spin down Fermi surfaces is needed. The two criteria described here can be fulfilled by tuning the exchange-splitting and/or by modifying the volume. Calculations on several elements; bcc V, bcc and fcc Mn, bcc Fe, bcc and fcc Co, and bcc and fcc Ni show that a non-collinear state can be stabilized provided that the criteria discussed above are met. More complex materials have also been analyzed in terms of these two criteria. The substitutional alloys TlCo2Se2-xSx are found in experiments to possess spin spiral structures for x = {0-1.5} and at a concentration x = 1.75 the alloys become ferromagnetic. As S takes the place of Se in the crystal structure the distance between the Co layers is reduced and the turn angle of the spin spiral becomes smaller until it totally vanishes at x = 1.75. This thesis show that the evolution of the magnetic structure in these alloys is the consequence of a modification of the distance between Co layers, which induces a change in the interlayer exchange coupling. Fermi surfaces have been analyzed in TbNi5 in order to determine nesting features which would be responsible for the magnetic spin spiral observed in this material. The electronic structure of CeRhIn5 is also reported in this thesis. Furthermore, the 3-k magnetic structure of UO2 was investigated and the crystal field levels were calculated. Transition metal systems such as Fe in the superconducting high-pressure hcp phase and in the fcc crystal structure were also studied. The results obtained for fcc Fe are in accordance with previous reports. However the paramagnetic state in hcp Fe is found to be more stable than the antiferromagnetic configurations discussed earlier in the literature as being favored in the volume range where the hcp phase is stable and superconductivity appears (~ 15 GPa). The complex non-collinear magnetic structure in Mn3IrSi was calculated and the results are found to be in good agreement with experiments.

The present thesis concerns method development and applications in the field of first principles electronic structure calculations. Augmented planewaves combine the simple planewaves with exact solutions of the Schrödinger equation for a spherical potential. This combination yields a very good set of basis functions for describing the electronic structure everywhere in a crystal potential. In the present work, developments of the original augmented planewave (APW) method are presented. It is shown that the exact APW eigenvalues can be found using information from the eigenvalues of the APW secular matrix. This provides a more efficient scheme to solve the APW eigenvalue problem, than the traditional evaluation of the secular determinant. Further, a new way of linearizing the APW method is presented and compared to the traditional linearized APW method (LAPW). Using a combination of the original APW basis functions and the so called local orbitals (lo), the APW+lo linearization is found to reproduce the results of the LAPW method, but already at a smaller basis set size. Another advantage of the new linearization is a faster convergence of forces, with respect to the basis set size, as compared to the LAPW method. The applications include studies of the non-collinear magnetic configuration in the fcc-based high-temperature phase of iron, γ-Fe. The system is found to be extremely sensitive to volume changes, as well as to a tetragonal distortion of the cubic unit cell. A continuum of degenerate spin spiral configurations, including the global energy minimum, are found for the undistorted crystal. The in-plane anisotropy of the ideal interface between a ferromagnetic layer of bcc Fe and the semiconducting ZnSe crystal is also investigated. In contrast to the four-fold symmetric arrangement of the atoms at the interface, the in-plane magnetic anisotropy displays a large uniaxiality. The calculated easy axes are in agreement with experiments for both Se and Zn terminated interfaces. In addition, calculations of the hyperfine parameters were performed for Li intercalated battery materials.

The development of interatomic potentials for magnetic transition metals, and particularly for iron, is difficult, yet it is also necessary for large-scale atomistic simulations of industrially important iron and steel alloys. The magnetism of iron is especially important as it is responsible for many of the element's unique physical properties -- its bcc ground state structure, its high-temperature phase transitions, and the mobility of its self-interstitial atom (SIA) defects. Yet an accurate description of itinerant magnetism within a real-space formalism is particularly challenging and existing interatomic potentials based on the Embedded Atom Method are suited only for studies of near-equilibrium ferritic iron, due to their restricted functional forms. For this work, the magnetic analytic Bond-Order Potential (BOP) method has been implemented in full to test the convergence properties in both collinear and non-collinear magnetic iron. The known problems with negative densities of states (DOS) are addressed by assessing various possible definitions for the bandwidth and by including the damping factors adapted from the Kernel Polynomial Method. A 9-moment approximation is found to be sufficient to reproduce the major structural energy differences observed in Density Functional Theory (DFT) and Tight Binding (TB) reference calculations, as well as the volume dependence of the atomic magnetic moments. The Bain path connecting bcc and fcc structures and the formation energy of mono- and divacancies are also described well at this level of approximation. Other quantities such as the high-spin/low-spin transition in fcc iron, the bcc elastic constants and the SIA formation energies converge more slowly towards the TB reference data. The theory of non-collinear magnetism within analytic BOP is extended as required for a practical implementation. The spin-rotational behaviour of the energy is shown to converge more slowly than the collinear bulk energy differences, and there are specific problems at low angles of rotation where the magnitude of the magnetic moment depends sensitively on the detailed structure of the local DOS. Issues of charge transfer in relation to magnetic defects are discussed, as well as inadequacies in the underlying d-electron TB model.