Search NASA⌕ Search

SEARCH · Search NASA

Results for “in silico”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Equilibrium Separation of Siloxanes in Metal–Organic Frameworks

We present an in silico assessment of metal–organic frameworks (MOFs) for the equilibrium separation of linear and cyclic siloxanes. Using a combination of configurational bias/continuous fractional component Monte Carlo (CB/CFCMC) simulations and Ideal Adsorbed Solution Theory (IAST), we investigated the adsorption of both equimolar and nonequimolar mixtures of linear and cyclic siloxanes in a selection of synthesizable MOFs with medium to large pore volumes. We showed that configurational entropy effects drive the preferential adsorption of linear siloxanes over cyclic siloxanes. Based on synthesizability metrics we identified ZIF-70 as a promising adsorbent for the separation of linear and cyclic siloxanes. Self-diffusivities of linear and cyclic siloxanes in ZIF-70 calculated from molecular dynamics simulations show that equilibrium can be reached on reasonable time scales. We explored vacuum-temperature swing adsorption (VTSA) as a potential method for the recovery of adsorbed linear siloxanes from ZIF-70, achieving a 75% reduction in adsorbed phase concentration. Additionally, we demonstrated that supercritical CO 2 offers an alternative desorption strategy by displacing adsorbed linear siloxanes in ZIF-70 at pressures above 200 bar driven by entropy effects.

Adsorption↗
DOE OSTI2020

Full-Scale Ab Initio Simulation of Magic-Angle-Spinning Dynamic Nuclear Polarization

Theoretical models aimed at describing magic-angle-spinning (MAS) dynamic nuclear polarization (DNP) NMR have great potential in facilitating the in silico design of DNP polarizing agents and formulations. These models must typically face a trade-off between the accuracy of a strict quantum mechanical description and the need for using realistically large spin systems, for instance, using phenomenological models. In this paper, we show that the use of aggressive state-space restrictions and an optimization strategy allows full-scale ab initio MAS-DNP simulations of spin systems containing thousands of nuclei. Our simulations are shown to reproduce experimental DNP enhancements quantitatively, including their MAS rate dependence, for both frozen solutions and solid materials. They also reveal the importance of a previously unrecognized structural feature found in some polarizing agents that helps minimize the sensitivity losses imposed by the spin diffusion barrier.

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

Controlling Structure and Properties of Vapor-Deposited Glasses of Organic Semiconductors: Recent Advances and Challenges

The past decade has seen great progress in manipulating the structure of vapor-deposited glasses of organic semiconductors. Upon varying the substrate temperature during deposition, glasses with a wide range of density and molecular orientation can be prepared from a given molecule. Here, we review recent studies that show the structure of vapor-deposited glasses can be tuned to significantly improve the external quantum efficiency and lifetime of organic light-emitting diodes (OLEDs). We highlight the ability of molecular simulations to reproduce experimentally observed structures, setting the stage for in silico design of vapor-deposited glasses in the coming decade. Finally, we identify research opportunities for improving the properties of organic semiconductors by controlling the structure of vapor-deposited glasses.

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

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Efficient Multiscale Optoelectronic Prediction for Conjugated Polymers

Conjugated polymers represent a high-potential material class due to the tunability of their optoelectronic properties via straightforward processing protocols. However, to in silico tailor these optoelectronic properties, one must compute the conformationally dependent electronic structure over mesoscopic length and time scales. This task represents a challenging multiscale computational problem in which quantum chemistry and atomistic or coarse-grained molecular dynamics must be integrated. Recently, we introduced a methodology, called electronic coarse graining (ECG), that utilizes artificial neural networks to compute ab initio quality electronic structures using only the reduced degrees of freedom of coarse-grained models. Here, we adapt ECG to variable-molecular weight conjugated oligomers by casting the problem in the framework of sequence detection. A machine-learning (ML) architecture employing one-dimensional convolutional neural networks and long short-term memory networks is utilized to compute ground-state orbital energies, charge density distributions, and optical spectra solely from the coarse-grained model's configurational degrees of freedom. Robust molecular weight transferability for ECG is established via a A-ML approach that leverages model electronic Hamiltonians for ground and excited-state property determination. The accuracy and transferability of the ECG methodology opens the door for scalable optoelectronic property prediction in conjugated polymers directly from coarse-grained degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Subnanometer Visualization of Spatially Varying Local Field Resonances that Drive Tip-Enhanced Optical Spectroscopy

