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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2016

Transfer of energy between a pair of molecules near a plasmonic core-shell nanoparticle: Tunability and sensing

Our model is applied to the calculation of interaction energy between a pair of dipolar molecules (point dipoles) in the vicinity of a nanoshell monomer with core-shell structure, based on the dipole quasi-electrostatic theory of classical electrodynamics and using the Drude and Maxwell-Garnett model. In other words, this work discusses the intermolecular energy transfer from a donor molecule to an acceptor molecule near a spherical nanoparticle that is important for practical applications like sensing. It is shown that the proximity of plasmonic nanoparticles can have a strong effect on the energy transfer between molecules. In addition to the influence of the size, composition, embedding medium, and the filling fraction of doped particles on the interaction energy, the contribution of the dipolar, quadrupolar, octupolar, hexadecapolar, triakontadipolar, and higher order multipole interactions is presented and analyzed. Briefly, we will show that it is possible to achieve enhanced energy transfer by manipulation of different parameters as mentioned above.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Blind test of density-functional-based methods on intermolecular interaction energies

In the past decade, a number of approaches have been developed to fix the failure of (semi)local density-functional theory (DFT) in describing intermolecular interactions. In this work, the performance of several such approaches with respect to highly accurate benchmarks is compared on a set of separation-dependent interaction energies for ten dimers. Since the benchmarks were unknown before the DFT-based results were collected, this comparison constitutes a blind test of these methods.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Diffusion quantum Monte Carlo calculations of SrFeO 3 and LaFeO 3

The equations of state, formation energy, and migration energy barrier of the oxygen vacancy in SrFeO 3 and LaFeO 3 were calculated in this paper with the diffusion quantum Monte Carlo (DMC) method. Calculations were also performed with various Density Functional Theory (DFT) approximations for comparison. DMC reproduces the measured cohesive energies of these materials with errors below 0.23(5) eV and the structural properties within 1% of the experimental values. The DMC formation energies of the oxygen vacancy in SrFeO 3 and LaFeO 3 under oxygen-rich conditions are 1.3(1) and 6.24(7) eV, respectively. Similar calculations with semi-local DFT approximations for LaFeO 3 yielded vacancy formation energies 1.5 eV lower. Comparison of charge density evaluated with DMC and DFT approximations shows that DFT tends to overdelocalize the electrons in defected SrFeO 3 and LaFeO 3 . Finally, calculations with DMC and local density approximation yield similar vacancy migration energy barriers, indicating that steric/electrostatic effects mainly determine migration barriers in these materials.

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

Polymorphic improvement of Stillinger-Weber potential for InGaN

A Stillinger-Weber potential is computationally very efficient for molecular dynamics simulations. Despite its simple mathematical form, the Stillinger-Weber potential can be easily parameterized to ensure that crystal structures with tetrahedral bond angles (e.g., diamond-cubic, zinc-blende, and wurtzite) are stable and have the lowest energy. As a result, the Stillinger-Weber potential has been widely used to study a variety of semiconductor elements and alloys. When studying an A-B binary system, however, the Stillinger-Weber potential is associated with two major drawbacks. First, it significantly overestimates the elastic constants of elements A and B, limiting its use for systems involving both compounds and elements (e.g., an A/AB multilayer). Second, it prescribes equal energy for zinc-blende and wurtzite crystals, limiting its use for compounds with large stacking fault energies. Here in this paper, we utilize the polymorphic potential style recently implemented in LAMMPS to develop a modified Stillinger-Weber potential for InGaN that overcomes these two problems.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Energy decomposition analysis for exciplexes using absolutely localized molecular orbitals

An energy decomposition analysis (EDA) scheme is developed for understanding the intermolecular interaction involving molecules in their excited states. The EDA utilizes absolutely localized molecular orbitals to define intermediate states and is compatible with excited state methods based on linear response theory such as configuration interaction singles and time-dependent density functional theory. The shift in excitation energy when an excited molecule interacts with the environment is decomposed into frozen, polarization, and charge transfer contributions, and the frozen term can be further separated into Pauli repulsion and electrostatics. These terms can then be added to their counterparts obtained from the ground state EDA to form a decomposition of the total interaction energy. The EDA scheme is applied to study a variety of systems, including some model systems to demonstrate the correct behavior of all the proposed energy components as well as more realistic systems such as hydrogen-bonding complexes (e.g., formamide-water, pyridine/pyrimidine-water) and halide (F - ,Cl - )-water clusters that involve charge-transfer-to-solvent excitations.

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

