Self‐consistent shielding in atoms and molecules

The orbital exponents of trial wave functions for simple systems can be found from the potential energy terms alone. Shielding of the nuclear charge by one electron on another is determined by the relative values of the nuclear–electron attraction and the electron–electron repulsion. For two electrons in the same orbital, the shielding is divided equally. For different orbitals, only the inner electron shields the outer. The systems tested are first‐row atoms, using Slater orbitals. It appears that if this approach can be generalized, it may not be necessary to calculate kinetic energies in chemical systems, since they will be determined by the orbital exponents. This would be useful if trial wave functions were not available, but trial electron density functions were. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ab initio MO studies of nuclear spin-spin coupling constants in CH4, SiH4, AlH 4 − and GeH4 systems

Summary Ab initio molecular orbital calculations of electron coupled nuclear spin-spin coupling constants are performed for CH₄, SiH₄, AlH ₄ ^– and GeH₄ systems using the SCF perturbation theory. Basis set dependence of the major contributing terms such as orbital diamagnetic, orbital paramagnetic, spin dipolar and Fermi contact terms are studied. The study also illustrates the relative importance of bond centred functions and nuclear centred polarization functions in predicting the directly bonded and geminal couplings in the systems selected. Basis sets having uncontracted cores functions and augmented with bond functions seem to predict most of these couplings fairly satisfactorily when compared to the experimental values.     

Local potentials in independent‐electron models

Density functional theory, exact for ground states, is commonly assumed to imply an independent‐electron model in which only local potential functions appear. It has recently been shown that several paradoxes in different aspects of the theory can be resolved if this locality hypothesis is abandoned. However, the locality hypothesis itself appears to be implied by rigorous variational theory. This conflict is discussed and resolved here. The resolution involves embedding the density functional theory in an orbital functional theory whose functional derivatives are not confined to normalized ground states. The orbital Euler–Lagrange equations of this extended theory in general contain effective potentials that are linear operators acting on orbital wave functions. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 384–388, 2001     

Orbital symmetry,orbital stability,and orbital pairing rules for organic reactions in the ground state

A generalization of the Hartree–Fock molecular orbital (MO) theory for treating diradical intermediates was explained pictorially by drawing molecular orbitals of diradical species such as ring-opened trimethylene. The generalized MO theory applied to elucidate electronic mechanisms of concerted, ionic, radical, and ion-radical reactions of organic reactants in the ground state. Generalized MO computations revealed the most essential characteristics of these reactions and mutal relationships between the worlds of Woodward–Hoffmann and Hughes–Ingold. Generalized MO studies supported our orbital symmetry, stability and pairing rules for concerted, ionic and radical reactions in the ground state, respectively. An extension of MO treatments to excited states reactions was briefly pointed out in relation to the density and spin correlation functions by the multireference CI wave functions.     

Analytical Hartree–Fock wave functions subject to cusp and asymptotic constraints: He to Xe,Li+ to Cs+, H− to I−

Analytical, variational approximations to Hartree–Fock wave functions are constructed for the ground states of all the neutral atoms from He to Xe, the cations from Li⁺ to Cs⁺, and the stable anions from H⁻ to I⁻. The wave functions are constrained so that each atomic orbital agrees well with the electron–nuclear cusp condition and has good long‐range behavior. Painstaking optimization of the exponents and principal quantum numbers of the Slater‐type basis functions allows us to reach this goal while obtaining total energies that, at worst, are a few microHartrees above the numerical Hartree–Fock limit values. The wave functions are freely available by anonymous ftp from okapi.chem.unb.ca or upon request to the authors. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 491–497, 1999     

An efficient algorithm for complete active space valence bond self‐consistent field calculation

This article describes a novel algorithm for the optimization of valence bond self‐consistent field (VBSCF) wave function for a complete active space (CAS), so‐called VBSCF(CAS). This was achieved by applying the strategies adopted in the optimization of CASSCF wave functions to VBSCF(CAS) wave functions, using an auxiliary orthogonal orbital set that generates the same configuration space as the original nonorthogonal orbital set. Theoretical analyses and test calculations show that the VBSCF(CAS) method shares the same computational scaling as CASSCF. The test calculations show the current capability of VBSCF method, which involves millions of VB structures. © 2012 Wiley Periodicals, Inc.     

Molecular MC–SCF calculations

A practical method for finding multi-configurational SCF wave functions is proposed. The basic equation is equivalent to the Brillouin theorem; comparison with the usual SCF equations obtained through effective hamiltonians gives an interpretation of the offdiagonal Lagrange multipliers. Numerical applications to Formaldehyde in a minimum Slater-type orbital basis with four different variational wave functions are reported. The molecular orbitals found in these calculations are localized on the chemical bonds. The largest contributions to the energy are obtained from π-π and dispersion-type σ-π correlation.     

