SEARCH · Engineering Papers
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.
The advances in theory and computing technology over the last decade have led to enormous progress in applying atomistic molecular dynamics (MD) methods to the characterization, prediction, and design of chemical, biological, and material systems,.
Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.
Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.
Through evolution, nature has presented a set of remarkable protein materials, including elastins, silks, keratins and collagens with superior mechanical performances that play crucial roles in mechanobiology. However, going beyond natural designs to discover proteins that meet specified mechanical properties remains challenging. Here, we report a generative model that predicts protein designs to meet complex nonlinear mechanical property-design objectives. Our model leverages deep knowledge on protein sequences from a pretrained protein language model and maps mechanical unfolding responses to create proteins. Via full-atom molecular simulations for direct validation, we demonstrate that the designed proteins are de novo, and fulfill the targeted mechanical properties, including unfolding energy and mechanical strength, as well as the detailed unfolding force-separation curves. Our model offers rapid pathways to explore the enormous mechanobiological protein sequence space unconstrained by biological synthesis, using mechanical features as the target to enable the discovery of protein materials with superior mechanical properties.
To move toward rational design of efficient organic light emitting diodes based on the radical idea of inverted singlet-triplet gap (INVEST) systems, we propose a set of novel quantum chemical approaches, predictive but low-cost, to unveil a set of structural-property relationships. We perform a computational study of a series of substituted molecules based on a small set of known INVEST molecules. Our study demonstrates a high degree of correlation between the intramolecular charge transfer and the singlet-triplet energy gap and hints towards the use of a quantitative estimate of charge transfer to predict and modulate these energy gaps. We aim to create a database of INVEST molecules that includes accurate benchmarks of singlet-triplet energy gaps. Furthermore, we aim to link structural features and molecular properties, enabling a control knob for rational design.
Understanding molecular motion in diffusion-driven complex environments is critical for designing sustainable materials and improving chemical processes. Here, we introduce a multidimensional nuclear magnetic resonance (NMR) method that captures how molecular populations exchange across different dynamic regimes. By extending the modulated gradient spin-echo technique to include frequency-frequency correlations, our approach reveals diffusion pathways that are otherwise obscured in heterogeneous systems. Implemented on a unilateral NMR magnet, the method eliminates gradient pulsing constraints and accesses dynamics in the kilohertz regime. We apply this technique to swelling and acid-catalyzed deconstruction of cross-linked and linear polymers to observe how structural heterogeneity evolves over time. By linking molecular motion to topology and chemical state, we extract physical metrics such as fractal surface dimensionality and reaction wavefront velocity, properties inaccessible with standard diffusion measurements. This work expands the capabilities of NMR for probing soft matter, with implications for polymer recycling and materials design.
The design rules for materials are clear for applications with a single objective. For most applications, however, there are often multiple, sometimes competing objectives where there is no single best material and the design rules change to finding the set of Pareto optimal materials. In this work, we leverage an active learning algorithm that directly uses the Pareto dominance relation to compute the set of Pareto optimal materials with desirable accuracy. We apply our algorithm to de novo polymer design with a prohibitively large search space. Using molecular simulations, we compute key descriptors for dispersant applications and drastically reduce the number of materials that need to be evaluated to reconstruct the Pareto front with a desired confidence. This work showcases how simulation and machine learning techniques can be coupled to discover materials within a design space that would be intractable using conventional screening approaches.
Directing the formation of nanoscale architectures from nanoparticles is one of the key challenges in designing nanomaterials with prescribed functions. Atomic systems, given their ability to form molecules and crystals via directional chemical bonds, provide an inspiration for establishing approaches where nanoparticles with designed anisotropic binding modalities can be assembled into nanoscale architectures. However, fabricating such nanoparticles has been challenging due to their small dimensions and limited ways for site-specific control of their surface. To this end, we present a Molecular Stamping (MOST) approach to pattern DNA-coated nanoparticles with molecules at the predefined positions on a nanoparticle surface. This patterning is realized by use of a rigid and coordinative DNA frame as a molecular stamping apparatus (MOST App). The MOST App transfers multiple types of molecular “inks”, DNA sequences, onto nanoparticle surface and fixes these molecular inks into place to form a designed pattern. After a nanoparticle release from MOST App, it possesses single-molecule patches that can provide anisotropic bonds with distinctive affinities. We further use these stamped nanoparticles to assemble prescribed clusters, whose structure is determined by the locations of patches. Using electron microscopy and tomographic methods, we investigate the efficiency of cluster formation and the resulting spatial arrangements of nanoparticles. The presented approach provides a single-molecule and spatiallydetermined control over nanoparticle functionalization for creating nanoparticles with designed placement of different molecules and for realizing a rational fabrication of nanomaterial architectures.
The article is an interview with Dr. Yi Liu. Dr. Liu is the facility director at the Molecular Foundry, Lawrence Berkeley National Laboratory. Trained as a supramolecular chemist and after a postdoctoral stay and working on click chemistry in Sharpless’s group he started his independent research career working on organic electronics, porous materials, and covalent organic frameworks (COFs). His aim is to develop nanostructured electronic materials through the design, synthesis and manipulation of tailor-made molecular constituents. He answers questions about research on covalent organic frameworks.
Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.
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.
Here, we present the theory, implementation, and benchmarking of a real-time time-dependent density functional theory (RT-TDDFT) module within the RMG code, designed to simulate the electronic response of molecular systems to external perturbations. Our method offers insights into nonequilibrium dynamics and excited states across a diverse range of systems, from small organic molecules to large metallic nanoparticles. Benchmarking results demonstrate excellent agreement with established TDDFT implementations and showcase the superior stability of our time integration algorithm, enabling long-term simulations with minimal energy drift. The scalability and efficiency of RMG on massively parallel architectures allow for simulations of complex systems, such as plasmonic nanoparticles with thousands of atoms. Future extensions, including nuclear and spin dynamics, will broaden the applicability of this RT-TDDFT implementation, providing a powerful toolset for studies of photoactive materials, nanoscale devices, and other systems where real-time electronic dynamics is essential.
Biological and heterogeneous catalysts for the electrochemical CO 2 reduction reaction (CO 2 RR) often exhibit a high degree of electronic delocalization that serves to minimize overpotential and maximize selectivity over the hydrogen evolution reaction (HER). Here, we report a molecular iron(II) system that captures this design concept in a homogeneous setting through the use of a redox non-innocent terpyridine-based pentapyridine ligand (tpyPY2Me). As a result of strong metal-ligand exchange coupling between the Fe(II) center and ligand, [Fe(tpyPY2Me)] 2+ exhibits redox behavior at potentials 640 mV more positive than the isostructural [Zn(tpyPY2Me)] 2+ analog containing the redox-inactive Zn(II) ion. This shift in redox potential is attributed to the requirement for both an open-shell metal ion and a redox non-innocent ligand. The metal-ligand cooperativity in [Fe(tpyPY2Me)] 2+ drives the electrochemical reduction of CO 2 to CO at low overpotentials with high selectivity for CO 2 RR (>90%) and turnover frequencies of 100,000 s -1 with no degradation over 20 h. The decrease in the thermodynamic barrier engendered by this coupling also enables homogeneous CO 2 reduction catalysis in water without compromising selectivity or rates. Synthesis of the two-electron reduction product, [Fe(tpyPY2Me)] o , and characterization by X-ray crystallography, Mössbauer spectroscopy, X-ray absorption spectroscopy (XAS), variable temperature NMR, and density functional theory (DFT) calculations, support assignment of an open-shell singlet electronic structure that maintains a formal Fe(II) oxidation state with a doubly reduced ligand system. Furthermore, this work provides a starting point for the design of systems that exploit metal-ligand cooperativity for electrocatalysis where the electrochemical potential of redox non-innocent ligands can be tuned through secondary metal-dependent interactions.
α 1 -adrenergic receptors (α 1 -ARs) play critical roles in the cardiovascular and nervous systems where they regulate blood pressure, cognition, and metabolism. However, the lack of specific agonists for all α 1 subtypes has limited our understanding of the physiological roles of different α 1 -AR subtypes, and led to the stagnancy in agonist-based drug development for these receptors. Here we report cryo-EM structures of α 1A -AR in complex with heterotrimeric G-proteins and either the endogenous common agonist epinephrine or the α 1A -AR-specific synthetic agonist A61603. These structures provide molecular insights into the mechanisms underlying the discrimination between α 1A -AR and α 1B -AR by A61603. Guided by the structures and corresponding molecular dynamics simulations, we engineer α 1A -AR mutants that are not responsive to A61603, and α 1B -AR mutants that can be potently activated by A61603. Together, these findings advance our understanding of the agonist specificity for α 1 -ARs at the molecular level, opening the possibility of rational design of subtype-specific agonists.
Abstract There has been considerable recent progress in designing new proteins using deep-learning methods 1–9 . Despite this progress, a general deep-learning framework for protein design that enables solution of a wide range of design challenges, including de novo binder design and design of higher-order symmetric architectures, has yet to be described. Diffusion models 10,11 have had considerable success in image and language generative modelling but limited success when applied to protein modelling, probably due to the complexity of protein backbone geometry and sequence–structure relationships. Here we show that by fine-tuning the RoseTTAFold structure prediction network on protein structure denoising tasks, we obtain a generative model of protein backbones that achieves outstanding performance on unconditional and topology-constrained protein monomer design, protein binder design, symmetric oligomer design, enzyme active site scaffolding and symmetric motif scaffolding for therapeutic and metal-binding protein design. We demonstrate the power and generality of the method, called RoseTTAFold diffusion (RFdiffusion), by experimentally characterizing the structures and functions of hundreds of designed symmetric assemblies, metal-binding proteins and protein binders. The accuracy of RFdiffusion is confirmed by the cryogenic electron microscopy structure of a designed binder in complex with influenza haemagglutinin that is nearly identical to the design model. In a manner analogous to networks that produce images from user-specified inputs, RFdiffusion enables the design of diverse functional proteins from simple molecular specifications.
The recognition of membrane separations as a vital technology platform for enhancing the efficiency of separation processes has been steadily increasing. Concurrently, 3D printing has emerged as an innovative approach to fabricating reverse osmosis membranes for water desalination and treatment purposes. This method provides a high degree of control over membrane chemistry and structural properties. In particular, when compared to traditional manufacturing techniques, 3D printing holds the potential to expedite customization, a feat that is typically achieved through conventional manufacturing methods but often involves numerous processes and significant costs. This review aims to present the current advancements in membrane manufacturing technology specifically tailored for water desalination purposes, with a particular focus on the development of 3D-printed membranes. A comprehensive analysis of recent progress in 3D-printed membranes is provided. However, conducting experimental work to investigate various influential factors while ensuring consistent results poses a significant challenge. To address this, we explore how membrane manufacturing processes and performance can be effectively pre-designed and guided through the use of molecular dynamics simulations. Finally, this review outlines the challenges faced and presents future perspectives to shed light on research directions for optimizing membrane manufacturing processes and achieving optimal membrane performance.
Here, to explore how anion substitution modifies the existing magnetism in strongly correlated oxides, we investigate local electronic states and magnetic ordering in nickel oxide (NiO) induced by substituting oxygen (O) with nitrogen (N). Each N introduces an additional N 2p hole and modifies the magnetic moment of a neighboring nickel (Ni) cation site, as the exchange interaction between this hole and the Ni $e_{g}$ electrons exceeds the Ni-O-Ni superexchange interaction. This leads to the formation of Ni-N-Ni centers consisting of five spins, without perturbing the antiferromagnetic NiO lattice. These centers are studied using density functional theory and confirmed through high-resolution spectroscopy on N-substituted NiO thin films grown by molecular beam epitaxy. This type of magnetic design may enable future advances in quantum technologies based on strongly correlated materials, such as quantum sensors and spin-based qubits.
Conducting long-term simulations of small molecules using Molecular Dynamics (MD) is crucial in drug design. However, traditional methods to accelerate the process including ASICs or GPUs, have limitations. ASICs are difficult to access, while publicly accessible GPU packages for MD have issues with more GPUs when processing small molecules, causing decreased performance. Conversely, FPGAs are renowned for their scalability and have been commonly used as accelerators for various applications, making them suitable for the task. Moreover, with the availability of FPGAs in the cloud, pharmaceutical developers can now access FPGA devices. The only missing piece is a design deployed for MD acceleration. Therefore, we present FASDA, the first FPGA-based MD accelerator available for community development. FASDA is designed to evaluate range-limited MD, which is the most resource-intensive and computationally demanding component in MD. It outperforms the state of-the-art GPU solution by 4.67x, significantly reducing the time for MD in drug design.