Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular conformation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

CoarsenConf: Equivariant Coarsening with Aggregated Attention for Molecular Conformer Generation

Molecular conformer generation (MCG) is an important task in cheminformatics and drug discovery. The ability to efficiently generate low-energy 3D structures can avoid expensive quantum mechanical simulations, leading to accelerated virtual screenings and enhanced structural exploration. Several generative models have been developed for MCG, but many struggle to consistently produce high-quality conformers for meaningful downstream applications. To address these issues, we introduce CoarsenConf, which coarse-grains molecular graphs based on torsional angles and integrates them into an SE(3)-equivariant hierarchical variational autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained atomic coordinates of subgraphs connected via rotatable bonds, creating a variable-length coarse-grained latent representation. Our model uses a novel aggregated attention mechanism to restore fine-grained coordinates from the coarse-grained latent representation, enabling efficient generation of accurate conformers. Furthermore, we evaluate the chemical and biochemical quality of our generated conformers on multiple downstream applications, including property prediction and large-scale oracle-based protein docking. Overall, CoarsenConf generates more accurate conformer ensembles compared to prior generative models.

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

Effect of Molecular Conformation on Intermolecular Interactions and Photovoltaic Performances of Giant Molecule Acceptors

The molecular conformation of giant molecule acceptors (GMAs) plays a significant role in regulating the intermolecular interactions and their photovoltaic performances in organic solar cells (OSCs). For the linear GMA GT‐l, the stronger homo‐molecular interaction causes its aggregation being weakly affected by the donor, thus forming an ordered molecular stacking and proper phase separation in its blend film. The star‐shaped GMA GT‐s‐based blend film shows a dominant hetero‐molecular interaction that suppresses the aggregation of the donor and acceptor, resulting in smaller phase separation and more uniform vertical phase distribution. While for another star‐shaped GMA GTs, the weakest hetero‐molecular interaction causes its blend film to form larger phase separation. Therefore, the GT‐l based OSC with PM6 as donor shows the highest charge mobilities, the fastest charge transfer (CT) process, reduced energy loss and less charge recombination, contributing to a higher power conversion efficiency (PCE) of 19.03%. Comparatively, the PCEs of the OSCs based on GTs and GT‐s are 18.05% and 17.58% respectively. Notably, all the three GMAs based OSCs show excellent thermal stability and long‐term storage stability. This study provides a facile strategy by tuning the linking unit and its connecting mode for designing highly efficient and stable organic photovoltaic materials.

Chemistry↗
DOE OSTI2023

Similarity Downselection: Finding the n Most Dissimilar Molecular Conformers for Reference-Free Metabolomics

Computational methods for creating in silico libraries of molecular descriptors (e.g., collision cross sections) are becoming increasingly prevalent due to the limited number of authentic reference materials available for traditional library building. These so-called “reference-free metabolomics” methods require sampling sets of molecular conformers in order to produce high accuracy property predictions. Due to the computational cost of the subsequent calculations for each conformer, there is a need to sample the most relevant subset and avoid repeating calculations on conformers that are nearly identical. The goal of this study is to introduce a heuristic method of finding the most dissimilar conformers from a larger population in order to help speed up reference-free calculation methods and maintain a high property prediction accuracy. Finding the set of the n items most dissimilar from each other out of a larger population becomes increasingly difficult and computationally expensive as either n or the population size grows large. Because there exists a pairwise relationship between each item and all other items in the population, finding the set of the n most dissimilar items is different than simply sorting an array of numbers. For instance, if you have a set of the most dissimilar n = 4 items, one or more of the items from n = 4 might not be in the set n = 5. An exact solution would have to search all possible combinations of size n in the population exhaustively. We present an open-source software called similarity downselection (SDS), written in Python and freely available on GitHub. SDS implements a heuristic algorithm for quickly finding the approximate set(s) of the n most dissimilar items. We benchmark SDS against a Monte Carlo method, which attempts to find the exact solution through repeated random sampling. We show that for SDS to find the set of n most dissimilar conformers, our method is not only orders of magnitude faster, but it is also more accurate than running Monte Carlo for 1,000,000 iterations, each searching for set sizes n = 3–7 out of a population of 50,000. We also benchmark SDS against the exact solution for example small populations, showing that SDS produces a solution close to the exact solution in these instances. Using theoretical approaches, we also demonstrate the constraints of the greedy algorithm and its efficacy as a ratio to the exact solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Hyperconjugation-controlled molecular conformation weakens lithium-ion solvation and stabilizes lithium metal anodes

