Simulations of solvation free energies and solubilities in supercritical solvents

Computer simulations are used to study solvation free energies and solubilities in supercritical solvents. Solvation free energies are calculated using the particle insertion method. The equilibrium solvent configurations required for these calculations are based on molecular dynamics simulations employing model solvent potentials previously tuned to reproduce liquid-vapor coexistence properties of the fluids Xe, C(2)H(6), CO(2), and CHF(3). Solutes are represented by all-atom potentials based on ab initio calculations and the OPLS-AA parameter set. Without any tuning of the intermolecular potentials, such calculations are found to reproduce the solvation free energies of a variety of typical solid solutes with an average accuracy of +/-2 kJmol. Further calculations on simple model solutes are also used to explore general aspects of solvation free energies in supercritical solvents. Comparisons of solutes in Lennard-Jones and hard-sphere representations of Xe show that solvation free energies and thus solubilities are not significantly influenced by solvent density fluctuations near the critical point. The solvation enthalpy and entropy do couple to these fluctuations and diverge similarly to solute partial molar volumes. Solvation free energies are also found to be little affected by the local density augmentation characteristic of the compressible regime. In contrast to solute-solvent interaction energies, which often provide a direct measure of local solvent densities, solvation free energies are remarkably insensitive to the presence of local density augmentation.     

A novel approach to solvation time scale in nonpolar solvents via instability of solvent density modes

We perform linear stability analysis of solvent density modes in the presence of nonpolar solute-solvent interaction in a nonpolar solvent. The dominant instability given by the maximum positive eigenvalue of the stability matrix provides the time scale of the solvent rearrangement around a solute. Our theory predicts two long time scales for both in normal nonpolar and supercritical fluids. We discuss the existing experimental results on nonpolar solvation dynamics in light of our prediction.     

The effects of solute-solvent electrostatic interactions on solvatochromic shifts in supercritical CO2

Solvent clustering around attractive solutes is an important feature of supercritical solvation. We examine here the effects of the local density enhancement on solvatochromic shifts in electronic absorption and emission spectra in supercritical CO2. We use molecular dynamics (MD) simulation to study the spectral line shifts for model diatomic solutes that become more polar upon electronic excitation. The electronic transition is modeled as either a change from a quadrupolar to a dipolar solute charge distribution or as an increase in the magnitude of the solute dipole. Our main focus is on the density dependence of the line shifts at 320 K, which corresponds to about 1.05 times the solvent critical temperature, Tc, but results for higher temperatures are also obtained in order to determine the behavior of the line shifts in the absence of local density enhancement. We find that the extent of local density enhancement at 1.05Tc is strongly correlated with solute-solvent electrostatic attraction and that the density dependence of the emission line shifts resembles the behavior of the effective local densities, rho(eff), obtained from the first-shell coordination numbers. The differences that are seen are shown to be due to solute-solvent orientational correlations which provide an additional source of enhancement for electrostatic solvation energies and spectral line shifts.     

Solvation of coumarin 153 in supercritical fluoroform

We present a study of local density augmentation around an attractive solute (i.e., giving rise to more attractive interaction with the solvent than solvent-solvent interactions) in supercritical fluoroform. This work is based on molecular dynamics simulations of coumarin 153 in supercritical fluoroform at densities both above and below the critical density, ranging from dilute gas-like to liquid-like, at a reduced temperature (T/T(c)) of 1.03. We focused on studying the structure of the solvation shell and the variation of the solute electronic absorption and emission shifts with density. Quantum calculations at the density functional theory (DFT) level were run on the solute in the ground state, and time-dependent DFT calculations were performed in the solute excited state in order to determine the solute-solvent potential parameters. The results obtained for the Stokes shift are in agreement with the experimental measurements. To evaluate local density augmentation from simulations, we used two different definitions, one based on the solvation number and the other derived from solvatochromic shifts. In the former case, the agreement with experimental results is good, while, in the latter case, better agreement is achieved by perturbatively including the induced-dipole contribution to the solvation energy.     

Free energy of solvation from molecular dynamics simulation applying Voronoi-Delaunay triangulation to the cavity creation

