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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Quantum and quasi-classical dynamics of the OH + CO → H + CO{sub 2} reaction on a new permutationally invariant neural network potential energy surface

A permutationally invariant global potential energy surface for the HOCO system is reported by fitting a larger number of high-level ab initio points using the newly proposed permutation invariant polynomial-neural network method. The small fitting error (∼5 meV) indicates a faithful representation of the potential energy surface over a large configuration space. Full-dimensional quantum and quasi-classical trajectory studies of the title reaction were performed on this potential energy surface. While the results suggest that the differences between this and an earlier neural network fits are small, discrepancies with state-to-state experimental data remain significant.

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

A nine-dimensional ab initio global potential energy surface for the H{sub 2}O{sup +} + H{sub 2} → H{sub 3}O{sup +} + H reaction

An accurate full-dimensional global potential energy surface (PES) is developed for the title reaction. While the long-range interactions in the reactant asymptote are represented by an analytical expression, the interaction region of the PES is fit to more than 81 000 of ab initio points at the UCCSD(T)-F12b/AVTZ level using the permutation invariant polynomial neural network approach. Fully symmetric with respect to permutation of all four hydrogen atoms, the PES provides a faithful representation of the ab initio points, with a root mean square error of 1.8 meV or 15 cm{sup −1}. The reaction path for this exoergic reaction features an attractive and barrierless entrance channel, a submerged saddle point, a shallow H{sub 4}O{sup +} well, and a barrierless exit channel. The rate coefficients for the title reaction and kinetic isotope effect have been determined on this PES using quasi-classical trajectories, and they are in good agreement with available experimental data. It is further shown that the H{sub 2}O{sup +} rotational enhancement of reactivity observed experimentally can be traced to the submerged saddle point. Using our recently proposed Sudden Vector Projection model, we demonstrate that a rotational degree of freedom of the H{sub 2}O{sup +} reactant is strongly coupled with the reaction coordinate at this saddle point, thus unraveling the origin of the pronounced mode specificity in this reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deprotonated sulfamic acid and its homodimers: Does sulfamic acid adopt zwitterion during cluster growth?

Here we present a joint experimental and computational study on the geometric and electronic structures of deprotonated sulfamic acid (SA) clusters [(SA) n –H] – (n = 1, 2) employing negative ion photoelectron spectroscopy and high-level ab initio calculations. The photoelectron spectra provide the vertical/adiabatic detachment energy (VDE/ADE) of the sulfamate anion (SM – ) H 2 N•SO 3 – at 4.85 ± 0.05 and 4.58 ± 0.08 eV, respectively, and the VDE and ADE of the SM – •SA dimer at 6.41 ± 0.05 and 5.87 ± 0.08 eV, respectively. The significantly increased electron binding energies of the dimer confirm the enhanced electronic stability upon the addition of one SA molecule. The CCSD(T)-predicted VDEs/ADEs agree excellently with the experimental data, confirming the identified structures as the most stable ones. Two types of dimer isomers possessing different hydrogen bonding (HB) motifs are identified, corresponding to SM – binding to a zwitterionic SA (SM – •SA z ) and a canonical SA (SM – •SA c ), respectively. Two N–H$\cdots$O HBs and one superior O–H$\cdots$O HB are formed in the lowest-lying SM – •SA c , while SM – •SA z has three moderate N–H$\cdots$O HBs, with the former being 4.71 kcal/mol more stable. Further theoretical analyses reveal that the binding strength advantage of SM – •SA c over SM – •SA z arises from its significant contributions of orbital interactions between fragments, illustrating that sulfamate strongly interacts with its parent SA acid and preferably chooses the canonical SA in the subsequent cluster formations. Given the prominent presence of SA, this study provides the first evidence that the canonical dimer model of sulfamic acid should exist as a superior configuration during cluster growth.

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

Quantum Monte Carlo calculations of weak transitions in A = 6 – 10 nuclei

{\it Ab initio} calculations of the Gamow-Teller (GT) matrix elements in the $\beta$ decays of $^6$He and $^{10}$C and electron captures in $^7$Be are carried out using both variational and Green's function Monte Carlo wave functions obtained from the Argonne $v_{18}$ two-nucleon and Illinois-7 three-nucleon interactions, and axial many-body currents derived from either meson-exchange phenomenology or chiral effective field theory. The agreement with experimental data is excellent for the electron captures in $^7$Be, while theory overestimates the $^6$He and $^{10}$C data by $\sim 2\%$ and $\sim 10\%$, respectively. We show that for these systems correlations in the nuclear wave functions are crucial to explain the data, while many-body currents increase by $\sim 2$--$3\%$ the one-body GT contributions. These findings suggest that the longstanding $g_A$-problem, {\it i.e.}, the systematic overprediction ($\sim 20 \%$ in $A\le 18$ nuclei) of GT matrix elements in shell-model calculations, may be resolved, at least partially, by correlation effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Ab initio Modeling of the Electronic Structures and Physical Properties of a‐Si 1− x Ge x O 2 Glass ( x = 0 to 1)

