Spin-Dependent operators in the spin-free quantum chemistry

The two-particle spatial density matrix components introduced by McWeeny are expressed in terms of the Fock coordinate wave function, which is constructed from an arbitrary function of N spatial coordinates. The integral relations for these components are verified. The necessary matrix elements of a standard representation of the S_N group are calculated.     

Low-cost evaluation of the exchange Fock matrix from Cholesky and density fitting representations of the electron repulsion integrals

The authors propose a new algorithm, "local K" (LK), for fast evaluation of the exchange Fock matrix in case the Cholesky decomposition of the electron repulsion integrals is used. The novelty lies in the fact that rigorous upper bounds to the contribution from each occupied orbital to the exchange Fock matrix are employed. By formulating these inequalities in terms of localized orbitals, the scaling of computing the exchange Fock matrix is reduced from quartic to quadratic with only negligible prescreening overhead and strict error control. Compared to the unscreened Cholesky algorithm, the computational saving is substantial for systems of medium and large sizes. By virtue of its general formulation, the LK algorithm can be used also within the class of methods that employ auxiliary basis set expansions for representing the electron repulsion integrals.     

Density matrix of crystalline systems. II. The influence of structural and computational parameters

The one-electron density matrix (DM ) of a number of crystalline systems—lithium, graphite, boron nitride, silicon, and beryllium—are considered here, as resulting from Hartree–Fock–SCF-LCAO calculations. The influence of structural and computational parameters is discussed. It is shown in particular why the structure of chemical bonds in semiconductors leads to an oscillating long-range behavior of the DM , similar to that observed in metals, where these oscillations are related to the very existence of a Fermi surface. Concerning computational parameters, the influence of the density of sampling k points and of basis set on the calculated DM is considered; it is shown that the choice of the basis set is not a very critical one as far as the DM range is concerned. Some critical aspects of the interrelation between DM range and exchange part of the Fock Hamiltonian are analyzed.     

A novel parallel algorithm for large-scale Fock matrix construction with small locally distributed memory architectures: RT parallel algorithm

We developed a novel parallel algorithm for large-scale Fock matrix calculation with small locally distributed memory architectures, and named it the "RT parallel algorithm." The RT parallel algorithm actively involves the concept of integral screening, which is indispensable for reduction of computing times with large-scale biological molecules. The primary characteristic of this algorithm is parallel efficiency, which is achieved by well-balanced reduction of both communicating and computing volume. Only the density matrix data necessary for Fock matrix calculations are communicated, and the data once communicated are reutilized for calculations as many times as possible. The RT parallel algorithm is a scalable method because required memory volume does not depend on the number of basis functions. This algorithm automatically includes a partial summing technique that is indispensable for maintaining computing accuracy, and can also include some conventional methods to reduce calculation times. In our analysis, the RT parallel algorithm had better performance than other methods for massively parallel processors. The RT parallel algorithm is most suitable for massively parallel and distributed Fock matrix calculations for large-scale biological molecules with more than thousands of basis functions.     

A comparison of Hartree—Fock,MP2, and DFT results for the HCN dimer and crystal

A number of hydrogen-bond related quantities—geometries, interaction energies, dipole moments, dipole moment derivatives, and harmonic vibrational frequencies—were calculated at the Hartree—Fock, MP2, and different DFT levels for the HCN dimer and the periodic HCN crystal. The crystal calculations were performed with the Hartree—Fock program CRYSTAL92, which routinely allows an a posteriori electron-correlation correction of the Hartree—Fock obtained lattice energy using different correlation-only functionals. Here, we have gone beyond this procedure by also calculating the electron-correlation energy correction during the structure optimization, i.e., after each CRYSTAL92 Hartree—Fock energy evaluation, the a posteriori density functional scheme was applied. In a similar manner, we optimized the crystal structure at the MP2 level, i.e., for each Hartree—Fock CRYSTAL92 energy evaluation, an MP2 correction was performed by summing the MP2 pair contributions from all HCN molecules within a specified cutoff distance. The crystal cell parameters are best reproduced at the Hartree—Fock and the nongradient-corrected HF + LDA and HF + VWN levels. The BSSE-corrected MP2 method and the HF + P91, HF + LDA, and HF + VWN methods give lattice energies in close agreement with the ZPE-corrected experimental lattice energy. The (HCN)₂ dimer properties are best reproduced at the MP2 level, at the gradient-corrected DFT levels, and with the B3LYP and BHHLYP methods. © 1996 John Wiley & Sons, Inc.     

