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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Oxidation Dynamics of Supported Catalytic Cu Clusters: Coupling to Fluxionality

Copper oxide nanoclusters have a wide range of catalytic applications, such as the selective oxidation of hydrocarbons. O 2 binding to the catalyst, activation, and release upon reagent oxidation are key events in these catalytic chemistries. Furthermore, these events are expected to be accompanied by significant structural changes of the Cu clusters, because O atoms integrate into the cluster, rather than bind to its surface. Topping the complexity of the problem, partially oxidized Cu clusters are known to exhibit strong fluxionality and feature diverse and interconverting structures and oxygen contents in conditions of oxidative dehydrogenation (ODH). Hence, a significant dynamic coupling between the “hot” O 2 molecule impacting the cluster at reaction temperatures and the cluster fluxionality can be expected. In this work, we focus on the dynamics of dioxygen integration into a partially oxidized Cu cluster supported on hydroxylated amorphous alumina–a system recently reported to be an exceptionally selective catalyst for cyclohexane ODH with very little CO produced, whose mechanistic underpinnings are of utmost interest. The statistics over a swarm of adsorption and scattering trajectories where O 2 hits various sites on the cluster at reaction temperature shows that the O 2 binding does not only follow the minimal energy paths. O 2 also rarely integrates into the cluster in a single step and instead first binds to a single Cu atom via either an η 1 -O 2 or an η 2 -O 2 mode. Surprisingly, this step often has a higher barrier than the subsequent O 2 integration and dissociation, which in turn take multiple steps and complete the oxidation process. Dynamic trajectories starting from the key transition state of integration of the adsorbed O 2 can also lead to different intermediate structures during or right after the dissociation, due to the energy released from the transition state and the thermal intracluster effects. From these activated O 2 chemisorbed structures, O 2 dissociation occurs with moderate barriers (~0.5 eV), producing multiple final oxidized Cu 4 O 4 states. Hence, a diversity of reaction profiles for the attack of supported Cu cluster by O 2 emerges due to the dynamic effects, with implications for mechanisms, kinetic models, and catalyst design principles.

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

Cabana: A Performance Portable Library for Particle-Based Simulations

Particle-based simulations are ubiquitous throughout many fields of computational science and engineering, spanning the atomistic level with molecular dynamics (MD), to mesoscale particle-in-cell (PIC) simulations for solid mechanics, device-scale modeling with PIC methods for plasma physics, and massive N-body cosmology simulations of galaxy structures, with many other methods in between (Hockney & Eastwood, 1989). While these methods use particles to represent significantly different entities with completely different physical models, many low-level details are shared including performant algorithms for short- and/or long-range particle interactions, multi-node particle communication patterns, and other data management tasks such as particle sorting and neighbor list construction. Cabana is a performance portable library for particle-based simulations, developed as part of the Co-Design Center for Particle Applications (CoPA) within the Exascale Computing Project (ECP) (Alexander et al., 2020). The CoPA project and its full development scope, including ECP partner applications, algorithm development, and similar software libraries for quantum MD, is described in (Mniszewski et al., 2021). Cabana uses the Kokkos library for on-node parallelism (Edwards et al., 2014; Trott et al., 2022), enabling simulation on multi-core CPU and GPU architectures, and MPI for GPU-aware, multi-node communication. Cabana provides particle simulation capabilities on almost all current Kokkos backends, including serial execution, OpenMP (including OpenMP-Target for GPUs), CUDA (NVIDIA GPUs), HIP (AMD GPUs), and SYCL (Intel GPUs), providing a clear path for the coming generation of accelerator-based exascale hardware. Cabana builds on Kokkos by providing new particle data structures and particle algorithms resulting in a similar execution policy-based, node-level programming model that is intended to be used in addition to the core Kokkos library within an application. Cabana is designed as an application and physics agnostic, but particle-specific toolkit which can either be used to generate a new application, or to be used as needed in existing applications at various levels of invasiveness including through interfaces that wrap user memory in existing data structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

First-Principles Grand-Canonical Simulations of Water Adsorption in Proton-Exchanged Zeolites

Water appears by design or as impurities in many important reactive systems. For those catalyzed by porous solid acids, such as widely used zeolites, experimentally quantifying the amount or elucidating the structure of adsorbed water clusters at reaction conditions is challenging, while computational studies (e.g., first-principles molecular dynamics simulations) often need to assume the loading to examine solvation effects. Furthermore, we perform first-principles grand-canonical simulations to predict water adsorption to H-ZSM-5 zeolites under specified experimental conditions. Presampling with inexpensive force fields and a pool-based parallelization algorithm are used to improve simulation efficiency, while molecular dynamics is used to sample configurations involving hydronium species. We observe that H + exchange dramatically increases the hydrophilicity of zeolite MFI and an appreciable amount of water is present at very low relative humidities. At all conditions examined, the zeolitic protons are found to dissociate readily from the surface basic sites and become mobile by participating in the hydrogen-bonded chains of adsorbed water molecules.

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

