Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 OSTI2020

Protein Molecular Dynamics Simulations with Approximate QM: What Can We Learn?

Classical force fields are essential for computer simulations of proteins and are typically parameterized to reproduce secondary and tertiary structure of isolated proteins. However, while protein-protein interactions are ubiquitous in nature, they are not considered in parameterization efforts and are far less understood than isolated proteins. A better characterization of intermolecular interactions is widely recognized as a key to revolutionizing drug and therapeutic developments with high-throughput computational screening. Urgently needed is a critical assessment of the performance of modern protein force fields against first-principles electronic structure methods and experiments. In a daring step toward this goal, we here describe a comparison of peptide folding dynamics as predicted by a molecular mechanics force field on the one hand and by an approximate electronic structure quantum mechanical (QM) method based on density-functional tight-binding (DFTB) on the other. We further compare the dynamics from straightforward DFTB simulations with a near-linear scaling version of DFTB for massively parallel computation based on the fragment molecular orbital (FMO-DFTB) method. We illustrate differences between the phenomenology of the folding dynamics from these three methods for a small model peptide, as well as charge polarization and dynamic fluctuations, point out possible correlations and implications for force field developers, and discuss the lessons learned that might become applicable to future predictive high-throughput computer screening for personalized neoantigen cancer therapy.

Irle, Stephan↗
DOE OSTI2024

Multiscale investigation of the mechanism of biomass deconstruction in the dimethyl isosorbide/water Co-solvent pretreatment system

In the context of promoting a circular bioeconomy, the development of green and efficient lignocellulosic biomass pretreatment technologies so as to realize high value-added biomass utilization is of intense interest. We demonstrated the potential of the bio-based green solvent dimethyl isosorbide (DMI) for the fractionation of Eucalyptus biomass with excellent performance. Here, to investigate the mechanisms involved in biomass fractionation, microimaging and microspectroscopic techniques were employed together with molecular dynamics (MD) simulation and COSMO-RS quantum chemical calculations to derive multiscale information. Both the microstructure and regional chemistry of the cell wall vary significantly with the volume ratio of DMI/H 2 O. The strongest effects were found at DMI/H 2 O = 9:1 and showed visible cell wall tearing cracks and cell wall deformation and collapse as well as the lowest values of cell wall thickness and circularity. From the MD simulations, lignin exhibits collapsed-like structure in pure H 2 O with low solvent accessibility surface area (SASA) and radius of gyration (R g ). In contrast, lignin in DMI/H 2 O shows extended structure with high SASA and solvent interactions dominated by van der Waals forces, with maximal contact in the 9:1 (v/v) system. Further, the COSMO-RS calculated sigma (σ-) potential suggests the intermolecular interactions in DMI and DMI/H 2 O co-solvent are weak, leading to stronger interaction with lignin and correspondingly higher lignin dissolution. The radial distribution functions and σ-potential all show that again DMI/H 2 O at 9:1 is an optimal volume ratio for high lignin dissolution. Finally, this study provides a solvent-ratio dependent mechanism for the action of polar aprotic solvents in the deconstruction of biomass.

09 BIOMASS FUELS↗
DOE OSTI2021

Predicting and Understanding Non-Covalent Interactions Using Novel Forms of Symmetry-Adapted Perturbation Theory

Although sometimes derided as “weak” interactions, non-covalent forces play a critical role in ligand binding and crystal packing and in determining the conformational landscape of flexible molecules. Symmetry-adapted perturbation theory (SAPT) provides a framework for accurate ab initio calculation of intermolecular interactions and furnishes a natural decomposition of the interaction energy into physically meaningful components: semiclassical electrostatics (rigorously obtained from monomer charge densities), Pauli or steric repulsion, induction (including both polarization and charge transfer), and dispersion. This decomposition helps to foster deeper understanding of non-covalent interactions and can be used to construct transferable, physics-based force fields. Separability of the SAPT interaction energy also provides the flexibility to construct composite methods, a feature that we exploit to improve the description of dispersion interactions. Furthermore, these are challenging to describe accurately because they arise from nonlocal electron correlation effects that appear for the first time at second order in perturbation theory but are not quantitatively described at that level.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of Long-Range Forces on the Transition States and Dynamics of NaCl Ion-Pair Dissociation in Water