Tuning the solvation structure of lithium ions via electrolyte engineering has proven effective for lithium metal (Li) anodes. Further advancement that bypasses the trial-and-error practice relies on the establishment of molecular design principles. Expanding the scope of our previous work on solvent fluorination, we report here an alternative design principle for non-fluorinated solvents, which potentially have reduced cost, environmental impact, and toxicity. By studying non-fluorinated ethers systematically, we found that the short-chain acetals favor the [gauche, gauche] molecular conformation due to hyperconjugation, which leads to weakened monodentate coordination with Li + . The dimethoxymethane electrolyte showed fast activation to >99% coulombic efficiency (CE) and high ionic conductivity of 8.03 mS cm -1 . The electrolyte performance was demonstrated in anode-free Cu$∥$LFP pouch cells at current densities up to 4 mA cm -2 (70 to 100 cycles) and thin-Li$∥$high-loading-LFP coin cells (200–300 cycles). Overall, we demonstrated and rationalized the improvement in Li metal cyclability by the acetal structure compared to ethylene glycol ethers. We expect further improvement in performance by tuning the acetal structure.

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

Towards interpretable Cryo-EM: disentangling latent spaces of molecular conformations

Molecules are essential building blocks of life and their different conformations (i.e., shapes) crucially determine the functional role that they play in living organisms. Cryogenic Electron Microscopy (cryo-EM) allows for acquisition of large image datasets of individual molecules. Recent advances in computational cryo-EM have made it possible to learn latent variable models of conformation landscapes. However, interpreting these latent spaces remains a challenge as their individual dimensions are often arbitrary. The key message of our work is that this interpretation challenge can be viewed as an Independent Component Analysis (ICA) problem where we seek models that have the property of identifiability. That means, they have an essentially unique solution, representing a conformational latent space that separates the different degrees of freedom a molecule is equipped with in nature. Thus, we aim to advance the computational field of cryo-EM beyond visualizations as we connect it with the theoretical framework of (nonlinear) ICA and discuss the need for identifiable models, improved metrics, and benchmarks. Moving forward, we propose future directions for enhancing the disentanglement of latent spaces in cryo-EM, refining evaluation metrics and exploring techniques that leverage physics-based decoders of biomolecular systems. Moreover, we discuss how future technological developments in time-resolved single particle imaging may enable the application of nonlinear ICA models that can discover the true conformation changes of molecules in nature. The pursuit of interpretable conformational latent spaces will empower researchers to unravel complex biological processes and facilitate targeted interventions. This has significant implications for drug discovery and structural biology more broadly. More generally, latent variable models are deployed widely across many scientific disciplines. Thus, the argument we present in this work has much broader applications in AI for science if we want to move from impressive nonlinear neural network models to mathematically grounded methods that can help us learn something new about nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

True molecular conformation and structure determination by three-dimensional electron diffraction of PAH by-products potentially useful for electronic applications

The true mol­ecular conformation and the crystal structure of benzo[e]di­naphtho­[2,3-a;1',2',3',4'-ghi]fluoranthene, 7,14-di­phenyl­naphtho­[1,2,3,4-cde]bis­anthene and 7,16-di­phenyl­naphtho­[1,2,3,4-cde]heli­anthrene were determined ab initio by 3D electron diffraction. All three mol­ecules are remarkable polycyclic aromatic hydro­carbons. The mol­ecular conformation of two of these com­pounds could not be determined via classical spectroscopic methods due to the large size of the mol­ecule and the occurrence of multiple and reciprocally connected aromatic rings. The mol­ecular structure of the third mol­ecule was previously considered provisional. These com­pounds were isolated as by-products in the synthesis of similar products and were at the same time nanocrystalline and available only in very limited amounts. 3D electron diffraction data, taken from submicrometric single crystals, allowed for direct ab initio structure solution and the unbiased determination of the inter­nal mol­ecular conformation. Detailed synthetic routes and spectroscopic analyses are also discussed. Based on many-body perturbation theory simulations, benzo[e]di­naphtho­[2,3-a;1',2',3',4'-ghi]fluoranthene may be a promising candidate for triplet–triplet annihilation and 7,14-di­phenyl­naphtho­[1,2,3,4-cde]bis­anthene may be a promising candidate for inter­molecular singlet fission in the solid state.

