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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Theoretical study on electromagnetically induced transparency in molecular aggregate models using quantum Liouville equation method

Electromagnetically induced transparency (EIT), which is known as an efficient control method of optical absorption property, is investigated using the polarizability spectra and population dynamics obtained by solving the quantum Liouville equation. In order to clarify the intermolecular interaction effect on EIT, we examine several molecular aggregate models composed of three-state monomers with the dipole-dipole coupling. On the basis of the present results, we discuss the applicability of EIT in molecular aggregate systems to a new type of optical switch.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Molecular simulation of fluids with non-identical intermolecular potentials: Thermodynamic properties of 10-5 + 12-6 Mie potential binary mixtures

General methods for combining interactions between particles characterised by non-identical intermolecular potentials are investigated. The combination methods are tested by performing molecular dynamics simulations to determine the pressure, energy, isochoric and isobaric heat capacities, thermal expansion coefficient, isothermal compressibility, Joule-Thomson coefficient, and speed of sound of 10-5 + 12-6 Mie potential binary mixtures. In addition to the two non-identical Mie potentials, mixtures are also studied with non-identical intermolecular parameters. The combination methods are compared with results obtained by simply averaging the Mie exponents. When either the energy or size parameters are non-identical, very significant differences emerge in the thermodynamic properties predicted by the alternative combination methods. The isobaric heat capacity is the thermodynamic property that is most affected by the relative magnitude of the intermolecular potential parameters and the method for combining non-identical potentials. Either the arithmetic or geometric combination of potentials provides a simple and effective way of performing simulations involving mixtures of components characterised by non-identical intermolecular potentials, which is independent of their functional form.

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

The infrared spectrum of the He–C{sub 2}D{sub 2} complex

Spectra of the helium-acetylene complex are elusive because this weakly bound system lies close to the free rotor limit. Previously, limited assignments of He–C{sub 2}D{sub 2} transitions in the R(0) region of the ν{sub 3} fundamental band (≈2440 cm{sup −1}) were published. Here, new He–C{sub 2}D{sub 2} infrared spectra of this band are obtained using a tunable optical parametric oscillator laser source to probe a pulsed supersonic slit jet expansion from a cooled nozzle, and the analysis is extended to the weaker and more difficult P(1) and R(1) regions. A term value approach is used to obtain a consistent set of “experimental” energy levels. These are compared directly with calculations using two recently reported ab initio intermolecular potential energy surfaces, which exhibit small but significant differences. Rovibrational energies for the He–C{sub 2}H{sub 2} complex are also calculated using both surfaces. A Coriolis model, useful for predicting spectral intensities, is used to interpret the energy level patterns, and a comparison with the isoelectronic complex He–CO is made.

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

Ground state analytical ab initio intermolecular potential for the Cl{sub 2}-water system

The chlorine/water interface is of crucial importance in the context of atmospheric chemistry. Modeling the structure and dynamics at this interface requires an accurate description of the interaction potential energy surfaces. We propose here an analytical intermolecular potential that reproduces the interaction between the Cl{sub 2} molecule and a water molecule. Our functional form is fitted to a set of high level ab initio data using the coupled-cluster single double (triple)/aug-cc-p-VTZ level of electronic structure theory for the Cl{sub 2} − H{sub 2}O complex. The potential fitted to reproduce the three minima structures of 1:1 complex is validated by the comparison of ab initio results of Cl{sub 2} interacting with an increasing number of water molecules. Finally, the model potential is used to study the physisorption of Cl{sub 2} on a perfectly ordered hexagonal ice slab. The calculated adsorption energy, in the range 0.27 eV, shows a good agreement with previous experimental results.

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

Microwave spectroscopy of the seeded binary and ternary clusters CO-(pH{sub 2}){sub 2}, CO-pH{sub 2}-He, CO-HD, and CO-(oD{sub 2}){sub N=1,2}

