Calculations on the vertical and adiabatic ionization energies of (H2S)2

After determining reliable procedures for calculating the ionization energy of H₂S using many-body perturbation theory and electron propagator theory, the same procedures are used to calculate the vertical ionization energy of the van der Waals molecule (H₂S)₂. The adiabatic ionization energy calculated for a dimer cation with a H₃S⁺...SH structure is in satisfactory agreement with photoionization threshold experiments. An alternative dimer cation structure with a three-electron S-S bond is also discussed.     

A note on the convergence of multiconfigurational many-body perturbation theory

The convergence of multiconfigurational many-body perturbation theory (MC MBPT ) is discussed in connection with the intruder state. Its convergence properties are first examined with a fictitious three-level system employing a Hermitian version of MC MBPT , which permits a general model space. It is then applied to the H₂—H₂ and N₂ systems. The results suggest that a more extensive model space is likely to embrace new intruder states and the space extension be executed with due caution.     

Calculation of vertical ionization potentials of H2O and Ne by many-body Rayleigh-Schrödinger perturbation theory

The vertical ionization potentials of H₂O and Ne are calculated by many-body Rayleigh-Schrödinger perturbation theory up to third order. The comparison of the present method with the other approaches is done.     

Electron correlation effects on magnetic properties of BH

The nuclear magnetic shielding and magnetizability tensors for the BH molecule are calculated by the coupled-Hartree–Fock method and many-body perturbation theory. As in the case of H₂, HF, and F₂, the second-order non-CHF diagrams make an inappreciable contribution.     

Perturbative calculation of intermolecular interactions in orthogonalized or biorthogonal basis sets

Two versions of a many-body perturbation theory for the computation of molecular interaction energies are investigated. The methods are based on the partitioning of the second-quantized form of the dimer Hamiltonian written either in the orthogonalized basis of the monomer MOs, or, alternatively, in the original non-orthogonal dimer basis set handling the overlap by the biorthogonal formalism. The zeroth-order Hamiltonian H ⁰ is the sum of effective monomer Fockians and the zeroth-order wave functions are exact eigenfunctions of H ⁰. Full antisymmetry is ensured by the second-quantized formalism. Basis set superposition error is accounted for by the counterpoise correction recipe. Results for He₂, (H₂)₂ and (H₂O)₂ indicate the reliability of the biorthogonal technique.     

Semiclassical Trajectory Studies of Reactive and Nonreactive Scattering of OH(A 2Σ+) by H2 Based on an Improved Full-Dimensional Ab Initio Diabatic Potential Energy Matrix

We present a new full-dimensional diabatic potential energy matrix (DPEM) for electronically nonadiabatic collisions of OH(A ²Σ⁺) with H₂, and we calculate the probabilities of electronically adiabatic inelastic collisions, nonreactive quenching, and reactive quenching to form H₂O+H. The DPEM was fitted using a many-body expansion with permutationally invariant polynomials in bond-order functions to represent the many-body part. The dynamics calculations were carried out with the fewest-switches with time uncertainty and stochastic decoherence (FSTU/SD) semiclassical trajectory method. We present results both for head-on collisions (impact parameter b equal to zero) and for a full range of impact parameters. The results are compared to experiment and to earlier FSTU/SD and quantum dynamics calculations with a previously published DPEM. The various theoretical results all agree that nonreactive quenching dominates reactive quenching, but there are quantitative differences between the two DPEMs and between the b=0 results and the all-b results, especially for the probability of reactive quenching.     

Calculation of the correlation energy of molecules by many-body rayleigh-schrödinger perturbation theory up to third order

The correlation energy of the H₂O, BH, HF and Ne systems is calculated by many-body Rayleigh-Schrödinger perturbation theory up to third order, and compared with CI calculations in exactly the same basis sets.     

Calculation of the electronic properties of metal dimers

The low-energy excitation spectra of metal dimers are determined by solving exactly a realistic many-body Hamiltonian with inter- and intra-atomic Coulomb interactions. Our results for Cu₂ and Ag₂ are in very good agreement with the excitation energies derived from recent photodetachment experiments. The characteristics of the many-body excited states in these clusters are briefly discussed.     

