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.     

Absorption and emission lineshapes and solvation dynamics of NO in supercritical Ar

We perform a theoretical study of electronic spectroscopy of dilute NO in supercritical Ar fluid. Absorption and emission lineshapes for the A(2)Sigma(+)<--X(2)Pi Rydberg transition of NO in argon have been previously measured and simulated, which yielded results for the NO/Ar ground- and excited-state pair potentials [Larregaray et al., Chem. Phys. 308, 13 (2005)]. Using these potentials, we have performed molecular dynamics simulations and theoretical statistical mechanical calculations of absorption and emission lineshapes and nonequilibrium solvation correlation functions for a wide range of solvent densities and temperatures. Theory was shown to be in good agreement with simulation. Linear response treatment of solvation dynamics was shown to break down at near-critical temperature due to dramatic change in the solute-solvent microstructure upon solute excitation to the Rydberg state and the concomitant increase of the solute size.     

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.     

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.     

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.     

Polar solvation dynamics in supercritical fluids: a mode-coupling treatment

A mode-coupling treatment of polar solvation dynamics in supercritical fluids is presented. The equilibrium solvation time correlation function for the solute fluctuating transition frequency is obtained from the mode-coupling theory method and from molecular-dynamics simulations. The theory is shown to be in good agreement with the simulation. The solvation time correlation function exhibits three distinct time scales, with rapid initial decay, followed by a recurrence at intermediate times, and a slowly decaying long-time tail. Our theoretical analysis shows that the short-time decay arises from the coupling of the solute energy gap to the solvent polarization modes, the recurrence at intermediate times is due to the energy modes, while the slow long-time decay reflects the coupling to the number density modes.     

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 dynamics in acetonitrile: a study incorporating solute electronic response and nuclear relaxation

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.     

Solvation dynamics of Hoechst 33258 in water: an equilibrium and nonequilibrium molecular dynamics study

Integrated within an appropriate theoretical framework, molecular dynamics (MD) simulations are a powerful tool to complement experimental studies of solvation dynamics. Together, experiment, theory, and simulation have provided substantial insight into the dynamic behavior of polar solvents. MD investigations of solvation dynamics are especially valuable when applied to the heterogeneous environments found in biological systems, where the calculated response of the environment to the electrostatic perturbation of the probe molecule can easily be decomposed by component (e.g., aqueous solvent, biomolecule, ions), greatly aiding the molecular-level interpretation of experiments. A comprehensive equilibrium and nonequilibrium MD study of the solvation dynamics of the fluorescent dye Hoechst 33258 (H33258) in aqueous solution is presented. Many fluorescent probes employed in experimental studies of solvation dynamics in biological systems, such as the DNA minor groove binder H33258, have inherently more conformational flexibility than prototypical fused-ring chromophores. The role of solute flexibility was investigated by developing a fully flexible force-field for the H33258 molecule and by simulating its solvation response. While the timescales for the total solvation response calculated using both rigid (0.16 and 1.3 ps) and flexible (0.17 and 1.4 ps) models of the probe closely matched the experimentally measured solvation response (0.2 and 1.2 ps), there were subtle differences in the response profiles, including the presence of significant oscillations for the flexible probe. A decomposition of the total response of the flexible probe revealed that the aqueous solvent was responsible for the overall decay, while the oscillations result from fluctuations in the electrostatic terms in the solute intramolecular potential energy. A comparison of equilibrium and nonequilibrium approaches for the calculation of the solvation response confirmed that the solvation dynamics of H33258 in water is well-described by linear response theory for both rigid and flexible models of the probe.     

Generalized Langevin theory on the dynamics of simple fluids under external fields

A theory on the time development of the density and current fields of simple fluids under an external field is formulated through the generalized Langevin formalism. The theory is applied to the linear solvation dynamics of a fixed solute regarding the solute as the external field on the solvent. The solute-solvent-solvent three-body correlation function is taken into account through the hypernetted-chain integral equation theory, and the time correlation function of the random force is approximated by that in the absence of the solute. The theoretical results are compared with those of molecular-dynamics (MD) simulation and the surrogate theory. As for the transient response of the density field, our theory is shown to be free from the artifact of the surrogate theory that the solvent can penetrate into the repulsive core of the solute during the relaxation. We have also found a large quantitative improvement of the solvation correlation function compared with the surrogate theory. In particular, the short-time part of the solvation correlation function is in almost perfect agreement with that from the MD simulation, reflecting that the short-time expansion of the theoretical solvation correlation function is exact up to t(2) with the exact three-body correlation function. A quantitative improvement is found in the long-time region, too. Our theory is also applied to the force-force time correlation function of a fixed solute, and similar improvement is obtained, which suggests that our present theory can be a basis to improve the mode-coupling theory on the solute diffusion.     

A molecular dynamics computer simulation study of room-temperature ionic liquids. II. Equilibrium and nonequilibrium solvation dynamics

