Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular chemistry”

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 145 records · Page 8

Structural insights into β-1,3-glucan cleavage by a glycoside hydrolase family

The fundamental and assorted roles of β-1,3-glucans in nature are underpinned on diverse chemistry and molecular structures, demanding sophisticated and intricate enzymatic systems for their processing. Herein, the selectivity and modes of action of a glycoside hydrolase family active on β-1,3-glucans were systematically investigated combining sequence similarity network, phylogeny, X-ray crystallography, enzyme kinetics, mutagenesis and molecular dynamics. This family exhibits a minimalist and versatile (α/β)-barrel scaffold, which can harbor distinguishing exo or endo modes of action, including an ancillary-binding site for the anchoring of triple-helical β-1,3-glucans. The substrate binding occurs via a hydrophobic knuckle complementary to the canonical curved conformation of β-1,3-glucans or through a substrate conformational change imposed by the active-site topology of some fungal enzymes. Together, these findings expand our understanding of the enzymatic arsenal of bacteria and fungi for the breakdown and modification of β-1,3-glucans, which can be exploited for biotechnological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unraveling the Heterogeneous but Ordered Microstructure of the Nonionic Deep Eutectic Solvent Formed by Lauric Acid and N -Methylacetamide

The nonionic deep eutectic solvent, formed by lauric acid (LA) and N-methylacetamide (NMA), has been shown to have a heterogeneous molecular structure in which the LA and NMA form nonpolar and polar domains, respectively. Previous vibrational spectroscopy experiments demonstrated that the ability of the LA domains to solvate compounds was limited to long carbon chains, whereas other nonpolar molecules, such as W(CO) 6 , were found to be solvated by both LA and NMA. These experiments were not fully compatible with the previously proposed micelle-like structure of the nonpolar domains of the LA-NMA DES. In this work, the modeling of the DES molecular structure is pursued using classical molecular dynamics simulations. The new classical model reproduces both the SAXS structural factors and the previously experimentally derived interaction map for these LA-NMA DESs. In addition, the simulation also shows that LA-NMA DESs form highly organized LA aggregates that are difficult to disorganize. Further evidence of the correct description provided by the newly derived model is obtained using a moderately polar probe: chloroform-d. Computations using the classical model have a good agreement with the solvation behavior of the probe derived from experiments, in which the location of the probe is found to be mostly within the polar domain of the DES. The computational model also demonstrates that the probe solvation is a consequence of the tightly packed LA structure, which causes nonpolar molecules to be located at the interphase of the DES nonpolar domains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zero-Field NMR and Millitesla-SLIC Spectra for >200 Molecules from Density Functional Theory and Spin Dynamics

NMR is usually performed at magnetic fields of 1 T and above to obtain sufficient sensitivity and spectral dispersion to identify chemicals based on chemical shifts and J couplings. At lower fields, the advent of hyperpolarization technologies and sensitive detectors can address sensitivity concerns. However, it remains disputed whether spectral signatures at zero and ultra-low fields are sufficient for chemical identification. Here, we report an all–electron DFT-based batch calculation of J-coupling constants, which are used to generate J coupling NMR spectra at zero field and 6.5 mT for over 200 small molecules. In the developed computational tool chain, we first used the all-electron FHI-aims code to calculate the molecular J couplings and chemical shifts. We then fed the calculated NMR parameters into the NMR simulation package SPINACH to simulate both heteronuclear J coupling spectra at zero-field, and homonuclear J coupling spectra as spin-lock induced crossing (SLIC) spectra at ultra-low field (6.5 mT). The resulting spectra demonstrate that zero and ultra-low field NMR spectra can represent unique identifiers of chemical structure for small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shaping the Water-Harvesting Behavior of Metal–Organic Frameworks Aided by Fine-Tuned GPT Models

Here, we construct a data set of metal–organic framework (MOF) linkers and employ a fine-tuned GPT assistant to propose MOF linker designs by mutating and modifying the existing linker structures. This strategy allows the GPT model to learn the intricate language of chemistry in molecular representations, thereby achieving an enhanced accuracy in generating linker structures compared with its base models. Aiming to highlight the significance of linker design strategies in advancing the discovery of water-harvesting MOFs, we conducted a systematic MOF variant expansion upon state-of-the-art MOF-303 utilizing a multidimensional approach that integrates linker extension with multivariate tuning strategies. We synthesized a series of isoreticular aluminum MOFs, termed Long-Arm MOFs (LAMOF-1 to LAMOF-10), featuring linkers that bear various combinations of heteroatoms in their five-membered ring moiety, replacing pyrazole with either thiophene, furan, or thiazole rings or a combination of two. Beyond their consistent and robust architecture, as demonstrated by permanent porosity and thermal stability, the LAMOF series offers a generalizable synthesis strategy. Importantly, these 10 LAMOFs establish new benchmarks for water uptake (up to 0.64 g g –1 ) and operational humidity ranges (between 13 and 53%), thereby expanding the diversity of water-harvesting MOFs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

