Search

With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. In the present work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn–Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer. [ABSTRACT FROM AUTHOR]

The many-body simulation of quantum systems is an active field of research that involves several different methods targeting various computing platforms. Many methods commonly employed, particularly coupled cluster methods, have been adapted to leverage the latest advances in modern high-performance computing. Selected configuration interaction (sCI) methods have seen extensive usage and development in recent years. However, the development of sCI methods targeting massively parallel resources has been explored only in a few research works. Here, we present a parallel, distributed memory implementation of the adaptive sampling configuration interaction approach (ASCI) for sCI. In particular, we will address the key concerns pertaining to the parallelization of the determinant search and selection, Hamiltonian formation, and the variational eigenvalue calculation for the ASCI method. Load balancing in the search step is achieved through the application of memory-efficient determinant constraints originally developed for the ASCI-PT2 method. The presented benchmarks demonstrate near optimal speedup for ASCI calculations of Cr2 (24e, 30o) with 106, 107, and 3 × 108 variational determinants on up to 16 384 CPUs. To the best of the authors' knowledge, this is the largest variational ASCI calculation to date. [ABSTRACT FROM AUTHOR]

For many computational chemistry packages, being able to efficiently and effectively scale across an exascale cluster is a heroic feat. Collective experience from the Department of Energy's Exascale Computing Project suggests that achieving exascale performance requires far more planning, design, and optimization than scaling to petascale. In many cases, entire rewrites of software are necessary to address fundamental algorithmic bottlenecks. This in turn requires a tremendous amount of resources and development time, resources that cannot reasonably be afforded by every computational science project. It thus becomes imperative that computational science transition to a more sustainable paradigm. Key to such a paradigm is modular software. While the importance of modular software is widely recognized, what is perhaps not so widely appreciated is the effort still required to leverage modular software in a sustainable manner. The present manuscript introduces PluginPlay, https://github.com/NWChemEx-Project/PluginPlay, an inversion-of-control framework designed to facilitate developing, maintaining, and sustaining modular scientific software packages. This manuscript focuses on the design aspects of PluginPlay and how they specifically influence the performance of the resulting package. Although, PluginPlay serves as the framework for the NWChemEx package, PluginPlay is not tied to NWChemEx or even computational chemistry. We thus anticipate PluginPlay to prove to be a generally useful tool for a number of computational science packages looking to transition to the exascale. [ABSTRACT FROM AUTHOR]

We develop a stochastic resolution of identity representation to the second-order Matsubara Green's function (sRI-GF2) theory. Using a stochastic resolution of the Coulomb integrals, the second order Born self-energy in GF2 is decoupled and reduced to matrix products/contractions, which reduces the computational cost from O(N5) to O(N3) (with N being the number of atomic orbitals). The current approach can be viewed as an extension to our previous work on stochastic resolution of identity second order Møller-Plesset perturbation theory [T. Y. Takeshita et al., J. Chem. Theory Comput. 13, 4605 (2017)] and offers an alternative to previous stochastic GF2 formulations [D. Neuhauser et al., J. Chem. Theory Comput. 13, 5396 (2017)]. We show that sRI-GF2 recovers the deterministic GF2 results for small systems, is computationally faster than deterministic GF2 for N > 80, and is a practical approach to describe weak correlations in systems with 103 electrons and more. [ABSTRACT FROM AUTHOR]

Data-driven prediction of molecular properties presents unique challenges to the design of machine learning methods concerning data structure/dimensionality, symmetry adaption, and confidence management. In this paper, we present a kernel-based pipeline that can learn and predict the atomization energy of molecules with high accuracy. The framework employs Gaussian process regression to perform predictions based on the similarity between molecules, which is computed using the marginalized graph kernel. To apply the marginalized graph kernel, a spatial adjacency rule is first employed to convert molecules into graphs whose vertices and edges are labeled by elements and interatomic distances, respectively. We then derive formulas for the efficient evaluation of the kernel. Specific functional components for the marginalized graph kernel are proposed, while the effects of the associated hyperparameters on accuracy and predictive confidence are examined. We show that the graph kernel is particularly suitable for predicting extensive properties because its convolutional structure coincides with that of the covariance formula between sums of random variables. Using an active learning procedure, we demonstrate that the proposed method can achieve a mean absolute error of 0.62 ± 0.01 kcal/mol using as few as 2000 training samples on the QM7 dataset. [ABSTRACT FROM AUTHOR]

