Spectral-product methods for electronic structure calculations

Progress is reported in development, implementation, and application of a spectral method for ab initio studies of the electronic structure of matter. In this approach, antisymmetry restrictions are enforced subsequent to construction of the many-electron Hamiltonian matrix in a complete orthonormal spectral-product basis. Transformation to a permutation-symmetry representation obtained from the eigenstates of the aggregate electron antisymmetrizer is seen to enforce the requirements of the Pauli principle ex post facto, and to eliminate the unphysical (non-Pauli) states spanned by the product representation. Results identical with conventional use of prior antisymmetrization of configurational state functions are obtained in applications to many-electron atoms. The development provides certain advantages over conventional methods for polyatomic molecules, and, in particular, facilitates incorporation of fragment information in the form of Hermitian matrix representatives of atomic and diatomic operators which include the non-local effects of overall electron antisymmetry. An exact atomic-pair expression is obtained in this way for polyatomic Hamiltonian matrices which avoids the ambiguities of previously described semi-empirical fragment-based methods for electronic structure calculations. Illustrative applications to the well-known low-lying doublet states of the H₃ molecule in a minimal-basis-set demonstrate that the eigensurfaces of the antisymmetrizer can anticipate the structures of the more familiar energy surfaces, including the seams of intersection common in high-symmetry molecular geometries. The calculated H₃ energy surfaces are found to be in good agreement with corresponding valence-bond results which include all three-center terms, and are in general accord with accurate values obtained employing conventional high-level computational-chemistry procedures. By avoiding the repeated evaluations of the many-centered one- and two-electron integrals required in construction of polyatomic Hamiltonian matrices in the antisymmetric basis states commonly employed in conventional calculations, and by performing the required atomic and atomic-pair calculations once and for all, the spectral-product approach may provide an alternative potentially efficient ab initio formalism suitable for computational studies of adiabatic potential energy surfaces more generally. Contribution to the Mark S. Gordon 65th Birthday Festschrift Issue.     

The trust-region self-consistent field method in Kohn-Sham density-functional theory

The trust-region self-consistent field (TRSCF) method is extended to the optimization of the Kohn-Sham energy. In the TRSCF method, both the Roothaan-Hall step and the density-subspace minimization step are replaced by trust-region optimizations of local approximations to the Kohn-Sham energy, leading to a controlled, monotonic convergence towards the optimized energy. Previously the TRSCF method has been developed for optimization of the Hartree-Fock energy, which is a simple quadratic function in the density matrix. However, since the Kohn-Sham energy is a nonquadratic function of the density matrix, the local energy functions must be generalized for use with the Kohn-Sham model. Such a generalization, which contains the Hartree-Fock model as a special case, is presented here. For comparison, a rederivation of the popular direct inversion in the iterative subspace (DIIS) algorithm is performed, demonstrating that the DIIS method may be viewed as a quasi-Newton method, explaining its fast local convergence. In the global region the convergence behavior of DIIS is less predictable. The related energy DIIS technique is also discussed and shown to be inappropriate for the optimization of the Kohn-Sham energy.     

Atomic spectral methods for molecular electronic structure calculations

Theoretical methods are reported for ab initio calculations of the adiabatic (Born-Oppenheimer) electronic wave functions and potential energy surfaces of molecules and other atomic aggregates. An outer product of complete sets of atomic eigenstates familiar from perturbation-theoretical treatments of long-range interactions is employed as a representational basis without prior enforcement of aggregate wave function antisymmetry. The nature and attributes of this atomic spectral-product basis are indicated, completeness proofs for representation of antisymmetric states provided, convergence of Schrodinger eigenstates in the basis established, and strategies for computational implemention of the theory described. A diabaticlike Hamiltonian matrix representative is obtained, which is additive in atomic-energy and pairwise-atomic interaction-energy matrices, providing a basis for molecular calculations in terms of the (Coulombic) interactions of the atomic constituents. The spectral-product basis is shown to contain the totally antisymmetric irreducible representation of the symmetric group of aggregate electron coordinate permutations once and only once, but to also span other (non-Pauli) symmetric group representations known to contain unphysical discrete states and associated continua in which the physically significant Schrodinger eigenstates are generally embedded. These unphysical representations are avoided by isolating the physical block of the Hamiltonian matrix with a unitary transformation obtained from the metric matrix of the explicitly antisymmetrized spectral-product basis. A formal proof of convergence is given in the limit of spectral closure to wave functions and energy surfaces obtained employing conventional prior antisymmetrization, but determined without repeated calculations of Hamiltonian matrix elements as integrals over explicitly antisymmetric aggregate basis states. Computational implementations of the theory employ efficient recursive methods which avoid explicit construction the metric matrix and do not require storage of the full Hamiltonian matrix to isolate the antisymmetric subspace of the spectral-product representation. Calculations of the lowest-lying singlet and triplet electronic states of the covalent electron pair bond (H(2)) illustrate the various theorems devised and demonstrate the degree of convergence achieved to values obtained employing conventional prior antisymmetrization. Concluding remarks place the atomic spectral-product development in the context of currently employed approaches for ab initio construction of adiabatic electronic eigenfunctions and potential energy surfaces, provide comparisons with earlier related approaches, and indicate prospects for more general applications of the method.     

