Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational molecular design”

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 325 records · Page 18

Impact of Hydrogen Bonds on CO 2 Binding in Eutectic Solvents: An Experimental and Computational Study toward Sorbent Design for CO 2 Capture

Choline-based amino acid ionic liquids with anions glycinate, β-alaninate, phenylalaninate, and prolinate were synthesized and mixed with ethylene glycol to form lower-viscosity benign eutectic solvents for CO 2 capture. The highest capacity measured was 0.7 moles of CO 2 per mole of ionic liquid (2 moles CO 2 per kg solvent) for a 1 to 2 mole ratio mixture of choline prolinate to ethylene glycol at 1 bar of CO 2 and 25 °C. Under 5000 ppm of CO 2 , half of this capacity was realized. Here, through a combined study of quantitative 13 C NMR spectroscopy, molecular dynamics simulations and density functional theory calculations, we show that hydrogen bonding in the eutectic solvent prevents proton-transfer between prolinate anions upon CO 2 absorption, which occurs in the absence of ethylene glycol and deactivates binding sites. Blocking this proton transfer leads to a higher binding capacity compared to neat choline prolinate. This work demonstrates the impact of hydrogen bonding on the CO 2 binding mechanism and energetics, as well as physical and thermal properties in eutectic solvents, thus addressing an unmet need and informing future studies on the development of benign sorbents for capturing CO 2 from dilute streams.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combining Physics and Machine Learning for the Next Generation of Molecular Simulation

Simulating molecules and atomic systems at quantum accuracy is a grand challenge for science in the 21 st century. Quantum-accurate simulations would enable the design of new medicines and the discovery of new materials. The defining problem in this challenge is that quantum calculations on large molecules, like proteins or DNA, are fundamentally impossible with current algorithms. In this work, we explore a range of different methods that aim to make large, quantum-accurate simulations possible. We show that using advanced classical models, we can accurately simulate ion channels, an important biomolecular system. We show how advanced classical models can be implemented in an exascale-ready software package. Lastly, we show how machine learning can learn the laws of quantum mechanics from data and enable quantum electronic structure calculations on thousands of atoms, a feat that is impossible for current algorithms. Altogether, this work shows that combining advances in physics models, computing, and machine learning, we are moving closer to the reality of accurately simulating our molecular world.

74 ATOMIC AND MOLECULAR PHYSICS↗

GFCCLib: Scalable and efficient coupled-cluster Green's function library for accurately tackling many-body electronic structure problems

Coupled-cluster Green’s function (GFCC) calculation has drawn much attention in the recent years for targeting the molecular and material electronic structure problems from a many-body perspective in a systematically improvable way. However, GFCC calculations on scientific computing clusters usually suffer from expensive higher di- mensional tensor contractions in the complex space, expensive inter-process communi- cation, and severe load imbalance, which limits it’s routine use for tackling electronic structure problems. Here we present a numerical library prototype that is specifically designed for large-scale GFCC calculations. The design of the library is focused on a systematically optimal computing strategy to improve its scalability and efficiency. The performance of the library is demonstrated by the relevant profiling analysis of running GFCC calculations on remote giant computing clusters. The capability of the library is highlighted by computing a wide near valence band of a fullerene C60 molecule for the first time at the GFCCSD level that shows excellent agreement with the experimental spectrum.

Peng, Bo↗

Ligand-Functionalized Polymer Membranes for Selective Ion Separations

Selective ion separations are central to technologies spanning water purification, resource recovery, and clean energy. Conventional polymer membranes, which rely on steric hindrance or Donnan exclusion, struggle to discriminate between chemically similar ions in high-ionic-strength environments. Ligand-functionalized membranes offer a transformative strategy by embedding molecular recognition directly into polymer matrices, enabling selective complexation and transport. Here, this Viewpoint highlights the structure–function relationships underlying ligand-mediated ion separation, emphasizing the interplay of dehydration penalties, ligand coordination, and nanoscale confinement. We discuss design principles, denticity, donor identity, rigidity, and spatial organization, alongside the permeability–selectivity trade-off, multicomponent effects, and stability challenges. Finally, we outline emerging strategies, from bioinspired ligands to computationally guided design, that chart a path toward next-generation membranes for precise and energy-efficient ion separations.