The free energy of solvation for a large number of representative solutes in various solvents has been calculated from the polarizable continuum model coupled to molecular dynamics computer simulation. A new algorithm based on the Voronoi-Delaunay triangulation of atom-atom contact points between the solute and the solvent molecules is presented for the estimation of the solvent-accessible surface surrounding the solute. The volume of the inscribed cavity is used to rescale the cavitational contribution to the solvation free energy for each atom of the solute atom within scaled particle theory. The computation of the electrostatic free energy of solvation is performed using the Voronoi-Delaunay surface around the solute as the boundary for the polarizable continuum model. Additional short-range contributions to the solvation free energy are included directly from the solute-solvent force field for the van der Waals-type interactions. Calculated solvation free energies for neutral molecules dissolved in benzene, water, CCl4, and octanol are compared with experimental data. We found an excellent correlation between the experimental and computed free energies of solvation for all the solvents. In addition, the employed algorithm for the cavity creation by Voronoi-Delaunay triangulation is compared with the GEPOL algorithm and is shown to predict more accurate free energies of solvation, especially in solvents composed by molecules with nonspherical molecular shapes.     

Solvent density effects on the solvation behavior and configurational structure of bare and passivated 38-atom gold nanoparticle in supercritical ethane

In exploring the effects of solvent density on the mode and the degree of solvation of the bare and passivated 38-atom gold particle in supercritical ethane, we have extended the molecular dynamics simulations of the system, reported previously,(34) to cover a range of isotherms in the T > T(c) regime, where T(c) is the critical temperature of the solvent. Consonant with our previous observations, the modes of solvation of the bare and the passivated particle, deduced from the radial distribution of the solvent about the metal core center of mass, are found to be vastly different from each other at all solvent densities: while the molecules solvating the bare particle form a well-defined, two-region layer around it, those solvating the passivated particle are loosely dispersed in the passivating layer. For the bare particle, the degree of solvation (vartheta) as a function of solvent density passes through a maximum occurring in the close vicinity of the critical point, consistent with our previous results and in agreement with Debenedetti's theoretical analysis,(22,23) which predicts a solvation enhancement effect in the critical region for systems where the unlike solvent/solute interaction is much stronger than the solvent/solvent interaction. Taking the degree of solvation (vartheta) as a measure of solvent quality, we have investigated how the solvent quality would vary along the solvent-density isotherms. In the solvent-density regime rho > rho(c), the solvent quality is found to be a decreasing function of the density as a result of progressive dominance of the excluded volume effect over the attractive particle/solvent interactions. The particle/solvent affinity is greatly reduced in the presence of the passivating layer, resulting in considerable shrinkage of the good-solvent-quality domain in the supercritical regime. The solvent environment and the presence of the passivating chains produce significant disorder in the equilibrium structure assumed by the nanoparticle core.     

用离散-连续模型计算NH2-,NH3和NH4+的溶剂化自由能

通过理论计算推测NH2-,NH3和NH4+在水溶液第一溶剂化层中与之直接作用的水分子分别为2,4和4个,并采用离散-连续模型计算了NH2-,NH3,NH3和NH4+在水溶液中的溶剂化自由能.结果表明,由于离散-连续模型在从头算水平考虑了溶质分子与第一溶剂化层溶剂分子之间的作用,能更准确地描述溶剂化作用.此外,采用更加符合溶液中真实情况的溶剂化构型,能得到更准确的溶剂化性质.     

Solvation in supercritical water

Solvation in supercritical water under equilibrium and nonequilibrium conditions is studied via molecular dynamics simulations. The influence of solute charge distributions and solvent density on the solvation structures and dynamics is examined with a diatomic probe solute molecule. It is found that the solvation structure varies dramatically with the solute dipole moment, especially in low-density water, in accord with many previous studies on ion solvation. This electrostrictive effect has important consequences for solvation dynamics. In the case of a nonequilibrium solvent relaxation, if there are sufficiently many water molecules close to the solute at the outset of the relaxation, the solvent response measured as a dynamic Stokes shift is almost completely governed by inertial rotations of these water molecules. By contrast, in the opposite case of a low local solvent density near the solute, not only rotations but also translations of water molecules play an important role in solvent relaxation dynamics. The applicability of a linear response is found to be significantly restricted at low water densities.     

Pyrazine in supercritical xenon: local number density defined by experiment and calculation

