New approach to free energy of solvation applying continuum models to molecular dynamics simulation

A new approach to the calculation of the free energy of solvation from trajectories obtained by molecular dynamics simulation is presented. The free energy of solvation is computed as the sum of three contributions originated at the cavitation of the solute by the solvent, the solute-solvent nonpolar (repulsion and dispersion) interactions, and the electrostatic solvation of the solute. The electrostatic term is calculated based on ideas developed for the broadly used continuum models, the cavitational contribution from the excluded volume by the Claverie-Pierotti model, and the Van der Waals term directly from the molecular dynamics simulation. The proposed model is tested for diluted aqueous solutions of simple molecules containing a variety of chemically important functions: methanol, methylamine, water, methanethiol, and dichloromethane. These solutions were treated by molecular dynamics simulations using SPC/E water and the OPLS force field for the organic molecules. Obtained free energies of solvation are in very good agreement with experimental data.     

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.     

Microscopic models for quantum mechanical calculations of chemical processes in solutions: LD/AMPAC and SCAAS/AMPAC calculations of solvation energies

Our previously developed approaches for integrating quantum mechanical molecular orbital methods with microscopic solvent models are refined and examined. These approaches consider the nonlinear solute–solvent coupling in a self-consistent way by incorporating the potential from the solvent dipoles in the solute Hamiltonian, while considering the polarization of the solvent by the potential from the solute charges. The solvent models used include the simplified Langevin Dipoles (LD) model and the much more expensive surface constrained All Atom Solvent (SCAAS) model, which is combined with a free energy pertubation (FEP) approach. Both methods are effectively integrated with the quantum mechanical AMPAC package and can be easily combined with other quantum mechanical programs. The advantages of the present approaches and their earlier versions over macroscopic reaction field models and supermolecular approaches are considered. A LD/MNDO study of solvated organic ions demonstrates that this model can yield reliable solvation energies, provided the quantum mechanical charges are scaled to have similar magnitudes to those obtained by high level ab initio methods. The incorporation of a field-dependent hydrophobic term in the LD free energy makes the present approach capable of evaluating the free energy of transfer of polar molecules from non polar solvents to aqueous solutions. The reliability of the LD approach is examined not only by evaluating a rather standard set of solvation energies of organic ions and polar molecules, but also by considering the stringent test case of sterically hindered hydrophobic ions. In this case, we compare the LD/MNDO solvation energies to the more rigorous FEP/SCAAS/MNDO solvation energies. Both methods are found to give similar results even in this challenging test case. The FEP/SCAAS/AMPAC method is incorporated into the current version of the program ENZYMIX. This option allows one to study chemical reactions in enzymes and in solutions using the MNDO and AM1 approximations. A special procedure that uses the EVB method as a reference potential for SCF MO calculations should help in improving the reliability of such studies.     

定标粒子理论在动力学介质效应上的应用

应用定标粒子溶液理论,导出了动力学介质效应公式,分析计算了乙酸乙酯皂化反应的介质效应,并与本文实验值相比较,结果,规律一致,数值相近。溶剂体系为水-乙醇和水-二甲亚砜。由于反应体系中包含有极性分子及离子,故在计算相互作用自由焓时,对定标粒子理论作适当补充,包括了色散能、诱导能、取向能和静电作用能。     

HIGH RESOLUTION SCINTILLATION SPECTRA OBTAINED WITH NANOSECOND PULSES OF 3 MeV ELECTRONS*

In studies involving an energy transfer process it is important to have definite knowledge of the excited state of the solute which is responsible for the emitted fluorescence that is measured. When solutions with low solute concentration are irradiated with ionizing particles, the initial excitations are produced in the solvent molecules (excited molecules or ions). The excited solvent molecules and ions can transfer energy to the solute producing excited solute molecules. The measured scintillations are the result of the return of the solute molecules to the ground state.     

Vibrational energy relaxation of a diatomic molecule in a room-temperature ionic liquid

