On the microscopic theory of polar solvation dynamics

A microscopic model of the time-resolved Stokes shift is developed. The model calculates the solvation dynamics by combining the atomic resolution of the solute structure with dipolar dynamics from the polarization structure factors of the homogeneous solvent. Calculations are made for coumarin 153 and quinoxaline optical dyes with atomic geometries and charge distributions taken from quantum calculations. Stokes shift dynamics is calculated and compared to experiment in high-temperature acetonitrile and methanol and in low-temperature 2-methyl-tetrahydrofurane using dielectric relaxation data from experiment.     

Theory of electron solvation in polar liquids: a continuum model

The solvation of electrons in polar liquids is analyzed on the basis of an extended continuum model. In addition to the long-range electron-dipole interaction two short-range interactions are introduced. Among others one accounts for interactions with groups capable of forming hydrogen bonds and the second for quadrupolar characteristics of the liquid molecules. Both are induced by the orientation of the molecular dipole. Applying the scaling method a proper reaction coordinate is introduced and the solvation dynamics are discussed for the electron in the electronic ground state and after excitation to the p-type excited state. The observed spectral evolution of the transient absorption spectra, after two photon excitations for electrons in water and in methanol, is well described by this theory. An analytic estimate for the nonradiative deactivation from the electronically excited solvated electron is found to be consistent with an observed lifetime of 50 fs for the electron in water. The theory predicts an about three times slower internal conversion in methanol as solvent in comparison with water.     

Spontaneous vesicle self-assembly: a mesoscopic view of membrane dynamics

Amphiphilic vesicles are ubiquitous in living cells and industrially interesting as drug delivery vehicles. Vesicle self-assembly proceeds rapidly from nanometer to micrometer length scales and is too fast to image experimentally but too slow for molecular dynamics simulations. Here, we use parallel dissipative particle dynamics (DPD) to follow spontaneous vesicle self-assembly for up to 445 μs with near-molecular resolution. The mean mass and radius of gyration of growing amphiphilic clusters obey power laws with exponents of 0.85 ± 0.03 and 0.41 ± 0.02, respectively. We show that DPD provides a computational window onto fluid dynamics on scales unreachable by other explicit-solvent simulations.     

Theory of solvation in polar nematics

We develop a linear response theory of solvation of ionic and dipolar solutes in anisotropic, axially symmetric polar solvents. The theory is applied to solvation in polar nematic liquid crystals. The formal theory constructs the solvation response function from projections of the solvent dipolar susceptibility on rotational invariants. These projections are obtained from Monte Carlo simulations of a fluid of dipolar spherocylinders which can exist both in the isotropic and nematic phases. Based on the properties of the solvent susceptibility from simulations and the formal solution, we have obtained a formula for the solvation free energy which incorporates the experimentally available properties of nematics and the length of correlation between the dipoles in the liquid crystal. The theory provides a quantitative framework for analyzing the steady-state and time-resolved optical spectra and makes several experimentally testable predictions. The equilibrium free energy of solvation, anisotropic in the nematic phase, is given by a quadratic function of cosine of the angle between the solute dipole and the solvent nematic director. The sign of solvation anisotropy is determined by the sign of dielectric anisotropy of the solvent: solvation anisotropy is negative in solvents with positive dielectric anisotropy and vice versa. The solvation free energy is discontinuous at the point of isotropic-nematic phase transition. The amplitude of this discontinuity is strongly affected by the size of the solute becoming less pronounced for larger solutes. The discontinuity itself and the magnitude of the splitting of the solvation free energy in the nematic phase are mostly affected by microscopic dipolar correlations in the nematic solvent. Illustrative calculations are presented for the equilibrium Stokes shift and the Stokes shift time correlation function of coumarin-153 in 4-n-pentyl-4'-cyanobiphenyl and 4,4-n-heptyl-cyanopiphenyl solvents as a function of temperature in both the nematic and isotropic phases.     

Probing polar solvation dynamics in proteins: a molecular dynamics simulation analysis

Measurements of time-resolved Stokes shifts on picosecond to nanosecond time scales have been used to probe the polar solvation dynamics of biological systems. Since it is difficult to decompose the measurements into protein and solvent contributions, computer simulations are useful to aid in understanding the details of the molecular behavior. Here we report the analysis of simulations of the electrostatic interactions of the rest of the protein and the solvent with 11 residues of the immunoglobulin binding domain B1 of protein G. It is shown that the polar solvation dynamics are position-dependent and highly heterogeneous. The contributions due to interactions with the protein and with the solvent are determined. The solvent contributions are found to vary from negligible after a few picoseconds to dominant on a scale of hundreds of picoseconds. The origin for the latter is found to involve coupled hydration and protein conformational dynamics. The resulting microscopic picture demonstrates that a wide range of possibilities have to be considered in the interpretation of time-resolved Stokes shift measurements.     