Calculation of molecular polarizabilities for large systems within the random phase approximation

We examine the static-field molecular polarizability from a sum over uncoupled Hartree—Fock states (SOS), the Tamm—Dancoff approximation (TDA), and the random phase approximation (RPA). An efficient algorithm for the inversion of the TDA or RPA matrix is outlined, which avoids matrix diagonalization and explicit construction of matrix elements over states, allowing for rapid calculation of molecular polarizabilities. The extension of the method is straightforward; third-order hyperpolarizability is developed as an example. Test cases are reported for molecules represented by an intermediate neglect of differential overlap (INDO) wavefunction.     

Ab initio studies on hydrogen bonded chains. I. Equilibrium geometry of the infinite,linear chain of hydrogen fluoride molecules

The idealized case of an infinite, linear chain of hydrogen fluoride molecules is studied at the Hartree—Fock level with the aid of the crystal orbital method. Extended gaussian basis sets have been used to compute the equilibrium structure and the stabilization energy (hydrogen bond energy) per HF molecule. It is demonstrated that near Hartree—Fock limit results for this model system account for a large part of the observed differences between isolated dimers in the gas phase and the infinite periodic crystal. For the infinite chain the following results were obtained: r_HF = 1.721 bohr, r_FF = 5.049 bohr and ΔE (hydrogen bond energy per HF) = 5.9 kcal/mole.     

Comparative study of the Hartree–Fock and the Hartree–Fock–Slater method

The Hartree–Fock method (standard Roothaan closed-shell HF –LCAO theory) and the Hartree–Fock–Slater method (restricted HFS –LCAO –DV method developed by Baerends and Ros) have been compared with emphasis on the respective one-electron equations and on the matrix elements of the respective Fock operators. Using the same STO basis in the two cases, the matrix elements of the Fock operators and of their separate one-electron, Coulomb, and exchange contributions have been calculated for the same orbitals and density of the ground state of the diatomic molecule ZnO. The effects of methodical (exchange potential) and numerical (DV method, density fit) differences between the HF and HFS methods on the various matrix elements have been analyzed. As expected the methodical effect prevails and is responsible for the higher (less negative) values of the matrix elements of the HFS Fock operator compared to those of the HF Fock operator. Numerical effects are observable also and are caused by the difference in integration procedures (DV method), not by the density fit.     

Exploitation of symmetry in periodic Self-Consistent-Field ab initio calculations: application to large three-dimensional compounds

Symmetry can dramatically reduce the computational cost (running time and memory allocation) of Self-Consistent-Field ab initio calculations for crystalline systems. Crucial for running time is use of symmetry in the evaluation of one- and two-electron integrals, diagonalization of the Fock matrix at selected points in reciprocal space, reconstruction of the density matrix. As regards memory allocation, full square matrices (overlap, Fock and density) in the Atomic Orbital (AO) basis are avoided and a direct transformation from the packed AO to the SACO (Symmetry Adapted Crystalline Orbital) basis is performed, so that the largest matrix to be handled has the size of the largest sub-block in the latter basis. We here illustrate the effectiveness of this scheme, following recent advancements in the CRYSTAL code, concerning memory allocation and direct basis set transformation. Quantitative examples are given for large unit cell systems, such as zeolites (all-silica faujasite and silicalite MFI) and garnets (pyrope). It is shown that the full SCF of 3D systems containing up to 576 atoms and 11136 Atomic Orbitals in the cell can be run with a hybrid functional on a single core PC with 500 MB RAM in about 8 h.     

Molecular symmetry. IV. The coupled perturbed Hartree–Fock method