Fast dynamics of a hydrogen-bonding glass forming liquid: Chemical exchange-induced spectral diffusion in 2D IR spectroscopy

Polarization-selective Two Dimensional Infrared (2D IR) and IR pump-probe spectroscopies have been performed on the hydrogen bonding glass forming liquid 2-biphenylmethanol doped with the long-lived vibrational probe phenylselenocyanate over a wide range of temperatures. The spectral diffusion seen in the 2D spectra was found to have a large polarization dependence, in large excess of what is predicted by standard theory. This anomaly was explained by decomposing the 2D spectra into hydrogen-bonding and non-bonding components, which exchange through large-angle orientational motion. By adapting chemical exchange theories, parameters for the component peaks were then calculated by fitting the polarization-dependent spectral diffusion and the pump-probe anisotropy. A model of highly heterogeneous exchange and orientational dynamics was used to explain the observed time dependences as a function of temperature on fast time scales. The experimental observations, the kinetic modeling, and physical arguments lead to the determination of the times for interconversion of slow dynamics structural domains to fast dynamics structural domains in the supercooled liquid as a function of temperature. The slow to fast domain interconversion times range from 40 ps at 355 K to 5000 ps at 270 K.

14 SOLAR ENERGY↗
DOE OSTI2019

One-directional thermal transport in densely aligned single-wall carbon nanotube films

Individual carbon nanotubes (CNTs) possess extremely high thermal conductivities. However, the thermal conductivities and their anisotropy of macroscopic assemblies of CNTs have so far remained small. Here, we report the results of directional thermal transport measurements on a nearly perfectly aligned CNT film fabricated via controlled vacuum filtration. We found the thermal conductivity to be 43 ± 2.2 W m -1 K -1 with a record-high thermal anisotropy of 500. From the temperature dependence of the thermal conductivity and its agreement with the atomistic phonon transport calculation, we conclude that the effect of intertube thermal resistance on heat conduction in the alignment direction is negligible because of the large contact area between CNTs. These observations thus represent ideal unidirectional thermal transport, i.e., the thermal conductivity of the film is determined solely by that of individual CNTs.

42 ENGINEERING↗
DOE OSTI2020

Examining the order-of-limits problem and lattice constant performance of the Tao–Mo functional

In their recent communication a semi-local density functional derived from the density matrix expansion of the exchange hole localised by a general coordinate transformation. We show that the order-of-limits problem present in the functional, dismissed as harmless in the original publication, causes severe errors in predicted phase transition pressures. We also show that the claim that lattice volume prediction accuracy exceeds that of existing similar functionals was based on comparison to reference data that misses anharmonic zero-point expansion and consequently overestimates accuracy. Here, by highlighting these omissions, we give a more accurate assessment of the Tao-Mo functional and show a possible direction for resolving the order-of-limits problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

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

Ligand structure and adsorption free energy of nanocrystals on solid substrates

In this paper, we present an investigation on the absorption of alkylthiolated nanocrystals on a solid substrate. We calculate adsorption free energies and report a number of effects induced by the substrate. Nearest neighbor distances and bonding free energies are significantly different than for a free floating case, there is a weakening of bonding free energies among nanocrystals, and the adsorption is manifestly anisotropic, i.e., stronger along certain directions of the nanocrystal core. We contend that this last result accounts for the Bain transition (fcc → bcc) observed in experimental results. We report the presence of vortices induced by the substrate, which explain the increased nearest neighbor distance among nanocrystals, which is in excellent quantitative agreement with experimental results and with the predictions of the Orbifold Topological Model. Implications for the assembly of nanostructures and future experiments are also discussed.

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

Phonon gas model for thermal conductivity of dense, strongly interacting liquids