The amorphous silica (a‐SiO 2 ) and germania (a‐GeO 2 ) have a wide range of applications in glass industry. Based on a previously constructed near‐perfect continuous random network model with 1296 atoms and periodic boundary conditions, we extend our study to amorphous Si 1− x G e x O 2 models of homogeneous random substitution of Si by Ge with x ranging from 0 to 1. We have calculated the structural, electronic, mechanical, and optical properties for the series by using the first‐principles density functional theory methods. The x ‐dependence of the variations in the properties is analyzed and critically compared with available experimental data. The mass density, volume, total bond order density, bulk mechanical properties, and refractive index are found to vary linearly as a function of x . For x = 0.5, we have also constructed six different kinds of particle immersion models to test the effect of inclusion of spherical particles of one glass of different sizes in the medium of the other glass on their physical properties. It is shown that particle sizes do affect the properties of particle immersion. Our calculations provide deep insight on the properties of mixture and nanocomposites of a‐SiO 2 and a‐GeO 2 glasses.

Baral, Khagendra↗
DOE OSTI2014

Ab initio calculation of excess properties of La{sub 1−x}(Ln,An){sub x}PO{sub 4} solid solutions

We used ab initio computational approach to predict the excess enthalpy of mixing and the corresponding regular/subregular model parameters for La{sub 1−x}Ln{sub x}PO{sub 4} (Ln=Ce,…, Tb) and La{sub 1−x}An{sub x}PO{sub 4} (An=Pu, Am and Cm) monazite-type solid solutions. We found that the regular model interaction parameter W computed for La{sub 1−x}Ln{sub x}PO{sub 4} solid solutions matches the few existing experimental data. Within the lanthanide series W increases quadratically with the volume mismatch between LaPO{sub 4} and LnPO{sub 4} endmembers (ΔV=V{sub LaPO{sub 4}}−V{sub LnPO{sub 4}}), so that W(kJ/mol)=0.618(ΔV(cm{sup 3}/mol)){sup 2}. We demonstrate that this relationship also fits the interaction parameters computed for La{sub 1−x}An{sub x}PO{sub 4} solid solutions. This shows that lanthanides can be used as surrogates for investigation of the thermodynamic mixing properties of actinide-bearing solid solutions. - Highlights: • The excess enthalpies of mixing for monazite-type solid solutions are computed. • The excess enthalpies increase with the endmembers volume mismatch. • The relationship derived for lanthanides is transferable to La{sub 1−x}An{sub x}PO{sub 4} systems.

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

Ab initio study of the adsorption, diffusion, and intercalation of alkali metal atoms on the (0001) surface of the topological insulator Bi{sub 2}Se{sub 3}

Ab initio study of the adsorption, diffusion, and intercalation of alkali metal adatoms on the (0001) step surface of the topological insulator Bi{sub 2}Se{sub 3} has been performed for the case of low coverage. The calculations of the activation energies of diffusion of adatoms on the surface and in van der Waals gaps near steps, as well as the estimate of diffusion lengths, have shown that efficient intercalation through steps is possible only for Li and Na. Data obtained for K, Rb, and Cs atoms indicate that their thermal desorption at high temperatures can occur before intercalation. The results have been discussed in the context of existing experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

On the electronic structure and spectrum of CuF2 and CuCl2

Results are reported from SDCI and CPF ab initio computations for the ground and low-lying states of CuF2 and CuCl2. A large Gaussian basis set with functions up to the g level is employed, and the results are presented in tables and characterized in detail. Particular attention is given to the bond distances and force constants for the ligand-field states and the bond distances for the charge-transfer states. Good agreement is obtained between the CPF and SDCI values and with published experimental data.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2025

Consistent thermodynamic properties for alicyclic components of jet fuels: Experimental data, estimation methods, and homologous series trends