Our knowledge of the electromagnetic fields that power modern nanoscale optical measurements, including (non)linear tip-enhanced Raman and photoluminescence, chiefly stems from numerical simulations. Aside from idealized in silico vs heterogeneous (nano)structures in the laboratory, challenges in quantitative descriptions of nanoscale light–matter interactions more generally stem from the very nature of the problem, which lies at the interface of classical and quantum theories. This is particularly the case in ultrahigh spatial resolution measurements that are sensitive to local optical field variations that take place on subnanometer length scales. Here this work approaches this challenge through extinction-based spectral nanoimaging experiments. We demonstrate <1 nm spatial resolution in hyperspectral extinction measurements that track spatially varying plasmon resonances. We describe the principles behind our experiments and highlight more general implications of our observations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Predicting Device Parameters for Dye-Sensitized Solar Cells from Electronic Structure Calculations to Reproduce Experiment.

Given that improvements to the power-conversion efficiency (PCE) of dye-sensitized solar cells (DSSCs) have slowed in recent years, a means to accurately predict the device parameters yielded by trial dyes in silico, without having to synthesize them, would be extremely valuable to speed up the design process. Currently, the best-performing methods of calculating device parameters rely on a set of experimentally determined kinetic coefficients. In practice, it is very difficult to measure these kinetic parameters accurately, limiting the overall accuracy of such predictive methods. This work proposes a model to obtain key parameters such as J(SC), V-OC, and PCE using only the results from density functional theory (DFT) and time-dependent DFT calculations, noting that rates of electron-transfer steps are ultimately linked to the electronic structure of the dye center dot center dot center dot TiO2 working electrode. Six organic DSSC dyes from dissimilar chemical classes (L0, L1, L2, WS-2, WS-92, and C281) were chosen to demonstrate the power of this approach. Their a priori known experimentally determined device performance metrics served to validate our predictions. The greatest absolute error in our predicted PCE values was 0.36% relative to the experiment, while the greatest fractional error was 0.042. This indicates that the proposed model offers a dramatic improvement on previous predictive methods for DSSC device parameters, both in accuracy and in consistency. Moreover, such a predictive model has great potential to be applied to other photovoltaic applications, further enabling the design of novel, highly efficient photoactive materials.

density functional theory↗
DOE OSTI2023

Active Learning Guided Computational Discovery of Plant-Based Redoxmers for Organic Nonaqueous Redox Flow Batteries

Organic nonaqueous redox flow batteries (ONRFBs) are promising energy storage devices due to their scalability and reliance on sourceable materials. However, finding suitable redox-active organic molecules (redoxmers) for these batteries remains a challenge. Using plant-based compounds as precursors for these redoxmers can decrease their costs and environmental toxicity. In this computational study, flavonoid molecules have been examined as potential redoxmers for ONRFBs. Flavone and isoflavone derivatives were selected as catholyte (positive charge carrier) and anolyte (negative charge carrier) molecules, respectively. To drive their redox potentials to the opposite extremes, in silico derivatization was performed using a novel algorithm to generate a library of > 40000 candidate molecules that penalizes overly complex structures. A multiobjective Bayesian optimization based active learning algorithm was then used to identify best redoxmer candidates in these search spaces. Furthermore, our study provides methodologies for molecular design and optimization of natural scaffolds and highlights the need of incorporating expert chemistry awareness of the natural products and the basic rules of synthetic chemistry in machine learning.

