Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

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↗
DOE OSTI2017

Initiation reactions in acetylene pyrolysis

In gas-phase combustion systems the interest in acetylene stems largely from its role in molecular weight growth processes. The consensus is that above 1500 K acetylene pyrolysis starts mainly with the homolytic fission of the C–H bond creating an ethynyl radical and an H atom. However, below ~1500 K this reaction is too slow to initiate the chain reaction. It has been hypothesized that instead of dissociation, self-reaction initiates this process. Nevertheless, rigorous theoretical or direct experimental evidence is lacking, to an extent that even the molecular mechanism is debated in the literature. In this work we use rigorous ab initio transition-state theory master equation methods to calculate pressure- and temperature-dependent rate coefficients for the association of two acetylene molecules and related reactions. We establish the role of vinylidene, the high-energy isomer of acetylene in this process, compare our results with available experimental data, and assess the competition between the first-order and second-order initiation steps. As a result, we also show the effect of the rapid isomerization among the participating wells and highlight the need for time-scale analysis when phenomenological rate coefficients are compared to observed time scales in certain experiments.

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

Resolving Heterogeneous Dynamics of Excess Protons in Aqueous Solution with Rate Theory

Rate theories have found great utility across the chemical sciences by providing a physically transparent way to analyze dynamical processes. Here we demonstrate the benefits of using transition state theory and Marcus theory to study the rate of proton transfer in HCl solutions. By using long ab initio molecular dynamics simulations we show that good agreement is obtained between these two different formulations of rate theory and how they can be used to study the pathways and life-time of proton transfer in aqueous solution. Since both rate theory formulations utilize identical sets of molecular data, this provides a self-consistent theoretical picture of the rates of aqueous phase proton transfer. Specifically, we isolate and quantify the rates of proton transfer, ion-pair dissociation, and solvent exchange in concentrated HCl solutions. Our analysis predicts a concentration dependence to both proton transfer and ion-pairing. Moreover, our estimate of the life-time for the Zundel species is 0.8 ps and 1.3 ps for 2M and 8M HCl, respectively. In this paper, we demonstrate that concentration effects, can indeed be quantified through the combination of state-of-the-art simulation and theory and provides a picture of the important correlations between the cation (hydronium) and the counter-ion in acid solutions.

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

PyCrystalField : software for calculation, analysis and fitting of crystal electric field Hamiltonians

PyCrystalField is a Python software package for calculating single-ion crystal electric field (CEF) Hamiltonians. This software can calculate a CEF Hamiltonian ab initio from a point charge model for any transition or rare earth ion in either the J basis or the LS basis, perform symmetry analysis to identify nonzero CEF parameters, calculate the energy spectrum and observables such as neutron spectrum and magnetization, and fit CEF Hamiltonians to any experimental data. Here, the theory, implementation and examples of its use are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Density functional theory (DFT) study of UF 6 hydrolysis: reaction pathways, spectroscopy, and chemical kinetics

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the gas-phase reaction mechanism and associated chemical kinetics are poorly understood. During the performance period we undertook development of a state-of-the-art ab initio gas-phase chemical kinetics simulation workflow to model the hydrolysis of uranium hexafluroride (UF 6 ). In doing so, we addressed several outstanding issues in the theoretical treatment of uranium-containing systems. At the outset it was unclear how to generate accurate estimates of kinetic and thermodynamic data for U-containing chemical reactions. Generation of such data has been made routine. Prior to our work, the literature associated with UF 6 hydrolysis were disparate and inaccurate. This body of work provides a modern and comprehensive theoretical assessment of the reaction mechanism, molecular clustering towards deposition, and chemical kinetics. New methodological implementations and software integrations resulting from this work are also highlighted. As much as possible, our predictions were validated against experimental data including particle morphologies, vibrational spectroscopy, atomization enthalpies, and kinetic rate constants. Nevertheless, we were unable to reconcile kinetic measurements with high-accuracy simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

High-Temperature Thermodynamics Modeling of Graphite

We present high-temperature thermodynamic properties for graphite from first-principles anharmonic theory. The ab initio electronic structure is obtained from density-functional theory coupled to a lattice dynamics method that is used to model anharmonic lattice vibrations. This combined approach produces free energies and specific heats for graphite that compare well with available experiments and results from models that empirically represent experimental data, such as CALPHAD. We show that anharmonic theory for the phonons is essential for accurate thermodynamic quantities above about 1000 K.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Influence of Cation Size on the Electronic and Magnetic Structure of M CoO{sub 3} (M = Bi, La) from Ab Initio Calculations