Experimental and theoretical binding energy spectra and momentum distributions for the valence orbitals of H2O

The binding energy spectra (10–46 eV) and momentum distributions of the valence orbitals of H₂O have been measured using a new high-sensitivity binary (e,2e) electron spectrometer employing position-sensitive detectors. The binding energy spectrum shows a previously unreported feature at = 27 eV which is shown to be associated with the (2a₁)^?1 ionization process. The region between 25 and 46 eV is compared with previous (e,2e) and X-ray photoelectron measurements as well as with several existing and new many-body calculations indicating a splitting of the 2a₁ ionization pole strength. In addition the separate momentum distributions of the three outer valence orbitals of H₂O have been obtained from deconvoluted binding energy spectra run at a series of azimuthal angles. The results, which show considerably improved signal-to-noise ratio over earlier measurements using single-channel instrumentation are compared with spherically averaged momentum distributions calculated with a variety of wavefunctions.     

[Ni(bpy)(mal)(H2O)3]·H2O and [Ni(4,4′-dmbpy)(mal)(H2O)3]·1.5H2O: syntheses,crystal structures,magnetic properties,and computational study of stacking interactions

[Ni(bpy)(mal)(H₂O)₃]·H₂O and [Ni(4,4′-dmbpy)(mal)(H₂O)₃]·1.5H₂O (mal = maleato; bpy = 2,2′-bipyridine; 4,4′-dmbpy = 4,4′-dimethyl-2,2′-bipyridine) exhibit molecular crystal structures. The Ni(II) central ions in both complexes are six-coordinate by one chelate bonded L–L ligand, three aqua ligands, and one position is occupied by a maleato oxygen donor atom. Hydrogen bonded ribbon-like supramolecular structural motifs are present in both studied complexes; these are linked by weaker C–H?O hydrogen bonds in [Ni(bpy)(mal)(H₂O)₃]·H₂O, whereas in [Ni(4,4′-dmbpy)(mal)(H₂O)₃]·1,5H₂O the hydrogen bonded ribbons are linked by O?H-O-H?O hydrogen bonds with the participation of the additional water solvate molecule positioned on the twofold axis. In both structures, π–π stacking interactions with different patterns in respective structures were found. The role of dispersion energy and many-body effects in the stabilization of bpy and 4,4′-dmbpy stacking interactions were investigated using methods of computational chemistry. Those confirm the dispersion-dominated stabilization of the 4,4′-dmbpy supramolecular chain-like structure, with only marginal impact of cooperativity effects. Thermal decompositions of both complexes start with dehydration. Magnetic susceptibility studies performed from 2 to 300 K revealed a dominant effect of the zero-field splitting of the Ni(II) ion, governing the low-temperature magnetic properties of both compounds.     

Time-dependent Wave Packet Quantum Scattering Study of Reaction S(3P)+H2→HS+H on a New ab initio Potential Energy Surface 3A'

A new potential energy surface is presented for the triplet state ³A' of the chemical reaction S(³P)+H₂ from a set of accurate ab initio data. The single point energies are computed using highly correlated complete active space self-consistent-field and multi-reference config-uration interaction wave functions with a basis set of aug-cc-pV5Z. We have fitted the full set of energy values using many-body expansion method with an Aguado-Paniagua function. Based on the new potential energy surface, we carry out the time-dependent wave packet scattering calculations over the collision energy range of 0.8～2.2 eV. Both the centrifugal-sudden approximation and Coriolis Coupling cross sections are obtained. In addition, the total reaction probabilities are calculated for the reactant H₂ initially in the vibrational states v=0～3 (j=0). It is found that initial vibrational excitation enhances the title reaction.     

Ordinary field-theoretic methods for self-consistent wave functions which describe bond formation and dissociation. II. Commutative coupling case