A new large-scale parallel multiconfigurational self-consistent field (MCSCF) implementation in the open-source NWChem computational chemistry code is presented. The generalized active space approach is used to partition large configuration interaction (CI) vectors and generate a sufficient number of batches that can be distributed to the available cores. Massively parallel CI calculations with large active spaces can be performed. The new parallel MCSCF implementation is tested for the chromium trimer and for an active space of 20 electrons in 20 orbitals, which can now routinely be performed. Unprecedented CI calculations with an active space of 22 electrons in 22 orbitals for the pentacene systems were performed and a single CI iteration calculation with an active space of 24 electrons in 24 orbitals for the chromium tetramer was possible. The chromium tetramer corresponds to a CI expansion of one trillion Slater determinants (914 058 513 424) and is the largest conventional CI calculation attempted up to date. [ABSTRACT FROM AUTHOR]

A combined theoretical and solid-state 17O nuclear magnetic resonance (NMR) study of the electronic structure of the uranyl ion UO22+ in (NH4)4UO2(CO3)3 and rutherfordine (UO2CO3) is presented, the former representing a system with a hydrogen-bonding environment around the uranyl oxygens and the latter exemplifying a uranyl environment without hydrogens. Relativistic density functional calculations reveal unique features of the U–O covalent bond, including the finding of 17O chemical shift anisotropies that are among the largest for oxygen ever reported (>1200 ppm). Computational results for the oxygen electric field gradient tensor are found to be consistently larger in magnitude than experimental solid-state 17O NMR measurements in a 7.05 T magnetic field indicate. A modified version of the Solomon theory of the two-spin echo amplitude for a spin-5/2 nucleus is developed and applied to the analysis of the 17O echo signal of U 17O22+. [ABSTRACT FROM AUTHOR]

Dynamic polarizabilities for open- and closed-shell molecules were obtained by using coupled-cluster (CC) linear response theory with full treatment of singles, doubles, and triples (CCSDT-LR) with large basis sets utilizing the NWChem software suite. By using four approximate CC methods in conjunction with augmented cc-pVNZ basis sets, we are able to evaluate the convergence in both many-electron and one-electron spaces. For systems with primarily dynamic correlation, the results for CC3 and CCSDT are almost indistinguishable. For systems with significant static correlation, the CC3 tends to overestimate the triples contribution, while the PS(T) approximation [J. Chem. Phys. 127, 164105 (2007)] produces mixed results that are heavily dependent on the accuracies provided by noniterative approaches used to correct the equation-of-motion CCSD excitation energies. Our results for open-shell systems show that the choice of reference (restricted open-shell Hartree–Fock versus unrestricted Hartree–Fock) can have a significant impact on the accuracy of polarizabilities. A simple extrapolation based on pentuple-zeta CCSD calculations and triple-zeta CCSDT calculations reproduces experimental results with good precision in most cases. [ABSTRACT FROM AUTHOR]

