Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular energetics”

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 OSTI2025

Second‐ and Third‐Order Elastic Constants of Inert and Energetic Molecular Crystals From Density Functional Theory

Complete tensors of the second- and third-order elastic constants of the organic molecular crystals acetaminophen, pentaerythritol tetranitrate (PETN), cyclotrimethylene trinitramine (RDX), cyclotetramethylene tetranitramine (HMX), 1,1-diamino-2,2-dinitroethylene (FOX-7), hexanitrohexaazaisowurtzitane (CL-20), and erythritol tetranitrate (ETN) have been calculated using dispersion-corrected density functional theory. The sets of second- and third-order elastic constants are expected to provide a more accurate and reliable description of the behavior of these materials under nonhydrostatic loads than pressure- and volume-dependent second-order elastic constants. The tensors of second-order constants have been compared with experimental data and/or other calculations when possible, and with the exception of results for CL-20 from Brillouin scattering experiments, we find good agreement. The calculated third-order elastic constants of PETN are in very good agreement with the subset of third-order constants derived from experimental wave speed measurements. The elastic anisotropies of the crystals have been estimated using the universal elastic anisotropy index, which shows that the crystals fall into three groups with low anisotropy (PETN, RDX, and CL-20), moderate anisotropy (acetaminophen, HMX, and ETN), and high elastic anisotropy (FOX-7).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantifying Energetic and Entropic Pathways in Molecular Systems

When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Lastly, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.

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

Orientation-Dependent Nanomechanical Behavior of Pentaerythritol Tetranitrate as Probed by Multiple Nanoindentation Tip Geometries

Nanoindentation can be leveraged to aid in the high fidelity modeling of dislocation mediated plasticity in pentaerythritol tetranitrate (PETN), an anisotropic energetic molecular crystal. Moreover, nanoindentation tip parameters such as tip geometry, size, and degree of acuity can be utilized to target anisotropic behavior. In this work, nanoindentation was conducted across a range of orientations on the (110) face of PETN to characterize resultant yield behavior, mechanical property measurements, and resultant slip behavior and fracture initiation. Three different indentation tips were utilized: a 3-sided pyramidal Berkovich tip, a 4-sided high aspect ratio Knoop tip, and a 90° conical tip. Ultimately, indenter tip radius was documented to impact yield behavior, whereas tip geometry affected larger scale processes such as slip, and tip acuity was the dominating factor that led to fracture. The axisymmetric conical tip, serving as a baseline, showed the least amount of variation in mechanical property measurements but also the largest distribution of maximum shear stress at which initial yielding occurred. Its high degree of acuity, however, was more prone to induce fracture at higher loads. The Knoop tip was shown to be suitable for average measurements, but also for elucidation of certain anisotropic features. A distinctly higher perceived hardness at 45° was measured with the Knoop tip, indicating less dislocation motion in that direction also observed in this work via scanning probe microscopy. Lastly, the commonly used Berkovich tip was a good compromise whereby it provided a representative volume element describing the average behavior of the material. These results can be utilized to target desired anisotropic behavior in a wider range of molecular crystals, as well as to inform theoretical considerations for dislocation mediated plasticity in PETN.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring the impact of the local environment on charge transfer states at molecular donor-acceptor heterojunctions

This program explored a broad range of charge transfer (CT) states in organic solar cells with the goal of understanding how molecular energetics, structure, and morphology influence CT state energy and dynamics, and how this in turn influences solar cell performance. We identified new types of delocalized CT states at ordered molecular heterojunctions, determined the factors that dictate CT state energetic disorder in bulk heterojunctions, determined that hot CT state dissociation is negligible for most systems, that the occupation of CT states is generally non-thermal in disordered heterojunctions, and that quasi-equilibrium generally does not hold in organic solar cells. These findings allowed us to demonstrate the lowest loss in potential of any organic solar cell to date and provide valuable guidance for designing higher efficiency organic solar cells in the future.

14 SOLAR ENERGY↗
DOE OSTI2024

Octupole moment driven free charge generation in partially chlorinated subphthalocyanine for planar heterojunction organic photodetectors

