NH#tex2html_wrap_inline2748# in ammonia and OH#tex2html_wrap_inline2750# in water next up previous
Next: Conclusions Up: Comparison with HO and Previous: NH in ammonia and

NH tex2html_wrap_inline1898 in ammonia and OH tex2html_wrap_inline1970 in water

By far, the most striking similarity between NH tex2html_wrap_inline1898 in ammonia and OH tex2html_wrap_inline1970 in water is the unusually high coordination number around the heavy atom. Based on the number of lone pairs around the heavy atom, a naive estimate of the number of hydrogen bonds would be 3 for OH tex2html_wrap_inline1970 in water and 2 for NH tex2html_wrap_inline1898 in ammonia. Indeed, gas phase quantum chemical calculations of OH tex2html_wrap_inline1970 with three water molecules confirm that threefold coordination of the oxygen atom is the lowest energy structure [73]. However, ab initio MD simulations of OH tex2html_wrap_inline1970 in water at 300 K [8, 9] showed that, in the liquid, structures in which four or five hydrogen bonds surround the oxygen are possible, and, in addition, the OH tex2html_wrap_inline1970 hydrogen can also form a weak hydrogen bond with nearby water molecules. The overall coordination number predicted from these simulations is approximately 5.3, which is in excellent agreement with the recently measured value from dieletric relaxation experiments of 5.5 tex2html_wrap_inline2766 0.5 [74]. Recent quantum chemical calculations of OH tex2html_wrap_inline1970 with four water molecules [73] show that threefold coordination of the OH tex2html_wrap_inline1970 oxygen with the fourth water molecule in the second solvation shell is the lowest energy structure but that fourfold coordination of the oxygen is only about 1.2 kcal/mol higher in energy and represents a local minimum on the potential energy surface. In the case of NH tex2html_wrap_inline1898 in ammonia, the coordination is also considerably higher than expected based on the number of lone pairs. In Sec. 5.1, it was seen that the coordination number of NH tex2html_wrap_inline1898 in ammonia is between 7 and 8. However, unlike the OH tex2html_wrap_inline1970 case, not all of this is due to hydrogen bonding (cf. Fig. 9) and the discussion of Fig. 5.1).

The ELF used in Sec. 4.3 and 5.4 can be used to provide an explanation of why the coordination of the anionic defects is so large. Figure 5.4 shows the ELF for NH tex2html_wrap_inline1898 in the gas phase and in some typical complexes. For comparison, the ELF for OH tex2html_wrap_inline1970 in the gas phase and in some typical complexes is shown in Fig. 6.2.

   figure701

Figure 19: The electron localization function (ELF) (cf. Eq. (9)) corresponding to OH tex2html_wrap_inline1970 in the gas phase (a), with a threefold coordinated first solvation shell (H tex2html_wrap_inline2002 O tex2html_wrap_inline2004 )(b), and with a fourfold coordinated first solvation shell (H tex2html_wrap_inline2006 O tex2html_wrap_inline2008 )(c). Because tex2html_wrap_inline1912 corresponds to perfect localization, the contour cutoff level is chosen to be 0.93 in the figure.

Plots such as these were also introduced in Ref. [76] in the context of fully quantum mechanical studies of OH tex2html_wrap_inline1970 solvation in water. The figure shows the dramatic delocalization of the three lone pairs of OH tex2html_wrap_inline1970 leading to the formation of a torus of nearly equal probability. This negatively charged ring allows hydrogens to form hydrogen bonds to the hydroxyl oxygen without needing to adhere to specific ``basins'' of negative attraction, thereby allowing for a higher coordination than would be expected based on the actual number of lone pairs. Similarly, the ELF for NH tex2html_wrap_inline1898 shows that half of such a ring is formed by the two lone pairs. This half ring also allows hydrogen bonding to the NH tex2html_wrap_inline1898 nitrogen to occur with no adherence to basins of attraction and leads to the structures observed in Fig. 5.1.

Several key differences exist between these ions as well. One clear difference is that the NH tex2html_wrap_inline1898 hydrogens do not form hydrogen bonds to neighboring ammonia molecules, while the OH tex2html_wrap_inline1970 hydrogen is capable of forming weak hydrogen bonds to solvent molecules [8, 9, 76]. The formation of such hydrogen bonds has been implicated as a crucial element in the structural diffusion of the OH tex2html_wrap_inline1970 defect in water [76]. The inability of NH tex2html_wrap_inline1898 hydrogens to form such weak hydrogen bonds, combined with the coordination number mismatch alluded to in the previous section, is the most likely explanation for the lack of proton transfer mediated structural diffusion of this anionic defect in ammonia.


next up previous
Next: Conclusions Up: Comparison with HO and Previous: NH in ammonia and

Mark Tuckerman
Mon Mar 26 01:25:53 EST 2001