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.
Dynamic intermolecular forces through hydrogen bonding between water and polymer chains are shown to enhance thermal transport in soft materials.
This work describes the use of C–H∙∙∙F–C contacts in the solid-state from the stator towards the rotator to fine-tune their internal motion, by constructing a set of interactions that generate close-fitting cavities in three supramolecular rotors 1–3I. The crystal structures of these rotors, determined by synchrotron radiation experiments at different temperatures, show the presence of such C–H∙∙∙F–C contacts between extended carbazole stators featuring fluorinated phenyl rings and the 1,4-diazabicyclo[2.2.2]octane (DABCO) rotator. According to the 2 H NMR results, using deuterated samples, and periodic density functional theory computations, the rotators experience fast angular displacements (preferentially 120° jumps) due to their low rotational activation energies (E a = 0.8–2.0 kcal mol -1 ). The higher rotational barrier for 1 (2.0 kcal mol -1 ) is associated with a larger number of weak C–H∙∙∙F–C contacts generated by the stators. This strategy offers the possibility to explore the correlation among weak intermolecular forces, cavity shape, and internal dynamics, which has strong implications in the design of future fine-tuned amphidynamic crystals.
Liquid emulsion droplet evaporation is of importance for various sensing and imaging applications. The liquid-to-gas phase transformation is typically triggered thermally or acoustically by low–boiling point liquids, or by inclusion of solid structures that pin the vapor/liquid contact line to facilitate heterogeneous nucleation. However, these approaches lack precise tunability in vaporization behavior. Here, we describe a previously unused approach to control vaporization behavior through an endoskeleton that can melt and blend into the liquid core to either enhance or disrupt cohesive intermolecular forces. This effect is demonstrated using perfluoropentane (C 5 F 12 ) droplets encapsulating a fluorocarbon (FC) or hydrocarbon (HC) endoskeleton. FC skeletons inhibit vaporization, whereas HC skeletons trigger vaporization near the rotator melting transition. Our findings highlight the importance of skeletal interfacial mixing for initiating droplet vaporization. Tuning molecular interactions between the endoskeleton and droplet phase is generalizable for achieving emulsion or other secondary phase transitions, in emulsions.
An anti-viral coating for the nose & throat is needed to quickly address the spread of COVID-19 infections and to aid future pandemics. Current nasal delivery systems are typically a 1-spray homogenous solution, which is convenient but may not be as efficient or durable as a multi-spray solution that contains a binding layer to anchor the anti-viral components. Many of the current nasal solutions being investigated to aid in pandemic response have solution-based mechanisms of action and are quickly cleared from the nose/mouth limiting their efficacy lifetime. A multi-spray approach utilizing strong intermolecular forces between polymeric materials and anti-viral agents to provide a robust biocompatible coating is expected to have increased physical and chemical properties to combat viral infection. This work evaluated which tailored biocompatible materials are needed to generate a layered system to combat viral infection. This layered system takes advantage of electrostatic interactions to bind anti-viral components, shown in Figure 1. Proof of concept success was demonstrated through cell toxicity studies and anti-viral assays on both the individual layer components and the complete layered system. Project work began with chemical modifications of Xantham Gum to hydrolyzing xanthan gum to give a negatively charged polymer with varying degrees of ionic character and investigate bonding efficiencies. However, this was abandoned in favor of comparing COTS available materials listed in Table 1 (Results section). Concentration and adhesion studies were performed on layers created using the listed polymeric materials and anti-viral agents. Finally, biocompatibility of materials, layers, and coating system were confirmed through cytotoxicity studies and the efficacy of the anti-viral properties tested with phi6 bacteriophage as a SARS-CoV2 surrogate.
Six hundred and seventy-five measurements of dynamic viscosity and density have been used to assess the prediction error of the Arrhenius blending rule for kinematic viscosity of hydrocarbon mixtures. Major trends within the data show that mixture complexity–binary to hundreds of components—and temperature are more important determinants of prediction error than differences in molecular size or hydrogen saturation between the components of the mixtures. Over the range evaluated, no correlation between prediction error and mole fractions was observed, suggesting the log of viscosity truly is linear in mole fraction, as indicated by the Arrhenius blending rule. Mixture complexity and temperature also impact molar volume and its prediction. However, a linear regression between the two model errors explains less than 20% of the observed variation, indicating that mixture viscosity and/or molar volume are not linear with respect to temperature and/or mixture complexity. Extensive discussion of the intermolecular forces and the geometric arrangement of molecules and vacancies in liquids, which ultimately determines its viscosity, is brought into context with the implicit approximations behind the Arrhenius blending rule. The complexity of this physics is not compatible with a simple algebraic correction to the model. However, sufficient data is now available to determine confidence intervals around the prediction of fuel viscosity based on its component mole fractions and viscosities. At -40°C, when all identified components are pure molecules the modeling error is 13.2% of the predicted (nominal) viscosity times the root mean square of the component mole fractions.
Here we introduce the synthesis and thermal properties of a series of sophisticated imidazolium ionenes with alternating amide-amide or amide-imide backbone functionality, and investigate the structural effects of mono(imidazolium) and unprecedented tris(imidazolium) ionic liquids (ILs) in these ionenes. The new set of poly(amide-amide) (PAA) and poly(amide-imide) (PAI) ionenes represent the intersection of conventional high-performance polymers with the ionene archetype–presenting polymers with alternating functional and ionic elements precisely sequenced along the backbone. The effects of polymer composition on the thermal properties and morphology were analyzed. Five distinct polymer backbones were synthesized and combined with a stoichiometric equivalent of the IL 1-benzyl-3-methylimidazolium bistriflimide ([Bnmim][Tf 2 N]), which were studied to probe the self-assembly, structuring, and contributions of intermolecular forces when IL is added. Furthermore, three polyamide (PA) or polyimide (PI) ionenes with simpler xylyl linkages were interfaced with [Bnmim][Tf 2 N] as well as a novel amide-linked tris(imidazolium) IL, to demonstrate the structural changes imparted by the inclusion of functional, ionic additives dispersed within the ionene matrix. This work highlights the possibilities for utilizing concepts from small molecules which exhibit supramolecular self-assembly to guide creative design and manipulate the structuring of ionenes.
The quest for converting lignin into high-value products has been continuously pursued in the past few decades. In its native form, lignin is a group of heterogeneous polymers comprised of phenylpropanoids. The major commercial lignin streams, including Kraft lignin, lignosulfonates, soda lignin and organosolv lignin, are produced from industrial processes including the paper and pulping industry and emerging lignocellulosic biorefineries. Although lignin has been viewed as a low-cost and renewable feedstock to replace petroleum-based materials, its utilization in polymeric materials has been suppressed due to the low reactivity and inherent physicochemical properties of lignin. Hence, various lignin modification strategies have been developed to overcome these problems. Herein, we review recent progress made in the utilization of functionalized lignins in commodity polymers including thermoset resins, blends/composites, grafted functionalized copolymers and carbon fiber precursors. In the synthesis of thermoset resins such as polyurethane, phenol-formaldehyde and epoxy, they are covalently incorporated into the polymer matrix, and the discussion is focused on chemical modifications improving the reactivity of technical lignins. In blends/composites, functionalization of technical lignins is based upon tuning the intermolecular forces between polymer components. In addition, grafted functional polymers have expanded the utilization of lignin-based copolymers to biomedical materials and value-added additives. Different modification approaches have also been applied to facilitate the application of lignin as carbon fiber precursors, heavy metal adsorbents and nanoparticles. These emerging fields will create new opportunities in cost-effectively integrating the lignin valorization into lignocellulosic biorefineries.
Abstract. Adsorption and desorption of gases on liquid or solid substrates are involved in multiphase processes and heterogeneous chemical reactions. The desorption energy (Edes0), which depends on the intermolecular forces between adsorbate and substrate, determines the residence time of chemical species at interfaces. We show how Edes0 and temperature influence the net uptake or release of gas species, the rates of surface–bulk exchange and surface or bulk reactions, and the equilibration timescales of gas–particle partitioning. Using literature data, we derive a parameterization to estimate Edes0 for a wide range of chemical species based on the molecular mass, polarizability, and oxygen-to-carbon ratio of the desorbing species independent of substrate-specific properties, which is possible because of the dominant role of the desorbing species' properties. Correlations between Edes0 and the enthalpies of vaporization and solvation are rooted in molecular interactions. The relation between Edes0 and desorption kinetics reflects the key role of interfacial exchange in multiphase processes. For small molecules and semi-volatile organics (VOC, IVOC, SVOC), Edes0 values around 10–100 kJ mol−1 correspond to desorption lifetimes around nanoseconds to days at room temperature. Even higher values up to years are obtained at low temperatures and for low volatile organic compounds (LVOC, ELVOC/ULVOC) relevant for secondary organic aerosols (SOA). Implications are discussed for SOA formation, gas–particle partitioning, organic phase changes, and indoor surface chemistry. We expect these insights to advance the mechanistic and kinetic understanding of multiphase processes in atmospheric and environmental physical chemistry, aerosol science, materials science, and chemical engineering.
Recent experiments suggest that vibrational strong coupling (VSC) may significantly modify ground-state chemical reactions and their rates even without external pumping. The intrinsic mechanism of this “vacuum-field catalysis” remains largely unclear. Generally, modifications of thermal reactions in the ground electronic states can be caused by equilibrium or non-equilibrium effects. The former are associated with modifications of the reactant equilibrium distribution as expressed by the transition state theory of chemical reaction rates, while the latter stem from the dynamics of reaching and leaving transition state configurations. In this work, we examine how VSC can affect chemical reactions rates in a cavity environment according to transition state theory. Our approach is to examine the effect of coupling to cavity mode(s) on the potential of mean force (PMF) associated with the reaction coordinate. Within the context of classical nuclei and classical photons and also assuming no charge overlap between molecules, we find that while the PMF can be affected by the cavity environment, this effect is negligible for the usual micron-length cavities used to examine VSC situations.
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.
Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.
In their recent communication a semi-local density functional derived from the density matrix expansion of the exchange hole localised by a general coordinate transformation. We show that the order-of-limits problem present in the functional, dismissed as harmless in the original publication, causes severe errors in predicted phase transition pressures. We also show that the claim that lattice volume prediction accuracy exceeds that of existing similar functionals was based on comparison to reference data that misses anharmonic zero-point expansion and consequently overestimates accuracy. Here, by highlighting these omissions, we give a more accurate assessment of the Tao-Mo functional and show a possible direction for resolving the order-of-limits problem.
In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.
In this paper, we present an investigation on the absorption of alkylthiolated nanocrystals on a solid substrate. We calculate adsorption free energies and report a number of effects induced by the substrate. Nearest neighbor distances and bonding free energies are significantly different than for a free floating case, there is a weakening of bonding free energies among nanocrystals, and the adsorption is manifestly anisotropic, i.e., stronger along certain directions of the nanocrystal core. We contend that this last result accounts for the Bain transition (fcc → bcc) observed in experimental results. We report the presence of vortices induced by the substrate, which explain the increased nearest neighbor distance among nanocrystals, which is in excellent quantitative agreement with experimental results and with the predictions of the Orbifold Topological Model. Implications for the assembly of nanostructures and future experiments are also discussed.
Developing predictive thermal property models for liquids based on microscopic principles has been elusive. The difficulty is that liquids have gas-like and solid-like attributes that are at odds when considering the frameworks of microscopic models: Models for gases are simple due to randomness and low density, whereas models for crystalline solids rely on symmetry and long-range order for easier calculation. The short-range order in liquids does, however, provide structure to neighboring molecules similar to amorphous solids, and there have been recent advances indicating that collective vibrational modes store heat in liquids. Models combining Debye approximations from solid-state physics and Frenkel’s theory of liquids can accurately predict the heat capacity of liquids. Phonon-like dispersions in liquids have also been widely observed in neutron scattering experiments. These developments motivate us to propose a model where high-frequency vibrational modes, which travel at the speed of sound and have a mean free path on the order of the average intermolecular distance, conduct heat in liquids. We use this liquid phonon gas model to calculate the thermal conductivity of liquids with varying intermolecular interaction energies from strongest to weakest—Coulomb, hydrogen-bonding, Keesom, and London dispersion energy. Generally, the model is more accurate as the intermolecular interaction energy and density of liquids increase. The calculated thermal conductivity of Coulombic-bound molten sodium nitrate and hydrogen-bonded water is within 1.46% and 2.98% of the experimentally measured values, respectively, across their entire temperature ranges. Further modal analysis of the velocity and the mean free path of collective vibrations could establish the liquid phonon gas model as an accurate model for weakly interacting liquids as well.
Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.
Materials with high dielectric constant, ε s , are desirable in a wide range of applications including energy storage and actuators. Recently, zwitterionic liquids have been reported to have the largest ε s of any liquid and, thus, have the potential to replace inorganic fillers to modulate the material ε s . Although the large ε s for zwitterionic liquids is attributed to their large molecular dipole, the role of chemical substituents attached to the zwitterion cation on ε s is not fully understood, which is necessary to enhance the performance of soft energy materials. Here, we report the impact of zwitterionic liquid cation chemical substituents on ε s (50 < ε s < 300 at room temperature). Dielectric relaxation spectroscopy reveals that molecular reorientation is the main contributor to the high ε s . The low Kirkwood factor g calculated for zwitterionic liquids (e.g., 0.1–0.2) suggests the tendency for the antiparallel zwitterion dipole alignment expected from the strong electrostatic intermolecular interactions. With octyl cation substituents, the g is decreased due to the formation of hydrophobic-rich domains that restrict molecular reorientation under applied electric fields. In contrast, when zwitterion cations are functionalized with ethylene oxide (EO) segments, g increases due to the EO segments interacting with the cations, allowing more zwitterion rotation in response to the applied field. The reported results suggest that high ε s zwitterionic liquids require a large molecular dipole, compositionally homogeneous liquids (e.g., no aggregation), a maximized zwitterion number density, and a high g, which is achievable by incorporating polar chemical substituents onto the zwitterion cations.
Molecular Dynamics (MD) simulations based on the Effective Fragment Potential (EFP) method are utilized to provide a comprehensive assessment of diffusion in liquid n-hexane. We decompose translational diffusion into components along and orthogonal to the long axis of the molecule. Rotational diffusion is decomposed into tumbling and spinning motions about this axis. Our analysis yields four corresponding diffusion coefficients which are related to diagonal entries in the complete 6 × 6 diffusion tensor accounting for the three rotational and three translational degrees of freedom and for the potential coupling between them. However, coupling between different degrees of freedom is expected to be minimal for a natural choice of the molecular body-fixed axis, so then off-diagonal entries in the tensor are negligible. This expectation is supported by a hydrodynamic analysis of the diffusion tensor which treats the liquid surrounding the molecule being tracked as a viscous continuum. Thus, the EFP MD analysis provides a comprehensive characterization of diffusion and also reveals expected shortcomings of the hydrodynamic treatment, particularly for rotational diffusion, when applied to neat liquids.