Graphically Contracted Function Construction with the Recursive Pairwise Merge Algorithm.

A new implementation of a recursive pairwise merge algorithm to construct a GCF from a list of CSF expansion coefficients is presented. The essential new feature is the preallocation of some work arrays used within the intermediate steps of the merge procedure. This results in roughly an order of magnitude improvement in overall efficiency and also approximately eliminates a factor of , the molecular orbital dimension, from the original implementation.

Quantification of the Ionic Character of Multiconfigurational Wave Functions: The Diagnostic.

The complete active space self-consistent field (CASSCF) method is a cornerstone in modern excited-state quantum chemistry providing the starting point for most common multireference computations. However, CASSCF, when used with a minimal active space, can produce significant errors (>2 eV) even for the excitation energies of simple hydrocarbons if the states of interest possess ionic character. After illustrating this problem in some detail, we present a diagnostic for ionic character, denoted as  , that is readily computed from the transition density.

The OpenMolcas : A Community-Driven Approach to Advancing Computational Chemistry.

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Observation of site-selective chemical bond changes via ultrafast chemical shifts.

The concomitant motion of electrons and nuclei on the femtosecond time scale marks the fate of chemical and biological processes. Here we demonstrate the ability to initiate and track the ultrafast electron rearrangement and chemical bond breaking site-specifically in real time for the carbon monoxide diatomic molecule. We employ a local resonant x-ray pump at the oxygen atom and probe the chemical shifts of the carbon core-electron binding energy.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry.

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g.

Representations of Shavitt Graphs Within the Graphical Unitary Group Approach.

The Shavitt graph is a visual representation of a distinct row table (DRT) within the graphical unitary group approach. The DRT is a compact representation of the entire configuration state function expansion space within a molecular electronic structure calculation. Each node of the graph is associated with an integer triple (a ,b ,c ).

The Representation and Parametrization of Orthogonal Matrices.

Four representations and parametrizations of orthogonal matrices Q ∈ R(m×n) in terms of the minimal number of essential parameters {φ} are discussed: the exponential representation, the Householder reflector representation, the Givens rotation representation, and the rational Cayley transform representation. Both square n = m and rectangular n < m situations are considered. Two separate kinds of parametrizations are considered: one in which the individual columns of Q are distinct, the Stiefel manifold, and the other in which only span(Q) is significant, the Grassmann manifold.

The multifacet graphically contracted function method. II. A general procedure for the parameterization of orthogonal matrices and its application to arc factors.

Practical algorithms are presented for the parameterization of orthogonal matrices Q ∈ R(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 multifacet graphically contracted function method. I. Formulation and implementation.

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, spin-density matrix construction, the computation of optimization gradients for single-state and state-averaged 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 full-CI Shavitt graphs, and it is observed that both the energy and the gradient computation scale as O(N(2)n(4)) for N electrons and n orbitals.

The multiradical character of one- and two-dimensional graphene nanoribbons.

When is an acene stable? The pronounced multiradical character of graphene nanoribbons of different size and shape was investigated with high-level multireference methods. Quantitative information based on the number of effectively unpaired electrons leads to specific estimates of the chemical stability of graphene nanostructures.

Evaluation of the spin-orbit interaction within the graphically contracted function method.

The graphically contracted function (GCF) method is extended to include an effective one-electron spin-orbit (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.

Computation of determinant expansion coefficients within the graphically contracted function method.

Most electronic structure methods express the wavefunction as an expansion of N-electron 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.

Hamiltonian matrix and reduced density matrix construction with nonlinear wave functions.

An efficient procedure to compute Hamiltonian matrix elements and reduced one- and two-particle density matrices for electronic wave functions using a new graphical-based nonlinear expansion form is presented. This method is based on spin eigenfunctions using the graphical unitary group approach (GUGA), and the wave function is expanded in a basis of product functions (each of which is equivalent to some linear combination of all of the configuration state functions), allowing application to closed- and open-shell systems and to ground and excited electronic states. In general, the effort required to construct an individual Hamiltonian matrix element between two product basis functions H(MN) = M|H|N scales as theta (beta n4) for a wave function expanded in n molecular orbitals.

A general nonlinear expansion form for electronic wave functions.

A new expansion form is presented for electronic wave functions. The wave function is a linear combination of product basis functions, and each product basis function in turn is formally equivalent to a linear combination of configuration state functions that comprise an underlying linear expansion space. The expansion coefficients that define the basis functions are nonlinear functions of a smaller number of variables.

Ab initio study of AmCl(+): f-f spectroscopy and chemical binding.

A valence full configuration interaction study with a polarized double-zeta quality basis set has been carried out for the lowest 49 electronic states of AmCl(+). The calculations use a pseudopotential treatment for the core electrons and incorporate a one-electron spin-orbit interaction operator. Electrons in the valence s, p, d, and f subshells were included in the active space.

Analytic evaluation of nonadiabatic coupling terms at the MR-CI level. II. Minima on the crossing seam: formaldehyde and the photodimerization of ethylene.

The method for the analytic calculation of the nonadiabatic coupling vector at the multireference configuration-interaction (MR-CI) level and its program implementation into the COLUMBUS program system described in the preceding paper [Lischka et al., J. Chem.

Analytic evaluation of nonadiabatic coupling terms at the MR-CI level. I. Formalism.

An efficient and general method for the analytic computation of the nonandiabatic coupling vector at the multireference configuration interaction (MR-CI) level is presented. This method is based on a previously developed formalism for analytic MR-CI gradients adapted to the use for the computation of nonadiabatic coupling terms. As was the case for the analytic energy gradients, very general, separate choices of invariant orbital subspaces at the multiconfiguration self-consistent field and MR-CI levels are possible, allowing flexible selections of MR-CI wave functions.

Reducing I/O costs for the eigenvalue procedure in large-scale configuration interaction calculations.

Several aspects of the matrix diagonalization method used for CI calculations in the COLUMBUS Program System are discussed, including a linear basis-contraction algorithm and the use of a nonorthogonal expansion basis. Both of these features significantly reduce the I/O requirements during the iterations.