Hartree–Bogoliubov–Valatin (HBV ) theory may be implemented with Lipkin Hamiltonians to obtain self-consistent BCS wave functions which describe bond formation and dissociation. These wave functions are in turn vacuua for Nambu's representation of Feynman–Dyson–Goldstone diagrammatic perturbation theory, and hence provide suitable references for the many-body treatment of correlation. Exact SCF solutions of the HBV equations are equivalent to special even-replacement MC –SCF solutions. The latter are similar to generalized valence bond theory, and require one Fock operator for each one-particle shell. The commutative coupling case of HBV theory is realized when the number-conserving renormalized one-body and number-nonconserving pairing operators commute. In this case, a set of orbital equations which involves a single Fock operator may be solved. Since this could prove to be a significant simplification for large systems, the commutative coupling and exact solutions are compared here for the fragmentation of H₂ and F₂. Results suggest that commutative coupling orbitals will be useful for the aforementioned many-body theory.     

Application of the many-body perturbation theory by using localized orbitals

Diagrammatic formulation of the many-body perturbation theory is investigated when both the occupied orbitals and the virtual ones are localized, i.e., they are unitary transforms of the canonical Hartree–Fock orbitals. All diagrams representing ground state correlation energy can be generated through fifth order. For cyclic polyenes C₆H₆ and C₁₀H₁₀ as model systems, the energy corrections are calculated in the Pariser–Parr–Pople approximation for a wide range of the coupling constant β^?1, through fourth order including some fifth order terms. The results are compared to those obtained by other methods: perturbation theory by using canonical orbitals and full CI. The effect of neglecting contributions from orbitals localized into neighboring sites is also studied.     

Near-degeneracy effects and efficiency of MBPT calculations of molecular properties. Polarizability curve of H2

The polarizability curve of H₂ is calculated by using the finite-field perturbation method. All self-consistency effects are accounted for at the HF level and many-body perturbation theory (MBPT) is used to evaluate the correlation contributions. Using a single HF determinant as a reference in MBPT calculations makes the near-degeneracy effects of essential importance on increasing the interatomic distance. Nevertheless, applying the MBPT scheme with appropriate fourth-order terms gives nearly exact values of both components of the polarizability tensor for interatomic distances up to ≈3.6 au.     

An investigation of the valence orbitals of water by high momentum resolution electron momentum spectroscopy

The momentum distributions of the valence orbitals for water well as the binding energy spectra in the region 10–45 eV have been reinvestigated with a high momentum resolution (≈0.1 a₀^?1 fwhm) binary (e.2e) spectrometer. The binding energy spectra show considerable satellite structure in the region > 25 eV which is consistent with theoretical predictions of final state configuration interaction (many-body effects) involving the (2a₁)^?1 hole state. An investigation of the momentum distribution in the satellite region confirms this assignment. This is in accord with a variety of recent theoretical studies and also consistent with earlier experiments. Differences suggested in earlier comparisons between theory and low momentum resolution experiments for the momentum distributions of the 1b₁ and 3a₁ orbitals have been verified. Several possible theoretical studies are suggested to investigate further this discrepancy between experiment and theory. Bonding effects and thenature of the molecular orbitals of H₂O in momentum space are also discussed.     

Looking at orbitals in the laboratory: The experimental investigation of molecular wavefunctions and binding energies by electron momentum spectroscopy

The experimental technique of electron momentum spectroscopy (EMS ) (i.e., binary (e, 2e) spectroscopy) is discussed together with typical examples of its applications over the past decade in the area of experimental quantum chemistry. Results interpreted within the framework of the plane wave impulse and the target Hartree—Fock approximations provide direct measurements of, spherically averaged, orbital electron momentum distributions. Results for a variety of atoms and small molecules are compared with calculations using a range of Fourier transformed SCF position space wavefunctions of varying sophistication. Measured momentum distributions (MD ) provide a “direct” view of orbitals. In addition to offering a sensitive experimental diagnostic for semiempirical molecular wavefunctions, the MD's provide a chemically significant, additional experimental constraint to the usual variational optimization of wavefunctions. The measured MD's clearly reflect well known characteristics of various chemical and physical properties. It appears that EMS and momentum space chemistry offer the promise of supplementary perspectives and new vistas in quantum chemistry, as suggested by Coulson more than 40 years ago. Binding energy spectra in the inner valence region reveal, in many cases, a major breakdown of the simple MO model for ionization in accord with the predictions of many-body calculations. Results are considered for atomic targets, including H and the noble gases. The measured momentum distribution for H₂ is also compared with results from Compton scattering. Results for H₂ and H are combined to provide a direct experimental assessment of the bond density in H₂, which is compared with calculations. The behavior of the outer valence MD ''s for small row two and row three hydride molecules such as H₂O and H₂S, NH₃, HF, and HCl are consistent with well known differences in chemical and physical behavior such as ligand-donor activity and hydrogen bonding. MD measurements for the outermost valence orbitals of HF, H₂O and NH₃ show significant differences from those calculated using even very high-quality wavefunctions. Measurements of MD's for outer σ_g orbitals of small polyatomic molecules such as CO₂, COS, CS₂, and CF₄ show clear evidence of mixed s and p character. It is apparent that EMS is a sensitive probe of details of electronic structure and electron motion in atoms and molecules.     