In this study, high-performance organic photodetectors are presented which utilize a pristine chlorinated subphthalocyanine photoactive layer. Optical and optoelectronic analyses indicate that the device photocurrent is primarily generated through direct charge generation within the chlorinated subphthalocyanine layer, rather than exciton separation at layer interfaces. Molecular modelling suggests that this direct charge generation is facilitated by chlorinated subphthalocyanine high octupole moment (-80 DÅ 2 ), which generates a 200 meV shift in molecular energetics. Increasing the thickness of chlorinated subphthalocyanine leads to faster response time, correlated with a decrease in trap density. Notably, photodetectors with a 50 nm thick chlorinated subphthalocyanine photoactive layer exhibit detectivities approaching 10 13 Jones, with a dark current below 10 -7 A cm -2 up to -5 V. Based on these findings, we conclude that high octupole moment molecular semiconductors are promising materials for high-performance organic photodetectors employing single-component photoactive layer.

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

Thermal conductivity tensor of γ and ε -hexanitrohexaazaisowurtzitane as a function of pressure and temperature

Using reverse non-equilibrium molecular dynamics simulations, we have determined the dependences on temperature and pressure of the thermal conductivity tensors for the monoclinic γ and ε polymorphs of hexanitrohexaazaisowurtzitane (HNIW or CL20). In this work, a recently developed non-reactive force field [X. Bidault and S. Chaudhuri, RSC Adv. 9, 39649–39661 (2019)], designed to study polymorphism and phase transitions in CL20, is employed. The effects of temperature and pressure are investigated between 200 and 500 K and up to 0.5 GPa for γ-CL20 and 2 GPa for ε-CL20. In order to obtain the full thermal conductivity tensor, $κ_{ij}$, for the monoclinic crystals, four distinct heat propagation directions are used. We find that $κ_{ij}$ for both polymorphs is more isotropic than for other energetic molecular crystals, including α- and γ-RDX, β-HMX, and PETN, with a maximum difference of 9.8% between orientations observed at 300 K and 0 GPa for γ-CL20 and a maximum difference of 4.8% for ε-CL20. The average thermal conductivity, $\bar{κ}$, of ε-CL20 is 6.4% larger than that of γ-CL20 at 300 K and 0 GPa. Analytic linear functions of the inverse temperature and the pressure are provided, which fit the data well and can be used to predict the thermal conductivity of both polymorphs for any orientation, pressure, and temperature in and around the fitting range. Our predictions agree reasonably well with the limited available experimental data, for which the polymorph type is unknown.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulation and Emulation of X-Ray Diffraction from Dynamic Compression Experiments

Many important aspects of the dynamic thermo-mechanical response of materials occur at the mesoscale, i.e. a physical scale of interactions smaller than what can be adequately described by homogenous behaviors, yet larger than the scale of the atomic lattice. Concurrent advancements in computational power, continuum theory, and experimental diagnostics are enabling unprecedented understanding of such interactions. However, we cannot develop a sufficient level of confidence in such mesoscale capability until the constitutive description of the underlying constituents is reliably representative of their actual physical behavior. Therefore, there is a strong need to combine experimental, modeling, and data-science techniques to validate models of the thermomechanical response of individual single crystals. One experimental diagnostic with high potential impact to shock physics and materials science is in-situ x-ray diffraction. This paper is primarily focused on simulation of x-ray diffraction in shock physics, but with an aim toward quantifying parametric uncertainty of simulation models. Here, we develop and demonstrate a data-science and model-driven approach to constrain the parameterization of continuum models of crystal lattice deformation associated with the shock response of crystalline materials. The framework is built around the connection between continuum hydrodynamic simulations of lattice deformation and a new Bragg diffraction simulation code, BarberShop. The dynamic deformation of a crystal lattice is modeled using the DiscoFlux model within an arbitrary Lagrangian-Eulerian hydrodynamic code, FLAG. These detailed continuum simulations of lattice deformation can be computationally slow, thus a statistical model is used to emulate the evolution of lattice deformation fields in time and across the considered model parameter space. Emulated lattice deformation fields can then be generated rapidly for any combination of physics model parameters. In turn, these fields can be fed into BarberShop to realize a rapid prediction of Bragg diffraction patterns associated with particular values of physics model parameters. The framework enables parameterization of the single crystal model to obtain Bragg diffraction patterns that most closely resemble a corresponding measurement. Furthermore, the framework naturally provides sensitivities of the lattice deformation to the physics parameters. We highlight the utility of this framework through the application to a synthetic closed-loop inverse problem leading to the parameterization of a single crystal material model. As a model problem, we consider the dynamic response of the energetic molecular crystal, cyclotrimethylenetrinitramine (or RDX), under dynamic compression induced by simulated flyer plate impact experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Maximum Entropy Theory of Multiscale Coarse-Graining via Matching Thermodynamic Forces: Application to a Molecular Crystal (TATB)

