Polarization functions are added in two steps to a split-valence extended gaussian basis set: d-type gaussians on the first row atoms C. N, O and F and p-type gaussians on hydrogen. The same d-exponent of 0.8 is found to be satisfactory for these four atoms and the hydrogen p-exponent of 1.1 is adequate in their hydrides. The energy lowering due to d functions is found to depend on the local symmetry around the heavy atom. For the particular basis used, the energy lowerings due to d functions for various environments around the heavy atom are tabulated. These bases are then applied to a set of molecules containing up to two heavy atoms to obtain their LCAO-MO-SCF energies. The mean absolute deviation between theory and experiment (where available) for heats of hydrogenation of closed shell species with two non-hydrogen atoms is 4 kcal/mole for the basis set with full polarization. Estimates of hydrogenation energy errors at the Hartree-Fock limit, based on available calculations, are given.

The influence of polarization functions on molecular orbital hydrogenation energies

We present an approach for deriving net atomic charges from ab initio quantum mechanical calculations using a least squares fit of the quantum mechanically calculated electrostatic potential to that of the partial charge model. Our computational approach is similar to those presented by Momany [J. Phys. Chem., 82, 592 (1978)], Smit, Derissen, and van Duijneveldt [Mol. Phys., 37, 521 (1979)], and Cox and Williams [J. Comput. Chem., 2, 304 (1981)], but differs in the approach to choosing the positions for evaluating the potential. In this article, we present applications to the molecules H2O, CH3OH, (CH3)2O, H2CO, NH3, (CH3O)2PO, deoxyribose, ribose, adenine, 9-CH3 adenine, thymine, 1-CH3 thymine, guanine, 9-CH3 guanine, cytosine, 1-CH3 cytosine, uracil, and 1-CH3 uracil. We also address the question of inclusion of “lone pairs,” their location and charge.

An approach to computing electrostatic charges for molecules

