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

Hoobas: A highly object-oriented builder for molecular dynamics

Polydispersity and random sequences are ubiquitous features of polymers, and molecular dynamics simulations can help elucidate the impact of disorder in polymer systems. However, currently available packages for building polymer topologies do not enable the user to include randomness in a straightforward fashion. Here, we introduce Hoobas, a molecular builder package that easily handles polydispersity using a prototype-builder design pattern. This enables fast and easy building of systems comprised of thousands of distinct objects. It is written in the Python programming language, which ensures compatibility with a wide range of molecular dynamics packages and tools, as well as easy integration into most workflows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Understanding Spin Coherence in Polyoxometalate-Based Molecular Qubits

This research advances the understanding of spin decoherence mechanisms and addresses the grand challenges in quantum information science of generating and stabilizing quantum states that may be manipulated with single-atom precision. It is accomplished using unique capabilities developed at PNNL for ion soft landing and spatially resolved vibrational spectroscopy of well-defined supported molecular qubit (MQ) arrays. It addresses knowledge gaps impeding the development of MQ-based quantum devices by answering two scientific questions: 1) Can we minimize spin-lattice relaxation rates by tuning the substrate-MQ interface? 2) What is the optimum number of spin centers that enables sufficient coherent signal generation prior to decoherence through inter-MQ and MQ-support interactions? The objective is to develop predictive design principles for prolonging spin coherence lifetimes in arrays of optically-addressable MQs for quantum computing applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

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

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

25 ENERGY STORAGE↗
DOE OSTI2021

Analytical gradients for molecular-orbital-based machine learning

We report molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.

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

Fast methods for multisite charge transfer processes. I. Constrained, state averaged CASSCF(1,n) and CASSCF(2n − 1,n) simulations

We design a dynamically weighted state-averaged constrained complete active space self-consistent field (DW-SA-cCASSCF) algorithm to treat electrons or holes moving between n molecular fragments (where n can be larger than 2). Within such a so-called eDSCn/hDSCn approach, we consider configurations that are mutually single excitations of each other, and we apply a generalized set of constraints to tailor the method for studying charge transfer problems. The constrained optimization problem is efficiently solved using a DIIS-SQP algorithm, thus maintaining computational efficiency. We demonstrate the method for a finite Su–Schrieffer–Heeger chain, successfully reproducing the expected exponential decay of diabatic couplings with distance. When combined with a gradient, the current extension immediately enables efficient nonadiabatic dynamics simulations of complex multi-state charge transfer processes.

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

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

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

First-principles-derived transport properties of Molten chloride salts

Molten salts have many applications ranging from a heat transfer medium in both generation IV nuclear reactor designs and the solar industry to thermal storage systems. While molten salts show promising properties for these applications, there still exists a knowledge gap for the transport properties of molten salts at elevated temperatures. This work uses ab initio Molecular Dynamics to investigate the transport properties of KCl, LiCl, KCl-LiCl eutectic, NaCl, MgCl 2 , and NaCl-MgCl 2 eutectic molten salt systems. Further, the properties presented here are the diffusion coefficient, viscosity, and isochoric heat capacity. These properties are compared to experimental data where available and other computational work in cases where no experimental data is available. This is the first work to explore timescales over 100 ps via AIMD for the determination of transport properties in molten salts.

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

w19_OMEC Scientific Highlight [Slides]

Scientific Achievement: Revealed the effect of graphite substrate on H binding / hydrogen evolution reaction (HER) activity of organic molecular electrocatalysts (OMECs); and demonstrated the effect of molecular modification on calculated HER activity descriptor. Significance and Impact: Understanding of HER processes is important for discovering low-cost OMECs for hydrogen production and energy application; design of OMEC structures remains challenging due to the unexplored redox properties of organic molecules; and density functional theory (DFT) has been proposed to provide needed insights into HER activity for improved materials design. Research Details: DFT calculations on the binding energy of HER intermediates are used for predicting reduction potentials and and calculated H adsorption energy as activity descriptor; and LANL Institutional Computing (IC) resources are essential for successful project execution given the size of structures involved in modeling molecule-support interaction.

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

