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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An automatically curated first-principles database of ferroelectrics

Ferroelectric materials have technological applications in information storage and electronic devices. The ferroelectric polar phase can be controlled with external fields, chemical substitution and size-effects in bulk and ultrathin film form, providing a platform for future technologies and for exploratory research. In this work, we integrate spin-polarized density functional theory (DFT) calculations, crystal structure databases, symmetry tools, workflow software, and a custom analysis toolkit to build a library of known, previously-proposed, and newly-proposed ferroelectric materials. With our automated workflow, we screen over 67,000 candidate materials from the Materials Project database to generate a dataset of 255 ferroelectric candidates, and propose 126 new ferroelectric materials. We benchmark our results against experimental data and previous first-principles results. The data provided includes atomic structures, output files, and DFT values of band gaps, energies, and the spontaneous polarization for each ferroelectric candidate. Finally, we contribute our workflow and analysis code to the open-source python packages atomate and pymatgen so others can conduct analogous symmetry driven searches for ferroelectrics and related phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Voltage-Dependent First-Principles Simulation of Insertion of Chloride Ions into Al/Al 2 O 3 Interfaces Using the Quantum Continuum Approximation

Experiments have shown that pitting corrosion can develop in aluminum surfaces at potentials > − 0.5 V relative to the standard hydrogen electrode (SHE). Until recently, the onset of pitting corrosion in aluminum has not been rigorously explored at an atomistic scale because of the difficulty of incorporating a voltage into density functional theory (DFT) calculations. We introduce the Quantum Continuum Approximation (QCA) which self-consistently couples explicit DFT calculations of the metal-insulator and insulator-solution interfaces to continuum Poisson-Boltzmann electrostatic distributions describing the bulk of the insulating region. By decreasing the number of atoms necessary to explicitly simulate with DFT by an order of magnitude, QCA makes the first-principles prediction of the voltage of realistic electrochemical interfaces feasible. After developing this technique, we apply QCA to predict the formation energy of chloride atoms inserting into oxygen vacancies in Al(111)/α-Al 2 O 3 (0001) interfaces as a function of applied voltage. We predict that chloride insertion is only favorable in systems with a grain boundary in the Al 2 O 3 for voltages > − 0.2 V (SHE). Here our results roughly agree with the experimentally demonstrated onset of corrosion, demonstrating QCA's utility in modeling realistic electrochemical systems at reasonable computational cost.

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

First-principles and cluster expansion study of the effect of magnetism on short-range order in Fe–Ni–Cr austenitic stainless steels

Short-range order (SRO), the regular and predictable arrangement of atoms over short distances, alters the mechanical properties of technologically relevant structural materials such as medium/high entropy alloys and austenitic stainless steels. In this study, we present a generalized spin cluster expansion (CE) model and show that magnetism is a primary factor influencing the level of SRO present in austenitic Fe-Ni-Cr alloys. The spin CE consists of a chemical cluster expansion combined with an Ising model for Fe-Ni-Cr austenitic alloys. It explicitly accounts for local magnetic exchange interactions, thereby capturing the effects of finite temperature magnetism on SRO. Model parameters are obtained by fitting to a first-principles data set comprising both chemically and magnetically diverse FCC configurations. The magnitude of the magnetic exchange interactions are found to be comparable to the chemical interactions. Compared to a conventional implicit magnetism CE built from only magnetic ground state configurations, the spin CE shows improved performance on several experimental benchmarks over a broad spectrum of compositions, particularly at higher temperatures due to the explicit treatment of magnetic disorder. We find that SRO is strongly influenced by alloy Cr content, since Cr atoms prefer to align antiferromagnetically with nearest neighbors but become magnetically frustrated with increasing Cr concentration. Using the spin CE, we predict that increasing the Cr concentration in typical austenitic stainless steels promotes the formation of SRO and increases order-disorder transition temperatures. Furthermore, this study underscores the significance of considering magnetic interactions explicitly when exploring the thermodynamic properties of complex transition metal alloys. It also highlights guidelines for customizing SRO through adjustments of alloy composition.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A First-principle Model for Polarization Swings during Reconnection-powered Flares

