Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer simulation”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Rapid Droplet Sampling Interface for Low-Volume, High-Throughput Mass Spectrometry Analysis

Here, we present a rapid droplet sampling interface (RDSI) electrospray ionization mass spectrometry (ESI-MS) system as a high-throughput, low-volume, noncontact, and minimal-carryover approach for characterization of liquids. Liquid characterization was achieved by combining droplet ejection with an open-face microflow capillary with a 2.5 μL/min continuous flow of carrier solvent. Through this implementation, single 0.3 nL droplets containing the analyte effectively mix with 4–8 nL of carrier solvent and create a combined electrospray plume. The carrier solvent continuously cleaned the system, eliminating carryover. A sampling rate of 5 Hz was achieved for droplets containing 1 μM propranolol or 5 μM leu-enkephalin with each droplet fully baseline-resolved (138 ± 32 ms baseline-to-baseline). Using a SCIEX API4000 mass spectrometer, a lower limit of quantification (LLOQ) of propranolol was 15 nM, corresponding to 1.16 fg of propranolol in the droplet, and was linear across 3 orders of magnitude. Further, quantitation could be achieved by adding an isotopically labeled internal standard, as done in conventional ESI. Signal transients were faster than the acquisition speed of the mass spectrometer, resulting in artificially high reproducibility of 15–30% RSD droplet-to-droplet. Analyte–solvent mixing ratios could be controlled by adjusting droplet positioning along the open-face capillary with an optimal position about 0.4 mm from the tip end. The range of analyte coverage was exemplified by measures of peptides and drugs in methanol, water, and buffer solutions. In a comparison to the Open Port Sampling Interface (OPSI) implemented on the same system, the RDSI had 78× greater sensitivity, 6× greater throughput and used significantly less carrier solvent.

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

Risk Assessment in a Chemical Laboratory Following an Explosive Incident Involving a Novel Diazonium Compound: Retrospective Analysis and Lessons Learned

Diazonium compounds are synthetically useful in the production of dyes and textiles, however they are highly explosive under dry conditions. Explosion prevention becomes more difficult when new diazonium compounds are synthesized, because while some syntheses include a counterion to increase their stability, this is not always a reliable method to prevent an explosive incident. Due to the uncertainty surrounding the explosiveness of different diazonium compounds, it is important to understand how to safely clean up after an incident and how to determine when it is safe to return a laboratory to typical operational use, particularly when the incident involves a novel compound where a standard does not exist for instrument calibration. Here, an explosive event is discussed involving the synthesis of 4-bromo-benzenediazonium-2-carboxylate. Following the explosive incident and 3-step cleanup, which involved a precautionary neutralization step, samples were collected from the fume hood where the incident occurred. Because the incident involved an unstable, novel compound that is not commercially available and was deemed unsafe to resynthesize for instrument calibration, we assessed the risk of further explosion by analyzing for the stable decomposition products. Mass spectrometry analysis confirmed that the residue in the fume hood contained 5-bromosalicylic acid, a decomposition product of 4-bromo-benzenediazonium-2-carboxylate. Samples were taken from multiple points in the fume hood and analyzed to estimate the spatial distribution of the decomposition product. Based on this analysis, we inferred that the primary decomposition product was far more abundant than residual energetic, indicating the energetic had been consumed or neutralized to a trace quantity where the risk of further explosion was low. Furthermore, the steps presented here─specifically, initial neutralization and then analyzing the spatial distribution of expected decomposition products to assess risk when a novel explosive material is detonated in a confined space─were our approach to assess further risk following an explosion due to a novel diazonium compound without the need for any further handling or resynthesis of the energetic. Here, we present our approach and critically analyze these steps by discussing retrospective lessons learned and alternative analytical approaches.

Computer simulations↗
DOE OSTI2024

Data Generation for Machine Learning Interatomic Potentials and Beyond

The field of data-driven chemistry is undergoing an evolution, driven by innovations in machine learning models for predicting molecular properties and behavior. Recent strides in ML-based interatomic potentials have paved the way for accurate modeling of diverse chemical and structural properties at the atomic level. The key determinant defining MLIP reliability remains the quality of the training data. A paramount challenge lies in constructing training sets that capture specific domains in the vast chemical and structural space. This Review navigates the intricate landscape of essential components and integrity of training data that ensure the extensibility and transferability of the resulting models. We delve into the details of active learning, discussing its various facets and implementations. We outline different types of uncertainty quantification applied to atomistic data acquisition and the correlations between estimated uncertainty and true error. The role of atomistic data samplers in generating diverse and informative structures is highlighted. Furthermore, we discuss data acquisition via modified and surrogate potential energy surfaces as an innovative approach to diversify training data. The Review also provides a list of publicly available data sets that cover essential domains of chemical space.

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