We study NaCl ion-pair dissociation in a dilute aqueous solution using computer simulations both for the full system with long-range Coulomb interactions and for a well-chosen reference system with short-range intermolecular interactions. Analyzing results using concepts from Local Molecular Field (LMF) theory and the recently proposed AI-based analysis tool “State predictive information bottleneck” (SPIB), we show that the system with short-range interactions can accurately reproduce the transition rate for the dissociation process, the dynamics for moving between the underlying metastable states, and the transition state ensemble. Contributions from long-range interactions can be largely neglected for these processes because long-range forces from the direct interionic Coulomb interactions are almost completely canceled (>90%) by those from solvent interactions over the length scale where the transition takes place. Thus, for this important monovalent ion-pair system, short-range forces alone are able to capture detailed consequences of the collective solvent motion, allowing the use of physically suggestive and computationally efficient short-range models for the dissociation event. Here, we believe that the framework here should be applicable to disentangling mechanisms for more complex processes such as multivalent ion disassociation, where previous work has suggested that long-range contributions may be more important.

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

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deconstructing the Local Intermolecular Ordering and Dynamics of Liquid Chloroform and Bromoform

Local intermolecular structure and dynamics of the polar molecular liquids chloroform and bromoform are studied by molecular dynamics simulation. Structural distribution functions, including 1- and 2-D pair correlations and dipole contour plots allow direct comparison and show agreement with recent analyses of diffraction experiments. Studies of the haloforms’ reorientational dynamics and longevity of structural features resulting from intermolecular interaction extend previous work toward deeper understanding of the factors controlling these features. Analyses of ensemble average structures and dynamical properties isolate mass, electrostatics, and steric packing as driving forces or contributing factors for the observed ordering and dynamics.

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

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Science & Technology - Other Topics↗
DOE OSTI2025

Structure and Dynamics of CO 2 at the Air–Water Interface from Classical and Neural Network Potentials

The accurate description of the structure and dynamics of CO 2 at the instantaneous air–water interface, along with the effects of surface fluctuations on the CO 2 -transport processes, is essential for the development of negative emission technologies aimed at minimizing climate change. In this study, we performed molecular dynamics simulations of CO 2 at the air–water interface using neural network potentials (NNPs) trained on ab initio data generated through density-functional-theory-based molecular dynamics simulations. We compared these results with classical force fields to assess their performance in modeling interfacial CO 2 behavior. Our findings revealed that the asymmetric interactions, coupled with thermal surface fluctuations at the air–water interface, significantly influence CO 2 transport into the aqueous phase. The simulations demonstrate that classical force fields underestimate both the free energy of CO 2 transport and the strength of its interactions at the interface compared with the neural network potentials. In conclusion, the free energy and the interfacial dynamics of CO 2 are primarily influenced by the distribution of water within the instantaneous interfacial water layer, responsible for creating an asymmetric intermolecular interaction environment within the interfacial region.

Ab initio molecular dynamics↗
DOE OSTI2022

A classical model for three-body interactions in aqueous ionic systems

We present a classical induction model to evaluate the three-body ion–water–water (I–W–W) and water–water–water (W–W–W) interactions in aqueous ionic systems. The classical description of the induction energy is based on electrostatic distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole–quadrupole on the O and H atoms of water. The monatomic ions were described by a point charge and a dipole–dipole polarizability, while for the polyatomic ions, distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole–quadrupole were used. The accuracy of the classical model is benchmarked against an accurate dataset of 936 (I–W–W) and 2184 (W–W–W) three-body terms for 13 different monatomic and polyatomic cation and anion systems. The classical model shows excellent agreement with the reference second order Moller–Plesset and coupled-cluster single double and perturbative triple [CCSD(T)] three-body energies. The Root-Mean-Square-Errors (RMSEs) for monatomic cations, monatomic anions, and polyatomic ions were 0.29, 0.25, and 0.12 kcal/mol, respectively. The corresponding RMSE for 1744 CCSD(T)/aVTZ three-body (W–W–W) energies, used to train MB-pol, was 0.12 kcal/mol. The accuracy of the proposed classical model demonstrates that the three-body term for aqueous ionic systems can be accurately modeled classically. This approach provides a fast, efficient, and as-accurate path toward modeling the three-body term in aqueous ionic systems that is fully transferable across systems with different ions without the need to fit to tens of thousands of ab initio calculations for each ion to extend existing many-body force fields to interactions between water and ions.

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

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

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

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

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

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

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

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

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

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

Fluids↗
DOE OSTI2022

Formation of the N≡N Triple Bond from Reductive Coupling of a Paramagnetic Diruthenium Nitrido Compound

