We use the color diffusion (CD) algorithm in nonequilibrium (accelerated) ab initio molecular dynamics simulations to determine Ti monovacancy jump frequencies in NaCl-structure titanium nitride (TiN), at temperatures ranging from 2200 to 3000 K. Our results showthat theCDmethod extended beyond the linear-fitting rate-versus-force regime [Sangiovanni et al., Phys. Rev. B 93, 094305 (2016)] can efficiently determine metal vacancy migration rates in TiN, despite the low mobilities of lattice defects in this type of ceramic compound. We propose a computational method based on gamma-distribution statistics, which provides unambiguous definition of nonequilibrium and equilibrium (extrapolated) vacancy jump rates with corresponding statistical uncertainties. The acceleration-factor achieved in our implementation of nonequilibrium molecular dynamics increases dramatically for decreasing temperatures from 500 for T close to the melting point T-m, up to 33 000 for T approximate to 0.7 T-m

The first-principles calculation of many material properties, in particular related to defects and disorder, starts with the relaxation of the atomic positions of the system under investigation. This procedure is routine for nonmagnetic and magnetically ordered materials. However, when it comes to magnetically disordered systems, in particular the paramagnetic phase of magnetic materials, it is not clear how the relaxation procedure should be performed or which geometry should be used. Here we propose a method for the structural relaxation of magnetic materials in the paramagnetic regime, in an adiabatic fast-magnetism approximation within the disordered local moment (DLM) picture in the framework of density functional theory. The method is straightforward to implement using any ab initio code that allows for structural relaxations. We illustrate the importance of considering the disordered magnetic state during lattice relaxations by calculating formation energies and geometries for an Fe vacancy and C insterstitial atom in body-centered cubic (bcc) Fe as well as bcc Fe1-xCrx random alloys in the paramagnetic state. In the vacancy case, the nearest neighbors to the vacancy relax toward the vacancy of 0.14 angstrom (-5% of the ideal bcc nearest-neighbor distance), which is twice as large as the relaxation in the ferromagnetic case. The vacancy formation energy calculated in the DLM state on these positions is 1.60 eV, which corresponds to a reduction of about 0.1 eV compared to the formation energy calculated using DLM but on ferromagnetic-relaxed positions. The carbon interstitial formation energy is found to be 0.41 eV when the DLM relaxed positions are used, as compared to 0.59 eV when the FM-relaxed positions are employed. For bcc Fe0.5Cr0.5 alloys, the mixing enthalpy is reduced by 5 meV/atom, or about 10%, when the DLM state relaxation is considered, as compared to positions relaxed in the ferromagnetic state.

Density functional theory is a standard model for condensed-matter theory and computational material science. The accuracy of density functional theory is limited by the accuracy of the employed approximation to the exchange-correlation functional. Recently, the so-called strongly constrained appropriately normed (SCAN) [Sun, Ruzsinszky, and Perdew, Phys. Rev. Lett. 115, 036402 (2015)] functional has received a lot of attention due to promising results for covalent, metallic, ionic, as well as hydrogen- and van der Waals-bonded systems alike. In this work, we focus on assessing the performance of the SCAN functional for itinerant magnets by calculating basic structural and magnetic properties of the transition metals Fe, Co, and Ni. We find that although structural properties of bcc-Fe seem to be in good agreement with experiment, SCAN performs worse than standard local and semilocal functionals for fcc-Ni and hcp-Co. In all three cases, the magnetic moment is significantly overestimated by SCAN, and the 3d states are shifted to lower energies, as compared to experiments.

Recent theoretical investigations [A. B. Belonoshko et aL Nat. Geosci. 10, 312 (2017)] revealed the occurrence of the concerted migration of several atoms in bcc Fe at inner-core temperatures and pressures. Here, we combine first-principles and semiempirical atomistic simulations to show that a diffusion mechanism analogous to the one predicted for bcc iron at extreme conditions is also operative and of relevance for the high-temperature bcc phase of pure Ti at ambient pressure. The mechanism entails a rapid collective movement of numerous (from two to dozens) neighbors along tangled closed-loop paths in defect-free crystal regions. We argue that this phenomenon closely resembles the diffusion behavior of superionics and liquid metals. Furthermore, we suggest that concerted migration is the atomistic manifestation of vanishingly small co-mode phonon frequencies previously detected via neutron scattering and the mechanism underlying anomalously large and markedly non-Arrhenius self-diffusivities characteristic of bcc Ti.