Developing predictive thermal property models for liquids based on microscopic principles has been elusive. The difficulty is that liquids have gas-like and solid-like attributes that are at odds when considering the frameworks of microscopic models: Models for gases are simple due to randomness and low density, whereas models for crystalline solids rely on symmetry and long-range order for easier calculation. The short-range order in liquids does, however, provide structure to neighboring molecules similar to amorphous solids, and there have been recent advances indicating that collective vibrational modes store heat in liquids. Models combining Debye approximations from solid-state physics and Frenkel’s theory of liquids can accurately predict the heat capacity of liquids. Phonon-like dispersions in liquids have also been widely observed in neutron scattering experiments. These developments motivate us to propose a model where high-frequency vibrational modes, which travel at the speed of sound and have a mean free path on the order of the average intermolecular distance, conduct heat in liquids. We use this liquid phonon gas model to calculate the thermal conductivity of liquids with varying intermolecular interaction energies from strongest to weakest—Coulomb, hydrogen-bonding, Keesom, and London dispersion energy. Generally, the model is more accurate as the intermolecular interaction energy and density of liquids increase. The calculated thermal conductivity of Coulombic-bound molten sodium nitrate and hydrogen-bonded water is within 1.46% and 2.98% of the experimentally measured values, respectively, across their entire temperature ranges. Further modal analysis of the velocity and the mean free path of collective vibrations could establish the liquid phonon gas model as an accurate model for weakly interacting liquids as well.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water

Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.

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

Effect of chemical substituents attached to the zwitterion cation on dielectric constant

Materials with high dielectric constant, ε s , are desirable in a wide range of applications including energy storage and actuators. Recently, zwitterionic liquids have been reported to have the largest ε s of any liquid and, thus, have the potential to replace inorganic fillers to modulate the material ε s . Although the large ε s for zwitterionic liquids is attributed to their large molecular dipole, the role of chemical substituents attached to the zwitterion cation on ε s is not fully understood, which is necessary to enhance the performance of soft energy materials. Here, we report the impact of zwitterionic liquid cation chemical substituents on ε s (50 < ε s < 300 at room temperature). Dielectric relaxation spectroscopy reveals that molecular reorientation is the main contributor to the high ε s . The low Kirkwood factor g calculated for zwitterionic liquids (e.g., 0.1–0.2) suggests the tendency for the antiparallel zwitterion dipole alignment expected from the strong electrostatic intermolecular interactions. With octyl cation substituents, the g is decreased due to the formation of hydrophobic-rich domains that restrict molecular reorientation under applied electric fields. In contrast, when zwitterion cations are functionalized with ethylene oxide (EO) segments, g increases due to the EO segments interacting with the cations, allowing more zwitterion rotation in response to the applied field. The reported results suggest that high ε s zwitterionic liquids require a large molecular dipole, compositionally homogeneous liquids (e.g., no aggregation), a maximized zwitterion number density, and a high g, which is achievable by incorporating polar chemical substituents onto the zwitterion cations.

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

Rotational and translational diffusion of liquid n-hexane: EFP-based molecular dynamics analysis

Molecular Dynamics (MD) simulations based on the Effective Fragment Potential (EFP) method are utilized to provide a comprehensive assessment of diffusion in liquid n-hexane. We decompose translational diffusion into components along and orthogonal to the long axis of the molecule. Rotational diffusion is decomposed into tumbling and spinning motions about this axis. Our analysis yields four corresponding diffusion coefficients which are related to diagonal entries in the complete 6 × 6 diffusion tensor accounting for the three rotational and three translational degrees of freedom and for the potential coupling between them. However, coupling between different degrees of freedom is expected to be minimal for a natural choice of the molecular body-fixed axis, so then off-diagonal entries in the tensor are negligible. This expectation is supported by a hydrodynamic analysis of the diffusion tensor which treats the liquid surrounding the molecule being tracked as a viscous continuum. Thus, the EFP MD analysis provides a comprehensive characterization of diffusion and also reveals expected shortcomings of the hydrodynamic treatment, particularly for rotational diffusion, when applied to neat liquids.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi’s golden rule rate

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [Li et al., J. Chem. Phys. 154, 094124 (2021)], we systematically study the mechanism of polariton relaxation for liquid CO 2 under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations confirm that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi’s golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. Finally, we also summarize recent progress on probing nonreactive VSC dynamics with CavMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