The MSCG/FM (multiscale coarse-graining via force-matching) approach is an efficient supervised machine learning method to develop microscopically informed coarse-grained (CG) models. Here we present a theory based on the principle of maximum entropy (PME) enveloping the existing MSCG/FM approaches. This theory views the MSCG/FM method as a special case of matching the thermodynamic forces from the extended ensemble described by the set of thermodynamic (relevant) system coordinates. This set may include CG coordinates, the stress tensor, applied external fields, and so forth, and may be characterized by nonequilibrium conditions. Following the presentation of the theory, we discuss the consistent matching of both bonded and nonbonded interactions. The proposed PME formulation is used as a starting point to extend the MSCG/FM method to the constant strain ensemble, which together with the explicit matching of the bonded forces is better suited for coarse-graining anisotropic media at a submolecular resolution. The theory is demonstrated by performing the fine coarse-graining of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a well-known insensitive molecular energetic material, which exhibits highly anisotropic mechanical properties.

1,3,5-triamino-2,4,6-trinitrobenzene↗
DOE OSTI2023

Structural elucidation of polydopamine facilitated by ionic liquid solvation

Minimal understanding of the formation mechanism and structure of polydopamine (pDA) and its natural analogue, eumelanin, impedes the practical application of these versatile polymers and limits our knowledge of the origin of melanoma. The lack of conclusive structural evidence stems from the insolubility of these materials, which has spawned significantly diverse suggestions of pDA's structure in the literature. We discovered that pDA is soluble in certain ionic liquids. Using these ionic liquids (ILs) as solvents, we present an experimental methodology to solvate pDA, enabling us to identify pDA's chemical structure. The resolved pDA structure consists of self-assembled supramolecular aggregates that contribute to the increasing complexity of the polymer. The underlying molecular energetics of pDA solvation and a macroscopic picture of the disruption of the aggregates using IL solvents have been investigated, along with studies of the aggregation mechanism in water.

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

Subsurface spectroscopy of heterogeneous materials using optical wavefront shaping

Plastic-bonded explosives and propellants consist of energetic molecular crystals embedded in a polymer matrix that also includes additives, such as taggants, plasticizers, grit, antioxidants, etc. Furthermore, this heterogeneous composition renders these materials optically opaque, limiting optical characterization techniques (e.g., Raman spectroscopy) to probing chemical reactions at the sur-face. However, many reactions of interest are believed to occur inside the material where current optical techniques can not probe. To address this challenge, we have developed a new optical technique that utilizes feedback assisted wavefront shaping (using spatial light modulators) to focus probe light inside a heterogeneous material, such that we can perform subsurface spectroscopy. Recently, we demonstrated this technique by monitoring subsurface photodegradation and thermal degradation of Eu-doped molecular crystals using both fluorescence and Raman spectroscopy. Based off these successes we are now looking into improvements to this technique to allow for greater time resolution, including replacing our spatial light modulators with fast digital mirror devices, using ns optical phase conjugation, and self-healing optical beams.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Predicting Melt Curves of Energetic Materials Using Molecular Models

