Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computing methods”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tutorial: Extracting entanglement signatures from neutron spectroscopy

This tutorial is a pedagogical introduction to recent methods of computing quantum spin entanglement witnesses from spectroscopy, with a special focus on neutron scattering on quantum spin systems. We offer a brief introduction to the concepts and equations, define a data analysis protocol, and discuss the interpretation of three entanglement witnesses: one-tangle, two-tangle, and Quantum Fisher Information. We also discuss practical experimental considerations, and give three examples of extracting entanglement witnesses from experimental data: Copper Nitrate, KCuF 3 , and NiPS 3 .

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Monte Carlo simulations for understanding the transport properties of metallic nanowires

We report the conductance calculation of metallic nanowires of different metals under the application of a tensile force. The elongation induced by the tensile force was simulated by means of the canonical Monte Carlo method combined with the embedded atom method to compute the energy of the system. We studied nanowires composed of Ag, Au, Pt and Cu. The conductance was calculated by means of the Landauer formula and the non-equilibrium Green function formalism. As previously verified experimentally and using other simulation techniques, the conductance is quantized. In addition to this, we have found that there is a strong correlation between the magnitude of the conductance and the width of the nanowire at its thinnest point. Conductance histograms were analyzed for all the different metals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

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

ab initio Sub-Mechanism Development for Cyclopentene Oxidation

To accurately predict low-temperature oxidation behavior, chemical kinetics mechanisms must contain complete reaction networks that include detailed consumption reactions of intermediates produced directly from hydroperoxyalkyl radicals, Q̇OOH, which undergo competing unimolecular reactions and bimolecular reactions with O2. Rates of chain-branching are governed by the flux between the two competing pathways, and inherently depend on temperature, pressure, and oxygen concentration. Neglect of consumption pathways for major oxidation intermediates leads to mechanism truncation error that is ameliorated by expanding the level of detail included in sub-mechanisms and employing ab initio methods for computing rates of elementary reactions and thermochemical properties of species involved. In the present work, an ab initio-derived sub-mechanism is developed using AutoMech to model the chemical kinetics of cyclopentene, a major product of cyclopentane oxidation. The ab initio sub-mechanism builds on a detailed mechanism developed using Reaction Mechanism Generator (RMG) for the specific purpose of determining the extent to which replacing cyclopentene-specific reactions and species with quantum chemical computations reduces model inaccuracies resulting from mechanism truncation error. In an effort to minimize interference from other reactions present during the formation of cyclopentene from cyclopentyl + O2, providing a narrower experimental scope, the model is compared against speciation measurements from jet-stirred reactor (JSR) experiments on cyclopentene oxidation. The experiments utilize vacuum ultraviolet-absorption spectroscopy and mass spectrometry for isomer-resolved speciation of intermediates at 835 Torr from 700 – 950 K. [O2]-dependent experiments were also conducted from 0.057 – 2.01 · 1018 molecules cm–3 at 825 K to examine the influence of oxygen on species profiles. Model predictions using the ab initio-revised mechanism yielded significant improvements in species profiles for both the temperature- and [O2]-dependent measurements, owing in part to increased rates of HOȮ and H2O2 production, which underscores the influence of theoretical calculations of reaction rates involving species produced from Ṙ + O2 such as cyclopentene.

AutoMech↗
DOE OSTI2020

Atomic resolution convergent beam electron diffraction analysis using convolutional neural networks