Integrating Contaminant Source Indicators, Water Quality Measures, and Ecotoxicity to Characterize Contaminant Mixtures and Per- and Polyfluoroalkyl Substance (PFAS) Variability in an Urban Watershed

Thousands of chemical contaminants threaten watersheds but are time and cost prohibitive to monitor. Identifying their sources, transport, and ecological risk is limited in heterogeneous urban watersheds. We present an integrative watershed approach using source-specific indicator compounds, common water quality measures, and ecotoxicity assays to examine the distribution of contaminant mixtures in an urbanized watershed. Indicator compound concentrations were temporally and spatially distributed for treated/untreated sewage (sucralose, artificial sweetener), road runoff (diphenyl-guanidine [DPG] and 6PPD-quinone [6PPD-Q], automobile tire additives), and lawncare runoff (aminomethanephosphonic acid (AMPA), major degradant of the herbicide glyphosate). Sucralose was predominately sourced from treated wastewater; measurable concentrations in tributaries indicated raw sewage inputs. DPG and 6PPD-Q concentrations correlated to road density during base flow and were elevated during stormflow. AMPA was measurable spring through fall, especially where lawns were dense. When specific sources dominated flow, water quality measures correlated with wastewater (sulfate, potassium, chloride, and sodium) and road runoff (chromium and lead) indicators. The limited behavioral toxicity observed in exposed zebrafish (Danio rerio) (18%) was not well explained by source-indicators. PFAS concentrations were highly variable spatially but not well explained by our source-specific indicator compounds. Here, more costly compound-specific monitoring may be necessary when multiple sources exist or when unexpected toxicity trends occur.

computer simulations↗
DOE OSTI2025

Multiple Sources of Riparian Wetland Suspended Solids during Episodic Rain Events: Influence on Uranium Transport

Suspended solids can be the primary vector for transporting contaminants in streams. The objective of this study was to determine whether changes in the properties of suspended solids during rain events impacted contaminant transport. Stream water was collected during five episodic events downstream from a U-contaminated wetland located in South Carolina, USA. The suspended particles were initially composed of Fe-flocs (particles formed in situ prior to the rain event) that had significantly greater Fe, Mn, organic-C, and U content than particles collected later during a sampling rain event. XANES and EXAFS revealed that U in the Fe-flocs was U(VI) and that it was not incorporated in a mineral structure but existed as inner- or outer-sphere adsorbed uranyl species associated with organic matter and Fe-oxides. The uranyl had an extraordinarily high affinity for the suspended solids, with solid to liquid U ratios of >72,000 (μg/kg)/(μg/L). After the initial flush of Fe-flocs, a greater fraction of the suspended solids had lower organic-C, Fe, Mn, and amorphous phases and were composed of more quartz, kaolinite, and gibbsite, resulting in lower U concentrations than those in the solids collected earlier in the rain event. This study highlights the importance of understanding suspended solids as transport vectors and their potential dynamic nature during rain events.

computer simulations↗
DOE OSTI2020

Consistency-Enhanced Evolution for Variable Selection Can Identify Key Chemical Information from Spectroscopic Data

In the last few decades, spectroscopic techniques such as near-infrared (NIR) spectroscopy have gained wide applications in several industries, such as the pharmaceutical, agricultural, oil, and gas industries. As a result, various soft sensors have been developed to predict sample properties from spectroscopic readings. Because the spectroscopic readings at different wavelengths, especially at the adjacent wavelengths, are highly correlated, it has been shown that variable selection could significantly improve a soft sensor’s prediction performance while reducing the model complexity. To improve the prediction performance, most variable selection methods focus on identifying the variables (i.e., wavelengths or wavelength segments) that are strongly correlated with the dependent variable. Although many successful applications have been reported, these variable selection methods do have their limitations. Specifically, the selected wavelengths sometimes show little connection to the chemical bounds or functional groups presenting in the sample. In addition, the selected variables can be quite sensitive to the choice of the training samples. In this work, we address these limitations from a different perspective: if a variable selection algorithm can identify the truly relevant input variables, it should consistently identify the same subset of variables regardless of the choice of the training samples. Therefore, we propose a variable selection method that aims to improve the consistency of variable selection resulting from different training samples. Furthermore, the new algorithm is termed consistency-enhanced evolution for variable selection (CEEVS). To demonstrate the performance and robustness of CEEVS, we compare the proposed method with three representative variable selection methods using five published NIR data sets. These case studies clearly demonstrate that by improving the variable selection consistency, we can not only achieve improved prediction performance, but also identify key chemical information from spectroscopic data.

