Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERMOLECULAR FORCES”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

High temperature spin dynamics in linear magnetic chains, molecular rings, and segments by nuclear magnetic resonance

We present the room temperature proton nuclear magnetic resonance (NMR) nuclear spin-lattice relaxation rate (NSLR) results in two 1D spin chains: the Heisenberg antiferromagnetic (AFM) Eu(hfac){sub 3}NITEt and the magnetically frustrated Gd(hfac){sub 3}NITEt. The NSLR as a function of external magnetic field can be interpreted very well in terms of high temperature spin dynamics dominated by a long time persistence of the decay of the two-spin correlation function due to the conservation of the total spin value for isotropic Heisenberg chains. The high temperature spin dynamics are also investigated in Heisenberg AFM molecular rings. In both Cr{sub 8} closed ring and in Cr{sub 7}Cd and Cr{sub 8}Zn open rings, i.e., model systems for a finite spin segment, an enhancement of the low frequency spectral density is found consistent with spin diffusion but the high cut-off frequency due to intermolecular anisotropic interactions prevents a detailed analysis of the spin diffusion regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Aqueous-Based Polyimine Functionalization of Cellulose Nanofibrils for Effective Drying and Polymer Composite Reinforcement

Here a simple strategy was developed to synthesize polyimine-coated cellulose nanofibrils (CNF) for effective CNF drying and composite reinforcement. The polyimine was synthesized in an aqueous medium using a selective hydrophilic and hydrophobic component that forces the polyimine to precipitate, which prevents the reverse imine reaction. The polyimine coating allowed the CNF to be easily oven-dried while maintaining a fibrillar morphology to provide mechanical reinforcement in poly(ethylene terephthalate glycol) (PETG) composites. In comparison, poor mechanical performance and a heterogeneous fracture surface were observed when the coated CNF were incorporated into poly(l-lactide) (PLA) composites. It is hypothesized that intermolecular aromatic–aromatic interactions are formed at the interface between fibers and the polymer matrix in the PETG system, while no such phenomena occur in the PLA system. Overall, this facile strategy to produce modified, easily dried CNF can be adapted to produce polyimine-coated CNF with unique functionalities that are useful in a range of applications.

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

Cavitation energies can outperform dispersion interactions

Amongst all intermolecular interactions that bring molecules, biomolecules, and even molecular fragments together, dispersion interactions are the omnipresent ones. Nevertheless, their contribution to solution-phase phenomena has recently become a matter of debate. In aqueous solution, in particular, they compete with the hydrophobic effect, which also drives associative processes, and when macrocyclic cavities are involved, the release of high-energy water has been identified as an extra driving force. We have now studied the binding of noble gases (He, Ne, Ar, Kr, Xe) with the macrocyclic receptor cucurbit[5]uril (CB5) in aqueous solution through a displacement strategy by using methane and ethane as 1H NMR probes. We dissected the hydration free energies of the noble gases into an attractive dispersive contribution and a repulsive one for cavity formation. This allowed us to fine-grain the components that contribute to host-guest binding and we reach the conclusion that the binding process in solution is driven by differential cavitation energies rather than dispersion interactions. This means that the free energy required to create a cavity to accept the noble gas inside CB5 is much lower (actually zero, since the cavity of the host is non-solvated) than that to create a similarly sized cavity in bulk water, and it is the recovery of the latter cavitation energy, which drives the overall process. The concept of differential cavitation energies merges seamlessly two previous approaches for predicting supramolecular thermodynamics, namely that of cohesive energy density (which applies to bulk solvent) and that of high-energy water (which becomes important for discrete concave binding sites). The dominance of differential cavitation energies in driving guest binding to poorly hydrated, concave binding sites (V < 100 Å3) has implications for the refinement of gas-storage materials and the understanding of biological receptors.

He, Suhang↗
DOE OSTI2014

Production of Hydrogen by Electrocatalysis: Making the H-H Bond by Combining Protons and Hydrides