Thermodynamic and Kinetic Activity Descriptors for the Catalytic Hydrogenation of Ketones

Activity descriptors are a powerful tool for the design of catalysts than can efficiently utilize H 2 with minimal energy losses. In this study, we develop the use of hydricity and H - self-exchange rates as thermodynamic and kinetic descriptors for the hydrogenation of ketones by molecular catalysts. Two complexes with known hydricity, HRh(dmpe) 2 and HCo(dmpe) 2 , were investigated for the catalytic hydrogenation of ketones under mild conditions (1.5 atm, 25 °C). The rhodium catalyst proved to be an efficient catalyst for a wide range of ketones, whereas the cobalt catalyst could only hydrogenate electron-deficient ketones. Using a combination of experiment and electronic structure theory, thermodynamic hydricity values were established for 46 alkoxide/ketone pairs in both MeCN and THF solvent. Through comparison of the hydricities of the catalysts and substrates, it was determined that catalysis was only observed for catalyst/ketone pairs with an exergonic H - transfer step. Mechanistic studies revealed that H - transfer was rate-limiting step for catalysis, allowing for the experimental and computation construction of linear free-energy relationships (LFERs) for H - transfer. Further analysis revealed the LFERs could be reproduced using Marcus theory, in which the H - self-exchange rates for the HRh/Rh + and ketone/alkoxide pairs were used to predict the experimentally measured catalytic barriers within 2 kcal mol -1 . Finally, these studies significantly expand the scope of catalytic reactions that can be analyzed with a thermodynamic hydricity descriptor and firmly establish Marcus theory as a valid approach to develop kinetic descriptors for designing catalysts for H - transfer reactions.

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

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗
DOE OSTI2025

Deconstructing Chirality: Probing Local and Nonlocal Effects in Azobenzene Derivatives with X-ray Circular Dichroism

Resolving molecular chirality on the atomic scale remains a critical challenge in chemistry. Conventional Optical Circular Dichroism spectroscopy often overlooks subtle and localized structural features. Here, we computationally investigate site-specific X-ray circular dichroism (XCD) across a series of trans-azobenzene derivatives to deconstruct and interpret chiroptical signals at the atomic level. Our modeling reveals that XCD is capable of distinguishing dichroic contributions arising from both a local chiral center and a global molecular twist, revealing their intricate interplay and potential for constructive or destructive interference. We show that sterically induced global distortions can dominate the XCD signal in some cases, even suppressing the response from the chiral center itself. This insight suggests a new molecular design principle for tuning the chiroptical activity, which we extend by proposing strategies to achieve unidirectional photoisomerization through steric gearing. Altogether, this work establishes a quantitative framework for engineering chiroptical responses, laying the foundation for the design of functional chiral systems utilizing principles of unidirectional molecular motor-like conformational dynamics.

X-ray spectroscopy↗
DOE OSTI2023

Ensuring thermodynamic consistency with invertible coarse-graining

Coarse-grained models are a core computational tool in theoretical chemistry and biophysics. A judicious choice of a coarse-grained model can yield physical insights by isolating the essential degrees of freedom that dictate the thermodynamic properties of a complex, condensed-phase system. The reduced complexity of the model typically leads to lower computational costs and more efficient sampling compared with atomistic models. Designing “good” coarse-grained models is an art. Generally, the mapping from fine-grained configurations to coarse-grained configurations itself is not optimized in any way; instead, the energy function associated with the mapped configurations is. In this work, we explore the consequences of optimizing the coarse-grained representation alongside its potential energy function. We use a graph machine learning framework to embed atomic configurations into a low-dimensional space to produce efficient representations of the original molecular system. Because the representation we obtain is no longer directly interpretable as a real-space representation of the atomic coordinates, we also introduce an inversion process and an associated thermodynamic consistency relation that allows us to rigorously sample fine-grained configurations conditioned on the coarse-grained sampling. We illustrate that this technique is robust, recovering the first two moments of the distribution of several observables in proteins such as chignolin and alanine dipeptide.

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

