Global Elucidation of Self-Consistent Field Solution Space Using Basin Hopping

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 Multireference Wave Function Description of Triplet-Triplet Energy Transfer Couplings

In this study, the use of self-consistent field quasi-diabats is investigated for calculation of triplet energy transfer diabatic coupling elements. It is proposed that self-consistent field quasi-diabats are particularly useful for studying energy transfer (EnT) processes because orbital relaxation in response to changes in electron configuration is implicitly built into the model. The conceptual model that is developed allows for the simultaneous evaluation of direct and charge-transfer mechanisms to establish the importance of the different possible EnT mechanisms. The method's performance is evaluated using two model systems: the ethylene dimer and ethylene with the methaniminium cation. While states that mediate the charge-transfer mechanism were found to be higher in energy than the states involved in the direct mechanism, the coupling elements that control the kinetics were found to be significantly larger in the charge-transfer mechanism. Subsequently, we discuss the advantage of the approach in the context of practical difficulties with the use of established approaches.

Excited States, Symmetry Breaking, and Unphysical Solutions in State-Specific CASSCF Theory

State-specific electronic structure theory provides a route toward balanced excited-state wave functions by exploiting higher-energy stationary points of the electronic energy. Multiconfigurational wave function approximations can describe both closed- and open-shell excited states and avoid the issues associated with state-averaged approaches. We investigate the existence of higher-energy solutions in complete active space self-consistent field (CASSCF) theory and characterize their topological properties. We demonstrate that state-specific approximations can provide accurate higher-energy excited states in H₂ (6-31G) with more compact active spaces than would be required in a state-averaged formalism. We then elucidate the unphysical stationary points, demonstrating that they arise from redundant orbitals when the active space is too large or symmetry breaking when the active space is too small. Furthermore, we investigate the singlet-triplet crossing in CH₂ (6-31G) and the avoided crossing in LiF (6-31G), revealing the severity of root flipping and demonstrating that state-specific solutions can behave quasi-diabatically or adiabatically. These results elucidate the complexity of the CASSCF energy landscape, highlighting the advantages and challenges of practical state-specific calculations.

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.

Excited States Calculations of MoS2@ZnO and WS2@ZnO Two-Dimensional Nanocomposites for Water-Splitting Applications

Transition metal dichalcogenide (TMD) MoS2 and WS2 monolayers (MLs) deposited atop of crystalline zinc oxide (ZnO) and graphene-like ZnO (g-ZnO) substrates have been investigated by means of density functional theory (DFT) using PBE and GLLBSC exchange-correlation functionals. In this work, the electronic structure and optical properties of studied hybrid nanomaterials are described in view of the influence of ZnO substrates thickness on the MoS2@ZnO and WS2@ZnO two-dimensional (2D) nanocomposites. The thicker ZnO substrate not only triggers the decrease of the imaginary part of dielectric function relatively to more thinner g-ZnO but also results in the less accumulated charge density in the vicinity of the Mo and W atoms at the conduction band minimum. Based on the results of our calculations, we predict that MoS2 and WS2 monolayers placed at g-ZnO substrate yield essential enhancement of the photoabsorption in the visible region of solar spectra and, thus, can be used as a promising catalyst for photo-driven water splitting applications.

Variational coupled cluster for ground and excited states

In single-reference coupled-cluster (CC) methods, one has to solve a set of non-linear polynomial equations in order to determine the so-called amplitudes that are then used to compute the energy and other properties. Although it is of common practice to converge to the (lowest-energy) ground-state solution, it is also possible, thanks to tailored algorithms, to access higher-energy roots of these equations that may or may not correspond to genuine excited states. Here, we explore the structure of the energy landscape of variational CC and we compare it with its (projected) traditional version in the case where the excitation operator is restricted to paired double excitations (pCCD). By investigating two model systems (the symmetric stretching of the linear H₄ molecule and the continuous deformation of the square H₄ molecule into a rectangular arrangement) in the presence of weak and strong correlations, the performance of variational pCCD (VpCCD) and traditional pCCD is gauged against their configuration interaction (CI) equivalent, known as doubly occupied CI, for reference Slater determinants made of ground- or excited-state Hartree-Fock orbitals or state-specific orbitals optimized directly at the VpCCD level. The influence of spatial symmetry breaking is also investigated.

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

Nonorthogonal approaches to electronic structure methods have recently received renewed attention, with the hope that new forms of nonorthogonal wavefunction Ansätze may circumvent the computational bottleneck of orthogonal-based methods. The basis in which nonorthogonal configuration interaction is performed defines the compactness of the wavefunction description and hence the efficiency of the method. Within a molecular orbital approach, nonorthogonal configuration interaction is defined by a "different orbitals for different configurations" picture, with different methods being defined by their choice of determinant basis functions. However, identification of a suitable determinant basis is complicated, in practice, by (i) exponential scaling of the determinant space from which a suitable basis must be extracted, (ii) possible linear dependencies in the determinant basis, and (iii) inconsistent behavior in the determinant basis, such as disappearing or coalescing solutions, as a result of external perturbations, such as geometry change. An approach that avoids the aforementioned issues is to allow for basis determinant optimization starting from an arbitrarily constructed initial determinant set. In this work, we derive the equations required for performing such an optimization, extending previous work by accounting for changes in the orthogonality level (defined as the dimension of the orbital overlap kernel between two determinants) as a result of orbital perturbations. The performance of the resulting wavefunction for studying avoided crossings and conical intersections where strong correlation plays an important role is examined.