36 MATERIALS SCIENCE↗
NASA NTRS1993

The molecular matching problem

Molecular chemistry contains many difficult optimization problems that have begun to attract the attention of optimizers in the Operations Research community. Problems including protein folding, molecular conformation, molecular similarity, and molecular matching have been addressed. Minimum energy conformations for simple molecular structures such as water clusters, Lennard-Jones microclusters, and short polypeptides have dominated the literature to date. However, a variety of interesting problems exist and we focus here on a molecular structure matching (MSM) problem.

Kincaid, Rex K.↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Engineering Structural Transitions in a Multilevel Molecular Switch via Intermolecular Coupling

Controlling molecular conformations with atomic precision is essential for advancing molecular functional electronics, as well as our understanding of molecular dynamics. While switching between bistable molecular conformers is common in nature, creating systems with multiple, addressable states remains synthetically challenging. Here, we demonstrate a bottom-up strategy in which intermolecular interactions give rise to multilevel functionality within a simple two-molecule assembly. Using low-temperature scanning tunneling microscopy, we show that a pyrrolidine dimer on Cu(100) exhibits six distinct adsorption conformations, exceeding the four expected from two independent bistable units. This unusual complexity arises from the interplay between intermolecular van der Waals attraction and steric repulsion, which reshapes the potential energy landscape and changes a single high-energy transition into a sequential two-step pathway. Each step is driven by low-energy inelastic electron excitations, achieving a switching efficiency an order of magnitude higher than that of the monomer. Here, by tuning the bias voltage and tip–molecule distance, we achieve deterministic control over multiple stable states, establishing a general design principle for on-demand engineering of collective molecular behavior and energy-efficient multilevel molecular devices.

Molecular interactions↗
DOE OSTI2025

Aromatic Amide Foldamers Show Conformation‐Dependent Electronic Properties

Electron transport in organic molecules and biomolecules is governed by electronic structure and molecular conformations. Despite recent progress, key challenges remain in understanding the role of intramolecular interactions and three-dimensional (3D) conformations on the electron transport behavior of organic molecules. In this work, the electronic properties of aromatic amide foldamers are characterized that organize into distinct 3D structures, including an extended secondary amide that adopts a trans-conformation and a folded N-methylated tertiary amide that adopts a cis-conformation. Results from single-molecule electronic experiments show that the extended secondary amide exhibits a fourfold enhancement in molecular conductance compared to the folded N-methylated tertiary amide, despite a longer contour length. The results show that extended amide molecules are governed by a through-bond electron transport mechanism, whereas folded amide molecules are dominated by through-space transport. Bulk spectroscopic characterization and density functional theory calculations further reveal that extended amides have a smaller HOMO–LUMO gap and larger transmission values compared to folded amides, consistent with single-molecule electronic experiments. Overall, this work shows that 3D molecular conformations significantly influence the electronic properties of single-molecule junctions.

electron transport↗
DOE OSTI2024

Shape-persistent ladder molecules exhibit nanogap-independent conductance in single-molecule junctions

Molecular electronic devices require precise control over the flow of current in single molecules. However, the electron transport properties of single molecules critically depend on dynamic molecular conformations in nanoscale junctions. Here, in this work, we report a unique strategy for controlling molecular conductance using shape-persistent molecules. Chemically diverse, charged ladder molecules, synthesized via a one-pot multicomponent ladderization strategy, show a molecular conductance (d[log( G/G 0 )]/d x ≈ -0.1 nm -1 ) that is nearly independent of junction displacement, in stark contrast to the nanogap-dependent conductance (d[log( G/G 0 )]/d x ≈ -7 nm -1 ) observed for non-ladder analogues. Ladder molecules show an unusually narrow distribution of molecular conductance during dynamic junction displacement, which is attributed to the shape-persistent backbone and restricted rotation of terminal anchor groups. These principles are further extended to a butterfly-like molecule, thereby demonstrating the strategy's generality for achieving gap-independent conductance. Overall, our work provides important avenues for controlling molecular conductance using shape-persistent molecules. Achieving robust and controllable conductance in single-molecule junctions is challenging due to the dynamic nature of molecular conformations that fluctuate over operational timescales. A strategy using shape-persistent molecules has now been developed that demonstrates nearly junction-displacement-independent conductance, providing a stable solution for single-molecule electronic properties.

