Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A molecular ground electronic state with an occupied 5g spinor—The superheavy (E125)F molecule

Fully relativistic calculations, primarily at the 4-component coupled-cluster singles and doubles with perturbative triples [CCSD(T)] level of theory with the Dirac–Coulomb (DC) Hamiltonian, have been carried out for the superheavy (E125)F molecule using large Gaussian basis sets. The electronic ground state is determined to have an [Og]8s 2 5g 1 6f 3 configuration on E125 with an Ω = 6 ground state and an 8p electron largely donated to F. A Mulliken population analysis indicates that the ground state is mainly ionic with a partial charge of +0.79 on E125 and a single sigma bond involving the F 2p and E125 8p spinors. The occupied g spinor is not involved in the bonding. With the largest basis set used in this work, the (0 K) dissociation energy was calculated at the DC-CCSD(T) level of theory to be 7.02 eV. Analogous calculations were also carried out for the E125 atom, both the neutral and its cation. The lowest energy electron configuration of E125 + , [Og]$8s$$^{2}_{1/2}$$5g$$^{1}_{7/2}$$6f$$^{3}_{5/2}$ with a J = 6 ground state, was found to be similar to that in (E125)F, while the neutral E125 atom has an [Og]$8s$$^{2}_{1/2}$$5g$$^{1}_{7/2}$$6f$$^{2}_{5/2}$$7d$$^{1}_{3/2}$$8p$$^{1}_{1/2}$ ground state electron configuration with a J = 17/2 ground state. The ionization energy (IE) of E125 is reported for the first time and is calculated to be 4.70 eV at the DC-CCSD(T) level of theory. Non-relativistic calculations were also carried out on the E125 atom and the (E125)F molecule. Here, the non-relativistic ground state of the E125 atom was calculated to have a 5g 5 ground state with an IE of just 3.4 eV. The net effect of relativity on (E125)F is to stabilize its bonding.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Reactivity of Sulfur Vacancy-Rich MoS 2 to Water Dissociation

Exposed Mo atoms on the surface of MoS 2 can catalyze certain useful chemical reactions, including the dissociation of water to produce hydrogen. However, a comprehensive understanding of water reactivity on defective MoS 2 surfaces remains elusive. Here, we use in situ near ambient pressure X-ray photoelectron spectroscopy (NAP-XPS) to investigate water dissociation reactions on MoS 2 surfaces before and after Ar + ion beam bombardment. To make the surfaces reactive, we treated them with Ar + ion beam sputtering first, which created exposed Mo sites. Using ultraviolet photoelectron spectroscopy (UPS) and density of states calculations conducted using density functional theory (DFT), we verified that stripping the surface of sulfur atoms creates metallic surface states that can catalyze water dissociation. At elevated H 2 O pressures, XPS measurements combined with DFT calculations suggested the presence of four distinct surface species from dissociated water. Specifically, we found that oxides and hydroxides are prominent at the surface, while chemisorbed and physisorbed H 2 O molecules are also present. In conclusion, this study provides new insights that reveal the prospects of surface-engineered MoS 2 as a catalyst for water dissociation.

08 HYDROGEN↗
DOE OSTI2020

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Solving the electronic structure problem for over 100000 atoms in real space

Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 ⁢H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.

0-dimensional systems↗
DOE OSTI2025

Atomic Energy Accuracy of Neural Network Potentials: Harnessing Pretraining and Transfer Learning

Machine learning-based interatomic potentials (MLIPs) have transformed the prediction of potential energy surfaces (PESs), achieving accuracy comparable to ab initio calculations. However, atomic energy predictions, often assumed to lack physical meaning, remain underexplored. In this study, we demonstrate that inaccuracies in atomic energy predictions reduce the robustness and transferability of Neural Network Potentials (NNPs) and atomic energy error can be masked in total energy predictions due to error cancellation. Here, we validate this finding using challenging configurations involving deformation and failure under tensile loading. By pretraining atomic energy predictions using empirical potentials and applying transfer learning with density functional theory (DFT) data, we achieve notable improvements in the accuracy of total energy, forces, and stress predictions. Furthermore, this approach enhances the robustness and transferability of NNPs, emphasizing the importance of atomic energy predictions in developing high-quality and reliable MLIPs.

Active Learning↗
DOE OSTI2024

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. Here, this study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure–temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mixed-Valent Copper Chalcogenides: Tuning Structures and Electronic Properties Using Multiple Anions