Investigating magnetic materials at realistic conditions with first-principles methods is a challenging task due to the interplay of vibrational and magnetic degrees of freedom. The most difficult contribution to include in simulations is represented by the longitudinal magnetic degrees of freedom [longitudinal spin fluctuation (LSF)] due to their inherent many-body nature; nonetheless, schemes that enable to take into account this effect on a semiclassical level have been proposed and employed in the investigation of magnetic systems. However, assessment of the effect of vibrations on LSF is lacking in the literature. For this reason, in this work we develop a supercell approach within the framework of constrained density functional theory to calculate self-consistently the size of local-environment-dependent magnetic moments in the paramagnetic, high-temperature state in the presence of lattice vibrations and for liquid Fe in different conditions. First, we consider the case of bcc Fe at the Curie temperature and ambient pressure. Then, we perform a similar analysis on bcc Fe at Earths inner-core conditions, and we find that LSFs stabilize nonzero moments which affect atomic forces and the electronic density of states of the system. Finally, we employ the present scheme on liquid Fe at the melting point at ambient pressure and at Earths outer-core conditions (p approximate to 200 GPa, T approximate to 6000 K). In both cases, we obtain local magnetic moments of sizes comparable to the solid-state counterparts.

Iron-based alloys are widely used as structural components in engineering applications. This calls for a fundamental understanding of their mechanical properties, including those of pure iron. Under operational temperatures the mechanical and magnetic properties will differ from those of ferromagnetic body-centered-cubic iron at 0 K. In this theoretical work we study the effect of disordered magnetism on the screw dislocation core structure and compare with results for the ordered ferromagnetic case. Dislocation cores control some local properties such as the choice of glide plane and the associated dislocation mobility. Changes in the magnetic state can lead to modifications in the structure of the core and affect dislocation mobility. In particular, we focus on the core properties of the 1/2 < 111 > screw dislocation in the paramagnetic state. Using the noncollinear disordered local moment approximation to address paramagnetism, we perform structural relaxations within density functional theory. We obtain the dislocation core structure for the easy and hard cores in the paramagnetic state, and compare them with their ferromagnetic counterparts. By averaging the energy of several disordered magnetic configurations, we obtain an energy difference between the easy- and hard-core configurations, with a lower, but statistically close, value than the one reported for the ferromagnetic case. The magnetic moment and atomic volume at the dislocation core differ between paramagnetic and ferromagnetic states, with possible consequences on the temperature dependence of defect-dislocation interactions.

Chromium-nitride based materials have shown unexpected promise as thermo-electric materials for, e.g., wasteheat harvesting. Here, CrN and (Cr,V)N thin films were deposited by reactive magnetron sputtering. Thermoelectric measurements of pure CrN thin films show a low electrical resistivity between 1.2 and 1.5 x 10(-3) Omega cm and very high values of the Seebeck coefficient and thermoelectric power factor, in the range between 370-430 mu V/K and 9-11 x 10(-3) W/mK(2), respectively. Alloying of CrN films with small amounts (less than 15 %) of vanadium results in cubic (Cr,V)N thin films. Vanadium decreases the electrical resistivity and yields powerfactor values in the same range as pure CrN. Density functional theory calculations of sub-stoichiometric CrN1-delta and (Cr,V)N1-delta show that nitrogen vacancies and vanadium substitution both cause n-type conductivity and features in the band structure typically correlated with a high Seebeck coefficient. The results suggest that slight variations in nitrogen and vanadium content affect the power factor and offers a means of tailoring the power factor and thermoelectric figure of merit.