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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Deep learning of dynamically responsive chemical Hamiltonians with semiempirical quantum mechanics

Conventional machine-learning (ML) models in computational chemistry learn to directly predict molecular properties using quantum chemistry only for reference data. While these heuristic ML methods show quantum-level accuracy with speeds several orders of magnitude faster than traditional quantum chemistry methods, they suffer from poor extensibility and transferability; i.e., their accuracy degrades on large or new chemical systems. Incorporating quantum chemistry frameworks into the ML models directly solves this problem. Here we take the structure of semiempirical quantum mechanics (SEQM) methods to construct dynamically responsive Hamiltonians. SEQM methods use empirical parameters fitted to experimental properties to construct reduced-order Hamiltonians, facilitating much faster calculations than ab initio methods but with compromised accuracy. By replacing these static parameters with machine-learned dynamic values inferred from the local environment, we greatly improve the accuracy of the SEQM methods. Trained on molecular energies and atomic forces, these dynamically generated Hamiltonian parameters show a strong correlation with atomic hybridization and bonding. Trained with only about 60,000 small organic molecular conformers, the resulting model retains interpretability, extensibility, and transferability when testing on much larger chemical systems and predicting various molecular properties. Overall, this work demonstrates the virtues of incorporating physics-based descriptions with ML to develop models that are simultaneously accurate, transferable, and interpretable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2021

A High Energy X-ray Diffraction Study of Amorphous Indomethacin

Amorphous pharmaceuticals often possess a wide range of molecular conformations and bonding arrangements. The x-ray pair distribution function (PDF) method is a powerful technique for the characterization of variations in both intra-molecular and inter-molecular packing arrangements. Here, the x-ray PDF of amorphous Indomethacin is shown to be particularly sensitive to the preferred orientations of the chlorobenzyl ring found in isomers in the crystalline state. In some cases, the chlorobenzyl ring has no preferred torsional angle in the amorphous form, while in others evidence of distinct isomer orientations are observed. Amorphous samples with no preferred torsion angles of the chlorobenzyl ring are found to favor enhanced inter-molecular hydrogen bonding, and this is reflected in the intensity of the first sharp diffraction peak. These significant variations in structure rule out amorphous Indomethacin as a possible standard for x-ray PDF measurements. At high humidity, time resolved PDF's for >40 h reveal water molecules forming hydrogen bonds with Indomethacin molecules. Finally, a simple linear hydrogen bond model indicates that water molecules in the wet amorphous form have similar hydrogen bond strengths to those found between Indomethacin dimers or chains in the dry amorphous form.

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

Vascular K ATP channel structural dynamics reveal regulatory mechanism by Mg-nucleotides

Vascular tone is dependent on smooth muscle K ATP channels comprising pore-forming Kir6.1 and regulatory SUR2B subunits, in which mutations cause Cantú syndrome. Unique among K ATP isoforms, they lack spontaneous activity and require Mg-nucleotides for activation. Structural mechanisms underlying these properties are unknown. Here, we determined cryogenic electron microscopy structures of vascular K ATP channels bound to inhibitory ATP and glibenclamide, which differ informatively from similarly determined pancreatic K ATP channel isoform (Kir6.2/SUR1). Unlike SUR1, SUR2B subunits adopt distinct rotational “propeller” and “quatrefoil” geometries surrounding their Kir6.1 core. The glutamate/aspartate-rich linker connecting the two halves of the SUR-ABC core is observed in a quatrefoil-like conformation. Molecular dynamics simulations reveal MgADP-dependent dynamic tripartite interactions between this linker, SUR2B, and Kir6.1. The structures captured implicate a progression of intermediate states between MgADP-free inactivated, and MgADP-bound activated conformations wherein the glutamate/aspartate-rich linker participates as mobile autoinhibitory domain, suggesting a conformational pathway toward K ATP channel activation.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2021

Conformer-specific photochemistry imaged in real space and time