Integrating Systems and Synthetic Biology to Understand and Engineer Microbiomes

Microbiomes are complex and ubiquitous networks of microorganisms whose seemingly limitless chemical transformations could be harnessed to benefit agriculture, medicine, and biotechnology. Here, the spatial and temporal changes in microbiome composition and function are influenced by a multitude of molecular and ecological factors. This complexity yields both versatility and challenges in designing synthetic microbiomes and perturbing natural microbiomes in controlled, predictable ways. In this review, we describe factors that give rise to emergent spatial and temporal microbiome properties and the meta-omics and computational modeling tools that can be used to understand microbiomes at the cellular and system levels. We also describe strategies for designing and engineering microbiomes to enhance or build novel functions. Throughout the review,we discuss key knowledge and technology gaps for elucidating the networks and deciphering key control points for microbiome engineering, and highlight examples where multiple omics and modeling approaches can be integrated to address these gaps.

42 ENGINEERING↗
DOE OSTI2025

Probing Condensed-Phase Structure and Dynamics in Hierarchical Zeolites and Nanosheets for Catalytic Upgradation of Biomass (Final Report)

Understanding complex reaction pathways in systems governed by multi-scale collective interactions across time and length scales remains a central scientific challenge. This project was guided by the hypothesis that the interplay among oligomers, solvents, and active sites can be tuned by a suitable choice of solvation environment and pore architecture in solid-acid catalysts to direct chemical transformations relevant to biomass conversion. Zeolites and zeolite nanosheets were used as model platforms, allowing for the interaction of macromolecules with the surface of the zeolite nanosheets and with smaller pores that host catalytically active sites. To investigate these coupled phenomena, we employ a multi-scale computational framework that integrates molecular-level descriptions with advanced sampling approaches to capture key physical and chemical interactions. Our work through this project improved fundamental understanding of how reactants and solid-acid catalysts interact in solvent-rich environments, thereby enabling the rational design of catalytic systems that upgrade biomass with enhanced selectivity and energy efficiency. In addition, the project developed advanced sampling methodologies critical for disentangling complex, reactive processes in multi-component catalytic environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

Algorithms for Haptic Rendering of 3D Objects

Algorithms have been developed to provide haptic rendering of three-dimensional (3D) objects in virtual (that is, computationally simulated) environments. The goal of haptic rendering is to generate tactual displays of the shapes, hardnesses, surface textures, and frictional properties of 3D objects in real time. Haptic rendering is a major element of the emerging field of computer haptics, which invites comparison with computer graphics. We have already seen various applications of computer haptics in the areas of medicine (surgical simulation, telemedicine, haptic user interfaces for blind people, and rehabilitation of patients with neurological disorders), entertainment (3D painting, character animation, morphing, and sculpting), mechanical design (path planning and assembly sequencing), and scientific visualization (geophysical data analysis and molecular manipulation).

Basdogan, Cagatay↗
DOE OSTI2025

Evaluating the Structural Response of Amphiphilic Monolayers to Environmental Stimuli

Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.

functional groups↗
DOE OSTI2026

Noncanonical folding of peptoid oligomers: Formation of a closed conformation in nonpolar solvent

Peptoids provide a versatile platform for foldamer design, yet their conformational behavior in low-dielectric media remains poorly understood. The structural characteristics of N-1-phenylethylglycine (Nspe) homo-oligomers were inves-tigated in chloroform, a solvent that mimics the interior of lipid bilayers, using nuclear magnetic resonance (NMR) spec-troscopy and molecular dynamics (MD) simulations. Nspe7 populated two distinct compact conformations, while Nspe10 adopted a single homogeneous conformation related to the previously reported Nspe9 threaded-loop structure. Integrated experimental and computational analysis reveals that these structures are stabilized by cooperative end-to-end intramo-lecular hydrogen bonding, cis-trans backbone isomerism, and hydrophobic side-chain shielding. The resulting structures minimize exposed polar surface area, demonstrating a closed conformation in the low-dielectric environment. These find-ings establish specific chain-length requirements for achieving well-defined closed conformations. This work provides insights into peptoid folding in nonpolar media, enabling rational design strategies for solvent-directed conformational switching systems.