25 ENERGY STORAGE↗
DOE OSTI2022

Metabolomics Analysis of Bacterial Pathogen Burkholderia thailandensis and Mammalian Host Cells in Co-culture

The Tier 1 HHS/USDA Select Agent Burkholderia pseudomallei is a bacterial pathogen that is highly virulent when introduced into the respiratory tract and intrinsically resistant to many antibiotics. Transcriptomic- and proteomic-based methodologies have been used to investigate mechanisms of virulence employed by B. pseudomallei and Burkholderia thailandensis, a convenient surrogate; however, analysis of the pathogen and host metabolomes during infection is lacking. Changes in the metabolites produced can be a result of altered gene expression and/or post-transcriptional processes. Thus, metabolomics complements transcriptomics and proteomics by providing a chemical readout of a biological phenotype, which serves as a snapshot of an organism’s physiological state. However, the poor signal from bacterial metabolites in the context of infection poses a challenge in their detection and robust annotation. In this work, we coupled mammalian cell culture-based metabolomics with feature-based molecular networking of mono- and co-cultures to annotate the pathogen’s secondary metabolome during infection of mammalian cells. These methods enabled us to identify several key secondary metabolites produced by B. thailandensis during infection of airway epithelial and macrophage cell lines. Additionally, the use of in silico approaches provided insights into shifts in host biochemical pathways relevant to defense against infection. Using chemical class enrichment analysis, for example, we identified changes in a number of host-derived compounds including immune lipids such as prostaglandins, which were detected exclusively upon pathogen challenge. Taken together, our findings indicate that co-culture of B. thailandensis with mammalian cells alters the metabolome of both pathogen and host and provides a new dimension of information for in-depth analysis of the host–pathogen interactions underlying Burkholderia infection.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Role of the Mineral in the Self-Healing of Cracks in Human Enamel

Human enamel is an incredibly resilient biological material, withstanding repeated daily stresses for decades. The mechanisms behind this resilience remain an open question, with recent studies demonstrating a crack-deflection mechanism contributing to enamel toughness and other studies detailing the roles of the organic matrix and remineralization. Here, we focus on the mineral and hypothesize that self-healing of cracks in enamel nanocrystals may be an additional mechanism acting to prevent catastrophic failure. To test this hypothesis, we used a molecular dynamics (MD) approach to compare the fracture behavior of hydroxyapatite (HAP) and calcite, the main minerals in human enamel and sea urchin teeth, respectively. In this work, we find that cracks heal under pressures typical of mastication by fusion of crystals in HAP but not in calcite, which is consistent with the resilience of HAP enamel that calcite teeth lack. Scanning transmission electron microscopy (STEM) images of structurally intact ("sound") human enamel show dashed-line nanocracks that resemble and therefore might be the cracks healed by fusion of crystals produced in silico . The fast, self-healing mechanism shown here is common in soft materials and ceramics but has not been observed in single crystalline materials at room temperature. The crack self-healing in sound enamel nanocrystals, therefore, is unique in the human body and unique in materials science, with potential applications in designing bioinspired materials.

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

Coarse-Grained Simulations for the Characterization and Optimization of Hybrid Protein–DNA Nanostructures

We present here the combination of experimental and computational modeling tools for the design and characterization of protein–DNA hybrid nanostructures. Our work incorporates several features in the design of these nanostructures: (1) modeling of the protein–DNA linker identity and length; (2) optimizing the design of protein–DNA cages to account for mechanical stresses; (3) probing the incorporation efficiency of protein–DNA conjugates into DNA nanostructures. The modeling tools were experimentally validated using structural characterization methods like cryo-TEM and AFM. Our method can be used for fitting low-resolution electron density maps when structural insights cannot be deciphered from experiments, as well as enable in-silico validation of nanostructured systems before their experimental realization. Finally, these tools will facilitate the design of complex hybrid protein–DNA nanostructures that seamlessly integrate the two different biomolecules.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Extraction of Furfural and Furfural/5-Hydroxymethylfurfural from Mixed Lignocellulosic Biomass-Derived Feedstocks