Construction of nitrogen–nitrogen triple bonds via homocoupling of metal nitrides is an important fundamental reaction relevant to a potential Nitrogen Economy. Here, we report that room temperature photolysis of Ru 2 (chp) 4 N 3 (chp– = 2-chloro-6-hydroxypyridinate) in CH 2 Cl 2 produces N 2 via reductive coupling of Ru2(chp)4N nitrido species. Computational analysis reveals that the nitride coupling transition state (TS) features an out-of-plane “zigzag” geometry instead of the anticipated planar zigzag TS. However, with intentional exclusion of dispersion correction, the planar zigzag TS geometry can also be found. Both the out-of-plane and planar zigzag TS geometries feature two important types of orbital interactions: (1) donor–acceptor interactions involving intermolecular donation of a nitride lone pair into an empty Ru–N π* orbital and (2) Ru–N π to Ru–N π* interactions derived from coupling of nitridyl radicals. The relative importance of these two interactions is quantified both at and after the TS. Our analysis shows that both interactions are important for the formation of the N–N σ bond, while radical coupling interactions dominate the formation of N–N π bonds. Comparison is made to isoelectronic Ru 2 -oxo compounds. Formation of an O–O bond via bimolecular oxo coupling is not observed experimentally and is calculated to have a much higher TS energy. As a result, the major difference between the nitrido and oxo systems stems from an extremely large driving force, ~–500 kJ/mol, for N–N coupling vs a more modest driving force for O–O coupling, –40 to –140 kJ/mol.

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

Competing Intermolecular and Molecule–Surface Interactions: Dipole–Dipole-Driven Patterns in Mixed Carborane Self-Assembled Monolayers

Carboranedithiol isomers adsorbing with opposite orientations of their dipoles on surfaces are self-assembled together to form mixed monolayers where both lateral dipole–dipole and lateral thiol–thiolate (S–H···S) interactions provide enhanced stability over single-component monolayers. Here, we demonstrate the first instance of the ability to map individual isomers in a mixed monolayer using the model system carboranedithiols on Au{111}. The addition of methyl groups to one isomer provides both an enhanced dipole moment and extra apparent height for differentiation via scanning tunneling microscopy (STM). Associated computational investigations rationalize favorable interactions of mixed pairs and the associated stability changes that arise from these interactions. Both STM images and Monte Carlo simulations yield similarly structured mixed monolayers, where approximately 10% of the molecules have reversed dipole moment orientations but no direct chemical attachment to the surface, leading to homogeneous monolayers with no apparent phase separation. Deprotonating the thiols by depositing the molecules under basic conditions eliminates the lateral S–H···S interactions while accentuating the dipole–dipole forces. The molecular system investigated is composed of isomeric molecules with opposite orientations of dipoles and identical surface packing, which enables the mapping of individual molecules within the mixed monolayers and enables analyses of the contributions of the relatively weak lateral interactions to the overall stability of the assemblies.

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

A new metric to control nucleation and grain size distribution in hybrid organic–inorganic perovskites by tuning the dielectric constant of the antisolvent

In perovskite research, there is a widely exploited but poorly explained phenomenon in which the addition of “antisolvents (ATS)” to precursor solutions results in higher-quality films. We explain the mechanism and driving force underlying an antisolvent-driven solvent extraction process. Density functional theory calculations uncover the defining effects of antisolvent choice on the extent of complexation between a lead salt and a methylammonium cation in solution. We experimentally validate the computational results using ultraviolet-visible spectroscopy and 207 Pb nuclear magnetic spectroscopy of methylammonium lead iodide solutions, containing both a processing solvent and an antisolvent. Furthermore, we uncover, and subsequently identify, the appearance of new species in solution as a result of the addition of the antisolvent. Furthermore, we observe that the choice of antisolvent has a substantial effect on the nature of the complexation of the methylammonium lead iodide (MAPbI 3 ) precursor species, whose origin we explain at an atomic level; specifically, the lower the dielectric of the antisolvent, the stronger the intermolecular binding energy between methylammonium cation (MA + ) cation and PbI 3 - plumbate, independent of the solvent or antisolvent interaction with the lead salt. Thin films were characterized using scanning electron microscopy; images of the films show how the addition of an antisolvent influences and, importantly, can be used to alter thin-film grain size. Grain size and distribution in thin films is reflected by the choice of antisolvent, promoting slower nucleation rates, a lower nucleation density, and hence larger final grain size.

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

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