Solvate systems

Be sure to remind students that these frequencies are gas phase data and arc thus not the same as the more-faniiliar solution spectra (we will treat solvated systems in Chapter 9). Even so, such gas phase calculations make excellent discovery-based exercises, For example, students may be asked to explain the substituent effects observed tising basic chemistry knowledge. 

As the plot of AE indicates, the energy difference between the two forms decreases in more polar solvents, and becomes nearly zero in acetonitrile. The left plot illustrates the fact that the IPCM model (at the B3LYP/6-31+G(d) level of theory) does a much better job of reproducing the observed solvent effect than the two Onsager SCRF models. In contrast, the Onsager model at the MP2 level treats the solvated systems more accurately than it does the gas phase system, leading to a poorer value for the solvent effect. ... 

The complex of protein, crystallographic water and the counter ions are treated as a fully solvated system. Two key developments were made to treat the system in a more realistic manner during molecular dynamics (i) water molecules were placed as a spherical shell with a radius 34 A from the center of the protein. The outer boundary of the spherical shell was defined by means of an artificial wall with a potential of the type W defined as ... 

This potential was developed to ensure that the molecules inside the sphere never escape and maintain a fully solvated system during molecular dynamics. Here, es, Rs, ew and Rw are the van der Waals constants for the solvent and the wall and rj is the distance between the molecule i and the center of the water sphere, Ro is the radius of the sphere. The quantities A, B and Rb are determined by imposing the condition that W and dW/dr, vanish at r, = Ro. The restraining potential W is set to zero for r, < R0. The van der Waals parameters Es, ew, Rs and Rw can also be specifically defined for different solvents. The constants Awaii and Cwan are computed using a well depth of es = ew = 0.1 kcal and the radius of Rs = Rw = 1.25 A. For the other set of simulations, especially for the hydride ion transfer, we applied periodic boundary conditions by using a spherical boundary shell of 10.0 A of TIP3P40 water to cover the edges of the protein. 

First, we remove the solvent and consider only the system of adsorbent and ligand molecules. We make this simplification not because solvent effects are unimportant or negligible. On the contrary, they are very important and sometimes can dominate the behavior of the systems. We do so because the development of the theory of cooperativity of a binding system in a solvent is extremely complex. One could quickly lose insight into the molecular mechanism of cooperativity simply because of notational complexity. On the other hand, as we shall demonstrate in subsequent chapters, one can study most of the aspects of the theory of cooperativity in unsolvated systems. What makes this study so useful, in spite of its irrelevance to real systems, is that the basic formalism is unchanged by introducing the solvent. The theoretical results obtained for the unsolvated system can be used almost unchanged, except for reinterpretation of the various parameters. We shall discuss solvated systems in Chapter 9. 

A high-level quantum chemical exploration of the Horner-Wandsworth-Emmons reaction has indicated that ring closure of the P—O bond (which favors formation of -product) is rate determining in the gas phase and that the C—C bond-forming addition step is rate determining in most solvated systems several effects that could account for the E/Z selectivities observed have been identified. 

Because of their unique characteristics, supercritical fluids have received a great deal of attention in a number of important scientific fields (1-14). Several reasons are given for choosing a supercritical fluid over another solvating system, but choice is governed generally by 1) the unique solvation and favorable mass transport properties (5) and 2) the ease with which the chemical potential can be varied simply by adjustment of the system pressure and/or temperature (13). 

For molecular systems in the vacuum, exact analytical derivatives of the total energy with respect to the nuclear coordinates are available [22] and lead to very efficient local optimization methods [23], The situation is more involved for solvated systems modelled within the implicit solvent framework. The total energy indeed contains reaction field contributions of the form ER(p,p ), which are not calculated analytically, but are replaced by numerical approximations Efp(p,p ), as described in Section 1.2.5. We assume from now on that both the interface Y and the charge distributions p and p depend on n real parameters (A, , A ). In the geometry optimization problem, the A, are the cartesian coordinates of the nuclei. There are several nonequivalent ways to construct approximations of the derivatives of the reaction field energy with respect to the parameters (A1 , A ) ... 

Density Functional Theory does not require specific modifications, in relation to the solvation terms [9], with respect to the Hartree-Fock formalism presented in the previous section. DFT also absorbs all the properties of the HF approach concerning the analytical derivatives of the free energy functional (see also the contribution by Cossi and Rega), and as a matter of fact continuum solvation methods coupled to DFT are becoming the routine approach for studies of solvated systems. 

This equation shows that vertical excitations in solvated systems are obtained as a sum of two terms, the difference in the excited and ground state energies in the presence of a frozen ground state solvent and a relaxation term determined by the mutual polarization... 

In this contribution we have presented some specific aspects of the quantum mechanical modelling of electronic transitions in solvated systems. In particular, attention has been focused on the ASC continuum models as in the last years they have become the most popular approach to include solvent effects in QM studies of absorption and emission phenomena. The main issues concerning these kinds of calculations, namely nonequilibrium effects and state-specific versus linear response formulations, have been presented and discussed within the most recent developments of modern continuum models. 