molecular electronics↗
DOE OSTI2023

Molecular Structure and Conformational Design of Donor‐Acceptor Conjugated Polymers to Enable Predictable Optoelectronic Property

Abstract Tuning the optoelectronic properties of donor‐acceptor conjugated polymers (D‐A CPs) is of great importance in designing various organic optoelectronic devices. However, there remains a critical challenge in precise control of bandgap through synthetic approach, since the chain conformation also affects molecular orbital energy levels. Here, D‐A CPs with different acceptor units are explored that show an opposite trend in energy band gaps with the increasing length of oligothiophene donor units. By investigating their chain conformation and molecular orbital energy, it is found that the molecular orbital energy alignment between donor and acceptor units plays a crucial role in dictating the final optical bandgap of D‐A CPs. For polymers with staggered orbital energy alignment, the higher HOMO with increasing oligothiophene length leads to a narrowing of the optical bandgap despite decreased chain rigidity. On the other hand, for polymers with sandwiched orbital energy alignment, the increased band gap with increasing oligothiophene length originates from the reduction of bandwidth due to more localized charge density distribution. Thus, this work provides a molecular understanding of the role of backbone building blocks on the chain conformation and bandgaps of D‐A CPs for organic optoelectronic devices through the conformation design and segment orbital energy alignment.

Cao, Zhiqiang↗
DOE OSTI2026

The Molecular Structures of Liquid and Glassy Nifedipine and Felodipine and Their Incorporation into PVP

Background: Amorphous drug formulations are commonly used to improve the solubility and bioavailability of poorly soluble molecular pharmaceuticals, yet less is known about their molecular conformations and local bonding interactions than their crystalline phases. Methods: High-energy X-ray diffraction structure factor measurements have been made on liquid and glassy nifedipine (NIF), felodipine (FEL), NIF 1:3 polyvinylpyrrolidone (PVP), and FEL 1:3 PVP wt.% mixtures. The corresponding X-ray pair distribution functions have been interpreted using empirical potential structure refinement using different models and density functional theory conformer calculations. Results: In both NIF and FEL, the NH···O inter-molecular hydrogen bonds between the pyridyl nitrogen and ester carbonyls are found to be considerably weaker than those observed in the crystalline polymorphs. For nifedipine, it is proposed that either inter-molecular NH…ON nitro bonds are present and/or a fraction (<20%) of conformational changes, with the aryl ring flipped, occur in the liquid state. For felodipine, the models indicate significant disorder associated with the methyl and ethyl side chains in the liquid state, with the main peak intensity at 3.0 Å arising from intra-molecular Cl-Cl atom pairs. When nifedipine molecules are incorporated into PVP, our models show they possess stronger NH···O bonds to the PVP polymer than felodipine molecules, which have stronger affinity for bonding to the polymer than to other felodipine molecules. Conclusions: The amorphous forms of both NIF and FEL show much weaker hydrogen bonding than found in their crystalline phases. Liquid NIF also exhibits configurations which are not observed in the crystal phases.

Amorphous Solid Dispersion↗
DOE OSTI2025

Effect of Crystallinity on Polymer Chain Compression in Polyethylene Dynamically Loaded with In-Situ X-Ray Diffraction

Abstract Polyethylene (PE) is a polymer widely used in commercial applications, where the properties of PE can be altered based on the molecular conformation. In this study, shock compression experiments with in-situ x-ray diffraction were conducted on two conformations of PE– low density polyethylene (LDPE) and ultra-high molecular weight polyethylene (UHMWPE). Encouragingly, X-ray diffraction patterns were easily observed in LDPE, which is only 25% crystalline. In agreement with a previous study on HDPE, LDPE exhibited a phase transition at 7.7 GPa. Lower pressure experiments on UHMWPE and LDPE exhibited polymer chain compression, but no occurrence of a phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2023

How Accurate Can Crystal Structure Predictions Be for High-Energy Molecular Crystals?

