Explicit anharmonicity, defined as the vibrational contribution beyond the quasiharmonic approximation, is qualitatively different between the group V and group VI bcc refractory elements. Group V elements show a small and mostly negative anharmonic entropy, whereas group VI elements have a large positive anharmonic entropy, strongly increasing with temperature. Here, we explain this difference utilizing highly accurate anharmonic free energies and entropies from ab initio calculations for Nb and Ta (group V), and Mo and W (group VI). The numerically calculated entropies are in agreement with prior experimental data. The difference in behavior between the two sets of elements arises not from their high-temperature behavior but rather from the 0K quasiharmonic reference state. We understand this by analyzing the 0K and the high-temperature phonon density of states and the electronic density of states. The qualitative difference disappears when the anharmonicity is instead referenced with a high-temperature effective harmonic potential. However, even for an optimized effective harmonic reference, the remaining effective anharmonicity is significant. The reason is that the anharmonicity in the bcc systems - carried by asymmetric distributions in the nearest neighbors - can never be accounted for by a harmonically restricted potential.

The microstructure of advanced high-strength steels often shows a sensitive dependence on alloying. For example, adding Cr to improve the corrosion resistance of medium-Mn steels also enhances the precipitation of carbides. The current study focuses on the behavior of H in such complex multicomponent carbides by employing different methodological strategies. We systematically analyze the impact of Cr, Mn, and Fe using density functional theory (DFT) for two prototype precipitate phases, M3C and M23C6, where M represents the metal sublattice. Our results show that the addition of these alloying elements yields strong nonmonotonic chemical trends for the H solubility. We identify magnetovolume effects as the origin for this behavior, which depend on the considered system, the sites occupied by H, and short- vs long-range interactions between H and the alloying elements. We further show that the H solubility is directly correlated with the occupation of its nearest-neighbor shells by Cr and Mn. Based on these insights, DFT data from H containing binary-metal carbides are used to design a ridge regression based model that predicts the solubility of H in the ternary-metal carbides (Fe-Cr-Mn-C).

Since its first emergence in 2004, the high-entropy alloy (HEA) concept has aimed at stabilizing single- or dual-phase multi-element solid solutions through high mixing entropy. Here, this strategy is changed and renders such massive solid solutions metastable, to trigger spinodal decomposition for improving the alloys’ magnetic properties. The motivation for starting from a HEA for this approach is to provide the chemical degrees of freedom required to tailor spinodal behavior using multiple components. The key idea is to form Fe-Co enriched regions which have an expanded volume (relative to unconstrained Fe-Co), due to coherency constraints imposed by the surrounding HEA matrix. As demonstrated by theory and experiments, this leads to improved magnetic properties of the decomposed alloy relative to the original solid solution matrix. In a prototype magnetic FeCoNiMnCu HEA, it is shown that the modulated structures, achieved by spinodal decomposition, lead to an increase of the Curie temperature by 48% and a simultaneous increase of magnetization by 70% at ambient temperature as compared to the homogenized single-phase reference alloy. The findings thus open a pathway for the development of advanced functional HEAs.