Basis set selections for relativistic self-consistent field calculations

Practical methods of generating reliable and economic basis sets for relativistic self-consistent fields (RSCF) calculations are developed. Large component basis sets are generated from constrained optimizations of exponents in the nonrelativistic atomic calculations for light atoms. For heavy atoms, large component basis sets for inner core orbitals are generated by fitting numerical atomic spinors of Dirac-Hartree-Fock calculations with appropriate number of Slater-type functions. Small component basis sets are obtained by using the kinetic balance condition and other computational criteria. With judicious selections of the basis sets, virtual orbitals in RSCF calculations become very similar to those in nonrelativistic calculations, implying that relativistic virtual orbitals can be used in electron correlation calculations in the same manner as the conventional nonrelativistic virtual orbitals. It is also evident that the Koopmans' theorem is also valid in RSCF results.     

Solution of self-consistent field electronic structure equations by a pseudospectral method

A new method for solving self-consistent electronic structure equations which offers orders of magnitude reductions in computation time and storage space with no obvious loss of accuracy is described and applied to the restricted Hartree-Fock equations for the neon atom. Application to polyatomic molecules is straightforward both conceptually and numerically.     

Hard-surface effects in polymer self-consistent field calculations

We have investigated several effects due to the confinement of polymer melts by impenetrable (hard) surfaces in the self-consistent field calculations. To adequately represent such confinement, the total (normalized) polymer segmental density (volume fraction) is usually constrained to an imposed profile that continuously decreases from 1 in the interior of confined melts to 0 at the surfaces over a short distance. The choice of this profile strongly influences the numerical performance of the self-consistent field calculations. In addition, for diblock copolymers A-B the hard-surface confinement has both energetic and entropic effects: On one hand, the decrease of polymer density from 1 reduces A-B repulsion and favors morphologies with more A-B interfaces near the surfaces. On the other hand, the enrichment of chain ends and depletion of middle segments near the surfaces favor parallel morphologies where chains orient mainly perpendicular to the surfaces. These two effects are comparable in magnitude, and for asymmetric diblock copolymers result in an entropic preference of a neutral surface for the shorter block as proposed previously [Q. Wang et al., Macromolecules 34, 3458 (2001)]. The hard-surface effects are weak in practice and thus manifested only when the surfaces are nearly neutral.     

Approaches to large-scale parallel self-consistent field calculations

The availability of massively parallel computers with high computation rates but limited memory and input/output bandwidth prompts the reevaluation of appropriate solution schemes for the self-consistent field (SCF) equations. Several algorithms are considered which exhibit between linear and quadratic convergence using various approximations to the orbital Hessian. A prototype is developed to understand the computational expense of each approach. The optimal choice is found to be a conjugate–gradient method preconditioned with a level-shifted approximation to the orbital Hessian. This is a compromise between efficiency, stability, and low memory usage. Sample benchmarks on two parallel supercomputers are also reported. © 1995 John Wiley & Sons, Inc.     

Wavelets for electronic structure calculations

Molecular electronic structure calculations have a multi‐scale character through the presence of a set of singularities corresponding to atomic nuclei, and thus there exists a potential to improve the efficiency of these calculations using fast wavelet transform techniques. We report on the development of a one dimensional prototype benchmark problem of sufficient complexity to capture the features of 3‐D problems that are being solved today in quantum electronics calculations. Theoretical estimates of decay across scales and spatial distribution of wavelet coefficients for the solutions of the 1‐D and 3‐D problems are derived and verified experimentally. Equivalence in a multi‐resolution context of the solutions of the 1‐D prototype and the 3‐D problem is established. This revised version was published online in July 2006 with corrections to the Cover Date.     

Monte Carlo free energy calculations using electronic structure methods

