Disorder in crystalline materials can take different forms and originate from different sources. In particular, temperature introduces disorder in any kind of material. This can be observed as the appearance of vacant lattice sites in an otherwise perfect crystal, or as a random distribution of different elements on the same lattice in an alloy; at the same time, if the material is magnetic, temperature induces disorder also on the magnetic degrees of freedom. In this thesis, different levels of disorder associated to structure and magnetism are investigated by means of density functional theory and thermodynamic models.

I start with diffusion of Ti vacancies in TiN, which is studied by means of nonequilibrium ab initio molecular dynamics using the color diffusion algorithm at different temperatures. The result is an Arrhenius behavior of Ti vacancy jump rates.

A method to perform structural relaxations in magnetic materials in their hightemperature paramagnetic phase is then developed based on the disordered local moments approach in order to study vacancies, interstitial atoms, and combinations of defects in paramagnetic bcc Fe and B1 CrN, as well as the mixing enthalpy of bcc Fe1−xCrx random alloys. A correction to the energetics of every system due to the relaxation in the disordered magnetic state is observed in all cases.

Not related to temperature and disorder, but very important for an accurate description of magnetic materials, is the choice of the exchange and correlation functional to be employed in the first principles calculations. We have investigated the performance of a recently developed meta-GGA functional, the strongly constrained and appropriately normed (SCAN) functional, in comparison with the more commonly used LDA and PBE on the ferromagnetic elemental solids bcc Fe, fcc Ni, and hcp Co, and SCAN it is found to give negligible improvements, if not a worsening, in the description of these materials.

Finally, the coupling between vibrational and magnetic degrees of freedom is discussed by reviewing the literature and proposing an investigation of the influence of vibrations on longitudinal spin fluctuations. These excitations are here studied by means of thermodynamic models based on Landau expansion of the energy in even powers of the magnitude of the local magnetic moments. We find that vibrational and magnetic disorder alter the energy landscapes as a function of moment size also in bcc Fe, which is often considered a Heisenberg system, inducing a more itinerant electron behavior.

The modeling of magnetic materials at finite temperatures is an ongoing challenge in the field of theoretical physics. This field has strongly benefited from the development of computational methods, which allow to predict material’s properties and explain physical effects on the atomic scale, and are now employed to direct the design of new materials. However, simulations need to be as accurate as possible to give reliable insights into solid-state phenomena, which means that, most desirably, all competing effects occurring in a system at realistic conditions should be included. This task is particularly difficult in the modeling of magnetic materials from first principles, due to the quantum nature of magnetism and its interplay with other phenomena related to the atomic degrees of freedom. The aim of this thesis is therefore to develop methods that enable the inclusion of magnetic effects in finite temperature simulations based on density functional theory (DFT), while considering on the same footing vibrational and structural degrees of freedom,with a particular focus on the high-temperature paramagnetic phase. The type of couplings investigated in this thesis can be separated in two big categories: interplay between magnetism and structure, and between magnetism and vibrations.

Regarding the former category, I have tried to shine some light on the effect of the paramagnetic state on atomic positions in a crystal in the presence of defects or for complicated systems, as opposed to the ordered magnetic state. To model the high-temperature paramagnetic phase of magnetic materials, the disordered local moment (DLM) approach is employed in the whole work. In this framework, I have developed a method to perform local lattice relaxations in the disordered magnetic state, which consists of a step-wise partial relaxation of the atomic positions, while changing the configuration of the magnetic moments at each step of the procedure. This method has been tested on point defects in paramagnetic bcc Fe, namely the single vacancy and, separately, the C interstitial in octahedral position, and on Fe_1-xCr_x alloys, finding non-negligible effects on formation energies. In addition, the feasibility of investigating extended defects like dislocations in the paramagnetic state with this method has also been proven by studying the screw dislocation in bcc Fe. The DLM-relaxation method has then been used to investigate intrinsic and extrinsic defects in CrN, an antiferromagnetic semiconductor studied for thermoelectric applications, found in the paramagnetic state at operating temperature, and a newly synthesized compound, Fe₃CO₇, which features a complicated crystal structure and unusual electronic properties, with possible important implications for the chemistry of Earth’s mantle.

The other focus of this thesis is the coupling between magnetism and lattice vibrations. As a pre-step to perform fully coupled atomistic spin dynamics-ab initio molecular dynamics (ASD-AIMD) simulations, I have first investigated the effect of vibrations on the so called longitudinal spin fluctuations, a mechanism occurring at finite temperatures and important for itinerant electron magnetic systems. I have developed a framework to investigate the dependence of the local moment’s energy landscapes on the instantaneous positions of the atoms, testing it on Fe at different temperature and pressure conditions. This study has laid the foundation to apply machine learning techniques to the prediction of the energy landscapes during an ASD-AIMD simulation. Finally, I have investigated the phase stability of Fe at ambient pressure from the theoretical Curie temperature up to its melting point with ASD-AIMD. This task is carried out by applying a pool of thermodynamic techniques to calculate free energy differences, and therefore I have defined a strategy to discern the thermodynamic equilibrium structure in magnetic materials in the high temperature paramagnetic phase based on first principles dynamical simulations. The methodologies developed and applied in this work constitute an improvement towards the simulation of magnetic materials accounting for the coupling of all effects, and the hope is to bridge a gap between theory and experiments.