General Protocol for the Accurate Prediction of Molecular 13 C/ 1 H NMR Chemical Shifts via Machine Learning Augmented DFT

An accurate prediction of NMR chemical shifts at affordable computational cost is very important for different types of structural assignments in experimental studies. Density functional theory (DFT) and gauge-including atomic orbital (GIAO) are two of the most popular computational methods for NMR calculation, yet they often fail to resolve ambiguities in structural assignments. In this work, we present a new method that uses machine learning (ML) techniques (DFT + ML) that significantly increases the accuracy of 13 C/ 1 H NMR chemical shift prediction for a variety of organic molecules. The input of the generalizable DFT + ML model contains two critical parts: one is a vector providing insights into chemical environments, which can be evaluated without knowing the exact geometry of the molecule; the other one is the DFT-calculated isotropic shielding constant. The DFT + ML model was trained with a data set containing 476 13 C and 270 1 H experimental chemical shifts. For the DFT methods used here, the root mean square deviations (RMSDs) for the errors between predicted and experimental 13 C/ 1 H chemical shifts can be as small as 2.10/0.18 ppm, which is much lower than those from simple DFT (5.54/0.25 ppm), or DFT + linear regression (LR) (4.77/0.23 ppm) approaches. It also has a smaller maximum absolute error than two previously proposed NMR-predicting ML models. The robustness of the DFT + ML model is tested on two classes of organic molecules (TIC10 and hyacinthacines), where the correct isomers were unambiguously assigned to the experimental ones. Overall, the DFT + ML model shows promise for structural assignments in a variety of systems, including stereoisomers, that are often challenging to determine experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

State-to-State Molecular Reactions in the Ultracold Regime (Final Scientific Report)

We have achieved the proposed goal to experimentally probe the AB + CD and AB + C types of reactions with state-to-state resolution, which we also compared to advanced theoretical calculations to help elucidate the role of quantum mechanics in the processes of bond breakage and formation. Our approach uses reactants that are prepared at ultracold temperatures (< 1µK) such that the quantum effects of translational motion are an important factor. Specific example reactions, including the potassium-rubidium metathesis reaction KRb + KRb → K 2 + Rb 2 as well as the atom exchange reaction Rb + KRb → Rb2 + K, are chosen because the technology of quantum internal and motional state control of these types of molecules is particularly advanced. The results for the entire funding period are fruitful. For the majority of this grant, we have constructed a one-of-the-kind quantum degenerate gas apparatus that integrates ion detection and velocity map imaging capabilities, allowing us to explore the KRb + KRb → K 2 + Rb 2 bimolecular reaction in detail. Specifically, we first verified such a reaction indeed proceed at ultracold temperatures by direct detection of reaction products. We then mapped out the complete product state distribution, which was compared to a state-counting model based on statistical theory. Our results show an overall agreement with the statistical state counting model, but also reveal several deviating state-pairs. An exact quantum calculation for molecule-molecule collisions, that is needed to understand these deviations, is however beyond the current state-of-the-art. Beside scrutinizing the reaction products, we also directly observe the reaction intermediate complex, which was quite a surprise to us. The intermediate complexes are long-lived and can interact with the inferred light that we use to trap the ultracold gas. After molecule-molecule collisions, we then explored the more theoretically tractable Rb + KRb reaction, which is endothermic. Surprisingly, we observed an exceedingly long-lived KRb$^*_2$ collisional complexes, with our experimentally measured complex lifetime deviating from conventional theoretical calculations by five orders of magnitude. This discrepancy has motivated many explorations of possible underlying causes, though no model yet captures this phenomenon completely. In the final year and the work that continues today, we extend upon these atom-molecule collision experiments to explore the origin of the long-lived KRb$^*_2$ complex lifetime and develop means to control the outcome of the reaction complex. The 5-year funded work advanced our understanding of chemical reactions at the lowest possible temperatures and at the same time opened up many new questions that are beyond our initial imaginations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insight into Pyroelectricity and Phase Transitions in Ferroelectrics from Nonequilibrium Approach: The Case of PbTiO 3

