Search

The B-Pt alloy system contains several Pt-rich phases exhibiting complex structures, many with partial site occupation. It also exhibits a deep (nearly 1000 degrees C) eutectic. We evaluate the ab-initio total energies of the crystalline solids to clarify the identity and character of the phases. Our work identifies inconsistencies in materials databases such as the Inorganic Crystallographic Structure Database and the ASM Phase Diagram Database, but our total energy calculations allow us to match up experimentally reported stable phases with specific structures. High temperature Gibbs free energy calculations in the liquid and solid states at the composition Pt$_2$B reveal that the depth of the eutectic arises from the low energy and high entropy of the liquid state., Comment: Proceedings of ISBB 2024

The AlCrTiV high entropy alloy undergoes an order-disorder transition from body centered cubic (Strukturbericht A2) at high temperatures to the CsCl structure (B2) at intermediate temperatures. We model this transition using first principles Monte Carlo/molecular dynamics simulations. Simulation results yield the temperature-dependent energy, entropy, heat capacity, occupancy fluctuations, and diffraction patterns. The contribution of chemical disorder to the entropy is calculated on the basis of point and pair cluster frequencies. The simulated structures exhibit compensated ferrimagnetism, and the Fermi level lies in a pseudogap. Sensitivity of structure and magnetism to the exchange-correlation functional is discussed, and neutron diffraction experiments are proposed to help resolve the true chemical order.

Vibrational properties of solids are key to determining stability, response and functionality. However, they are challenging to computationally predict at Ab-Initio accuracy, even for elemental systems. Ab-Initio methods for modeling atomic interactions are limited in the system sizes and simulation times that can be achieved. Due to these limitations, Machine Learning Interatomic Potentials (MLIPs) are gaining popularity and success as a faster, more scalable approach for modeling atomic interactions, potentially at Ab-Initio accuracy. Even with faster potentials, methodologies for predicting entropy, free energy and vibrational properties vary in accuracy, cost and difficulty to implement. Using the Covariance of Atomic Displacements (CAD) to predict entropy, free energy and finite-temperature phonon dispersions is a promising approach but thorough benchmarking has been hampered by the cost of Ab-Initio methods for sampling. In this work, we use a MLIP and the CAD to characterize the convergence of the predicted properties and determine optimal sampling strategies. We focus on solid lithium at zero pressure, showing that the MLIP-CAD approach reproduces experimental entropy, phonon dispersions and the martensitic transition while also comparing to more established methods.

Atomic structures of Al-Co-Cu decagonal quasicrystals (QCs) are investigated using empirical oscillating pair potentials (EOPP) in molecular dynamic (MD) simulations that we enhance by Monte Carlo (MC) swapping of chemical species and replica exchange. Predicted structures exhibit planar decagonal tilin g patterns and are periodic along the perpendicular direction. We then recalculate the energies of promising structures using first-principles density functional theory (DFT), along with energies of competing phases. We find that our $\tau$-inflated sequence of QC approximants are energetically unstable a t low temperature by at least 3 meV/atom. Extending our study to finite temperatures by calculating harmonic vibrational entropy, as well as anharmonic contributions that include chemical species swaps and tile flips, our results suggest that the quasicrystal phase is entropically stabilized at temperatur es in the range 600-800K and above. It decomposes into ordinary (though complex) crystal phases at low temperatures, including a partially disordered B2-type phase. We discuss the influence of density and composition on QC phase stability; we compare the structural differences between Co-rich and Cu-rich quasicrystals; and we analyze the role of entropy in stabilizing the quasicrystal, concluding with a discussion of the possible existence of "high entropy" quasicrystals.

The In-Sn binary alloy system exhibits several unusual features that challenge crystallographic and thermodynamic expectations. We combine first principles total energy calculation with simple thermodynamic modeling to address two key points. First, we evaluate energies along the Bain path to interpret the discontinuous transition between the phases $\alpha$-In (Pearson type tI2) and $\beta$--In$_3$Sn (also Pearson type tI2) that are identical in symmetry. Second, we demonstrate that the solid solution phases $\beta$-In$_3$Sn and $\gamma$-InSn$_4$ (Pearson type hP1) exist at high temperatures only, and they exhibit eutectoid decompositions at low temperatures.