Vibrational energy relaxation (VER) dynamics of a diatomic solute in ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate (EMI(+)PF(6) (-)) are studied via equilibrium and nonequilibrium molecular dynamics simulations. The time scale for VER is found to decrease markedly with the increasing solute dipole moment, consonant with many previous studies in polar solvents. A detailed analysis of nonequilibrium results shows that for a dipolar solute, dissipation of an excess solute vibrational energy occurs almost exclusively via the Lennard-Jones interactions between the solute and solvent, while an oscillatory energy exchange between the two is mainly controlled by their electrostatic interactions. Regardless of the anharmonicity of the solute vibrational potential, VER becomes accelerated as the initial vibrational energy increases. This is attributed primarily to the enhancement in variations of the solvent force on the solute bond, induced by large-amplitude solute vibrations. One interesting finding is that if a time variable scaled with the initial excitation energy is employed, dissipation dynamics of the excess vibrational energy of the dipolar solute tend to show a universal behavior irrespective of its initial vibrational state. Comparison with water and acetonitrile shows that overall characteristics of VER in EMI(+)PF(6) (-) are similar to those in acetonitrile, while relaxation in water is much faster than the two. It is also found that the Landau-Teller theory predictions for VER time scale obtained via equilibrium simulations of the solvent force autocorrelation function are in reasonable agreement with the nonequilibrium results.     

The importance of excluded solvent volume effects in computing hydration free energies

Continuum dielectric methods such as the Born equation have been widely used to compute the electrostatic component of the solvation free energy, DeltaG(solv)(elec), because they do not need to include solvent molecules explicitly and are thus far less costly compared to molecular simulations. All of these methods can be derived from Gauss Law of Maxwell's equations, which yields an analytical solution for the solvation free energy, DeltaG(Born), when the solute is spherical. However, in Maxwell's equations, the solvent is assumed to be a structureless continuum, whereas in reality, the near-solute solvent molecules are highly structured unlike far-solute bulk solvent. Since we have recently reformulated Gauss Law of Maxwell's equations to incorporate the near-solute solvent structure by considering excluded solvent volume effects, we have used it in this work to derive an analytical solution for the hydration free energy of an ion. In contrast to continuum solvent models, which assume that the normalized induced solvent electric dipole density P(n) is constant, P(n) mimics that observed from simulations. The analytical formula for the ionic hydration free energy shows that the Born radius, which has been used as an adjustable parameter to fit experimental hydration free energies, is no longer ill defined but is related to the radius and polarizability of the water molecule, the hydration number, and the first peak position of the solute-solvent radial distribution function. The resulting DeltaG(solv)(elec) values are shown to be close to the respective experimental numbers.     

Simulations of time-dependent fluorescence in nano-confined solvents

The time-dependent fluorescence of a model diatomic molecule with a charge-transfer electronic transition in confined solvents has been simulated. The effect of confining the solvent is examined by comparing results for solutions contained within hydrophobic spherical cavities of varying size (radii of 10-20 angstroms). In previous work [J. Chem. Phys. 118, 6618 (2002)] it was found that the solute position in the cavity critically affects the absorption and fluorescence spectra and their dependence on cavity size. Here we examine the effect of cavity size on the time-dependent fluorescence, a common experimental probe of solvent dynamics. The present results confirm a prediction that motion of the solute in the cavity after excitation can be important in the time-dependent fluorescence. The effects of solvent density are also considered. The results are discussed in the context of interpreting time-dependent fluorescence measurements of confined solvent systems.     

Energy components of aqueous solution: Insight from hybrid QM/MM simulations using a polarizable solvent model

An energy decomposition method is present in statistical Monte Carlo simulations of aqueous solutions of a series of organic solutes, making use of a hybrid quantum mechanical and polarizable molecular mechanical (QM/MM-PIPF) approach. In the hybrid QM/MM-PIPF method, the mutual solute–solvent polarization effect is specifically considered through a coupled iterative procedure that ensures the convergence of solvent induced dipoles and the solute wave function. It should be noted that the method is an approximate approach without specifically considering the electronic correlation effect between solute and solvent electrons, and energetic results have not been verified by free energy calculations. Nevertheless, the energy decomposition analysis provides insight into the details of the molecular polarization effect. Qualitative trends of the energy components from such analyses provide guidance in the understanding of the nature of intermolecular interactions in biomolecular systems, whereas quantitative results on specific terms may be utilized to derive empirical, yet computationally more efficient, force fields. Polarization effects are found to be significant, which contribute 10% to 23% to the total solute–solvent interaction energies. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 :1061–1071, 1997     

The Role of Computational Chemistry in the Experimental Determination of the Dipole Moment of Molecules in Solution

The methods for the experimental determination of electric dipole moment of molecules in solution from measurements of dielectric permittivity and refractive index are traditionally based on the classical Onsager model. In this model the molecular solute is approximated as a simple polarizable point dipole inside a spherical or ellipsoidal cavity of a dielectric medium representing the solvent. However, the inadequacies of the model resulting from the assumption of a simple shape of the cavity, for the evaluation of the cavity field effect, and from the uncertainty of the polarizability of the molecular solute influences the results and hampers the comparison with the electric dipole moments computed from quantum chemical solvation models. In this article we propose a new method for the experimental determination of the electric dipole moment in solution in which information from the Polarizable Continuum Model calculations are used in place of the Onsager model. The new method overcomes the limitations of this latter model regarding both the cavity field effect and the polarizability of the molecular solutes, and thus allows a coherent comparison between experimental and computed dipole moments of solvated molecules. © 2019 Wiley Periodicals, Inc.     

