Mechanisms of solvation dynamics of polyatomic solutes in polar and nondipolar solvents: A simulation study
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523 (United States)
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802 (United States)
Molecular dynamics (MD) simulations of a benzenelike solute in acetonitrile and CO{sub 2} (298 K and 52.18thinspcm{sup 3}/mol) are used to investigate the molecular basis of solvation dynamics in polar and nondipolar solvents. The solvation response to various charge rearrangements within the benzene solute are simulated in order to mimic the type of electrostatic solvation observed in typical experimental systems. From equilibrium MD simulations the solvation time correlation function [TCF; C(t)] and the corresponding solvation velocity TCF [G(t)] are used to study the mechanisms underlying time-dependent solvation within the linear response limit. Decomposition of G(t) into contributions from rotational and translational solvent velocities reveals that the relative mix of these two types of motion is quite similar in the two solvents but is strongly dependent on the multipolar order ({ital m}) of the solute perturbation. The contribution of translational solvent motions to both the short and long time dynamics of C(t) increases from about 10{percent} for a monopolar perturbation (m=0; i.e., a change in net charge) to about 40{percent} for a perturbation of octopolar (m=3) symmetry. Decomposition of both C(t) and G(t) into single-molecule and molecular-pair contributions shows that the collective nature of the solvation response depends markedly on the charge symmetry of both the solvent molecule{close_quote}s charge distribution and the solute perturbation. In the nondipolar solvent CO{sub 2} neither C(t) nor G(t) differ significantly from their single-molecule counterparts{emdash}collective effects are therefore of little consequence to solvation in this solvent. However, in the highly dipolar solvent acetonitrile pair contributions to C(t) greatly suppress the magnitude of the solvation response and as a consequence greatly increase the speed of the response over what it would be in their absence. The importance of these intermolecular correlations in acetonitrile decreases substantially with {ital m}, such that the {open_quotes}suppression factors{close_quotes} ({alpha}{sub s}) vary from {approximately}9 for m=0 to {approximately}2 for m=3. The intermolecular correlations of primary importance in acetonitrile are of a static rather than a dynamic nature (i.e., pair effects on G(t) are of only secondary importance). This feature makes it possible to employ several approximate relationships to relate the collective dynamics of solvation in polar fluids to simpler single-solvent molecule dynamics. {copyright} {ital 1998 American Institute of Physics.}
- OSTI ID:
- 639166
- Journal Information:
- Journal of Chemical Physics, Vol. 109, Issue 8; Other Information: PBD: Aug 1998
- Country of Publication:
- United States
- Language:
- English
Similar Records
Three pulse photon echo studies of nondipolar solvation: Comparison with a viscoelastic model
Modeling solvation contributions to conformational free energy changes of biomolecules using a potential of mean force expansion