Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational 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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Heterogenized Molecular Catalysts: Vibrational Sum-Frequency Spectroscopic, Electrochemical, and Theoretical Investigations

Rhenium and manganese bipyridyl tricarbonyl complexes have attracted intense interest for their promising applications in photocatalytic and electrocatalytic CO 2 reduction in both homogeneous and heterogenized systems. To date, there have been extensive studies on immobilizing Re catalysts on solid surfaces for higher catalytic efficiency, reduced catalyst loading, and convenient product separation. However, in order for the heterogenized molecular catalysts to achieve the combination of the best aspects of homogeneous and heterogeneous catalysts, it is essential to understand the fundamental physicochemical properties of such heterogeneous systems, such as surface-bound structures of Re/Mn catalysts, substrate–adsorbate interactions, and photoinduced or electric-field-induced effects on Re/Mn catalysts. For example, the surface may act to (un)block substrates, (un)trap charges, (de)stabilize particular intermediates (and thus affect scaling relations), and shift potentials in different directions, just as protein environments do. The close collaboration between the Lian, Batista, and Kubiak groups has resulted in an integrated approach to investigate how the semiconductor or metal surface affects the properties of the attached catalyst. Synthetic strategies to achieve stable and controlled attachment of Re/Mn molecular catalysts have been developed. Steady-state, time-resolved, and electrochemical vibrational sum-frequency generation (SFG) spectroscopic studies have provided insight into the effects of interfacial structures, ultrafast vibrational energy relaxation, and electric field on the Re/Mn catalysts, respectively. Various computational methods utilizing density functional theory (DFT) have been developed and applied to determine the molecular orientation by direct comparison to spectroscopy, unravel vibrational energy relaxation mechanisms, and quantify the interfacial electric field strength of the Re/Mn catalyst systems. This Account starts with a discussion of the recent progress in determining the surface-bound structures of Re catalysts on semiconductor and Au surfaces by a combined vibrational SFG and DFT study. The effects of crystal facet, length of anchoring ligands, and doping of the semiconductor on the bound structures of Re catalysts and of the substrate itself are discussed. This is followed by a summary of the progress in understanding the vibrational relaxation (VR) dynamics of Re catalysts covalently adsorbed on semiconductor and metal surfaces. The VR processes of Re catalysts on TiO 2 films and TiO 2 single crystals and a Re catalyst tethered on Au, particularly the role of electron–hole pair (EHP)-induced coupling on the VR of the Re catalyst bound on Au, are discussed. The Account also summarizes recent studies in quantifying the electric field strength experienced by the catalytically active site of the Re/Mn catalyst bound on a Au electrode based on a combined electrochemical SFG and DFT study of the Stark tuning of the CO stretching modes of these catalysts. Finally, future research directions on surface-immobilized molecular catalyst systems are discussed.

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

Characterization of the In Vivo Deuteration of Native Phospholipids by Mass Spectrometry Yields Guidelines for Their Regiospecific Customization

Customization of deuterated biomolecules is vital for many advanced biological experiments including neutron scattering. However, because it is challenging to control the proportion and regiospecificity of deuterium incorporation in live systems, often only two or three synthetic lipids are mixed together to form simplistic model membranes. This limits the applicability and biological accuracy of the results generated with these synthetic membranes. Despite some limited prior examination of deuterating Escherichia coli lipids in vivo, this approach has not been widely implemented. In this report an extensive mass spectrometry-based profiling of E. coli phospholipid deuteration states with several different growth media was performed, and a computational method to describe deuterium distributions with a one-number summary is introduced. The deuteration states of 36 lipid species were quantitatively profiled in 15 different growth conditions, and tandem mass spectrometry was used to reveal deuterium localization. Regressions were employed to enable the prediction of lipid deuteration for untested conditions. Small-angle neutron scattering was performed on select deuterated lipid samples, which validated the deuteration states calculated from the mass spectral data. Based on these experiments, guidelines for the design of specifically deuterated phospholipids are described. This unlocks even greater capabilities from neutron-based techniques, enabling experiments that were formerly impossible.

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

Increasing the Scale of the Mass Spectrometry Query Language Compendium with Explainable AI