Using the first-principle DFT method and LSDA+ U calculations, we have performed a detailed investigation on the structural and magnetic ground states of strongly correlated M CoO{sub 3} (M= Bi, La) perovskites. Magneticphase stabilities, structural behaviors and electronic properties of both materials have been investigated. LSDA+ U calculations show that the antiferromagnetic (AFM)-C phase in the tetragonal La/BiCoO{sub 3} is energetically more favorable than the other configurations, and both rhombohedral materials are nonmagnetic. All structural data of tetragonal BiCoO{sub 3} and rhombohedral LaCoO{sub 3} compounds are found in excellent agreement with the experiment for U{sub eff}= 6 eV. We have predicted also band gaps and magnetic moments with a good accuracy via the band structure calculations. From the total and partial densities of states, we have shown the important role of O p Co d and Bi/La s-O p hybridizations in both systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ab initio calculation of the β-decay spectrum of 6 He

Here, we calculate the β spectrum in the decay of 6 He using quantum Monte Carlo methods with nuclear interactions derived from chiral effective field theory and consistent weak vector and axial currents. We work at second order in the multipole expansion, retaining terms suppressed by O(q 2 /m$^{2}_{π}$), where q denotes low-energy scales such as the reaction’s Q value or the electron energy, and m π is the pion mass. We go beyond the impulse approximation by including the effects of two-body vector and axial currents. We estimate the theoretical error on the spectrum by using four potential models in the Norfolk family of local two- and three-nucleon interactions, which have different cutoffs, fit two-nucleon data up to different energies, and use different observables to determine the couplings in the three-body force. We find the theoretical uncertainty on the β spectrum, normalized by the total rate, to be well below the permille level, and to receive contributions of comparable size from first- and second-order corrections in the multipole expansion. We consider corrections to the β decay spectrum induced by beyond-standard-model charged-current interactions in the standard model effective field theory, with and without sterile neutrinos, and discuss the sensitivity of the next generation of experiments to these interactions.

6 ≤ A ≤ 19↗
DOE OSTI2023

AL-ASMR: Active Learning of Atomistic Surrogate Models for Rare Events

Atomistic simulation with artificial intelligence (AI) is an emerging tool for understanding materials' properties and behaviors and predicting novel materials with optimized/targeted properties. Neural network potentials (NNPs) are outstanding in this field as they have shown a comparable accuracy to ab initio electronic structure calculations for reproducing potential energy surfaces while being several orders of magnitude faster. However, such NNPs can perform poorly outside of their training domain and typically fail catastrophically in the prediction of rare events in molecular dynamics (MD) simulations. For effective AL loops to distinguish the informative data from enhanced sampled configurations, we developed a decision engine by configurational similarity and uncertainty quantification (UQ) with data augmentation.

Jung, Gang Seob↗
NASA NTRS1986

Electron-nitrogen molecular collisions in high-temperature nonequilibrium air

Ab initio calculation of vibrational excitation and deexcitation cross sections of N2 by low energy resonant electron impact have been carried out. The calculation includes initial target state in v = 0-12 and final target state in v = 0-17. The calculated cross sections are found to agree with available experimental data to within reported experimental error. Effect of target rotation has been investigated. Comparison of vibrational excitation cross sections at J = 0, 50, and 150 shows that the major effect of target rotation is to lower the resonance energy. Vibrational excitation and deexcitation rate coefficients are tabulated for electron temperatures of 1100-55,000 K. The present tabulation supercedes an earlier compilation by the Naval Research Laboratory that was derived from experimental data with a factor-of-two error. The tabulated rate coefficients are suited for flowfield modeling of aeroassisted orbital transfer vehicles.

Huo, W. M.↗
DOE OSTI2020

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

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

Large exchange-driven intrinsic circular dichroism of a chiral 2D hybrid perovskite