Final Report

This final progress report describes the development of the CANDO Software package at ThirdLaw Technologies LLC (TLT) towards a software platform that will be used by its subsidiary company ThirdLaw Molecular LLC (TLM) for the development of therapeutic Spiroligomer™ molecules. ThirdLaw Molecular LLC is 95% owned by ThirdLaw Technologies LLC. TLM is developing a new therapeutics platform based on Spiroligomer™ molecules and is currently funded by the Department of Defense, Defense Threat Reduction Agency (DOD/DTRA W15QKN209C004). Spiroligomer™ molecules are complex, highly-preorganized macromolecules with programmable three-dimensional shapes. Spiroligomer™ molecules are currently synthesized by TLM in the form of large (>880,000 member) DNA-encoded libraries and screened against dangerous bacteria and viruses as part of the DOD/DTRA project. The CANDO software package is being developed to enable the computational design of Spiroligomer™ molecules for the development of therapeutic and diagnostic molecules. Towards this end, over the period of the grant we implemented software algorithms and a user interface within CANDO to enable the computational design of Spiroligomer™ molecules.

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

Replication Data for: Deconstructing Chirality: Probing Local and Non-local Effects in Azobenzene Derivatives with X-ray Circular Dichroism

Resolving molecular chirality at the atomic scale remains a critical challenge in chemistry. Conventional Optical Circular Dichroism spectroscopy often overlooks subtle and localized structural features. Here, we computationally investigate site-specific X-ray Circular Dichroism (XCD) across a series of trans-azobenzene derivatives to deconstruct and interpret chiroptical signals at the atomic level. Our modeling reveals that XCD is capable of distinguishing dichroic contributions arising from both a local chiral center and global molecular twist, revealing their intricate interplay and potential for constructive or destructive interference. We show that sterically-induced global distortions can dominate the XCD signal in some cases, even suppressing the response from the chiral center itself. This insight suggests a new molecular design principle for tuning chiroptical activity, which we extend by proposing strategies to achieve unidirectional photoisomerization through steric gearing. Altogether, this work establishes a quantitative framework for engineering chiroptical responses, laying the foundation for the design of functional chiral systems utilizing principles of unidirectional molecular motor-like conformational dynamics.

Chemistry↗
DOE OSTI2022

Predictive coarse-grained (CG) modeling of morphologies in polymer nanocomposites with specific and directional intermolecular interactions (Final Report)

The overarching goal of the proposed work was to develop predictive models for investigating structure and dynamics in soft materials with chemistries that have specific and directional molecular interactions. The motivation behind studying materials with specific and directional interactions lies in the many desirable features these interactions provide when designing novel soft materials. Soft materials with specific and direction interactions (such as hydrogen bonds or H-bonds) can have a) thermally reversible phase behavior with different functions with varying temperature, b) precisely tuned nanostructure with desirable geometries that afford unique physical properties (e.g., color response, mechanical properties) and c) well-mixed/blended morphologies that are useful for variety of applications in energy field (e.g., materials for batteries require use of blended polymers where one domain gives superior mechanical properties and one domain promotes electrical conduction). Notably, biology makes extensive use of specific and directional interactions, in many cases based on H-bonds, to construct materials with precisely defined architectures and properties. Engineered soft materials with precisely tuned nanostructures and improved processiblity through thermoresponsive phase behavior are useful in numerous applications that are relevant to the Department of Energy (DOE) including high efficiency electronic devices, light-weight high-strength composite materials for batteries and fuel cells, and polymer membranes for separations, etc. While past computational studies have been tremendously useful in understanding molecular phenomena and guiding synthesis of new macromolecular soft materials for a wide variety of applications, the inability to capture small scale specific and directional interactions alongside macromolecular length and time scales represented a key limitation of most studies to date. Our work in this project addressed this grand challenge in computational materials chemistry, i.e., ability to model the anisotropic, directional, and specific interactions that govern the behavior of many macromolecular soft matter systems of interest, thus, has greatly expanded the predictive potential of simulations. Specifically the key outcomes were: successful development of new coarse-grained (CG) polymer models to study generic and specific polymer chemistries in which hydrogen-bonding interactions are dominant. These CG models were then used in molecular simulations to study structure and thermodynamics in polymer nanocomposites and blends; some studies were conducted in collaboration with experimentalists. We also published a perspective and a viewpoint which included some of the work we completed in this DOE project; we believe these perspective and viewpoint articles guide other researchers in the soft materials community on how to extend the computational approaches and models we have developed for the purposes of their studies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Autonomous continuous flow reactor synthesis for scalable atom-precision

