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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling Intercalation Chemistry with Multiredox Reactions by Sparse Lattice Models in Disordered Rocksalt Cathodes

Modern battery materials can contain many elements with substantial site disorder, and their configurational state has been shown to be critical for their performance. The intercalation voltage profile is a critical parameter to evaluate the performance of energy storage. The application of commonly used cluster expansion techniques to model the intercalation thermodynamics of such systems ab initio is challenged by the combinatorial increase in configurational degrees of freedom as the number of species grows. Such challenges necessitate the efficient generation of lattice models without overfitting and proper sampling of the configurational space under the requirement of charge balance in ionic systems. In this work, we introduce a combined approach that addresses these challenges by (1) constructing a robust cluster expansion Hamiltonian using the sparse regression technique, including -norm regularization and structural hierarchy; and (2) implementing semigrand-canonical Monte Carlo to sample charge-balanced ionic configurations using the table-exchange method and an ensemble average approach. These techniques are applied to a disordered rocksalt oxyfluoride (LMNOF) that is part of a family of promising earth-abundant cathode materials. The simulated voltage profile is found to be in good agreement with experimental data and particularly provides a clear demonstration of the and oxygen contributions to the redox potential as a function of content.

25 ENERGY STORAGE↗
DOE OSTI2017

First-principles study on the electronic, optical and thermodynamic properties of ABO 3 (A = La,Sr, B = Fe,Co) perovskites

The electronic, optical and thermodynamic properties of ABO 3 (A = La,Sr, B = Fe,Co) perovskites are investigated using first-principles calculations. The obtained results indicate that SrCoO 3 and SrFeO 3 are metals, while LaCoO 3 and LaFeO 3 are insulators and all of them exhibit strong hybridization of the Fe/Co-3d and O-2p orbitals. By correlating the energy band structures with the peaks of the imaginary part of the dielectric function, we obtained the origin of each electron excitation to provide information about the active bands for the corresponding optical transitions observed in the experiment. Moreover, the Debye temperatures θ D obtained from the phonon frequencies are comparable to the available data. In conclusion, the thermodynamic properties of the Helmholtz free energy F , entropy S , and constant-volume heat capacity C v are investigated based on the phonon spectra.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Light nuclei with semilocal momentum-space regularized chiral interactions up to third order

In this work, we present a systematic investigation of few-nucleon systems and light nuclei using the current Low Energy Nuclear Physics International Collaboration interactions comprising semilocal momentum-space regularized two- and three-nucleon forces up to third chiral order (N 2 LO). Following our earlier study utilizing the coordinate-space regularized interactions, the two low-energy constants entering the three-body force are determined from the triton binding energy and the differential cross-section minimum in elastic nucleon-deuteron scattering. Predictions are made for selected observables in elastic nucleon-deuteron scattering and in the deuteron breakup reactions, for properties of the A = 3 and A = 4 nuclei, and for spectra of p-shell nuclei up to A = 16. A comprehensive error analysis is performed including an estimation of correlated truncation uncertainties for nuclear spectra. The obtained predictions are generally found to agree with experimental data within errors. Similarly to the coordinate-space regularized chiral interactions at the same order, a systematic overbinding of heavier nuclei is observed, which sets in for A~10 and increases with A.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

First-principles Hubbard U approach for small molecule binding in metal-organic frameworks

We apply first-principles approaches with Hubbard U corrections for calculation of small molecule binding energetics to open-shell transition metal atoms in metal-organic frameworks (MOFs). Using density functional theory with van der Waals dispersion-corrected functionals, we determine Hubbard U values ab initio through an established linear response procedure for M-MOF-74, for a number of different metal centers (M = Ti, V, Cr, Mn, Fe, Co, Ni, and Cu). While our ab initio U values differ from those used in previous work, we show that they result in lattice parameters and electronic contributions to CO{sub 2}-MOF binding energies that lead to excellent agreement with experiments and previous results, yielding lattice parameters within 3%. In addition, U-dependent calculations for an example system, Co-MOF-74, suggest that the CO{sub 2} binding energy grows monotonically with the value of Hubbard U, with the binding energy shifting 4 kJ/mol (or 0.041 eV) over the range of U = 0-5.4 eV. These results provide insight into an approximate but computationally efficient means for calculation of small molecule binding energies to open-shell transition metal atoms in MOFs and suggest that the approach can be predictive with good accuracy, independent of the cations used and the availability of experimental data.

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

Revisiting Selenium Interactions with Pyrite: From Adsorption to Coprecipitation