42 ENGINEERING↗
DOE OSTI2023

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge or artificial intelligence driven, and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary- structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models, which contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly, while refining them against the density maps. Here, we introduce such an adaptive decision making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, and the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two protein systems, Adenylate Kinase and Carbon Monoxide Dehydrogenase. For the test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2 - 3 Å, map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that is inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from Chimpanzee adenovirus. The overhead for decision making remaining low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install the R-MDFF on different computing environments, from local Linux based workstations to High Performance Computing (HPC) environments.

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

QuantumScents: Quantum-Mechanical Properties for 3.5k Olfactory Molecules

Quantitative structure–odor relationships are critically important for studies related to the function of olfaction. Current literature data sets contain expert-labeled molecules but lack feature data. This paper introduces QuantumScents, a quantum mechanics augmented derivative of the Leffingwell data set. QuantumScents contains 3.5k structurally and chemically diverse molecules ranging from 2 to 30 heavy atoms (CNOS) and their corresponding 3D coordinates, total PBE0 energy, molecular dipole moment, and per-atom Hirshfeld charges, dipoles, and ratios. The authors demonstrate that Hirshfeld charges and ratios contain sufficient information to perform molecular classification by training a Message Passing Neural Network with chemprop (Heid, E.; et al. ChemRxiv, 2023, DOI: 10.26434/chemrxiv-2023-3zcfl) to predict scent labels. Finally, the QuantumScents data set is freely available on Zenodo along with the authors’ code, example models, and data set generation workflow (https://zenodo.org/doi/10.5281/zenodo.8239853).

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

Rapid Adaptation of Chemical Named Entity Recognition Using Few-Shot Learning and LLM Distillation

Named entity recognition (NER) has been widely used in chemical text mining for the automatic identification and extraction of chemical entities. However, existing chemical NER systems primarily focus on scenarios with abundant training data, requiring significant human effort on annotations. This poses challenges for applications in the chemical field, such as catalysis, where many advancements have traditionally relied on trial-and-error investigations and incremental adjustment of variables. This hinders catalysis science and technology progress in addressing emerging energy and environmental crises. In this work, we propose a few-shot NER model that can quickly adapt to extract new types of chemical entities by using only a limited number of annotated examples. Our model employs a metric-learning approach to transfer entity similarity knowledge from high-resource chemical domains (with abundant annotations) to enable effective entity recognition in low-resource specialized domains (limited annotation). We validate the effectiveness of our model on a few-shot chemical NER benchmark built based on six existing chemical NER data sets. Experiments show that the proposed few-shot NER model can achieve reasonable performance with only 5 examples per entity type and shows consistent improvement as the number of examples increases. Furthermore, we demonstrate how the proposed model can be trained with large language model (LLM) annotated data, opening a new pathway for rapid adaptation of NER systems. Furthermore, our approach leverages the knowledge broadness of large language models for chemistry while distilling this knowledge into a lightweight model suitable for efficient and in-house use.

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

Multitask Machine Learning of Collective Variables for Enhanced Sampling of Rare Events

Computing accurate reaction rates is a central challenge in computational chemistry and biology because of the high cost of free energy estimation with unbiased molecular dynamics. In this work, a data-driven machine learning algorithm is devised to learn collective variables with a multitask neural network, where a common upstream part reduces the high dimensionality of atomic configurations to a low dimensional latent space and separate downstream parts map the latent space to predictions of basin class labels and potential energies. Here, the resulting latent space is shown to be an effective low-dimensional representation, capturing the reaction progress and guiding effective umbrella sampling to obtain accurate free energy landscapes. This approach is successfully applied to model systems including a 5D Müller Brown model, a 5D three-well model, the alanine dipeptide in vacuum, and an Au(110) surface reconstruction unit reaction. It enables automated dimensionality reduction for energy controlled reactions in complex systems, offers a unified and data-efficient framework that can be trained with limited data, and outperforms single-task learning approaches, including autoencoders.

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

Coarse-Grained Density Functional Theory Predictions via Deep Kernel Learning

Scalable electronic predictions are critical for soft materials design. Recently, the Electronic Coarse-Graining (ECG) method was introduced to renormalize all-atom quantum chemical (QC) predictions to coarse-grained (CG) resolutions using deep neural networks (DNNs). While DNNs can learn complex representations that prove challenging for kernel-based methods, they are susceptible to overfitting and the overconfidence of uncertainty estimations. Here, we develop ECG within a GPU-accelerated Deep Kernel Learning (DKL) framework to enable CG QC predictions using range-separated hybrid density functional theory (DFT), obtaining a 107 speedup relative to naive all-atom QC. By treating the predicted electronic properties as random Gaussian Processes, DKL incorporates CG mapping degeneracy by learning the distribution of electronic energies as a function of CG configuration. DKL-ECG accurately reproduces molecular orbital energies from range-separated DFT while facilitating efficient training via active learning using the uncertainties provided by DKL. Further, we show that while active learning algorithms enable efficient sampling of a more diverse configurational space relative to random sampling, all explored query methods exhibit comparable performance for the examined system. We attribute this result to the significant overlap of the feature space and output property distributions across multiple temperatures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Comparison of Equilibrium and Nonequilibrium Approaches for Relative Binding Free Energy Predictions

Alchemical relative binding free energy calculations have recently found important applications in drug optimization. A series of congeneric compounds are generated from a preidentified lead compound, and their relative binding affinities to a protein are assessed in order to optimize candidate drugs. While methods based on equilibrium thermodynamics have been extensively studied, an approach based on nonequilibrium methods has recently been reported together with claims of its superiority. However, these claims pay insufficient attention to the basis and reliability of both methods. Here we report a comparative study of the two approaches across a large data set, comprising more than 500 ligand transformations spanning in excess of 300 ligands binding to a set of 14 diverse protein targets. Ensemble methods are essential to quantify the uncertainty in these calculations, not only for the reasons already established in the equilibrium approach but also to ensure that the nonequilibrium calculations reside within their domain of validity. If and only if ensemble methods are applied, we find that the nonequilibrium method can achieve accuracy and precision comparable to those of the equilibrium approach. Compared to the equilibrium method, the nonequilibrium approach can reduce computational costs but introduces higher computational complexity and longer wall clock times. There are, however, cases where the standard length of a nonequilibrium transition is not sufficient, necessitating a complete rerun of the entire set of transitions. This significantly increases the computational cost and proves to be highly inconvenient during large-scale applications. Our findings provide a key set of recommendations that should be adopted for the reliable implementation of nonequilibrium approaches to relative binding free energy calculations in ligand-protein systems.

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

Permutationally Invariant Networks for Enhanced Sampling (PINES): Discovery of Multimolecular and Solvent-Inclusive Collective Variables

The typically rugged nature of molecular free energy landscapes can frustrate efficient sampling of the thermodynamically relevant phase space due to the presence of high free energy barriers. Enhanced sampling techniques can improve phase space exploration by accelerating sampling along particular collective variables (CVs). A number of techniques exist for data-driven discovery of CVs parameterizing the important large scale motions of the system. A challenge to CV discovery is learning CVs invariant to symmetries of the molecular system, frequently rigid translation, rigid rotation, and permutational relabeling of identical particles. Of these, permutational invariance has proved a persistent challenge in frustrating the data-driven discovery of multi-molecular CVs in systems of self-assembling particles and solvent-inclusive CVs for solvated systems. In this work, we integrate Permutation Invariant Vector (PIV) featurizations with autoencoding neural networks to learn nonlinear CVs invariant to translation, rotation, and permutation, and perform interleaved rounds of CV discovery and enhanced sampling to iteratively expand sampling of configurational phase space and obtain converged CVs and free energy landscapes. Here, we demonstrate the Permutationally Invariant Network for Enhanced Sampling (PINES) approach in applications to the self-assembly of a 13-atom Argon cluster, association/dissociation of a NaCl ion pair in water, and hydrophobic collapse of a C 45 H 92 n-pentatetracontane polymer chain. We make the approach freely available as a new module within the PLUMED2 enhanced sampling libraries.

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

Efficient 3D Molecular Design with an E(3) Invariant Transformer VAE

This work introduces a three-dimensional (3D) invariant graph-to-string transformer variational autoencoders (VAE) (Vagrant) for generating molecules with accurate density functional theory (DFT)-level properties. Vagrant learns to model the joint probability distribution of a 3D molecular structure and its properties by encoding molecular structures into a 3D-aware latent space. Directed navigation through this latent space implicitly optimizes the 3D structure of a molecule, and the latent embedding can be used to condition a generative transformer to predict the candidate structure as a one-dimensional (1D) sequence. Additionally, we introduce two novel sampling methods that exploit the latent characteristics of a VAE to improve performance. Here, we show that our method outperforms comparable 3D autoregressive and diffusion methods for predicting quantum chemical property values of novel molecules in terms of both sample quality and computational efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