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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

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

Linearly dispersive bands at the onset of correlations in K 𝑥 ⁢C 60 films

Molecular crystals are a flexible platform to induce novel electronic phases. Due to the weak forces between molecules, intermolecular distances can be varied over larger ranges than interatomic distances in atomic crystals. On the other hand, the hopping terms are generally small, which results in narrow bands, strong correlations, and heavy electrons. Here, by growing K x C 60 fullerides on hexagonal layered Bi 2 Se 3 , we show that upon doping the series undergoes a Mott transition from a molecular insulator to a correlated metal and an in-gap state evolves into highly dispersive Dirac-like fermions at half filling, where superconductivity occurs. This picture challenges the commonly accepted description of the low-energy quasiparticles as appearing from a gradual electron doping of the conduction states and suggests an intriguing parallel with the more famous family of the cuprate superconductors. More in general, it indicates that molecular crystals offer a viable route to engineer electron-electron interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1974

Interaction between polymer constituents and the structure of biopolymers

The paper reviews the current status of methods for calculating intermolecular interactions between biopolymer units. The nature of forces contributing to the various domains of intermolecular separations is investigated, and various approximations applicable in the respective regions are examined. The predictive value of current theory is tested by establishing a connection with macroscopic properties and comparing the theoretical predicted values with those derived from experimental data. This has led to the introduction of a statistical model describing DNA.

Rein, R.↗
DOE OSTI2021

Neat, Simple, and Wrong: Debunking Electrostatic Fallacies Regarding Noncovalent Interactions

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.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Acoustic Mismatch Model for Thermal Contact Conductance of Van Der Waals Contacts Under Static Force

We report that Van der Waals interfaces play a major role in technology today. Thermal transport in material systems with van der Waals interfaces is mainly limited by the contact conductance. Although the effects of static force, such as pressure or the electrostatic part of hydrogen bonds, on the thermal contact conductance of van der Waals interfaces have been examined in a few studies, the focus was either on numerical simulation using techniques such as molecular dynamics or on experimental investigation. In this article, an analytical model of thermal contact conductance that accounts for the effects of static force and adhesion energy is presented. The application of static forces is found to cause a decrease in the intermolecular distance, which leads to increased interatomic forces across the interfaces and thus higher thermal conductance. The model is in good agreement with experimental data on the effect of pressure on thermal conductance collected by Gotsmann and Lantz (Nature Materials, Vol. 12, p. 59–65, 2012).

36 MATERIALS SCIENCE↗
DOE OSTI2018

The Water Dimer II: Theoretical Investigations

As the archetype of aqueous hydrogen bonding, the water dimer has been extensively studied by both theory and experiment for nearly seven decades. Here in this article, we present a detailed chronological review of the theoretical advances made using electronic structure methods to address the structure, hydrogen bonding and vibrational spectroscopy of the water dimer, as well as the role of its potential energy surface in the development of classical force fields to describe intermolecular interaction in clusters and the condensed phases of water.

08 HYDROGEN↗
DOE OSTI2018

Perspective: Ab initio force field methods derived from quantum mechanics

It is often desirable to accurately and efficiently model the behavior of large molecular systems in the condensed phase (thousands to tens of thousands of atoms) over long time scales (from nanoseconds to milliseconds). In these cases, ab initio methods are difficult due to the increasing computational cost with the number of electrons. A more computationally attractive alternative is to perform the simulations at the atomic level using a parameterized function to model the electronic energy. Many empirical force fields have been developed for this purpose. However, the functions that are used to model interatomic and intermolecular interactions contain many fitted parameters obtained from selected model systems, and such classical force fields cannot properly simulate important electronic effects. Furthermore, while such force fields are computationally affordable, they are not reliable when applied to systems that differ significantly from those used in their parameterization. They also cannot provide the information necessary to analyze the interactions that occur in the system, making the systematic improvement of the functional forms that are used difficult. Ab initio force field methods aim to combine the merits of both types of methods. The ideal ab initio force fields are built on first principles and require no fitted parameters. Ab initio force field methods surveyed in this perspective are based on fragmentation approaches and intermolecular perturbation theory. This perspective summarizes their theoretical foundation, key components in their formulation, and discusses key aspects of these methods such as accuracy and formal computational cost. The ab initio force fields considered here were developed for different targets, and this perspective also aims to provide a balanced presentation of their strengths and shortcomings. Finally, this perspective suggests some future directions for this actively developing area.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF–Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born–Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. In this study, using tetrahydrofuran (THF)–water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D 8 THF to examine the molecular origins of the isotope-dependent loss of the THF–water miscibility gap.

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