With new instrumentation design, robotics, and in-operando hyphenated analytical tool automation, the intelligent discovery of synthesis pathways is becoming feasible. It can potentially bridge the gap for the scale-up of new materials. In this article, we review current progress and describe a new system that uses an autonomous continuous flow chemistry framework to translate high-quality lead molecules and materials to quantities that can meet scalability demands. At the core is a continuous flow synthesis platform that can design its viable synthesis pathway to a particular molecule or material and then autonomously carry it out. This is realized by integrating: (1) A workflow/architecture for multimode chemical/materials characterization in-line. The in-line characterization modes are NMR, ESR, IR, Raman, UV-Vis, GC-MS, and HPLC, along with ex-situ modes for X-Ray and neutron scattering; (2) Integration for feedback/analysis/data storage of the control variables; (3) A core software stack that includes deep learning and reinforcement learning alongside quantum chemistry and molecular dynamics; (4) On-demand compute architectures that parse calculations to compute resources needed which include light-weight edge, mid-level edge (NVIDA DGX-2), and high-performance computing. We demonstrate preliminary results on how this autonomous reactor system can enhance our ability to deliver deuterated materials, copolymers, and site-substituted molecules.

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

In Silico Discovery of Multistep Chemistry Initiated by a Conical Intersection: The Challenging Case of Donor–Acceptor Stenhouse Adducts

Detailed mechanistic understanding of multistep chemical reactions triggered by internal conversion via a conical intersection is a challenging task that emphasizes limitations in theoretical and experimental techniques. We present a discovery-based, hypothesis-free computational approach based on first-principles molecular dynamics to discover and refine the switching mechanism of donor–acceptor Stenhouse adducts (DASAs). We simulate the photochemical experiment in silico, following the “hot” ground state dynamics for 10 ps after photoexcitation. Using state-of-the-art graphical processing units-enabled electronic structure calculations we performed in total ~2 ns of nonadiabatic ab initio molecular dynamics discovering (a) critical intermediates that are involved in the open-to-closed transformation, (b) several competing pathways which lower the overall switching yield, and (c) key elements for future design strategies. Finally, our dynamics describe the natural evolution of both the nuclear and electronic degrees of freedom that govern the interconversion between DASA ground-state intermediates, exposing significant elements for future design strategies of molecular switches.

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

Mechanical Roles of Polysaccharide Assembly and Interactions in Plant Cell Walls

Plants synthesize polysaccharide-based primary cell walls that possess unique microstructures and mechanical properties to accommodate plant growth and provide protection. Here, it remains challenging to assess the role of polysaccharide organization and interactions in the mechanical behavior of primary cell walls owing to their complex microstructure and highly nonlinear mechanical responses. Employing a coarse-grained molecular dynamics model developed for onion epidermal walls, this work explores the conditions under which polysaccharide assembly and interactions might play a significant role in primary cell wall mechanics. Cellulose–cellulose adhesion plays a dominant role in the wall load-bearing capacity, but when cellulose–cellulose adhesion was disrupted computationally, cellulose–xyloglucan adhesion could influence the wall load-bearing capacity. Contrary to the common concept that xyloglucans mechanically tether well-separated cellulose microfibrils, xyloglucans functioned in this case as interfibrillar adhesives capable of transmitting tensile forces between cellulose microfibrils. Our findings may inform design criteria of new materials inspired by plant cell walls.

59 BASIC BIOLOGICAL SCIENCES↗