Alkylcycloalkanes represent a significant fraction of jet fuel components. An evaluation of their thermodynamic properties, enthalpies of formation in liquid and gas phases and enthalpies of vaporization, was conducted. A combination of available experimental data, up-to-date group-contribution methods, high-level quantum-chemical calculations, and homologous series trends was used to identify outliers and to recommend the most reliable values. The group-contribution approach was found to work well for the enthalpies of vaporization. Its performance for the enthalpies of formation in the liquid and gas phases was found to be substantially less effective, especially considering notable differences in this property among stereoisomers. Computationally affordable high-level ab initio results and homologous series trend analysis appeared more reliable. In conclusion, the recommended property values for 212 individual compounds and their isomeric mixtures were provided.

09 BIOMASS FUELS↗
DOE OSTI2018

Pressure induced short-range structural changes in supercooled liquid Ge 2 Sb 2 Te 5

Phase-change material such as Ge 2 Sb 2 Te 5 is usually utilized to store data due to the pronounced contrast in optical and electrical properties between crystalline and amorphous phases. As the density differs in the two phases, it is necessary to explore the influence of pressure on the structures of Ge 2 Sb 2 Te 5 , especially for the supercooled liquid which is an inevitable state in the formation of the two phases. The short-range structures in supercooled liquid Ge 2 Sb 2 Te 5 under compression have been investigated by using ab initio molecular dynamics simulation. The supercooled liquid eventually changes to a solid with an increase in pressure. During the process, tetrahedrons decrease slightly, revealing that tetrahedral structures are insensitive to the pressure. Octahedrons increase as the pressure is less than 3.7 GPa and then decrease, suggesting that a moderate pressure can promote the formation of octahedrons. The body-centered-cubic, face-centered-cubic and hexagonal closed-packed structures are observed at 9.9 GPa and then increase gradually. Additionally, it is noticed that Sb- and Te-centered clusters prefer Ge-centered clusters to forming the high-coordinated short-range structures during the compression process. Here our research makes us aware of the effects of pressure on short-range structures in supercooled liquid, which is of great importance in the application of phase-change materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Temperature-Dependent Properties of Molten Li 2 BeF 4 Salt Using Ab Initio Molecular Dynamics

Molten lithium tetrafluoroberyllate (Li 2 BeF 4 ) salt, also known as FLiBe, with a 2:1 mixture of LiF and BeF 2 is being proposed as a coolant and solvent in advanced nuclear reactor designs, such as the molten salt reactor or the fluoride salt cooled high-temperature reactor. We present the results on the structure and properties of FLiBe over a wide range of temperatures, 0–2000 K, from high-throughput ab initio molecular dynamics simulation using a supercell model of 504 atoms. The variations in the local structures of solid and liquid FLiBe with temperature are discussed in terms of a pair distribution function, coordination number, and bond angle distribution. The temperature-dependent electronic structure and optical and mechanical properties of FLiBe are calculated. The optical and mechanical property results are reported for the first time. The results above and below the melting temperature (~732 K) are compared with the experimental data and with data for crystalline FLiBe. The electronic structure and interatomic bonding results are discussed in correlation with the mechanical strength. A novel concept of total bond order density (TBOD), an important quantum mechanical parameter, is used to characterize the internal cohesion and strength in the simulated models. The results show a variation in the rate of change in properties in solid and liquid phases with anomalous behavior across the melting region. The observed trend is the decrease in mechanical strength, band gap, and TBOD in a nonlinear fashion as a function of temperature. The refractive index shows a surprising minimum at 850 K, among the tested temperatures, which lies above the melting point. These findings provide a new platform to understand the interplay between the temperature-dependent structures and properties of FLiBe salt.

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

Ab initio predictions for polarized deuterium-tritium thermonuclear fusion

The fusion of deuterium (D) with tritium (T) is the most promising of the reactions that could power thermonuclear reactors of the future. It may lead to even more efficient energy generation if obtained in a polarized state, that is with the spin of the reactants aligned. Here in this paper, we report first-principles predictions of the polarized DT fusion using nuclear forces from effective field theory. By employing the ab initio no-core shell model with continuum reaction method to solve the quantum mechanical five-nucleon problem, we accurately determine the enhanced fusion rate and angular distribution of the emitted neutron and 4 He. Our calculations demonstrate in detail the small contribution of anisotropies, placing on a firmer footing the understanding of the rate of DT fusion in a polarized plasma. In the future, analogous calculations could be used to obtain accurate values for other, more uncertain thermonuclear reaction data critical to nuclear science applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermophysical properties of FLiBe using moment tensor potentials