Moment equations for partial filling capillary electrophoresis

Moment equations were developed for partial filling CE systems, in which solute dissolution phenomena by spherical molecular assemblies or intermolecular interactions take place. Because experimental conditions of partial filling CE are divided into five categories on the basis of the magnitude relationship between the migration velocity of solute molecules and that of molecular assemblies or ligand molecules, the moment equations were systematically developed for each case by using the Einstein equation for diffusion and the random walk model. In order to demonstrate the effectiveness of the moment equations, they were applied to the analysis of partial filling CE behavior, which is correlated with dissolution phenomena of small solute molecules into spherical molecular assemblies as specific examples. Simulation results only in the case that the migration velocity of solute molecules is faster than that of molecular assemblies were represented in this paper. Detailed explanations about the derivation procedure of the moment equations and the simulation results in other cases can be found in the Supporting Information. The moment equations are theoretical bases for applying partial filling CE to the study on solute permeation kinetics at the interface of spherical molecular assemblies and on reaction kinetics of intermolecular interactions.     

Solute-solvent interaction in liquids. I. Long-range forces and vibrational solvent shifts

The model of solute— solvent interaction based on dipole-induced-dipole forces (Kirkwood-Bauer-Magat) has been generalized, yielding an expression for the energy as a function of solute position and orientalion within a spherical cavity in a dielectric medium. An analogous relation has been derived for the dispersion energy. Barriers to rotation of the solute molecule and shifts in its vibrational frequency are calculated as functions of cavity radius and eccentricity for the case of dilute solutions of HCl in CCl₄. It is found that the effect of dispersion forces on the vibrational frequency of HCl is two-to-three times more important than the traditional dipole-induced-dipole contribution.     

Toward a new approach for determination of solute's charge distribution to analyze interatomic electrostatic interactions in quantum mechanical/molecular mechanical simulations

We present an alternative approach to determine "density-dependent property"-derived charges for molecules in the condensed phase. In the case of a solution, it is essential to take into consideration the electron polarization of molecules in the active site of this system. The solute and solvent molecules in this site have to be described by a quantum mechanical technique and the others are allowed to be treated by a molecular mechanical method (QM/MM scheme). For calculations based on this scheme, using the forces and interaction energy as density-dependent property our charges from interaction energy and forces (CHIEF) approach can provide the atom-centered charges on the solute atoms. These charges reproduce well the electrostatic potentials around the solvent molecules and present properly the picture of the electron density of the QM subsystem in the solution system. Thus, the CHIEF charges can be considered as the atomic charges under the conditions of the QM/MM simulation, and then enable one to analyze electrostatic interactions between atoms in the QM and MM regions. This approach would give a view of the QM nuclei and electrons different from the conventional methods.     

Rotational-diffusion anomalies in dye solutions from transient-dichroism experiments

Subnanosecond transient-dichroism experiments have been performed to investigate the rotational diffusion of dyes in solution. Dyes and solvents were chosen in a way to obtain information on the influence of size, shape and hydrogen-bonding abilities either of the solute or the solvent molecules. One finds slow orientational relaxation of di-anionic xanthene dyes in alcohols, while oblate cationic dyes rotate faster in spite of their comparable size. The rotational diffusion times for alcohol solutions exceed the theoretical values predicted by the Debye-Einstein model except for prolate molecules. For a solute molecule with internal mobility the rotational diffusion exhibits a partial slip behaviour. It is shown that the deviations from the Debye-Einstein model are restricted to alcohols since for other solvents either with or without strong hydrogen-bonding abilities the experimental values agree with the hydrodynamic model including the stick-boundary condition. Experiments on erythrosine B reveal the influence of size and shape of the attached solvent molecules.     

Modeling absorption spectra of molecules in solution

The presence of solvent tunes many properties of a molecule, such as its ground and excited state geometry, dipole moment, excitation energy, and absorption spectrum. Because the energy of the system will vary depending on the solvent configuration, explicit solute–solvent interactions are key to understanding solution-phase reactivity and spectroscopy, simulating accurate inhomogeneous broadening, and predicting absorption spectra. In this tutorial review, we give an overview of factors to consider when modeling excited states of molecules interacting with explicit solvent. We provide practical guidelines for sampling solute–solvent configurations, choosing a solvent model, performing the excited state electronic structure calculations, and computing spectral lineshapes. We also present our recent results combining the vertical excitation energies computed from an ensemble of solute–solvent configurations with the vibronic spectra obtained from a small number of frozen solvent configurations, resulting in improved simulation of absorption spectra for molecules in solution.     

