Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First principles”

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 1,009 records · Page 56

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ab initio studies of the impact of the Debye-Waller factor on the structural and dynamical properties of amorphous semiconductors: The case of a-Si

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (a-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in a-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts – originating from isolated and vacancy-induced clustered defects – and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in a-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of a-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of a-Si at low temperatures.

36 MATERIALS SCIENCE↗

Electronic structure and lattice dynamics of 1⁢𝑇−VSe 2 : Origin of the three-dimensional charge density wave

To characterize in detail the charge density wave (CDW) transition of 1⁢𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇⁢ CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇⁢ CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Molecular Simulation of Lithium Carbonate Reactive Vapor–Liquid Equilibria Using a Deep Potential Model

We developed a first-principles machine learning model for the reactive vapor–liquid phase behavior of molten Li 2 CO 3 . The model was trained on ab initio electronic density functional theory data using the Deep Potential (DP) methodology, and its accuracy was evaluated by comparing model predictions of density and viscosity to experimental measurements. Direct coexistence simulations with the DP model over time scales of tens of nanoseconds were used to observe equilibrium dissociation of Li 2 CO 3 into CO 2 residing primarily in the vapor phase and Li 2 O which remains dissolved in the liquid. The simulations covered a range of temperatures, overall system sizes, and vapor-to-liquid volume ratios. Results were analyzed in terms of the observed chemical composition of the liquid and vapor phases, product structure, and CO 2 partial pressures. In addition, we calculated equilibrium constants for the dissociation reaction by assuming ideal-solution behavior for the liquid. As expected on the basis of thermodynamic arguments and prior experiments for this system, the observed partial pressure of CO 2 in the gas phase depends on both the temperature and the ratio of vapor to liquid volumes, while the calculated equilibrium constants only depend on temperature. DP model predictions for the equilibrium constant of the reaction are generally consistent with the available experimental measurements. Furthermore, the present study establishes the validity of the DP methodology for the description of reactive, multiphase equilibria from first principles, with possible applications to many other systems of scientific and technological interest even in the absence of relevant experimental measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anisotropic magnetism and electronic structure of trigonal EuAl 2 Ge 2 single crystals

Understanding the interplay between magnetic and electronic degrees of freedom is of profound recent interest in different Eu-based magnetic topological materials. Here, we studied the magnetic and electronic properties of the layered Zintl-phase compound EuAl 2 Ge 2 crystallizing in the trigonal CaAl 2 Si 2 -type structure. We report zero-field neutron diffraction, temperature T - and magnetic-field H-dependent magnetic susceptibility χ(T, H), isothermal magnetization M(T, H), heat capacity C p (T, H), and electrical resistivity ρ(T, H) measurements, together with T-dependent angle-resolved photoemission spectroscopy (ARPES) measurements complemented with first-principle calculations. EuAl 2 Ge 2 undergoes second-order A-type antiferromagnetic (AFM) ordering below T N = 27.5(5) K, with the Eu moments (Eu 2+ , S = 7/2) aligned ferromagnetically in the ab plane while these layers are stacked antiferromagnetically along the c axis. The critical fields at which all moments become parallel to the field are 37.5(5) and 52.5(5) kOe for H ∥ ab and H ∥ c, respectively. The H = 0 magnetic structure consists of trigonal AFM domains associated with ab-plane magnetic anisotropy and a field-induced reorientation of the Eu spins in the domains is also evident at T = 2 K below the critical field H c1 = 2.5(1) kOe. The ρ(T) measurements reveal metallic behavior transforming into a slight resistivity increase on cooling towards T N . A pronounced loss of spin-disorder scattering is observed below T N . The ARPES results show that EuAl 2 Ge 2 is metallic both above and below T N , and the Fermi surface is anisotropic with two hole pockets at the zone center and one small electron pocket at each M point. In the AFM phase, we directly observe folded bands in ARPES due to the doubling of the magnetic unit cell along the c axis with an enhancement of quasiparticle weight due to the complex change in the coupling between the magnetic moments and itinerant electrons on cooling below T N . The observed electronic structure is well reproduced by first-principle calculations, which also predict the presence of nontrivial electronic states near the Fermi level in the AFM phase with Z 2 topological numbers 1;(000)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Impact of Sr-Containing Secondary Phases on Oxide Conductivity in Solid-Oxide Electrolyzer Cells

Solid-oxide electrolyzer cells (SOECs) based on a yttria-stabilized zirconia (YSZ) oxide electrolyte produce hydrogen from water with the assistance of excess thermal energy; however, Sr diffusion within the Gd-doped CeO 2 (GDC) barrier layer during processing or operation can lead to the formation of unwanted secondary phases such as SrO and SrZrO 3 . Here, to establish and compare the degree of impact of these phases on SOEC performance, we conduct first-principles calculations to study their bulk oxide conductivities and compare them to that of the YSZ electrolyte. We find that SrO has a low conductivity arising from the poor mobility and low concentration of mobile oxygen vacancies, and its presence in SOECs should therefore be avoided. SrZrO 3 also has a lower oxide conductivity than YSZ; however, this discrepancy is primarily due to lower vacancy concentrations rather than low mobility. We find that sufficient levels of Y-doping on the Zr site can increase oxygen vacancy concentrations in SrZrO 3 to achieve an oxide ionic conductivity on par with that of YSZ, thereby mitigating any potential deleterious effect on transport performance. Energy-dispersive X-ray spectroscopy confirms that Y is the most common minority element present in SrZrO 3 forming near the GDC–YSZ interface, alleviating concerns regarding the impact of SrZrO 3 on device performance. These results from our combined computational–experimental analysis can inform future engineering strategies designed to limit the detrimental effects of Sr-induced secondary phase formation on SOEC performance.

08 HYDROGEN↗

(NH 3 (CH 2 ) 7 NH 3 ) 2 Sn 3 I 10 , a Vacancy-Ordered Three-Dimensional Tin(II) Perovskite-Derived Semiconductor

Ordering vacancies in hybrid Sn(II) halide semiconductors provides a strategy for preventing uncontrolled oxidation and formation of mobile holes. In this study, we report the structure and optical and electronic properties of (NH 3 (CH 2 ) 7 NH 3 ) 2 Sn 3 I 10 , a vacancy-ordered perovskite derivative with three-dimensional inorganic connectivity. The crystal structure resembles that of a Dion–Jacobson layered perovskite derivative, but with [SnI 5 ] square pyramids bridging the layers. UV–vis diffuse reflectance spectroscopy reveals a sharp onset of light absorption at 1.86(1) eV with the photoluminescence emission maximum at 1.90(1) eV. However, the maximum excitation occurs from 3.42 to 3.81 eV (325 to 370 nm), revealing a significant Stokes shift of 1.3 eV. The electronic properties determined from dark and time-resolved microwave conductivity measurements reveal a minimum carrier mobility of 4.3 × 10 –2 cm 2 V –1 s –1 and a maximum carrier density of 5.96 × 10 16 cm –3 , a uniquely low value for a hybrid Sn(II) halide semiconductor. The transport behavior in combination with first-principles calculations of the electronic band structure and dielectric permittivity suggest polaron-mediated electronic transport, yet the photogenerated carriers have a fast and fluence-dependent nonradiative recombination rate, suggestive of localized “defect-like” states at the band edge. The observed photoluminescence is most consistent with single-ion-like behavior of an asymmetric Sn(II) environment. Together, these results suggest that defect ordering presents a strategy for the reduction of mobile charge carriers at equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HfTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-titanium (Hf-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗

NbTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-titanium (Nb-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP), https://github.com/ORNL/ScalableWorkflow_VASP_Calculations. Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗

Localized High Concentration Electrolytes for High Voltage Lithium-Metal Batteries: Correlation between the Electrolyte Composition and Its Reductive/Oxidative Stability

We demonstrate a first-principles screening methodology as an effective tool to explore electrolyte formulations for the new generation of high energy density rechargeable batteries. We study the liquid structure and electronic properties in dilute electrolytes, high concentration electrolytes (HCE), and localized high concentration electrolytes (LHCE), with focus on electrolyte formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethyl carbonate (DMC), and bis(2,2,2-trifluoroethyl) ether (BTFE) as a diluent. We describe the solvation complexes in the dilute electrolyte and explore structural changes triggered by the increase in lithium salt concentration for HCEs and the diluent effects in LHCEs. In HCE formulations, there is a 4-fold coordination environment of lithium-ions as in the dilute electrolyte, but the number of lithium-ion interactions with O atoms from FSI- anions dominates. In these solutions, the ability of the FSI- anions to interact with multiple lithium-ions allows complex 3D network formation and influences the reductive/oxidative behavior of the electrolyte. Interestingly, in LHCEs, the BTFE diluent molecules do not change the 3D solution structure when diluting the HCE formulation from 5.49 to 3.83 M. However, there is a composition threshold where the structural and electronic behavior may change. We show that diluting the HCE electrolyte with BTFE down to 1.77 M breaks the three-dimensional solution structure into an island-like solvation complex. We relate these structural changes to the electronic properties of the electrolytes finding a causal relationship between the reductive/oxidative behavior and the lithium-oxygen interaction mechanisms in the solvated complexes. The coordination with lithium-ions lowers the electrolyte LUMO and HOMO levels: the higher is the number of interactions with lithium-ions, the more likely the solvent molecule, FSI- anion, or diluent molecule is to be reduced and the less likely it is to become oxidized. The evolution of the solvated ion structure in HCE and LHCE suggests a close connection to a corresponding change in the lithium-ion transport mechanisms for these electrolytes.

anions, molecules, oxygen, solvents, electrolytes↗

A Multi-Scale Computational Platform for Predictive Modeling of Corrosion in Al-Steel Joints (Final Report)

The research team proposed to develop innovative multi-scale models to predict corrosion and the resulting mechanical performances in aluminum-steel joints. The methods of joining considered are resistance spot welding, self-piercing riveting, and rivet-welding, all suitable for mass production applications. The multi-scale models integrate high throughput first-principle calculations based on density functional theory (DFT), high throughput calculation of phase diagrams (CALPHAD) modeling, and finite element method (FEM) simulations. These models are to be validated through laboratory experiments. Furthermore, the models are available as open source so as to enable scientists and engineers in the community to adapt and contribute to the development and application. The approaches rely on the research team’s extensive experience on the prediction of properties of individual phases at finite temperatures and variable compositions through DFT calculations, and our broad expertise on dissimilar material joining and their corrosion. The proposed computational framework enables high throughput computations for improved predictions of corrosion and the associated mechanical performance in dissimilar material joints, resulting in significant reduction in computational time needed by the current state-of-the-art methods. With the participation of researchers from three universities, an auto manufacturer, two manufacturing technology/equipment suppliers, and a software developer/vendor, the interdisciplinary research team applies the technical development on both phase-based modeling and laboratory experiments into the automobile body joining processes for validation and technology demonstration. The global cost of corrosion was estimated at about 3.4% of the global GDP in 2013. By using available corrosion control practices, it is estimated a saving between 15-35% of the cost of corrosion. In the U.S., more than $276 billion is spent repairing corrosion damage. Prediction of the corrosion and its impact on performance of the dissimilar material joints is critical for reducing the massive number of the current corrosion-based recalls for automobiles. Thus, the project goal is to develop models to enable predictive maintenance and end-of-life planning of multi-metal joints with risk of corrosion under different conditions such as exposure to high temperatures in summer and salt solutions in winter, quantified through its pH. An academia-industry consortium led by the University of Michigan and including Pennsylvania State University, University of Illinois Urbana-Champaign, University of Georgia, General Motors Company, Livermore Software Technology Corporation, and Optimal Process Technologies, LLC. created multi-scale models for prediction of corrosion in aluminum-steel joint structures such of them used in vehicle subassemblies – chassis and transmission systems. Starting from the first principle calculations, the team developed mathematical and data-driven models to predict the metallic components, which are formed during joining of two metals, for example aluminum and steel - a lightweight multilateral system which is currently used in more than 60% car bodies. These models were used for simulating chemical reactions that are happening when the joining metallic components are exposed to high temperatures and different pH values. The team was able to predict how the corrosion installs on the metallic components and how they lead to a sudden failure of components in cars. Newly developed machine learning algorithms combining Science, Technology, Engineering and Math disciplines, advanced finite element simulation and experimental validations have been integrated in a platform for prediction of the corrosion evolution and prediction the failure of joints under mechanical loadings and fatigue. Moreover, based on machine learning and inverse analysis, the team proposed solutions for designing new metallic alloys less susceptible to corrosion when joining multi-material assembles. An average of 4% error compared with experiments was achieved for the most common joints that are used in vehicle subassemblies.

36 MATERIALS SCIENCE↗

Developing single-site Pt catalysts for the preferential oxidation of CO: A surface science and first principles-guided approach

Here, we report a comprehensive study combining surface science, Density Functional Theory (DFT) calculations, and catalyst synthesis, characterization, and testing to investigate the preferential oxidation of CO in the presence of H 2 over single-site Pt 1 /Cu x O catalysts. Surface science studies show that while Pt 1 /Cu x O model surfaces enable low-temperature CO oxidation via a Mars-van Krevelen mechanism, there was no evidence for H 2 activation or oxidation. DFT-based calculations explain these results and demonstrate that the H 2 oxidation barrier is high as compared to H 2 desorption from Pt 1 /Cu x O. Inspired by these model catalyst studies, nanoporous Pt 1 /Cu x O catalysts were synthesized and demonstrated to be active and highly selective for the preferential oxidation of CO. This work highlights the potential of combined surface science, theory, and catalyst studies for identifying new catalytic materials, which in this case led to the development of a promising single-site Pt 1 /Cu x O catalyst for the preferential oxidation of CO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetism and finite-temperature effects in UZr 2 : A density functional theory analysis

The structure and the thermophysical properties of d-UZr 2 are investigated using 0 K density functional theory and ab initio molecular dynamics (AIMD). Modeling the true paramagnetic state of this intermetallic compound has been challenging using first-principles calculations. For the first time, we find that the generalized gradient approximation method without applying an on-site Coulomb interaction term (Hubbard U) can result in a ground state that is antiferromagnetic (AFM). We believe that this weak AFM ground state is the closest to the real paramagnetic state. We found that structure optimization at finite temperatures using AIMD is necessary to achieve this ground state instead of the non-magnetic and ferromagnetic states previously reported in the literature. Our findings indicate that applying the Hubbard U on uranium f-orbitals in this metallic system is unnecessary and not recommended, as it leads to a large overestimation of the volume and introduces an unphysical strong spin polarization. Finally, our approach results in atomic volume, thermal expansion, and heat capacities that have strong agreement with experiments.

36 MATERIALS SCIENCE↗

Molecular beam epitaxy and crystal structure of majority a -plane-oriented and substrate-strained Mn 3 Sn thin films grown directly on sapphire (0001)

The Kagome antiferromagnet Mn 3 Sn has garnered much attention due to the presence of exciting properties such as anomalous Hall and Nernst effects. This paper discusses the synthesis of crystalline Mn 3 Sn thin films, prepared on Al 2 O 3 (0001) substrates at 453 ± 5°C using molecular beam epitaxy. The growth is monitored in situ using reflection high energy electron diffraction and measured ex situ using x-ray diffraction, Rutherford back-scattering, and cross-sectional scanning transmission electron microscopy. Our analysis shows the in-plane lattice constants of a 1,M = 4.117 ± 0.027 Å and a 2,M = 4.943 ± 0.033 Å, which is a very unexpected result when compared to the bulk a-plane Mn 3 Sn. This indicates a strain in the film and makes it challenging to provide a straightforward explanation. In an effort to explain our results, we discuss two possible orientation relationships between the Mn 3 Sn films and the sapphire substrates. Samples prepared under these conditions appear to have smooth surfaces locally, but overall the film has a 3D island morphology. First-principles calculations provide atomic models of the Mn 3 Sn (112¯0) lattice on Al 2 O 3 (0001) high symmetry sites, indicating that the L3-R90° is the most stable configuration. A detailed discussion of the experimental data and theoretical results, as well as strain effects, is provided.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiative capture rates at deep defects from electronic structure calculations

We present a methodology to calculate radiative carrier capture coefficients at deep defects in semiconductors and insulators from first principles. Electronic structure and lattice relaxations are accurately described with hybrid density functional theory. Calculations of capture coefficients provide an additional validation of the accuracy of these functionals in dealing with localized defect states. We also discuss the validity of the Condon approximation, showing that even in the event of large lattice relaxations the approximation is accurate. We test the method on GaAs:VGa-TeAs and GaN:CN, for which reliable experiments are available, and demonstrate very good agreement with measured capture coefficients.

36 MATERIALS SCIENCE↗

Predicting the Phase Stability of Multicomponent High-Entropy Compounds

A generic method to estimate the relative feasibility of formation of high-entropy compounds in a single phase, directly from first principles, is developed. As a first step, the relative formation abilities of 56 multicomponent, AO, oxides were evaluated. These were constructed from five cation combinations chosen from A = {Ca, Co, Cu, Fe, Mg, Mn, Ni, Zn}. Candidates for multicomponent oxides are predicted from descriptors related to the enthalpy and configurational entropy obtained from the mixing enthalpies of two-component oxides. The utility of this approach is evaluated by comparing the predicted combinations with the experimentally realized entropy-stabilized oxide, (MgCoCuNiZn)O. In the second step, Monte Carlo simulations are utilized to investigate the phase composition and local ionic segregation as a function of temperature. This approach allows for the evaluation of potential secondary phases, thereby making realistic predictions of novel multicomponent compounds that can be synthesized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