Computational methodology is proposed to simulate time dependent temperature variations via thermostat approach. Here, the methodology is applied to study pyroelectricity under nonequilibrium conditions in prototypical ferroelectric PbTiO 3 . It is found that at room temperature the nonequilibrium effects in pyroelectricity begin manifesting above 1.0 THz and originate from the overlap with the intrinsic soft mode frequency. Therefore, the proposed computational methodology can be used to predict equilibrium pyroelectric response by choosing a frequency of temperature variation below the intrinsic polar mode frequencies for the material. Probing the dynamics of ferroelectric phase transition with the proposed methodology revealed that it takes about a nanosecond for the computational supercell to undergo phase transition. It is concluded that computational techniques capable of simulating picoseconds can be used to compute pyroelectric response from nonequilibrium dynamical approach while those that can reach nanoseconds are suitable for simulating phase transitions in ferroelectrics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Static and Dynamic Gradient Based Directional Transportation of Neutral Molecules in Swollen Polymer Films

Abstract Materials which selectively transport molecules offer powerful opportunities for concentrating and separating chemical agents. Here, utilizing static and dynamic chemical gradients, transport of molecules within swollen crosslinked polymers is demonstrated. Using an ≈200 μm static hydroxyl to hexyl gradient, the neutral ambipolar nerve agent surrogate diethyl (cyanomethyl)phosphonate (DECP) is directionally transported and concentrated 60‐fold within 4 hours. To accelerate transport kinetics, a dynamic gradient (a “travelling wave”) is utilized. Here, the non‐polar dye pyrene was transported. The dynamic gradient is generated by an ion exchange process triggered by the localized introduction of an aqueous NaCl solution, which converts the gel from hydrophobic to hydrophilic. As the hydrophilic region expands, associated water enters the gel, and pyrene is pushed ahead of the expansion front. The dynamic gradient provides about 10‐fold faster transport kinetics than the static gradient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cobalt Carbonyl Clusters Enable Independent Control of Two Proton Transfer Rates in the Mechanism for Hydrogen Evolution