Symmetry methods employed in the ab initio polyatomic program HONDO are extended to the coupled perturbed Hartree–Fock (CPHF) formalism, a key step in the analytical computation of energy first derivatives for configuration interaction (CI) wavefunctions, and energy second derivatives for Hartree–Fock (HF) wavefunctions. One possible computational strategy is to construct Fock-like matrices for each nuclear coordinate in which the one- and two-electron integrals of the usual Fock matrix are replaced by the integral first derivatives. “Skeleton” matrices are constructed from the unique blocks of electron-repulsion integral derivatives. The correct matrices are generated by applying a symmetrization operator. The analysis is valid for many wavefunctions, including closed- or open-shell spin-restricted and spin-unrestricted HF wavefunctions. To illustrate the method, we compare the computer time required for setting up the coupled perturbed HF equations for eclipsed ethane using D_3h symmetry point group and various subgroups of D_3h. Computational times are roughly inversely proportional to the order of the point group.     

Towards an order-N DFT method

One of the most important steps in a Kohn-Sham (KS) type density functional theory calculation is the construction of the matrix of the KS operator (the “Fock” matrix). It is desirable to develop an algorithm for this step that scales linearly with system size. We discuss attempts to achieve linear scaling for the calculation of the matrix elements of the exchange-correlation and Coulomb potentials within a particular implementation (the Amsterdam density functional, ADF, code) of the KS method. In the ADF scheme the matrix elements are completely determined by 3D numerical integration, the value of the potentials in each grid point being determined with the help of an auxiliary function representation of the electronic density. Nearly linear scaling for building the total Fock matrix is demonstrated for systems of intermediate size (in the order of 1000 atoms). For larger systems further development is desirable for the treatment of the Coulomb potential. Received: 30 March 1998 / Accepted: 6 July 1998 / Published online: 15 September 1998     

Self-consistent molecular hartree—fock—slater calculations. IV. On electron densities,spectroscopic constants and proton affinities of some small molecules

The use of the Xα exchange approximation in calculations on small molecules is studied. Electron densities are very similar to Hartree—Fock densities, as judged from density difference maps. The statistical total energy, E_Xα, is used in order to calculate R_eB_e, ω₃ and D_e of a series of diatomic molecules. The agreement with experiment is again similar to that in Hartree—Fock calculations. Proton affinities can also be calculated very well. The Hartree—Fock—Slater and Hartree—Fock models show on the whole very analogous behaviour. These results are obtained by using accurate, unapproximated, potentials and densities.     

Simple derivation of conditions for instability in the Hartree–Fock and projected Hartree–Fock schemes

The conditions for instability of solutions of Hartree–Fock and projected Hartree–Fock equations are derived in a form involving finite real symmetric matrices. These conditions are also expressed in terms of the Fock–Dirac density matrix, both at the spin–orbital and at the orbital level. The particular variations which give rise to the so-called singlet and triplet instabilities are described.     

SCF theory of intermolecular interactions without basis set superposition error

Special SCF LCAO MO type equations are derived, permitting “supermolecule” calculations for intermolecular interactions, excluding basis set superposition error (BSSE) from the beginning on the basis of the “chemical Hamiltonian approach”. (No additional “monomer” calculations are necessary to correct for BSSE.) The formalism excluding the BSSE results in a non-Hermitean Fock matrix; an algorithm is proposed to obtain the required molecular orbitals, in which no integral transformation is needed.     

Parallel computation of the Moller–Plesset second-order contribution to the electronic correlation energy

We report herein, the implementation of a second-order Moller–Plesset perturbation theory (MP2) program on the IBM LCAP parallel supercomputers. The LCAP systems comprise IBM 308X hosts and 10 FPS-X64 attached processing units (APs). The APs are interconnected by a 512 Mbyte shared memory which allows rapid interprocessor communication. All the computationally demanding steps of the MP2 procedure execute efficiently in parallel. Parallel computation of two-electron integrals is accomplished by distributing the loop over shell blocks among the APs. Parallel Fock matrix formation is achieved by having each AP evaluate the contribution of its own integral sublist to the total Fock matrix. The contributions are added together on the host, and the sum diagonalized either on the host or on a single AP. The parallel implementations of the integral transformation and the MP2 calculation are less straightforward. In each case, the use of the shared memory is essential for an efficient implementation. Details of the implementations and performance data are given.     