Two types of convolutional neural network (CNN) models, a discrete classification network and a continuous regression network, were trained to determine local sample thickness from convergent beam diffraction (CBED) patterns of SrTiO 3 collected in a scanning transmission electron microscope (STEM) at atomic column resolution. Acquisition of atomic resolution CBED patterns for this purpose requires careful balancing of CBED feature size in pixels, acquisition speed, and detector dynamic range. The training datasets were derived from multislice simulations, which must be convolved with incoherent source broadening. Sample thicknesses were also determined using quantitative high-angle annular dark-field (HAADF) STEM images acquired simultaneously. The regression CNN performed well on sample thinner than 35 nm, with 70% of the CNN results within 1 nm of HAADF thickness, and 1.0 nm overall root mean square error between the two measurements. The classification CNN was trained for a thicknesses up to 100 nm and yielded 66% of CNN results within one classification increment of 2 nm of HAADF thickness. Our approach depends on methods from computer vision including transfer learning and image augmentation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computing the viscous effect in early-time drop impact dynamics

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. ( Phys. Rev. Lett. , vol. 112, no. 13, 2014 b , p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.

42 ENGINEERING↗
DOE OSTI2021

Identification of Known and Novel Monomers for Poly(hydroxyurethanes) from Biobased Materials

Sourcing polymers from biobased materials is desirable for a transition to a more sustainable world. However, finding new monomers sourced from biobased molecules requires extensive experimental effort, due to the fact that most biomolecules lack the functional groups or properties needed to serve as a monomer in a particular polymeric system without chemical modification. In this work, we demonstrate a computational screening method to obtain monomers of a non-isocyanate polyurethane (NIPU), polyhydroxyurethane (PHU), from biomolecules. Our method generates candidate monomers from a list of 15 selected biobased starting molecules by performing chemical transformations that are standard to chemical catalysis and organic synthesis on these starting molecules and then screening the results for molecules with moieties that can be used as monomers for a PHU. Product molecules from two generations of reaction, a total number of 20,057, were analyzed automatically to identify molecules that have two or more epoxy or cyclic carbonate functional groups as candidates for further derivatization. These 21 candidate molecules were then subjected to another two generations of reaction, producing a set of 17,913 additional molecules. Further scrutiny of this set for desired functional groups characterizing PHUs and querying against the PubChem database was carried out, and 15 known molecules that are commercially available were identified, with the desirable features that they are predicted to be derived from biobased starting molecules in four or fewer steps. Moreover, molecule space generated from this process was primarily comprised of molecules that were not known to PubChem, indicating the framework also has the ability to generate novel candidate monomers that originate from biobased starting points.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Discovery of Novel Rhizoctonia solani DHFR Inhibitors as Fungicides Using Virtual Screening

Dihydrofolate reductase (DHFR) is an essential enzyme in the folate pathway and has been recognized as a well-known target for antibacterial and antifungal drugs. We discovered eight compounds from the ZINC database using virtual screening to inhibit Rhizoctonia solani (R. solani), a fungal pathogen in crops. These compounds were evaluated with in vitro assays for enzymatic and antifungal activity. Among these, compound Hit8 is the most active R. solani DHFR inhibitor, with the IC 50 of 10.2 μM. The selectivity of inhibition is 22.3 against human DHFR with the IC 50 of 227.7 μM. Moreover, Hit8 has higher antifungal activity against R. solani (EC 50 of 38.2 mg L –1 ) compared with validamycin A (EC 50 of 67.6 mg L –1 ), a well-documented fungicide. These results suggest that Hit8 may be a potential fungicide. Finally, our study exemplifies a computer-aided method to discover novel inhibitors that could target plant pathogenic fungi.

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

Toward DMC Accuracy Across Chemical Space with Scalable Δ-QML

In the past decade, quantum diffusion Monte Carlo (DMC) has been demonstrated to successfully predict the energetics and properties of a wide range of molecules and solids by numerically solving the electronic many-body Schrödinger equation. With O(N 3 ) scaling with the number of electrons N, DMC has the potential to be a reference method for larger systems that are not accessible to more traditional methods such as CCSD(T). Assessing the accuracy of DMC for smaller molecules becomes the stepping stone in making the method a reference for larger systems. We show that when coupled with quantum machine learning (QML)-based surrogate methods, the computational burden can be alleviated such that quantum Monte Carlo (QMC) shows clear potential to undergird the formation of high-quality descriptions across chemical space. We discuss three crucial approximations necessary to accomplish this: the fixed-node approximation, universal and accurate references for chemical bond dissociation energies, and scalable minimal amons-set-based QML (AQML) models. Numerical evidence presented includes converged DMC results for over 1000 small organic molecules with up to five heavy atoms used as amons and 50 medium-sized organic molecules with nine heavy atoms to validate the AQML predictions. Finally, numerical evidence collected for Δ-AQML models suggests that already modestly sized QMC training data sets of amons suffice to predict total energies with near chemical accuracy throughout chemical space.

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

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

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

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

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

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer ice-philicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface–ice contact at the expense of surface–water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Furthermore, our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

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

Adsorbate Partition Functions via Phase Space Integration: Quantifying the Effect of Translational Anharmonicity on Thermodynamic Properties

A new method for computing anharmonic thermophysical properties for adsorbates on metal surfaces is presented. Classical Monte Carlo phase space integration is performed to calculate the partition function for the motion of a hydrogen atom on Cu(111). Here, a minima-preserving neural network potential energy surface is used within the integration routine. Two different sampling schema for generating the training data are presented, and two different density functionals are used. The results are benchmarked against direct state counting results by using discrete variable representation. The phase space integration results are in excellent quantitative agreement with the benchmark results. Additionally, both the discrete variable representation and the phase space integration results confirm that the motion of H on Cu(111) is highly anharmonic. The results were applied to calculate the free energy of dissociative adsorption of H 2 and the resulting Langmuir isotherms at 400, 800, and 1200 K in a partial pressure range of 0–1 bar. It shows that the anharmonic effects lead to significantly higher predicted surface site fractions of hydrogen.

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

Accurate and Efficient Prediction of p K w in Aqueous Electrolytes Using Local Electrostatic Potentials

The K w =10 -pK w in aqueous electrolytes exhibit significant variations as a function of ion concentration. This behavior is not well-described by the Pitzer model, particularly at high concentrations. Here, we provide a molecular interpretation of the concentration dependent trends in p K w and develop a new method that accurately predicts the dissociation constant based upon the local electrostatic potential of the water oxygen that is determined by its nearest neighbors. The method is computationally efficient, relying only upon sampling of the local configurations (via molecular dynamics) and a small set of experimentally measured p K w . It is extensible to a range of salt compositions and the molecular understanding of the local potentials provides new routes toward the development of electrolytes with tailored physicochemical properties like p K w .

anions↗
DOE OSTI2022

Trends in Siting of Metals in Heterometallic Nd–Yb Metal–Organic Frameworks and Molecular Crystals

Several studies suggest that metal ordering within metal–organic frameworks (MOFs) is important for understanding how MOFs behave in relevant applications; however, these siting trends can be difficult to determine experimentally. To garner insight into the energetic driving forces that may lead to nonrandom ordering within heterometallic MOFs, we employ density functional theory (DFT) calculations on several bimetallic metal–organic crystals composed of Nd and Yb metal atoms. We also investigate the metal siting trends for a newly synthesized MOF. Our DFT-based energy of mixing results suggest that Nd will likely occupy sites with greater access to electronegative atoms and that local homometallic domains within a mixed-metal Nd–Yb system are favored. We also explore the use of less computationally extensive methods such as classical force fields and cluster expansion models to understand their feasibility for large system sizes. Here, this study highlights the impact of metal ordering on the energetic stability of heterometallic MOFs and crystal structures.

energy of mixing↗
DOE OSTI2022

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

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

Atomic Resolution Cryogenic 4D-STEM Imaging via Robust Distortion Correction

Cryogenic four-dimensional scanning transmission electron microscopy (4D-STEM) imaging is a useful technique for studying quantum materials and their interfaces by simultaneously probing charge, lattice, spin, and chemistry on the atomic scale with the sample held at temperatures ranging from room to cryogenic. However, its applications are currently limited by the instabilities of cryo-stages and electronics. To overcome this challenge, we develop an algorithm to effectively correct the complex distortions present in atomic resolution cryogenic 4D-STEM data sets. This method uses nonrigid registration to identify localized distortions in a 4D-STEM and relate them to an undistorted experimental STEM image, followed by a series of affine transformations for distortion corrections. This method allows a minimum loss of information in both reciprocal and real spaces, enabling the reconstruction of sample information from 4D-STEM data sets. This method is computationally cheap, fast, and applicable for on-the-fly data analysis in future in situ cryogenic 4D-STEM experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