Liquid state entropy formulas based on configurational probability distributions are examined for Lennard-Jones fluids across a range temperatures and densities. These formulas are based on expansions of the entropy in series of $n$-body distribution functions. We focus on two special cases. One, which we term the ``perfect gas'' series, starts with the entropy of an ideal gas; the other, which we term the ``dense liquid series'' removes a many-body contribution from the ideal gas entropy and reallocates it among the subsequent $n$-body terms. We show that the perfect gas series is most accurate at low density, while the dense liquid series is most accurate at high density. We propose empirical interpolation methods that are capable of connecting the two series and giving consistent predictions in most situations.

Refractory high-entropy alloys, RHEAs, are promising high-temperature structural materials. Their large compositional space poses great design challenges for phase control and high strength-ductility synergy. The present research pioneers using integrated high-throughput machine learning with Monte Carlo simulations to effectively navigate phase-selection and mechanical-properties predictions, developing aluminum-containing RHEAs in single-phase ordered B2 alloys demonstrating both high strength and ductility. These aluminum-containing RHEAs achieve remarkable mechanical properties, including compressive yield strengths up to 1.6 GPa, fracture strains exceeding 50 percent, and significant high-temperature strength retention. They also demonstrate a tensile yield strength of 1.1 GPa with a tension ductility of 6.3 percent. Besides, we identify a valence-electron-count domain for alloy brittleness with the explanation from density-functional theory and provide crucial insights into elements' influence on atomic ordering and mechanical performance. The work sets forth a strategic blueprint for high-throughput alloy design and reveals fundamental principles that govern the mechanical properties of advanced structural alloys.

Refractory metals exhibit high strength at high temperature, but often lack ductility. Multiprinciple element alloys such as high entropy alloys offer the potential to improve ductility while maintaining strength, but we don't know $a-priori$ what compositions will be suitable. A number of measures have been proposed to predict the ductility of metals, notably the Pugh ratio, the Rice-Thomson D-parameter, among others. Here we examine direct $ab-initio$ simulation of deformation under tensile strain, and we apply this to a variety of Nb- and Mo-based binary alloys and to several quaternary alloy systems. Our results exhibit peak stresses for elastic deformation, beyond which defects such as lattice slip, stacking faults, transformation, and twinning, relieve the stress. The peak stress grows strongly with increasing valence electron count. Correlations are examined among several physical properties, including the above-mentioned ductility parameters., Comment: 22 pages, 13 figures, accepted for publication in Physical Review Materials

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with large number of atoms. In this work, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate first-principles residual resistivity of large systems. This has been implemented in the open-source code lsms. We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (KKR-CPA). We discuss future applications of this method to complex systems where other methods are not applicable., Comment: 16 pages, 5 figures

Predicting quasicrystal structures is a multifaceted problem that can involve predicting a previously unknown phase, predicting the structure of an experimentally observed phase, or predicting the thermodynamic stability of a given structure. We survey the history and current state of these prediction efforts with a focus on methods that have improved our understanding of the structure and stability of known metallic quasicrystal phases. Advances in the structural modeling of quasicrystals, along with first principles total energy calculation and statistical mechanical methods that enable the calculation of quasicrystal thermodynamic stability, are illustrated by means of cited examples of recent work.

This paper reports formation enthalpies of phases in the Al-Mn-Pd ternary alloy system as calculated from first principles using electronic density functional theory. We consider all crystal structures as reported in the assessed phase diagrams of the ternary and its binary alloy subsystems (Al-Mn, Al-Pd, and Mn-Pd), as well as additional reported or hypothetical structures. Icosahedral and decagonal quasicrystalline approximants are among the structures that we predict to be stable, or nearly so. Our results suggest the need for careful experimental reexamination of phase stability in each of the alloy systems, in tandem with further efforts to refine crystallographic and ab-initio structures.

We introduce a simplified method to model the interatomic interactions of high entropy alloys based on a lookup table of cluster energies. These interactions are employed in replica exchange Monte Carlo simulations with histogram analysis to obtain thermodynamic properties across a broad temperature range. Kikuchi's Cluster Variation Method entropy formalism and high temperature series expansions are applied to directly calculate entropy from statistics on short- and long-range chemical order, and we discuss the convergence of the entropy as clusters of differing size are included.