Fluoride salts are prospective materials for applications in some next-generation nuclear reactors and their thermophysical properties at various conditions are of interest. Experimental measurement of the properties of these salts is often difficult and, in some cases, unfeasible due to challenges from high temperatures, impurity control, and corrosivity. Therefore, accurate theoretical methods are needed for fluoride salt property prediction. In this work, we used moment tensor potentials (MTP) to approximate the potential energy surface of eutectic FLiBe (66.6% LiF – 33.3% BeF2) predicted by the ab initio (DFT-D3) method. Here, we then used the developed potential and molecular dynamics to obtain several thermophysical properties of FLiBe, including radial distribution functions, density, self-diffusion coefficients, thermal expansion, specific heat capacity, bulk modulus, viscosity, and thermal conductivity. Our results show that the MTP potential approximates the potential energy surface accurately and the overall approach yields very good agreement with experimental values. The converged fitting can be obtained with less than 600 configurations generated from DFT calculations, which data can be generated in just 1200 core hours on today's typical processors. The MTP potential is faster than many machine learning potentials and about one order of magnitude slower than widely used empirical molten salt potentials such as Tosi/Fumi.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atomic Energy Accuracy of Neural Network Potentials: Harnessing Pretraining and Transfer Learning

Machine learning-based interatomic potentials (MLIPs) have transformed the prediction of potential energy surfaces (PESs), achieving accuracy comparable to ab initio calculations. However, atomic energy predictions, often assumed to lack physical meaning, remain underexplored. In this study, we demonstrate that inaccuracies in atomic energy predictions reduce the robustness and transferability of Neural Network Potentials (NNPs) and atomic energy error can be masked in total energy predictions due to error cancellation. Here, we validate this finding using challenging configurations involving deformation and failure under tensile loading. By pretraining atomic energy predictions using empirical potentials and applying transfer learning with density functional theory (DFT) data, we achieve notable improvements in the accuracy of total energy, forces, and stress predictions. Furthermore, this approach enhances the robustness and transferability of NNPs, emphasizing the importance of atomic energy predictions in developing high-quality and reliable MLIPs.

Active Learning↗
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↗
NASA NTRS1985

Conformational properties, torsional potential, and vibrational force field for methacryloyl fluoride - An ab initio investigation

The structure, torsional potentials, vibrational spectra, and harmonic force fields for s-cis and s-trans isomers of methacryloyl fluoride are examined to understand the conformational properties of the molecules and their relationship to macroscopic polymer properties. The structure is found to be in good agreement with experiment. It is shown by calculations that the energy difference between the cis and the transisomers is less than 1 kcal/mol at both the split valence and the split valence polarized levels, with the trans form favored. Analysis of the torsional potentials indicates that a rigid rotor model provides a reasonable description of the motion of the COF group in the molecule. The torsional barrier to interconvert the s-trans to the s-cis form is found to be 7.0 kcal/mol. A fit of the data to a three-term Fourier series shows that it is possible to reproduce the experimentally derived barrier, even though a direct determination indicates that the barrier is higher.

Laskowski, B. C.↗
DOE OSTI2021

Persistent Homology Metrics Reveal Quantum Fluctuations and Reactive Atoms in Path Integral Dynamics

Nuclear quantum effects (NQEs) are known to impact a number of features associated with chemical reactivity and physicochemical properties, particularly for light atoms and at low temperatures. In the imaginary time path integral formalism, each atom is mapped onto a “ring polymer” whose spread is related to the quantum mechanical uncertainty in the particle’s position, i.e., its thermal wavelength. A number of metrics have previously been used to investigate and characterize this spread and explain effects arising from quantum delocalization, zero-point energy, and tunneling. Many of these shape metrics consider just the instantaneous structure of the ring polymers. However, given the significant interest in methods such as centroid molecular dynamics and ring polymer molecular dynamics that link the molecular dynamics of these ring polymers to real time properties, there exists significant opportunity to exploit metrics that also allow for the study of the fluctuations of the atom delocalization in time. Here we consider the ring polymer delocalization from the perspective of computational topology, specifically persistent homology, which describes the 3-dimensional arrangement of point cloud data, (i.e. atomic positions). We employ the Betti sequence probability distribution to define the ensemble of shapes adopted by the ring polymer. The Wasserstein distances of Betti sequences adjacent in time are used to characterize fluctuations in shape, where the Fourier transform and associated principal components provides added information differentiating atoms with different NQEs based on their dynamic properties. We demonstrate this methodology on two representative systems, a glassy system consisting of two atom types with dramatically different de Broglie thermal wavelengths, and ab initio molecular dynamics simulation of an aqueous 4 M HCl solution where the H-atoms are differentiated based on their participation in proton transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