ions↗

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗

Assembly and Repair of the Photosystem II Reaction Center

This project investigated the biochemical and biophysical mechanisms governing assembly and repair of Photosystem II (PSII), the membrane protein complex responsible for solar-driven water oxidation in oxygenic photosynthesis. The work focused on three integrated areas: (1) protein–protein interactions that facilitate PSII assembly in cyanobacterial biogenesis centers, (2) the chemical mechanism of photo-assembly of the Mn 4 CaO 5 oxygen-evolving complex (OEC), and (3) mechanisms that target PSII reaction centers for degradation and repair in photosynthetic organisms. Using electron paramagnetic resonance spectroscopy, protein biochemistry, molecular genetics, quantitative mass spectrometry, and computational modeling, the project demonstrated that proton release events limit early steps of OEC assembly and that chloride and calcium ions facilitate Mn oxidation and intermediate stabilization. Complementary studies identified chaperone recruitment mechanisms in cyanobacterial PSII biogenesis centers and translation and protease factors involved in PSII turnover in Chlamydomonas. Together, these results establish proton management and coordinated protein quality control as central design principles in PSII assembly and repair and provide mechanistic insight relevant to biological and artificial photosynthetic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coordination-Dependent Chemical Reactivity of TFSI Anions at a Mg Metal Interface

Charge transfer across the electrode–electrolyte interface is a highly complex and convoluted process involving diverse solvated species with varying structures and compositions. Despite recent advances in in situ and operando interfacial analysis, molecular specific reactivity of solvated species is inaccessible due to a lack of precise control over the interfacial constituents and/or an unclear understanding of their spectroscopic fingerprints. However, such molecular-specific understanding is critical to the rational design of energy-efficient solid–electrolyte interphase layers. We have employed ion soft landing, a versatile and highly controlled method, to prepare well-defined interfaces assembled with selected ions, either as solvated species or as bare ions, with distinguishing molecular precision. Equipped with precise control over interfacial composition, we employed in situ multimodal spectroscopic characterization to unravel the molecular specific reactivity of Mg solvated species comprising (i.e., bis(trifluoromethanesulfonyl)imide, TFSI – ) anions and solvent molecules (i.e., dimethoxyethane, DME/G1) on a Mg metal surface relevant to multivalent Mg batteries. In situ multimodal spectroscopic characterization revealed higher reactivity of the undercoordinated solvated species [Mg-TFSI-G1] + compared to the fully coordinated [Mg-TFSI-(G1) 2 ] + species or even the bare TFSI – . These results were corroborated by the computed reaction pathways and energy barriers for decomposition of the TFSI – within Mg solvated species relative to bare TFSI – . Finally, we evaluated the TFSI reactivity under electrochemical conditions using Mg(TFSI) 2 –DME-based phase-separated electrolytes representing different solvated constituents. Based on our multimodal study, we report a detailed understanding of TFSI – decomposition processes as part of coordinated solvated species at a Mg-metal anode that will aid the rational design of improved sustainable electrochemical energy technologies.

25 ENERGY STORAGE↗

Optimal decision-making in high-throughput virtual screening pipelines

Screening large pools of molecular candidates to identify those with specific design criteria or targeted properties is demanding in various science and engineering domains. While a high-throughput virtual screening (HTVS) pipeline can provide efficient means to achieving this goal, its design and operation often rely on experts' intuition, potentially resulting in suboptimal performance. In this paper, we fill this critical gap by presenting a systematic framework that can maximize the return on computational investment (ROCI) of such HTVS campaigns. Based on various scenarios, we empirically validate the proposed framework and demonstrate its potential to accelerate scientific discoveries through optimal computational campaigns, especially in the context of virtual screening.

97 MATHEMATICS AND COMPUTING↗

Bottom-up design of actinide materials from molecular clusters: Demonstration of a general-purpose simulation capability leveraging machine-learned atomic potentials

