Water molecules method

The occurrence of excited-state proton transfer during the lifetime of the excited state depends on the relative rates of de-excitation and proton transfer. The general equations will be presented first, but only for the most extensively studied case where the excited-state process is proton ejection (pK < pK) the proton donor is thus an acid, AH, and the proton acceptor is a water molecule. Methods for the determination of pK are then described and finally, the various cases of pH dependence of the absorption and fluorescence spectra are examined. 

The anhydrous chloride is prepared by standard methods. It is readily soluble in water to give a blue-green solution from which the blue hydrated salt CuClj. 2H2O can be crystallised here, two water molecules replace two of the planar chlorine ligands in the structure given above. Addition of dilute hydrochloric acid to copper(II) hydroxide or carbonate also gives a blue-green solution of the chloride CuClj but addition of concentrated hydrochloric acid (or any source of chloride ion) produces a yellow solution due to formation of chloro-copper(ll) complexes (see below). 

Fig. 2. The time evolution of the total energy of four water molecules (potential-energy details are given in [48]) as propagated by the symplectic Verlet method (solid) and the nonsymplectic fourth-order Runge-Kutta method (dashed pattern) for Newtonian dynamics at two timestep values.

In Table 1 the CPU time required by the two methods (LFV and SISM) for 1000 MD integration steps computed on an HP 735 workstation are compared for the same model system, a box of 50 water molecules, respectively. The computation cost per integration step is approximately the same for both methods so that th< syieed up of the SISM over the LFV algorithm is deter-... 

The correct treatment of boundaries and boundary effects is crucial to simulation methods because it enables macroscopic properties to be calculated from simulations using relatively small numbers of particles. The importance of boundary effects can be illustrated by considering the following simple example. Suppose we have a cube of volume 1 litre which is filled with water at room temperature. The cube contains approximately 3.3 X 10 molecules. Interactions with the walls can extend up to 10 molecular diameters into the fluid. The diameter of the water molecule is approximately 2.8 A and so the number of water molecules that are interacting with the boundary is about 2 x 10. So only about one in 1.5 million water molecules is influenced by interactions with the walls of the container. The number of particles in a Monte Carlo or molecular dynamics simulation is far fewer than 10 -10 and is frequently less than 1000. In a system of 1000 water molecules most, if not all of them, would be within the influence of the walls of the boundary. Clecirly, a simulation of 1000 water molecules in a vessel would not be an appropriate way to derive bulk properties. The alternative is to dispense with the container altogether. Now, approximately three-quarters of the molecules would be at the surface of the sample rather than being in the bulk. Such a situation would be relevcUit to studies of liquid drops, but not to studies of bulk phenomena. 

The use of QM-MD as opposed to QM-MM minimization techniques is computationally intensive and thus precluded the use of an ab initio or density functional method for the quantum region. This study was performed with an AMi Hamiltonian, and the first step of the dephosphorylation reaction was studied (see Fig. 4). Because of the important role that phosphorus has in biological systems [62], phosphatase reactions have been studied extensively [63]. From experimental data it is believed that Cys-i2 and Asp-i29 residues are involved in the first step of the dephosphorylation reaction of BPTP [64,65]. Alaliambra et al. [30] included the side chains of the phosphorylated tyrosine, Cys-i2, and Asp-i 29 in the quantum region, with link atoms used at the quantum/classical boundaries. In this study the protein was not truncated and was surrounded with a 24 A radius sphere of water molecules. Stochastic boundary methods were applied [66]. 

However, theories that are based on a basis set expansion do have a serious limitation with respect to the number of electrons. Even if one considers the rapid development of computer technology, it will be virtually impossible to treat by the MO method a small system of a size typical of classical molecular simulation, say 1000 water molecules. A logical solution to such a problem would be to employ a hybrid approach in which a chemical species of interest is handled by quantum chemistry while the solvent is treated classically. 

In the studies by Skipper et al. the number of water layers (and thus molecules) was fixed on the basis of experimental evidence consequently, the stable states or degrees of swelhng were presumed. Quite differently, Karaborni et al. [44] determined, by means of a combination of GCMC and MD, the number of water molecules directly from a series of simulations in which the distance between montmorillonite planes was varied systematically. They observed that swelling proceeded from the dry state through the formation of one, three, and then five layers of water. This is very different from the usually beheved hydration sequence from one layer to two, then to three layers, and so on, which has been intrinsically assumed by Skipper and coworkers. The authors conclude that the complex swelling behavior accounts for many of the experimental facts. This work demonstrates impressively the power of the grand canonical simulation method. 

We ll look at a simple example of the latter method here, converting the water molecule structure saved in Brookhaven Protein Data Bank (PDB) format. 

The above method is unsatisfactory when hydration takes place at two alternative sites in the molecule, although one hydrate is usually present in only a very small proportion, at equilibrium. Which oxo compound is preferentially formed in such a case depends on the rates of oxidation at the different sites and on the rate of isomerization of the water molecule from one position to the other, hence this method does not indicate which is the thermodynamically more stable hydrate. 

---