A significant bottleneck in metabolomics data interpretation is the effective use of domain knowledge to assign structural information based on fragmentation patterns. The mass spectrometry query language (MassQL) aims to make this process accessible and applicable across multiple analysis platforms. While advanced computational methods are capable of predicting compound structures from fragmentation data, AI/ML approaches often rely on complex, opaque criteria that are difficult to interpret or modify. As a result, their predictive patterns cannot be readily translated into human-readable rules, such as those used in MassQL. Here, in this study, we introduce ChemEcho, a machine learning embedding method that converts tandem mass spectrometry data into sparse feature vectors containing peak and neutral mass subformulae to enhance explainable AI/ML-based methods. An advantage of this approach is that decision trees trained using these feature vectors can be directly translated to MassQL. Using a battery of decision trees trained using ChemEcho embeddings to predict molecular attributes, we generated over 1500 MassQL queries for 765 molecular features and evaluated their precision and recall. From these queries, the 50 highest-performing queries were integrated into the MassQL compendium. This set of generated MassQL queries included environmentally and biologically relevant classes such as PFAS and molecules containing phosphate or sulfate substructures. To illustrate the impact these queries would have on a typical metabolomics experiment, these MassQL queries were applied to a public metabolomics data set─resulting in a marked increase in the structural information derived from tandem mass spectra. Access and reuse of these queries is expected to enhance structural annotation in untargeted experiments, leading to more specific claims and advancing many applications in metabolomics.