First principles molecular dynamics simulations of the hydration shells surrounding UO22+ ions are reported for temperatures near 300 K. Most of the simulations were done with 64 solvating water molecules (22 ps). Simulations with 122 water molecules (9 ps) were also carried out. The hydration structure predicted from the simulations was found to agree with very well-known results from x-ray data. The average U==O bond length was found to be 1.77 Å. The first hydration shell contained five trigonally coordinated water molecules that were equatorially oriented about the O–U–O axis with the hydrogen atoms oriented away from the uranium atom. The five waters in the first shell were located at an average distance of 2.44 Å (2.46 Å, 122 water simulation). The second hydration shell was composed of distinct equatorial and apical regions resulting in a peak in the U–O radial distribution function at 4.59 Å. The equatorial second shell contained ten water molecules hydrogen bonded to the five first shell molecules. Above and below the UO22+ ion, the water molecules were found to be significantly less structured. In these apical regions, water molecules were found to sporadically hydrogen bond to the oxygen atoms of the UO22+, oriented in such a way as to have their protons pointed toward the cation. While the number of apical waters varied greatly, an average of five to six waters was found in this region. Many water transfers into and out of the equatorial and apical second solvation shells were observed to occur on a picosecond time scale via dissociative mechanisms. Beyond these shells, the bonding pattern substantially returned to the tetrahedral structure of bulk water. [ABSTRACT FROM AUTHOR]

The electronic structure and spectroscopic properties (Re, ωe, ωexe, βe, and Te) of the ground state and the 22 lowest excited states of chlorine molecule were studied within a four-component relativistic framework using the MOLFDIR program package. The potential energy curves of all possible 23 covalent states were calculated using relativistic complete open shell configuration interaction approach. In addition, four component multireference configuration interaction with single and double excitation calculations were performed in order to infer the effects due to dynamical correlation in vertical excitations. The calculated properties are in good agreement with the available experimental data. [ABSTRACT FROM AUTHOR]

This paper discusses practical scheme for correcting the linear response coupled cluster with singles and doubles (CCSD) equations by shifting the poles corresponding to the equation-of-motion CCSD excitation energies by adding noniterative corrections due to triples. A simple criterion is derived for the excited states to be corrected in the spectral resolution of similarity transformed Hamiltonian on the CCSD level. Benchmark calculations were performed to compare the accuracies of static and dynamic polarizabilities obtained in this way with the CC3 and CCSDT counterparts. [ABSTRACT FROM AUTHOR]

Coupled-cluster theory with single and double excitations is applied to the calculation of optical properties of large polyaromatic hydrocarbons. Dipole polarizabilities are reported for benzene, pyrene, and the oligoacenes sequence n=2–6. Dynamic polarizabilities were calculated on polyacences as large as pentacene for a single frequency and for benzene and pyrene at many frequencies. The basis set effect was studied for benzene using a variety of basis sets in the Pople [Theor. Chim. Acta 28, 213 (1973)] and Dunning [J. Chem. Phys. 90, 1007 (1989)] families up to aug-cc-pVQZ and the Sadlej pVTZ basis [Collect. Czech. Chem. Commun. 53, 1995 (1998)], which was used exclusively for the largest molecules. Geometries were optimized using HF, B3LYP, PBE0, and MP2 and compared to experiment to measure method dependence and the possible role of bond-length alternation. Finally, the polarizability results were compared to four common density functionals (B3LYP, BLYP, PBE0, PBE). [ABSTRACT FROM AUTHOR]

The self-consistent charge density-functional tight-binding (SCC-DFTB) method is employed for studying various molecular properties of small fullerenes: C28, C60, and C70. The computed bond distances, vibrational infrared and Raman spectra, vibrational densities of states, and electronic densities of states are compared with experiment (where available) and density-functional theory (DFT) calculations using various basis sets. The presented DFT benchmark calculations using the correlation-consistent polarized valence triple zeta basis set are at present the most extensive calculations on harmonic frequencies of these species. Possible limitations of the SCC-DFTB method for the prediction of molecular vibrational and optical properties are discussed. The presented results suggest that SCC-DFTB is a computationally feasible and reliable method for predicting vibrational and electronic properties of such carbon nanostructures comparable in accuracy with small to medium size basis set DFT calculations at the computational cost of standard semiempirical methods. [ABSTRACT FROM AUTHOR]