We present a combined in silico and experimental study on the extraction of furfural and 5-hydroxymethylfurfural (HMF) in aqueous–organic biphasic systems. We predict the liquid–liquid equilibria and furfural/HMF partition coefficients of over 2200 water-organic biphasic systems using the multiscale COSMO-RS model and measure experimentally single-component (furfural) and mixture (furfural and HMF) partition coefficients at room and dehydration reaction-relevant temperatures in 28 solvents. We find the experimental data to be within a factor of 2 from the COSMO-RS predictions. Even though furfural and HMF have chemical similarity, the slight differences in molecular structure render the separation of furfural easier by the supply of more solvents of higher partition coefficient for extraction. We leverage this molecular difference and experimentally demonstrate that with an additional extraction step, using dichloromethane or toluene, we can selectively extract furfural from furfural-HMF mixtures, which can coexist in lignocellulosic biomass dehydration products, despite their partition coefficients being generally correlated. We complement solvent selection criteria for biphasic lignocellulosic biomass processes with a simple mass balance extraction model for determining volume ratios in multistage extraction. Lastly, the molecular nature of the preferential furfural extraction is rationalized using COSMO-RS σ-profile analysis.

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

Mapping Composition Evolution through Synthesis, Purification, and Depolymerization of Random Heteropolymers

Random heteropolymers (RHPs) consisting of three or more comonomers have been routinely used to synthesize functional materials. While increasing the monomer variety diversifies the side-chain chemistry, this substantially expands the sequence space and leads to ensemble-level sequence heterogeneity. Most studies have relied on monomer composition and simulated sequences to design RHPs, but the questions remain unanswered regarding heterogeneities within each RHP ensemble and how closely these simulated sequences reflect the experimental outcomes. Here, we quantitatively mapped out the evolution of monomer compositions in four-monomer-based RHPs throughout a design-synthesis-purification-depolymerization process. By adopting a Jaacks method, we first determined 12 reactivity ratios directly from quaternary methacrylate RAFT copolymerization experiments to account for the influences of competitive monomer addition and the reversible activation/deactivation equilibria. The reliability of in silico analysis was affirmed by a quantitative agreement (<4% difference) between the simulated RHP compositions and the experimental results. Furthermore, we mapped out the conformation distribution within each ensemble in different solvents as a function of monomer chemistry, composition, and segmental characteristics via high-throughput computation based on self-consistent field theory (SCFT). These comprehensive studies confirmed monomer composition as a viable design parameter to engineer RHP-based functional materials as long as the reactivity ratios are accurately determined and the livingness of RHP synthesis is ensured.

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

The Factors Governing Metal Dependence of an Emergent Superfamily of Bimetallic Oxygenases

Metalloenzyme superfamilies are typically defined by their protein scaffolds and active sites. Owing to the high tunability of protein structures, members of a single superfamily can catalyze diverse reactions with the same metallocofactor. Some superfamilies, such as amidohydrolase-related dinuclear oxygenases (AROs), display further versatility by utilizing multiple metallocofactors. We have shown that certain AROs catalyze monooxygenation reactions with diiron, dimanganese, and/or mixed manganese−iron cofactors, but the molecular factors governing the selection of a particular cofactor remain unknown, and the extent of this superfamily in biology is unclear. Here, we report bioinformatic analyses that expand the ARO superfamily to approximately 17,000 unique UniProt sequences, far exceeding the number of previously characterized enzymes. Through the integration of structural, spectroscopic, and thermodynamic analyses of representative proteins with a bioinformatic pipeline that identifies key secondary- and tertiary-sphere residues, we can predict in silico the metal preference for the majority of reported ARO sequences. These annotations were validated via the characterization of multiple new AROs, including ones implicated in key oxidative steps of natural product biosyntheses. This study establishes the key structure−function relationships governing metal preferences in AROs and highlights their vastly underappreciated role in myriad biological processes.