Oh, Jinyoung↗
DOE OSTI2025

Dynamic separation of gases using microsieves

Separation of light weight molecules, such as nitrogen, argon, and oxygen, from heavier compounds can have significant impacts on energy capture, environmental monitoring, or isotopic applications. Large-scale gas separation techniques, like gas centrifugation and membrane mitigation, can be problematic as they impart tremendous energy and induce high mechanical stress onto the instrumentation. Microsieves, also known as micronozzles or microfunnels, are developed to create physical barriers to separate specific isotopes and gases. Separation is achieved using a converging and diverging micronozzle to impose supersonic gas flow around a curved wall, and it has been used for the separation of heavy actinide isotopes in low weight gas as well as separation of low weight gas compositions of nitrogen and argon back in 1900s. However, systematic reviews of this unique technology are lacking. The application of the Laval style nozzle, which has a converging/diverging entrance fundamental to the micronozzle, is included in this review due to its importance in industrial applications in uranium (U) isotope refinement. Using advanced computational fluid dynamic (CFD) simulations, the extent of gas separation can be modelled. Herein, we first examine the literature and survey recent advances on fabrication techniques for creating curved micronozzles, methods and separation principles used to design devices. Furthermore, we then follow with highlights of CFD simulations applied to evaluate the separation effects using microsieves. Finally, identification of the gap and recommendation for future development and applications are suggested for using intrinsic molecular features and fluidic dynamics in formulating separation strategies.

30 Microfluidics↗
DOE OSTI2023

Transferring a Molecular Foundation Model for Polymer Property Predictions

Transformer-based large language models have remarkable potential to accelerate design optimization for applications such as drug development and material discovery. Self-supervised pretraining of transformer models requires large-scale data sets, which are often sparsely populated in topical areas such as polymer science. Further, state-of-the-art approaches for polymers conduct data augmentation to generate additional samples but unavoidably incur extra computational costs. In contrast, large-scale open-source data sets are available for small molecules and provide a potential solution to data scarcity through transfer learning. In this work, we show that using transformers pretrained on small molecules and fine-tuned on polymer properties achieves comparable accuracy to those trained on augmented polymer data sets for a series of benchmark prediction tasks.

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

Advanced Modeling and Process-Materials Co-Optimization Strategies for Swing Adsorption Based Gas Separations

This project devised a computational framework for simultaneously co-optimizing pressure swing adsorption process designs along with the sorbent materials (specifically, metal-organic frameworks) to be employed in the associated packed bed columns. The materials optimization aspect involved search over a design space that can describe the material’s molecular structure, while the process optimization aspect considered various process degrees of freedom for steps arising in various cycle configurations. This framework was demonstrated on the separation of nitrogen and carbon dioxide, which arises ubiquitously in a multitude of post-combustion carbon capture and “blue” hydrogen production applications. Our results led to metal-organic framework molecular descriptor choices that are predicted to outperform standard structures used in practice, providing guidance for future metal-organic framework synthesis efforts.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS1995

Application of Direct Simulation Monte Carlo to Satellite Contamination Studies

A novel method is presented to estimate contaminant levels around spacecraft and satellites of arbitrarily complex geometry. The method uses a three-dimensional direct simulation Monte Carlo algorithm to characterize the contaminant cloud surrounding the space platform, and a computer-assisted design preprocessor to define the space-platform geometry. The method is applied to the Upper Atmosphere Research Satellite to estimate the contaminant flux incident on the optics of the halogen occultation experiment (HALOE) telescope. Results are presented in terms of contaminant cloud structure, molecular velocity distribution at HALOE aperture, and code performance.

Rault, Didier F. G.↗
NASA NTRS1997

Peptide Folding and Translocation Across the Water-Membrane Interface