Here, in this work, the solid–liquid coexistence curves of classical fully flexible atomistic models of α-RDX and β-HMX were calculated using thermodynamically rigorous methodologies that identify where the free energy difference between the phases is zero. The free energy difference between each phase at a given state point was computed using the pseudosupercritical path (PSCP) method, and Gibbs–Helmholtz integration was used to evaluate the solid–liquid free energy difference as a function of temperature. This procedure was repeated for several pressures to determine points along the coexistence curve, which were then fit to the Simon–Glatzel functional form. While effective, this method is computationally expensive. An alternative approach is to compute the melting point at a single pressure via the PSCP method, and then use the Gibbs–Duhem integration technique to trace out the coexistence curve in a more computationally economical manner. Both approaches were used to determine the coexistence curve of α-RDX. The Gibbs–Duhem integration method was shown to generate a melt curve that is in good agreement with the PSCP-derived melt curve, while only costing ~10 % of the computational resources used for the PSCP method. For α-RDX, the predicted melting temperature increases significantly more for a given increase in pressure when compared to available experimental data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Biased Spin-State Energetics of Fe(II) Molecular Complexes within Density-Functional Theory and the Linear-Response Hubbard U Correction

The spin-state energetics of six Fe(II) molecular complexes are computed using the linear-response Hubbard U approach within DFT. The adiabatic energy differences, ΔE H-L , between the high-spin (S = 2) and the low-spin (S = 0) states are computed and compared with accurate-coupled cluster-corrected CASPT2 results. We show that DFT+U fails in correctly capturing the ground state for strong-field ligands yielding ΔE H-L that are almost constant throughout the molecular series. This bias toward high spin together with the metal/ligand charge transfer upon U correction are here quantified and explained using molecular orbital diagrams involving both σ- and π-bonding interactions. Furthermore, with increasing ligand-field strengths this bias also increases owing to the stronger molecular character of the metal/ligand Kohn–Sham orbitals thus resulting in large deviations from the reference larger than 4 eV. Smaller values of U can be employed to mitigate this effect and recover the right energetics.

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

Dynamical electron correlation and the chemical bond. I. Covalent bonds in AH and AF (A = B–F)

Dynamical electron correlation has a major impact on the computed values of molecular properties and the energetics of molecular processes. This study focused on the effect of dynamical electron correlation on the spectroscopic constants ($R$ e , $ω$ e , $D$ e ), and potential energy curves, Δ$E$($R$), of the covalently bound AH and AF molecules, A = B–F. The changes in the spectroscopic constants (Δ$R$ e , Δ$ω$ e , Δ$D$ e ) caused by dynamical correlation are erratic and, at times, even surprising. These changes can be understood based on the dependence of the dynamical electron correlation energies of the AH and AF molecules as a function of the bond distance, i.e., Δ$E$ DEC ($R$). At large $R$, the magnitude of Δ$E$ DEC ($R$) increases nearly exponentially with decreasing $R$, but this increase slows as $R$ continues to decrease and, in many cases, even reverses at very short $R$. The changes in Δ$E$ DEC ($R$) in the region around $R$ e were as unexpected as they were surprising, e.g., distinct minima and maxima were found in the curves of Δ$E$ DEC ($R$) for the most polar molecules. The variations in Δ$E$ DEC ($R$) for $R$ ≲ $R$ e are directly correlated with major changes in the electronic structure of the molecules as revealed by a detailed analysis of the spin-coupled generalized valence bond wave function. In conclusion, the results reported here indicate that we have much to learn about the nature of dynamical electron correlation and its effect on chemical bonds and molecular properties and processes.

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

Supercritical fluids behave as complex networks

Abstract Supercritical fluids play a key role in environmental, geological, and celestial processes, and are of great importance to many scientific and engineering applications. They exhibit strong variations in thermodynamic response functions, which has been hypothesized to stem from the microstructural behavior. However, a direct connection between thermodynamic conditions and the microstructural behavior, as described by molecular clusters, remains an outstanding issue. By utilizing a first-principles-based criterion and self-similarity analysis, we identify energetically localized molecular clusters whose size distribution and connectivity exhibit self-similarity in the extended supercritical phase space. We show that the structural response of these clusters follows a complex network behavior whose dynamics arises from the energetics of isotropic molecular interactions. Furthermore, we demonstrate that a hidden variable network model can accurately describe the structural and dynamical response of supercritical fluids. These results highlight the need for constitutive models and provide a basis to relate the fluid microstructure to thermodynamic response functions.