Generation of hydrogen by reduction of two protons by two electrons can be catalysed by molecular electrocatalysts. Determination of the thermodynamic driving force for elimination of H2 from molecular complexes is important for the rational design of molecular electrocatalysts, and allows the design of metal complexes of abundant, inexpensive metals rather than precious metals (“Cheap Metals for Noble Tasks”). The rate of H2 evolution can be dramatically accelerated by incorporating pendant amines into diphosphine ligands. These pendant amines in the second coordination sphere function as protons relays, accelerating intramolecular and intermolecular proton transfer reactions. The thermodynamics of hydride transfer from metal hydrides and the acidity of protonated pendant amines (pKa of N-H) contribute to the thermodynamics of elimination of H2; both of the hydricity and acidity can be systematically varied by changing the substituents on the ligands. A series of Ni(II) electrocatalysts with pendant amines have been developed. In addition to the thermochemical considerations, the catalytic rate is strongly influenced by the ability to deliver protons to the correct location of the pendant amine. Protonation of the amine endo to the metal leads to the N-H being positioned appropriately to favor rapid heterocoupling with the M-H. Designing ligands that include proton relays that are properly positioned and thermodynamically tuned is a key principle for molecular electrocatalysts for H2 production as well as for other multi-proton, multi-electron reactions important for energy conversions. The research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences. Pacific Northwest National Laboratory is operated by Battelle for DOE.

catalysis↗
DOE OSTI2016

FOB-SH: Fragment orbital-based surface hopping for charge carrier transport in organic and biological molecules and materials

We introduce a fragment orbital-based fewest switches surface hopping method, FOB-SH, designed to efficiently simulate charge carrier transport in strongly fluctuating condensed phase systems such as organic semiconductors and biomolecules. The charge carrier wavefunction is expanded and the electronic Hamiltonian constructed in a set of singly occupied molecular orbitals of the molecular sites that mediate the charge transfer. Diagonal elements of the electronic Hamiltonian (site energies) are obtained from a force field, whereas the off-diagonal or electronic coupling matrix elements are obtained using our recently developed analytic overlap method. We derive a general expression for the exact forces on the adiabatic ground and excited electronic state surfaces from the nuclear gradients of the charge localized electronic states. Applications to electron hole transfer in a model ethylene dimer and through a chain of ten model ethylenes validate our implementation and demonstrate its computational efficiency. On the larger system, we calculate the qualitative behaviour of charge mobility with change in temperature T for different regimes of the intermolecular electronic coupling. For small couplings, FOB-SH predicts a crossover from a thermally activated regime at low temperatures to a band-like transport regime at higher temperatures. For higher electronic couplings, the thermally activated regime disappears and the mobility decreases according to a power law. This is interpreted by a gradual loss in probability for resonance between the sites as the temperature increases. The polaron hopping model solved for the same system gives a qualitatively different result and underestimates the mobility decay at higher temperatures. Taken together, the FOB-SH methodology introduced here shows promise for a realistic investigation of charge carrier transport in complex organic, aqueous, and biological systems.

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

Intermolecular symmetry-adapted perturbation theory study of large organic complexes

Binding energies for the complexes of the S12L database by Grimme [Chem. Eur. J. 18, 9955 (2012)] were calculated using intermolecular symmetry-adapted perturbation theory combined with a density-functional theory description of the interacting molecules. The individual interaction energy decompositions revealed no particular change in the stabilisation pattern as compared to smaller dimer systems at equilibrium structures. This demonstrates that, to some extent, the qualitative description of the interaction of small dimer systems may be extrapolated to larger systems, a method that is widely used in force-fields in which the total interaction energy is decomposed into atom-atom contributions. A comparison of the binding energies with accurate experimental reference values from Grimme, the latter including thermodynamic corrections from semiempirical calculations, has shown a fairly good agreement to within the error range of the reference binding energies.

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

Liquid chloroform structure from computer simulation with a full ab initio intermolecular interaction potential