Computer simulation of cellulose solvation in polar solvents

The solvation of cellulose molecules in water and N-methylmorpholine N-oxide has been investigated by the molecular-dynamics method. An analysis of simulation results yields the conformational characteristics of cellulose molecules in these solvents: the mean-square distance between polymer chain ends and the radial distribution function of monomer units. The radial distribution functions of oxygen atoms of solvents relative to protons of carbohydrate molecules are estimated. This step makes it possible to draw conclusions about the number and character of hydrogen bonds and the local structure of solvate shells.     

Trimethylglycine complexes with carboxylic acids and HF: solvation by a polar aprotic solvent

A series of strong H-bonded complexes of trimethylglycine, also known as betaine, with acetic, chloroacetic, dichloroacetic, trifluoroacetic and hydrofluoric acids as well as the homo-conjugated cation of betaine with trifluoroacetate as the counteranion were investigated by low-temperature (120-160 K) liquid-state NMR spectroscopy using CDF(3)/CDF(2)Cl mixture as the solvent. The temperature dependencies of (1)H NMR chemical shifts are analyzed in terms of the solvent-solute interactions. The experimental data are explained assuming the combined action of two main effects. Firstly, the solvent ordering around the negatively charged OHX region of the complex (X = O, F) at low temperatures, which leads to a contraction and symmetrisation of the H-bond; this effect dominates for the homo-conjugated cation of betaine. Secondly, at low temperatures structures with a larger dipole moment are preferentially stabilized, an effect which dominates for the neutral betaine-acid complexes. The way this second contribution affects the H-bond geometry seems to depend on the proton position. For the Be(+)COO(-)···HOOCCH(3) complex (Be = (CH(3))(3)NCH(2)-) the proton displaces towards the hydrogen bond center (H-bond symmetrisation, O···O contraction). In contrast, for the Be(+)COOH···(-)OOCCF(3) complex the proton shifts further away from the center, closer to the betaine moiety (H-bond asymmetrisation, O···O elongation). Hydrogen bond geometries and their changes upon lowering the temperature were estimated using previously published H-bond correlations.     

Shape effect on nanoparticle solvation: a comparison of morphometric thermodynamics and microscopic theories

Conventional wisdom for controlling the nanoparticle size and shape during synthesis is that particle growth favors the direction of a facet with the highest surface energy. However, the particle solvation free energy, which dictates the particle stability and growth, depends not only on the surface area and surface free energy but also on other geometric measures such as the solvent excluded volume and the surface curvature and their affiliated thermodynamic properties. In this work, we study the geometrical effects on the solvation free energies of nonspherical nanoparticles using morphometric thermodynamics and density functional theories. For idealized systems that account for only molecular excluded-volume interactions, morphometric thermodynamics yields a reliable solvation free energy when the particle size is significantly larger than the solvent correlation length. However, noticeable deviations can be identified in comparison to the microscopic theories for predicting the solvation free energies of small nanoparticles. This conclusion also holds for predicting the potential of mean force underlying the colloidal "key-and-lock" interactions. Complementary to the microscopic theories, morphometric thermodynamics requires negligible computational cost, therefore making it very appealing for a broad range of practical applications.     

Molecular density functional theory of solvation: from polar solvents to water

A classical density functional theory approach to solvation in molecular solvent is presented. The solvation properties of an arbitrary solute in a given solvent, both described by a molecular force field, can be obtained by minimization of a position and orientation-dependent free-energy density functional. In the homogeneous reference fluid approximation, limited to two-body correlations, the unknown excess term of the functional approximated by the angular-dependent direct correlation function of the pure solvent. We show that this function can be extracted from a preliminary MD simulation of the pure solvent by computing the angular-dependent pair distribution function and solving subsequently the molecular Ornstein-Zernike equation using a discrete angular representation. The corresponding functional can then be minimized in the presence of an arbitrary solute on a three-dimensional cubic grid for positions and Gauss-Legendre angular grid for orientations to provide the solvation structure and free-energy. This two-step procedure is proved to be much more efficient than direct molecular dynamics simulations combined to thermodynamic integration schemes. The approach is shown to be relevant and accurate for prototype polar solvents such as the Stockmayer solvent or acetonitrile. For water, although correct for neutral or moderately charged solute, it tends to underestimate the tetrahedral solvation structure around H-bonded solutes, such as spherical ions. This can be corrected by introducing suitable three-body correlation terms that restore both an accurate hydration structure and a satisfactory energetics.     