42 ENGINEERING↗
DOE OSTI2025

Quantum Monte Carlo Benchmarking of Molecular Adsorption on Graphene-Supported Single Pt Atom

The precise understanding of adsorption energetics and molecular geometry at catalytic sites is fundamental for advancing catalysis, particularly under the constraints of resource efficiency and environmental sustainability. Here, this study benchmarks the performance of density functional theory (DFT) calculations against diffusion Monte Carlo (DMC) calculations for adsorption properties of small gas molecules relevant to CO oxidation—namely O 2 , CO, CO 2 , and atomic oxygen—on a single Pt atom supported by pristine graphene. Our findings reveal that DMC calculations provide a significantly different landscape of adsorption energetics compared to DFT results. Notably, DFT predicts different lowest-energy configurations and spin states, particularly for O 2 , which suggests potential discrepancies in predicting the catalytic behavior. Furthermore, this study identifies the critical issue of CO poisoning, highlighted by the large disparity between the DMC adsorption energies of O 2 (−1.23(2) eV) and CO (−3.37(1) eV), which can inhibit the catalytic process. These results emphasize the necessity for more sophisticated computational approaches in catalysis research, aiming to refine the prediction accuracy of reaction mechanisms and to enhance the design of more effective catalysts.

Ahn, Jeonghwan [University of Illinois at Urbana-C↗
DOE OSTI2023

Progressive evolutions of non-fused ring electron acceptors toward efficient organic solar cells with improved photocurrent and reduced energy loss

Resolving competitive counterbalance between high external quantum efficiency (EQE) response and low energy loss is critical for achieving efficient organic solar cells (OSCs), including those based on low-cost non-fused-ring electron acceptor (NFREA). In this work, we performed a comparable study by utilizing four NFREAs with progressive evolutions in molecular structures, energetic levels, spectral coverages, crystallinity, etc, to pair with polymer donor J52. By regulating the terminal and side chains of NFREAs, the enhanced intramolecular charge transfer (ICT) effect allows O-PC-EH possessing the broader spectral coverage for higher photocurrent, the reduced energetic offset enables O-PC-EH-based device owning higher luminescence property for lower energy loss, the strengthened crystallinity endows O-PC-EH-based blend better charge transport property for less charge recombination. With those comprehensive managements from molecule to blend properties, the J52:O-PC-EH-based OSCs achieved the best efficiency of 14.4% (certified: 14.1%) with benefits in all the three device parameters. Here, this work reveals the importance of fine-tuning molecular structures of NFREAs according to paired polymer donor for maximizing the device performances.

14 SOLAR ENERGY↗
DOE OSTI2021

Nonintuitive Surface Self-Assembly of Functionalized Molecules on Ag(111)

The fabrication of nanomaterials involves selfordering processes of functional molecules on inorganic surfaces. To obtain specific molecular arrangements, a common strategy is to equip molecules with functional groups. However, focusing on the functional groups alone does not provide a comprehensive picture. Especially at interfaces, processes that govern self-ordering are complex and involve various physical and chemical effects, often leading to unexpected structures, as we showcase here on the example of a homologous series of quinones on Ag(111). Naively, one could expect that such quinones, which all bear the same functionalization, form similar motifs. In salient contrast, our joint theoretical and experimental study shows that profoundly different structures are formed. Using a machine-learning-based structure search algorithm, we find that this is due to a shift of the balance of three antagonizing driving forces: adsorbate-substrate interactions governing adsorption sites, adsorbate-adsorbate interactions favoring close packing, and steric hindrance inhibiting certain otherwise energetically beneficial molecular arrangements. The theoretical structures show excellent agreement with our experimental characterizations of the organic/ inorganic interfaces, both for the unit cell sizes and the orientations of the molecules within. The nonintuitive interplay of similarly important interaction mechanisms will continue to be a challenging aspect for the design of functional interfaces. With a detailed examination of all driving forces, we are, however, still able to devise a design principle for self-assembly of functionalized molecules.

structure prediction↗