These effective properties represent the main result of the theoretical formulation of NLO phenomena for solvated systems [2-6] as they describe the response of the solute in terms of the macroscopic field in the surrounding medium and thus they may be directly related to the macroscopic properties determined from experiments. 

The present contribution considers general electronic states of solvated molecules and is not limited to closed shell molecular compounds. For closed shell molecular systems, methods utilizing closed shell coupled-cluster electronic structure and closed shell density density functional theory for the electronic structure of the solvated system have appeared in the literature [54-67],... 

The second chapter ends with two overviews by Stephens Devlin and by Hug on the theoretical and the physical aspects of two vibrational optical activity spectroscopies (VCD and VROA, respectively). In both overviews the emphasis is more on their basic formalism and the gas-phase quantum chemical calculations than on the analysis of solvent effects. For these spectroscopies, in fact, both the formulation of continuum solvation models and their applications to realistic solvated systems are still in their infancy. 

The radicals that have been investigated in the gas phase are listed by Jacox (1984, 1988, 1990), Miller (Foster and Miller 1989), and Huber (Huber and Herzberg 1979). Solvated systems have also been studied and have been reviewed by Heaven (Heaven 1992) and Lester (Giancarlo et al. 1994). The cluster chemistry studies, as such, are discussed below for benzyl (C6H5CH2 ) and methyl radicals. Work in progress on small radicals, carbenes, and nitrenes will be briefly mentioned at the end of the section. 

In order to characterise ion solvation processes, gas phase studies can be performed providing detailed information about individual interactions. These studies can explore changes in some properties between the complexes in the gas phase and the solvated systems in the liquid phase. Theoretical methods can thus provide valuable complementary information not accessible to experimental approaches, both in the characterization of the complexes and in the specific mechanism of the relevant interactions. 

By contrast, the description given by a continuum description does not require any knowledge of the solvent configuration around the solute as a structureless continuum dielectric is introduced instead. The response of such a dielectric to the presence of the solute is determined by its macroscopic properties (namely the dielectric constant and the refractive index) and thus it will be implicitly averaged. Contrary to what happens in a QM/MM approach, here a single calculation on a given solute contained within the continuum dielectric will be sufficient to get the correct picture of the solvated system. 

In the previous sections we have briefly summarized the basic theory of QM/MM and QM/continuum methods showing their differences and similarities, now we can move on to describe their applications to the calculation of molecular response properties and the related spectroscopies for a generic solvated system. 

The QM theory of chemical shielding was originally developed many years ago [22,23], but only later have ab initio methods and density functional theories (DFT) been reliably used for the prediction of NMR properties of isolated molecular systems, and finally of solvated systems. The latter step has been achieved by extending the gas-phase theoretical methods to continuum solvation models (see Ref. [11] for a sufficiently updated list of papers). 

An Application to Solvated Systems N Nuclear Shieldings of Diazines... 

In this section we shall present and compare different computational strategies one can adopt to simulate the effect of the environment on spectroscopic properties of solvated systems. In particular, as a representative example, we shall summarize the results of two studies [30,31] we have published in the last years on the environment effects on the nitrogen nuclear shieldings of a specific class of molecular systems containing sp2-type nitrogens three diazines, also known as pyridazine (1,2-diazine), pyrimidine (1,3-diazine) and pyrazine (1,4-diazine). 

All calculations on the isolated and the PCM solvated systems have been performed using the Gaussian code [39] while the DPM calculations have been performed using the development version of the Dalton Quantum Chemistry Program [40],... 

In Figure 7-11 the charge densities related to the KS eigenstates near the Fermi level of the hydrated C6o are shown. The HOMO of the solvated system... 

The examples of applications given in the latter part of this chapter will show that even at the present state of the art and technology, many solvated ions could be treated with sufficient quality to obtain reliable results not only for structural details and species distributions, but also for the aforementioned ultrafast dynamical processes determining the chemical behaviour of such solvates. On the other hand, the latest improvements of the simulation methodology have opened a straightforward access to the treatment of other arbitrary solvated systems as computational capabilities increase. Therefore, simulation methods are not only becoming a valuable research field of their own, but also an essential supplement - if not prerequisite -for the interpretation of experimental investigations of solvates. 

The earliest applications of quantum chemistry were targeted toward molecules in the gas phase. This was in part due to an interest in such an environment which allows researchers to focus on the intrinsic properties of the molecule of interest. But also, since much of practical chemistry takes place in solution of some sort, this environment presents an important avenue of inquiry as well. However, such solvated systems were typically out of reach of quantum calculations. For one thing, the inclusion of a number of solvent molecules into the calculations would commonly take the system beyond the capabilities of computers at the time. Secondly, solvent... 

---