Ab-initio crystal structure prediction depends on accurate calculation of the energies of competing structures. Many DFT codes are available that utilize different approaches to solve the Kohn-Sham equation. We evaluate the consistency of three software packages (WIEN2k, VASP and MuST) that utilize three different methods (FL-APW, plane-wave pseudopotential and the KKR-Green's Function methods) by comparing the relative total energies obtained for a set of BCC and FCC binary metallic alloys. We focus on the impact of choices such as muffin-tin $vs.$ full-potential, angular momentum cutoff and other important KKR parameters. Different alloy systems prove more or less sensitive to these choices, and we explain the differences through study of the angular variation of their potentials. Our results can provide guidance in the application of KKR as a total energy method for structure prediction., Comment: 18 pages, 9 figures

The vibrational entropy of a solid at finite temperature is investigated from the perspective of information theory. Ab initio molecular dynamics (AIMD) simulations generate ensembles of atomic configurations at finite temperature from which we obtain the $N$-body distribution of atomic displacements, $\rho_N$. We calculate the information-theoretic entropy from the expectation value of $\ln{\rho_N}$. At a first level of approximation, treating individual atomic displacements independently, our method may be applied using Debye-Waller B-factors, allowing diffraction experiments to obtain an upper bound on the thermodynamic entropy. At the next level of approximation we correct the overestimation through inclusion of displacement covariances. We apply this approach to elemental body-centered cubic sodium and face-centered cubic aluminum, showing good agreement with experimental values above the Debye temperatures of the metals. Below the Debye temperatures we extract an effective vibrational density of states from eigenvalues of the covariance matrix, and then evaluate the entropy quantum mechanically, again yielding good agreement with experiment down to low temperatures. Our method readily generalizes to complex solids, as we demonstrate for a high entropy alloy. Further, our method applies in cases where the quasiharmonic approximation fails, as we demonstrate by calculating the HCP/BCC transition in Ti.

In this paper we obtain formulas for the distribution of the left-most up-spin in the Heisenberg-Ising spin-1/2 chain with anisotropy parameter $\Delta$, also known as the XXZ spin-1/2 chain, on the one-dimensional lattice $\mathbb{Z}$ with domain wall initial conditions. We use the Bethe Ansatz to solve the Schr$\"o$dinger equation and a recent antisymmetrization identity of Cantini, Colomo, and Pronko (arXiv:1906.07636) to simplify the marginal distribution of the left-most up-spin. In the $\Delta=0$ case, the distribution $F_2$ arises. In the $\Delta \neq 0$ case, we propose a conjectural series expansion type formula based on a saddle point analysis. The conjectural formula turns out to be a Fredholm series expansion in the $\Delta \rightarrow 0$ limit and recovers the result for $\Delta = 0$., Comment: 41 pages, 4 figures. Version 2 corrects a typo in notation for equation (98)

Application of isotropic pressure or uniaxial strain alters the elastic properties of materials; sufficiently large strains can drive structural transformations. Linear elasticity describes stability against infinitesimal strains, while nonlinear elasticity describes the response to finite deformations. It was previously shown that uniaxial strain along [100] drives refractory metals and alloys towards mechanical instabilities. These include an extensional instability, and a symmetry-breaking orthorhombic distortion caused by a Jahn-Teller-Peierls instability that splays the cubic lattice vectors. Here, we analyze these transitions in depth. Eigenvalues and eigenvectors of the Wallace tensor identify and classify linear instabilities in the presence of strain. We show that both instabilities are discontinuous, leading to discrete jumps in the lattice parameters. We provide physical intuition for the instabilities by analyzing the changes in first principles energy, stress, bond lengths and angles upon application of strain. Electronic band structure calculations show differential occupation of bonding and anti-bonding orbitals, driven by the changing bond lengths and leading to the structural transformations. Strain thresholds for these instabilities depend on the valence electron count., Comment: 31 pages, 17 figures, 1 table, accepted to Physical Review Materials

The Kubo-Greenwood equation, in combination with the first-principles Korringa-Kohn-Rostoker Coherent Potential Approximation (KKR-CPA) can be used to calculate the DC residual resistivity of random alloys at T = 0 K. We implemented this method in a multiple scattering theory based ab initio package, MuST, and applied it to the ab initio study of the residual resistivity of the high entropy alloy Al$_x$CoCrFeNi as a function of $x$. The calculated resistivities are compared with experimental data. We also predict the residual resistivity of refractory high entropy alloy MoNbTaV$_x$W. The calculated resistivity trends are also explained using theoretical arguments., Comment: 18 pages, 9 figures (not including Supplementary Material). The article has been submitted to Applied Physics Letters

Comparison of free energies between different phases and different compositions underlies the prediction of alloy phase diagrams. To allow direct comparison, consistent reference points for the energies or enthalpies are required, and the entropy must be placed on an absolute scale, yielding absolute free energies. Here we derive absolute free energies of liquids from ab-initio molecular dynamics (AIMD) by combining the directly simulated enthalpies with an entropy derived from simulated densities and pair correlation functions. As an example of the power of this method we calculate the phase diagrams of two binary alkali metal alloys, Li-Na and K-Na, revealing a critical point and liquid-liquid phase separation in the former case, and a deep eutectic in the latter. Good agreement with experimental data demonstrates the power of this simple method.

Energy efficiency is motivating the search for new high-temperature metals. Some new body-centered-cubic random multicomponent "high entropy alloys (HEAs)" based on refractory elements (Cr-Mo-Nb-Ta-V-W-Hf-Ti-Zr) possess exceptional strengths at high temperatures but the physical origins of this outstanding behavior are not known. Here we show, using integrated neutron-diffraction (ND), high-resolution transmission electron microscopy (HRTEM), and theory, that the high strength and strength retention of a NbTaVTi alloy and a new high-strength/low-density CrMoNbV alloy are attributable to edge dislocations. This is surprising because plastic-flow in BCC elemental metals and dilute alloys is universally accepted to be controlled by screw dislocations. We use the insight and theory to perform a computationally-guided search over $10^7$ BCC HEAs and identify over $10^6$ possible ultra-strong high-temperature alloy compositions for future exploration.

Solid solutions occur when multiple chemical species share sites of a common crystal lattice. Although the single site occupation is random, chemical interaction preferences bias the occupation probabilities of neighboring sites, and this bias reduced the entropy of mixing below its ideal value. Sufficiently strong bias leads to symmetry-breaking phase transitions. We apply the cluster variation method to explore solid solutions on body centered cubic lattices in the context of two specific compounds that exhibit opposite ordering trends. Employing density functional theory to model the energetics, we show that CuZn exhibits an order-disorder transition to the CsCl prototype structure, while AlLi instead takes the NaTl prototype structure, and we evaluate their temperature-dependent order parameters, correlations and entropies., Comment: 13 pages, 8 figures

The single-site Korringa-Kohn-Rostoker Coherent Potential Approximation (KKR-CPA) ignores short range ordering present in disordered metallic systems. In this paper, we establish a new technique to fix this shortcoming by embedding an averaged cluster that displays chemical short range order (SRO). The degree of SRO can be tuned by externally defined order parameters. This averaged cluster can be embedded in the single site CPA medium, or a self-consistently obtained effective medium that contains SRO information. The validity of this method is demonstrated by applying it to two alloy systems - the CuZn body centered cubic (BCC) solid solution, and AlCrTiV, a four-element BCC high entropy alloy. A comparison between the non-self-consistent and self-consistent modes is also provided for the two above mentioned systems. We make the code available on the internet. Planned extensions to this work are discussed., Comment: 9 pages, 7 figures; minor revisions based on reviewer comments; added a section on self-consistent CA-CPA

Icosahedral quasicrystals spontaneously form from the melt in simulations of Al--Cu--Fe alloys. We model the interatomic interactions using oscillating pair potentials tuned to the specific alloy system based on a database of density functional theory (DFT)-derived energies and forces. Favored interatomic separations align with the geometry of icosahedral motifs that overlap to create face-centered icosahedral order on a hierarchy of length scales. Molecular dynamics simulations, supplemented with Monte Carlo steps to swap chemical species, efficiently sample the configuration space of our models, which reach up to 9846 atoms. Exchanging temperatures of independent trajectories (replica exchange) allows us to achieve thermal equilibrium at low temperatures. By optimizing structure and composition we create structures whose DFT energies reach to within $\sim$2 meV/atom of the energies of competing crystal phases. Free energies obtained by adding contributions due to harmonic and anharmonic vibrations, chemical substitution disorder, phasons, and electronic excitations, show that the quasicrystal becomes stable against competing phases at temperatures above 600K. The average structure can be described succinctly as a cut through atomic surfaces in six-dimensional space that reveal specific patterns of preferred chemical occupancy. Atomic surface regions of mixed chemical occupation demonstrate the proliferation of phason fluctuations, which can be observed in real space through the formation, dissolution and reformation of large scale icosahedral motifs -- a picture that is hidden from diffraction refinements due to averaging over the disorder and consequent loss of information concerning occupancy correlations.

The structure of the SiC(000-1) surface, the C-face of the {0001} SiC surfaces, is studied as a function of temperature and of pressure in a gaseous environment of disilane (Si2H6). Various surface reconstructions are observed, both with and without the presence of an overlying graphene layer (which spontaneously forms at sufficiently high temperatures). Based on cross-sectional scanning transmission electron microscopy measurements, the interface structure that forms in the presence of the graphene is found to contain 1.4 - 1.7 monolayers (ML) of Si, a somewhat counter-intuitive result since, when the graphene forms, the system is actually under C-rich conditions. Using ab initio thermodynamics, it is demonstrated that there exists a class of Si-rich surfaces containing about 1.3 ML of Si that are stable on the surface (even under C-rich conditions) at temperatures above about 400 K. The structures that thus form consist of Si adatoms atop a Si adlayer on the C-face of SiC, with or without the presence of overlying graphene., Comment: 19 pages with 8 figures, along with 6 pages of Supplemental Material; v2 adds refs. 56,59-62 all being prior work, along with several paragraphs in the Discussion or Summary sections of the manuscript discussing this prior work; v3 corrects several typographical errors

Experiments on rare-earth filled skutterudites demonstrate an intriguing array of thermodynamic, transport and superconducting properties, and bring to fore theoretical challenges posed by f-electron systems. First principle calculations based density functional theory and its extensions for strongly correlated systems such as the Hubbard U correction, provide valuable information about electronic structure that can be used to understand experiments. We present a comprehensive study of the electronic structure and Fermi surface of a series of rare earth filled skutterudites, RPt4Ge12 (where R = La, Ce, Pr), aimed at shedding light on: consequences of progressive increase of f-orbital occupancy in the series; the effects of the Hubbard parameter U; the Fermi surfaces, band structures and densities of states. The calculated Fermi surfaces may be relevant to the question of multi-band versus single-band superconductivity. Computed densities of states qualitatively explain the available resonant photoemission spectroscopy experiments, and (together with available specific heat measurements) provide estimates of the effective masses. We also show the existence of pseudogaps in the total density of states which may be relevant for the thermoelectric properties of these systems., Comment: 17 pages, 11 figures

The quantum spin Hall (QSH) state was recently demonstrated in monolayers of the transition metal dichalcogenide 1T'-WTe$_2$ and is characterized by a band gap in the two-dimensional (2D) interior and helical one-dimensional (1D) edge states. Inducing superconductivity in the helical edge states would result in a 1D topological superconductor, a highly sought-after state of matter. In the present study, we use a novel dry-transfer flip technique to place atomically-thin layers of WTe$_2$ on a van der Waals superconductor, NbSe$_2$. Using scanning tunneling microscopy and spectroscopy (STM/STS), we demonstrate atomically clean surfaces and interfaces and the presence of a proximity-induced superconducting gap in the WTe$_2$ for thicknesses from a monolayer up to 7 crystalline layers. At the edge of the WTe$_2$ monolayer, we show that the superconducting gap coexists with the characteristic spectroscopic signature of the QSH edge state. Taken together, these observations provide conclusive evidence for proximity-induced superconductivity in the QSH edge state in WTe$_2$, a crucial step towards realizing 1D topological superconductivity and Majorana bound states in this van der Waals material platform.