Picosecond observations of electron solvation in dilute polar fluids

Electron solvation has been studied in dilute polar fluids in order to quantify the role of the fluid in the proposed mechanisms of electron trapping and solvation. In a series of dilute alcohol-alkane systems, the picosecond evolution of the absorption spectrum is shown to be a sensitive function of the local liquid structure and dynamics. A solvation mechanism is outlined which correlates the absorption and mobility data from neat and dilute polar fluids.     

Spectroscopic studies of anion complexes and clusters: a microscopic approach to understanding anion solvation

Recent infrared spectroscopic studies of negatively charged clusters in the gas phase have furnished new information on non-covalent bonds between anions and neutral molecules, and provided fresh perspectives on the microscopic details of anion solvation. We describe the central spectroscopic techniques employed for obtaining infrared spectra of mass-selected solvated anions in the gas phase, and illustrate recent progress by describing studies of simple halide-H2 dimers, and larger clusters in which up to 9 C2H2 molecules are attached to a Cl- anion.     

Dynamics of polar solvation: Route to single exponential relaxation via translational diffusion

A microscopic theoretical calculation of time-dependent solvation energy shows that the solvation of an ion or a dipole is dominated by a single relaxation time if the translational contribution to relaxation is significant. Contribution No. 545 from the Solid State and Structural Chemistry Unit     

COSMO-RS: A novel view to physiological solvation and partition questions

Both, dielectric continuum solvation models as well as surface or group based methods using polarity and lipophilicity parameters have been proven to be useful tools for the analysis of solvation and partition questions. For the first time, COSMO-RS provides an integrated theory, which combines the aspects of continuum solvation and surface interactions, and which ends up with chemical potentials of molecules in almost arbitrary solvents and mixtures. Due to its sound theoretical basis, COSMO-RS does not only provide a new quantitative access to solvation and partition properties in well defined solvents, but it also opens a novel view and gives a better understanding of the general problem of solvation. Finally, this allows for a generalisation of COSMO-RS to sophisticatedphysiological partition problems involving as complex phases as blood, brain, or cell membranes. The use of COSMO-RS for drug discovery and design is demonstrated by applications to blood-brain partition coefficients, and water solubility.     

Solvent reorganization entropy of electron transfer in polar solvents

We report the results of molecular dynamics simulations of the solvent reorganization energy of intramolecular electron transfer in a charge-transfer molecule dissolved in water and acetonitrile at varying temperatures. The simulations confirm the prediction of microscopic solvation theories of a positive reorganization entropy in polar solvents. The results of simulations are analyzed in terms of the splitting of the reorganization entropy into the contributions from the solute-solvent interaction and from the alteration of the solvent structure induced by the solute. These two contributions mutually cancel each other, resulting in the reorganization entropy amounting to only a fraction of each component.     

The dielectric constant of a polar liquid by the simulation of microscopic drops

We report a new method of computing the static dielectric constant for a polar liquid using a natural system, a drop of liquid in its own vapour. Results for the Stockmayer potential agree with those computed for a homogeneous “cyclic” liquid with Ewald corrections. The method is quite general.     

A REMPI study of solvation of small Hg n clusters by polar molecules

Free Hg _n (DME) _m clusters (where DME=dimethyl-ether,n=1, 2, 3,m=1÷5) formed in a supersonic expansion were studied by the REMPI (Resonance-Enhanced Multi-Photon Ionization) technique. A large decrease of ionization energies due to solvation of Hg _n clusters is observed. Preliminary results are discussed in terms of different equilibrium configurations of the electronic ground, excited and ionic states of clusters.     

Inference of the solvation energy parameters of amino acids using maximum entropy approach

We present a novel technique, based on the principle of maximum entropy, for deriving the solvation energy parameters of amino acids from the knowledge of the solvent accessible areas in experimentally determined native state structures as well as high quality decoys of proteins. We present the results of detailed studies and analyze the correlations of the solvation energy parameters with the standard hydrophobic scale. We study the ability of the inferred parameters to discriminate between the native state structures of proteins and their decoy conformations.