Water wave functions

Drukker, K., Hammes-Schiffer, S. An analytical derivation of MC-SCF vibrational wave functions for the quantum dynamical simulation of multiple proton transfer reactions Initial application to protonated water chains. J. Chem. Phys. 107 (1997) 363-374. 

In most solids, the sound speed is an increasing function of pressure, and it is that property that causes a compression wave to steepen into a shock. The situation is similar to a shallow water wave, whose velocity increases with depth. As the wave approaches shore, a small wavelet on the trailing, deeper part of the wave moves faster, and eventually overtakes similar disturbances on the front part of the wave. Eventually, the water wave becomes gravitationally unstable and overturns. 

Electron Nuclear Dynamics (48) departs from a variational form where the state vector is both explicitly and implicitly time-dependent. A coherent state formulation for electron and nuclear motion is given and the relevant parameters are determined as functions of time from the Euler equations that define the stationary point of the functional. Yngve and his group have currently implemented the method for a determinantal electronic wave function and products of wave packets for the nuclei in the limit of zero width, a "classical" limit. Results are coming forth protons on methane (49), diatoms in laser fields (50), protons on water (51), and charge transfer (52) between oxygen and protons. 

In a recent paper. Mo and Gao [5] used a sophisticated computational method [block-localized wave function energy decomposition (BLW-ED)] to decompose the total interaction energy between two prototypical ionic systems, acetate and meth-ylammonium ions, and water into permanent electrostatic (including Pauli exclusion), electronic polarization and charge-transfer contributions. Furthermore, the use of quantum mechanics also enabled them to account for the charge flow between the species involved in the interaction. Their calculations (Table 12.2) demonstrated that the permanent electrostatic interaction energy dominates solute-solvent interactions, as expected in the presence of ion species (76.1 and 84.6% for acetate and methylammonium ions, respectively) and showed the active involvement of solvent molecules in the interaction, even with a small but evident flow of electrons (Eig. 12.3). Evidently, by changing the solvent, different results could be obtained. 

Overall, the accuracy reached by the hybrid functionals tested is quite pleasing, in particular if we keep in mind the favorable scaling of computational demands and the obviously much lower basis set requirements. The binding energy of the water molecule computed with the B3LYP and the mPWIPW functional miss the lower experimental error bar by only 0.1 kcal/mol and approach quite closely the best available conventional wave function based data. 

Secondly, correlations in the initial state can lead to experimental orbital momentum densities significantly different from the calculated Hartree-Fock ones. Figure 3 shows such a case for the outermost orbital of water, showing how electron-electron correlations enhance the density at low momentum. Since low momentum components correspond in the main to large r components in coordinate space, the importance of correlations to the chemically interesting long range part of the wave function is evident. 

The possible wave functions for the molecular orbitals for molecules are those constructed from the irreducible representations of the groups giving the symmetry of the molecule. These are readily found in the character table for the appropriate point group. For water, which has the point group C2 , the character table (see Table 5.4) shows that only A1 A2, B1 and B2 representations occur for a molecule having C2 symmetry. 

How well can continuum solvation models distinguish changes in one or another of these solvent properties This is illustrated in Table 2, which compares solvation energies for three representative solutes in eight test solvents. Three of the test solvents are those shown in Table 1, one is water, and the other four were selected to provide useful comparisons on the basis of their solvent descriptors, which are shown in Table 3. Notice that all four solvents in Table 3 have no acidity, which makes them more suitable, in this respect, than 1-octanol or chloroform for modeling biomembranes. Table 2 shows that the SM5.2R model, with gas-phase geometries and semiempirical molecular orbital theory for the wave function, does very well indeed in reproducing all the trends in the data. 

---