Molecular crystals have shallow potential energy landscapes, with multiple local minima separated by very small differences in total energy. Predicting molecular packing and molecular conformation in the crystal generally requires ab initio methods of high accuracy, especially when polymorphs are involved. We used dispersion-corrected density functional theory (DFT-D) to assess the capabilities of an evolutionary algorithm (EA) for the crystal structure prediction (CSP) of well-known but challenging high-energy molecular crystals (HMX, RDX, CL-20, and FOX-7). While providing the EA with the experimental conformation of the molecule quickly re-discovers the experimental packing, it is more realistic to start instead from a naïve, flat, or neutral initial conformation, which reflects the limited experimental knowledge we generally have in the computational design of molecular crystals. By doing so, and using fully flexible molecules in fully variable unit cells, we show that the experimental structures can be predicted in fewer than 20 generations. Nonetheless, one must be aware that some molecular crystals have naturally hindered evolutions, requiring as many attempts as there are space groups of interest to predict their structures, and some may require the accuracy of all-electron calculations to discriminate between closely ranked structures. To save resources in this computationally demanding process, we showed that a hybrid xTB/DFT-D approach could be considered in a subsequent study to push the limits of CSP beyond 200+ atoms and for cocrystals.

42 ENGINEERING↗
DOE OSTI2023

Electronic coarse-graining of long conjugated molecules: Case study of non-fullerene acceptors

By considering only one electronic state per molecule, charge transport models of molecular solids neglect intramolecular charge transfer. This approximation excludes materials with quasi-degenerate spatially separated frontier orbitals, such as non-fullerene acceptors (NFAs) and symmetric thermally activated delayed fluorescence emitters. By analyzing the electronic structure of room-temperature molecular conformers of a prototypical NFA, ITIC-4F, we conclude that the electron is localized on one of the two acceptor blocks with the mean intramolecular transfer integral of 120 meV, which is comparable with intermolecular couplings. Therefore, the minimal basis for acceptor–donor–acceptor (A–D–A) molecules consists of two molecular orbitals localized on the acceptor blocks. This basis is robust even with respect to geometry distortions in an amorphous solid, in contrast to the basis of two lowest unoccupied canonical molecular orbitals withstanding only thermal fluctuations in a crystal. The charge carrier mobility can be underestimated by a factor of two when using single site approximation for A–D–A molecules in their typical crystalline packings.

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

Ca-dimers, solvent layering, and dominant electrochemically active species in Ca(BH4)2 in THF

Abstract Divalent ions (Mg, Ca, and Zn) are being considered as competitive, safe, and earth-abundant alternatives to Li-ion electrochemistry, but present challenges for stable cycling due to undesirable interfacial phenomena. We explore the formation of electroactive species in the electrolyte Ca(BH 4 ) 2 ∣THF using molecular dynamics coupled with a continuum model of bulk and interfacial speciation. Free-energy analysis and unsupervised learning indicate a majority population of neutral Ca dimers and monomers with diverse molecular conformations and an order of magnitude lower concentration of the primary electroactive charged species – the monocation, $${\rm{CaBH}}_{4}^{+}$$ CaBH 4 + – produced via disproportionation of neutral complexes. Dense layering of THF molecules within ~1 nm of the electrode surface strongly modulates local electrolyte species populations. A dramatic increase in monocation population in this interfacial zone is induced at negative bias. We see no evidence for electrochemical activity of fully-solvated Ca 2+ . The consequences for performance are discussed in light of this molecular-scale insight.

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

Dynamic Interfacial Architectures: Cruciferin‐Stabilized Oil/Water Interfaces for Sustainable Emulsions

Stabilizing oil-water interfaces in emulsions by plant-based proteins provides sustainable and tunable ways for designing emulsions with specific properties, for food, healthcare, and pharmaceuticals. Cruciferin, a protein from rapeseed, has great potential as green emulsifier, but details about its structure and mobility at oil-water interfaces are largely unknown. Here, these properties are studied with small angle neutron and x-ray scattering, and neutron spin echo spectroscopy, analyzed by atomistic modelling of scattering curves and coarse-grained modelling, to gain insight into interface coverage, and molecular conformation and mobility at the interface. Cruciferin assumes trimeric conformations at the interface, as in solution, but with its protrusions from the central core of the subunits (“arms”) more compressed. Interfacial mobility is only marginally lower than in solution, indicating the arms still transiently extend and preserve a network, for the first time revealing the mechanism how cruciferin forms highly elastic 2d gel-like oil-water interfaces, as observed in macroscopic rheology. The high interfacial mobility may help in self-repairing non-stabilized interfacial fractions, reducing coalescence. These findings provide a deeper molecular level understanding of proteins at oil-water interfaces, which can stimulate development of new plant-based emulsion products, and contribute to the global protein transition.

cruciferin↗