Toward our goal of using supercritical fluids to study solvent effects in physical and chemical phenomena, we develop a method to spatially define the solvent local number density at the solute in the highly compressible regime of a supercritical fluid. Experimentally, the red shift of the pyrazine n-pi* electronic transition was measured at high dilution in supercritical xenon as a function of pressure from 0 to approximately 24 MPa at two temperatures: one (293.2 K) close to the critical temperature and the other (333.2 K) remote. Computationally, several representative stationary points were located on the potential surfaces for pyrazine and 1, 2, 3, and 4 xenons at the MP2/6-311++G(d,p)/aug-cc-pVTZ-PP level. The vertical n-pi* ((1)B(3u)) transition energies were computed for these geometries using a TDDFT/B3LYP/DGDZVP method. The combination of experiment and quantum chemical computation allows prediction of supercritical xenon bulk densities at which the pyrazine primary solvation shell contains an average of 1, 2, 3, and 4 xenon molecules. These density predictions were achieved by graphical superposition of calculated shifts on the experimental shift versus density curves for 293.2 and 333.2 K. Predicted bulk densities are 0.50, 0.91, 1.85, and 2.50 g cm(-3) for average pyrazine primary solvation shell occupancy by 1, 2, 3, and 4 xenons at 293.2 K. Predicted bulk densities are 0.65, 1.20, 1.85, and 2.50 g cm(-3) for average pyrazine primary solvation shell occupancy by 1, 2, 3, and 4 xenons at 333.2 K. These predictions were evaluated with classical Lennard-Jones molecular dynamics simulations designed to replicate experimental conditions at the two temperatures. The average xenon number within 5.0 A of the pyrazine center-of-mass at the predicted densities is 1.3, 2.1, 3.0, and 4.0 at both simulation temperatures. Our three-component method-absorbance measurement, quantum chemical prediction, and evaluation of prediction with classical molecular dynamics simulation-therefore has a high degree of internal consistency for a system in which the intermolecular interactions are dominated by dispersion forces.     

Solvent effects in density functional calculations of uracil and cytosine tautomerism

The effects of the solvent on the tautomeric equilibria of cytosine and uracil are studied using Onsager's reaction field model in the framework of density functional theory. In this model, the solute molecule is placed in a spherical cavity of radius ao immersed in a continuous medium with a dielectric constant e, treating the solute-solvent electrostatic interactions at the dipole level. The cavity radius is evaluated by determining the molecular volume with a quantum mechanical approach. The solvent effect brings significant changes in the geometrical parameters of some cytosine tautomers, but only very small changes to those of the uracil tautomers. Our results are in good agreement with available experimental results and confirm that the polarization of the solute by the continuum has important effects on the absolute and relative solvation energies. Frequency shifts and intensity variations in the infrared spectra due to the presence of the solvent are also presented. © 1995 John Wiley & Sons, Inc.     

The effects of solute-solvent electrostatic interactions on solvation dynamics in supercritical CO2

We present here the results of molecular-dynamics simulation of solvation dynamics in supercritical CO(2) at a temperature of about 1.05T(c), where T(c) is the critical temperature, and at a series of densities ranging from 0.4 to 2.0 of the critical density rho(c). We focus on electrostatic solvation dynamics, representing the electronic excitation of the chromophore as a change in its charge distribution from a quadrupolar-symmetry ground state to a dipolar excited state. Two perturbations are considered, corresponding to different magnitudes of solute excited-state dipoles, denoted as d5 and d8. The d8 solute is more attractive, leading to a larger enhancement in CO(2) clustering upon solute electronic excitation. This has a large impact on solvation dynamics, especially at densities below rho(c). At these densities, solvation dynamics is much slower for the d8 than for the d5 solute. For both solutes, solvation dynamics becomes faster at densities above rho(c) at which solvent clustering diminishes. We show that the slowest solvation time scale is associated with solvent clustering and we relate it to solute-solvent mutual translational diffusion and the extent of change in effective local density resulting from solute electronic excitation.     

Free-energy analysis of the electron-density fluctuation in the quantum-mechanical/molecular-mechanical simulation combined with the theory of energy representation