Interactions of selenium (Se), a trace element bio-essential at low concentrations but highly toxic at high concentrations, with the most abundant sulfide mineral in the Earth's crust, namely pyrite, was investigated over a wide range of time scales. At the nanosecond scale, selenate Se(VI)O 4 2– adsorption onto the net pyrite surface is shown by ab-initio computations to proceed via the formation of a chemical bond between an oxyanion oxygen atom and a surface Fe atom, weakening the other Se-O bonds and reducing Se atom oxidation state. At the hour-to-day scale, adsorption and coprecipitation of selenate Se(VI)O 4 2– and selenite, Se(IV)O 3 2– , were investigated through wet chemical batch experiments at various pH values at different sulfide concentrations. Selenium removal from solution is slower and weaker for selenate than for selenite. After 24h, only 10% of selenate, against 60% of selenite (up to 100% in the presence of sulfide), is removed by the pyrite surface. Independently of its original oxidation state, adsorbed Se is completely reduced to elemental trigonal selenium via adsorption, precipitation or coprecipitation, as shown by XANES spectroscopy. Here, our EXAFS results, compared to published data on Se-rich pyrite, show a Se to S substitution within the pyrite structure. The reductive coprecipitation mechanism of selenium with pyrite represents valuable new insights for improving our understanding of modern and ancient biogeochemical cycles involving Se. In addition, several industries can benefit from direct applications of our findings, such as water treatment, green technologies and sustainable mining.

Adsorption↗
DOE OSTI2014

Depolarizing collisions with hydrogen: Neutral and singly ionized alkaline earths