The molecular mechanics-based importance sampling function (MMBIF) algorithm [R. Iftimie, D. Salahub, D. Wei, and J. Schofield, J. Chem. Phys. 113, 4852 (2000)] is extended to incorporate semiempirical electronic structure methods in the secondary Markov chain, creating a fully quantum mechanical Monte Carlo sampling method for simulations of reactive chemical systems which, unlike the MMBIF algorithm, does not require the generation of a system-specific force field. The algorithm is applied to calculating the potential of mean force for the isomerization reaction of HCN using thermodynamic integration. Constraints are implemented in the sampling using a modification of the SHAKE algorithm, including that of a fixed, arbitrary reaction coordinate. Simulation results show that sampling efficiency with the semiempirical secondary potential is often comparable in quality to force fields constructed using the methods suggested in the original MMBIF work. The semiempirical based importance sampling method presented here is a useful alternative to MMBIF sampling as it can be applied to systems for which no suitable MM force field can be constructed.     

NDDO fragment self-consistent field approximation for large electronic systems

A semi-empirical NDDO method, generalized from a similar scheme at the CNDO/2 level developed previously, is presented to treat very large molecules. The extended molecular system is divided into a relatively small subsystem where substantial chemical changes take place and an environment remaining more-or-less unperturbed during the process. Expanding the wave function on an atomic hybrid basis an SCF procedure is performed for the subsystem in the field of the iteratively determined electronic distribution of the environment. A computer program has been written for the IBM RISC System/6000 530 computer and several test calculations were done for a variety of large classical molecules, like substituted aliphatic hydrocarbons, water oligomers, and a heptapeptide. Protonation energies, proton transfer potential curves, rotational barriers, atomic net charges, and HOMO and LUMO energies, as computed by the exact version of the NDDO method, are fairly well reproduced by our approximation if the subsystem is appropriately defined. © 1992 by John Wiley & Sons, Inc.     

Convergence optimization of restricted open-shell self-consistent field calculations

Different self-consistent field (SCF) iteration schemes for open-shell systems are discussed. After a brief summary of the well-known level shifting and damping procedure, we describe the quadratically convergent SCF (QCSCF) approach based on the gradient and the Hessian matrix in a space of orbital rotation parameters. An analytical expression for the latter is derived for the general many-shell case. Starting from the expression for the energy change obtained by the QCSCF method, we then present a simplified direct procedure avoiding matrix diagonalization but also the difficulties of the QCSCF method in handling the Hessian matrix. Numerical calculations on some open-shell systems involving transition-metal complexes show that this method leads to rapid and reliable convergence of the iteration process in cases where the usual SCF procedure of iterative diagonalization tends to diverge. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 62: 617–637, 1997     

Molecular-orbital self-consistent field calculations of quinoline and its derivatives

The electronic spectra of the thiol, anionic, protonated, and mesoionic forms of 8-mercaptoquinoline and the mesoionic form of 5-mercaptoquinoline and the 5-mercapto-N-methylquinolinium ion were calculated within the CNDO (complete neglect of differential overlap) approximation. The inclusion of the d orbitals of the sulfur atom in the AO basis has virtually no effect on the energies of the electron transitions of these compounds.See [4] for communication I.Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 11, pp. 1530–1534, November, 1978.     

Pre-processing two-electron integrals for efficient utilization in many-electron self-consistent field calculations

A procedure for pre-processing a random list of two-electron integrals is described. Construction of Fock-matrices with the resulting PK-file is shown to be three times faster than use of a file containing the random two-electron integrals, indices, and a branching code.     

Linear-scaling implementation of molecular electronic self-consistent field theory

A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field (SCF) theories is presented and illustrated with applications to molecules consisting of more than 1000 atoms. The diagonalization bottleneck of traditional SCF methods is avoided by carrying out a minimization of the Roothaan-Hall (RH) energy function and solving the Newton equations using the preconditioned conjugate-gradient (PCG) method. For rapid PCG convergence, the Lowdin orthogonal atomic orbital basis is used. The resulting linear-scaling trust-region Roothaan-Hall (LS-TRRH) method works by the introduction of a level-shift parameter in the RH Newton equations. A great advantage of the LS-TRRH method is that the optimal level shift can be determined at no extra cost, ensuring fast and robust convergence of both the SCF iterations and the level-shifted Newton equations. For density averaging, the authors use the trust-region density-subspace minimization (TRDSM) method, which, unlike the traditional direct inversion in the iterative subspace (DIIS) scheme, is firmly based on the principle of energy minimization. When combined with a linear-scaling evaluation of the Fock/Kohn-Sham matrix (including a boxed fitting of the electron density), LS-TRRH and TRDSM methods constitute the linear-scaling trust-region SCF (LS-TRSCF) method. The LS-TRSCF method compares favorably with the traditional SCF/DIIS scheme, converging smoothly and reliably in cases where the latter method fails. In one case where the LS-TRSCF method converges smoothly to a minimum, the SCF/DIIS method converges to a saddle point.     