Actinide thin-film coatings such as uranium dioxide (UO 2 ) play an important role in nuclear reactors and other mission-relevant applications, but realization of their potential requires a deep fundamental understanding of the chemical vapor deposition (CVD) processes used for their growth. The slow experimental progress can be attributed, in part, to the standard safety guidelines associated with handling uranium byproducts, which are often corrosive, toxic, and radioactive. Accurate simulation techniques, when used in concert with experiment, can improve laboratory safety, material durability, and deliverable timeframes. However, state-of-the-art computational methods are either insufficiently accurate or intractably expensive. To remedy this situation, in this project we suggested a machine-learning (ML) accelerated workflow for simulating molecular clustering toward deposition. As a benchmark test case, we considered molecular clustering in steam and assessed independent components of our workflow by comparing with measured thermodynamic properties of water. After analyzing each component individually and finding no fundamental barrier to realization of the workflow, we attempted to integrate the ML component, a Sandia-developed tool called FitSNAP. As this was the first application of FitSNAP to atoms and molecules in the gas phase at Sandia, the method required more fitting data than was originally anticipated. Systematic improvements were made by including in the fit data diatomic potentials, molecular single-bond-breaking curves, and symmetry-constrained intermolecular potentials. We concluded that our strategy provides a feasible pathway toward modeling CVD and related processes, but that extensive training data must be generated before it can be of practical use.

36 MATERIALS SCIENCE↗

Accelerating Materials Discovery: Artificial Intelligence for Sustainable, High-Performance Polymers

PolyID enables the discovery of polymers with advanced performance and greater sustainability while reducing material development timelines. The material design space is immense and cannot be reasonably probed using an Edisionian approach. High-throughput property prediction, enabled by artificial intelligence provides a hypothesis driven approach for down selection of candidate polymers to pursue experimentally. To aid experimentalists in the down selection of material targets this high-throughput, machine learning-based tool is capable of predicting polymer properties simply from molecular structures. Currently, transport, thermal, and mechanical properties across 7 polymer class (polyamides, polyesters, polycarbonates, polyimides, polyolefins, polyacrylates, and polyurethanes) can be predicted, and the PolyID platform has been flexibly designed so new materials and properties can be added.

artificial intelligence↗

The Development of Catalysts for Upgrading of Pyrolysis Vapor for Refinery Feedstocks and Intermediates (CRADA Final Report)

Catalytic fast pyrolysis (CFP) is a versatile technology platform to convert biomass into fungible hydrocarbon transportation fuels and chemical co-products. Key technical barriers to reaching this goal include increasing the product yields and achieving the desired fuel properties for gasoline, diesel, and jet range fuels or blendstocks that would be suitable for introduction into existing refinery unit operations. Overcoming these barriers will require durable catalysts that are effective at upgrading and stabilizing biomass pyrolysis vapors. Towards these goals, this CRADA leveraged NREL experience as a leader in biomass pyrolysis research and Johnson Matthey's (JM) experience as a leader in the production of advanced catalytic materials. The scope spanned CFP catalyst development, characterization, multi-scale reaction testing, and computational modeling. CRADA benefits to DOE, Participant, and U.S. Taxpayer: Assists laboratory in achieving programmatic scope, Uses the laboratory’s core competencies. The purpose of this CRADA was to develop and deploy catalysts for biomass CFP to help achieve cost-competitive biofuels and bio-based products. This was accomplished through a close collaboration between biomass conversion researchers at NREL and catalyst development researchers at JM. Summary of Research Results: Focus Area 1. Foundational research on catalytic conversion and deactivation: Key interactions between pyrolysis vapors and heterogeneous catalysts were probed through catalyst characterization, model compound reaction testing, and atomistic-scale computational modeling. Catalyst development focused on multifunctional materials, which include zeolites, oxides, carbides, and nitrides. Computational modeling identified reaction mechanisms and elucidated surface chemistry to test hypotheses regarding mechanisms of deoxygenation, coupling, cracking, dehydration, coke formation, hydrogen transfer, and aromatic ring reactions. This information was used to design multifunctional catalysts to increase product yields, control product selectivity, and reduce deactivation during CFP and downstream processing steps. The results served to increase fundamental understanding of key catalyst attributes and durability features for the upgrading of biomass pyrolysis vapors. Model compound experiments confirmed the importance of metal-acid bifunctionality for the deoxygenation of lignin-derived phenolic species under hydrodeoxygenation conditions. This insight led to the development of catalysts such as Pt/TiO2 and Mo2C, which were confirmed as high-performing materials during subsequent bench-scale experiments using biomass-derived pyrolysis vapors. This focus area also led to the identification of important catalyst deactivation mechanisms associated with the deposition of inorganic contaminants such as potassium. The molecular-level insight from model compound experiments and computational modeling, shown in Figure 1, informed the development of regeneration procedures that have been shown to be effective for restoration of > 90% of initial catalyst activity. This understanding has subsequently been translated to other catalyst systems, including zeolite materials that can be operated without requirements for co-fed hydrogen.

