Complete one‐ and two‐center partitioning scheme for the total energy in the Hartree–Fock theory

We proposed a complete calculation scheme for attributing the total energy by the Hartree–Fock theory to atoms (E_A) and the region between two atoms (E_AB). It was pointed out that the conventional method using the Fock matrix includes a large amount of mutual contamination in both E_A and E_AB. The new scheme was derived from the basic expression of the total energy. Calculated results by the new scheme satisfy the theoretical requirements. The scaling effect on partitioned energies was also examined. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 35–46, 1999     

Near-Dirac–Hartree–Fock results for first-row atoms calculated with GTO basis sets

Relativistic basis sets for first-row atoms have been constructed by using the near-Hartree–Fock (nonrelativistic) eigenvectors calculated by Partridge. These bases generate results of near-Dirac–Hartree–Fock quality. Relativistic total and orbital energies, relativistic corrections to the total energy, and magnetic interaction energies for the first-row atoms have been presented. The smallest Gaussian expansions (13s8 p expansions) yield Dirac–Hartree–Fock total energies accurate through six significant digits, while the largest expansions (18s13p expansions) give these energies accurate through seven significant digits. These results are more accurate than some of the results reported earlier, particularly for the open-shell atoms, indicating that the basis employed is reasonably economical for relativistic calculations. © 1995 John Wiley & Sons, Inc.     

Perturbation treatment of the Hartree–Fock equations of 1s2p3s 4Pθ state of the lithium isoelectronic sequence

The first order Hartree–Fock equations of the 1s2p3s ⁴P⁰ state of the three-electron atomic systems have been solved exactly. These solutions are used to evaluate Hartree–Fock energy up to third order with high accuracy. The third order Hartree–Fock energies for Li to Ne⁷⁺ are compared with those derived from experiment and other theoretical calculations.     

Examination of several density functionals in numerical Kohn–Sham calculations for atoms

We studied several exchange‐only and exchange–correlation energy density functionals in numerical, i.e., basis‐set‐free, nonrelativistic Kohn–Sham calculations for closed‐shell ¹S states of atoms and atomic ions with N electrons, where 2≤N≤120. Accurate total energies are presented to serve as reference data for algebraic approaches, as do the numerical Hartree–Fock results, which are also provided. Gradient‐corrected exchange‐only functionals considerably improve the total energies obtained from the usual local density approximation, when compared to the Hartree–Fock results. Such an improvement due to gradient corrections is not seen in general for highest orbital energies, neither for exchange‐only results (to be compared with Hartree–Fock results), nor for exchange–correlation results (to be compared with experimental ionization energies). © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 227–241, 2001     

Contraction of the well-tempered Gaussian basis sets: The first-row diatomic molecules

The well-tempered Gaussian basis sets (14s 10p) for atoms from lithium to neon were contracted and used in restricted Hartree–Fock calculations on 13 systems: Li₂(Σ), B₂(Σ), C₂(Σ), N₂(Σ), O₂(Σ), F₂(Σ), Ne₂(Σ), LiF(Σ), BeO(Σ), BF(Σ), CN^?(Σ), CO(Σ), and NO⁺(Σ). Spectroscopic constants (R_e, ω_e, ω_ex_e, B_e, α_e, and k_e) and one-electron properties (dipole, quadrupole, and octupole moments at the center of mass and electric field, electric field gradient, potential, and electron density at the nuclei) were evaluated and compared with the Hartree–Fock results. The largest contracted basis set (7_s6_p3_d) gives results very close to the Hartree–Fock values; the remaining differences are attributed to the absence of the f functions in the present basis sets. For Ne₂, the interaction energy was calculated; the magnitude of the basis-set superposition error was found to be very small (less than 3 μE_h at 2.8 a₀ and less than 2 μE_h at 5.0 a₀).     

Large-Z and -N dependence of atomic energies from renormalization of the large-dimension limit