We report the Fourier transform microwave spectra of the a-type J = 1-0 transitions of the binary and ternary CO-(pH{sub 2}){sub 2}, CO-pH{sub 2}-He, CO-HD, and CO-(oD{sub 2}){sub N=1,2} clusters. In addition to the normal isotopologue of CO for all clusters, we observed the transitions of the minor isotopologues, {sup 13}C{sup 16}O, {sup 12}C{sup 18}O, and {sup 13}C{sup 18}O, for CO-(pH{sub 2}){sub 2} and CO-pH{sub 2}-He. All transitions lie within 335 MHz of the experimentally or theoretically predicted values. In comparison to previously reported infrared spectra [Moroni et al., J. Chem. Phys. 122, 094314 (2005)], we are able to tentatively determine the vibrational shift for CO-pH{sub 2}-He, in addition to its b-type J = 1-0 transition frequency. The a-type frequency of CO-pH{sub 2}-He is similar to that of CO-He{sub 2} [Surin et al., Phys. Rev. Lett. 101, 233401 (2008)], suggesting that the pH{sub 2} molecule has a strong localizing effect on the He density. Perturbation theory analysis of CO-oD{sub 2} reveals that it is approximately T-shaped, with an anisotropy of the intermolecular potential amounting to ∼9 cm{sup −1}.

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

Spin-unrestricted random-phase approximation with range separation: Benchmark on atomization energies and reaction barrier heights

We consider several spin-unrestricted random-phase approximation (RPA) variants for calculating correlation energies, with and without range separation, and test them on datasets of atomization energies and reaction barrier heights. We show that range separation greatly improves the accuracy of all RPA variants for these properties. Moreover, we show that a RPA variant with exchange, hereafter referred to as RPAx-SO2, first proposed by Szabo and Ostlund [J. Chem. Phys. 67, 4351 (1977)] in a spin-restricted closed-shell formalism, and extended here to a spin-unrestricted formalism, provides on average the most accurate range-separated RPA variant for atomization energies and reaction barrier heights. Since this range-separated RPAx-SO2 method had already been shown to be among the most accurate range-separated RPA variants for weak intermolecular interactions [J. Toulouse et al., J. Chem. Phys. 135, 084119 (2011)], this works confirms range-separated RPAx-SO2 as a promising method for general chemical applications.

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

Similarity-transformed perturbation theory on top of truncated local coupled cluster solutions: Theory and applications to intermolecular interactions

Your correspondents develop and apply fully nonorthogonal, local-reference perturbation theories describing non-covalent interactions. Our formulations are based on a Löwdin partitioning of the similarity-transformed Hamiltonian into a zeroth-order intramonomer piece (taking local CCSD solutions as its zeroth-order eigenfunction) plus a first-order piece coupling the fragments. If considerations are limited to a single molecule, the proposed intermolecular similarity-transformed perturbation theory represents a frozen-orbital variant of the “(2)”-type theories shown to be competitive with CCSD(T) and of similar cost if all terms are retained. Different restrictions on the zeroth- and first-order amplitudes are explored in the context of large-computation tractability and elucidation of non-local effects in the space of singles and doubles. To accurately approximate CCSD intermolecular interaction energies, a quadratically growing number of variables must be included at zeroth-order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Reliable prediction of three-body intermolecular interactions using dispersion-corrected second-order Møller-Plesset perturbation theory

Three-body and higher intermolecular interactions can play an important role in molecular condensed phases. Recent benchmark calculations found problematic behavior for many widely used density functional approximations in treating 3-body intermolecular interactions. Here, we demonstrate that the combination of second-order Møller-Plesset (MP2) perturbation theory plus short-range damped Axilrod-Teller-Muto (ATM) dispersion accurately describes 3-body interactions with reasonable computational cost. The empirical damping function used in the ATM dispersion term compensates both for the absence of higher-order dispersion contributions beyond the triple-dipole ATM term and non-additive short-range exchange terms which arise in third-order perturbation theory and beyond. Empirical damping enables this simple model to out-perform a non-expanded coupled Kohn-Sham dispersion correction for 3-body intermolecular dispersion. The MP2 plus ATM dispersion model approaches the accuracy of O(N{sup 6}) methods like MP2.5 or even spin-component-scaled coupled cluster models for 3-body intermolecular interactions with only O(N{sup 5}) computational cost.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Pressure-induced structural changes in NH{sub 4}Br