We report the first harmonic vibrational spectra for each of the lowest lying isomers within the four major families of minima of (H2O)20, namely, the dodecahedron, fused cubes, face-sharing pentagonal prisms, and edge-sharing pentagonal prisms. These were obtained at the second-order Mo\ller–Plesset perturbation level of theory (MP2) with the augmented correlation consistent basis set of double zeta quality (aug-cc-pVDZ) at the corresponding minimum energy geometries. The computed infrared (IR) spectra are the first ones obtained from first principles for these clusters. They were found to contain spectral features, which can be directly mapped onto the distinctive spectroscopic signatures of their constituent tetramer, pentamer, and octamer fragments. The dodecahedron spectra show the richest structure in the OH stretching region and are associated with the most redshifted OH vibrations with respect to the monomer. The lowest lying edge-sharing pentagonal prism isomer displays intense IR active vibrations that are redshifted by ∼600 cm-1 with respect to the water monomer. Furthermore the most redshifted, IR-active OH stretching vibrations for all four networks correspond to hydrogen bonded OH groups, which exhibit the following two common characteristics: (i) they belong to fragments which have a “free” OH stretch and (ii) they act as donors to a neighboring water molecule along a “dimerlike” (strong) hydrogen bond. The zero-point energy corrected MP2/CBS (complete basis set) limit binding energies D0 for the four isomers are -163.1 kcal/mol (edge-sharing pentagonal prism), -160.1 kcal/mol (face-sharing pentagonal prism), -157.5 kcal/mol (fused cubes), and -148.1 kcal/mol (dodecahedron). [ABSTRACT FROM AUTHOR]

Coupled-cluster methods including through and up to the connected single, double, triple, and quadruple substitutions have been derived and implemented automatically for sequential and parallel executions by an algebraic and symbolic manipulation program TCE (TENSOR CONTRACTION ENGINE) for use in conjunction with a one-component third-order Douglas–Kroll approximation for relativistic corrections. A combination of the converging electron-correlation methods, the accurate relativistic reference wave functions, and the use of systematic basis sets tailored to the relativistic approximation has been shown to predict the experimental singlet–triplet separations within 0.02 eV (0.5 kcal/mol) for five triatomic hydrides (CH[sub 2], NH[sub 2][sup +], SiH[sub 2], PH[sub 2][sup +], and AsH[sub 2][sup +]), the experimental bond lengths (r[sub e] or r[sub 0]) within 0.002 Å, rotational constants (B[sub e] or B[sub 0]) within 0.02 cm[sup -1], vibration–rotation constants (α[sub e]) within 0.01 cm[sup -1], centrifugal distortion constants (D[sub e]) within 2%, harmonic vibration frequencies (ω[sub e]) within 8 cm[sup -1] (0.4%), anharmonic vibrational constants (xω[sub e]) within 2 cm[sup -1], and dissociation energies (D[sub 0][sup 0]) within 0.02 eV (0.4 kcal/mol) for twenty diatomic hydrides (BH, CH, NH, OH, FH, AlH, SiH, PH, SH, ClH, GaH, GeH, AsH, SeH, BrH, InH, SnH, SbH, TeH, and IH) containing main-group elements across the second through fifth rows of the periodic table. In these calculations, spin–orbit effects on dissociation energies, which were assumed to be additive, were estimated from the measured spin–orbit coupling constants of atoms and diatomic molecules, and an electronic energy in the complete-basis-set, complete-electron-correlation limit has been extrapolated in two ways to verify the robustness of the results: One assuming Gaussian-exponential dependence of total energies on double through quadruple ζ basis sets and the other assuming n[sup -3] dependence of correlation energies on double through quintuple ζ basis sets. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Atomization energies at 0 K and heats of formation at 298 K were obtained for a collection of small halogenated molecules from coupled cluster theory including noniterative, quasiperturbative triple excitations calculations with large basis sets (up through augmented septuple zeta quality in some cases). In order to achieve near chemical accuracy (± 1 kcal/mol) in the thermodynamic properties, we adopted a composite theoretical approach which incorporated estimated complete basis set binding energies based on frozen core coupled cluster theory energies and (up to) five corrections: (1) a core/valence correction; (2) a Douglas-Kroll-Hess scalar relativistic correction; (3) a first-order atomic spin-orbit correction; (4) a second-order spin-orbit correction for heavy elements; and (5) an approximate correction to account for the remaining correlation energy. The last of these corrections is based on a recently proposed approximation to full configuration interaction via a continued fraction approximant for coupled cluster theory [CCSD(T)-cf]. Failure to consider corrections (1) to (4) can introduce errors significantly in excess of the target accuracy of ± 1 kcal/mol. Although some cancellation of error may occur if one or more of these corrections is omitted, such a situation is by no means universal and cannot be relied upon for high accuracy. The accuracy of the Douglas-Kroll-Hess approach was calibrated against both new and previously published four-component Dirac Coulomb results at the coupled cluster level of theory. In addition, vibrational zero-point energies were computed at the coupled cluster level of theory for those polyatomic systems lacking an experimental anharmonic value. [ABSTRACT FROM AUTHOR]