We show that magnetic reconnection in a magnetically dominated fast-cooling plasma can naturally produce bright flares accompanied by rotations in the synchrotron polarization vector. With particle-in-cell simulations of reconnection, we find that flares are powered by efficient particle acceleration at the interface of merging magnetic flux ropes, or “plasmoids.” The accelerated particles stream through the post-merger plasmoid toward the observer, thus progressively illuminating regions with varying plane-of-sky field direction, and so leading to a rotation in the observed polarization vector. Our results provide evidence for magnetic reconnection as the physical cause of highenergy flares from the relativistic jets of blazars (which recent observations have shown to be frequently associated with polarization rotations), and provide a first-principle physical mechanism for such flares.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

High-accuracy first-principles-based rate coefficients for the reaction of OH and CH 3 OOH

Here, the ˙OH-initiated oxidation of methyl hydroperoxide, which plays an important role in the atmospheric chemistry of methane, was theoretically characterized using high-accuracy composite amHEAT-345(Q) coupled-cluster calculations followed by a two-dimensional E,J resolved master equation analysis. The reaction is found to proceed through two distinct hydrogen-bonded pre-reactive complexes leading to two product channels, in accord with the experimental observations: (i) ˙OH + CH 3 OOH → CH 3 OO˙ + H 2 O with a yield of 0.8 ± 0.1, and (ii) ˙OH + CH 3 OOH → HCHO + ˙OH + H 2 O with a yield of 0.2 ± 0.1. The calculated reaction enthalpies are within 0.2 kcal mol -1 of the benchmark ATcT values. Overall thermal rate coefficients obtained from first principles are found to be in the low-pressure limit at atmospheric pressure; the total rate coefficient can be expressed over the T = 200–450 K range as k(T) = 5.0 × 10 -12 × T -0.152 × exp(287/T) cm 3 s -1 , strongly supporting the experimental results of Vaghjiani and Ravishankara, with which this expression agrees within ca. 15%. The current results show that (i) is the principal reaction channel and support the view that, due to its inherently fast transformations, CH 3 OOH is an important redistribution species for HO x ˙ radicals in the Earth's atmosphere.

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

Combined First-Principles and Experimental Investigation into the Reactivity of Codeposited Chromium–Carbon under Pressure

High-pressure synthesis in the diamond anvil cell suffers from the lack of a general approach for the control of precursor stoichiometry and homogeneity. Here, we present results from a new method we have developed that uses magnetron cosputtering to prepare stoichiometrically precise and atomically mixed amorphous films of Cr:C. Laser-heated diamond anvil cell experiments carried out on a flake of this sample at pressures between 13.5 and 24.3 GPa lead to the observation of Cr 3 C (Pnma) over the entire pressure range–in good agreement with our in-house theoretical predictions–but also reveal two other metastable phases that were not expected: a novel monoclinic chromium carbide phase and the NaCl-type CrC (Fm3̅m) phase. The unexpected stability of CrC is investigated by using first-principles methods, revealing a large stabilizing effect tied to substoichiometry at the carbon site. These results offer an important case study into the current limitations of crystal structure prediction methods with regard to phase complexity and bolster the growing need for advanced theoretical approaches that can more completely survey experimentally unexplored phase space.

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

Structure, Magnetism, and First-Principles Modeling of the Na 0.5 La 0.5 RuO 3 Perovskite

High-purity polycrystalline Na 0.5 La 0.5 RuO 3 was synthesized by a solid-state method, and its properties were studied by magnetic susceptibility, heat capacity, and resistivity measurements. In this study, we find it to be an orthorhombic perovskite, in contrast to an earlier report, with random La/Na mixing. With a Curie–Weiss temperature of –231 K and an effective moment of 2.74 μB/mol Ru, there is no magnetic ordering down to 1.8 K. A broad hump at 1.4 K in the heat capacity, however, indicates the presence of a glassy magnetic transition, which we attribute to the influence of the random distribution of Na and La on the perovskite A-sites. Comparison to CaRuO 3 , a structurally ordered ruthenate perovskite with a similar A-site ionic radius and magnetic properties, is presented. First-principles calculations indicate that the Na–La distribution determines the local magnetic exchange interactions between Ru ions, favoring either antiferromagnetic or ferromagnetic coupling when the local environment is Na- or La-rich. Thus, our data and analysis suggest that mixing cations with different charges and sizes on the A-site in this perovskite results in magnetic frustration through a balance of local magnetic exchange interactions.

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