The ability of small peptides to organize at aqueous interfaces was examined by performing a series of large-scale, molecular dynamics computer simulations of several peptides composed of two amino acids, nonpolar leucine (L) and polar glutamine (Q). The peptides differed in size and sequence of the amino acids. Studies on dipeptides LL, LQ, QL and QQ were extended to two heptamers, LQQLLQL and LQLQLQL, designed to maximize interfacial stability of an alpha-helix and a beta-strand, respectively, by exposing polar side chains to water and nonpolar side chains to a nonpolar phase. Finally, a transition of an undecamer, composed entirely of leucine residues, from a disordered structure in water to an alpha-helix in a nonpolar phase representing the interior of the membrane was investigated. Complete folding of a peptide in solution was accomplished for the first time in computer simulations. The simulations revealed several basic principles governing the sequence-dependent organization of peptides at interfaces. Short peptides tend to accumulate at interfaces and acquire ordered structures, providing that they have a proper sequence of polar and nonpolar amino acids. The dominant factor determining the interfacial structure of peptides is the hydrophobic effect, which is manifested at aqueous interfaces as a tendency for polar and nonpolar groups of the solute to segregate into the aqueous and nonpolar phases, respectively. If peptides consist of nonpolar residue's only, they become inserted into the nonpolar phase. As demonstrated by the example of the leucine undecamer, such peptides fold into an alpha-helix as they partition into the nonpolar medium. The folding proceeds through an intermediate, called 3-10-helix, which remains in equilibrium with the alpha-helix. Once in the nonpolar environment, the peptides can readily change their orientation with respect to the interface from parallel to perpendicular, especially in response to local electric fields. The ability of nonpolar peptides to modify both the structure and orientation with respect to the interface from parallel to perpendicular, especially in response to local electric fields. The ability of nonpolar peptides to modify both the structure and orientation with changing external conditions may have provided a simple mechanism of transmitting signals from the environment to the interior of a cell.

Pohorille, Andrew↗
NASA NTRS2012

Lattice Thermal Conductivity from Atomistic Simulations: ZrB2 and HfB2

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗
NASA NTRS2012

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗
DOE OSTI2021

Molecular design and fabrication of PIM-1/polyphosphazene blend membranes with high performance for CO 2 /N 2 separation

New polymeric blend membranes for CO 2 separation were synthesized based on insights from molecular dynamics simulation. A molecular-level structure-property relationship in polymers of intrinsic microporosity (PIM) based blend membranes, was investigated in detail computationally. Calculated local density profiles and energy of interaction of the blend membranes, composed of PIM-1 and various polyphosphazenes, showed that using the polyphosphazene with a higher concentration of ether side chains can improve the compatibility with PIM-1. Furthermore, based on the findings of computational studies, blend membranes were experimentally fabricated from PIM-1 and polyphosphazenes with various polyether side chain concentrations. Polyether concentration in polyphosphazenes was correlated with the film properties and gas transport performance of the blend membranes. Blend membranes showed very high CO 2 permeability (3100-5300 barrer) and improved CO 2 /N 2 selectivity (24-28), outperforming all other PIM-based blend membranes reported to date. Moreover, the CO 2 permeability performance of the blend membranes was tested 566 hours under real post-combustion flue gas from a coal-fired power plant, including CO 2 , N 2 , H 2 O, O 2 , SO x and NO x .

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

Deep Learning Approaches for Predicting the Surface Tension of Ionic Liquids

Ionic liquids (ILs) are a novel class of solvents that have attracted significant attention due to their unique and tunable properties. Among their physiochemical characteristics, surface tension plays a critical role in various industrial applications including electrolytes, heat transfer fluids, and separation processes. However, because of the exploratory nature of IL design and the vast combinatorial space of possible anion–cation pairs, the experimental determination of these properties is often impractical, being both time-consuming and costly. To overcome these challenges, computational approaches are increasingly employed to develop accurate predictive models that can accelerate IL discovery and design. In this study, we present two deep learning (DL) models for predicting the surface tension of ILs across a broad temperature range at a constant pressure. The models use simplified molecular input line entry system, SMILES, representations of ILs to extract molecular features as inputs. Both DL models demonstrate excellent agreement with experimental data, achieving an R 2 value of 0.990 and a root-mean-square error of 0.792 mN/m. In conclusion, these results offer valuable insights for the rapid screening and rational design of ILs with tailored surface tension values.

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

Competitive Pi-Stacking and H-Bond Piling Increase Solubility of Heterocyclic Redoxmers

Redoxmers are organic molecules that carry electric charge in flow batteries. In many instances, they consist of heteroaromatic moieties modified with appended groups to prevent stacking of the planar cores and increase solubility in liquid electrolytes. This higher solubility is desired as it potentially allows achieving greater energy density in the battery. However, the present synthetic strategies often yield bulky molecules with low molarity even when they are neat and still lower molarity in liquid solutions. Fortunately, there are exceptions to this rule. Here, we examine one well-studied redoxmer, 2,1,3-benzothiadiazole, which has solubility similar to 5.7 M in acetonitrile at 25 degrees C. We demonstrate computationally and prove experimentally that the competition between two packing motifs, face-to-face it-stacking and random N-H bond piling, introduces frustration that confounds nucleation in crowded solutions. Our findings and examples from related systems suggest a complementary strategy for the molecular design of redoxmers for high energy density redox flow cells.

25 ENERGY STORAGE↗