We have calculated the intermolecular interaction energies of the chloroform dimer in 12 orientations using the second-order Møller-Plesset perturbation theory. Single point energies of important geometries were calibrated by the coupled cluster with single and double and perturbative triple excitation method. Dunning's correlation consistent basis sets up to aug-cc-pVQZ have been employed in extrapolating the interaction energies to the complete basis set limit values. With the ab initio potential data we constructed a 5-site force field model for molecular dynamics simulations. We compared the simulation results with recent experiments and obtained quantitative agreements for the detailed atomwise radial distribution functions. Our results were also consistent with previous results using empirical force fields with polarization effects. Moreover, the calculated diffusion coefficients reproduced the experimental data over a wide range of thermodynamic conditions. To the best of our knowledge, this is the first ab initio force field which is capable of competing with existing empirical force fields for liquid chloroform.

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

Communication: Contrasting effects of glycerol and DMSO on lipid membrane surface hydration dynamics and forces

Glycerol and dimethyl sulfoxide (DMSO) are commonly used cryoprotectants in cellular systems, but due to the challenges of measuring the properties of surface-bound solvent, fundamental questions remain regarding the concentration, interactions, and conformation of these solutes at lipid membrane surfaces. We measured the surface water diffusivity at gel-phase dipalmitoylphosphatidylcholine (DPPC) bilayer surfaces in aqueous solutions containing ≤7.5 mol. % of DMSO or glycerol using Overhauser dynamic nuclear polarization. We found that glycerol similarly affects the diffusivity of water near the bilayer surface and that in the bulk solution (within 20%), while DMSO substantially increases the diffusivity of surface water relative to bulk water. We compare these measurements of water dynamics with those of equilibrium forces between DPPC bilayers in the same solvent mixtures. DMSO greatly decreases the range and magnitude of the repulsive forces between the bilayers, whereas glycerol increases it. We propose that the differences in hydrogen bonding capability of the two solutes leads DMSO to dehydrate the lipid head groups, while glycerol affects surface hydration only as much as it affects the bulk water properties. The results suggest that the mechanism of the two most common cryoprotectants must be fundamentally different: in the case of DMSO by decoupling the solvent from the lipid surface, and in the case of glycerol by altering the hydrogen bond structure and intermolecular cohesion of the global solvent, as manifested by increased solvent viscosity.

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

Non-covalent interactions of nitrous oxide with aromatic compounds: Spectroscopic and computational evidence for the formation of 1:1 complexes

We present the first study of intermolecular interactions between nitrous oxide (N{sub 2}O) and three representative aromatic compounds (ACs): phenol, cresol, and toluene. The infrared spectroscopic experiments were performed in a Ne matrix and were supported by high-level quantum chemical calculations. Comparisons of the calculated and experimental vibrational spectra provide direct identification and characterization of the 1:1 N{sub 2}O-AC complexes. Our results show that N{sub 2}O is capable of forming non-covalently bonded complexes with ACs. Complex formation is dominated by dispersion forces, and the interaction energies are relatively low (about −3 kcal mol{sup −1}); however, the complexes are clearly detected by frequency shifts of the characteristic bands. These results suggest that N{sub 2}O can be bound to the amino-acid residues tyrosine or phenylalanine in the form of π complexes.

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

Communication: THz absorption spectrum of the CO{sub 2}–H{sub 2}O complex: Observation and assignment of intermolecular van der Waals vibrations

Terahertz absorption spectra have been recorded for the weakly bound CO{sub 2}–H{sub 2}O complex embedded in cryogenic neon matrices at 2.8 K. The three high-frequency van der Waals vibrational transitions associated with out-of-plane wagging, in-plane rocking, and torsional motion of the isotopic H{sub 2}O subunit have been assigned and provide crucial observables for benchmark theoretical descriptions of this systems’ flat intermolecular potential energy surface. A (semi)-empirical value for the zero-point energy of 273 ± 15 cm{sup −1} from the class of intermolecular van der Waals vibrations is proposed and the combination with high-level quantum chemical calculations provides a value of 726 ± 15 cm{sup −1} for the dissociation energy D{sub 0}.

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