Tradeoff between Intermolecular Interaction and Backbone Disorder by High Molecular Dipole Block for Improving Blend Morphology of Polymer Solar Cells

Constructing terpolymer has attracted increasing attention as a strategy to improve the performance of polymer solar cell. Terpolymer usually offers an opportunity to lower the frontier molecular orbital energy level, introduces additional absorption band and sometimes optimizes the morphology of the active blend. Generally, the additional segment in terpolymer backbone inevitably introduces backbone disorder, which causes entropy rises. However, selecting a suitable dipole unit introduces extra driving forces for crystallization by enhancing intermolecular interactions. This provides a handy knob for tradeoff between intermolecular interaction and backbone disorder, thus regulating the blend morphology. In this study, a high dipole and electron-deficient group of pyrrolo[3,4-f ]benzotriazole-5,7(6H)-dione (TzBI) is incorporated into the high-performance donor polymer and a series of terpolymers with different content of TzBI are designed. As expected, the morphology is optimized gradually for improving charge generation and charge transport, also suppressing charge recombination. The champion device with 10% TzBI exhibited a power conversion efficiency (PCE) of 18.36%, which is 5% increase compared to the controlled device. This study presents a charming terpolymer strategy by highly electron-deficient and high dipole segment to realize a tradeoff between intermolecular interaction and backbone disorder, facilitating the optimization of morphology and elevation of fill factor and device efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Simple Model for the Benzene Hexafluorobenzene Interaction

While the experimental intermolecular distance distribution functions of pure benzene and pure hexafluorobenzene are well described by transferable all-atom force fields, the interaction between the two molecules (in a 1:1 mixture) is not well simulated. We demonstrate that the parameters of the transferable force fields are adequate to describe the intermolecular distance distribution if the charges are replaced by a set of charges that are not located at the atoms. Here, the simplest model that well describes the experimental distance distribution, between benzene and hexafluorobenzene, is that of a single ellipsoid for each molecule, representing the van der Waals interactions, and a set of three point charges (on the axis perpendicular to the arene plane) which give the same quadrupole moment as do the all atom charges from the transferable force fields.

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

Tension amplification in tethered layers of bottle-brush polymers

In this paper, molecular dynamics simulations of a coarse-grained bead–spring model have been used to study the effects of molecular crowding on the accumulation of tension in the backbone of bottle-brush polymers tethered to a flat substrate. The number of bottle-brushes per unit surface area, Σ, as well as the lengths of the bottle-brush backbones N bb (50 ≤ N bb ≤ 200) and side chains N sc (50 ≤ N sc ≤ 200) were varied to determine how the dimensions and degree of crowding of bottle-brushes give rise to bond tension amplification along the backbone, especially near the substrate. From these simulations, we have identified three separate regimes of tension. For low Σ, the tension is due solely to intramolecular interactions and is dominated by the side chain repulsion that governs the lateral brush dimensions. With increasing Σ, the interactions between bottle-brush polymers induce compression of the side chains, transmitting increasing tension to the backbone. For large Σ, intermolecular side chain repulsion increases, forcing side chain extension and reorientation in the direction normal to the surface and transmitting considerable tension to the backbone.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Statistical mechanics of transport processes in active fluids: Equations of hydrodynamics

We present the equations of hydrodynamics including mass, linear momentum, angular momentum, and energy are derived by coarse-graining the microscopic equations of motion for systems consisting of rotary dumbbells driven by internal torques. In deriving the balance of linear momentum, we find that the symmetry of the stress tensor is broken due to the presence of non-zero torques on individual particles. The broken symmetry of the stress tensor induces internal spin in the fluid and leads us to consider the balance of internal angular momentum in addition to the usual moment of momentum. In the absence of spin, the moment of momentum is the same as the total angular momentum. In deriving the form of the balance of total angular momentum, we find the microscopic expressions for the couple stress tensor that drives the spin field. We show that the couple stress contains contributions from both intermolecular interactions and the active forces. The presence of spin leads to the idea of balance of moment of inertia due to the constant exchange of particles in a small neighborhood around a macroscopic point. We derive the associated balance of moment of inertia at the macroscale and identify the moment of inertia flux that induces its transport. Lastly, we obtain the balances of total and internal energy of the active fluid and identify the sources of heat and heat fluxes in the system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2005

Electronic Transport through Self Assembled Thiol Molecules: Effect of Monolayer Order, Dynamics and Temperature