We report angle dispersive X-ray diffraction (XRD) measurements and Raman spectroscopy on NH{sub 4}Br up to 70.0 GPa at room temperature. Three thermodynamically stable phases (phases II, IV, and V) are confirmed and a new possible phase (phase VI) of P2{sub 1}/m symmetry is proposed whose structure was established from Rietveld refinement of synchrotron XRD data for the first time. The phase sequence observed in NH{sub 4}Br is in accordance with phase II → IV → V → VI. Phase V transforms into phase VI at about 57.8 GPa with a huge volume reduction of 30%. Still, the intramolecular distances are analyzed to better understand the nature of structures. The H–H interactions become markedly more important as the N–Br distances are compacted, which is probably the reason of the kink of symmetric stretching band (ν{sub 1}) at the transition pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Graphene as a hexagonal 2-lattice: Evaluation of the in-plane material constants for the linear theory. A multiscale approach

Continuum modeling of free-standing graphene monolayer, viewed as a two dimensional 2-lattice, requires specification of the components of the shift vector that acts as an auxiliary variable. If only in-plane motions are considered, the energy depends on an in-plane strain measure and the shift vector. The assumption of geometrical and material linearity leads to quadratic energy terms with respect to the shift vector, the strain tensor, and their combinations. Graphene's hexagonal symmetry reduces the number of independent moduli then to four. We evaluate these four material parameters using molecular calculations and the adaptive intermolecular reactive empirical bond order potential and compare them with standard linear elastic constitutive modeling. The results of our calculations show that the predicted values are in reasonable agreement with those obtained solely from our molecular calculations as well as those from the literature. To the best of our knowledge, this is the first attempt to measure mechanical properties when graphene is modeled as a hexagonal 2-lattice. This work targets at the continuum scale when the insight measurements come from finer scales using atomistic simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Water droplet excess free energy determined by cluster mitosis using guided molecular dynamics

Atmospheric aerosols play a vital role in affecting climate by influencing the properties and lifetimes of clouds and precipitation. Understanding the underlying microscopic mechanisms involved in the nucleation of aerosol droplets from the vapour phase is therefore of great interest. One key thermodynamic quantity in nucleation is the excess free energy of cluster formation relative to that of the saturated vapour. In our current study, the excess free energy is extracted for clusters of pure water modelled with the TIP4P/2005 intermolecular potential using a method based on nonequilibrium molecular dynamics and the Jarzynski relation. The change in free energy associated with the “mitosis” or division of a cluster of N water molecules into two N/2 sub-clusters is evaluated. This methodology is an extension of the disassembly procedure used recently to calculate the excess free energy of argon clusters [H. Y. Tang and I. J. Ford, Phys. Rev. E 91, 023308 (2015)]. Our findings are compared to the corresponding excess free energies obtained from classical nucleation theory (CNT) as well as internally consistent classical theory (ICCT). The values of the excess free energy that we obtain with the mitosis method are consistent with CNT for large cluster sizes but for the smallest clusters, the results tend towards ICCT; for intermediate sized clusters, we obtain values between the ICCT and CNT predictions. Furthermore, the curvature-dependent surface tension which can be obtained by regarding the clusters as spherical droplets of bulk density is found to be a monotonically increasing function of cluster size for the studied range. The data are compared to other values reported in the literature, agreeing qualitatively with some but disagreeing with the values determined by Joswiak et al. [J. Phys. Chem. Lett. 4, 4267 (2013)] using a biased mitosis approach; an assessment of the differences is the main motivation for our current study.

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

Simulations of the THz spectrum of liquid water incorporating the effects of intermolecular charge fluxes through hydrogen bonds

The intensity of the band at ∼200 cm{sup −1} (∼6 THz) in the Terahertz spectrum of liquid water mainly arises from the modulations of the extent of intermolecular charge transfer through hydrogen bonds, called intermolecular charge fluxes, occurring upon molecular translations along the O…H hydrogen bonds. To obtain reasonable spectral profiles from simulations, it is necessary to effectively incorporate the effects of those intermolecular charge fluxes, but apparently it is not possible by simple classical molecular dynamics simulations with fixed atomic partial charges even if they are amended by molecular induced dipoles due to intermolecular electrostatic interactions. The present paper shows how we can do reasonably correct spectral simulations, without resorting to fully ab initio molecular dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