By combining Hartree–Fock results for nonrelativistic ground-state energies of N-electron atoms with analytic expressions for the large-dimension limit, we have obtained a simple renormalization procedure. For neutral atoms, this yields energies typically threefold more accurate than the Hartree–Fock approximation. Here, we examine the dependence on Z and N of the renormalized energies E(N, Z) for atoms and cations over the range Z, N = 2 → 290. We find that this gives for large Z = N an expansion of the same form as the Thomas–Fermi statistical model, E → Z^7/2(C₀ + C₁Z^?1/3 + C₂Z^?2/3 + C₃Z^?3/3 + ?), with similar values of the coefficients for the three leading terms. Use of the renormalized large-D limit enables us to derive three further terms. This provides an analogous expansion for the correlation energy of the form δE δZ^4/3(δC₃ + δC₅Z^?2/3 + δC₆Z^?3/3 + ?); comparison with accurate values of δE available for the range Z ? 36 indicates the mean error is only about 10%. Oscillatory terms in E and δE are also evaluated. © 1994 John Wiley & Sons, Inc.     

Kinetics of the reactions of cyclopropane derivatives. V. The gas-phase unimolecular reactions of 1 α-chloro-2α,3α-dimethylcyclopropane and 1α-chloro-2α,3β-dimethylcyclopropane

Thermolysis of the “all-cis” compound 1α-chloro-2α,3α-dimethylcyclopropane (A) at 550–607 K and 6–115 torr is a first-order homogeneous non-radical-chain process giving penta-1,3-diene (PD) and HCl as products. The Arrhenius parameters are log₁₀A(sec^?1) = 13.92 ± 0.08 and E = 199.6 ± 0.9 kJ/mol. The isomer with trans-methyl groups, 1α-chloro-2α,3β-dimethylcyclopropane (B) reacts by two parallel first-order processes giving as observed products trans-4-chloropent-2-ene (4CP) and PD + HCl, with log₁₀A(sec^?1) = 14.6 and 13.8, respectively, and E = 199.5 and 190.2 kJ/mol, respectively. The 4CP undergoes secondary decomposition to PD + HCl (as investigated previously). Comparison of the results for compounds (A) and (B) with those for other gas-phase and solution reactions leads to the conclusion that the gas-phase thermolyses proceed by rate-determining ring opening to form olefins which may decompose further by thermal or chemically activated reactions, and that the ring opening is a semiionic electrocyclic reaction in which alkyl groups in the 2,3-positions trans to the migrating chlorine semianion move apart, with appropriate consequences for the rate of reaction and the stereochemistry of the products.     

Coulomb-Hole–Hartree–Fock functional

We have extended to molecules a density functional previously parametrized for atomic computations. The Coulomb-hole–Hartree–Fock functional, introduced by Clementi in 1963, estimates the dynamical correlation energy by the computations of a Hartree–Fock-type single-determinant wave function, where the Hartree–Fock potential was augmented with an effective potential term, related to a hard Coulomb hole enclosing each electron. The method was later revisited by S. Chakravorty and E. Clementi [Phys. Rev. A 39 , 2290 (1989)], where a Yukawa-type soft Coulomb hole replaced the previous hard hole; atomic correlation energies, computed for atoms with Z = 2 to Z = 54 as well as for a number of excited states, validated the method. In this article, we parametrized a function, which controls the width of the soft Coulomb hole, by fitting the first and second atomic ionization potentials of the atoms with 1 ? Z ? 18. The parametrization has been preliminarily validated by computing the dissociation energy for a number of molecules. A few-determinant version of the Coulomb-hole–Hartree–Fock method, necessary to account for the nondynamic correlation corrections, is briefly discussed. © 1994 John Wiley & Sons, Inc.     

Variation of parameters in Becke‐3 hybrid exchange‐correlation functional

We have investigated the consequences of varying the three parameters in Becke's hybrid exchange‐correlation functional, which includes five contributions: Hartree–Fock exchange, local exchange, Becke's gradient exchange correction, local correlation, and some form of gradient correlation correction. Our primary focus was upon obtaining orbital energies with magnitudes that are reasonable approximations to the electronic ionization potentials; however, we also looked at the effects on molecular geometries and atomization enthalpies. A total of 12 parameter combinations was considered for each of three different gradient correlation corrections: the Lee–Yang–Parr, the Perdew‐86, and the Perdew–Wang 91. Five molecules were included in the study: HCN, N₂, N₂O, F₂O, and H₂O. For comparison, a Hartree–Fock calculation was also carried out for each of these. The 6‐31+G** basis set was used throughout this work. We found that the ionization potential estimates can be greatly improved (to much better than Hartree–Fock levels) by increasing the Hartree–Fock exchange contribution at the expense of local exchange. In itself, this also introduces major errors in the atomization enthalpies. However, this can be largely or even completely counteracted by reducing or eliminating the role of the gradient exchange correction. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 227–238, 2000     