Harwood, Thomas V. [USDOE Joint Genome Institute (↗
DOE OSTI2018

Toward First-Principles Design of Organic Nonlinear Optical Materials: Crystal Structure Prediction and Halogen Bonding Impact on Hyperpolarizabilities of 2-Iodo-3-hydroxypyridine

Computational methods can potentially accelerate development of more efficient organic materials for second harmonic generation. Here, we test the method that includes the evolutionary algorithm for predicting crystal structure and prognosis of nonlinear optical properties based on the predicted structure. For this test, we selected 2-iodo-3-hydroxypyridine, which exhibits second harmonic generation intensity comparable to that of urea. We performed global minimization of the lattice energy and found the experimental structure when many-body dispersion correction is added to the density functional theory values. We analyzed geometric preferences of the halogen bonding in predicted virtual polymorphs. We also found linear correlation between the lengths of the iodine–iodine halogen bonds and calculated second order susceptibilities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dynamics of Hydroxyl Anions Promotes Lithium Ion Conduction in Antiperovskite Li 2 OHCl

Li 2 OHCl is an exemplar of the antiperovskite family of ionic conductors, for which high ionic conductivities have been reported, but in which the atomic-level mechanism of ion migration is unclear. The stable phase is both crystallographically defective and disordered, having similar to 1/3 of the Li sites vacant, while the presence of the OH(- )anion introduces the possibility of rotational disorder that may be coupled to cation migration. In this work, complementary experimental and computational methods are applied to understand the relationship between the crystal chemistry and ionic conductivity in Li 2 OHCl , which undergoes an orthorhombic to cubic phase transition near 311 K (approximate to 38 degrees C) and coincides with the more than a factor of 10 change in ionic conductivity (from 1.2 x 10 -5 mS/cm at 37 degrees C to 1.4 x 10 -3 mS/cm at 39 degrees C). X ray and neutron experiments conducted over the temperature range 20-200 degrees C, including diffraction, quasi-elastic neutron scattering (QENS), the maximum entropy method (MEM) analysis, and ab initio molecular dynamics (AIMD) simulations, together show conclusively that the high lithium ion conductivity of cubic Li2OHCl is correlated to "paddlewheel" rotation of the dynamic OH - anion. The present results suggest that in antiperovskites and derivative structures a high cation vacancy concentration combined with the presence of disordered molecular anions can lead to high cation mobility.

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

Acid Gas Stability of Zeolitic Imidazolate Frameworks: Generalized Kinetic and Thermodynamic Characteristics

Here, acid gases such as SO 2 and CO 2 are present in many environments in which the use of nanoporous metal-organic frameworks (MOFs) is envisaged. Among metal-organic frameworks, zeolitic imidazolate frameworks (ZIFs) have been extensively explored as membranes or adsorbents. However, there is little systematic knowledge of the effects of acid gas exposure on the structure of ZIFs, in particular the mechanistic aspects of ZIF degradation by acid gases as well as the effects of ZIF crystal topology and linker composition on their stability. Here we present a generalized and quantitative investigation of the kinetic and thermodynamic acid gas stability of a diverse range of ZIF materials. The stability of 16 ZIFs (of SOD, RHO, ANA, and GME topologies) under different environments – humid air, liquid water, and acid gases CO 2 and SO 2 (dry, humid, and aqueous) – are investigated by a suite of experimental and computational methods. The kinetics of ZIF degradation under exposure to humid SO 2 is studied in detail, and effective rate constants for acid gas degradation of ZIFs are reported for the first time. Remarkably, the kinetics of degradation of the diverse ZIFs correlate strongly with the linker pKa and ZIF water adsorption in a manner contrary to that expected from previous predictions in the literature. Furthermore, we find that the material ZIF-71 (RHO topology) shows much higher stability relative to the other ZIFs in humid SO 2 and CO 2 environments.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Soft X-ray Absorption Spectroscopy of Liquids and Solutions.

X-ray absorption spectroscopy (XAS) is an electronic absorption technique for which the initial state is a deeply buried core level. The photon energies corresponding to such transitions are governed primarily by the binding energies of the initial state. Because the binding energies of core electrons vary significantly among atomic species, this makes XAS an element-selective spectroscopy. Proper interpretation of XA spectra can provide detailed information on the local chemical and geometric environment of the target atom. The introduction of liquid microjet and flow cell technologies into XAS experiments has enabled the general study of liquid samples. Liquids studied to date include water, alcohols, and solutions with relevance to biology and energy technology. This Review summarizes the experimental techniques employed in XAS studies of liquid samples and computational methods used for interpretation of the resulting spectra and summarizes salient experiments and results obtained in the XAS investigations of liquids.

Smith, Jacob W↗
DOE OSTI2020

In Silico Guidance for In Vitro Androgen and Glucocorticoid Receptor ToxCast Assays

Molecular initiating events (MIEs) are key events in adverse outcome pathways (AOPs) that link molecular chemistry to target biology. As they are based in chemistry, these interactions are excellent targets for computational chemistry approaches to in silico modelling. In this work, we aim to link ligand chemical structure to MIEs for androgen receptor (AR) and glucocorticoid receptor (GR) binding using ToxCast data. This has been done using an automated computational algorithm to perform maximal common substructure searches on chemical binders for each target from the ToxCast dataset. The models developed show a high level of accuracy, correctly assigning 87.20% of AR binders and 96.81% of GR binders in a 25% test set using holdout cross-validation. The 2D structural alerts developed can be used as in silico models to predict these MIEs, and as guidance for in vitro ToxCast assays to confirm hits. These models can target such experimental work, reducing the number of assays to be performed to gain required toxicological insight. Development of these models has also allowed some structural alerts to be identified as predictors for agonist or antagonist behavior at the receptor target. This work represents a first step in using computational methods to guide and target experimental approaches.

Allen, Timothy H.↗
DOE OSTI2020

Perturbation of 1 J C,F Coupling in Carbon–Fluorine Bonds on Coordination to Lewis Acids: A Structural, Spectroscopic, and Computational Study

A lithiated m-terphenyl ligand bearing fluorine atoms at the ortho positions of the flanking aryl rings was synthesized and characterized using single crystal X-ray diffraction, variable-temperature multinuclear NMR spectroscopy, and computational methods. Here, changes in 1 J C,F on coordination to lithium as a spectroscopic observable parametrizing the strength of the C–F···Li interaction are described, and a general, qualitative relationship between C–F bond lengths, Δ 1 J C,F values, and the extent of C–F bond activation as a result of Lewis acid coordination is proposed.

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

Four-Coordinate Fe N 2 and Imido Complexes Supported by a Hemilabile NNC Heteroscorpionate Ligand

Inspired by mechanistic proposals for N 2 reduction at the nitrogenase FeMo cofactor, we report herein a new, strongly σ-donating heteroscorpionate ligand featuring two weak-field pyrazoles and an alkyl donor. This ligand supports four-coordinate Fe(I)-N 2 , Fe(II)-Cl, and Fe(III)-imido complexes, which we have characterized using a variety of spectroscopic and computational methods. Structural and quantum mechanical analysis reveal the nature of the Fe–C bonds to be essentially invariant between the complexes, with conversion between the (formally) low-valent Fe-N 2 and high-valent Fe-imido complexes mediated by pyrazole hemilability. This presents a useful strategy for substrate reduction at such low-coordinate centers and suggests a mechanism by which FeMoco might accommodate the binding of both π-acidic and π-basic nitrogenous substrates.

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

General Protocol for the Accurate Prediction of Molecular 13 C/ 1 H NMR Chemical Shifts via Machine Learning Augmented DFT

An accurate prediction of NMR chemical shifts at affordable computational cost is very important for different types of structural assignments in experimental studies. Density functional theory (DFT) and gauge-including atomic orbital (GIAO) are two of the most popular computational methods for NMR calculation, yet they often fail to resolve ambiguities in structural assignments. In this work, we present a new method that uses machine learning (ML) techniques (DFT + ML) that significantly increases the accuracy of 13 C/ 1 H NMR chemical shift prediction for a variety of organic molecules. The input of the generalizable DFT + ML model contains two critical parts: one is a vector providing insights into chemical environments, which can be evaluated without knowing the exact geometry of the molecule; the other one is the DFT-calculated isotropic shielding constant. The DFT + ML model was trained with a data set containing 476 13 C and 270 1 H experimental chemical shifts. For the DFT methods used here, the root mean square deviations (RMSDs) for the errors between predicted and experimental 13 C/ 1 H chemical shifts can be as small as 2.10/0.18 ppm, which is much lower than those from simple DFT (5.54/0.25 ppm), or DFT + linear regression (LR) (4.77/0.23 ppm) approaches. It also has a smaller maximum absolute error than two previously proposed NMR-predicting ML models. The robustness of the DFT + ML model is tested on two classes of organic molecules (TIC10 and hyacinthacines), where the correct isomers were unambiguously assigned to the experimental ones. Overall, the DFT + ML model shows promise for structural assignments in a variety of systems, including stereoisomers, that are often challenging to determine experimentally.

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

Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction

Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO 2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. Here, we show that the sCI framework enables accurate prediction of band gaps in LiCoO 2 , closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t 2g and e g orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.

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

qFit-ligand Reveals Widespread Conformational Heterogeneity of Drug-Like Molecules in X-Ray Electron Density Maps

Here, proteins and ligands sample a conformational ensemble that governs molecular recognition, activity, and dissociation. In structure-based drug design, access to this conformational ensemble is critical to understand the balance between entropy and enthalpy in lead optimization. However, ligand conformational heterogeneity is currently severely underreported in crystal structures in the Protein Data Bank, owing in part to a lack of automated and unbiased procedures to model an ensemble of protein–ligand states into X-ray data. Here, we designed a computational method, qFit-ligand , to automatically resolve conformationally averaged ligand heterogeneity in crystal structures, and applied it to a large set of protein receptor–ligand complexes. In an analysis of the cancer related BRD4 domain, we found that up to 29% of protein crystal structures bound with drug-like molecules present evidence of unmodeled, averaged, relatively isoenergetic conformations in ligand–receptor interactions. In many retrospective cases, these alternate conformations were adventitiously exploited to guide compound design, resulting in improved potency or selectivity. Combining qFit-ligand with high-throughput screening or multitemperature crystallography could therefore augment the structure-based drug design toolbox.

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

Experimental and Theoretical Study of Oxolan-3-one Thermal Decomposition

The thermal decomposition of oxolan-3-one, a common component of the bio-oil formed during biomass pyrolysis, has been studied using ab initio calculations and experiments employing pulsed gas-phase pyrolysis with matrix isolation FTIR product detection. Four pathways for unimolecular decomposition were predicted using computational methods. The dominant reaction channel led to carbon monoxide, formaldehyde, and ethylene, all of which were observed experimentally. Here, the other channels led to an assortment of products including ketene, water, propyne, and acetylene, which were all confirmed in the matrix-isolation FTIR spectra. There is also evidence for the production of substituted ketenes in pyrolysis, most likely hydroxyketene and methylketene.

09 BIOMASS FUELS↗
DOE OSTI2020

Realistic Ion Dynamics through Charge Renormalization in Nonaqueous Electrolytes

While many practically important electrolytes contain lithium ions, interactions of these ions are particularly difficult to probe experimentally because of their small X-ray and neutron scattering cross sections and large neutron absorption cross sections. Molecular dynamics (MD) is a powerful tool for understanding the properties of nonaqueous electrolyte solutions from the atomic level, but the accuracy of this computational method crucially depends on the physics built into the classical force field. In this study, we demonstrate that several force fields for lithium bistriflimide (LiTFSI) in acetonitrile yield a solution structure that is consistent with the neutron scattering experiments, yet these models produce dramatically different ion dynamics in solution. Such glaring discrepancies indicate that inadequate representation of long-range interactions leads to excessive ionic association and ion-pair clustering. We show that reasonable agreement with the experimental observations can be achieved by renormalization of the ion charges using a “titration” method suggested herewith. This simple modification produces realistic concentration dependencies for ionic diffusion and conductivity in <2 M solutions, without loss in quality for simulation of the structure.

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

Concentration-Dependent Solvation Structure and Dynamics of Aqueous Sulfuric Acid Using Multinuclear NMR and DFT

Sulfuric acid is a ubiquitous compound for industrial processes and aqueous sulfate solutions also play a critical role in electrolytes for many prominent battery chemistries. While the thermodynamic literature for it is quite well-developed, comprehensive studies of the solvation structure, and particularly molecular-scale dynamical and transport properties, are less available. This study applies a multinuclear NMR approach to the elucidation of the solvation structure and dynamics over a wide temperature (-10 – 50?C) and concentration (0 – 18 M) range, combining 17O shift, linewidth, and T1 relaxation measurements, 33S shift and linewidth measurements, and 1H PFG-NMR measurements of the proton self-diffusivity. In conjunction, these results indicate a crossover between two regimes of solvation structure and dynamics, occurring above the concentration associated with the deep eutectic point (~ 4.5 M), with the high-concentration regime dominated by strong water-sulfate correlation. This description was borne out in detail by the activation energy trends with increasing concentration derived from the relaxation of both the H2O/H3O+ and H2SO4/HSO4-/SO42- 17O resonances and the 1H self-diffusivity. However, the 17O chemical shift difference between the H2O/H3O+ and H2SO4/HSO4-/SO42- resonances across the entire temperature range is nevertheless strikingly linear. A computational approach coupling MD simulations and DFT NMR shift calculations to reproduce this trend is presented, which will be the subject of further development. This combination of multinuclear, dynamical NMR and computational methods, and the results furnished by this study, will provide a platform for future studies on battery electrolytes where aqueous sulfate chemistry plays a central role in the solution structure.

Bazak, Jonathan D.↗
DOE OSTI2022

Structural and Electronic Effect Driven Distortions in Visible Light Absorbing Polar Materials A Ta 2 V 2 O 11 ( A = Sr, Pb)

Complex vanadates of tantalum(V), such as ATa 2 V 2 O 11 (A = Sr, Pb), are rare and underrated materials, which have potential application domains that could be substantially expanded, mitigating the existing controversy on their atomic and electronic organization. Herein, we present a thorough structural examination combining synchrotron powder X-ray diffraction-aided distortion mode analysis with computational methods to study hettotypes of SrTa 2 V 2 O 11 (STVO) and PbTa 2 V 2 O 11 (PTVO). Being distinct from the perovskite family due to the presence of [VO 4 ] groups, both compounds are polar dielectric materials with certain similarities to SBT and PBT Aurivillius phases. Applying the model of anions of metallic matrices to the analysis of electron localization functions calculated on top of as-established equilibrium structures helps retrace the effects in the Sr and Pb surroundings on the respective crystal packings of STVO and PTVO.

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

Impact of Disorder and Defects on Bond Energies in Highly Fluorinated Graphene

Fluorinated graphite materials are of interest for an assortment of applications and can be synthesized under a variety of synthetic conditions from many different types of carbon. Due to such variations, structural disorders in the form of defects and polymorphism are often present. Here, we investigate the impact of local structural variations on the C–F bond dissociation energies (BDEs) in carbon-based fluoride materials using density functional theory (DFT) computational methods. Employing fluorographene (FG) cluster models, we determine the impact of different C–F bonding configurations in the core of each platelet on the equilibrium BDEs for each C–F bond. The introduction of structural disorder decreases the first C–F BDE by approximately 1 eV compared to the canonical arrangement of axial C–F bonds ordered as in a network of cyclohexane “chairs”. Variability of calculated BDEs among the different polymorphs decreases upon subsequent F removal. Common structural tendencies of the adiabatic defluorination pathways for each polymorph are identified. In conclusion, our analysis suggests that at F/C ratios near 1.0, disorder in the local structure can play a significant role in the energetics of the initial carbon fluoride defluorination and that the influence of this configurational disorder diminishes with decreasing F/C ratios.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