A first principles framework to predict the transient performance of latent heat thermal energy storage

Thermal energy storage (TES) is increasingly recognized as an essential component of efficient Combined Heat and Power (CHP), Concentrated Solar Power (CSP), Heating Ventilation and Air Conditioning (HVAC), and refrigeration as it reduces peak demand while helping to manage intermittent availability of energy (e.g., from solar or wind). Latent Heat Thermal Energy Storage (LHTES) is a viable option because of its high energy storage density. Parametric analysis of LHTES in terms of dimensionless numbers is highly desired as a tool to model LHTES systems. One approach is to develop a model equation so as to minimize the error between the model and data obtained from experiments or simulations. While this approach can produce an accurate correlation applicable within the range of data used for its creation, it does not provide physical understanding of the rate-limiting process controlling the transient behavior of the device. In this paper we present an alternative approach whereby the potential rate-limiting processes are identified from first principles and then the key process is determined as a function of time as a LHTES device is charged. For example, in a simple geometry, the melt-fraction can be expected to vary linearly in time if the heat transfer rate is limited by natural convection of the phase changing material and we show it scales with the PCM Grashof number as $Gr^1_p$ and PCM Prandtl number as $Pr_p^{(1/3)}$. On the other hand, if surface area of solid PCM limits the heat transfer rate, the melt fraction increases asymptotically to reach full melting. The existence of these linear and asymptotic regions and the $Gr^1_pP r^{1/3}_p$ shape of the melt fraction curve is verified using our database of 64 simulations. Of practical importance in designing LHTES devices is the melt fraction at which the heat transfer rate ceases to be limited by convection, after which the heat storage rate deteriorates. For our geometry, this is found to be about 90%. This test case of our methodology shows the value of our approach, that predicting heat storage rate based on the rate-limiting physical phenomenon as a function of time is an effective approach to modeling LHTES devices.

25 ENERGY STORAGE↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

First-principles calculation of lattice distortion, electronic structure, and bonding properties of GeTe-based and PbSe-based high-entropy chalcogenides

The massive amount of wasted heat energy from industry has pushed the development of thermoelectric (TE) materials that directly convert heat into electricity to a new level of concern. Recently, multicomponent alloys such as GeTe-based and PbSe-based high-entropy (HE) chalcogenides have attracted a great deal of attention due to their potential application as TE materials. The nature of the interatomic bonding, lattice distortion (LD), and the electronic structure in this class of materials is not fully understood. Herein, we report a comprehensive computational investigation of nine GeTe-based HE alloys with eight metallic elements (Ag, Pb, Sb, Bi, Cu, Cd, Mn, and Sn) with large supercells of 1080 atoms each; seven PbSe-based HE solid solutions: Pb0.99−ySb0.012SnySe1−2xTexSx (x = 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, with y = 0) with supercells of 1000 atoms each; and five Pb0.99−ySb0.012SnySe1−2xTexSx (y = 0.05, 0.1, 0.15, 0.2, 0.25 with x = 0.25) solid solutions with supercells of 1000 atoms each. All these HE models are theoretically investigated for the first time. The electronic structure, interatomic bonding, charge transfer, and lattice distortion (LD) are investigated by first-principles calculations based on density functional theory. Multicomponent HE alloys can cause a significant LD, which affects their mechanical, thermal, and TE properties. The calculations for the GeTe-based HE chalcogenides showed that they are semiconductors with a narrow bandgap, except for m8, which has a semi-metallic characteristic, and this makes them good candidates for TE applications. For most of these models, the Fermi level shifts upward and locates deeply in the conduction bands, resulting in the enhancement of the electrical conductivity (σ). The bonding properties showed that most bonds in m5 are more dispersed, indicating highest LD and lower lattice thermal conductivity. For PbSe-based HE solid solutions, the LD calculations showed that the models Pb0.99Sb0.012Se0.5Te0.25S0.25 and Pb0.89Sb0.012Sn0.1Se0.5Te0.25S0.25 have the higher LD, and thus a lower lattice thermal conductivity. Such investigations are in high demand since it enables us to design new HE chalcogenides for TE applications. We use the novel concept of total bond order density as a single quantum mechanical metric to characterize the internal cohesion of these HE alloys and correlate with calculated properties, especially the mechanical properties. This work provides a solid database for HE chalcogenides and a road map for many potential applications. Moreover, the computational procedure we developed can be used to design new HE chalcogenides for specific TE applications.