Liu, Chang [University of California, Berkeley, CA↗
DOE OSTI2024

Advancing the Prediction of MS/MS Spectra Using Machine Learning

Tandem mass spectrometry (MS/MS) is an important tool for the identification of small molecules and metabolites where resultant spectra are most commonly identified by matching them with spectra in MS/MS reference libraries. While popular, this strategy is limited by the contents of existing reference libraries. In response to this limitation, various methods are being developed for the in silico generation of spectra to augment existing libraries. Recently, machine learning and deep learning techniques have been applied to predict spectra with greater speed and accuracy. Here, in this work, we investigate the challenges these algorithms face in achieving fast and accurate predictions on a wide range of small molecules. The challenges are often amplified by the use of generic machine learning benchmarking tactics, which lead to misleading accuracy scores. Curating data sets, only predicting spectra for sufficiently high collision energies, and working more closely with experimental mass spectrometrists are recommended strategies to improve overall prediction accuracy in this nuanced field.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Host population diversity as a driver of viral infection cycle in wild populations of green sulfur bacteria with long standing virus-host interactions

Abstract Temperate phages are viruses of bacteria that can establish two types of infection: a lysogenic infection in which the virus replicates with the host cell without producing virions, and a lytic infection where the host cell is eventually destroyed, and new virions are released. While both lytic and lysogenic infections are routinely observed in the environment, the ecological and evolutionary processes regulating these viral dynamics are still not well understood, especially for uncultivated virus-host pairs. Here, we characterized the long-term dynamics of uncultivated viruses infecting green sulfur bacteria (GSB) in a model freshwater lake (Trout Bog Lake, TBL). As no GSB virus has been formally described yet, we first used two complementary approaches to identify new GSB viruses from TBL; one in vitro based on flow cytometry cell sorting, the other in silico based on CRISPR spacer sequences. We then took advantage of existing TBL metagenomes covering the 2005–2018 period to examine the interactions between GSB and their viruses across years and seasons. From our data, GSB populations in TBL were constantly associated with at least 2-8 viruses each, including both lytic and temperate phages. The dominant GSB population in particular was consistently associated with two prophages with a nearly 100% infection rate for >10 years. We illustrate with a theoretical model that such an interaction can be stable given a low, but persistent, level of prophage induction in low-diversity host populations. Overall, our data suggest that lytic and lysogenic viruses can readily co-infect the same host population, and that host strain-level diversity might be an important factor controlling virus-host dynamics including lytic/lysogeny switch.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

PplD is a de-N-acetylase of the cell wall linkage unit of streptococcal rhamnopolysaccharides

The cell wall of the human bacterial pathogen Group A Streptococcus (GAS) consists of peptidoglycan decorated with the Lancefield group A carbohydrate (GAC). GAC is a promising target for the development of GAS vaccines. In this study, employing chemical, compositional, and NMR methods, we show that GAC is attached to peptidoglycan via glucosamine 1-phosphate. This structural feature makes the GAC-peptidoglycan linkage highly sensitive to cleavage by nitrous acid and resistant to mild acid conditions. Using this characteristic of the GAS cell wall, we identify PplD as a protein required for deacetylation of linkage N-acetylglucosamine (GlcNAc). X-ray structural analysis indicates that PplD performs catalysis via a modified acid/base mechanism. Genetic surveys in silico together with functional analysis indicate that PplD homologs deacetylate the polysaccharide linkage in many streptococcal species. We further demonstrate that introduction of positive charges to the cell wall by GlcNAc deacetylation protects GAS against host cationic antimicrobial proteins.

59 BASIC BIOLOGICAL SCIENCES↗