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 127 records · Page 7

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

A combined experimental and theoretical study of the supramolecular self-assembly of Cu(II) malonate complex assisted by various weak forces and water dimer

A Cu(II) malonate complex with formula [Cu(C{sub 3}H{sub 2}O{sub 4})(C{sub 6}H{sub 8}N{sub 2})(H{sub 2}O)]{sub 2}·4H{sub 2}O (1) [C{sub 6}H{sub 8}N{sub 2}=2-picolylamine, C{sub 3}H{sub 2}O{sub 4}{sup 2−}=malonate dianion] has been synthesized by mixing the reactants in their stoichiometric proportion and its crystal structure has been determined by single-crystal X-ray diffraction. In 1, monomeric neutral metal malonate units [Cu(C{sub 3}H{sub 2}O{sub 4})(C{sub 6}H{sub 8}N{sub 2})(H{sub 2}O)] are interlinked with each other through hydrogen bonds, weak lone pair⋯π and cuprophilic interactions to generate supramolecular dimers, which in turn further associated through hydrogen bonding to form infinite 1D chains. Water dimers, through series of hydrogen bonds and weak π–stacking forces are found to be responsible for interconnection of 1D chains, which resulted in a 3D network. A density functional (DFT) study of the energetic features of several noncovalent interactions observed in the solid state have been analyzed and characterized using Bader's theory of “atoms-in-molecules”. We also present here Hirshfeld surface analysis to investigate the close intermolecular contacts. - Graphical Abstract: Interplay of weak forces like hydrogen bonding, lone pair⋯π, Cu⋯Cu and π–stacking interactions leading to the formation of supramolecular network in [Cu(C{sub 3}H{sub 2}O{sub 4})(C{sub 6}H{sub 8}N{sub 2})(H{sub 2}O)]{sub 2}·4H{sub 2}O complex. - Highlights: • A complex of Cu(II) with malonate and 2-picolylamine is synthesized and X-ray characterized. • We report a density functional study of the energetic features of several noncovalent interactions • We perform Hirshfeld surface analysis to investigate the close intermolecular contacts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of intermolecular interactions on the infrared complex indices of refraction for binary liquid mixtures

This paper investigates the accuracy of deriving the composite optical constants of binary mixtures from only the complex indices of refraction of the neat materials. These optical constants enable the reflectance spectra of the binary mixtures to be modeled for multiple scenarios (e.g., different substrates, thicknesses, volume ratios), which is important for contact and standoff chemical detection. Using volume fractions, each mixture's complex index of refraction was approximated via three different mixing rules. To explore the impact of intermolecular interactions, these predictions are tested by experimental measurements for two representative sets of binary mixtures: (1) tributyl phosphate combined with n-dodecane, a non-polar medium, to represent mixtures which primarily interact via dispersion forces and (2) tributyl phosphate and 1-butanol to represent mixtures with polar functional groups that can also interact via dipole–dipole interactions, including hydrogen bonding. The residuals and the root-mean-square error between the experimental and calculated index values are computed and demonstrate that for miscible liquids in which the average geometry of the cross-interactions can be considered isotropic (e.g., dispersion), the refractive indices of the mixtures can be modeled using composite n and k values derived from volume fractions of the neat liquids. Conversely, in spectral regions where the geometry of the cross-interactions is more restricted and anisotropic (e.g., hydrogen bonding), the calculated n and k values vary from the measured values. The impact of these interactions on the reflectance spectra are then compared by modeling a thin film of the binary mixtures on an aluminum substrate using both the measured and the mathematically computed indices of refraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