Evidences of nonideal mixing in poly(ethylene glycol)/organic solvent mixtures by Brillouin scattering

The concentration dependence of the hypersonic properties of solutions of poly(ethylene glycol) of mean molecular mass 600 g/mol (PEG600) in benzene and toluene has been investigated by Brillouin scattering. The two solvents are very similar in structure and chemical properties, but while benzene is nonpolar, toluene possess a modest dipole. In both solvents a high-frequency relaxation process has been observed at high concentrations which has been assigned to conformational rearrangements of the polymeric chains, triggered by reorientation of the side groups. In both cases, the concentration dependence of the adiabatic compressibility deviates significantly from linearity, indicating the existence of nonideal mixing phenomena driven by aggregation processes taking place in the systems. However, there is no temperature dependence for solutions of PEG600 in benzene; on the contrary, the results obtained for solutions of PEG600 in toluene are noticeably dependent on the temperature. The comparison of the experimental data with the results of previous experiments on similar systems allows a general picture for weakly interacting mixtures of hydrogen-bonded systems and organic solvents to be developed. In particular, in the presence of a nonpolar solvent molecule the local structure of the mixture is dominated by solute self-association processes and any resulting solute-solvent correlation is barely induced by excluded volume effects. At high enough dilution the self-aggregation of solute molecules produces a variety of new local topologies that cannot be observed in bulk solute, and as a consequence, the concentration evolution of the system is too rich to be described in terms of a linear combination of a few components over the whole concentration range. The situation seems to be simpler for the polar toluene solvent molecules, where a three-component model seems able to fit the experimental concentration dependence of the hypersonic velocity. This result is interpreted to imply that the interaction between the solvent dipoles and the active sites of the solute produces a relatively stable heterocoordination, while the relevance of self-association is partially reduced.     

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

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

Solvent effect on the absorption spectra of coumarin 120 in water: A combined quantum mechanical and molecular mechanical study

The solvent effect on the absorption spectra of coumarin 120 (C120) in water was studied utilizing the combined quantum mechanical∕molecular mechanical (QM∕MM) method. In molecular dynamics (MD) simulation, a new sampling scheme was introduced to provide enough samples for both solute and solvent molecules to obtain the average physical properties of the molecules in solution. We sampled the structure of the solute and solvent molecules separately. First, we executed a QM∕MM MD simulation, where we sampled the solute molecule in solution. Next, we chose random solute structures from this simulation and performed classical MD simulation for each chosen solute structure with its geometry fixed. This new scheme allowed us to sample the solute molecule quantum mechanically and sample many solvent structures classically. Excitation energy calculations using the selected samples were carried out by the generalized multiconfigurational perturbation theory. We succeeded in constructing the absorption spectra and realizing the red shift of the absorption spectra found in polar solvents. To understand the motion of C120 in water, we carried out principal component analysis and found that the motion of the methyl group made the largest contribution and the motion of the amino group the second largest. The solvent effect on the absorption spectrum was studied by decomposing it in two components: the effect from the distortion of the solute molecule and the field effect from the solvent molecules. The solvent effect from the solvent molecules shows large contribution to the solvent shift of the peak of the absorption spectrum, while the solvent effect from the solute molecule shows no contribution. The solvent effect from the solute molecule mainly contributes to the broadening of the absorption spectrum. In the solvent effect, the variation in C-C bond length has the largest contribution on the absorption spectrum from the solute molecule. For the solvent effect on the absorption spectrum from the solvent molecules, the solvent structure around the amino group of C120 plays the key role.     

Effect of attractions on the structure of polymer solutions confined between surfaces: a density functional approach

A density functional theory is presented to study the effect of attractions on the structure of polymer solutions confined between surfaces. The polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres, both having Yukawa-type attractions and the mixture being confined between attractive Yukawa-type surfaces. The present theory treats the ideal gas free energy functional exactly and uses weighted density approximation for the hard chain and hard sphere contributions to the excess free energy functional. The attractive interactions are calculated using the direct correlation function obtained from the polymer reference interaction site model theory along with the mean spherical approximation closure. The theoretical predictions on the density profiles of the polymer and the solvent molecules are found to agree quite well with the Monte Carlo simulation results for varying densities, chain lengths, wall separations, and different sets of interaction potentials.