Three-dimensional numerical integration for electronic structure calculations

Two three-dimensional numerical schemes are presented for molecular integrands such as matrix alements of one-electron operators occuring in the Fock operator and expectation values of one-electron operators describing molecular properties. The schemes are based on a judicious partitioning of space so that product-Gauss integration rules can be used in each region. Convergence with the number of integration points is such that very high accuracy (8–10 digits) may be obtained with obtained with a modest number of points. The use of point group symmetry to reduce the required number of points is discussed. Examples are given for overlap, nuclear potential, and electric field gradient integrals.     

The trust-region self-consistent field method: towards a black-box optimization in Hartree-Fock and Kohn-Sham theories

The trust-region self-consistent field (TRSCF) method is presented for optimizing the total energy E(SCF) of Hartree-Fock theory and Kohn-Sham density-functional theory. In the TRSCF method, both the Fock/Kohn-Sham matrix diagonalization step to obtain a new density matrix and the step to determine the optimal density matrix in the subspace of the density matrices of the preceding diagonalization steps have been improved. The improvements follow from the recognition that local models to E(SCF) may be introduced by carrying out a Taylor expansion of the energy about the current density matrix. At the point of expansion, the local models have the same gradient as E(SCF) but only an approximate Hessian. The local models are therefore valid only in a restricted region-the trust region-and steps can only be taken with confidence within this region. By restricting the steps of the TRSCF model to be inside the trust region, a monotonic and significant reduction of the total energy is ensured in each iteration of the TRSCF method. Examples are given where the TRSCF method converges monotonically and smoothly, but where the standard DIIS method diverges.     

Accurate ab initio density fitting for multiconfigurational self-consistent field methods

Using Cholesky decomposition and density fitting to approximate the electron repulsion integrals, an implementation of the complete active space self-consistent field (CASSCF) method suitable for large-scale applications is presented. Sample calculations on benzene, diaquo-tetra-mu-acetato-dicopper(II), and diuraniumendofullerene demonstrate that the Cholesky and density fitting approximations allow larger basis sets and larger systems to be treated at the CASSCF level of theory with controllable accuracy. While strict error control is an inherent property of the Cholesky approximation, errors arising from the density fitting approach are managed by using a recently proposed class of auxiliary basis sets constructed from Cholesky decomposition of the atomic electron repulsion integrals.     

Nonvariational time-dependent multiconfiguration self-consistent field equations for electronic dynamics in laser-driven molecules

A time-dependent multiconfiguration self-consistent field (TDMCSCF) scheme is developed to describe the time-resolved electron dynamics of a laser-driven many-electron atomic or molecular system, starting directly from the time-dependent Schrodinger equation for the system. This nonvariational formulation aims at the full exploitations of concepts, tools, and facilities of existing, well-developed quantum chemical MCSCF codes. The theory uses, in particular, a unitary representation of time-dependent configuration mixings and orbital transformations. Within a short-time, or adiabatic approximation, the TDMCSCF scheme amounts to a second-order split-operator algorithm involving generically the two noncommuting one-electron and two-electron parts of the time-dependent electronic Hamiltonian. We implement the scheme to calculate the laser-induced dynamics of the two-electron H2 molecule described within a minimal basis, and show how electron correlation is affected by the interaction of the molecule with a strong laser field.     

Fast vibrational self-consistent field calculations through a reduced mode-mode coupling scheme

We present a new methodology to perform fast correlation-corrected vibrational self-consistent field (CC-VSCF) calculations using ab initio potential energy points calculated on the fly. Our method is based on the replacement of all-electron basis sets with a pseudo-potential basis for heavy atoms, and on an efficient reduction of the number of pair-coupling elements used in the CC-VSCF procedure. The method is applied to several test systems: H2O, NH3, and CH4, where it proves to be efficient, providing a speedup factor of 2 compared to a standard CC-VSCF calculation. We also apply our technique to the simulation of the vibrational spectrum of ethane and show that very accurate results can be obtained with a substantial speedup for this system.