Uncovering the mechanism of selective stabilization of high-energy diastereoisomers via inclusion

Abstract Supramolecular systems may be used to stabilize otherwise unstable isomers to find alternative synthetic pathways. It has been reported that cucurbit[8]uril can stabilize trans -I and trans -II Cu $$^{{\textrm{II}}}$$ II cyclam, whereas trans -III is the only non-substituted trans Cu $$^{{\textrm{II}}}$$ II cyclam diastereoisomer found outside of the host molecule experimentally. Quantum chemistry methods can provide valuable insight into the intermolecular interactions involved in these inclusion complexes. All five possible trans diastereoisomers of Cu $$^{{\textrm{II}}}$$ II cyclam were studied within the host molecule to calculate the interaction energy and free energy of association for each complex. The relative free energies of the five free cyclams confirm that trans -I and trans -II are the most energetically accessible diastereoisomers from the initial trans -III starting point. Energy decomposition analysis was used to identify the attractive and repulsive interactions between cyclam and cucurbit[8]uril and showed that trans -II encounters repulsive forces almost three times greater than trans -I, which may explain the 7:3 ratio of trans -I to trans -II within cucurbit[8]uril that occurs experimentally. Optimized complex geometries with trans -III, IV, and V show that the cyclams protrude out of cucurbit[8]uril, whereas trans -I and trans -II become more encapsulated and elongate the host, suggesting that the position of the cyclam is extremely important when forming non-covalent interactions. Our results agree with the experimental findings and provide greater insight into why the most stable isolated cyclam diastereoisomer, trans -III, does not form a complex.

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

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The mechanism of chemisorption of hydrogen atom on graphene: Insights from the reaction force and reaction electronic flux

At the PBE-D3/cc-pVDZ level of theory, the hydrogen chemisorption on graphene was analyzed using the reaction force and reaction electronic flux (REF) theories in combination with electron population analysis. It was found that chemisorption energy barrier is mainly dominated by structural work (∼73%) associated to the substrate reconstruction whereas the electronic work is the greatest contribution of the reverse energy barrier (∼67%) in the desorption process. Moreover, REF shows that hydrogen chemisorption is driven by charge transfer processes through four electronic events taking place as H approaches the adsorbent surface: (a) intramolecular charge transfer in the adsorbent surface; (b) surface reconstruction; (c) substrate magnetization and adsorbent carbon atom develops a sp{sup 3} hybridization to form the σC-H bond; and (d) spontaneous intermolecular charge transfer to reach the final chemisorbed state.

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

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
DOE OSTI2016

Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions

An accurate treatment of the long-range electron correlation energy, including van der Waals (vdW) or dispersion interactions, is essential for describing the structure, dynamics, and function of a wide variety of systems. Among the most accurate models for including dispersion into density functional theory (DFT) is the range-separated many-body dispersion (MBD) method [A. Ambrosetti et al., J. Chem. Phys., 2014, 140, 18A508], in which the correlation energy is modeled at short-range by a semi-local density functional and at long-range by a model system of coupled quantum harmonic oscillators. In this work, we develop analytical gradients of the MBD energy with respect to nuclear coordinates, including all implicit coordinate dependencies arising from the partitioning of the charge density into Hirshfeld effective volumes. To demonstrate the efficiency and accuracy of these MBD gradients for geometry optimizations of systems with intermolecular and intramolecular interactions, we optimized conformers of the benzene dimer and isolated small peptides with aromatic side-chains. We find excellent agreement with the wavefunction theory reference geometries of these systems (at a fraction of the computational cost) and find that MBD consistently outperforms the popular TS and D3(BJ) dispersion corrections. To demonstrate the performance of the MBD model on a larger system with supramolecular interactions, we optimized the C 60 @C 60 H 28 buckyball catcher host–guest complex. In our analysis, we also find that neglecting the implicit nuclear coordinate dependence arising from the charge density partitioning, as has been done in prior numerical treatments, leads to an unacceptable error in the MBD forces, with relative errors of ~20% (on average) that can extend well beyond 100%.

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