09 BIOMASS FUELS↗

Biochemical Process Modeling and Simulation (BPMS)

The Biochemical Process Modeling and Simulation project aims to reduce the cost and time of research by applying theory, modeling, and simulation to the most relevant bottlenecks in the biochemical process. We use molecular modeling, quantum mechanics, metabolic modeling, fluid dynamics, and reaction-diffusion methods in close collaboration with pretreatment, hydrolysis, upgrading, and TEA. The project's outcomes are increased yields and efficiency of the biochemical process, added value to products, and reduced price of fuels by specifically targeting catalytic efficiency, reactor design, enzyme efficiency, and microbial design. We work closely with experimental projects to identify problems and iterate with experiments to find and refine solutions. By working with experimentalists, we decide on problems that can be solved with simulation that could otherwise not be solved or would take too long with experiment alone to reach BETO's targets. Over the years, we have produced solutions that have resulted in determining the most likely fatty-acid derivative for passive transport out of bacteria that upgrade biomass, and we have also designed enzyme mutations for enhanced lignin upgrading. Metabolic models have been developed to tune the activity of 2,3 butanediol production for the 2030 target. A computational method to deliver understanding of how complex omics data can be interpreted in the metabolic pathways of organisms used in the Agile Biofoundry. We have found methods to overcome specific barriers and continue to develop those methods. Our reactor studies have guided the design of both the microbes and reactors for aerobic and micro-aerobic production at all scales and have been instrumental in improving the accuracy of techno-economic analysis models. This project is essential in the process of selecting the final processes for 2030 SAF production targets. More specifically, recently, we have: 1) Predicted the strength of the basic structural interactions in commodity plastics to provide guidance for plastics upcycling strategies. 2) Developed computational tool to improve the characterization of lignin-derived compounds 3) Developed new methodologies to enable Machine Learning-based Directed Evolution for protein engineering. 4) Developed Machine Learning methods to predict protein promiscuity and mutations to further improve microbial and enzymatic driven processes and demonstrated the utility of ML approaches to engineering proteins from sparse experimental datasets. 5) Developed new methods to enable high-fidelity simulation of aerobic fermentation at industrial scale and resolving mismatch of time scales through subcycling/operator splitting 7) Identified the difficulty in preventing local high-oxygen conditions in industrial bubble columns, which leads to less-desirable acetoin production, suggesting future research directions in alternative reactor configurations (e.g loop reactors, shallow-channel reactors).

BIOMASS FUELS↗

Disentangling cation effects on ion mobility and structure in ionic liquid electrolytes

Ionic liquids (ILs) are low-temperature molten salts, where ion transport is primarily governed by ion–ion interactions. Yet, the impact of organic IL cations on critical electrolyte properties such as ion dissociation and overall transport behavior in lithium-salt-doped ILs remains poorly understood. Moreover, despite their critical role in designing IL-based electrolytes for energy storage applications, ion–ion interactions and ion-specific transport under an applied electrical potential are seldom quantified, largely due to the unique experimental and computational challenges involved. Herein, we compare transport properties obtained using 1 H, 7 Li, and 19 F pulsed-field gradient nuclear magnetic resonance (NMR) and electrophoretic NMR (eNMR) with those measured by electrochemical impedance spectroscopy. Non-equilibrium molecular dynamics (MD) simulations and eNMR confirm the presence of negatively charged [Li(TFSI) n ] (1−n) aggregates that migrate towards the positive electrode, resulting in negative lithium transference numbers. Equilibrium MD simulations reveal a vehicular Li ion transport mechanism facilitated by long-lived aggregates with Li + cations strongly bound to multiple TFSI − anions. Finally, we observe an inverse relationship between the apparent charge of the TFSI − anion in the neat IL, which is dictated by the IL cation, and Li + transport in the salt-doped systems. This highlights the opportunity to tune electrolyte performance by tailoring cation chemistry.