A generalized formula for the energies of alternant molecular orbitals. II. Heteronuclear molecules

Using the method of alternant molecular orbitals (AMO ), it is shown that the energies of AMOS (E_kσ) for an arbitrary heteronuclear alternant system, having a singlet ground state, are connected with the energies of MOS (e_{k(k )}) obtained by means of the conventional Hartree–Fock (HF ) method (SCF -LCAO -MO -PPP ) via the formula: In the general case, the determination of the correlation corrections δ_i,kσ is connected with the solving of a complicated system of integral equations, which is considerably simplified if the Hubbard approximation is accepted for the electron interaction. The energy spectrum of a chain with two atoms in the elementary cell (AB)_n is considered as an example. It is shown that if nontrivial solutions exist (δ_i,kσ ≠ 0), the correlation correction for AMOS of different spin are different (δ_i,kσ ≠ δ_i,kβ), from which it follows, that the width of the energy gap ΔE_∞ for AMOS with different spin is different: ΔE_∞,α ≠ ΔE_∞,β.     

Adapted Gaussian basis sets for atoms Cs to Lr based on the generator coordinate Hartree–Fock method

We have applied a discretized version of the generator coordinate Hartree–Fock method to generate adapted Gaussian basis sets for atoms Cs (Z=55) to Lr (Z=103). Our Hartree–Fock total energy results, for all atoms studied, are better than the corresponding Hartree–Fock energy results attained with previous Gaussian basis sets. For the atoms Cs to Lr we have obtained an energy value within the accuracy of 10⁻⁴ to 10⁻³ hartree when compared with the corresponding numerical Hartree–Fock total energy results. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 858–865, 1998     

Absorption spectrum of manganian andalusite: Cluster calculation by an ab initio method

The ground state and the first few excited states of an MnO₆^9? cluster are calculated in the unrestricted Hartree–Fock model. The state ordering is ⁵B_{1 g}, ⁵A_{1 g}, ⁵B_{2 g}, and ⁵E_g as can be expected from simpler models. Consistent with the results by the same method for copper complexes, we obtain d–d transition energies about one half or less of the experimental energies. The charge transfer spectrum is subject to a large spin polarization in the sense that the lowest charge transfer state (⁵E_u) has five unpaired spins on Mn.     

Ab initio studies of F−(H2O)n and Cl−(H2O)n clusters for n = 1, 2

Ab initio Hartree–Fock calculations are performed on hydrates of the F^? and Cl^? ions using 6-31G, 6-31G**, and 6-21G basis sets. Geometries and binding energies are obtained. An estimate of the correlation energy is provided by an MP2/6-31G (Møller-Plesset second-order perturbation) calculation. Comparisons are made between the Cl^?(SO₂) and the Cl^?(H₂O) complexes.     

Core polarization and relativistic effects competition in the first ionization potentials for some systems in the Cu,Ag, and Au isoelectronic sequences

Single-configuration relativistic Hartree–Fock values of the first ionization potentials for Cu through Kr⁷⁺, Ag through I⁶⁺, and Au through Pb³⁺ are computed in “frozen” and “relaxed core” approximations with and without allowance for core polarization. Effects of polarization of the atomic core by the valence electron are included by introducing a polarization potential in the one-electron Hamiltonian of the valence electron. The core polarization potential depends on two parameters, the static dipole polarizability of the core α and the cut-off radius r₀, which are chosen independently of the ionization potential data. It is demonstrated that by including the core polarization potential with α and r₀ parameters, which are simply chosen instead of being empirically fitted, it is still possible to account, on the average, for at least 70% of the discrepancy between the single-configuration relativistic Hartree–Fock ionization potentials and the experiment, a discrepancy usually ascribed to the contribution of valence-core electron correlations, and to bring the theoretical ionization potentials to an average agreement with experiment of around 1%. It can be concluded from this study that for low and medium Z elements the core polarization dominates for neutral systems or systems in low ionization stages, whereas for highly ionized systems the relativistic effects prevail. For heavy elements, however, the core polarization influence is comparable to the relativistic one only for neutral systems, whereas for ions the relativistic effects are overwhelmingly predominant.     