The origins of intra- and inter-molecular vibrational couplings: A case study of H{sub 2}O-Ar on full and reduced-dimensional potential energy surface

The origin and strength of intra- and inter-molecular vibrational coupling is difficult to probe by direct experimental observations. However, explicitly including or not including some specific intramolecular vibrational modes to study intermolecular interaction provides a precise theoretical way to examine the effects of anharmonic coupling between modes. In this work, a full-dimension intra- and inter-molecular ab initio potential energy surface (PES) for H{sub 2}O–Ar, which explicitly incorporates interdependence on the intramolecular (Q{sub 1}, Q{sub 2}, Q{sub 3}) normal-mode coordinates of the H{sub 2}O monomer, has been calculated. In addition, four analytic vibrational-quantum-state-specific PESs are obtained by least-squares fitting vibrationally averaged interaction energies for the (v{sub 1}, v{sub 2}, v{sub 3}) = (0, 0, 0), (0, 0, 1), (1, 0, 0), (0, 1, 0) states of H{sub 2}O to the three-dimensional Morse/long-range potential function. Each vibrationally averaged PES fitted to 442 points has root-mean-square (rms) deviation smaller than 0.15 cm{sup −1}, and required only 58 parameters. With the 3D PESs of H{sub 2}O–Ar dimer system, we employed the combined radial discrete variable representation/angular finite basis representation method and Lanczos algorithm to calculate rovibrational energy levels. This showed that the resulting vibrationally averaged PESs provide good representations of the experimental infrared data, with rms discrepancies smaller than 0.02 cm{sup −1} for all three rotational branches of the asymmetric stretch fundamental transitions. The infrared band origin shifts associated with three fundamental bands of H{sub 2}O in H{sub 2}O–Ar complex are predicted for the first time and are found to be in good agreement with the (extrapolated) experimental values. Upon introduction of additional intramolecular degrees of freedom into the intermolecular potential energy surface, there is clear spectroscopic evidence of intra- and intermolecular vibrational couplings.

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

Communication: Accurate higher-order van der Waals coefficients between molecules from a model dynamic multipole polarizability

Due to the absence of the long-range van der Waals (vdW) interaction, conventional density functional theory (DFT) often fails in the description of molecular complexes and solids. In recent years, considerable progress has been made in the development of the vdW correction. However, the vdW correction based on the leading-order coefficient C{sub 6} alone can only achieve limited accuracy, while accurate modeling of higher-order coefficients remains a formidable task, due to the strong non-additivity effect. Here, we apply a model dynamic multipole polarizability within a modified single-frequency approximation to calculate C{sub 8} and C{sub 10} between small molecules. We find that the higher-order vdW coefficients from this model can achieve remarkable accuracy, with mean absolute relative deviations of 5% for C{sub 8} and 7% for C{sub 10}. Inclusion of accurate higher-order contributions in the vdW correction will effectively enhance the predictive power of DFT in condensed matter physics and quantum chemistry.

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

Carbon dioxide in an ionic liquid: Structural and rotational dynamics