Abstract In two-dimensional chiral metal-halide perovskites, chiral organic spacers endow structural and optical chirality to the metal-halide sublattice, enabling exquisite control of light, charge, and electron spin. The chiroptical properties of metal-halide perovskites have been measured by transmissive circular dichroism spectroscopy, which necessitates thin-film samples. Here, by developing a reflection-based approach, we characterize the intrinsic, circular polarization-dependent complex refractive index for a prototypical two-dimensional chiral lead-bromide perovskite and report large circular dichroism for single crystals. Comparison with ab initio theory reveals the large circular dichroism arises from the inorganic sublattice rather than the chiral ligand and is an excitonic phenomenon driven by electron-hole exchange interactions, which breaks the degeneracy of transitions between Rashba-Dresselhaus-split bands, resulting in a Cotton effect. Our study suggests that previous data for spin-coated films largely underestimate the optical chirality and provides quantitative insights into the intrinsic optical properties of chiral perovskites for chiroptical and spintronic applications.

25 ENERGY STORAGE↗
DOE OSTI2014

An ab initio study of the electronic structure of the boron oxide neutral (BO), cationic (BO{sup +}), and anionic (BO{sup −}) species

The BO neutral, cationic, and anionic molecular species have been painstakingly studied through multireference configuration interaction and single reference coupled cluster methods employing basis sets of quintuple cardinality. Potential energy curves have been constructed for 38 (BO), 37 (BO{sup +}), and 12 (BO{sup −}) states and the usual molecular parameters have been extracted most of which are in very good agreement with the scarce experimental data. Numerous avoided crossings appear on more or less all of the studied states of the neutral and cationic species challenging the validity of the Born Oppenheimer approximation. Finally, all excited states of the anionic system lie above the ground state of the neutral BO system and are therefore resonances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
NASA NTRS2011

Assessing the Age of an Asteroid's Surface with Data from the International Rosetta Mission

Rosetta is an international mission led by the European Space Agency (ESA) with key support and instrumentation from the National Aeronautics and Space Administration (NASA). Rosetta is currently on a ten-year mission to catch comet 67P/Churyumov-Gerasimenko (C-G); throughout its voyage, the spacecraft has performed flybys of two main belt asteroids (MBA): Steins and Lutetia. Data on the physical, chemical, and geological properties of these asteroids are currently being processed and analyzed. Accurate interpretation of such data is fundamental in the success of Rosetta's mission and overall objectives. Post-flyby data analyses strive to correlate the size, shape, volume, and rotational rate of Lutetia, in addition to interpreting its multi-color imagining, albedo, and spectral mapping. Although advancements in science have contributed to the examination of celestial bodies, methods to analyze asteroids remain largely empirical, not semi-empirical, nor ab initio. This study aims to interpret and document the scientific methods currently utilized in the characterization of asteroid (21) Lutetia in order to render these processes and methods accessible to the public. Examples include a standardized technique for assessing the age of an asteroid surface, complete with clickable reference maps, methodology of grouping surface characteristics together, and a standardized power law equation for the age. Other examples include determining the density of an object. Context for what both density and age mean is a bi-product of this study. Results of the study will aid in the development of pedagogical material on asteroids for public use, and in creation of an academic database for selected targets that might be used as a reference.

asteriods↗
DOE OSTI2019

Shock compression of niobium from first-principles

The equation of state (EOS) of bulk niobium (Nb) was investigated within the framework of density functional theory, with Mermin’s generalization to finite temperatures. The shock Hugoniot for fully-dense and porous Nb was obtained from canonical ab initio molecular dynamics simulations with Erpenbeck’s approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled along isotherms between 300 and 4000 K, for densities ranging from ρ = 5.5 to 12 g/cm 3 . Results from simulations compare favorably with room-temperature multianvil and diamond anvil cell data for fully-dense Nb samples and with a recent tabulated SESAME EOS. Here, the results of this study indicate that, for the application of weak and intermediate shocks, the tabular EOS models are expected to give reliable predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0−1 and only five systems (0.07%) with an absolute error larger than 0.01 a0−1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Chemistry↗
DOE OSTI2021

Database of ab initio L-edge X-ray absorption near edge structure

Abstract The L-edge X-ray Absorption Near Edge Structure (XANES) is widely used in the characterization of transition metal compounds. Here, we report the development of a database of computed L-edge XANES using the multiple scattering theory-based FEFF9 code. The initial release of the database contains more than 140,000 L-edge spectra for more than 22,000 structures generated using a high-throughput computational workflow. The data is disseminated through the Materials Project and addresses a critical need for L-edge XANES spectra among the research community.

42 ENGINEERING↗