SEARCH · Search NASA
Results for “Molecular Structure”
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.
HamLib: A library of Hamiltonians for benchmarking quantum algorithms and hardware
In order to characterize and benchmark computational hardware, software, and algorithms, it is essential to have many problem instances on-hand. This is no less true for quantum computation, where a large collection of real-world problem instances would allow for benchmarking studies that in turn help to improve both algorithms and hardware designs. To this end, here we present a large dataset of qubit-based quantum Hamiltonians. The dataset, called HamLib (for Hamiltonian Library), is freely available online and contains problem sizes ranging from 2 to 1000 qubits. HamLib includes problem instances of the Heisenberg model, Fermi-Hubbard model, Bose-Hubbard model, molecular electronic structure, molecular vibrational structure, MaxCut, Max- k -SAT, Max- k -Cut, QMaxCut, and the traveling salesperson problem. The goals of this effort are (a) to save researchers time by eliminating the need to prepare problem instances and map them to qubit representations, (b) to allow for more thorough tests of new algorithms and hardware, and (c) to allow for reproducibility and standardization across research studies.
Effect of layer bending on montmorillonite hydration and structure from molecular simulation
Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm 2 system that included dehydrated and hydrated regions in the same interlayer.
Peak2Patch: High-Fidelity Functional Group Identification through Attention-Based Fusion of Infrared and Mass Spectra
Identifying molecular structure based on spectroscopic readings is a key task in a variety of chemical and biological applications. Common spectroscopy techniques, such as Infrared (IR) Spectroscopy and Mass Spectrometry (MS), provide detailed information on the structure of molecular compounds but nonetheless require expert-level knowledge to decode. Machine learning has emerged as a potential solution for automating structure prediction from chemical spectra; however, current approaches generally focus on single sensor modalities, neglecting to leverage the complementary information contained within differing spectra. In this paper, we introduce Peak2Patch, a novel approach to fusion-enhanced prediction of functional groups from IR and mass spectra. First, we perform a detailed comparison of backbone networks for encoding both sparse mass spectra and dense IR spectra and demonstrate the superior performance of transformer neural networks over current state-of-the-art convolutional neural networks. Second, we evaluate three broad categories of fusion: early (raw feature), middle (deep feature), and late (decision) fusion, demonstrating the potential of a deep feature fusion-based approach. Lastly, we present Peak2Patch, our attention-based fusion scheme, which leverages cross-attention to mix features between encoded tokens of the two modalities. We validate our approach on a publicly available multimodal spectroscopic data set of 790k simulated molecules, demonstrating a large improvement in functional group prediction over both the previous state-of-the-art and our own strong single-modal baselines.
Beyond real: alternative unitary cluster Jastrow models for molecular electronic structure calculations on near-term quantum computers
Near-term quantum devices require wavefunction ansätze that are expressive while also of shallow circuit depth in order to both accurately and efficiently simulate molecular electronic structure. While the unitary coupled cluster ansatz (e.g., UCCSD) has become a standard, the high gate count associated with the implementation of this limits its feasibility on noisy intermediate-scale quantum (NISQ) hardware. k -Fold unitary cluster Jastrow (uCJ) ansätze mitigate this challenge by providing O( kN 2 ) circuit scaling and favorable linear depth circuit implementation. Previous work has focused on the real orbitalrotation (Re-uCJ) variant of uCJ, which allows an exact (Trotter-free) implementation. Here we extend and generalize the k -fold uCJ framework by introducing two new variants, Im-uCJ and g-uCJ, which incorporate imaginary and fully complex orbital rotation operators, respectively. Similar to Re-uCJ, both of the new variants achieve quadratic gate-count scaling. Our results focus on the simplest k = 1 model, and show that the uCJ models frequently maintain energy errors within chemical accuracy (∼1 kcal mol −1 ). Both g-uCJ and Im-uCJ are more expressive in terms of capturing electron correlation and are also more accurate than the earlier Re-uCJ ansatz. We further show that Im-uCJ and g-uCJ circuits can also be implemented exactly, without any Trotter decomposition. Numerical tests using k = 1 on H 2 , H 3 + , Be 2 , C 2 H 4 , C 2 H 6 and C 6 H 6 in various basis sets confirm the practical feasibility of these shallow Jastrow-based ansätze for applications on near-term quantum hardware.
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.
Molecular Design Principles for Photosystem I-Based Biohybrid Solar Fuel Catalysts
Direct solar-to-chemical conversion offers a compelling route to clean, dispatchable energy. Photosystem I (PSI), an evolutionarily optimized light-driven oxidoreductase, can be repurposed for solar-fuel production by coupling its photochemistry to catalytic interfaces. However, the molecular determinants that govern productive electron transfer to abiotic catalysts remain poorly understood. Here, we present molecular structures of active PSI-Pt nanoparticle (PtNP) biohybrids that reveal how protein architecture controls catalyst access, binding geometry, and photocatalytic efficiency. Removal of stromal subunits exposes the electron transfer chain and enables PtNP binding proximal to the F X cluster, demonstrating that steric occlusion limits access to native acceptor regions in PSI. In contrast, in trimeric PSI, PtNPs bind at multiple sites per monomer, but only a subset are positioned within electron transfer distance of terminal cofactors, resulting in a heterogeneous population of productive and nonproductive configurations. Structural analyses and molecular dynamics simulations define the interface topology, electrostatics, and cofactor-to-nanoparticle distances that govern catalyst binding and electron transfer. These results establish that catalytic inefficiency arises not only from intrinsic electron transfer constraints but also from the distribution of binding geometries imposed by the protein scaffold. Together, these findings provide a molecular framework linking protein structure to biohybrid function and define design principles for engineering PSI-based solar fuel systems and protein-nanomaterial interfaces for light-driven catalysis.
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.
Chemical Diversity of Oligomers in Biomass Fast Pyrolysis Oils, Part 2: Heavy Lignin-Derived Molecules and Highly Dehydrated Sugars from Dichloromethane-Insoluble Pyrolytic Lignin
Here, this paper, the second in a series, investigates the water-insoluble, dichloromethane (DCM)-insoluble residue of BTG pyrolysis oil, also known in the literature as high-molecular-weight pyrolytic lignin (HMW-PL). HMW-PL accounts for 9.2 wt % of the original bio-oil, which is known to increase upon aging and contribute to coke formation during bio-oil upgrading. Understanding its composition is crucial for managing pyrolysis oil during the processing. The HMW-PL was fractionated using a silica gel column with solvents of increasing polarity: ethyl acetate (EA), acetone (AC), isopropanol (ISO), and methanol (MeOH), yielding 49.2, 19.9, 8.2, and 8.1 wt %, respectively. Each subfraction was characterized by UV-fluorescence, FTIR, HSQC NMR, and (−)APCI-FT-Orbitrap MS. Analysis revealed that the EA subfraction primarily contained aromatic dimers (C 19 –C 20 ) and trimers (C 17 –C 13 ), while the MeOH subfraction was rich in oxygenated aliphatic monomers (C 8 –C 16 ) likely derived from sugar dehydration products. The resulting fractions were further separated via preparatory HPLC. Fifty-five candidate molecular structures were proposed based on the Orbitrap MS data, supported by UV-fluorescence, FTIR, and NMR results. This fractionation strategy defined four distinct subfractions, enabling the proposal of surrogate molecular structures and advancing the molecular-level understanding of the HMW-PL fraction in the pyrolysis bio-oil.
Analytic Nuclear Gradients Including Oriented External Electric Fields in a Molecule-Fixed Frame
Electric-field-assisted chemistry has attracted much attention in recent years, particularly in the context of oriented external electric fields for controlling molecular structure and reactivity. Such fields have been explored in a wide range of applications, including switching materials, nanoparticles, controllable catalysts, medicines, and clinical therapies. However, the determination of fixed fields in the laboratory frame becomes ineffective for flexible molecules, as conformational changes can significantly alter the relative orientation between the applied field and molecular structure. In this work, we propose two molecular reference frames─the principal axis frame and the local reference frame─to define oriented electric fields within the molecular framework. These coordinate systems powerfully eliminate ambiguities in the relative orientation between the applied field and the molecule. Analytic nuclear gradients in the presence of external electric fields are derived and implemented, with an initial application to field-dependent geometry optimizations of cis - and trans -formanilide. Analysis of the resulting field-induced equilibrium structures reveals distinct structural responses, validating the accuracy and robustness of the proposed formalism. The analytic gradient framework enables systematic investigations of molecular properties and reactivity under arbitrarily oriented electric fields, opening new opportunities for computational modeling and rational design in electric-field-controlled chemistry.
Influence of thermal treatment on structure and catalytic performance of ceria-zirconia supported copper oxide (CuO x /Ce y Zr 1-y O 2 ) catalysts for CO oxidation
Copper oxide (CuO x ) supported on ceria-zirconia (Ce y Zr 1-y O 2 , y = 1.0, 0.5, 0.0) catalysts were investigated to elucidate the effects of thermal treatment on their physicochemical properties and catalytic performance in carbon monoxide (CO) oxidation. Here, the catalysts were synthesized via a one-pot chemical vapor deposition (OP-CVD) method at 700˚C and 900˚C with controlled Cu loading. Characterization techniques, including synchrotron X-ray diffraction (S-XRD), Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), inductively coupled plasma spectroscopy (ICP), and N 2 adsorption-desorption, were implemented to probe the crystalline structure, molecular and electronic structure, oxygen vacancies, specific surface area (SSA) and metal loading. CO oxidation was chosen as a model reaction to explore the structure-catalytic performance relationship. A ∼100% CO conversion was achieved at < 150˚C, particularly with the CuO x /CeO 2 catalyst calcined at 700˚C. In contrast, calcination at 900˚C caused a ∼90% decrease in SSA and a ∼24% increase in T 50 . Activity tests revealed that increasing ZrO 2 content lowered CO oxidation activity despite generating more defect sites. In-situ measurement of the 700 °C calcined samples revealed the presence of stable and unstable defects in CuO x /Ce 0.5 Zr 0.5 O 2 and CeO 2 respectively, which play a key role in the activity of the catalysts. The results highlight that catalytic performance is closely related to the SSA. Furthermore, an optimum calcination temperature favor significant oxygen vacancy formation with required CuO x -support interactions, enhancing redox properties and catalytic performance.
Marine Algae Polysaccharides: An Overview of Characterization Techniques for Structural and Molecular Elucidation
Polysaccharides make up a large portion of the organic material from and in marine organisms. However, their structural characterization is often overlooked due to their complexity. With many high-value applications and unique bioactivities resulting from the polysaccharides’ complex and heterogeneous structures, dedicated analytical efforts become important to achieve structural elucidation. Because algae represent the largest marine resource of polysaccharides, the majority of the discussion is focused on well-known algae-based hydrocolloid polymers. The native environment of marine polysaccharides presents challenges to many conventional analytical techniques necessitating novel methodologies. We aim to deliver a review of the current state of the art in polysaccharide characterization, focused on capabilities as well as limitations in the context of marine environments. This review covers the extraction and isolation of marine polysaccharides, in addition to characterizations from monosaccharides to secondary and tertiary structures, highlighting a suite of analytical techniques.
Electronic structure and magnetic tendencies of trilayer La 4 Ni 3 O 10 under pressure: Structural transition, molecular orbitals, and layer differentiation
Motivated by the recent observation of superconductivity in the pressurized trilayer Ruddlesden-Popper (RP) nickelate La 4 Ni 3 O 10 , we explore its structural, electronic, and magnetic properties as a function of hydrostatic pressure from first-principles calculations. We find that an orthorhombic (monoclinic)-to-tetragonal transition under pressure takes place concomitantly with the onset of superconductivity. The electronic structure of La 4 Ni 3 O 10 can be understood using a molecular trimer basis wherein n molecular subbands arise as the d z 2 orbitals hybridize strongly along the c axis within the trilayer. The magnetic tendencies indicate that the ground state at ambient pressure is formed by nonmagnetic inner planes and stripe-ordered outer planes that are antiferromagnetically coupled along the c axis, resulting in an unusual ↑, 0, ↓ stacking that is consistent with the spin density wave model previously suggested by neutron diffraction. Such a state is destabilized at the pressure where superconductivity arises. Despite the presence of d z 2 states at the Fermi level, the d x 2 –y 2 orbitals also play a key role in the electronic structure of La4Ni3O10. Finally, this active role of the d x 2 –y 2 states in the low-energy physics of the trilayer RP nickelate, together with the distinct electronic behavior of the inner and outer planes, resembles the physics of multilayer cuprates.
Twins in rotational spectroscopy: Does a rotational spectrum uniquely identify a molecule?
Rotational spectroscopy is the most accurate method for determining structures of molecules in the gas phase. It is often assumed that a rotational spectrum is a unique “fingerprint” of a molecule. The availability of large molecular databases and the development of artificial intelligence methods for spectroscopy make the testing of this assumption timely. In this paper, we pose the determination of molecular structures from rotational spectra as an inverse problem. Within this framework, we adopt a funnel-based approach to search for molecular twins, which are two or more molecules, which have similar rotational spectra but distinctly different molecular structures. Here we demonstrate that there are twins within standard levels of computational accuracy by generating rotational constants for many molecules from several large molecular databases, indicating that the inverse problem is ill-posed. However, some twins can be distinguished by increasing the accuracy of the theoretical methods or by performing additional experiments.
Dashboard for Visualizing Molecular Property Prediction Machine Learning Results
This is a dashboard for exploring the results of machine learning models for predicting molecular properties from molecular structure. It includes tools for: 1. Modifying molecules to observe the change in predicted properties 2. Exploring the relationship between molecular structure and predicted properties 3. Recommending structurally similar molecules with improved properties 4. Exploring the impact of data subsampling on model performance metrics
Navigating Large Chemical Spaces Using Graph Theory and Integer Programming
Navigating and analyzing large chemical spaces are necessary to accelerate the design and discovery of new molecules and chemical processes. In this work, we introduce a computational framework that integrates graph theory and integer programming to enable the efficient navigation of large chemical spaces. Our framework represents the chemical space as a graph, wherein nodes represent molecules and edges represent the degree of similarity or connectivity based on domain-specific information. Using the graph representation, we identify representative molecules by computing the so-called minimum dominating set (MDS), which in our context is the minimum set of molecules that is connected to all other molecules. We present a suite of solution strategies for the MDS problem including heuristic and rigorous integer programming (IP) approaches. We show that these approaches allow us to capture physicochemical properties and domain-specific logic and constraints, facilitating the identification of molecules with the target properties. We demonstrate the effectiveness of the proposed approach by navigating the chemical space of per- and polyfluoroalkyl substances (PFAS); this comprises approximately 15,000 molecular structures. We compare our framework against traditional dimensionality reduction and clustering methods such as t-SNE and K-means clustering.
Two datasets are better than one: method of double moments for 3D reconstruction in cryo-EM
Cryo-electron microscopy is a powerful imaging technique for reconstructing three-dimensional molecular structures from noisy tomographic projection images of randomly oriented particles. We introduce a new data fusion framework, termed the method of double moments, which reconstructs molecular structures from two instances of the second-order moment of projection images obtained under distinct orientation distributions: one uniform, the other non-uniform and unknown. We prove that these moments generically uniquely determine the underlying structure, up to a global rotation and reflection, and we develop a convex-relaxation-based algorithm that achieves accurate recovery using only second-order statistics. Our results demonstrate the advantage of collecting and modeling multiple datasets under different experimental conditions, illustrating that leveraging dataset diversity can substantially enhance reconstruction quality in computational imaging tasks.
Molecular and structural characterization of a Bacillus cereus strain producing an anthrax-like capsule
Bacillus cereus is a ubiquitous Gram-positive, spore-forming, rod-shaped saprophytic bacterium, occasionally reported to cause food-borne illnesses. However, instances of B. cereus strains harboring anthrax toxin and capsule genes have elevated certain strains as formidable pathogens and biothreats. This study focuses on the genomic analysis and the structural characterization of capsular material produced by the virulent B. cereus PATH2418 strain, isolated from the wound of a traumatic open fracture patient. The genome was sequenced using Nanopore MinION sequencing, revealing a chromosome of 5,270,283 bp and three plasmids. One plasmid, pATH1, was found to encode an operon for the biosynthesis of a bacterial capsule. This operon had sequence homology to the Bacillus anthracis capBCADE operon, which encodes the poly-γ-D-glutamate (PDGA) capsule. The capsule production in B. cereus PATH2418 was influenced by temperature and CO 2 levels. Structural analysis of the capsular material using a combined approach of nuclear magnetic resonance (NMR) and high-performance liquid chromatography (HPLC) techniques confirmed the presence of a high-molecular-weight poly-γ-glutamate capsule, with an enantiomeric composition of approximately 67% D-glutamic acid and 33% L-glutamic acid, matching that of B. anthracis.