Conformation-dependent dynamics play an important role in organic chemistry syntheses such as electrocyclic reactions, as well as in biological processes such as protein folding. However, current time-resolved experimental methods struggle to distinguish conformers from each other, and conformational isomerism is usually analyzed through reactant and product distributions. Using a combination of mega–electron volt ultrafast electron diffraction and quantum wave packet simulations, Champenois et al. directly followed the photochemical electrocyclic ring opening of the molecule α-phellandrene with femtosecond time resolution and confirmed that the transformation of a specific molecular conformer follows the famous Woodward-Hoffmann rules. The proposed method is potentially a powerful tool to follow conformer specificity in various organic and biological systems in real time.

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

Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene

In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S 2 /S 3 (D 2d ) → S 1 (D 8h ) → S o (D 4h ) → S o (D 2d ) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. Here, the obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.

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↗
DOE OSTI2025

Fluorinated polythiophenes with ester side chains for boosting the V OC and efficiency in non-fullerene polymer solar cells

As the power conversion efficiency (PCE) of organic photovoltaics (OPVs) approximates the 19 % threshold, wide bandgap (WBG) polythiophenes (PTs) have gained increasing attention due to their simple chemical structures and ease of synthesis, making them promising candidates for large-scale production. However, the benchmark polymer poly(3-hexylthiophene) (P3HT) is limited by its high-lying HOMO energy level, which restricts the open-circuit voltage (V OC ) in solar cells. Here, in this study, we introduce a novel series of PT derivatives (PDC8-T, PDC8-3T, PDC16-3T, and PDC16-3T-2F) featuring ester side chains designed to fine-tune electronic properties through a streamlined three-step synthesis. Additionally, we incorporated a π-spacer to reduce steric hindrance and elongated alkyl side chains to improve solubility and processability. Compared to P3HT, these PT derivatives demonstrate a significant reduction in HOMO energy levels, lowering by approximately 0.3–0.4 eV. Among them, PDC16-3T-2F—with fluorine atom substitution—achieves the lowest HOMO energy level, induces a coplanar molecular conformation, and enhances polymer aggregation behavior. We evaluated these PT derivatives in inverted non-fullerene bulk-heterojunction (NFA BHJ) OPVs. The PDC8-3T device showed a relatively low PCE of 0.69 %, with a V OC of 0.76 V, a short-circuit current density (J SC ) of 3.32 mA/cm 2 , and a fill factor (FF) of 27.3 %. In contrast, the PDC16-3T-2F device achieved an impressive PCE of 7.21 %, with a V OC of 0.85 V, a J SC of 14.60 mA/cm 2 , and an FF of 58.4 %. This remarkable improvement is attributed to the fluorine substitution, which not only enhances molecular orientation but also downshifts the HOMO energy level and further boosts the V OC . Hence, these molecular design strategies have led to a fibrillar bicontinuous interpenetrating network with optimal nanoscale phase separation within the active layer.

14 SOLAR ENERGY↗
DOE OSTI2023

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

A Three-in-One Hybrid Strategy for High-Performance Semiconducting Polymers Processed from Anisole

The development of semiconducting polymers with good processability in green solvents and competitive electrical performance is essential for realizing sustainable large-scale manufacturing and commercialization of organic electronics. A major obstacle is the processability-performance dichotomy that is dictated by the lack of ideal building blocks with balanced polarity, solubility, electronic structures, and molecular conformation. Herein, through the integration of donor, quinoid and acceptor units, an unprecedented building block, namely TQBT, is introduced for constructing a serial of conjugated polymers. The TQBT, distinct in non-symmetric structure and high dipole moment, imparts enhanced solubility in anisole—a green solvent—to the polymer TQBT-T. Furthermore, PTQBT-T possess a highly rigid and planar backbone owing to the nearly coplanar geometry and quinoidal nature of TQBT, resulting in strong aggregation in solution and localized aggregates in film. Remarkably, PTQBT-T films spuncast from anisole exhibit a hole mobility of 2.30 cm 2 V -1 s -1 , which is record high for green solvent-processable semiconducting polymers via spin-coating, together with commendable operational and storage stability. The hybrid building block emerges as a pioneering electroactive unit, shedding light on future design strategies in high-performance semiconducting polymers compatible with green processing and marking a significant stride towards ecofriendly organic electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Regiospecific Incorporation of Acetylene Linker in High‐Electron Mobility Dimerized Acceptors for Organic Solar Cells with High Efficiency (18.8%) and Long 1‐Sun Lifetime (> 5000 h)