We present the charge transport and tunneling conductance of self assembled organic thiol molecules and discuss the influence of order and dynamics in the monolayer on the transport behavior and the effect of temperature. Conjugated thiol molecular wires and organometals such as terpyridine metal complexes provide a new platform for molecular electronic devices and we study their self assembly on Au(111) substrates by the scanning tunneling microscope. Determining the organization of the molecule and the ability to control the nature of its interface with the substrate is important for reliable performance of the molecular electronic devices. By concurrent scanning tunneling microscopy and spectroscopy studies on SAMs formed from oligo (phenelyne ethynelyne) monolayers with and without molecular order, we show that packing and order determine the response of a self assembled monolayer (SAM) to competing interactions. Molecular resolution STM imaging in vacuum shows that the OPES adopt an imcommensurate SAM structure on Au(111) with a rectangular unit cell. Tunneling spectroscopic measurements were performed on the SAM as a function of junction resistance. STS results show that the I-Vs are non linear and asymmetric due to the inherent asymmetry in the molecular structure, with larger currents at negative sample biases. The asymmetry increases with increasing junction resistance due to the asymmetry in the coupling to the leads. This is brought out clearly in the differential conductance, which also shows a gap at the Fermi level. We also studied the effect of order and dynamics in the monolayer on the charge transport and found that competing forces between the electric field, intermolecular interactions, tip-molecule physisorption and substrate-molecule chemisorption impact the transport measurements and its reliability and that the presence of molecular order is very important for reproducible transport measurements. Thus while developing new electronic platforms based on molecules, it is important to have a good control of the molecule-substrate interface, for the devices to perform reliably. While such a control would minimize fluctuations and dynamics in the ensemble, the real challenge is to develop device architectures that are tolerant to fluctuations, since they cannot be totally eliminated in these low dimensional soft systems. Results of temperature dependent STS measurements will also be discussed.

Dholakia, Geetha↗
DOE OSTI2016

The interplay of covalency, hydrogen bonding, and dispersion leads to a long range chiral network: The example of 2-butanol

The assembly of complex structures in nature is driven by an interplay between several intermolecular interactions, from strong covalent bonds to weaker dispersion forces. Understanding and ultimately controlling the self-assembly of materials requires extensive study of how these forces drive local nanoscale interactions and how larger structures evolve. Surface-based self-assembly is particularly amenable to modeling and measuring these interactions in well-defined systems. This study focuses on 2-butanol, the simplest aliphatic chiral alcohol. 2-butanol has recently been shown to have interesting properties as a chiral modifier of surface chemistry; however, its mode of action is not fully understood and a microscopic understanding of the role non-covalent interactions play in its adsorption and assembly on surfaces is lacking. In order to probe its surface properties, we employed high-resolution scanning tunneling microscopy and density functional theory (DFT) simulations. We found a surprisingly rich degree of enantiospecific adsorption, association, chiral cluster growth and ultimately long range, highly ordered chiral templating. Firstly, the chiral molecules acquire a second chiral center when adsorbed to the surface via dative bonding of one of the oxygen atom lone pairs. This interaction is controlled via the molecule’s intrinsic chiral center leading to monomers of like chirality, at both chiral centers, adsorbed on the surface. The monomers then associate into tetramers via a cyclical network of hydrogen bonds with an opposite chirality at the oxygen atom. The evolution of these square units is surprising given that the underlying surface has a hexagonal symmetry. Our DFT calculations, however, reveal that the tetramers are stable entities that are able to associate with each other by weaker van der Waals interactions and tessellate in an extended square network. This network of homochiral square pores grows to cover the whole Au(111) surface. Our data reveal that the chirality of a simple alcohol can be transferred to its surface binding geometry, drive the directionality of hydrogen-bonded networks and ultimately extended structure. Furthermore, this study provides the first microscopic insight into the surface properties of this important chiral modifier and provides a well-defined system for studying the network’s enantioselective interaction with other molecules.

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

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

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

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

Kinetic effects on double hysteresis in spin crossover molecular magnets analyzed with first order reversal curve diagram technique

In this paper, we analyze two types of hysteresis in spin crossover molecular magnets compounds in the framework of the First Order Reversal Curve (FORC) method. The switching between the two stable states in these compounds is accompanied by hysteresis phenomena if the intermolecular interactions are higher than a threshold. We have measured the static thermal hysteresis (TH) and the kinetic light induced thermal hysteresis (LITH) major loops and FORCs for the polycrystalline Fe(II) spin crossover compound [Fe{sub 1−x}Zn{sub x}(bbtr){sub 3}](ClO{sub 4}){sub 2} (bbtr = 1,4-di(1,2,3-triazol-1-yl)butane), either in a pure state (x = 0) or doped with Zn ions (x = 0.33) considering different sweeping rates. Here, we use this method not only to infer the domains distribution but also to disentangle between kinetic and static components of the LITH and to estimate the changes in the intermolecular interactions introduced by dopants. We also determined the qualitative relationship between FORC distributions measured for TH and LITH.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