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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Interaction and dynamics of ionic liquids based on choline and amino acid anions

The combination of amino acid anions with the choline cation gives origin to a new and potentially important class of organic ionic liquids that might represent a viable and bio-compatible alternative with respect to the traditional ones. We present here a detailed study of the bulk phase of the prototype system composed of the simplest amino acid (alanine) anion and the choline cation, based on ab initio and classical molecular dynamics. Theoretical findings have been validated by comparing with accurate experimental X-ray diffraction data and infrared spectra. We find that hydrogen bonding (HB) features in these systems are crucial in establishing their local geometric structure. We have also found that these HBs once formed are persistent and that the proton resides exclusively on the choline cation. In addition, we show that a classical force field description for this particular ionic liquid can be accurately performed by using a slightly modified version of the generalized AMBER force field.

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

FiXR : a framework to reconstruct fiber cross-sections from X-ray fiber diffraction experiments

Ab initio reconstruction methods have revolutionized the capabilities of small-angle X-ray scattering (SAXS), allowing the data-driven discovery of previously unknown molecular conformations, exploiting optimization heuristics and assumptions behind the composition of globular molecules. While these methods have been successful for the analysis of small particles, their impact on fibrillar assemblies has been more limited. The micrometre-range size of these assemblies and the complex interaction of their periodicities in their scattering profiles indicate that the discovery of fibril structures from SAXS measurements requires novel approaches beyond extending existing tools for molecular discovery. In this work, it is proposed to use SAXS measurements, together with diffraction theory, to infer the electron distribution of the average cross-section of a fiber. This cross-section is modeled as a discrete electron density with continuous support, allowing representations beyond binary distributions. Additional constraints, such as non-negativity or smoothness/connectedness, can also be added to the framework. The proposed approach is tested using simulated SAXS data from amyloid β fibril models and using measured data of Tobacco mosaic virus from SAXS experiments, recovering the geometry and density of the cross-sections in all cases. The approach is further tested by analyzing SAXS data from different amyloid β fibril assemblies, with results that are in agreement with previously proposed models from cryo-EM measurements. The limitations of the proposed method, together with an analysis of the robustness of the method and the combination with different experimental sources, are also discussed.

Biochemistry & Molecular Biology↗
DOE OSTI2021