One-Parameter double-zeta atomic functions for the hartree–fock energies of helium isoelectronic sequence

The double-zeta atomic functions are characterized by the nuclear charge z of the two-electron atomic system. The Hartree–Fock total energies and the corresponding orbital energies are calculated using various atomic wave functions for the helium isoelectronic sequence. The expectation values rⁿ of various wave functions are also examined. It is found that the accuracy of our one-parameter double-zeta functions corresponds to the accuracy of the usual five-parameter double-zeta functions.     

Matrix elements for Ĵ2 and Ĵz operators over explicitly correlated Cartesian Gaussian functions

General formalism for evaluation of multiparticle integrals involving J?² and J?_z operators over explicitly correlated Cartesian Gaussian functions is presented. The integrals are expressed in terms of the general overlap integrals. An explicitly correlated Cartesian Gaussian function is a product of spherical orbital Gaussian functions, powers of the Cartesian coordinates of the particle, and exponential Gaussian factors, which depend on interparticular distances. This development is relevant to both adiabatic and nonadiabatic calculations of energy and properties of multiparticle systems. © 1995 John Wiley & Sons, Inc.     

Density Functional Theory Investigations on Vibrational Spectra,Molecular Structure,and Properties of the L‐Serine,L‐Cysteine,and L‐Aspartic Acid Molecules

In this paper, we present a thorough investigation of the conformational space to characterize all possible gas‐phase structures of the neutral L‐serine, L‐cysteine, and L‐aspartic acid molecules. A total of 120 trial structures were generated for L‐aspartic acid and 96 trial structures for L‐serine and L‐cysteine by combining all internal single‐bond rotamers. Various combinations of the Hartree–Fock and density functional theory/B3LYP methods with different bases were used to optimize all possible trial structures. The theoretical studies on the structure, harmonic vibrational spectra, and molecular properties of these amino acids are presented. The assignments of the calculated wave numbers resulting from potential energy distributions were performed using the VEDA 4 program to allow a good interpretation of the theoretical vibrational spectra of the title compounds. The fundamental harmonic frequencies were found to be in good agreement with data in the literature. A natural bond orbital analysis was performed to investigate the charge delocalization throughout the molecules for the three test compounds. Moreover, an extensive discussion of the highest occupied molecular orbital–lowest unoccupied molecular orbital energy gap as well as other related molecular properties are reported.     

Elimination, in electronic structure calculations, of redundant orbital products

We propose a direct method for reducing the dimension of the space of orbital products that occur, for example, in the calculation of time dependent density functional theory linear response and in Hedin's GW approximation to the electron propagator. We do this by defining, within the linear space of orbital products, a subspace of dominant directions that are associated with a certain eigenvalue problem. These directions span the entire linear space of products with an error that decreases approximately exponentially with their number. Our procedure works best for atomic orbitals of finite range and it avoids the use of extra sets of auxiliary fit functions.     

A general,recursive, and open‐ended response code

We present a new implementation of a recent open‐ended response theory formulation for time‐ and perturbation‐dependent basis sets (Thorvaldsen et al., J. Chem. Phys. 2008, 129, 214108) at the Hartree–Fock and density functional levels of theory. A novel feature of the new implementation is the use of recursive programming techniques, making it possible to write highly compact code for the analytic calculation of any response property at any valid choice of rule for the order of perturbation at which to include perturbed density matrices. The formalism is expressed in terms of the density matrix in the atomic orbital basis, allowing the recursive scheme presented here to be used in linear‐scaling formulations of response theory as well as with two‐ and four‐component relativistic wave functions. To demonstrate the new code, we present calculations of the third geometrical derivatives of the frequency‐dependent second hyperpolarizability for HSOH at the Hartree–Fock level of theory, a seventh‐order energy derivative involving basis sets that are both time and perturbation dependent. © 2014 Wiley Periodicals, Inc.     

Configuration interaction by the method of bonded functions

An approach is described to construct a complete and unique list of bonded functions for arbitrary spin states from any number of reference functions, any degree of orbital substitution, and any orbital pool.     

Analytical optimization of exponent values in protonic and deuteronic Gaussian‐type functions by elimination of translational and rotational motions from multi‐component molecular orbital scheme

To optimize the exponent values in protonic and deuteronic Gaussian‐type functions (GTF) by the elimination of translational and rotational motions, we have proposed the new scheme of an analytical gradient formula with respect to the exponent values in the multi‐component molecular orbital scheme, which can take into account the quantum effects of protons and deuterons, under the Hartree‐Fock level of theory. Numerical assessment of H₂ and D₂ molecules confirms that there is a clear difference between distributions of protonic and deuteronic orbitals following the elimination of translational and rotational motions. In particular, the d‐type GTF in the protonic orbital drastically improves the total energy. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

