Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

Hydration of the Bisulfate Ion: Atmospheric Implications

Using molecular dynamics configurational sampling combined with ab initio energy calculations, we determined the low energy isomers of the bisulfate hydrates. We calculated the CCSD(T) complete basis set (CBS) binding electronic and Gibbs free energies for 53 low energy isomers of HSO4–(H2O)n=1–6 and derived the thermodynamics of adding waters sequentially to the bisulfate ion and its hydrates. Comparing the HSO4–/H2O system to the neutral H2SO4/H2O cluster, water binds more strongly to the anion than it does to the neutral molecules. The difference in the binding thermodynamics of HSO4–/H2O and H2SO4/H2O systems decreases with increasing number of waters. The thermodynamics for the formation of HSO4–(H2O)n=1–5 is favorable at 298.15 K, and that of HSO4–(H2O)n=1–6 is favorable for T < 273.15 K. The HSO4– ion is almost always hydrated at temperatures and relative humidity values encountered in the troposphere. Because the bisulfate ion binds more strongly to sulfuric acid than it does to water, it is expected to play a role in ion-induced nucleation by forming a strong complex with sulfuric acid and water, thus facilitating the formation of a critical nucleus.

Characterization of HgCl2 Tridentate Amine Complexes by X-ray Crystallography, NMR and ESI-MS

Two new HgCl2 complexes of tridentate nitrogen ligands were characterized by X-ray crystallography, proton NMR spectroscopy and ESI-MS. The five-coordinate complex [Hg(BMPA)Cl-2] (1) (BMPA = bis(2-pyridylmethyl)amine) crystallized from acetonitrile/m-xylene by slow evaporation in the monoclinic space group P2(1)/n with a = 8.3896(8) , b = 12.8020(13) , c = 13.3526(13) , alpha = 90A degrees, beta A = 90.480(2)A degrees, gamma A = 90A degrees and z = 4. The square pyramidal structure (tau = 0.009) has approximate C (s) symmetry. Despite comparable Hg-N bond lengths in 1, inversion of the central nitrogen was rapid on the chemical shift time scale in dilute solution except at very low temperatures. The related complex [Hg(BEPA)Cl-2] (2) (BEPA = bis(2-{pyrid-2-yl}ethyl)amine) crystallized from acetonitrile/ethyl acetate/hexanes by slow diffusion in the orthorhombic space group Pnma with a = 13.424(3) , b = 14.854(3) , c = 8.118(2) , alpha = 90A degrees, beta A = 90A degrees, gamma A = 90A degrees and z = 4. The mixed geometry structure (tau = 0.56) also has crystallographic mirror symmetry as well as C (s) point group symmetry. In dilute acetonitrile solution, 1 was stable while 2 slowly converted to a more thermodynamically stable complex.