Abstract The commercialization of organic solar cells (OSCs) requires both high power conversion efficiency (PCE) and long‐term stability. However, the lifetime of the OSCs containing small‐molecule acceptors (SMA) should be significantly enhanced. In this study, a series of planarity‐controlled is developed, high electron mobility dimerized SMAs (DSMAs) and realize OSCs with high‐performance (PCE = 18.8%) and high‐stability ( t 80% lifetime = 5380 h under 1‐Sun illumination). An acetylene linker with a planar triple bond is designed for dimerization of SMA units to afford the high backbone planarity necessary to achieve high crystallinity and electron mobility. To further engineer the molecular conformation and electron mobility of the DSMAs, different regioisomers of a Y‐based SMA are dimerized to yield three regioisomerically distinct DSMAs, denoted as DYA‐I, DYA‐IO, and DYA‐O, respectively. It is found that the crystallinity, electron mobility, and glass transition temperature of the DSMAs gradually increase in the order of DYA‐O, DYA‐IO, and DYA‐I, which, in turn, enhance the PCE and device stability of the resulting OSCs; DYA‐O (PCE = 16.45% and t 80% lifetime = 3337 h) < DYA‐IO (PCE = 17.54% and t 80% lifetime = 4255 h) < DYA‐I (PCE = 18.83% and t 80% lifetime = 5380 h).

Chemistry↗
DOE OSTI2025

Pressure–Modulated Luminescence Enhancement and Quenching in a Hydrogen–Bonded Organic Framework

Light emission in the solid state is central for illumination, sensing, and imaging applications. Unlike luminescence in dilute solutions, where the excited states are unimolecular in nature, intermolecular interaction plays a significant role in the quantum yield of solid-state luminophores, manifested as competing aggregation-caused quenching (ACQ) and aggregation-induced enhancement (AIE). Both effects are extensively studied in various systems; however, it remains unclear how their competition depends on molecular conformation and intermolecular stacking. Here the direct observation of pressure-modulated AIE-ACQ competition in a crystalline hydrogen-bonded organic framework (HOF) is reported. Using in situ spectroscopies and computational modeling, the intramolecular vibration and intermolecular π–π stacking directly responsible for the non-radiative decay of the excited state are identified. The extent of these two contributions is modulated by hydrostatic pressure and guest molecules in the HOF pores. Furthermore, this work demonstrates a physically neat model system to understand and control solid-state luminescence, and a potential material platform for piezoluminescent sensing.

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

Determination of the chemical shift tensor anisotropy and asymmetry of strongly dipolar coupled protons under fast MAS

Orientationally-dependent interactions such as dipolar coupling, quadrupolar coupling, and chemical shift anisotropy (CSA) contain a wealth of spatial information that can be used to elucidate molecular conformations and dynamics. To determine the sign of the chemical shift tensor anisotropy parameter (δ aniso ), both the |m| = 1 and |m| = 2 components of the CSA need to be symmetry allowed, while the recoupling of the |m| = 1 term is accompanied with the reintroduction of homonuclear dipolar coupling components. Therefore, previously suggested sequences which solely recouple the |m| = 2 term cannot determine the sign a 1 H's δ aniso in a densely-coupled network. Here, we demonstrate the CSA recoupling of strongly dipolar coupled 1 H spins using the $Cn^1_n (90_{0}360_{180}540_{0}360_{180}90_{0})$ sequence. This pulse scheme recouples both the |m| = 1 and |m| = 2 CSA terms but the scaling factors for the homonuclear dipolar coupling terms are zeroed. Consequently, the sequence is sensitive to the sign of δ aniso but is not influenced by homonuclear dipolar interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