Li-ion batteries↗

Architector 2.0: Expanded Capabilities for Metal Complex Engineering

Automated three-dimensional molecular construction from two-dimensional graph representations is critical to high-throughput discovery eIorts. Software capabilities in this area have accelerated research across fields ranging from protein design and drug discovery to transition metal catalyst development. When Architector was first introduced, it uniquely enabled high-throughput, chemically relevant three-dimensional construction of f-element complexes. Since its introduction, Architector has been applied in large-scale computational campaigns, targeted studies in critical mineral extraction, and artificial intelligence-driven discovery eIorts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HIResist: a database of HIV-1 resistance to broadly neutralizing antibodies

Changing the course of the human immunodeficiency virus type I (HIV-1) pandemic is a high public health priority with approximately 39 million people currently living with HIV-1 (PLWH) and about 1.5 million new infections annually worldwide. Broadly neutralizing antibodies (bnAbs) typically target highly conserved sites on the HIV-1 envelope glycoproteins (Envs), which mediate viral entry, and block the infection of diverse HIV-1 strains. But different mechanisms of HIV-1 resistance to bnAbs prevent robust application of bnAbs for therapeutic and preventive interventions. Here we report the development of a new database that provides data and computational tools to aid the discovery of resistant features and may assist in analysis of HIV-1 resistance to bnAbs. Bioinformatic tools allow identification of specific patterns in Env sequences of resistant strains and development of strategies to elucidate the mechanisms of HIV-1 escape; comparison of resistant and sensitive HIV-1 strains for each bnAb; identification of resistance and sensitivity signatures associated with specific bnAbs or groups of bnAbs; and visualization of antibody pairs on cross-sensitivity plots. The database has been designed with a particular focus on user-friendly and interactive interface. Our database is a valuable resource for the scientific community and provides opportunities to investigate patterns of HIV-1 resistance and to develop new approaches aimed to overcome HIV-1 resistance to bnAbs.

60 APPLIED LIFE SCIENCES↗

Electron-nuclear decoupling at a spin clock transition

The ability to design quantum systems that decouple from environmental noise sources is highly desirable for development of quantum technologies with optimal coherence. The chemical tunability of electronic states in magnetic molecules combined with advanced electron spin resonance techniques provides excellent opportunities to address this problem. Indeed, so-called clock transitions have been shown to protect molecular spin qubits from magnetic noise, giving rise to significantly enhanced coherence. Here we conduct a spectroscopic and computational investigation of this physics, focusing on the role of the nuclear bath. Away from the clock transition, linear coupling to the nuclear degrees of freedom causes a modulation and decay of electronic coherence, as quantified via electron spin echo signals generated experimentally and in silico. Meanwhile, the effective hyperfine interaction vanishes at the clock transition, resulting in electron-nuclear decoupling and an absence of quantum information leakage to the nuclear bath, providing opportunities to characterize other decoherence sources.

74 ATOMIC AND MOLECULAR PHYSICS↗

Logic-based analysis of gene expression data predicts association between TNF, TGFB1 and EGF pathways in basal-like breast cancer

For breast cancer, clinically important subtypes are well characterized at the molecular level in terms of gene expression profiles. In addition, signaling pathways in breast cancer have been extensively studied as therapeutic targets due to their roles in tumor growth and metastasis. However, it is challenging to put signaling pathways and gene expression profiles together to characterize biological mechanisms of breast cancer subtypes since many signaling events result from post-translational modifications, rather than gene expression differences. We designed a logic-based computational framework to explain the differences in gene expression profiles among breast cancer subtypes using Pathway Logic and transcriptional network information. Pathway Logic is a rewriting-logic-based formal system for modeling biological pathways including post-translational modifications. Our method demonstrated its utility by constructing subtype-specific path from key receptors (TNFR, TGFBR1 and EGFR) to key transcription factor (TF) regulators (RELA, ATF2, SMAD3 and ELK1) and identifying potential association between pathways via TFs in basal-specific paths, which could provide a novel insight on aggressive breast cancer subtypes. Lastly, codes and results are available at http://epigenomics.snu.ac.kr/PL/.

59 BASIC BIOLOGICAL SCIENCES↗

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