A series of mixed-anion copper chalcogenides have been prepared using solid-state methods including ACu 4.2 TeS 2 (A = K, Rb, Cs), which adopt the KCu 4 S 3 structure type. The mixedanion motif has an expanded sublattice, relative to KCu 4 S 3 , that can accommodate additional Cu atoms at its interstitial sites that are unoccupied in the parent structure. The variable temperature transport shows that the materials are p-type metals with carrier densities on the order of 10 21 cm -3 and room-temperature electrical conductivity as high as 4000 S cm -1 . Band structures calculated using density functional theory corroborate the experimental data and indicate that the interstitial Cu atoms lower the carrier concentration and increase the Fermi level of the materials. The layered structure has Te and S atoms occupying unique sites within the ACu 4.2 TeS 2 structure, where relatively hard S 2- anions prefer Wyckoff positions where they can form energetically favorable acid-base interactions with hard alkali cations. The phenomenon is observed in the related K 3 Cu 8 Te 2 S 4 system, which also has a fully ordered layered structure. Here, we believe that the report provides new chemical guidelines for targeting ordered multianion structures, as well as a unique method for tuning the electronic properties of metallic chalcogenides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Coupled Experimental–Theoretical Characterization of a Carbon Electrode in Vanadium Redox Flow Batteries using X-ray Absorption Spectroscopy

Vanadium redox flow batteries (VRFBs) have emerged as promising solutions for stationary grid energy storage due to their high efficiency, scalability, safety, near room-temperature operation conditions, and the ability to independently size power and energy capacities. The performance of VRFBs heavily relies on the redox couple reactions of V 2+ /V 3+ and VO 2+ /VO 2 + on carbon electrodes. Therefore, a thorough understanding of the surface functionality of carbon electrodes and their propensity for degradation during electrochemical cycles is crucial for designing VRFBs with extended lifespans. In this study, we present a coupled experimental–theoretical approach based on carbon K edge X-ray absorption spectroscopy (XAS) to characterize carbon electrodes prepared under different conditions and identify relevant functional groups that contribute to unique spectroscopic features. Atomic models were created to represent functional groups, such as hydroxyl, carboxyl, methyl, and aldehyde, bonded to carbon atoms in either sp 2 or sp 3 environments. The interactions between functionalized carbon and various solvated vanadium complexes were modeled using density functional theory. A library of carbon K-edge XAS spectra was generated for distinct carbon atoms in different functional groups, both before and after interacting with solvated vanadium complexes. Here we demonstrate how these simulated spectra can be used to deconvolve ex situ experimental spectra measured from carbon electrodes and to track changes in the electrode composition following immersion in different electrolytes or extended cycling within a functional VRFB. By doing so, we identify the active species present on the carbon electrodes, which play a crucial role in determining their electrochemical performance.

25 ENERGY STORAGE↗
DOE OSTI2024

Generalized rate theory for void and bubble swelling and its application to δ -plutonium

A rate theory for void and bubble swelling is derived that allows both vacancies and self-interstitial atoms to be generated by thermal activation at all sinks. In addition, they can also be produced by displacement damage from external and internal radiation. This generalized rate theory (GRT) is applied to swelling of gallium-stabilized δ-plutonium in which α-decay causes the displacement damage. Since the helium atoms produced also become trapped in vacancies, a distinction is made between empty and occupied vacancies. The growth of helium bubbles observed by transmission electron microscopy in weapons-grade and in material enriched with Pu238 is analyzed, using different values for the formation energy of self-interstitial atoms (SIA) and two different sets of relaxation volumes for the vacancy and for the SIA. One set allows preferential capture of SIA at dislocations, while the other set gives equal preference to both vacancy and SIA. It is found that the helium bubble diameters observed are in better agreement with GRT predictions if no preferential capture occurs at dislocations. Therefore, helium bubbles in δ-plutonium will not evolve into voids. Furthermore, the helium density within the bubbles remains sufficiently high to cause thermal emission of SIA. Based on a helium density between two to three helium atoms per vacant site, the sum of formation and migration energies must be around 2.0 eV for SIA in δ-plutonium.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

The role of atomic carbon in directing electrochemical CO (2) reduction to multicarbon products

Electrochemical reduction of carbon-dioxide/carbon-monoxide (CO (2) R) to fuels and chemicals presents an attractive approach for sustainable chemical synthesis, but it also poses a serious challenge in catalysis. Understanding the key aspects that guide CO (2) R towards value-added multicarbon (C 2+ ) products is imperative in designing an efficient catalyst. Herein, we identify the critical steps toward C 2 products on copper through a combination of energetics from density functional theory and micro-kinetic modeling. We elucidate the importance of atomic carbon in directing C 2+ selectivity and how it introduces surface structural sensitivity on copper catalysts. Overall, this insight enables us to propose two simple thermodynamic descriptors that effectively identify C 2+ selectivity on metal catalysts beyond copper and hence it defines an intelligible protocol to screen for materials that selectively catalyze CO (2) to C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗
DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
DOE OSTI2021

Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal–Organic Frameworks

Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO 2 . Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal–Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. In conclusion, a strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Tunable tunnel barriers in a semiconductor via ionization of individual atoms