Abstract Metal carbonyl clusters (MCC's) are atomically defined nanomaterials which can be characterized using the precise tools of molecular (electro)chemistry. HER mechanisms involve two proton transfer (PT) steps in the catalytic cycle. For HER catalyzed by [Co 11 C 2 (CO) 23 ] 3− ( 1 3− ), cyclic voltammetry measurements were used to determine the rate for PT1 as k PT1 =3×10 8 M −1 s −1 , whereas the rate for PT2 is k PT2 =3.7×10 3 M −1 s −1 . The fast, diffusion‐limited rate for PT1 is consistent with a previous report describing [Co 13 C 2 (CO) 24 ] 4− ( 2 4− ), with PT=2.9×10 9 M −1 s −1 . In both cases rate enhancement in PT1 is promoted by the many Co−Co bonds on the surface of the MCC that serve as PT sites: a statistical enhancement of rate akin to the effects of proton relays. In contrast, k PT2 varies and is five orders of magnitude slower for 1 2− compared with 2 4− . Thus, MCC's and nanomaterials offer an opportunity to enhance the rate for PT1 while maintaining thermochemical or kinetic control of PT2.

Pattanayak, Santanu↗

Sequence Design of Random Heteropolymers as Protein Mimics

Random heteropolymers (RHPs) have been computationally designed and experimentally shown to recapitulate protein-like phase behavior and function. However, unlike proteins, RHP sequences are only statistically defined and cannot be sequenced. Recent developments in reversible-deactivation radical polymerization allowed simulated polymer sequences based on the well-established Mayo–Lewis equation to more accurately reflect ground-truth sequences that are experimentally synthesized. This led to opportunities to perform bioinformatics-inspired analysis on simulated sequences to guide the design, synthesis, and interpretation of RHPs. We compared batches on the order of 10000 simulated RHP sequences that vary by synthetically controllable and measurable RHP characteristics such as chemical heterogeneity and average degree of polymerization. Our analysis spans across 3 levels: segments along a single chain, sequences within a batch, and batch-averaged statistics. We discuss simulator fidelity and highlight the importance of robust segment definition. Examples are presented that demonstrate the use of simulated sequence analysis for in-silico iterative design to mimic protein hydrophobic/hydrophilic segment distributions in RHPs and compare RHP and protein sequence segments to explain experimental results of RHPs that mimic protein function. To facilitate the community use of this workflow, the simulator and analysis modules have been made available through an open source toolkit, the RHPapp.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of Charge Sequence Pattern and Lysine-to-Arginine Substitution on the Structural Stability of Bioinspired Polyampholytes

A comprehensive study focusing on the combined influence of the charge sequence pattern and the type of positively charged amino acids on the formation of secondary structures in sequence-specific polyampholytes is presented. The sequences of interest consisting exclusively of ionizable amino acids (lysine, K; arginine, R; and glutamic acid, E) are (EKEK) 5 , (EKKE) 5 , (ERER) 5 , (ERRE) 5 , and (EKER) 5 . The stability of the secondary structure was examined at three pH values in the presence of urea and NaCl. The results presented here underscore the combined prominent effects of the charge sequence pattern and the type of positively charged monomers on secondary structure formation. Additionally, (ERRE) 5 readily aggregated across a wide range of pH. In contrast, sequences with the same charge pattern, (EKKE) 5 , as well as the sequences with the equivalent amino acid content, (ERER) 5 , exhibited no aggregate formation under equivalent pH and concentration conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discovery of Energy Storage Molecular Materials Using Quantum Chemistry-Guided Multiobjective Bayesian Optimization

Redox flow batteries (RFBs) are a promising technology for stationary energy storage applications due to their flexible design, scalability, and low cost. In RFBs, energy is carried in flowable redoxactive materials (redoxmers) which are stored externally and pumped to the cell during operation. Further improvements in the energy density of RFBs necessitates redoxmer designs with wider redox potential windows and higher solubility. Additionally, designing redoxmers with a fluorescence-enabled self-reporting functionality allows monitoring of the state of health of RFBs. To accelerate the discovery of redoxmers with desired properties, state-of-the- art machine learning ( ML) methods, such as multiobjective Bayesian optimization (MBO), are useful. Here, we first employed density functional theory calculations to generate a database of reduction potentials, solvation free energies, and absorption wavelengths for 1400 redoxmer molecules based on a 2,1,3-benzothiadiazole (BzNSN) core structure. From the computed properties, we identified 22 Pareto-optimal molecules that represent best trade-off among all of the desired properties. We further utilized these data to develop and benchmark an MBO approach to identify candidates quickly and efficiently with multiple targeted properties. With MBO, optimal candidates from the 1400-molecule data set can be identified at least 15 times more efficiently compared to the brute force or random selection approach. Importantly, we utilized this approach for discovering promising redoxmers from an unseen database of 1 million BzNSN-based molecules, where we discovered 16 new Pareto-optimal molecules with significant improvements in properties over the initial 1400 molecules. Furthermore, we anticipate that this active learning technique is general and can be utilized for the discovery of any class of functional materials that satisfies multiple desired property criteria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Guide to Nonaqueous Electrochemistry of f-Element Complexes

Electrochemistry is a powerful tool for assessing and understanding the redox chemistry of molecular complexes. Cyclic voltammetry enables the f-element community to study molecules in unusually high or low oxidation states, which potentially pose important broad-scope questions of electronic structure. In the pursuit of boundary-pushing compounds, reactive air- and moisture-sensitive species are often encountered, which can be challenging to characterize, especially when they are chemically incompatible with certain solvents, electrolytes, or electrodes or when their potentials lie outside of common electrochemical windows. Nonaqueous solvents and pseudo reference electrodes complicate many of the standard practices in acquiring high-quality and reproducible electrochemical data. This guide presents a detailed discussion of selecting appropriate cell conditions and referencing and addresses metrics for evaluating electrochemical and chemical reversibility. These methodological approaches have been extended to best practices for the electrochemical analysis of radioactive transuranic complexes.

Electrodes↗

Accurate Machine Learning for Predicting the Viscosities of Deep Eutectic Solvents

Deep eutectic solvents (DESs) are emerging as environmentally friendly designer solvents for mass transport and heat transfer processes in industrial applications; however, the lack of accurate tools to predict and thus control their viscosities under both a range of environmental factors and formulations hinders their general application. While DESs may serve as designer solvents, with nearly unlimited combinations, this unfortunately makes it experimentally infeasible to comprehensively measure the viscosities of all DESs of potential industrial interest. To assist in the design of DESs, we have developed several new machine learning (ML) models that accurately and rapidly predict the viscosities of a diverse group of DESs at different temperatures and molar ratios using, to date, one of the most comprehensive data sets containing the properties of over 670 DESs over a wide range of temperatures (278.15–385.25 K). Three ML models, including support vector regression (SVR), feed forward neural networks (FFNNs), and categorical boosting (CatBoost), were developed to predict DES viscosity as a function of temperature and molar ratio and contrasted with multilinear and two-factor polynomial regression baselines. Further, quantum chemistry-based, COSMO-RS-derived sigma profile (σ-profile) features were used as inputs for the ML models. The CatBoost model is excellent at externally predicting DES viscosity, as indicated by high R 2 (0.99) and low root-mean-square-error (RMSE) and average absolute relative deviations (AARD) (5.22%) values for the testing data sets, and 98% of the data points lie within the 15% of AARD deviations. Furthermore, SHapley additive explanation (SHAP) analysis was employed to interpret the ML results and rationalize the viscosity predictions. The result is an ML approach that accurately predicts viscosity and will aid in accelerating the design of appropriate DESs for industrial applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

COCOMO2: A Coarse-Grained Model for Interacting Folded and Disordered Proteins

Biomolecular interactions are essential in many biological processes, including complex formation and phase separation processes. Coarse-grained computational models are especially valuable for studying such processes via simulation. Here, we present COCOMO2, an updated residue-based coarse-grained model that extends its applicability from intrinsically disordered peptides to folded proteins. This is accomplished with the introduction of a surface exposure scaling factor, which adjusts interaction strengths based on solvent accessibility, to enable the more realistic modeling of interactions involving folded domains without additional computational costs. COCOMO2 was parametrized directly with solubility and phase separation data to improve its performance on predicting concentration-dependent phase separation for a broader range of biomolecular systems compared to the original version. COCOMO2 enables new applications including the study of condensates that involve IDPs together with folded domains and the study of complex assembly processes. COCOMO2 also provides an expanded foundation for the development of multiscale approaches for modeling biomolecular interactions that span from residue-level to atomistic resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vibrational Quantum-State-Controlled Reactivity in the O 2 + + C 3 H 4 Reaction

Quantum-state-controlled reactivity is a long-standing goal in the field of physical chemistry. In this work, we explore the vibrational-state-dependent behavior of the ion–molecule reaction between O 2 + in distinct vibrational states and two isomers of C 3 H 4 , allene (H 2 C 3 H 2 ) and propyne (H 3 C 3 H). While most products are formed regardless of the vibrational state of O 2 + , the branching ratios are influenced by vibrational excitation, and a new product, C 2 O + , appears exclusively in the excited-state reactions. This selective formation of C 2 O + demonstrates that vibrational excitation can effectively activate a reaction pathway, providing direct evidence of quantum-state control in the reactivity. These results represent an important step toward the goal of quantum-state-controlled chemistry in molecular systems.

Chemical reactions↗

Amphiphilic Baskets for Supramolecular Nanoarchitectures at Interfaces: Inverted Monolayer Formation on Water

Interfacial chemistry of molecular baskets remains poorly understood despite their promise for supramolecular applications of detection and sequestration of toxic molecules including those of illicit drugs, organophosphorus compounds, and anticancer agents. We present a fundamental investigation of the interfacial behavior of three amphiphilic supramolecular baskets (ASB 4, 8, and 12), having increasingly longer yet linear alkyl chains at the top of their bowl-shaped cavity. The studies were completed at the air−water interface to elucidate surface activity, interfacial stability, self-assembly, and monolayer organization that drive inverted monolayer formation, in which the molecular arms orient toward the aqueous phase in a configuration opposite to that typically observed for lipids. Herein, surface pressure−area isotherms of ASB 4, 8, 12, deposited on a water surface, were performed in tandem with nonequilibrium relaxation experiments to quantify surface activity, thermodynamic stability, and monolayer compressibility of the baskets’ monolayer assembly. Brewster angle microscopy enabled direct visualization of morphological evolution, aggregation, and packing at the interface. We show that systematic extension of the hydrocarbon arms, from four to 12 methylene groups, progressively modifies intermolecular packing, drives distinct two-dimensional aggregation pathways, and increases number densities at the air−water interface. Atomistic molecular dynamics simulations corroborate many of these experimentally observed trends and provide mechanistic detail on the cooperative roles of basket topology and interfacial concentration in regulating the hydration structure and dynamics within the cavities generated by surface-adsorbing baskets, consistent with observed variations in surface activity and packing. Our results establish how the topology of these unique supramolecules and their concentration govern interfacial organization and offer a rational framework for designing amphiphiles with predictable behavior at soft interfaces.

Hydrocarbons↗