The molecular dynamics (MD) simulation study of solvation structure and free energetics in 1-ethyl-3-methylimidazolium chloride and 1-ethyl-3-methylimidazolium hexafluorophosphate using a probe solute in the preceding article [Y. Shim, M. Y. Choi and H. J. Kim, J. Chem. Phys. 122, 044510 (2005)] is extended to investigate dynamic properties of these liquids. Solvent fluctuation dynamics near equilibrium are studied via MD and associated time-dependent friction is analyzed via the generalized Langevin equation. Nonequilibrium solvent relaxation following an instantaneous change in the solute charge distribution and accompanying solvent structure reorganization are also investigated. Both equilibrium and nonequilibrium solvation dynamics are characterized by at least two vastly different time scales--a subpicosecond inertial regime followed by a slow diffusive regime. Solvent regions contributing to the subpicosecond nonequilibrium relaxation are found to vary significantly with initial solvation configurations, especially near the solute. If the solvent density near the solute is sufficiently high at the outset of the relaxation, subpicosecond dynamics are mainly governed by the motions of a few ions close to the solute. By contrast, in the case of a low local density, solvent ions located not only close to but also relatively far from the solute participate in the subpicosecond relaxation. Despite this difference, linear response holds reasonably well in both ionic liquids.     

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.     

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.     

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.     

Solvent density inhomogeneities and solvation free energies in supercritical diatomic fluids: a density functional approach

Density functional theory is used to explore the solvation properties of a spherical solute immersed in a supercritical diatomic fluid. The solute is modeled as a hard core Yukawa particle surrounded by a diatomic Lennard-Jones fluid represented by two fused tangent spheres using an interaction site approximation. The authors' approach is particularly suitable for thoroughly exploring the effect of different interaction parameters, such as solute-solvent interaction strength and range, solvent-solvent long-range interactions, and particle size, on the local solvent structure and the solvation free energy under supercritical conditions. Their results indicate that the behavior of the local coordination number in homonuclear diatomic fluids follows trends similar to those reported in previous studies for monatomic fluids. The local density augmentation is particularly sensitive to changes in solute size and is affected to a lesser degree by variations in the solute-solvent interaction strength and range. The associated solvation free energies exhibit a nonmonotonous behavior as a function of density for systems with weak solute-solvent interactions. The authors' results suggest that solute-solvent interaction anisotropies have a major influence on the nature and extent of local solvent density inhomogeneities and on the value of the solvation free energies in supercritical solutions of heteronuclear molecules.     

Transient cavities and the excess chemical potentials of hard-spheroid solutes in dipolar hard-sphere solvents

Monte Carlo computer simulations are used to study transient cavities and the solvation of hard-spheroid solutes in dipolar hard-sphere solvents. The probability distribution of spheroidal cavities in the solvent is shown to be well described by a Gaussian function, and the variations of fit parameters with cavity elongation and solvent properties are analyzed. The excess chemical potentials of hard-spheroid solutes with aspect ratios x in the range of 15< or =x< or =5, and with volumes between 1 and 20 times that of a solvent molecule, are presented. It is shown that for a given molecular volume and solvent dipole moment (or temperature) a spherical solute has the lowest excess chemical potential and hence the highest solubility, while a prolate solute with aspect ratio x should be more soluble than an oblate solute with aspect ratio 1x. For a given solute molecule, the excess chemical potential increases with increasing temperature; this same trend can be observed in hydrophobic solvation. A scaled-particle theory based on the solvent equation of state and a fitted solute-solvent interfacial tension shows excellent agreement with the simulation results over the whole range of solute elongations and volumes considered. An information-theoretic model based on the solvent density and radial distribution function is less successful, being accurate only for small solute volumes and low solvent densities.     

丙酮分子n→π%*跃迁光谱的平均溶剂静电势/分子动力学方法研究

采用量子力学/分子动力学方法研究了具体溶剂分子结构对溶质光谱行为的静电影响. 通过拟合溶质所处外电场和引入溶剂分子极化率, 考虑了溶质溶剂分子之间的相互极化效应, 得到合理的溶质和溶剂分子的电荷分布. 经过严格推导发现, 在传统的显溶剂模型中, 平衡和非平衡溶剂化能表达式均未考虑溶剂分子永久偶极弹簧能的贡献. 因此, 在正确计算永久偶极弹簧能的基础上, 重新建立了溶剂化能的表达式和新的吸收/发射光谱移动公式. 采用修改后的ASEP/MD程序, 计算得到了与实验值比较吻合的丙酮在水溶液中n→π*跃迁的光谱移动值, 验证了新公式的合理性.     

A case against large cluster-induced nonequilibrium solvent effects in supercritical water

Using a partially compressible continuum solvation model, we have shown that solvent compression in just the first two solvation shells (or thereabouts) is all that is required to gain the bulk of the compression-induced enhancement to the solvation energy of ions in supercritical water. This result is found to hold even when the direct, equilibrium solvent-solute cluster involves well over a hundred solvent molecules. We argue that, for charge variation reactions in supercritical water, the observed short-range behavior of the compression-induced solvation free energy precludes the existence of any anomalously large nonequilibrium solvent effects which might be expected on the basis of the very large size of the equilibrium clusters. Received: 8 January 1997 / Accepted: 17 January 1997