The relationship is investigated for QM/MM (quantum-mechanical/molecular-mechanical) systems between the fluctuations of the electronic state of the QM subsystem and of the solvation effect due to the QM-MM interaction. The free-energy change due to the electron-density fluctuation around its average is highlighted, and is evaluated through an approximate functional formulated in terms of distribution functions of the many-body coupling (pairwise non-additive) part of the QM-MM interaction energy. A set of QM/MM simulations are conducted in MM water solvent for QM water solute in ambient and supercritical conditions and for QM glycine solute in the neutral and zwitterionic forms. The variation of the electronic distortion energy of the QM solute in the course of QM/MM simulation is then shown to be compensated by the corresponding variation of the free energy of solvation. The solvation free energy conditioned by the electronic distortion energy is further analyzed with its components. It is found that the many-body contribution is essentially equal between the free energy and the average sum of solute-solvent interaction energy.     

Ultrafast solvent response upon a change of the solute size in non-polar supercritical fluids

Non-polar solvation dynamics has been investigated using steady-state absorption and emission spectroscopy of the NO A ²Σ⁺(3sσ) Rydberg state in fluid Ar over a wide range of densities spanning the supercritical regime. Equilibrium molecular dynamics simulations were implemented to derive a new isotropic NO A(3sσ)–Ar pair potential which was further used to investigate the role of local density enhancements on the solvation process by non-equilibrium molecular dynamics simulations. These density inhomogeneities were found to have no influence on the solvation dynamics. Furthermore, the latter was shown to take place in a strongly non-linear regime, especially at low temperatures. This process results from the dramatic change of solute–solvent short range interaction associated with the large solute size change upon excitation to the Rydberg state.     

Solvation model based on order parameters and a fast sampling method for the calculation of the solvation free energies of peptides

An analytical solvation model is proposed as a function of an order parameter, which represents the local arrangement of water molecules in the first solvation shell of peptide atoms. The model is combined with a fast sampling method, rotational isomeric state Monte Carlo, to sample efficiently the torsional degrees of freedom on a peptide backbone. This order parameter solvation model is shown to reproduce without ad hoc fitting parameters the solvation free energies of single amino acids and tripeptides with slightly better accuracy than the generalized Born model but with several orders of magnitude improvement in efficiency. This method is a potential candidate for efficiently and accurately tackling some important issues in biophysical chemistry that are related to solvation, for example, protein folding, ligand binding, etc. Our results also present fundamental new insights into solvation. Specifically, the local water geometry, represented in this work by a properly defined order parameter, carries the majority, if not all, of the energetic information of solvation, including solute-solvent interactions and solvent reorganization in the presence of the solute.     

A novel algorithm for creating coarse-grained, density dependent implicit solvent models

Implicit solvent simulations are those in which solvent molecules are not explicitly simulated, and the solute-solute interaction potential is modified to compensate for the implicit solvent effect. Implicit solvation is well known in Brownian dynamics of dilute solutions but offers promise to speed up many other types of molecular simulations as well, including studies of proteins and colloids where the local density can vary considerably. This work examines implicit solvent potentials within a more general coarse-graining framework. While a pairwise potential between solute sites is relatively simple and ubiquitous, an additional parametrization based on the local solute concentration has the possibility to increase the accuracy of the simulations with only a marginal increase in computational cost. We describe here a method in which the radial distribution function and excess chemical potential of solute insertion for a system of Lennard-Jones particles are first measured in a fully explicit, all-particle simulation, and then reproduced across a range of solute particle densities in an implicit solvent simulation.     

Molecular dynamics study on hydrophobic effects in aqueous urea solutions.

Hydrophobic effects in aqueous urea were analyzed by molecular dynamics simulations. The contribution of solvents to the potential of mean force between two methane molecules was calculated by using molecular dynamics simulations and was compared with the solubility data of hydrocarbons in aqueous urea. Both the simulation results and the solubility data indicated that urea stabilizes methane-methane association. The stabilization was due to increasing the solvation free energies of small hydrocarbons such as methane by addition of urea. The solvation free energies of larger hydrocarbons, on the other hand, are decreased by addition of urea. This effect of the solute size on hydrophobic free energies in aqueous urea was also analyzed by using molecular dynamics simulations by means of division of the solvation process into two parts: the cavity formation and the introduction of the solute-solvent attractive interactions. In the cavity formation, urea increased hydrophobic free energies, and in the introduction of the solute-solvent attractive interactions, urea decreased hydrophobic free energies. The influence of urea on hydrophobic free energies was determined by the balance of effects of the two parts of the solvation process.