Toward transferable interatomic van der Waals interactions without electrons: The role of multipole electrostatics and many-body dispersion

We estimate polarizabilities of atoms in molecules without electron density, using a Voronoi tesselation approach instead of conventional density partitioning schemes. The resulting atomic dispersion coefficients are calculated, as well as many-body dispersion effects on intermolecular potential energies. We also estimate contributions from multipole electrostatics and compare them to dispersion. We assess the performance of the resulting intermolecular interaction model from dispersion and electrostatics for more than 1300 neutral and charged, small organic molecular dimers. Applications to water clusters, the benzene crystal, the anti-cancer drug ellipticine—intercalated between two Watson-Crick DNA base pairs, as well as six macro-molecular host-guest complexes highlight the potential of this method and help to identify points of future improvement. The mean absolute error made by the combination of static electrostatics with many-body dispersion reduces at larger distances, while it plateaus for two-body dispersion, in conflict with the common assumption that the simple 1/R{sup 6} correction will yield proper dissociative tails. Overall, the method achieves an accuracy well within conventional molecular force fields while exhibiting a simple parametrization protocol.

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

Scanning tunneling microscopy/spectroscopy of picene thin films formed on Ag(111)

Using ultrahigh-vacuum low-temperature scanning tunneling microscopy and spectroscopy combined with first principles density functional theory calculations, we have investigated structural and electronic properties of pristine and potassium (K)-deposited picene thin films formed in situ on a Ag(111) substrate. At low coverages, the molecules are uniformly distributed with the long axis aligned along the [112{sup ¯}] direction of the substrate. At higher coverages, ordered structures composed of monolayer molecules are observed, one of which is a monolayer with tilted and flat-lying molecules resembling a (11{sup ¯}0) plane of the bulk crystalline picene. Between the molecules and the substrate, the van der Waals interaction is dominant with negligible hybridization between their electronic states; a conclusion that contrasts with the chemisorption exhibited by pentacene molecules on the same substrate. We also observed a monolayer picene thin film in which all molecules were standing to form an intermolecular π stacking. Two-dimensional delocalized electronic states are found on the K-deposited π stacking structure.

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

Characteristics of energy exchange between inter- and intramolecular degrees of freedom in crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) with implications for coarse-grained simulations of shock waves in polyatomic molecular crystals

In this report, we characterize the kinetics and dynamics of energy exchange between intramolecular and intermolecular degrees of freedom (DoF) in crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). All-atom molecular dynamics (MD) simulations are used to obtain predictions for relaxation from certain limiting initial distributions of energy between the intra- and intermolecular DoF. The results are used to parameterize a coarse-grained Dissipative Particle Dynamics at constant Energy (DPDE) model for TATB. Each TATB molecule in the DPDE model is represented as an all-atom, rigid-molecule mesoparticle, with explicit external (molecular translational and rotational) DoF and coarse-grained implicit internal (vibrational) DoF. In addition to conserving linear and angular momentum, the DPDE equations of motion conserve the total system energy provided that particles can exchange energy between their external and internal DoF. The internal temperature of a TATB molecule is calculated using an internal equation of state, which we develop here, and the temperatures of the external and internal DoF are coupled using a fluctuation-dissipation relation. The DPDE force expression requires specification of the input parameter σ that determines the rate at which energy is exchanged between external and internal DoF. We adjusted σ based on the predictions for relaxation processes obtained from MD simulations. The parameterized DPDE model was employed in large-scale simulations of shock compression of TATB. We show that the rate of energy exchange governed by σ can significantly influence the transient behavior of the system behind the shock.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