Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Ghost states and surface structures of the charge density wave kagome metal ScV 6 Sn 6

In this work, we investigate the high-temperature phase of the kagome metal ScV 6 Sn 6 using scanning tunneling microscopy/spectroscopy (STM/S) and density functional theory calculations. STM topographic images of the cleaved sample reveal two distinct surface terminations: flat islands with Sn termination and trenches terminated by kagome layers with Sn as the outermost atomic layer. STS measurements on the Sn-terminated and kagome-terminated surfaces show significant differences, in particular the presence of large density of states near the Fermi level in the former case. Our first-principles calculations reveal that the charge density on the kagome-terminated surface gives rise to “ghost states” which show intensity away from surface atoms, arising due to hybridization of orbitals above the surface. These states can obscure the intrinsic properties of the surface, potentially leading to misattribution of the surface termination. This underscores the need for careful interpretation in STM studies, especially when discerning surface states of localized states. Understanding the surface structure of this versatile quantum material provides essential information for interpreting surface-sensitive experiments, tailoring material properties, engineering interfaces, and controlling stability and reactivity. This knowledge paves the way for further exploration and potential applications of kagome lattice materials in various fields, including quantum computing, topological physics, and advanced electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Magnetic-field-induced nontrivial electronic state in the Kondo-lattice semimetal CeSb

Synergic effect of electronic correlation and spin-orbit coupling is an emerging topic in topological materials. Central to this rapidly developing area are the prototypes of strongly correlated heavy-fermion systems. Recently, some Ce-based compounds are proposed to host intriguing topological nature, among which the electronic properties of CeSb are still under debate. Herein, we report a comprehensive study combining magnetic and electronic transport measurements, and electronic band-structure calculations of this compound to identify its topological nature. Quantum oscillations are clearly observed in both magnetization and magnetoresistance at high fields, from which one pocket with a nontrivial Berry phase is recognized. Angular-dependent magnetoresistance shows that this pocket is elongated in nature and corresponds to the electron pocket as observed in LaBi. Nontrivial electronic structure of CeSb is further confirmed by first-principle calculations, which arises from spin splitting in the fully polarized ferromagnetic state. These features indicate that magnetic field can induce nontrivial topological electronic states in this prototypical Kondo semimetal.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stable Cuprous Hydroxide Nanostructures by Organic Ligand Functionalization

Abstract Copper compounds have been extensively investigated for diverse applications. However, studies of cuprous hydroxide (CuOH) have been scarce due to structural metastability. Herein, a facile, wet‐chemistry procedure is reported for the preparation of stable CuOH nanostructures via deliberate functionalization with select organic ligands, such as acetylene and mercapto derivatives. The resulting nanostructures are found to exhibit a nanoribbon morphology consisting of small nanocrystals embedded within a largely amorphous nanosheet‐like scaffold. The acetylene derivatives are found to anchor onto the CuOH forming CuC linkages, whereas CuS interfacial bonds are formed with the mercapto ligands. Effective electronic coupling occurs at the ligand‐core interface in the former, in contrast to mostly non‐conjugated interfacial bonds in the latter, as manifested in spectroscopic measurements and confirmed in theoretical studies based on first principles calculations. Notably, the acetylene‐capped CuOH nanostructures exhibit markedly enhanced photodynamic activity in the inhibition of bacteria growth, as compared to the mercapto‐capped counterparts due to a reduced material bandgap and effective photocatalytic generation of reactive oxygen species. Results from this study demonstrate that deliberate structural engineering with select organic ligands is an effective strategy in the stabilization and functionalization of CuOH nanostructures, a critical first step in exploring their diverse applications.

Liu, Qiming↗
DOE OSTI2020

Ca 2 MnO 4 structural path: Following the negative thermal expansion at the local scale

The oxygen octahedral rotations in Ca 2 MnO 4 , the first member of the CaO(CaMnO 3 )n Ruddlesden-Popper family, is probed through a set of complementary techniques, including temperature-dependent neutron and x-ray diffraction, combined with local probe studies and ab initio calculations. Here we demonstrate the enhancement of the uniaxial negative thermal expansion coefficient from -1.26 ± 0.25 to -21 ± 1.8 ppm/K at the second order I4 1 /acd to I4/mmm structural phase transition, providing direct evidence for the corkscrew atomic mechanism. We establish, also, that the predicted I4/mmm high symmetry is attained around 1050 K. At lower temperatures, within the 10–1000 K temperature range, our first-principles calculations and detailed analysis of the Ca local environment reveals that the reported Aba2 structural phase, coexisting with the I4 1 /acd one, cannot describe correctly this compound. On the other hand, our data allow for the coexistence of the locally identical I4 1 /acd and Acam structural phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