A theoretical investigation of the excited states of OCLO radical, cation, and anion using the CASSCF/CASPT2 method

Using the complete active space self-consistent field method with a large atomic natural orbital basis set, 10, 13, and 9 electronic states of the OClO radical, OClO(+) cation, and OClO(-) anion were calculated, respectively. Taking the further correlation effects into account, the second-order perturbation (CASPT2) calculations were carried out for the energetic calibration. The photoelectron spectroscopy of the OClO radical and OClO(-) anion were extensively studied in the both case of the adiabatic and vertical ionization energies. The calculated results presented the relatively complete assignment of the photoelectron bands of the experiments for OClO and its anion. Furthermore, the Rydberg states of the OClO radical were investigated by using multiconfigurational CASPT2 (MS-CASPT2) theory under the basis set of large atomic natural orbital functions augmented with an adapted 1s1p1d Rydberg functions that have specially been built for this study. Sixteen Rydberg states were obtained and the results were consistent with the experimental results.     

Solution of the Hartree–Fock problem by expansion onto nested bases

A method for solving the Hartree–Fock problem in a finite basis set is derived, which permits each orbital to be expanded in a different basis. If the basis set for each orbital ?_i contains the basis functions for the preceding orbitals, ?_i?1, ?_i?2,… ?₁, then the ?_i form an orthonormal set. One advantage over the standard Hartree–Fock method is that a different long range behavior for each orbital, as for example is required in the Hartree–Fock-Slater method, can be forced. A calculation on the ground state of beryllium is performed using the nested procedure. Very little energy is lost because of nesting, and the node in the 1s orbital disappears.     

Photoionization cross sections of H2(+) and H2 with complex Gaussian-type basis functions optimized for the frequency-dependent polarizabilities

Within the framework of the complex basis function method, the photoionization cross sections of H(2)(+) and H(2) were calculated based on the variational principle for the frequency-dependent polarizabilities. In these calculations, complex orbital exponents of Gaussian-type basis functions for the final state continuum wavefunctions were fully optimized for each photon energy with the numerical Newton-Raphson method. In most cases, the use of only one or two complex Gaussian-type basis functions was enough to obtain excellent agreement with previous high precision calculations and available experimental results. However, there were a few cases, in which the use of complex basis functions having various angular momentum quantum numbers was crucial to obtain the accurate results. The behavior of the complex orbital exponents as a function of photon energy was discussed in relation to the scaling relation and the effective charge for photoelectron. The success of this method implies the effectiveness of the optimization of orbital exponents to reduce the number of basis functions and shows the possibility to calculate photoionization cross sections of general molecules using only Gaussian-type basis functions.     

Relative Stability of closo‐closo,closo‐nido,and nido‐nido Macropolyhedral Boranes: The Role of Orbital Compatibility

The concept of orbital compatibility is used to explain the relative energies of different macropolyhedral structural patterns such as closo‐closo, closo‐nido, and nido‐nido. A large polyhedral borane condenses preferentially with a smaller polyhedron owing to orbital compatibility. Calculations carried out at the B3LYP/6‐31G* level show that the macropolyhedron closo(12)‐closo(6) is the most preferred structural pattern among the face‐sharing closo‐closo systems. The relative stabilities of four‐shared‐atom closo‐closo, three‐shared‐atom closo‐closo, three‐shared‐atom closo‐nido, edge‐sharing closo‐nido, and edge‐sharing nido‐nido structures are in accordance with the difference in the number of vertices of the individual polyhedra of the macropolyhedra. When the difference in the number of vertices of the individual polyhedra is large, the stability of the macropolyhedra is also large. Calculations further show that the orbital compatibility plays an important role in deciding the stability of the macropolyhedral boranes with more than two polyhedral units. The dependence of the orbital compatibility on the relative stability of the macropolyhedron varies with other factors such as inherent stability of the individual polyhedron and steric factors.     

HF and MP2 calculations on CN−, N2, AlF,SiO, PN,SC, ClB,and P2 using correlated molecular wave functions

Contracted basis sets of double zeta valence quality plus polarization functions (DZP) and augmented DZP basis sets, which were recently constructed for the first‐ and second‐row atoms, are applied to study the electronic ground states of the diatomic molecules CN^?, N₂, AlF, SiO, PN, SC, ClB, and P₂. At the Hartree–Fock (HF) and/or Møller–Plesset second‐order (MP2) levels, total and molecular orbital energies, dissociation energies, bond lengths, harmonic vibrational frequencies, and dipole moments are calculated and compared with available experimental data and with the results obtained from correlation consistent polarized valence basis sets of Dunning's group. For N₂, calculations of polarizabilities at the HF and MP2 levels with the sets presented above are also done and compared with results reported in the literature. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006