Correlation energies in open shell systems. Comparison of CEPA,PNO-CI and perturbation treatments based on the restricted Roothaan-Hartree-Fock formalism

A set of simple molecules in closed and open-shell ground states is treated by the three techniques mentioned in the title, using the same geometries and basis sets (DZ + P). It is found that for nearly all molecules treated in this study (exceptions are H₂ and CH₃) consistently about 98% of the CEPA valence shell correlation energy is obtained by third-order many-body Ray-leigh-Schrödinger perturbation theory (MB-RSPT). The CEPA and MB-RSPT results for reaction energies and barrier heights for some simple reactions differ by 0 to 30 kJ/mol, the CEPA results being in most cases closer to experiment than MB-RSPT, while CI results are much less reliable as long as CI is limited to singly and doubly substituted configurations only.     

Revisiting small clusters of water molecules

The geometries and relative energies of small clusters of water molecules, (H₂O)_n with 4 ⩽ n ⩽ 8, are reported. For each value of n we have considered the conformations corresponding to the lowest-energy minimum and those in nearby relative minima. Thus we report on six tetramers, four pentamers, six hexanlers, four heptamers, and eigth octamers. The geometrical conformations have been obtained using the Metropolis Monte Carlo method as a minimization technique, where the interaction energy is computed with the MCY potential plus three- and four-body corrections previously discussed. All the reported structures for a given cluster size are found to be close in energy. For the lowest conformation the geometry was optimized with ab initio SCF computations using energy gradients. Our results are compared with previous theoretical studies. We discuss the convergence of the interaction potential for liquid water when expressed in terms of a many-body series expansion.     

Monte Carlo study of electron correlation functions for small molecules

A study is made of electron-electron correlation functions for use in trial wave functions for small molecules. New forms are proposed that have only a few variational parameters, and these parameters have physical meanings that are easily discerned. Total energies for H₂, LiH and Li₂ computed using these correlation functions are presented, and comparison is made with previous forms, including the Jastrow-Pade form often used in Monte Carlo studies. We further treat the possibility that correlation depends not only on the separation of a pair of electrons but also on the location of the electron pair relative to the nuclei — indicative of a density-dependent or many body correlation effect. Our results indicate that such a many-body correlation effect is weakly present.This work was supported by the Director, Office of Energy Research, Office of Basic Energy Sciences, Chemical Sciences Division of the U.S. Department of Energy under Contract No. DE-AC03-76SF00098     

Correlation entropy of the H2 molecule

Correlation of a quantum many-body state makes the one-particle density matrix nonidempotent. Therefore, the Shannon entropy of the natural occupation numbers measures the correlation strength on the one-particle level. Here, it is shown how this general idea of a correlation entropy must be adapted for two-electron systems in view of conservation laws which mix Slater determinants even in the noninteracting limit. Results are presented for the correlation entropy s of H₂ as a function of the nucleus-nucleus separation R. In the ground state, the entropy of the spatial factor of the wave function maximizes 1.7 bohr beyond the Coulson-Fischer separation. The role of the correlation entropy in density functional theory is also discussed. © 1997 John Wiley & Sons, Inc.