𝜂 and 𝜂′ Production in 𝐽/𝜓 Radiative Decays from Quantum Chromodynamics

We present a first principles calculation within quantum chromodynamics (QCD) of the radiative decays of the 𝐽/𝜓 into the light pseudoscalar mesons 𝜂 and 𝜂′. Within a lattice computation we obtain the transition form factors as a function of photon virtuality from the timelike region, accessible experimentally via the “Dalitz” decay 𝐽/𝜓 →𝑒 + ⁢𝑒 − ⁢𝜂 (′) , through to the real photon point corresponding to 𝐽/𝜓 →𝛾⁢𝜂 (′) . This is the first calculation in lattice QCD with two (heavier than physical) degenerate flavors of light quark and a heavier strange quark, in which the 𝜂′ appears as the first excited state with pseudoscalar isoscalar quantum numbers. We access it reliably by using variationally optimized operators , the use of which also improves the purity of the 𝐽/𝜓 and 𝜂 signals, reducing systematic uncertainties. High quality results at a large number of kinematic points are obtained in a typically noisy disconnected process by using a novel correlator averaging procedure. Our results show the expected enhanced production of the 𝜂′ over the 𝜂 in this process, and suggest that the demonstrated lattice technology is suitable for future calculations considering processes in which light meson resonances are produced.

form factors↗
DOE OSTI2020

Piezoelectric Electron-Phonon Interaction from Ab Initio Dynamical Quadrupoles: Impact on Charge Transport in Wurtzite GaN

First-principles calculations of $\textit{e}$-ph interactions are becoming a pillar of electronic structure theory. However, the current approach is incomplete. The piezoelectric (PE) $\textit{e}$-ph interaction, a long-range scattering mechanism due to acoustic phonons in noncentrosymmetric polar materials, is not accurately described at present. Current calculations include short-range $\textit{e}$-ph interactions (obtained by interpolation) and the dipolelike Frölich long-range coupling in polar materials, but lack important quadrupole effects for acoustic modes and PE materials. Here we derive and compute the long-range $\textit{e}$-ph interaction due to dynamical quadrupoles, and apply this framework to investigate $\textit{e}$-ph interactions and the carrier mobility in the PE material wurtzite GaN. In this work, we show that the quadrupole contribution is essential to obtain accurate $\textit{e}$-ph matrix elements for acoustic modes and to compute PE scattering. Our work resolves the outstanding problem of correctly computing $\textit{e}$-ph interactions for acoustic modes from first principles, and enables studies of $\textit{e}$-ph coupling and charge transport in PE materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole–dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives. The embedding-potential parameters used to describe the environment in the PE model are derived through first-principles calculations, thus allowing a wide variety of systems to be modeled, including solvents, proteins, and other large and complex molecular environments. Here, we present proof-of-principle calculations of IR and Raman spectra of acetone in different solvents. This work is an important step toward calculating accurate vibrational spectra of molecules embedded in realistic environments.

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

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗
NASA NTRS2021

Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures

Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.

Thermodynamics↗
DOE OSTI2024

Derivation of low-energy Hamiltonians for heavy-fermion materials

Here, by utilizing a multiorbital periodic Anderson model with parameters obtained from ab initio band structure calculations, combined with degenerate perturbation theory, we derive effective Kondo-Heisenberg and spin Hamiltonians that capture the interaction among the effective magnetic moments. This derivation encompasses fluctuations via both nonmagnetic 4⁢𝑓 0 and magnetic 4⁢𝑓 2 virtual states, and its accuracy is confirmed through comparison with experimental data obtained from CeIn 3 . The significant agreement observed between experimental results and theoretical predictions underscores the potential of deriving minimal models from first-principles calculations for achieving a quantitative description of 4⁢𝑓 materials. Moreover, our microscopic derivation unveils the underlying origin of anisotropy in the exchange interaction between Kramers doublets, shedding light on the conditions under which this anisotropy may be weak compared to the isotropic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Intrinsic reactivity and dissolution characteristics of tetracalcium aluminoferrite