IA Basic Concepts.- 1 The Importance of Hydrogen Bonds.- 1.1 Historical Perspective.- 1.2 The Importance of Hydrogen Bonds in Biological Structure and Function.- 1.3 The Role of the Water Molecules.- 1.4 Significance of Small Molecule Crystal Structural Studies.- 1.5 The Structural Approach.- 2 Definitions and Concepts.- 2.1 Definition of the Hydrogen Bond - Strong and Weak Bonds.- 2.2 Hydrogen-Bond Configurations: Two- and Three-Center Hydrogen Bonds Bifurcated and Tandem Bonds.- 2.3 Hydrogen Bonds Are Very Different from Covalent Bonds.- 2.4 The van der Waals Radii Cut-Off Criterion Is Not Useful.- 2.5 The Concept of the Hydrogen-Bond Structure.- 2.6 The Importance of ? and ? Cooperativity.- 2.7 Homo-, Anti- and Heterodromic Patterns.- 2.8 Hydrogen Bond Flip-Flop Disorder: Conformational and Configurational.- 2.9 Proton-Deficient Hydrogen Bonds.- 2.10 The Excluded Region.- 2.11 The Hydrophobic Effect.- 3 Experimental Studies of Hydrogen Bonding.- 3.1 Infrared Spectroscopy and Gas Electron Diffraction.- 3.2 X-Ray and Neutron Crystal Structure Analysis.- 3.3 Treatment of Hydrogen Atoms in Neutron Diffraction Studies.- 3.4 Charge Density and Hydrogen-Bond Energies.- 3.5 Neutron Powder Diffraction.- 3.6 Solid State NMR Spectroscopy.- 4 Theoretical Calculations of Hydrogen-Bond Geometries.- 4.1 Calculating Hydrogen-Bond Geometries.- 4.2 Ab-Initio Molecular Orbital Methods.- 4.3 Application to Hydrogen-Bonded Complexes.- 4.4 Semi-Empirical Molecular Orbital Methods.- 4.5 Empirical Force Field or Molecular Mechanics Methods.- 5 Effect of Hydrogen Bonding on Molecular Structure.- IB Hydrogen-Bond Geometry.- 6 The Importance of Small Molecule Structural Studies.- 6.1 Problems Associated with the Hydrogen-Bond Geometry.- 6.2 The Hydrogen Bond Can Be Described Statistically.- 6.3 The Problems of Measuring Hydrogen-Bond Lengths and Angles in Small Molecule Crystal Structures.- 7 Metrical Aspects of Two-Center Hydrogen Bonds.- 7.1 The Metrical Properties of O-H *** O Hydrogen Bonds.- 7.1.1 Very Strong and Strong OH *** O Hydrogen Bonds Occur with Oxyanions, Acid Salts, Acid Hydrates, and Carboxylic Acids.- 7.1.2 OH *** O Hydrogen Bonds in the Ices and High Hydrates.- 7.1.3 Carbohydrates Provide the Best Data for OH ... O Hydrogen Bonds: Evidence for the Cooperative Effect.- 7.2 N-H *** O Hydrogen Bonds.- 7.3 N-H *** N Hydrogen Bonds.- 7.4 O-H *** N Hydrogen Bonds.- 7.5 Sequences in Lengths of Two-Center Hydrogen Bonds.- 7.6 H/D Isotope Effect.- 8 Metrical Aspects of Three- and Four-Center Hydrogen Bonds.- 8.1 Three-Center Hydrogen Bonds.- 8.2 Four-Center Hydrogen Bonds.- 9 Intramolecular Hydrogen Bonds.- 10 Weak Hydrogen-Bonding Interactions Formed by C-H Groups as Donors and Aromatic Rings as Acceptors.- 11 Halides and Halogen Atoms as Hydrogen-Bond Acceptors.- 12 Hydrogen-Bond Acceptor Geometries.- II Hydrogen Bonding in Small Biological Molecules.- 13 Hydrogen Bonding in Carbohydrates.- 13.1 Sugar Alcohols (Alditols) as Model Cooperative Hydrogen-Bonded Structures.- 13.2 Influence of Hydrogen Bonding on Configuration and Conformation in Cyclic Monosaccharides.- 13.3 Rules to Describe Hydrogen-Bonding Patterns in Monosaccharides.- 13.4 The Water Molecules Link Hydrogen-Bond Chains into Nets in the Hydrated Monosaccharide Crystal Structures.- 13.5 The Disaccharide Crystal Structures Provide an Important Source of Data About Hydrogen-Bonding Patterns in Polysaccharides.- 13.6 Hydrogen Bonding in the Tri- and Tetrasaccharides Is More Complex and Less Well Defined.- 13.7 The Hydrogen Bonding in Polysaccharide Fiber Structures Is Poorly Defined.- 14 Hydrogen Bonding in Amino Acids and Peptides: Predominance of Zwitterions.- 15 Purines and Pyrimidines.- 15.1 Bases Are Planar and Each Contains Several Different Hydrogen-Bonding Donor and Acceptor Groups.- 15.2 Many Tautomeric Forms Are Feasible But Not Observed.- 15.3 ?-Bond Cooperativity Enhances Hydrogen-Bonding Forces.- 15.4 General, Non-Base-Pairing Hydrogen Bonds.- 16 Base Pairing in the Purine and Pyrimidine Crystal Structures.- 16.1 Base-Pair Configurations with Purine and Pyrimidine Homo-Association.- 16.2 Base-Pair Configurations with Purine-Pyrimidine Hetero-Association: the Watson-Crick Base-Pairs.- 16.3 Base Pairs Can Combine to Form Triplets and Quadruplets.- 17 Hydrogen Bonding in the Crystal Structures of the Nucleosides and Nucleotides.- 17.1 Conformational and Hydrogen-Bonding Characteristics of the Nucleosides and Nucleotides.- 17.2 A Selection of Cyclic Hydrogen-Bonding Patterns Formed in Nucleoside and Nucleotide Crystal Structures.- 17.3 General Hydrogen-Bonding Patterns in Nucleoside and Nucleotide Crystal Structures.- III Hydrogen Bonding in Biological Macromolecules.- 18 O-H *** O Hydrogen Bonding in Crystal Structures of Cyclic and Linear Oligoamyloses: Cyclodextrins, Maltotriose, and Maltohexaose.- 18.1 The Cyclodextrins and Their Inclusion Complexes.- 18.2 Crystal Packing Patterns of Cyclodextrins Are Determined by Hydrogen Bonding.- 18.3 Cyclodextrins as Model Compounds to Study Hydrogen-Bonding Networks.- 18.4 Cooperative, Homodromic, and Antidromic Hydrogen-Bonding Patterns in the ?-Cyclodextrin Hydrates.- 18.5 Homodromic and Antidromic O-H *** O Hydrogen-Bonding Systems Analyzed Theoretically.- 18.6 Intramolecular Hydrogen Bonds in the ?-Cyclodextrin Molecule are Variable - the Induced-Fit Hypothesis.- 18.7 Flip-Flop Hydrogen Bonds in ?-Cyclodextrin * 11 H2O.- 18.8 From Flip-Flop Disorder to Ordered Homodromic Arrangements at Low lbmperature: The Importance of the Cooperative Effect.- 18.9 Maltohexaose Polyiodide and Maltotriose - Double and Single Left-Handed Helices With and Without Intramolecular O(2) *** O(3?) Hydrogen Bonds.- 19 Hydrogen Bonding in Proteins.- 19.1 Geometry of Secondary-Structure Elements: Helix, Pleated Sheet, and Turn.- 19.2 Hydrogen-Bond Analysis in Protein Crystal Structures.- 19.3 Hydrogen-Bonding Patterns in the Secondary Structure Elements.- 19.4 Hydrogen-Bonding Patterns Involving Side-Chains.- 19.5 Internal Water Molecules as Integral Part of Protein Structures.- 19.6 Metrical Analysis of Hydrogen Bonds in Proteins.- 19.7 Nonsecondary-Structure Hydrogen-Bond Geometry Between Main-Chains, Side-Chains and Water Molecules.- 19.8 Three-Center (Bifurcated) Bonds in Proteins.- 19.9 Neutron Diffraction Studies on Proteins Give Insight into Local Hydrogen-Bonding Flexibility.- 19.10 Site-Directed Mutagenesis Gives New Insight into Protein Thermal Stability and Strength of Hydrogen Bonds.- 20 The Role of Hydrogen Bonding in the Structure and Function of the Nucleic Acids.- 20.1 Hydrogen Bonding in Nucleic Acids is Essential for Life.- 20.2 The Structure of DNA and RNA Double Helices is Determined by Watson-Crick Base-Pair Geometry.- 20.3 Systematic and Accidental Base-Pair Mismatches: "Wobbling" and Mutations.- 20.4 Noncomplementary Base Pairs Have a Structural Role in tRNA.- 20.5 Homopolynucleotide Complexes Are Stabilized by a Variety of Base-Base Hydrogen Bonds - Three-Center (Bifurcated) Hydrogen Bonds in A-Tracts.- 20.6 Specific Protein-Nucleic Acid Recognition Involves Hydrogen Bonding.- IV Hydrogen Bonding by the Water Molecule.- 21 Hydrogen-Bonding Patterns in Water, Ices, the Hydrate Inclusion Compounds, and the Hydrate Layer Structures.- 21.1 Liquid Water and the Ices.- 21.2 The Hydrate Inclusion Compounds.- 21.3 Hydrate Layer Structures.- 22 Hydrates of Small Biological Molecules: Carbohydrates, Amino Acids, Peptides, Purines, Pyrimidines, Nucleosides and Nucleotides.- 23 Hydration of Proteins.- 23.1 Characterization of "Bound Water" at Protein Surfaces - the First Hydration Shell.- 23.2 Sites of Hydration in Proteins.- 23.3 Metrics of Water Hydrogen Bonding to Proteins.- 23.4 Ordered Water Molecules at Protein Surfaces - Clusters and Pentagons.- 24 Hydration of Nucleic Acids.- 24.1 Two Water Layers Around the DNA Double Helix.- 24.2 Crystallographically Determined Hydration Sites in A-, B-, Z-DNA. A Statistical Analysis.- 24.3 Hydration Motifs in Double Helical Nucleic Acids.- 24.3.1 Sequence-Independent Motifs.- 24.3.2 Sequence-Dependent Motifs.- 24.4 DNA Hydration and Structural Transitions Are Correlated: Some Hypotheses.- 25 The Role of Three-Center Hydrogen Bonds in the Dynamics of Hydration and of Structure Transition.- References.- Refcodes.

