Orbital optimization in nonorthogonal multiconfigurational self-consistent field applied to the study of conical intersections and avoided crossings

Selected Nonorthogonal Configuration Interaction with Compressed Single and Double Excitations

Addressing both dynamic and static correlations accurately is a primary goal in electronic structure theory. Nonorthogonal configuration interaction (NOCI) is a versatile tool for treating static correlation, offering chemical insights by combining diverse reference states. Nevertheless, achieving quantitative accuracy requires the inclusion of the missing dynamic correlation. This work introduces a framework for compressing orthogonal single and double excitations into a NOCI of a much smaller dimension. This compression is repeated with each Slater determinant in a reference NOCI, resulting in another NOCI that includes all of its single and double excitations (NOCISD), effectively recovering the missing dynamic correlations from the reference. This compressed NOCISD is further refined through a selection process using metric and energy tests (SNOCISD). We validate the effectiveness of SNOCISD through its application to the dissociation of the nitrogen molecule and the hole-doped two-dimensional Hubbard model at various interaction strengths.

Correlated pair ansatz with a binary tree structure

We develop an efficient algorithm to implement the recently introduced binary tree state (BTS) ansatz on a classical computer. BTS allows a simple approximation to permanents arising from the computationally intractable antisymmetric product of interacting geminals and respects size-consistency. We show how to compute BTS overlap and reduced density matrices efficiently. We also explore two routes for developing correlated BTS approaches: Jastrow coupled cluster on BTS and linear combinations of BT states. The resulting methods show great promise in benchmark applications to the reduced Bardeen-Cooper-Schrieffer Hamiltonian and the one-dimensional XXZ Heisenberg Hamiltonian.

Time propagation of electronic wavefunctions using nonorthogonal determinant expansions

The use of truncated configuration interaction in real-time time-dependent simulations of electron dynamics provides a balance of computational cost and accuracy, while avoiding some of the failures associated with real-time time-dependent density functional theory. However, low-order truncated configuration interaction also has limitations, such as overestimation of polarizability in configuration interaction singles, even when perturbative doubles are included. Increasing the size of the determinant expansion may not be computationally feasible, and so, in this work, we investigate the use of nonorthogonality in the determinant expansion to establish the extent to which higher-order substitutions can be recovered, providing an improved description of electron dynamics. Model systems are investigated to quantify the extent to which different methods accurately reproduce the (hyper)polarizability, including the high-harmonic generation spectrum of H2, water, and butadiene.

Nonorthogonal Active Space Decomposition of Wave Functions with Multiple Correlation Mechanisms

The nonorthogonal active space decomposition (NO-ASD) methodology is proposed for describing systems containing multiple correlation mechanisms. NO-ASD partitions the wave function by a correlation mechanism, such that the interactions between different correlation mechanisms are treated with an effective Hamiltonian approach, while interactions between correlated orbitals in the same correlation mechanism are treated explicitly. As a result, the determinant expansion scales polynomially with the number of correlation mechanisms rather than exponentially, which significantly reduces the factorial scaling associated with the size of the correlated orbital space. Despite the nonorthogonal framework of NO-ASD, the approach can take advantage of computational efficient matrix element evaluation when performing nonorthogonal coupling of orthogonal determinant expansions. In this work, we introduce and examine the NO-ASD approach in comparison to complete active space methods to establish how the NO-ASD approach reduces the problem dimensionality and the extent to which it affects the amount of correlation energy recovered. Calculations are performed on ozone, nickel-acetylene, and isomers of μ-oxo dicopper ammonia.

Role of Exact Exchange in Difference Projected Double-Hybrid Density Functional Theory for Treatment of Local, Charge Transfer, and Rydberg Excitations

Difference approaches to the study of excited states have undergone a renaissance in recent years, with the development of a plethora of algorithms for locating self-consistent field approximations to excited states. Density functional theory is likely to offer the best balance of cost and accuracy for difference approaches, and yet there has been little investigation of how the parametrization of density functional approximations affects performance. In this work, we aim to explore the role of the global Hartree-Fock exchange parameter in tuning accuracy of different excitation types within the framework of the recently introduced difference projected double-hybrid density functional theory approach and contrast the performance with conventional time-dependent double-hybrid density functional theory. Difference projected double-hybrid density functional theory was demonstrated to give vertical excitation energies with average error and standard deviation with respect to multireference perturbation theory comparable to more expensive linear-response coupled cluster approaches ( J. Chem. Phys.2020, 153, 074103). However, despite benchmarking of local excitations, there has been no investigation of the methods performance for charge transfer or Rydberg excitations. In this work we report a new benchmark of charge transfer, Rydberg, and local excited state vertical excitation energies and examine how the exact Hartree-Fock exchange affects the benchmark performance to provide a deeper understanding of how projection and nonlocal correlation balance differing sources of error in the ground and excited states.

Multistate Density Functional Theory of Excited States

We report a rigorous formulation of density functional theory for excited states, providing a theoretical foundation for a multistate density functional theory. We prove the existence of a Hamiltonian matrix functional H[D] of the multistate matrix density D(r) in the subspace spanned by the lowest N eigenstates. Here, D(r) is an N-dimensional matrix of state densities and transition densities. Then, a variational principle of the multistate subspace energy is established, whose minimization yields both the energies and densities of the individual N eigenstates. Furthermore, we prove that the N-dimensional matrix density D(r) can be sufficiently represented by N² nonorthogonal Slater determinants, based on which an interacting active space is introduced for practical calculations. This work establishes that the ground and excited states can be treated on an equal footing in density functional theory.

GronOR: Scalable and Accelerated Nonorthogonal Configuration Interaction for Molecular Fragment Wave Functions

GronOR is a program package for nonorthogonal configuration interaction calculations. Electronic wave functions are constructed in terms of antisymmetrized products of multiconfiguration molecular fragment wave functions. The computational complexity of the nonorthogonal methodologies implemented in GronOR applied to large molecular assemblies requires a design that takes full advantage of massively parallel supercomputer architectures and accelerator technologies. This work describes the implementation strategy and resulting performance characteristics. In addition to parallelization and acceleration, the software development strategy includes aspects of fault resiliency and heterogeneous computing. The program was designed for large-scale supercomputers but also runs effectively on small clusters and workstations for small molecular systems. GronOR is available as open source to the scientific community.