Tetracalcium aluminoferrite is one of the main phases of Portland cement, the dissolution characteristics of which however have not been fully revealed. This work adopted first-principles calculations and Mössbauer spectroscopy as well as X-ray diffraction analysis to clarify the intrinsic reactivity of different species of tetracalcium aluminoferrite crystals. Ca ions are the most reactive species of tetracalcium aluminoferrite and Al ions are slightly more reactive than Fe ions. Moreover, the tetrahedral Al and Fe ions are more reactive than the octahedral Al and Fe ions. These indicate dissolution characteristics of tetracalcium aluminoferrite crystals, namely, Ca ions dissolve first and the tetrahedral species dissolve faster than the octahedral species. After the initial dissolution, the dissolved ions nucleate and grow on the defective surface of C{sub 4}AF and the C{sub 4}AF residues are fully hydroxylated. The further hydration thereafter becomes the complex dissolution with the nucleation and growth of hydration products.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA's exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA's Entry Systems Modeling (ESM) Project, which serves as the Agency's only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specific models for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures.

Barnhardt, Michael D.↗
NASA NTRS2021

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA’s exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA’s Entry Systems Modeling (ESM) Project, which serves as the Agency’s only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specificmodels for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures. Results for both configurations have demonstrated the ability to land within the 50-meter precision requirement demanded by Mars human exploration missions, while also reducing propellant requirements by enabling more efficient control through entry and descent. Ongoing research in GN&C is developing mechanical specifications for the systems and establishing engineering feasibility.

aerothermodynamics↗
DOE OSTI2022

Experimental and theoretical study of hole scattering in RF sputtered p-type Cu 2 O thin films

One of the key materials of interest for p-type oxide semiconductor thin film electronics is cuprous oxide (Cu 2 O), due to its relatively high hole mobility. Here, in this work, we use experiments, analytical models, and density functional theory calculations to study the scattering mechanisms that determine Hall mobility in two Cu 2 O samples. First, we examine a polycrystalline Cu 2 O thin film deposited by RF magnetron sputtering, and second, a single-crystalline Cu2O bulk substrate. Temperature-dependent Hall measurements indicate that neutral impurity and grain boundary scattering are dominant for the polycrystalline Cu 2 O thin film, while phonon scattering is dominant for single-crystalline Cu 2 O. Our first-principles calculations show that the room-temperature intrinsic hole mobility of Cu 2 O is 106 cm 2 V –1 s –1 , indicating the great promise of the material for p-type electronic devices. This intrinsic mobility is limited by phonon scattering, with the most dominant scattering modes having phonon energies of 88.4 and 17.1 meV. These results indicate that the key pathways to increase the hole mobility in Cu 2 O thin films are by reducing the impurity concentration and by increasing grain size. Our work thus sets the stage for the future development of high performance Cu 2 O-based p-type thin film transistors.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Theoretical studies of spectroscopic problems of importance for atmospheric radiation measurements

Many of the instruments used to deduce the physical parameters of the Earth's atmosphere necessary for climate studies or for pollution monitoring (for instance, temperature versus pressure or number densities of trace molecules) rely on the existence of accurate spectroscopic data and an understanding of the physical processes responsible for the absorption or emission of radiation. During the summer, research was either continued or begun on three distinct problems: (1) an improved theoretical framework for the calculation of the far-wing absorption of allowed spectral lines; (2) a refinement of the calculation of the collision-induced fundamental spectrum of N2; and (3) an investigation of possible line-mixing effects in the fundamental spectrum of CH4. Progress in these three areas is summarized below. During the past few years, we have developed a theoretical framework for the calculation of the absorption of radiation by the far wings of spectral lines. Such absorption due to water vapor plays a crucial role in the greenhouse effect as well as limiting the retrieval of temperature profiles from satellite data. Several improvements in the theory have been made and the results are being prepared for publication. Last year we published results for the theoretical calculation of the absorption of radiation due to the dipoles induced during binary collisions of N2 molecules using independently measured molecular parameters; the results were in reasonable agreement with experimental data. However, recent measurements have revealed new fine structure that has been attributed to line-mixing effects. We do not think that this is correct, rather that the structure results from short-range anisotropic dipoles. We are in the process of including this refinement in our theoretical calculation in order to compare with the new experimental data. Subtle changes in the spectra of CH4 measured by researchers at Langley have also been attributed to line-mixing effects. By analyzing the same spectral lines we have attempted to verify or rule out possible line-mixing mechanisms. Due to the complexity and richness of the spectrum of this highly symmetric molecule, as well as the small magnitude of the effects, a detailed first-principle calculation of the mixing is a difficult problem. Before such a program is undertaken it is important to glean as much information as possible concerning the possible mechanisms by a systematic analysis of the existing data.