Depolarizing collisions are elastic or quasielastic collisions that equalize the populations and destroy the coherence between the magnetic sublevels of atomic levels. In astrophysical plasmas, the main depolarizing collider is neutral hydrogen. We consider depolarizing rates on the lowest levels of neutral and singly ionized alkali earths Mg I, Sr I, Ba I, Mg II, Ca II, and Ba II, due to collisions with H°. We compute ab initio potential curves of the atom-H° system and solve the quantum mechanical dynamics. From the scattering amplitudes, we calculate the depolarizing rates for Maxwellian distributions of colliders at temperatures T ≤ 10,000 K. A comparative analysis of our results and previous calculations in the literature is completed. We discuss the effect of these rates on the formation of scattering polarization patterns of resonant lines of alkali earths in the solar atmosphere, and their effect on Hanle effect diagnostics of solar magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbation estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10, 472 +/- 5, and 2051 +/- 10 cm(exp -1) for omega(sub 1), omega(sub 2), and omega(sub 3), respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental nu(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗
DOE OSTI2016

Thermophysical properties of liquid Ni around the melting temperature from molecular dynamics simulation

Thermophysical properties of liquid nickel (Ni) around the melting temperature are investigated by means of classical molecular dynamics (MD) simulation, using three different embedded atom method potentials to model the interactions between the Ni atoms. Melting temperature, enthalpy, static structure factor, self-diffusion coefficient, shear viscosity, and thermal diffusivity are compared to recent experimental results. Using ab initio MD simulation, we also determine the static structure factor and the mean-squared displacement at the experimental melting point. For most of the properties, excellent agreement is found between experiment and simulation, provided the comparison relative to the corresponding melting temperature. We discuss the validity of the Hansen-Verlet criterion for the static structure factor as well as the Stokes-Einstein relation between self-diffusion coefficient and shear viscosity. The thermal diffusivity is extracted from the autocorrelation function of a wavenumber-dependent temperature fluctuation variable.

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

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2500+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with an emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with experimental data obtained from X-ray diffraction, infrared spectroscopy and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Mobility of singly-charged lanthanide cations in rare gases: Theoretical assessment of the state specificity

High quality, ab initio calculations are reported for the potential energy curves governing the interactions of four singly-charged lanthanide ions (Yb{sup +}, Eu{sup +}, Lu{sup +}, and Gd{sup +}) with the rare gases (RG = He–Xe). Scalar-relativistic coupled cluster calculations are used for the first three S-state ions, but for Gd{sup +}({sup 10}D°) it is necessary to take the interaction anisotropy into account with the help of the multi-reference technique. The potential energy curves are used to determine the ion mobility and other transport properties describing the motion of the ions through the dilute RG, both as functions of the temperature, T, in the low-field limit, and at fixed T as functions of the ratio of the electrostatic field strength to the gas number density, E/N. The calculated mobilities are in good agreement with the very limited experimental data that have become available recently. The calculations show a pronounced dependence of the transport properties on the electronic configuration of the ion, as well as a significant effect of the spin-orbit coupling on the transport properties of the Gd{sup +} ion, and predict that state-specific mobilities could be detectable in Gd{sup +}–RG experiments.

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

Thermoelectric properties of IV–VI-based heterostructures and superlattices

Doping in a manner that introduces anisotropy in order to reduce thermal conductivity is a significant focus in thermoelectric research today. By solving the semiclassical Boltzmann transport equations in the constant scattering time (τ) approximation, in conjunction with ab initio electronic structure calculations, within Density Functional Theory, we compare the Seebeck coefficient (S) and figure of merit (ZT) of bulk PbTe to PbTe/SnTe/PbTe heterostructures and PbTe doping superlattices (SLs) with periodically doped planes. Bismuth and Thallium were used as the n- and p-type impurities, respectively. The effects of carrier concentration are considered via chemical potential variation in a rigid band approximation. The impurity bands near the Fermi level in the electronic structure of PbTe SLs are of Tl s- and Bi p-character, and this feature is independent of the doping concentration or the distance between impurity planes. We observe the impurity bands to have a metallic nature in the directions perpendicular to the doping planes, yet no improvement on the values of ZT is found when compared to bulk PbTe. For the PbTe/SnTe/PbTe heterostructures, the calculated S presents good agreement with recent experimental data, and an anisotropic behavior is observed for low carrier concentrations (n<10{sup 18} cm{sup −3}). A large value of ZT{sub ||} (parallel to the growth direction) of 3.0 is predicted for n=4.7×10{sup 18} cm{sup −3} and T=700 K, whereas ZT{sub p} (perpendicular to the growth direction) is found to peak at 1.5 for n=1.7×10{sup 17} cm{sup −3}. Both electrical conductivity enhancement and thermal conductivity reduction are analyzed. - Graphical abstract: Figure of merit for PbTe/SnTe/PbTe heterostructure along the [0 0 1] direction, P.D. Borges, J.E. Petersen, L. Scolfaro, H.W. Leite Alves, T.H. Myers, Improved thermoelectric properties of IV–VI-based heterostructures and superlattices. - Highlights: • Thermoelectric properties of IV–VI-based heterostructures and superlattices. • High figure of merit is predicted for the PbTe/SnTe/PbTe heterostructure. • Nanotechnology has an important role for the development of thermoelectric devices.

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

Nonlinear ionization dynamics of hot dense plasma observed in a laser-plasma amplifier

Understanding the behaviour of matter under conditions of extreme temperature, pressure, density and electromagnetic fields has profound effects on our understanding of cosmologic objects and the formation of the universe. Lacking direct access to such objects, our interpretation of observed data mainly relies on theoretical models. However, such models, which need to encompass nuclear physics, atomic physics and plasma physics over a huge dynamic range in the dimensions of energy and time, can only provide reliable information if we can benchmark them to experiments under well-defined laboratory conditions. Due to the plethora of effects occurring in this kind of highly excited matter, characterizing isolated dynamics or obtaining direct insight remains challenging. High-density plasmas are turbulent and opaque for radiation below the plasma frequency and allow only near-surface insight into ionization processes with visible wavelengths. Here, the output of a high-harmonic seeded laser-plasma amplifier using eight-fold ionized krypton as the gain medium operating at a 32.8 nm wavelength is ptychographically imaged. A complex-valued wavefront is observed in the extreme ultraviolet (XUV) beam with high resolution. Ab initio spatio-temporal Maxwell–Bloch simulations show excellent agreement with the experimental observations, revealing overionization of krypton in the plasma channel due to nonlinear laser-plasma interactions, successfully validating this four-dimensional multiscale model. This constitutes the first experimental observation of the laser ion abundance reshaping a laser-plasma amplifier. The presented approach shows the possibility of directly modelling light-plasma interactions in extreme conditions, such as those present during the early times of the universe, with direct experimental verification.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Contribution of the entropy on the thermodynamic equilibrium of vacancies in nickel

The equilibrium vacancy concentration in nickel was determined from ab initio calculations performed with both generalized gradient approximation and local density approximation up to the melting point. We focus the study on the vacancy formation entropy expressed as a sum of a vibration and an electronic contribution, which were determined from the vibration modes and the electronic densities of states. Applying a method based on the quasi-harmonic approximation, the temperature dependence of the defect formation energy and entropy were calculated. We show that the vibrations of the first shell of atoms around the defect are predominant to the vibration formation entropy. On the other hand, the electronic formation entropy is very sensitive to the exchange-correlation potential used for the calculations. Finally, the vacancy concentration is computed at finite temperature with the calculated values for the defect formation energy and entropy. In order to reconcile point-defects concentration obtained with our calculations and experimental data, we conducted complementary calorimetric measurements of the vacancy concentration in the 1073–1273 K temperature range. Close agreement between theory and experiments at high temperature is achieved if the calculations are performed with the generalized gradient approximation and both vibration and electronic contributions to the formation entropy are taken into account.

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

Computational Spectroscopy of Polycyclic Aromatic Hydrocarbons In Support of Laboratory Astrophysics

Polycyclic aromatic hydrocarbons (PAHs) are strong candidates for the molecular carriers of the unidentified infrared bands (UIR) and the diffuse interstellar bands (DIBs). In order to test the PAH hypothesis, we have systematically measured the vibronic spectra of a number of jet-cooled neutral and ionized PAHs in the near ultraviolet (UV) to visible spectral ranges using the cavity ring-down spectroscopy. To support this experimental effort, we have carried out theoretical studies of the spectra obtained in our measurements. Ab initio and (time-dependent) density.functiona1 theory calculations are performed to obtain the geometries, energetics, vibrational frequencies, transition dipole moments, and normal coordinates of these PAH molecules. Franck-Condon (FC) calculations and/or vibronic calculations are then performed using the calculated normal coordinates and vibrational frequencies to simulate the vibronic spectra. It is found that vibronic interactions in these conjugated pi electron systems are often strong enough to cause significant deviations from the Born-Oppenheimer (BO) approximation. For vibronic transitions that are well described by the BO approximation, the vibronic band profiles are simulated by calculating the rotational structure of the vibronic transitions. Vibronic oscillator strength factors are calculated in the frame of the FC approximation from the electronic transition dipole moments and the FC factors. This computational effort together with our experimental measurements provides, for the first time, powerful tools for comparison with space-based data and, hence, a powerful approach to understand the spectroscopy of interstellar PAH analogs and the nature of the UIR and DIBs.

Tan, Xiaofeng↗
DOE OSTI2018

An experimental, theoretical, and modeling study of the ignition behavior of cyclopentanone

The ignition delay times of cyclopentanone in air were measured using a high pressure shock tube (HPST) and a rapid compression machine (RCM) over the temperature range of 794–1368 K at P = 15 and 30 bar and at equivalence ratios of 0.5, 1.0 and 2.0. To provide more insight into the oxidation of cyclopentanone, CO time-histories during cyclopentanone oxidation in a shock tube at high temperatures and various pressures were also measured. In addition, quantum chemistry calculations have been performed to calculate the reaction rates for the olefin +HO 2 elimination reactions of fuel peroxyl radical decomposition reactions, which were suggested as critical reaction pathways for the oxidation of cyclopentanone in previous studies. Based on these experimental and theoretical investigations, a detailed kinetic model has been developed and validated using the experimental data. Here, the model has satisfactorily reproduced the ignition delay times in the RCM and shock tube, and CO histories in the shock tube over the wide range of temperature, pressure and equivalence ratio. Rate of production and sensitivity analyses were performed to determine the important reaction pathways and critical reactions that affect the predicted reactivity of cyclopentanone at the condition investigated.

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

When more data hurts: Optimizing data coverage while mitigating diversity-induced underfitting in an ultrafast machine-learned potential

Machine-learned interatomic potentials (MLIPs) are becoming an essential tool in materials modeling. However, optimizing the generation of training data used to parametrize the MLIPs remains a significant challenge. This is because MLIPs can fail when encountering local environments too different from those present in the training data. The difficulty of determining a priori the environments that will be encountered during molecular dynamics simulation necessitates diverse, high-quality training data. Here, this study investigates how training data diversity affects the performance of MLIPs using the Ultra-Fast force field (UF 3 ) to model amorphous silicon nitride. We employ expert and autonomously generated data to create the training data and fit four force field variants to subsets of the data. Our findings reveal a critical balance in training data diversity: insufficient diversity hinders generalization, while excessive diversity can exceed the MLIP's learning capacity, reducing simulation accuracy. Specifically, we found that the UF 3 variant trained on a subset of the training data, in which nitrogen-rich structures were removed, offered vastly better prediction and simulation accuracy than any other variant. By comparing these UF 3 variants, we highlight the nuanced requirements for creating accurate MLIPs, emphasizing the importance of application-specific training data to achieve optimal performance in modeling complex material behaviors.

ab initio molecular dynamics↗
DOE OSTI2015

Solving coiled-coil protein structures

With the availability of more than 100,000 entries stored in the Protein Data Bank (PDB) that can be used as search models, molecular replacement (MR) is currently the most popular method of solving crystal structures of macromolecules. Significant methodological efforts have been directed in recent years towards making this approach more powerful and practical. This resulted in the creation of several computer programs, highly automated and user friendly, that are able to successfully solve many structures even by researchers who, although interested in structures of biomolecules, are not very experienced in crystallography.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