Physics↗
DOE OSTI2022

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

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

Quaternary i-MAX Phases (Mo 2/3 RE 1/3 ) 2 AlC (RE: Dy, Tb, Er): Experimental Characterization and First-Principles Insights into their Fundamental Properties

Rare earth (RE)-based materials have unique electronic, magnetic, and optical properties, leading to the recent discovery of atomically layered solids with the chemical formula (M' 2/3 RE 1/3 ) 2 AlC, which have since garnered significant attention in the scientific community. This study aims to synthesize, characterize, and investigate the structural and thermal stability of the RE i-MAX phases. We prepared i-MAX phases using molybdenum (Mo) as M′ and RE elements as Dy, Tb, and Er, namely (Mo 2/3 Dy 1/3 ) 2 AlC, (Mo 2/3 Tb 1/3 ) 2 AlC, and (Mo 2/3 Er 1/3 ) 2 AlC. Structural characterization through x-ray diffraction (XRD) and Raman spectroscopy confirms the formation of the RE-based i-MAX phase, along with the presence of minor impurity phases in the alloys. Thermogravimetric analysis (TGA) conducted up to 1000°C under ambient conditions reveals that the i-MAX phases remain thermally stable up to approximately 450°C, beyond which oxidation leads to a noticeable weight gain in all samples. Differential scanning calorimetry (DSC) measurements during heating and cooling cycles show endothermic and exothermic peaks for (Mo 2/3 Dy 1/3 ) 2 AlC i-MAX in the 410–420°C range, indicating a temperature-induced minor atomic arrangement. In contrast, these peaks are absent in the Tb- and Er-based i-MAX phases. These findings offer valuable insights into the thermal behavior and stability of these i-MAX phases under thermal stress, contributing to a deeper understanding of their unique properties. Furthermore, first-principles density functional theory (DFT) calculations were performed to investigate the electronic and optical properties of the i-MAX phases. The results reveal their metallic nature, with pronounced contributions from Mo and RE elements near the Fermi level and within the conduction band.

Rare earth↗
DOE OSTI2023

Optical activity of solids from first principles

Within the framework of independent particle approximation, the optical activity tensor of solids is formulated as from different contributions: the magnetic dipole, electric quadrupole, and band dispersion terms. The first two terms have similar counterparts in the theory of finite systems, whereas the last term is unique for crystals. Further, the magnetic dipole and electric quadrupole transition moments are calculated with a sum-over-states formulation. We apply the formulation to calculate and analyze the optical rotation of elemental tellurium and the circular dichroism of the (6,4) carbon nanotube. Decomposed optical activity into different contributions is discussed. The calculated spectra agree well with experiments. As a showcase of achiral crystals, we calculate the optical activity of wurtzite GaN.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

A first principles approach for determining solute effect on electrical resistivity of Aluminum solid solution

There are no general approaches for predicting the effect of solutes on electrical resistivity changes in aluminum (Al) solid solutions. Here, we introduce a first principles-based approach for predicting temperature-dependent electrical resistivity of Al solid solutions. A combination of electronic structure methods was applied to calculate electronic scattering effects on chemical disorder and temperature induced atomic vibrations or phonons. Our calculations are validated by contrasting experimental resistivity measurements for two elements in Al solid solution: at 0.3 at. % concentration, Sn (an sp-valence solute) increases the residual resistivity by 0.2 μΩ·m, whereas Zr (a d-valence solute) causes a 2.5 μΩ·m increase. This trend may deviate for solute elements in Al that form semiconductors in pure state, like Si. Our findings provide an understanding of the role of valence electron character of solute atoms in determining electrical resistivity of Al alloys.

Samolyuk, German [ORNL] (ORCID:0000000168778255)↗