Hydrogen Bonding in Biological Structures

Ab initio LCAO–MO–SCF calculation for H2CO···H2O is carried out with a minimal Slater basis set. The most stable conformation has an O···H distance of 1.89 A with <C=O···H=− 64° and a stabilization energy of 3.5 kcal/mole, about a half of that for H2O···H2O. Nonlinear and π hydrogen bonds, H2CO···2H2O and the O···H–C hydrogen bond in H2O···HCHO, are also studied. An energy decomposition scheme is proposed and applied to H2CO···H2O and H2O···H2O. In the latter the electrostatic energy 8.0 kcal/mole, the exchange repulsion − 9.9 kcal/mole, the polarization and dispersion energy 0.3 kcal/mole, and the delocalization energy 8.2 kcal/mole are in good agreement with Coulson's estimates.

Molecular Orbital Studies of Hydrogen Bonds. III. C=O···H–O Hydrogen Bond in H2CO···H2O and H2CO···2H2O

Molecular dynamics calculations on a classical model for liquid water have been carried out at mass density 1 g/cm3 and at four temperatures. The effective pair potential employed is based on a four‐charge model for each molecule and represents a modification of the prior ``BNS'' interaction. Results for molecular structure and thermodynamic properties indicate that the modification improves the fidelity of the molecular dynamics simulation. In particular, the present version leads to a density maximum near 27°C for the liquid in coexistence with its vapor and to molecular distribution functions in better agreement with x‐ray scattering experiments.

Improved simulation of liquid water by molecular dynamics

Accurate SCF calculations have been carried out to investigate the potential of interaction for pairs and triplets of water molecules. The most stable pair configuration involves a linear hydrogen bond of length ROO = 3.00 A and strength 4.72 kcal/mole. Three‐molecule nonadditivities are large in magnitude and vary in sign according to the hydrogen‐bond pattern involved. In both aqueous liquids and solids, the net trimer nonadditivity effect amounts to increased binding energy, decreased neighbor distance, and slightly enhanced tendency toward perfect tetrahedral coordination symmetry. The nonadditivity furthermore is inconsistent with the phenomenology of simple mutual electrostatic polarization between neighboring molecules.

D. Hankins

Papers

Water Molecule Interactions

Hydrogen-bond energy nonadditivity in water☆

Erratum: Water molecule interactions