The multifacet graphically contracted function method. I. Formulation and implementation
Abstract
The basic formulation for the multifacet generalization of the graphically contracted function (MFGCF) electronic structure method is presented. The analysis includes the discussion of linear dependency and redundancy of the arc factor parameters, the computation of reduced density matrices, Hamiltonian matrix construction, spindensity matrix construction, the computation of optimization gradients for singlestate and stateaveraged calculations, graphical wave function analysis, and the efficient computation of configuration state function and Slater determinant expansion coefficients. Timings are given for Hamiltonian matrix element and analytic optimization gradient computations for a range of model problems for fullCI Shavitt graphs, and it is observed that both the energy and the gradient computation scale as O(N{sup 2}n{sup 4}) for N electrons and n orbitals. The important arithmetic operations are within dense matrixmatrix product computational kernels, resulting in a computationally efficient procedure. An initial implementation of the method is used to present applications to several challenging chemical systems, including N{sub 2} dissociation, cubic H{sub 8} dissociation, the symmetric dissociation of H{sub 2}O, and the insertion of Be into H{sub 2}. The results are compared to the exact fullCI values and also to those of the previous singlefacet GCF expansion form.
 Authors:
 Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439 (United States)
 Department of Chemistry and Biochemistry, Gonzaga University, 502 E. Boone Ave., Spokane, Washington 992580102 (United States)
 Publication Date:
 OSTI Identifier:
 22420011
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Chemical Physics; Journal Volume: 141; Journal Issue: 6; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; DENSITY MATRIX; DISSOCIATION; ELECTRONIC STRUCTURE; ELECTRONS; HAMILTONIANS; KERNELS; SLATER METHOD; SPIN
Citation Formats
Shepard, Ron, Brozell, Scott R., and Gidofalvi, Gergely. The multifacet graphically contracted function method. I. Formulation and implementation. United States: N. p., 2014.
Web. doi:10.1063/1.4890734.
Shepard, Ron, Brozell, Scott R., & Gidofalvi, Gergely. The multifacet graphically contracted function method. I. Formulation and implementation. United States. doi:10.1063/1.4890734.
Shepard, Ron, Brozell, Scott R., and Gidofalvi, Gergely. Thu .
"The multifacet graphically contracted function method. I. Formulation and implementation". United States.
doi:10.1063/1.4890734.
@article{osti_22420011,
title = {The multifacet graphically contracted function method. I. Formulation and implementation},
author = {Shepard, Ron and Brozell, Scott R. and Gidofalvi, Gergely},
abstractNote = {The basic formulation for the multifacet generalization of the graphically contracted function (MFGCF) electronic structure method is presented. The analysis includes the discussion of linear dependency and redundancy of the arc factor parameters, the computation of reduced density matrices, Hamiltonian matrix construction, spindensity matrix construction, the computation of optimization gradients for singlestate and stateaveraged calculations, graphical wave function analysis, and the efficient computation of configuration state function and Slater determinant expansion coefficients. Timings are given for Hamiltonian matrix element and analytic optimization gradient computations for a range of model problems for fullCI Shavitt graphs, and it is observed that both the energy and the gradient computation scale as O(N{sup 2}n{sup 4}) for N electrons and n orbitals. The important arithmetic operations are within dense matrixmatrix product computational kernels, resulting in a computationally efficient procedure. An initial implementation of the method is used to present applications to several challenging chemical systems, including N{sub 2} dissociation, cubic H{sub 8} dissociation, the symmetric dissociation of H{sub 2}O, and the insertion of Be into H{sub 2}. The results are compared to the exact fullCI values and also to those of the previous singlefacet GCF expansion form.},
doi = {10.1063/1.4890734},
journal = {Journal of Chemical Physics},
number = 6,
volume = 141,
place = {United States},
year = {Thu Aug 14 00:00:00 EDT 2014},
month = {Thu Aug 14 00:00:00 EDT 2014}
}

Practical algorithms are presented for the parameterization of orthogonal matrices Q ∈ R {sup m×n} in terms of the minimal number of essential parameters (φ). Both square n = m and rectangular n < m situations are examined. Two separate kinds of parameterizations are considered, one in which the individual columns of Q are distinct, and the other in which only Span(Q) is significant. The latter is relevant to chemical applications such as the representation of the arc factors in the multifacet graphically contracted function method and the representation of orbital coefficients in SCF and DFT methods. The parameterizations aremore »

Evaluation of the spinorbit interaction within the graphically contracted function method.
The graphically contracted function (GCF) method is extended to include an effective oneelectron spinorbit (SO) operator in the Hamiltonian matrix construction. Our initial implementation is based on a multiheaded Shavitt graph approach that allows for the efficient simultaneous computation of entire blocks of Hamiltonian matrix elements. Two algorithms are implemented. The SOGCF method expands the spinorbit wave function in the basis of GCFs and results in a Hamiltonian matrix of dimension N{sub dim} = N{sub a} ((S{sub max} + 1){sup 2} ? S{sub min}{sup 2}). N{sub a} is the number of sets of nonlinear arc factor parameters, and S{sub min}more » 
The evaluation of spindensity matrices within the graphically contracted function method.
An efficient algorithm is presented to compute spindensity matrices from wave functions expanded in a basis of graphically contracted functions (GCF). The GCFs are based on the graphical unitary group approach (GUGA), which is a 'spinfree' formulation of the electronic wave function. The spindensity matrix elements are computed from oneparticle and twoparticle chargedensity matrix elements. The recursive algorithm allows the spindensity matrix to be computed with O(N{sub 2}{sup GCF}{omega}n{sup 2}) total effort where N{sub GCF} is the dimension of the GCF basis and n is the dimension of the orbital basis. The scale factor {omega} depends on the number ofmore » 
Computation of determinant expansion coefficients within the graphically contracted function method.
Most electronic structure methods express the wavefunction as an expansion of Nelectron basis functions that are chosen to be either Slater determinants or configuration state functions. Although the expansion coefficient of a single determinant may be readily computed from configuration state function coefficients for small wavefunction expansions, traditional algorithms are impractical for systems with a large number of electrons and spatial orbitals. In this work, we describe an efficient algorithm for the evaluation of a single determinant expansion coefficient for wavefunctions expanded as a linear combination of graphically contracted functions. Each graphically contracted function has significant multiconfigurational character and dependsmore »