Ionic liquids (ILs), which have widely tunable structural motifs and intermolecular interactions with solutes, have been proposed as possible carbon capture media. In order to inform the choice of an optimal ionic liquid system, it can be useful to understand the details of dynamics and interactions on fundamental time scales (femtoseconds to picoseconds) of dissolved gases, particularly carbon dioxide (CO 2 ), within the complex solvation structures present in these uniquely organized materials. The rotational and local structural fluctuation dynamics of CO 2 in the room temperature ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EmimNTf2) were investigated by using ultrafast infrared spectroscopy to interrogate the CO 2 asymmetric stretch. Polarization-selective pump probe measurements yielded the orientational correlation function of the CO 2 vibrational transition dipole. It was found that reorientation of the carbon dioxide occurs on 3 time scales: 0.91 ± 0.03, 8.3 ± 0.1, 54 ± 1 ps. The initial two are attributed to restricted wobbling motions originating from a gating of CO 2 motions by the IL cations and anions. The final (slowest) decay corresponds to complete orientational randomization. Two-dimensional infrared vibrational echo (2D IR) spectroscopy provided information on structural rearrangements, which cause spectral diffusion, through the time dependence of the 2D line shape. Analysis of the time-dependent 2D IR spectra yields the frequency-frequency correlation function (FFCF). Polarization-selective 2D IR experiments conducted on the CO 2 asymmetric stretch in the parallel- and perpendicular-pumped geometries yield significantly different FFCFs due to a phenomenon known as reorientation-induced spectral diffusion (RISD), revealing strong vector interactions with the liquid structures that evolve slowly on the (independently measured) rotation time scales. To separate the RISD contribution to the FFCF from the structural spectral diffusion contribution, the previously developed first order Stark effect RISD model is reformulated to describe the second order (quadratic) Stark effect—the first order Stark effect vanishes because CO 2 does not have a permanent dipole moment. Through this analysis, we characterize the structural fluctuations of CO 2 in the ionic liquid solvation environment, which separate into magnitude-only and combined magnitude and directional correlations of the liquid’s time dependent electric field. This new methodology will enable highly incisive comparisons between CO 2 dynamics in a variety of ionic liquid systems.

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

Study of interaction in silica glass via model potential approach

Silica is one of the most commonly encountered substances in daily life and in electronics industry. Crystalline SiO{sub 2} (in several forms: quartz, cristobalite, tridymite) is an important constituent of many minerals and gemstones, both in pure form and mixed with related oxides. Cohesive energy of amorphous SiO{sub 2} has been investigated via intermolecular potentials i.e weak Van der Waals interaction and Morse type short-range interaction. We suggest a simple atom-atom based Van der Waals as well as Morse potential to find cohesive energy of glass. It has been found that the study of silica structure using two different model potentials is significantly different. Van der Waals potential is too weak (P.E =0.142eV/molecule) to describe the interaction between silica molecules. Morse potential is a strong potential, earlier given for intramolecular bonding, but if applied for intermolecular bonding, it gives a value of P.E (=−21.92eV/molecule) to appropriately describe the structure of silica.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Communication: Charge-population based dispersion interactions for molecules and materials

We introduce a system-independent method to derive effective atomic C{sub 6} coefficients and polarizabilities in molecules and materials purely from charge population analysis. This enables the use of dispersion-correction schemes in electronic structure calculations without recourse to electron-density partitioning schemes and expands their applicability to semi-empirical methods and tight-binding Hamiltonians. We show that the accuracy of our method is en par with established electron-density partitioning based approaches in describing intermolecular C{sub 6} coefficients as well as dispersion energies of weakly bound molecular dimers, organic crystals, and supramolecular complexes. We showcase the utility of our approach by incorporation of the recently developed many-body dispersion method [Tkatchenko et al., Phys. Rev. Lett. 108, 236402 (2012)] into the semi-empirical density functional tight-binding method and propose the latter as a viable technique to study hybrid organic-inorganic interfaces.

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

Vibrational spectroscopic, structural and nonlinear optical activity studies on 6-aminonicotinamide: A DFT approach

The conformational analysis was carried out for 6-aminonicotinamide (ANA) using potential energy surface scan method and the most stable optimized conformer was predicted. The theoretical vibrational frequencies were calculated for the optimized geometry using DFT/B3LYP cc-pVQZ basis set by Gaussian 09 Program. The vibrational frequencies were assigned on the basis of potential energy distribution calculation using VEDA 4.0 program. The Mulliken atomic charge values were calculated. In the Frontier molecular orbitals analysis, the molecular reactivity, kinetic stability, intermolecular charge transfer studies and the related molecular properties were calculated. The ultraviolet-visible spectrum was simulated for both in the gas phase and liquid phase (ethanol) and the π to π* electronic transition was predicted. The nonlinear optical (NLO) activity was studied by means of the first order hyperpolarizability value, which was 8.61 times greater than the urea and the natural bond orbital analysis was also performed to confirm the NLO activity of the molecule. Hence, the ANA molecule is a promising candidate for the NLO materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