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 451 records · Page 25

In Silico Screening of CO 2 –Dipeptide Interactions for Bioinspired Carbon Capture

Carbon capture, sequestration and utilization offers a viable solution for reducing the total amount of atmospheric CO 2 concentrations. On an industrial scale, amine-based solvents are extensively employed for CO 2 capture through chemisorption. Nevertheless, this method is marked by the high cost associated with solvent regeneration, high vapor pressure, and the corrosive and toxic attributes of by-products, such as nitrosamines. An alternative approach is the biomimicry of sustainable materials that have strong affinity and selectivity for CO 2 . Bioinspired approaches, such as those based on naturally occurring amino acids, have been proposed for direct air capture methodologies. In this study, we present a database consisting of 960 dipeptide molecular structures, composed of the 20 naturally occurring amino acids. Furthermore, those structures were analyzed with a novel computational workflow presented in this work that considers certain interaction sites that determine CO 2 affinity. Density functional theory (DFT) and symmetry-adapted perturbation theory (SAPT) computations were performed for the calculation of CO 2 interaction energies, which allowed to limit our search space to 400 unique dipeptide structures. Using this computational workflow, we provide statistical insights into dipeptides and their affinity for CO 2 binding, as well as design principles that can further enhance CO 2 capture through cooperative binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Challenges for density functional theory: calculation of CO adsorption on electrocatalytically relevant metals

Density Functional Theory (DFT) is currently the most tractable choice of theoretical model used to understand the mechanistic pathways for electrocatalytic processes such as CO 2 or CO reduction. In this work, we assess the performance of two DFT functionals designed specifically to describe surface interactions, RTPSS and RPBE, as well as two popular meta-GGA functionals, SCAN and B97M-rV, that have not been a priori optimized for better interfacial properties. We assess all four functionals against available experimental data for prediction of bulk and bare surface properties on four electrocatalytically relevant metals, Au, Ag, Cu, and Pt, and for binding CO to surfaces of these metals. To partially mitigate issues such as thermal and anharmonic corrections associated with comparing computations with experiments, molecular benchmarks against high level quantum chemistry are reported for CO complexes with Au, Ag, Cu and Pt atoms, as well as the CO–water complex and the water dimer. Overall, we find that the surface modified RPBE functional performs reliably for many of the benchmarks examined here, and the meta-GGA functionals also show promising results. Specifically B97M-rV predicts the correct site preference for CO binding on Ag and Au (the only functional tested here to do so), while RTPSS performs well for surface relaxations and binding of CO on Pt and Cu.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exchange, interpretation, and database-search of ion mobility spectra supported by data format JCAMP-DX

To assist peak assignment in ion mobility spectrometry it is important to have quality reference data. The reference collection should be stored in a database system which is capable of being searched using spectral or substance information. We propose to build such a database customized for ion mobility spectra. To start off with it is important to quickly reach a critical mass of data in the collection. We wish to obtain as many spectra combined with their IMS parameters as possible. Spectra suppliers will be rewarded for their participation with access to the database. To make the data exchange between users and system administration possible, it is important to define a file format specially made for the requirements of ion mobility spectra. The format should be computer readable and flexible enough for extensive comments to be included. In this document we propose a data exchange format, and we would like you to give comments on it. For the international data exchange it is important, to have a standard data exchange format. We propose to base the definition of this format on the JCAMP-DX protocol, which was developed for the exchange of infrared spectra. This standard made by the Joint Committee on Atomic and Molecular Physical Data is of a flexible design. The aim of this paper is to adopt JCAMP-DX to the special requirements of ion mobility spectra.

Baumback, J. I.↗

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↗

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↗

Design and Off-Design Performance of 100 kWe-Class Brayton Power Conversion Systems

The NASA Glenn Research Center in-house computer model Closed Cycle Engine Program (CCEP) was used to explore the design trade space and off-design performance characteristics of 100 kWe-class recuperated Closed Brayton Cycle (CBC) power conversion systems. Input variables for a potential design point included the number of operating units (1, 2, 4), cycle peak pressure (0.5, 1, 2 MPa), and turbo-alternator shaft speed (30, 45, 60 kRPM). The design point analysis assumed a fixed turbine inlet temperature (1150 K), compressor inlet temperature (400 K), helium-xenon working-fluid molecular weight (40 g/mol), compressor pressure ratio (2.0), recuperator effectiveness (0.95), and a Sodium-Potassium (NaK) pumped-loop radiator. The design point options were compared on the basis of thermal input power, radiator area, and mass. For a nominal design point with defined Brayton components and radiator area, off-design cases were examined by reducing turbine inlet temperature (as low as 900 K), reducing shaft speed (as low as 50 percent of nominal), and circulating a percentage (up to 20 percent) of the compressor exit flow back to the gas cooler. The off-design examination sought approaches to reduce thermal input power without freezing the radiator.

Johnson, Paul K.↗

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↗

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↗

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↗

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↗

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↗

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↗

A Synthesis of Fluid Dynamics and Quantum Chemistry for the Design of Nanoelectronics

In 1959, during a famous lecture entitled "There's Plenty of Room at the Bottom", Richard Feynman focused on the startling technical possibilities that would exist at the limit of miniaturization, that being atomically precise devices with dimensions in the nanometer range. A nanometer is both a convenient unit of length for medium to large sized molecules, and the root of the name of the new interdisciplinary field of "nanotechnology". Essentially, "nanoelectronics" denotes the goal of shrinking electronic devices, such as diodes and transistors, as well as integrated circuits of such devices that can perform logical operations, down to dimensions in the range of 100 nanometers. The thirty-year hiatus in the development of nanotechnology can figuratively be seen as a period of waiting for the bottom-up and atomically precise construction skills of synthetic chemistry to meet the top-down reductionist aspirations of device physics. The sub-nanometer domain of nineteenth-century classical chemistry has steadily grown, and state-of-the-art supramolecular chemistry can achieve atomic precision in non-repeating molecular assemblies of the size desired for nanotechnology. For nanoelectronics in particular, a basic understanding of the electron transport properties of molecules must also be developed. Quantum chemistry provides powerful computational methods that can accurately predict the properties of small to medium sized molecules on a desktop workstation, and those of large molecules if one has access to a supercomputer. Of the many properties of a molecule that quantum chemistry routinely predicts, the ability to carry a current is one that had not even been considered until recently. "Currently", there is a controversy over just how to define this key property. Reminiscent of the situation in high-Tc superconductivity, much of the difficulty arises from the different models that are used to simplify the complex electronic structure of real materials. A model-independent approach has been proposed, that sacrifices the plentiful molecular orbitals and Bloch functions of conventional approaches, for a single three-dimensional observable quantity, the electron momentum density Pi(sub rho). This quantity is simply the probability of any electron having momentum rho, multiplied by the total number of electrons in the system (the position of the electron is uncertain). We have explored the utility of this new approach in providing a fundamental understanding of the electron transport properties of molecules that have provi been nominated as candidates for components in the design of nanoelectronics; phenylene-ethynylene oligomers. Some of the molecular systems that have been studied are sketched below.

MacDougall, Preston J.↗

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↗

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↗

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↗

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↗

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↗