An approximate frozen core approach to valence only molecular calculations

A simple valence electron-only theory based on an approximate frozen core approach and an exact core-valence strong orthogonality condition is developed for atomic and molecular systems. A unique reduced basis is introduced in which both core and valence orbitals are expanded. The core representation is roughly approximated, and the valence orbital overlap with the corresponding all-electron reference functions is nearly exact. The size of the reduced basis in terms of primitive functions is practically the same as that adopted by effective core potential methods in which the valence orbitals have the correct nodal properties. Results obtained with the present approach are presented for LiO, BeO and CaO molecules, and compared with the corresponding all-electron frozen core calculations. In addition, a detailed investigation on Li _n Be clusters (n=1,..., 6) is carried out.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

Medium-size Gaussian basis sets for hydrogen through argon

Summary Medium-sized Gaussian basis sets are reoptimized for the ground states of the atoms from hydrogen through argon. The composition of these basis sets is (4s), (5s), and (6s) for H and He, (9s5p) and (12s7p) for the atoms Li to Ne, and (12s8p) and (12s9p) for the atoms Na to Ar. Basis sets for the² P states of Li and Na, and the³ P states of Be and Mg are also constructed since they are useful in molecular calculations. In all cases, our energies are lower than those obtained previously with Gaussian basis sets of the same size.     

Basis Set Orbital Relaxation in Atomic and Molecular Hydrogen Systems

Orbital relaxation (OR) amounts to variation of the orbital exponents in hydrogen molecules and ions relative to the exponents of the isolated atom; it is represented as the sum of the one- and two-center contributions depending on the effective atomic charge and on the presence of other atoms in the molecule. The procedure for isolating the contributions of the exponent includes treatment of the OR of hydrogen in a special set of neutral and charged atoms and molecules with certain multiplicities of their electronic states. Within the framework of the spin-unrestricted Hartree-Fock method, we found and discussed the optimal values of the exponents of the basis orbitals of hydrogen atoms and molecules using the minimal split valence-shell basis set, the basis set that includes the polarization function, and the expanded set of grouped natural orbitals. A simple energy model is suggested for OR. Expressions are derived for evaluating the exponents of the relaxed orbitals in hydrogen-containing systems.Original Russian Text Copyright © 2004 by A. I. Ermakov, A. E. Merkulov, A. A. Svechnikova, and V. V. Belousov__________Translated from Zhurnal Strukturnoi Khimii, Vol. 45, No. 6, pp. 973–978, November–December, 2004.     

On the use of local basis sets for localized molecular orbitals

Two procedures are discussed for the direct variational optimization of localized molecular orbitals which are expanded in local subsets of the molecular basis set. It is shown that a Newton-Raphson approach is more efficient than an iterative diagonalization scheme. The effect of the basis-set truncation on the quality ofab-initio SCF results is investigated for Be, Li₂, HF, H₂O, NH₃, CH₄ and C₂H₆.     

Analogy between real irreducible tensor operator and molecular two-particle operator

. Molecular matrix elements of a physical operator are expanded in terms of polycentric matrix elements in the atomic basis by multiplying each by a geometrical factor. The number of terms in the expansion can be minimized by using molecular symmetry. We have shown that irreducible tensor operators can be used to imitate the actual physical operators. The matrix elements of irreducible tensor operators are easily computed by choosing rational irreducible tensor operators and irreducible bases. A set of geometrical factors generated from the expansion of the matrix elements of irreducible tensor operator can be transferred to the expansion of the matrix elements of the physical operator to compute the molecular matrix elements of the physical operator. Two scalar product operators are employed to simulate molecular two-particle operators. Thus two equivalent approaches to generating the geometrical factors are provided, where real irreducible tensor sets with real bases are used. Received: 3 September 1996 / Accepted: 19 December 1996     

Complete graph conjecture for inner-core electrons: homogeneous index case