New theoretical predictions for the static and frequency dependent polarizabilities of C60 are reported. Using the linear response coupled cluster approach with singles and doubles and a basis set especially designed to treat the molecular properties in external electric field, we obtained 82.20 and 83.62 Å3 for static and dynamic (λ=1064 nm) polarizabilities. These numbers are in a good agreement with experimentally inferred data of 76.5±8 and 79±4 Å3 [R. Antoine et al., J. Chem. Phys.110, 9771 (1999); A. Ballard et al., J. Chem. Phys.113, 5732 (2000)]. The reported results were obtained with the highest wave function-based level of theory ever applied to the C60 system. [ABSTRACT FROM AUTHOR]

We present an EOMCCSD (equation of motion coupled cluster with singles and doubles) study of excited states of the small [UO2]2+ and [UO2]+ model systems as well as the larger UVIO2(saldien) complex. In addition, the triples contribution within the EOMCCSDT and CR-EOMCCSD(T) (completely renormalized EOMCCSD with non-iterative triples) approaches for the [UO2]2+ and [UO2]+ systems as well as the active-space variant of the CR-EOMCCSD(T) method-CR-EOMCCSd(t)-for the UVIO2(saldien) molecule are investigated. The coupled cluster data were employed as benchmark to choose the 'best' appropriate exchange-correlation functional for subsequent time-dependent density functional (TD-DFT) studies on the transition energies for closed-shell species. Furthermore, the influence of the saldien ligands on the electronic structure and excitation energies of the [UO2]+ molecule is discussed. The electronic excitations as well as their oscillator dipole strengths modeled with TD-DFT approach using the CAM-B3LYP exchange-correlation functional for the [UVO2(saldien)]- with explicit inclusion of two dimethyl sulfoxide molecules are in good agreement with the experimental data of Takao et al. [Inorg. Chem. 49, 2349 (2010)]. [ABSTRACT FROM AUTHOR]