The interface method is a well established approach for predicting melting points of materials using interatomic potentials. However, applying the interface method is tedious and involves significant human intervention. The whole procedure involves several successive tasks: estimate a rough melting point, set up the interface structure, run molecular dynamic calculations and analyze the data. Loop calculations are necessary if the predicted melting point is different from the estimated one by more than a certain convergence criterion, or if full melting/solidification occurs. In this case monitoring the solid–liquid phase transition in the interface structure becomes critical. As different initial random seeds for the molecular dynamic simulations within the interface method induce slightly different melting points, a few ten or hundred interface method calculations with different random seeds are necessary for performing a statistical analysis on these melting points. Considering all these technical details, the work load for manually executing and combining the various involved scripts and programs quickly becomes prohibitive. To simplify and automatize the whole procedure, we have implemented the interface method into pyiron (http://pyiron.org). Our fully automatized procedure allows to efficiently and precisely predict melting points of stable unaries represented by arbitrary potentials with only two user-specified parameters (interatomic potential file and element). For metastable or dynamically unstable unary phases, the crystal structure needs to be provided as an additional parameter. We have applied our automatized approach on fcc Al, Ni, dynamically unstable bcc Ti and hcp Mg and employed a large set of available interatomic potentials. Melting points for classical interatomic potentials of these metals have been obtained with a numerical precision well below 1 K.

The lattice dynamics in magnetic materials, such as Fe depends on the degree of disorder of the atomic magnetic moments and the time scale of spin fluctuations. Using first-principles methods, we have studied this effect by determining the force constant matrix in two limits: (i) When spin fluctuations are much faster than the atom vibrations, their combined impact is captured by a spin-space averaged force constant matrix, (ii) when individual spin fluctuations are sufficiently slow to scatter the phonon modes, the itinerant coherent potential approximation with spin-pair resolved force constants (i.e., φ↑↑,φ↓↓, and φ↑↓) is employed in this paper. The physical consequences for the vibrational spectral functions are analyzed by systematically modifying the input parameters (magnetization and ratio of force constants betweens atoms with equal and opposite spin directions) and by deriving them for the prototype material system bcc Fe from first-principles calculations. In the paramagnetic regime, the two limits yield identical phonon spectra. Below the Curie temperature, however, there are regions in the parametric study that show qualitative differences, including a broadening of the phonon peaks. For bcc Fe, however, the quantitative modifications of phonon frequencies turn out to be small.

We reveal the impact of magnetic ordering on stacking fault energy (SFE) and its influence on the deformation mechanisms and mechanical properties in a class of nonequiatomic quinary Mn-containing compositional complex alloys or high entropy alloys (HEAs). By combining ab initio simulation and experimental validation, we demonstrate magnetic ordering as an important factor in the activation and transition of deformation modes from planar dislocation slip to TWIP (twinning-induced plasticity) and/or TRIP (transformation-induced plasticity). A wide compositional space of Cr20MnxFeyCo20Niz(x+y+z=60, at. %) was probed by density-functional theory calculations to search for potential alloys displaying the TWIP/TRIP effects. Three selected promising HEA compositions with varying Mn concentrations were metallurgically synthesized, processed, and probed for microstructure, deformation mechanism, and mechanical property evaluation. The differences in the deformation modes of the probed HEAs are interpreted in terms of the computed SFEs and their dependence on the predicted magnetic state, as revealed by ab initio calculations and validated by explicit magnetic measurements. It is found that the Mn content plays a key role in the stabilization of antiferromagnetic configurations which strongly impact the SFEs and eventually lead to the prevalent deformation behavior.

Concentrated solid solutions including the class of high entropy alloys (HEAs) have attracted enormous attention recently. Among these alloys a recently developed face-centered cubic (fcc) equiatomic VCoNi alloy revealed extraordinary high yield strength, exceeding previous high-strength fcc CrCoNi and FeCoNiCrMn alloys. Significant lattice distortions had been reported in the VCoNi solid solution. There is, however, a lack of knowledge about potential short-range order (SRO) and its implications for most of these alloys. We performed first-principles calculations and Monte Carlo simulations to compute the degree of SRO for fcc VCoNi, namely, by utilizing the coherent-potential approximation in combination with the generalized perturbation method as well as the supercell method in combination with recently developed machine-learned potentials. We analyze the chemical SRO parameters as well as the impact on other properties such as relaxation energies and lattice distortions.

To facilitate the understanding of Invar effects and design of FeNiCo-base Invar alloys characterized by low thermal expansion coefficient (TEC), we investigated the magnetic and thermal expansion behavior of an equiatomic prototype medium entropy alloy FeNiCo and a non-equiatomic (super Invar)Fe₆₃Ni₃₂Co₅ (at. %)reference alloy by means of experiments and ab initio calculations. Both alloys consist of a single face-centered cubic phase with fully recrystallized microstructure in the homogenized state. Large spontaneous volume magnetostriction is observed in both alloys below their respective Curie temperatures. The Invar effect in the non-equiatomic Fe₆₃Ni₃₂Co₅ alloy is of step-type with nearly zero TEC over a wide temperature range (from room temperature to 120 °C)below its Curie temperature. The equiatomic FeNiCo alloy shows a peak-type Invar effect in a very narrow temperature range (from ∼675 °C to ∼730 °C)with relatively low TECs. The equiatomic FeNiCo alloy shows both higher saturation magnetization and Curie temperature than the non-equiatomic Fe₆₃Ni₃₂Co₅ alloy. The relationships among magnetic behavior, spontaneous volume magnetostriction and Invar effects for a wider array of metallic alloys are discussed mainly based on Masumoto's rule combined with Wohlfarth's itinerant electron theory. An Invar alloy search map is constructed based on the present results and available literature data to visualize the relationships among saturation magnetization, Curie temperature and thermal expansion coefficient for a wide range of Invar alloys. Based on this treasure map a design route for further developments of new Invar alloys by tuning their magnetic properties is discussed.

Severe lattice distortion is a core effect in the design of multiprincipal element alloys with the aim to enhance yield strength, a key indicator in structural engineering. Yet, the yield strength values of medium- and high-entropy alloys investigated so far do not substantially exceed those of conventional alloys owing to the insufficient utilization of lattice distortion. Here it is shown that a simple VCoNi equiatomic medium-entropy alloy exhibits a near 1 GPa yield strength and good ductility, outperforming conventional solid-solution alloys. It is demonstrated that a wide fluctuation of the atomic bond distances in such alloys, i.e., severe lattice distortion, improves both yield stress and its sensitivity to grain size. In addition, the dislocation-mediated plasticity effectively enhances the strength–ductility relationship by generating nanosized dislocation substructures due to massive pinning. The results demonstrate that severe lattice distortion is a key property for identifying extra-strong materials for structural engineering applications.

Recently, high-entropy alloys (HEAs) have attracted wide attention due to their extraordinary materials properties. A main challenge in identifying new HEAs is the lack of efficient approaches for exploring their huge compositional space. Ab initio calculations have emerged as a powerful approach that complements experiment. However, for multicomponent alloys existing approaches suffer from the chemical complexity involved. In this work we propose a method for studying HEAs computationally. Our approach is based on the application of machine-learning potentials based on ab initio data in combination with Monte Carlo simulations. The high efficiency and performance of the approach are demonstrated on the prototype bcc NbMoTaW HEA. The approach is employed to study phase stability, phase transitions, and chemical short-range order. The importance of including local relaxation effects is revealed: they significantly stabilize single-phase formation of bcc NbMoTaW down to room temperature. Finally, a so-far unknown mechanism that drives chemical order due to atomic relaxation at ambient temperatures is discovered.

We study the mutual coupling of spin fluctuations and lattice vibrations in paramagnetic CrN by combining atomistic spin dynamics and ab initio molecular dynamics. The two degrees of freedom are dynamically coupled, leading to nonadiabatic effects. Those effects suppress the phonon lifetimes at low temperature compared to an adiabatic approach. The dynamic coupling identified here provides an explanation for the experimentally observed unexpected temperature dependence of the thermal conductivity of magnetic semiconductors above the magnetic ordering temperature.

Medium and high entropy alloys (MEAs and HEAs) based on 3d transition metals, such as face-centered cubic (fcc) CrCoNi and CrMnFeCoNi alloys, reveal remarkable mechanical properties. The stacking fault energy (SFE) is one of the key ingredients that controls the underlying deformation mechanism and hence the mechanical performance of materials. Previous experiments and simulations have therefore been devoted to determining the SFEs of various MEAs and HEAs. The impact of local chemical environment in the vicinity of the stacking faults is, however, still not fully understood. In this work, we investigate the impact of the compositional fluctuations in the vicinity of stacking faults for two prototype fcc MEAs and HEAs, namely CrCoNi and CrMnFeCoNi by employing first-principles calculations. Depending on the chemical composition close to the stacking fault, the intrinsic SFEs vary in the range of more than 150 mJ/m² for both the alloys, which indicates the presence of a strong driving force to promote particular types of chemical segregations towards the intrinsic stacking faults in MEAs and HEAs. Furthermore, the dependence of the intrinsic SFEs on local chemical fluctuations reveals a highly non-linear behavior, resulting in a non-trivial interplay of local chemical fluctuations and SFEs. This sheds new light on the importance of controlling chemical fluctuations via tuning, e.g., the annealing condition to obtain the desired mechanical properties for MEAs and HEAs.

A first-principles-based computational tool for simulating phonons of magnetic random solid solutions including thermal magnetic fluctuations is developed. The method takes fluctuations of force constants due to magnetic excitations as well as due to chemical disorder into account. The developed approach correctly predicts the experimentally observed unusual phonon hardening of a transverse acoustic mode in Fe-Pd an Fe-Pt Invar alloys with increasing temperature. This peculiar behavior, which cannot be explained within a conventional harmonic picture, turns out to be a consequence of thermal magnetic fluctuations. The proposed methodology can be straightforwardly applied to a wide range of materials to reveal new insights into physical behaviors and to design materials through computation, which were not accessible so far.

Quantum-mechanical calculations are used to determine the temperature dependence of the Gibbs energy of vacancy formation in nickel. Existing data reveal a discrepancy between the high-temperature estimates from experiments and low-temperature approximations from density functional theory. Our finite-temperature calculations - which include the effects of magnetism and fully interacting phonon vibrations - demonstrate that this discrepancy is mostly caused by the previously neglected explicit anharmonic contribution.

The temperature-dependent intrinsic stacking fault Gibbs energy is computed based on highly converged density-functional-theory (DFT) calculations for the three prototype face-centered cubic metals Al, Cu, and Ni. All relevant temperature-dependent contributions are considered including electronic, vibrational, magnetic, and explicit anharmonic Gibbs energy contributions as well as coupling terms employing state-of-the-art statistical sampling techniques. Particular emphasis is put on a careful comparison of different theoretical concepts to derive the stacking fault energy such as the axial-next-nearest-neighbor-Ising (ANNNI) model or the vacuum-slab approach. Our theoretical results are compared with an extensive set of previous theoretical and experimental data. Large uncertainties in the experimental data highlight the necessity of complementary parameter-free calculations. Specifically, the temperature dependence is experimentally unknown and poorly described by thermodynamic databases. Whereas calphad derived data shows an increase of the stacking fault energy with temperature for two of the systems (Cu and Ni), our results predict a decrease for all studied systems. For Ni, the temperature induced change is in fact so strong that in the temperature interval relevant for super-alloy applications the stacking fault energy falls below one third of the low temperature value. Such large changes clearly call for a revision of the stacking fault energy when modeling or designing alloys based on such elements.

We introduce a new class of high-entropy alloys (HEAs), i.e., quinary (five-component) dual-phase (DP) HEAs revealing transformation-induced plasticity (TRIP), designed by using a quantum mechanically based and experimentally validated approach. Ab initio simulations of thermodynamic phase stabilities of Co₂₀Cr₂₀Fe_40-xMn₂₀Ni_x (x = 0–20 at. %) HEAs were performed to screen for promising compositions showing the TRIP-DP effect. The theoretical predictions reveal several promising alloys, which have been cast and systematically characterized with respect to their room temperature phase constituents, microstructures, element distributions and compositional homogeneity, tensile properties and deformation mechanisms. The study demonstrates the strength of ab initio calculations to predict the behavior of multi-component HEAs on the macroscopic scale from the atomistic level. As a prototype example a non-equiatomic Co₂₀Cr₂₀Fe₃₄Mn₂₀Ni₆ HEA, selected based on our ab initio simulations, reveals the TRIP-DP effect and hence exhibits higher tensile strength and strain-hardening ability compared to the corresponding equiatomic CoCrFeMnNi alloy.

The spin-space averaging formalism is applied to compute atomic forces and phonon spectra for magnetically excited states of fcc nickel. Transverse and longitudinal magnetic fluctuations are taken into account by a combination of magnetic special quasi random structures and constrained spin-density-functional theory. It turns out that for fcc Ni interatomic force constants and phonon spectra are almost unaffected by both kinds of spin fluctuations. Given the computational expense to simulate coupled magnetic and atomic fluctuations, this insight facilitates computational modeling of magnetic alloys such as Ni-based superalloys.