The complete graph conjecture that encodes the inner-core electrons of atoms with principal quantum number n >or= 2 with complete graphs, and especially with odd complete graphs, is discussed. This conjecture is used to derive new values for the molecular connectivity and pseudoconnectivity basis indices of hydrogen-suppressed chemical pseudographs. For atoms with n = 2 the new values derived with this conjecture are coincident with the old ones. The modeling ability of the new homogeneous basis indices, and of the higher-order terms, is tested and compared with previous modeling studies, which are centered on basis indices that are either based on quantum concepts or partially based on this new conjecture for the inner-core electrons. Two similar algorithms have been proposed with this conjecture, and they parallel the two "quantum" algorithms put forward by molecular connectivity for atoms with n > 2. Nine properties of five classes of compounds have been tested: the molecular polarizabilities of a class of organic compounds, the dipole moment, molar refraction, boiling points, ionization energies, and parachor of a series of halomethanes, the lattice enthalpy of metal halides, the rates of hydrogen abstraction of chlorofluorocarbons, and the pED(50) of phenylalkylamines. The two tested algorithms based on the odd complete graph conjecture give rise to a highly interesting model of the nine properties, and three of them can even be modeled by the same set of basis indices. Interesting is the role of some basis indices all along the model.     

Complete basis set ab initio computational study of ionization potential, electron affinity and the C-F bond dissociation energy for perfluorinated methane derivatives

A computational study of perfluorinated methane derivatives was performed with complete basis set ab initio methods. The total energies for their neutral, cation, and anionic states were computed. From these values, the energy gaps between different electronic states, ionization potentials, electron affinities, and C-F bond dissociation energies were calculated. The computed values are compared with experimental data and the reliability of complete basis set ab initio methods is discussed. New values for C-F bond dissociation energies are suggested. Received: 12 January 1998 / Accepted: 2 April 1998 / Published online: 29 July 1998     

Medium-size polarized basis sets for high-level-correlated calculations of molecular electric properties

Summary The first-order polarized basis sets for the use in high-level-correlated investigations of molecular electric properties have been generated for Pb, Bi, Po, and At. The performance of the standard [10.17.14.5/13.11.8.2] and extended [20.17.14.9/13.11.8.4] basis sets has been examined in nonrelativistic and quasirelativistic calculations for atoms and simple closed-shell hydrides. The relativistic contributions to electric dipole properties of those systems have been evaluated by using the recently developed quasirelativistic scheme. The predicted dipole polarizability of At is in good agreement with the results of other relativistic calculations. The calculated quasirelativistic dipole moments of BiH₃ (–0.499 a.u.), PoH₂ (–0.207 a.u.), and AtH (+0.036 a.u.) involve a significant relativistic contribution which amounts to —0.230 a.u., –0.177 a.u., and –0.097 a.u., respectively. The basis set details append this paper. They are also available as a part of the basis set library of the MOLCAS system.     

Ab initio and DFT studies on hydrolyses of phosphorus halides

Hydrolyses of phosphorus halides, (RO)(2)POX where R = H or Me and X = F or Cl, in the gas phase and in the reaction field have been investigated theoretically with ab initio and the density functional theory (DFT). The free energy of activation in the reaction field was also estimated using the Onsager method with a correction of entropy change and basis set superposition error (BSSE). The reaction of (MeO)(2)POF proceeds through a path with bifunctional catalysis regardless of the medium, but the reaction of (MeO)(2)POCl proceeds through bifunctional and general base catalysis in the gas phase and in water, respectively. The estimated free energy barrier of 23 kcal/mol for the hydrolysis of (MeO)(2)POF is in good agreement with the experimental values of 24 kcal/mol, and relative barrier of 3 kcal/mol to the (MeO)(2)POCl is also in good agreement with the experimental values of 5 kcal/mol of diisopropyl phosphorus halides ((Pr(i)O)(2)POX, X = F and Cl).     

X-H (X = C, N, O, Si, P, S) 键离解能计算的密度泛函方法研究

使用了不同密度泛函方法计算X-H (X = C, N, O, Si, P, S) 键离解能，并分析不同密度泛函方法的计算精度。研究发现大多数密度泛函方法包括B3LYP, B3P86, B3PW91, G96LYP, PBE1PBE，和BH&HLYP都明显低估键离解能13-25 kJ/mol。该现象与是否使用无限基组无关，因为即使使用无限基组键离解能仍然被低估。因此密度泛函方法不适合用于键离解能的估算。其中B3P86方法的偏差最小。进一步分析表明，使用限制性开壳层计算并无任何优势，在大多数情况下非限制性开壳层计算实际上比限制性开壳层计算要好。最后，我们发现了密度泛函方法对键离解能的低估是系统的，因此建议利用校准后的UDFT/6-311++G(d, p)方法计算化学键离解能。