Continuum Medium Model of Non-Equilibrium Solvation and Solvent Effect for Ultra-Fast Process
QA Zhu,KX Fu,XY Li
DOI: https://doi.org/10.3321/j.issn:0251-0790.2006.02.014
2006-01-01
Abstract:In the 1950s of the 20th century, the non-equilibrium solvation theory for ultra-fast processes such as electron transfer and light absorption/emission was paid particular attention to. A number of scientists made efforts to study this area and various models, which give reasonable qualitative descriptions for solvent reorganization energy in electron transfer and spectral shift in solution, were developed within the framework of continuous medium theory. However, in a series of publications by the authors ( see for example J. Comput. Chem. 2004, 25: 500; J. Comput. Chem. 2004, 25: 835; J. Comput. Chem. 2005, 26: 399; Chinese Sci. Bull. 2003, 48: 965; J. Mole. Struct. -Theochem 2005, 715: 157), it was clarified that the expression of the nonequilibrium electrostatic free energy, which is at the dominant position of non-equilibrium solvation and serves as the basis of various models, was incorrectly formulated. In this work, the authors argue that reversible charging work integration delta W = integral(V)Phi delta rho dV was inappropriately applied to an irreversible path linking the equilibrium state and the non-equilibrium one in the past. Because the step from the equilibrium state to the non-equilibrium one is factually thermodynamically irreversible, the conventional expression for non-equilibrium free energy, G(2)(non) (M) = (1/2) integral(V)(rho(2)Phi(non)(2) + rho(2)Phi(eq)(1) - rho(1)Phi(non)(2)) dV that was deduced in different ways, unreasonable. Here the authors derive the non-equilibrium free energy to a quite different form of G(2)(non) = (1/2) integral(V)rho(2)Phi(non)(2) dV according to Jackson integral formula, dG = (1/2) integral(V)(Phi delta rho + rho delta Phi)dV. Such a difference throws doubts to the v models including the famous Marcus two-sphere model for solvent reorganization energy of electron transfer and the Lippert-Mataga equation for spectral shift. By introducing the concept of "spring energy" arising from medium polarizations, the energy constitution of the non-equilibrium state is highlighted. For a solute-solvent system, the authors separate the total electrostatic energy into different components: the self-energies of solute charge and polarized charge, the interaction energy between them and the "spring energy" of the solvent polarization. With detailed reasoning and derivation, our formula for non-equilibrium free energy can be reached through different ways. Based on the new expression of non-equilibrium free energy, the generalized form for solvent reorganization energy, lambda(BV) = (1/4) integral(V)Delta rho(Delta phi(op) - Delta phi) dV, has been attained. A new two-sphere model for solvent reorganization energy is proved to have the form of lambda(av) = (1/2)Delta q(2)(1/epsilon(op) - 1/epsilon(s)) (1/(2r(D)) + 1/(2r(A)) 1/d). Compared with Marcus' expression, this new formula estimates the solvent reorganization,,only one half of the latter. This difference provides a pretty explanation for why Marcus' theory often energy over estimated the solvent reorganization energy by a factor about two in the past.With the single-sphere model and point dipole approximation, the authors argue that the total spectral shift should look like Delta h nu(total) = ( 1/2) R-slow (mu 1 - & m;(m))(2)(,) and this is also one half of the Lippert-Mataga result. The novel expressions for the spectral shifts for individual absorption and emission have also been given. Finally, a numerical algorism for the solution of Poisson equation is presented and the total non-equilibrium solvation energy is deduced to a quite different and much more compact form as Delta F-2(non) = - (Psi(gas)(2) vertical bar H-0 vertical bar Psi(gas)(2)) + (1/2) Sigma(i)Sigma(j)[q(2,fast) (i) + q(1,slow) (i) ]Z(j)/(vertical bar r(i) - R-j vertical bar) when compared with the most recently developed expression by other authors. As an application, the numerical algorism incorporated with COSMO was applied to a model system, and the solvent reorganization energy is found in excellent agreement with the experimental fitting, while the conventional theories always estimate twice this quantity.