Electronic structure calculations and direct chemical dynamics simulations are used to study the formation and decomposition of dioxetane on its ground state singlet potential energy surface. The stationary points for 1O2 + C2H4, the singlet ·O-O-CH2-CH2· biradical, the transition state (TS) connecting this biradical with dioxetane, and the two transition states and gauche ·O-CH2-CH2-O· biradical connecting dioxetane with the formaldehyde product molecules are investigated at different levels of electronic structure theory including UB3LYP, UMP2, MRMP2, and CASSCF and a range of basis sets. The UB3LYP/6-31G* method was found to give representative energies for the reactive system and was used as a model for the simulations. UB3LYP/6-31G* direct dynamics trajectories were initiated at the TS connecting the ·O-O-CH2-CH2· biradical and dioxetane by sampling the TS's vibrational energy levels, and rotational and reaction coordinate energies, with Boltzmann distributions at 300, 1000, and 1500 K. This corresponds to the transition state theory model for trajectories that pass the TS. The trajectories were directed randomly towards both the biradical and dioxetane. A small fraction of the trajectories directed towards the biradical recrossed the TS and formed dioxetane. The remainder formed 1O2 + C2H4 and of these ∼ 40% went directly from the TS to 1O2 + C2H4 without getting trapped and forming an intermediate in the ·O-O-CH2-CH2· biradical potential energy minimum, a non-statistical result. The dioxetane molecules which are formed dissociate to two formaldehyde molecules with a rate constant two orders of magnitude smaller than that predicted by Rice-Ramsperger-Kassel-Marcus theory. The reaction dynamics from dioxetane to the formaldehyde molecules do not follow the intrinsic reaction coordinate or involve trapping in the gauche ·O-CH2-CH2-O· biradical potential energy minimum. Important non-statistical dynamics are exhibited for this reactive system. [ABSTRACT FROM AUTHOR]

New theoretical predictions for the static and frequency dependent polarizabilities of C60 are reported. Using the linear response coupled cluster approach with singles and doubles and a basis set especially designed to treat the molecular properties in external electric field, we obtained 82.20 and 83.62 Å3 for static and dynamic (λ=1064 nm) polarizabilities. These numbers are in a good agreement with experimentally inferred data of 76.5±8 and 79±4 Å3 [R. Antoine et al., J. Chem. Phys.110, 9771 (1999); A. Ballard et al., J. Chem. Phys.113, 5732 (2000)]. The reported results were obtained with the highest wave function-based level of theory ever applied to the C60 system. [ABSTRACT FROM AUTHOR]

The self-consistent charge density-functional tight-binding (SCC-DFTB) method is employed for studying various molecular properties of small fullerenes: C28, C60, and C70. The computed bond distances, vibrational infrared and Raman spectra, vibrational densities of states, and electronic densities of states are compared with experiment (where available) and density-functional theory (DFT) calculations using various basis sets. The presented DFT benchmark calculations using the correlation-consistent polarized valence triple zeta basis set are at present the most extensive calculations on harmonic frequencies of these species. Possible limitations of the SCC-DFTB method for the prediction of molecular vibrational and optical properties are discussed. The presented results suggest that SCC-DFTB is a computationally feasible and reliable method for predicting vibrational and electronic properties of such carbon nanostructures comparable in accuracy with small to medium size basis set DFT calculations at the computational cost of standard semiempirical methods. [ABSTRACT FROM AUTHOR]

A combined theoretical and solid-state (17)O nuclear magnetic resonance (NMR) study of the electronic structure of the uranyl ion UO(2)(2+) in (NH(4))(4)UO(2)(CO(3))(3) and rutherfordine (UO(2)CO(3)) is presented, the former representing a system with a hydrogen-bonding environment around the uranyl oxygens and the latter exemplifying a uranyl environment without hydrogens. Relativistic density functional calculations reveal unique features of the U-O covalent bond, including the finding of (17)O chemical shift anisotropies that are among the largest for oxygen ever reported (>1200 ppm). Computational results for the oxygen electric field gradient tensor are found to be consistently larger in magnitude than experimental solid-state (17)O NMR measurements in a 7.05 T magnetic field indicate. A modified version of the Solomon theory of the two-spin echo amplitude for a spin-5/2 nucleus is developed and applied to the analysis of the (17)O echo signal of U (17)O(2)(2+).