Tipping, Richard H.↗
DOE OSTI2022

Multiscale modeling of hydrogenolysis of ethane and propane on Ru(0001): Implications for plastics recycling

Plastic waste presents an environmental threat. Chemical recycling via hydrogenolysis can convert plastic waste into waxes, lubricants, and fuels. Among catalysts, Ru stands out for its superior activity and selectivity. The chemistry of light alkane hydrogenolysis can help understanding plastics deconstruction. Here, we perform first-principles calculations, develop descriptor-based relations, and conduct microkinetic modeling and analysis on ethane and propane. Predictions are in excellent agreement with experimental data. We identify a similar cracking pattern for both hydrocarbons entailing a deeply dehydrogenated species with the removal of four hydrogen atoms: CHCH*+* → 2CH* for ethane and CH 3 CCH*+2* → CH 3 C* + CH* for propane. We find that the rate-determining step is the C-C cracking for ethane and the first dehydrogenation from the terminal carbon (CH 3 CH 2 CH 3 *+* → CH 3 CH 2 CH 2 *+H*) for propane. The vinyl species CH 2 CH* produced from propane cracking is responsible for whether a single or multiple cracking events occur and effectively controls the selectivity. Specifically, ethane formation in propane hydrogenolysis is suppressed at elevated temperatures due to over-cracking via multiple (two here) cracking events being preferred over hydrogenation and desorption of ethane from the catalyst. Our workflow and models provide a baseline for future studies on heavier hydrocarbons. Insights into recent experimental studies of polyethylene over Ru-based catalysts are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Bound-free Spectra for Diatomic Molecules

It is now recognized that prediction of radiative heating of entering space craft requires explicit treatment of the radiation field from the infrared (IR) to the vacuum ultra violet (VUV). While at low temperatures and longer wavelengths, molecular radiation is well described by bound-bound transitions, in the short wavelength, high temperature regime, bound-free transitions can play an important role. In this work we describe first principles calculations we have carried out for bound-bound and bound-free transitions in N2, O2, C2, CO, CN, NO, and N2+. Compared to bound ]bound transitions, bound-free transitions have several particularities that make them different to deal with. These include more complicated line shapes and a dependence of emission intensity on both bound state diatomic and atomic concentrations. These will be discussed in detail below. The general procedure we used was the same for all species. The first step is to generate potential energy curves, transition moments, and coupling matrix elements by carrying out ab initio electronic structure calculations. These calculations are expensive, and thus approximations need to be made in order to make the calculations tractable. The only practical method we have to carry out these calculations is the internally contracted multi-reference configuration interaction (icMRCI) method as implemented in the program suite Molpro. This is a widely used method for these kinds of calculations, and is capable of generating very accurate results. With this method, we must first of choose which electrons to correlate, the one-electron basis to use, and then how to generate the molecular orbitals.

Schwenke, David W.↗
DOE OSTI2020

Linking electronic structure calculations to generalized stacking fault energies in multicomponent alloys

Abstract The generalized stacking fault energy is a key ingredient to mesoscale models of dislocations. Here we develop an approach to quantify the dependence of generalized stacking fault energies on the degree of chemical disorder in multicomponent alloys. We introduce the notion of a “configurationally-resolved planar fault” (CRPF) energy and extend the cluster expansion method from alloy theory to express the CRPF as a function of chemical occupation variables of sites surrounding the fault. We apply the approach to explore the composition and temperature dependence of the unstable stacking fault energy (USF) in binary Mo–Nb alloys. First-principles calculations are used to parameterize a formation energy and CRPF cluster expansion. Monte Carlo simulations show that the distribution of USF energies is significantly affected by chemical composition and temperature. The formalism is broadly applicable to arbitrary crystal structures and alloy chemistries and will enable the development of rigorous models for deformation mechanisms in high-entropy alloys.

Natarajan, Anirudh Raju (ORCID:0000000202529480)↗