Ab initio Hartree–Fock and configuration-interaction treatment of the interaction between two nickel atoms

The interaction between two nickel atoms in the configurations (3d)⁸(4s)² and (3d)⁹ (4s)¹ has been calculated using ab initio methods (Hartree–Fock and configuration interaction). The results of the calculations compare favorably with the optical spectrum. The discrepancy between the calculated and the experimental dissociation energy is discussed, and a new estimate of the dissociation energy is given. The configuration-interaction calculations show that the interaction between the two nickel atoms is of a very complex nature. In spite of this the binding can be interpreted in a simple way. The bond is minly due to the 4sσ_g molecular orbital while the 3d orbitals of the two nuclei are exchange coupled.     

Radical-initiated homo- and copolymerization of α-methylene-γ-butyrolactone

The kinetics of α-methylene-γ-butyrolactone (α-MBL) homopolymerization was investigated in N,N-dimethylformamide (DMF) with azobis(isobutyronitrile) as initiator. The rate of polymerization (R_p) was expresed by R_p = k[AIBN]^0.54[α-MBL]^1.1 and the overall activation energy was calculated as 76.1 kJ/mol. Kinetic constants for α-MBL polymerization were obtained as follows: k_p/k_t^1/2 = 0.161 L^1/2 mol^?1/2·s^?1/2; 2fk_d = 2.18 × 10^?5 s^?1. The relative reactivity ratios of α-MBL(M₂) copolymerization with styrene (r₁ = 0.14, r₂ = 0.87) were obtained. Applying the Q–e scheme led to Q = 2.2 and e = 0.65. These Q and e values for α-MBL are higher than those for MMA     

Mξ, Mαβ, Mγ and Mm X-ray production cross-sections for elements with 71⩽z⩽92 at 5.96 keV photon energy

The M_ξ, M_αβ, M_γ and M_m X-ray production (XRP) cross-sections have been measured for the elements with 71⩽Z⩽92 at 5.96 keV incident photon energy satisfying E_M1<E_inc<E_L3, where E_M1(L3) is the M₁(L₃) subshell binding energy. These XRP cross-sections have been calculated using photoionization cross-sections based on the relativistic Dirac–Hartree–Slater (RDHS) model with three sets of X-ray emission rates, fluorescence, Coster–Kronig and super Coster–Kronig yields based on (i) the non-relativistic Hartree–Slater (NRHS) potential model, (ii) the RDHS model and (iii) the relativistic Dirac–Fock (RDF) model. For the third set, the M_i (i=1–5) subshell fluorescence yields have been calculated using the RDF model-based X-ray emission rates and total widths reevaluated to incorporate the RDF model-based radiative widths. The measured cross-sections have been compared with the calculated values to check the applicability of the physical parameters based on different models.     

Comprehensive understanding of size‐, shape‐, and composition‐dependent polarizabilities of SimCn (m,n = 1–4) clusters

We performed a comprehensive study of the size‐, shape‐, and composition‐dependent polarizabilities of Si_mC_n (m, n = 1–4) clusters on the basis of the density‐functional‐based coupled perturbed Hartree–Fock calculations. We found better correlations between the polarizabilities and both the binding energies (E_b) and change in charge distribution (Δq) than the energy gaps. The α values exhibit overall decreasing and increasing trends with increases in the E_b and Δq values, respectively. For isomers with the same E_b values and different polarizabilities, Δq can well explain the difference in polarizabilities. The π‐electron delocalization effect is the best factor for understanding the shape‐dependence. For a given m/n value, the linear clusters have an obviously larger polarizability than both the prolate and compact clusters, irrespective of the cluster size. We fit a quantitative expression [α = A ? (A ? B) × exp(?k(m/n))] to describe the composition‐dependent polarizabilities. © 2012 Wiley Periodicals, Inc.     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000