The self-consistent charge density-functional tight-binding (SCC-DFTB) method is employed for studying various molecular properties of small fullerenes: C(28), C(60), and C(70). The computed bond distances, vibrational infrared and Raman spectra, vibrational densities of states, and electronic densities of states are compared with experiment (where available) and density-functional theory (DFT) calculations using various basis sets. The presented DFT benchmark calculations using the correlation-consistent polarized valence triple zeta basis set are at present the most extensive calculations on harmonic frequencies of these species. Possible limitations of the SCC-DFTB method for the prediction of molecular vibrational and optical properties are discussed. The presented results suggest that SCC-DFTB is a computationally feasible and reliable method for predicting vibrational and electronic properties of such carbon nanostructures comparable in accuracy with small to medium size basis set DFT calculations at the computational cost of standard semiempirical methods.

We report the first harmonic vibrational spectra for each of the lowest lying isomers within the four major families of minima of (H2O)20, namely, the dodecahedron, fused cubes, face-sharing pentagonal prisms, and edge-sharing pentagonal prisms. These were obtained at the second-order Moller-Plesset perturbation level of theory (MP2) with the augmented correlation consistent basis set of double zeta quality (aug-cc-pVDZ) at the corresponding minimum energy geometries. The computed infrared (IR) spectra are the first ones obtained from first principles for these clusters. They were found to contain spectral features, which can be directly mapped onto the distinctive spectroscopic signatures of their constituent tetramer, pentamer, and octamer fragments. The dodecahedron spectra show the richest structure in the OH stretching region and are associated with the most redshifted OH vibrations with respect to the monomer. The lowest lying edge-sharing pentagonal prism isomer displays intense IR active vibrations that are redshifted by approximately 600 cm(-1) with respect to the water monomer. Furthermore the most redshifted, IR-active OH stretching vibrations for all four networks correspond to hydrogen bonded OH groups, which exhibit the following two common characteristics: (i) they belong to fragments which have a "free" OH stretch and (ii) they act as donors to a neighboring water molecule along a "dimerlike" (strong) hydrogen bond. The zero-point energy corrected MP2/CBS (complete basis set) limit binding energies D(0) for the four isomers are -163.1 kcal/mol (edge-sharing pentagonal prism), -160.1 kcal/mol (face-sharing pentagonal prism), -157.5 kcal/mol (fused cubes), and -148.1 kcal/mol (dodecahedron).

Coupled-cluster methods including through and up to the connected single, double, triple, and quadruple substitutions have been derived and implemented automatically for sequential and parallel executions by an algebraic and symbolic manipulation program TCE (TENSOR CONTRACTION ENGINE) for use in conjunction with a one-component third-order Douglas-Kroll approximation for relativistic corrections. A combination of the converging electron-correlation methods, the accurate relativistic reference wave functions, and the use of systematic basis sets tailored to the relativistic approximation has been shown to predict the experimental singlet-triplet separations within 0.02 eV (0.5 kcal/mol) for five triatomic hydrides (CH2, NH2+, SiH2, PH2+, and AsH2+), the experimental bond lengths (re or r0) within 0.002 angstroms, rotational constants (Be or B0) within 0.02 cm(-1), vibration-rotation constants (alphae) within 0.01 cm(-1), centrifugal distortion constants (De) within 2%, harmonic vibration frequencies (omegae) within 8 cm(-1) (0.4%), anharmonic vibrational constants (xomegae) within 2 cm(-1), and dissociation energies (D0(0)) within 0.02 eV (0.4 kcal/mol) for twenty diatomic hydrides (BH, CH, NH, OH, FH, AlH, SiH, PH, SH, ClH, GaH, GeH, AsH, SeH, BrH, InH, SnH, SbH, TeH, and IH) containing main-group elements across the second through fifth rows of the periodic table. In these calculations, spin-orbit effects on dissociation energies, which were assumed to be additive, were estimated from the measured spin-orbit coupling constants of atoms and diatomic molecules, and an electronic energy in the complete-basis-set, complete-electron-correlation limit has been extrapolated in two ways to verify the robustness of the results: One assuming Gaussian-exponential dependence of total energies on double through quadruple zeta basis sets and the other assuming n(-3) dependence of correlation energies on double through quintuple zeta basis sets., (Copyright 2004 American Institute of Physics)