SEARCH · Engineering Papers
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.
The vibrational and electronic properties of supported polyoxometalates determine their suitability for potential applications as molecular qubits.
Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.
Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe–Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.
The rise in energy demands and the deleterious environmental issues related to fossil fuels has led to a surge of interest in hydrogen as a "green" alternative. Hydrogen's extraordinary energy density makes it a potential energy and economic "power"-house. Significant research has been dedicated to materials-based hydrogen storage. One area, liquid organic hydrogen carriers (LOHC) is of substantial interest for the reversible transportation of hydrogen from production to end-use facilities. There are challenges associated with this technology including the dependency on precious metal-based catalysts. Recent work in frustrated Lewis pair (FLP) catalysis demonstrates promise for addressing these challenges. Here this review is focused on assessing recent literature on the utilization of intermolecular FLP main group catalysts for improved hydrogenation/dehydrogenation of various substrates including potential LOHC complexes. This review will present an overview of FLPs, highlight potential hydrogen storage applications, and propose areas where knowledge gaps exist that require further investigations.
Biomass-derived bio-oils are emerging as sustainable, low-carbon alternatives for construction materials, particularly in innovative pavement applications, 'BioPave'. Here, this study evaluates bio-oils from various sources, including algae and wastewater biomass, revealing variations in carbon, nitrogen, and sulfur contents, along with viscosity. We focused on their application in outdoor construction, assessing their resistance to thermal and solar radiation, and moisture. The study also examines their interaction with siliceous surfaces and their influence on asphalt adhesion in BioPave systems under moisture conditions. Laboratory experiments and computational modeling demonstrated molecular composition significantly influenced bio-oils' responses to thermal and UV exposure. Bio-oils rich in polar groups showed hardening upon thermal treatment, while those with higher concentrations of saturated aliphatics remained more stable. Furthermore, our findings highlight the significance of dosage control in maximizing beneficial effects of bio-oils in intermolecular interactions at bitumen-aggregate interface in BioPave applications. These findings offer insights into the potential of bio-oils in sustainable construction and the need for dosage control to optimize pavement performance.
Conductive polymers are promising materials as active elements for energy storage and conversion devices due to mixed ion–electron conduction. The ion diffusion coefficient is a relative measure of the efficacy of ion transport, allowing for comparison between materials and electrochemical conditions. In this work, diffusion coefficients of hexafluorophosphate (PF6–) counterions in poly(3-alkylthiophene) (P3AT) materials are measured as a function of both side-chain length and microstructure using electrochemical impedance spectroscopy (EIS). For semi-crystalline films, the diffusion coefficient is found to be anomalous and nearly independent of applied electrochemical potential. Here, the anomalous behavior of diffusion indicates that spin casting yields compact films with an enthalpic barrier to ion transport, attributed to ionic trapping. Diffusion coefficient values ~10 –11 cm 2 s –1 were measured for all films, indicating interchain spacing, in the absence of strong intermolecular interactions with the electrolyte, is not a viable design strategy to control ion transport. For the prototypical system of poly(3-hexylthiophene), we observe almost no potential dependence in ion transport for regioregular and regiorandom films of comparable molecular weight, with both exhibiting anomalous diffusion. Alternatively, changing the microstructure of poly(3-hexylthiophene) to a mostly amorphous, ion-imprinted structure yields ~500× increase in the diffusion coefficient to ~2 × 10 –8 cm 2 s –1 at 0.8 V vs. Ag/Ag + with behavior closer to ordinary diffusion. Collectively, these results indicate new insight into ion transport in conductive polymers, where ionic trapping effects can be mitigated through electrodeposition protocols over post-synthesis processing (i.e. spin coating).
The structural investigation of amorphous pharmaceuticals is of paramount importance in comprehending their physicochemical stability. However, it has remained a relatively underexplored realm primarily due to the limited availability of high-resolution analytical tools. Here in this study, we utilized the combined power of X-ray pair distribution functions (PDFs) and solid-state nuclear magnetic resonance (ssNMR) techniques to probe the molecular packing of amorphous posaconazole and its amorphous solid dispersion at the molecular level. Leveraging synchrotron X-ray PDF data and employing the empirical potential structure refinement (EPSR) methodology, we unraveled the existence of a rigid conformation and discerned short-range intermolecular C–F contacts within amorphous posaconazole. Encouragingly, our ssNMR 19 F– 13 C distance measurements offered corroborative evidence supporting these findings. Furthermore, employing principal component analysis on the X-ray PDF and ssNMR data sets enabled us to gain invaluable insights into the chemical nature of the intermolecular interactions governing the drug–polymer interplay. These outcomes not only furnish crucial structural insights facilitating the comprehension of the underlying mechanisms governing the physicochemical stability but also underscore the efficacy of synergistically harnessing X-ray PDF and ssNMR techniques, complemented by robust modeling strategies, to achieve a high-resolution exploration of amorphous structures.
By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. Here, we find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.
Singlet fission (SF) is the process of converting an excited singlet to a pair of excited triplets. Harvesting two charges from a single photon has the potential to increase photovoltaic device efficiencies. Acenes, such as tetracene and pentacene, are model molecules for studying SF. Despite SF being an endoergic process for tetracene and exoergic for pentacene, both acenes exhibit near unity SF quantum efficiencies, raising questions about how tetracene can overcome the energy barrier. Here, we use recently developed instrumentation to measure inelastic neutron scattering (INS) while optically exciting the model molecules using two different excitation energies. The spectroscopic results reveal intermolecular structural relaxation due to the presence of a triplet excited state. Here, the structural dynamics of the combined excited state molecule and surrounding tetracene molecules are further studied using time-dependent density functional theory (TD-DFT), which shows that the singlet and triplet levels shift due to the excited state geometry, reducing the uphill energy barrier for SF to within kT.
Octopole moments for methane and carbon tetrafluoride derived from second virial coefficients temperature dependence, using spherical shell potential model
Long-range intermolecular interaction plays an essential role in (ultra)cold collisions. For monomers with degenerate electronic states, an accurate description of the long-range interaction can be challenging due to the potential breakdown of the Born–Oppenheimer approximation. In this work, we developed a general method to construct long-range interaction potential energy surfaces (IPESs) between arbitrary molecules using the degenerate perturbation theory of Schrödinger. To further simplify the application of this method, a group-theoretical method is introduced to identify nonzero components of electrostatic properties of monomers, especially for the nonlinear molecules. To demonstrate this approach, we report the long-range IPESs of the O( 3 P)–OH( 2 Π) and OH( 2 Π)–OH( 2 Π) systems. The excellent agreement with ab initio calculations in large separations, even near electronic degeneracies, confirms the accuracy of our approach. The method is expected to lay the groundwork for the accurate investigation of (ultra)cold collisional dynamics involving open-shell species.
Type 3 Porous Liquids (PLs) are a new class of materials that offer transformative potential for gas capture and utilization. These PLs are comprised of a bulky solvent with a suspended empty nanoporous material that enables selective and high-capacity capture from dilute and complex gas steams. The relationship between the nanoporous material, the structure of the solvent molecules, and the time dependence of solvent adsorption into the nanoporous material governs long term adsorption properties. Herein, molecular dynamics (MD) calculations evaluated the solvent dynamics of nine neat solvents with a broad range of chemical identities and experimental density measurements probed the time-dependent adsorption of these solvents into the ZIF-8 metal-organic frameworks (MOFs) to form PL dispersions over three weeks. Two of the nine dispersion, using glyceryl triacetate and 2’-hydroxyacetopheneone as solvents, presented sufficiently low solvent infiltration—characterized by solvent adsorption less than 40% of the ZIF-8 pore volume—to be viable as PLs. The experimental solvent adsorption data for ZIF-8 dispersed in the four aromatic solvents (acetophenone, methylbenzoate, 2-isopropylphenol, and 2’-hydroxyacetophenone) was fit to a kinetic model to quantify rate of solvent adsorption into ZIF-8. The rates of adsorption of these four similar-sized solvents demonstrated that solvent adsorption is slowed and limited by solvent clustering dynamics that are, in turn, driven by intermolecular hydrogen bonding. Ultimately, experimental solvent sorption measurements indicated that glyceryl triacetate and 2’-hydroxyacetophenone with ZIF-8 formed PLs with the most stable porosity and the greatest potential for gas capture. In conclusion, the combined experimental-computational materials exploration framework used here reveals a novel relationship between molecular scale solvent dynamics and macroscale solvent adsorption kinetics critical for the discovery and rapid evaluation of Type 3 PL compositions.
Direct air capture (DAC) materials enable the removal of CO 2 from the atmosphere, but improving their efficiency requires a detailed understanding of the intermolecular interactions that govern CO 2 sorption and release. Here, we present an experimental electron density study of methylglyoxal-bis(iminoguanidine) (MGBIG), a promising DAC material, using high-resolution X-ray and neutron diffraction data combined with quantum crystallographic analysis. This approach bridges theoretical and experimental data by quantifying electron density distributions and revealing how hydrogen bonds stabilize CO 2 -derived carbonate phases and may influence the desorption behavior. We identify distinct hydrogen-bonding environments in two crystalline carbonate phases: P1, a transient kinetic product, and P3, a thermodynamically stable phase. Multipolar refinement and electrostatic potential and multipole moment calculations precisely map electron density distributions, revealing key hydrogen bonds involved in CO 2 capture. Topological analysis of electron density highlights a cooperative hydrogen-bonding network in the thermodynamically favored P3 phase, where enhanced electron density delocalization and water-mediated interactions contribute to a more stable lattice. Energetic analyses confirm that stronger hydrogen bonding networks enhance the stability of P3 with a binding energy of −607.0 kJ/mol and greater lattice stability (−847.3 kJ/mol) compared to P1 (−302.5 and −571.0 kJ/mol, respectively). Electrostatic potential maps further illustrate polarization patterns that may influence the stability of the binding of CO 2 and release conditions. These findings establish a direct experimental framework for linking electron density distributions to intermolecular interactions in DAC materials, providing a rational design strategy for optimizing sorbents with improved CO 2 capture efficiency and reduced energy demands.
Multipole moments such as charge, dipole, and quadrupole are often invoked to rationalize intermolecular phenomena, but a low-order multipole expansion is rarely a valid description of electrostatics at the length scales that characterize nonbonded interactions. This is illustrated by examining several common misunderstandings rooted in erroneous electrostatic arguments. First, the notion that steric repulsion originates in Coulomb interactions is easily disproved by dissecting the interaction potential for Ar2. Second, the Hunter-Sanders model of π–π interactions, which is based on quadrupolar electrostatics, is shown to have no basis in accurate calculations. Third, curved “buckybowls” exhibit unusually large dipole moments, but these are ancillary to the forces that control their intermolecular interactions, as illustrated by two examples involving corannulene. Finally, the assumption that interactions between water and small anions are dictated by the dipole moment of H 2 O is shown to be false in the case of binary halide–water complexes. Furthermore, these examples present a compelling case that electrostatic explanations based on low-order multipole moments are very often counterfactual for nonbonded interactions at close range and should not be taken seriously in the absence of additional justification.
Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.
Phosphine ligands with pendant amines have played an important role in the development of molecular electrocatalyst for H2 oxidation and production and for other potential energy storage and utilization reactions, such as CO2 and O2 reduction and NH3 oxidation. The pendant amines facilitate substrate binding, heterolytic bond cleavage and formation, intra- and intermolecular proton transfers, the formation of dihydrogen and hydrogen bonds, and provide a physical mechanism for coupling electron and proton transfer reactions. The structures of the phosphine ligands determine the ability of the pendant amines to participate in these different catalytic steps and the stability of the electrocatalysts as they are cycled between different oxidation states. These different features are illustrated by a detailed comparison of Mn, Fe, and Ni electrocatalysts for H2 oxidation that contain diphosphine ligands with pendant amines.
The eco-friendly processing of conjugated polymer binder for lithium-ion batteries demands improved polymer solubility by introducing functional moieties, while this strategy will concurrently sacrifice polymer conductivity. Employing the polyfluorene-based binder poly(2,7-9,9 (di(oxy-2,5,8-trioxadecane))fluorene) (PFO), soluble in water-ethanol mixtures, a novel approach is presented to solve this trade-off, which features integration of aqueous solution processing with subsequent controlled thermal-induced cleavage of solubilizing side chains, to produce hierarchically ordered structures (HOS). The thermal processing can enhance the intermolecular π–π stacking of polyfluorene backbone for better electrochemical performance. Notably, HOS-PFO demonstrated a substantial 6–7 orders of magnitude enhancement in electronic conductivity, showcasing its potential as a functional binder for lithium-ion batteries. As an illustration, HOS-PFO protected SiOx anodes, utilizing in situ side chain decomposition of PFO surrounding SiOx particles after aqueous processing are fabricated. HOS-PFO contributed to the stable cycling and high-capacity retention of practical SiOx anodes (3.0 mAh cm -2 ), without the use of any conducting carbon additives or fluorinated electrolyte additives. It is proposed that this technique represents a universal approach for fabricating electrodes with conjugated polymer binders from aqueous solutions without compromising conductivity.
The essential role of a well-defined hydrogen-bond network in achieving chemically reversible multiproton translocations triggered by one-electron electrochemical oxidation/reduction is investigated by using pyridylbenzimidazole–phenol models. The two molecular architectures designed for these studies differ with respect to the position of the N atom on the pyridyl ring. In one of the structures, a hydrogen-bond network extends uninterrupted across the molecule from the phenol to the pyridyl group. Experimental and theoretical evidence indicates that an overall chemically reversible two-protoncoupled electron-transfer process (E2PT) takes place upon electrochemical oxidation of the phenol. This E2PT process yields the pyridinium cation and is observed regardless of the cyclic voltammogram scan rate. In contrast, when the hydrogen-bond network is disrupted, as seen in the isomer, at high scan rates (~1000 mV s –1 ) a chemically reversible process is observed with an E 1/2 characteristic of a one-proton-coupled electron-transfer process (E1PT). At slow cyclic voltammetric scan rates (<1000 mV s –1 ) oxidation of the phenol results in an overall chemically irreversible two-proton-coupled electron-transfer process in which the second proton-transfer step yields the pyridinium cation detected by infrared spectroelectrochemistry. In this case, we postulate an initial intramolecular proton-coupled electron-transfer step yielding the E1PT product followed by a slow, likely intermolecular chemical step involving a second proton transfer to give the E2PT product. Insights into the electrochemical behavior of these systems are provided by theoretical calculations of the electrostatic potentials and electric fields at the site of the transferring protons for the forward and reverse processes. Furthermore, this work addresses a fundamental design principle for constructing molecular wires where protons are translocated over varied distances by a Grotthuss-type mechanism.