Restricted hartree–fock method instability

The explicit forms of the spin density matrix variations which are responsible for the external instability in the restricted Hartree–Fock (RHF ) method are found. The RHF open shell instability matrix which guarantees the distorted wave function to belong to the same space as the initial one is derived for arbitrary spin.     

On interelectronic magnetic and retardation effects within relativistic hartree-fock theory

Two forms and uses of interelectronic beyond-Coulomb interactions are discussed within the molecular analytical relativistic Hartree—Fock formalism. A completed version of the Fock matrix is derived. The matrix elements of the magnetic and the retardation interaction are evaluated by the use of gaussian basis sets.     

Hartree—Fock—Slater functions as basis for one-electron perturbation calculations for molecules: I. Calculation of ground state electric dipole polarizabilities

The analysis of the decoupling of Hartree—Fock—Slater SCF perturbation equations for an external field is undertaken. The points of departure from the corresponding Hartree—Fock perturbation equations are stressed. Both formal and numerical results suggest that the fully uncoupled Hartree—Fock—Slater expression is a less drastic approximation than the same Hartree—Fock one. The uncoupled expression for the ground state electric dipole polarizability is calculated for CO, N₂, ethylene, acethylene and trans-butadiene in the dipole length—dipole length, dipole velocity—dipole length and dipole velocity—dipole velocity alternative formulations with an ab initio Hartree—Fock—Slater SCF basis set. The results compare well with other non-empirical results and the dipole velocity-dipole length results are in remarkably good agreement with experiments.     

Extrapolation of real-space quantum chemical calculations from finite-size super cells to the ideal infinite system. II. Application to one-dimensional polymers

Calculations of finite atomic clusters based on the Hartree–Fock self-consistent field theory are modified to model more closely the ideal behavior of the infinite system. The density matrix of the standard finite supercell calculation is extrapolated to an infinite supercell so that it contains information from a continuum of k points in the first Brillouin zone. This modification is incorporated into the self-consistency loop of the MOPAC quantum chemistry program and leads to improved results compared to a standard finite supercell calculation. Heats of formation, bond lenghts, and electronic properties converge more quickly to the correct ground-state values. For polyacetylene, we obtain a reduced bond-length alternation of Δr = 0.084 Å, which is in agreement with more sophisticated calculations containing electron correlation effects. © 1994 John Wiley & Sons, Inc.     

The two-dimensional band structure of (polyphthalocyaninato)Ni(II)

The two-dimensional (2D) band structure of (polyphthalocyaninato)Ni(II), Ni(ppc), has been analyzed by a self-consistent field (SCF ) Hartree–Fock (HF ) crystal orbital (CO ) formalism based on an INDO (intermediate neglect of differential overlap) type Hamiltonian. The calculated HF band gap of Ni(ppc) amounts to 0.24 eV. The highest filled band is a ringlike a_1u combination (D_4h symmetry label) localized at the carbon sites of the organic fragment. Remarkable hybridization in the valence band leads to the considerable band width Δ?_v of 2.92 eV. This value is close to the Δ?_v numbers which are conventionally encountered in one-dimensional metallomacrocycles. The effective width of the states in Ni(ppc) is 13.8 eV. In graphite a net π interval of 13.0 eV is predicted by the present CO formalism; i.e., the energetic distribution of the π electrons is roughly comparable in both 2D solids. The Ni 3d states in Ni(ppc) are far below the Fermi level which is calculated at ?4.9 eV; they are predicted between ?12.2 and ?16.4 eV in the mean-field approximation. Quasi-particle corrections lead to a significant shift of these strongly metal-centered states. Important electronic structure properties of Ni(ppc) are compared with those of 1D metallomacrocycles with similar molecular stoichiometry. The total density of states distribution of Ni(ppc) has been fragmented into projected (ligand π and σ, Ni 3d) contributions in order to allow for a transparent interpretation of the 2D band structure.