individual adatoms deposited on an InSb(110) surface. The adatoms can be reproducibly dropped off from the STM tip by voltage pulses, and impact tunneling into the surface by up to ~100×. The spatial extent and magnitude of the tunneling effect are widely tunable by imaging conditions such as bias voltage, set current and photoillumination. We attribute the effect to occupation of a (+/0) charge transition level, and switching of the associated adatom-induced band bending. The effect in STM topographic images is well reproduced by transport modeling of filling and emptying rates as a function of the tip position. STM atomic contrast and tunneling spectra are in good agreement with density functional theory calculations for In adatoms. The adatom ionization effect can extend to distances greater than 50 nm away, which we attribute to the low concentration and low binding energy of the residual donors in the undoped InSb crystal. These studies demonstrate how individual atoms can be used to sensitively control current flow in nanoscale devices.

Physics↗
DOE OSTI2023

Modeling–Experiment–Theory Analysis of Reactions Initiated from Cl + Methyl Formate

Methyl formate (MF; CH 3 OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF—CH 3 OCO and CH 2 OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl 2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl 2 , HCl, CO 2 , CH 3 , CH 3 Cl, CH 2 O, and CH 2 ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400–750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH 3 OCO (R1a) and Cl + MF → HCl + CH 2 OCHO (R1b): k 1a,theory = 6.71 × 10 –15 ·T 1.14 ·exp(—606/T) cm 3 /molecule·s; k 1b,theory = 4.67 × 10 –18 ·T 2.21 ·exp(—245/T) cm 3 /molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH 3 OCO and CH 2 OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH 3 + CO 2 (R3) and CH 2 O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k 3,theory (P = ∞) = 2.94 × 10 9 ·T 1.21 ·exp(—6209/T) s –1 and k 5,theory (P = ∞) = 8.45 × 10 8 ·T 1.39 ·exp(—15132/T) s –1 over T = 300–1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH 2 OCHO + Cl, CH 2 OCHO + Cl 2 , and CH 3 + Cl 2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Three-body resonances in pionless effective field theory

We investigate the appearance of resonances in three-body systems using pionless effective field theory at leading order with two complementary methods. The Faddeev equation is analytically continued to the unphysical sheet adjacent to the positive real energy axis using a contour rotation. We consider both the three-boson system and the three-neutron system. For the former, we calculate the trajectory of Borromean three-body Efimov states turning into resonances as they cross the three-body threshold. For the latter, we find no sign of three-body resonances or virtual states at leading order. Furthermore, this result is validated by exploring the level structure of three-body states in a finite volume approach.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Accurate prediction of short-range order and its effect on thermodynamic, structural, and electronic properties of disordered alloys: exemplified in archetypical Cu 3 Au

Electronic-structure methods based on density-functional theory (DFT) were used to quantify the effect of chemical short-range order (SRO) on thermodynamic, structural, and electronic properties of archetypal face-centered-cubic (fcc) Cu3Au alloy. We showed that SRO can be tuned to alter bonding and lattice dynamics (i.e., phonons) and detail how these properties are changed with SRO. Thermodynamically favorable SRO significantly improved the phase stability of fcc Cu3Au from -0.0343 eV-atom -1 to –0.0682 eV-atom -1 . We used our DFT-based linear-response theory to predict SRO and its electronic origin, and accurately estimate the observed transition temperature, ordering instability (L1 2 ), and Warren-Cowley SRO parameters, in agreement with experiments. The accurate prediction of real-space SRO gives an edge over computationally and resource intensive approaches such as monte-carlo methods or experiments, which will enable large scale molecular dynamic simulations by providing supercells with optimized SRO. Here we also analyzed phonon dispersion and estimated the vibrational entropy change (from 9kB at 300 K to 6kB at 100 K) in fcc Cu3Au. We established from SRO analysis that exclusion of chemical interactions may lead to a skewed view of true properties in chemically complex alloys. The first-principles methods described in this work are generally applicable to any arbitrary solid-solution alloys, including multi-principal-element alloys, therefore, holds promise for designing technologically useful materials.

36 MATERIALS SCIENCE↗