ATOMIC Simulations and Experimental Data for CaCO3 Mixtures

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [1] and are used to develop a statistical characterization of matrix effects. Our overall suite of simulations includes contains several sets of simulations: training and validation sets of simulations for three and four element mixtures of calcium, carbon, oxygen, and nitrogen (included to account for atmosphere) along with simulations of the individual elements. The 4-element simulations include the mixture of all four elements mentioned and for each of the four individual elements. The 3-element simulations include output for the mixture of calcium, carbon, oxygen and for these three individual elements. The training data were produced using a 600-run design, shown in Figure 1, that varies input parameters temperature (T), electron density (Ne), and proportion of the elements calcium, carbon, oxygen, and nitrogen (Ca; C; O; N) for the 4 element output. The 3-element output includes all parameters except for the proportion of nitrogen. The element proportions (all the variables but T and Ne) sum to one and are unused in the single-element simulations. The validation data was produced with a 80-run design shown in Figure 2. The training and validation simulation outputs for the 4-element simulations for the mixture and for the single element calcium are shown as sample simulations in Figures 3 and 4 respectively. The simulations produce spectra over a range of 190nm - 950nm that roughly mimics the range collected by the SciAps Z-300 LIBS instrument that was used for the experimental data. The measured spectra for a CaCO3 (which may include contribution from Earth's atmosphere) in the experiment is shown in in Figure 5. All files are kept in directories whose names indicate the elemental composition (CaCO3, Ca, C, O, or N), number of elements (3 or 4), and purpose (training, which is not labeled in the file name, or validation) with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file. [1] J Colgan, EJ Judge, DP Kilcrease, and JE Barefield II. Ab-initio modeling of an iron laser-induced plasma: Comparison between theoretical and experimental atomic emission spectra. Spectrochimica Acta Part B: Atomic Spectroscopy, 97:65{73}, 2014.

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

Interstellar Isomers: The Importance of Bonding Energy Differences

We present strong detections of methyl cyanide (CH3CN), vinyl cyanide (CH2CHCN), ethyl cyanide (CH3CH2CN) and cyanodiacetylene (HC4CN) molecules with the Green Bank Telescope (GBT) toward the Sgr B2(N) molecular cloud. Attempts to detect the corresponding isocyanide isomers were only successful in the case of methyl isocyanide (CH3NC) for its J(sub K) = 1(sub 0) - 0(sub 0) transition, which is the first interstellar report of this line. To determine the spatial distribution of CH3NC, we used archival Berkeley-Illinois-Maryland Association (BIMA) array data for the J(sub K) = 1(sub 0) - 0(sub 0) transitions but no emission was detected. From ab initio calculations, the bonding energy difference between the cyanide and isocyanide molecules is greater than 8500 per centimeter (greater than 12,000 K). Thus, cyanides are the more stable isomers and would likely be formed more preferentially over their isocyanide counterparts. That we detect CH3NC emission with a single antenna (Gaussian beamsize(omega(sub B))=1723 arcsec(sup 2)) but not with an interferometer (omega(sub b)=192 arcsec(sup 2)), strongly suggests that CH3NC has a widespread spatial distribution toward the Sgr B2(N) region. Other investigators have shown that CH3CN is present both in the LMH hot core of Sgr B2(N) and in the surrounding medium, while we have shown that CH3NC appears to be deficient in the LMH hot core. Thus, largescale, non-thermal processes in the surrounding medium may account for the conversion of CH3CN to CH3NC while the LMH hot core, which is dominated by thermal processes, does not produce a significant amount of CH3NC. Ice analog experiments by other investigators have shown that radiation bombardment of CH3CN can produce CH3NC, thus supporting our observations. We conclude that isomers separated by such large bonding energy differences are distributed in different interstellar environments, making the evaluation of column density ratios between such isomers irrelevant unless it can be independently shown that these species are co-spatial.

Remijan, Anthony J.↗
DOE OSTI2024

Evaluation of Thermal Scattering Law for U 3 Si 2 Using Temperature-Dependent Phonon Density of States

The thermal scattering law for uranium silicide, a candidate material for high-density nuclear fuel, was evaluated from the temperature-dependent phonon density of states obtained using the quasi-harmonic approximation. For uranium silicide model systems of various volumes near the equilibrium volume, Hellmann-Feynman forces were calculated by performing ab initio lattice dynamics simulations based on density functional theory. Subsequently, the quasi-harmonic approximation was applied by interpolating the results of the harmonic approximation for each volume system. As a sensitivity test, the phonon density of states was calculated for the structures with temperature-dependent volumes at 298 K and 1000 K, obtained by the quasi-harmonic approximation, and the corresponding thermal neutron scattering cross section was evaluated. The temperature-dependent volume and heat capacity of uranium silicide obtained by the quasi-harmonic approximation showed reasonable agreement with experimental data while lattice thermal expansion showed deviation from experimental observations. Thermal scattering cross sections evaluated from 0 K and temperature-dependent phonon density of states exhibited differences of up to 4.1% at 1000 K. This study examines the impact of temperature on the thermal scattering law, employing ab initio lattice dynamics simulations through the quasi-harmonic approximation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Local cation ordering in compositionally complex Ruddlesden–Popper n = 1 oxides

The Ruddlesden–Popper (RP) layered perovskite structure is of great interest due to its inherent tunability, and the emergence and growth of the compositionally complex oxide (CCO) concept endows the RP family with further possibilities. Here, a comprehensive assessment of thermodynamic stabilization, local order/disorder, and lattice distortion was performed in the first two reported examples of lanthanum-deficient La n+1 B n O 3n+1 (n = 1, B = Mg, Co, Ni, Cu, Zn) obtained via various processing conditions. Chemical short-range order (CSRO) at the B-site and the controllable excess interstitial oxygen (δ) in RP-CCOs are uncovered by neutron pair distribution function analysis. Reverse Monte Carlo analysis of the data, Metropolis Monte Carlo simulations, and extended x-ray absorption fine structure analysis implies a modest degree of magnetic element segregation on the local scale. Further, ab initio molecular dynamics simulations results obtained from special quasirandom structure disagree with experimentally observed CSRO but confirm Jahn–Teller distortion of CuO 6 octahedra. These findings highlight potential opportunities to control local order/disorder and excess interstitial oxygen in layered RP-CCOs and demonstrate a high degree of freedom for tailoring application-specific properties. They also suggest a need for expansion of theoretical and data modeling approaches in order to meet the innate challenges of CCO and related high-entropy phases.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Quantum molecular dynamics simulation of shock-wave experiments in aluminum

We present quantum molecular dynamics calculations of principal, porous, and double shock Hugoniots, release isentropes, and sound velocity behind the shock front for aluminum. A comprehensive analysis of available shock-wave data is performed; the agreement and discrepancies of simulation results with measurements are discussed. Special attention is paid to the melting region of aluminum along the principal Hugoniot; the boundaries of the melting zone are estimated using the self-diffusion coefficient. Also, we make a comparison with a high-quality multiphase equation of state for aluminum. Independent semiempirical and first-principle models are very close to each other in caloric variables (pressure, density, particle velocity, etc.) but the equation of state gives higher temperature on the principal Hugoniot and release isentropes than ab initio calculations. Thus, the quantum molecular dynamics method can be used for calibration of semiempirical equations of state in case of lack of experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Adjusted Graphite Thermal Scattering Data Applied to Fluoride-Salt-Cooled Reactors

Thermal scattering data is a key component of fission reactor simulation. While most materials can be reasonably modeled with a free gas approximation, the moderators in thermal systems require a more rigorous treatment to sufficiently capture their specific neutron interaction physics. Graphite does not typically feature as a principle moderator in power reactors in service today, but remains an attractive prospect for a number of new reactor concepts currently in development. In comparison to light water, a material nearly ubiquitous in thermal reactors, graphite offers a better neutron economy by way of a small absorption cross section, high thermodynamic stability and can be incorporated into the structure of the core itself. Graphite in its purest form is a crystalline material composed of stacked layers of graphene, a hexagonal planar lattice of carbon atoms. In-plane, each carbon atom is covalently bonded to 3 neighbouring carbon atoms, creating a strong bond, while between planes only a weak van der Waals force exists. Pure graphite has a density of 2.26 g/cm{sup 3}. However, reactor-grade contains a significant amount of carbon binder, leading to a porous bulk medium with randomly oriented micro crystals. Densities of 1.5 - 1.9 g/cm{sup 3} are typical for reactor-grade graphite, corresponding to porosities on the order of 10-30%. This deviation from pure crystalline graphite should be accounted for in the thermal scattering data. When an incident neutron is in or near the thermal region, its de Broglie wavelength is the same order as the molecular bond lengths, leading to elastic scattering. At similar energies, the dynamics of the moderator lattice (in the case of solid moderators) can cause inelastic scattering. For these moderators, special cross sections should be calculated and used in reactor simulations. To generate this thermal scattering data, parameters pertaining to the properties and structure of the material in question are fed into a nuclear data processing code, such as NJOY. Within NJOY is a module called LEAPR which is used to generate the thermal scattering law S(α,β), which is in turn is fed into another module THERMR which is used to generate thermal cross sections that factor in the molecular-scale scattering phenomena. The core of the thermal scattering law is the phonon spectrum ρ(β) which describes the modes of excitation of the material lattice. An investigation into the phonon frequency spectrum was conducted at North Carolina State University (NCSU) for reactor-grade graphite which accounts for porosity. This phonon spectrum was calculated using the ab initio density functional theory (DFT) code VASP. The inelastic cross section calculated using this spectrum was found to have good agreement with relevant experimental data. The resulting spectra differed from the benchmark spectrum produced by Young and Koppel, the same data used in current nuclear data evaluations. The latter uses the lattice dynamics of perfect crystalline graphite with four force constants fitted to the model from thermodynamics data from reactor-grade (porous) graphite, which yields a reasonable approximation to the phonon spectrum of reactor-grade graphite.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Determination of the binding energies of the np Rydberg states of H{sub 2}, HD, and D{sub 2} from high-resolution spectroscopic data by multichannel quantum-defect theory

Multichannel quantum-defect theory (MQDT) is used to calculate the electron binding energies of np Rydberg states of H{sub 2}, HD, and D{sub 2} around n = 60 at an accuracy of better than 0.5 MHz. The theory includes the effects of rovibronic channel interactions and the hyperfine structure, and has been extended to the calculation of the asymmetric hyperfine structure of Rydberg states of a heteronuclear diatomic molecule (HD). Starting values for the eigenquantum-defect parameters of MQDT were extracted from ab initio potential-energy functions for the low-lying p Rydberg states of molecular hydrogen and subsequently refined in a global weighted fit to available experimental data on the singlet and triplet Rydberg states of H{sub 2} and D{sub 2}. The electron binding energies of high-np Rydberg states derived in this work represent important quantities for future determinations of the adiabatic ionization energies of H{sub 2}, HD, and D{sub 2} at sub-MHz accuracy.

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

Reference Correlation for the Viscosity of Carbon Dioxide

A comprehensive database of experimental and computed data for the viscosity of carbon dioxide (CO 2 ) was compiled and a new reference correlation was developed. Literature results based on an ab initio potential energy surface were the foundation of the correlation of the viscosity in the limit of zero density in the temperature range from 100 K to 2000 K. Guided symbolic regression was employed to obtain a new functional form that extrapolates correctly to T → 0 K and to 10 000 K. Coordinated measurements at low density made it possible to implement the temperature dependence of the Rainwater-Friend theory in the linear-in-density viscosity term. The residual viscosity could be formulated with a scaling term ργ/T the significance of which was confirmed by symbolic regression. The final viscosity correlation covers temperatures from 100 K to 2000 K for gaseous CO 2 , and from 220 K to 700 K with pressures along the melting line up to 8000 MPa for compressed and supercritical liquid states. The data representation is more accurate than with the previous correlations, and the covered pressure and temperature range is significantly extended. The critical enhancement of the viscosity of CO 2 is included in the new correlation.

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

Magnetic properties of Ni 40+x Mn 39-x Sn 21 (x = 0, 2, 4, 6 and 8 at.%) Heusler alloys

The low electron concentration region (e/a < 7.75) of the magnetic phase diagram of the off-stoichiometric Ni–Mn–Sn Heusler alloys was investigated in detail by DSC and magnetization measurements of the Ni 40+x Mn 39-x Sn 21 (x = 0, 2, 4, 6 and 8 at.%) alloys. The alloys show a stable austenitic phase without any martensitic transformation down to 5 K even after heat treatment. The Curie temperature exhibits a broad maximum over a large composition range. The evolution of the magnetic moment with the electron concentration fits the data of previous studies and confirms the peak-like dependence in the extended range of e/a values predicted by ab initio calculations. Finally, the explored part of the moment versus e/a curve can be explained in terms of a localized magnetic moment model and full atomic order in the alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Convergent Protocols for Computing Protein–Ligand Interaction Energies Using Fragment-Based Quantum Chemistry

Fragment-based quantum chemistry methods offer a way to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein–ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost. We also demonstrate that the HF-3c results are illustrative of trends obtained with density functional theory in basis sets up to augmented quadruple-ζ quality. Strategic deployment of fragmentation facilitates the use of converged biomolecular model systems alongside high-quality electronic structure methods and basis sets, bringing ab initio quantum chemistry to systems of hitherto unimaginable size. This will be useful for generation of high-quality training data for machine learning applications.

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

Shock melting method to determine melting curve by molecular dynamics: Cu, Pd, and Al

A melting simulation method, the shock melting (SM) method, is proposed and proved to be able to determine the melting curves of materials accurately and efficiently. The SM method, which is based on the multi-scale shock technique, determines melting curves by preheating and/or prepressurizing materials before shock. This strategy was extensively verified using both classical and ab initio molecular dynamics (MD). First, the SM method yielded the same satisfactory melting curve of Cu with only 360 atoms using classical MD, compared to the results from the Z-method and the two-phase coexistence method. Then, it also produced a satisfactory melting curve of Pd with only 756 atoms. Finally, the SM method combined with ab initio MD cheaply achieved a good melting curve of Al with only 180 atoms, which agrees well with the experimental data and the calculated results from other methods. It turned out that the SM method is an alternative efficient method for calculating the melting curves of materials.

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

New interaction potentials for alkali and alkaline-earth aluminosilicate glasses

Here, we apply a recently developed optimization scheme to obtain effective potentials for alkali and alkaline-earth aluminosilicate glasses that contain lithium, sodium, potassium, or calcium as modifiers. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by means of ab initio molecular dynamics simulations and density and elastic modulus of glass at room temperature from experiments. The new interaction potentials are able to reproduce reliably the structure and various mechanical and vibrational properties over a wide range of compositions for binary silicates. Previously, we tested these potentials for various ternary systems and find that they are transferable and can be mixed, thus allowing us to reproduce and predict the structure and properties of multicomponent glasses.

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

Ab Initio Quantum Monte Carlo Simulation of the Warm Dense Electron Gas in the Thermodynamic Limit

Here we perform ab initio quantum Monte Carlo (QMC) simulations of the warm dense uniform electron gas in the thermodynamic limit. By combining QMC data with the linear response theory, we are able to remove finite-size errors from the potential energy over the substantial parts of the warm dense regime, overcoming the deficiencies of the existing finite-size corrections by Brown et al. [Phys. Rev. Lett. 110, 146405 (2013)]. Extensive new QMC results for up to N = 1000 electrons enable us to compute the potential energy V and the exchange-correlation free energy F xc of the macroscopic electron gas with an unprecedented accuracy of | Δ V | / | V | , | Δ F xc | / | F | xc ~ 10 $-$3 . Finally, a comparison of our new data to the recent parametrization of F xc by Karasiev et al. [Phys. Rev. Lett. 112, 076403 (2014)] reveals significant deviations to the latter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

ζ-Glycine: insight into the mechanism of a polymorphic phase transition

Glycine is the simplest and most polymorphic amino acid, with five phases having been structurally characterized at atmospheric or high pressure. A sixth form, the elusive ζ phase, was discovered over a decade ago as a short-lived intermediate which formed as the high-pressure ϵ phase transformed to the γ form on decompression. However, its structure has remained unsolved. We now report the structure of the ζ phase, which was trapped at 100 K enabling neutron powder diffraction data to be obtained. The structure was solved using the results of a crystal structure prediction procedure based on fully ab initio energy calculations combined with a genetic algorithm for searching phase space. We show that the fate of ζ-glycine depends on its thermal history: although at room temperature it transforms back to the γ phase, warming the sample from 100 K to room temperature yielded β-glycine, the least stable of the known ambient-pressure polymorphs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermodynamics of Plutonium Monocarbide from Anharmonic and Relativistic Theory

Thermodynamics of plutonium monocarbide is studied from first-principles theory that includes relativistic electronic structure and anharmonic lattice vibrations. Density-functional theory (DFT) is expanded to include orbital-orbital coupling in addition to the relativistic spin-orbit interaction for the electronic structure and it is combined with anharmonic, temperature dependent, lattice dynamics derived from the self-consistent ab initio lattice dynamics (SCAILD) method. The obtained thermodynamics are compared to results from simpler quasi-harmonic theory and experimental data. Formation enthalpy, specific heat, and Gibbs energy calculated from the anharmonic model are validated by direct comparison with a calculation of phase diagram (CALPHAD) assessment of PuC and sub-stochiometric PuC0.896. Overall, the theory reproduces CALPHAD results and measured data for PuC rather well, but the comparison is hampered by the sub-stoichiometric nature of plutonium monocarbide. It was found that a bare theoretical approach that ignores spin-orbit and orbital-orbital coupling (orbital polarization) of the plutonium 5f electrons promotes too soft phonons and Gibbs energies that are incompatible with that of the CALPHAD assessment of the experimental data. The investigation of PuC suggests that the electronic structure is accurately described by plutonium 5f electrons as “band like” and delocalized, but correlate through spin polarization, orbital polarization, and spin-orbit coupling, in analogy to previous findings for plutonium metal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