iOI: An Iterative Orbital Interaction Approach for Solving the Self-Consistent Field Problem

Method and Implementation of Projected Hybrid Orbitals for Treating Multiple Covalent Bonds in Combined QM/MM Calculations

The projected hybrid orbital (PHO) method for treatment of multiple boundary atoms introduces a novel solution for handling the covalent connection between quantum mechanical (QM) and molecular mechanical (MM) regions in QM/MM calculations. By projecting the QM basis, typically adequately large for computational accuracy, onto a secondary minimal basis set on the boundary atom, it preserves electronic interactions between the two regions without system-specific parameters. Applicable in both ab initio wave function theory and density functional theory, PHO has been shown to maintain structural and electronic integrity across various systems. It offers key advantages over traditional QM/MM methods, such as avoiding modification of MM charge distribution and energy corrections, while being flexible for biochemical systems and potentially broader QM/QM embedding frameworks.

Spin-Adapted Externally Contracted Multireference Configuration Interaction Method Based on Selected Reference Configurations

As one kind of approximation of the full configuration interaction solution, the selected configuration interaction (sCI) methods have been shown to be valuable for large active spaces. However, the inclusion of dynamic correlation beyond large active spaces is necessary for more quantitative results. Since the sCI wave function can provide a compact reference for multireference methods, previously, we proposed an externally contracted multireference configuration interaction method using the sCI reference reconstructed from the density matrix renormalization group wave function [J. Chem. Theory Comput. 2018, 14, 4747-4755]. The DMRG2sCI-EC-MRCI method is promising for dealing with more than 30 active orbitals and large basis sets. However, it suffers from two drawbacks: spin contamination and low efficiency when using Slater determinant bases. To solve these problems, in this work, we adopt configuration state function bases and introduce a new algorithm based on the hybrid of tree structure for convenient configuration space management and the graphical unitary group approach for efficient matrix element calculation. The test calculation of naphthalene shows that the spin-adapted version could achieve a speed-up of 6.0 compared with the previous version based on the Slater determinant. Examples of dinuclear copper(II) compound as well as Ln(III) and An(III) complexes show that the sCI-EC-MRCI can give quantitatively accurate results by including dynamic correlation over sCI for systems with large active spaces and basis sets.

Extending multi-layer energy-based fragment method for excited-state calculations of large covalently bonded fragment systems

Recently, we developed a low-scaling Multi-Layer Energy-Based Fragment (MLEBF) method for accurate excited-state calculations and nonadiabatic dynamics simulations of nonbonded fragment systems. In this work, we extend the MLEBF method to treat covalently bonded fragment ones. The main idea is cutting a target system into many fragments according to chemical properties. Fragments with dangling bonds are first saturated by chemical groups; then, saturated fragments, together with the original fragments without dangling bonds, are grouped into different layers. The accurate total energy expression is formulated with the many-body energy expansion theory, in combination with the inclusion-exclusion principle that is used to delete the contribution of chemical groups introduced to saturate dangling bonds. Specifically, in a two-layer MLEBF model, the photochemically active and inert layers are calculated with high-level and efficient electronic structure methods, respectively. Intralayer and interlayer energies can be truncated at the two- or three-body interaction level. Subsequently, through several systems, including neutral and charged covalently bonded fragment systems, we demonstrate that MLEBF can provide accurate ground- and excited-state energies and gradients. Finally, we realize the structure, conical intersection, and path optimizations by combining our MLEBF program with commercial and free packages, e.g., ASE and SciPy. These developments make MLEBF a practical and reliable tool for studying complex photochemical and photophysical processes of large nonbonded and bonded fragment systems.

Tuning Catalyst-Free Photocontrolled Polymerization by Substitution: A Quantitative and Qualitative Interpretation

Catalyst-free photocontrolled reversible addition-fragmentation chain transfer (RAFT) polymerization avoids the side effects of photocatalysts but has the accompanying slow kinetics, thereby warranting more efficient photolysis and faster chain transfer. To understand the underlying mechanisms, both quantitative and qualitative interpretations are needed. Such a goal can be achieved by the iCAS (imposed automatic selection and localization of complete active spaces) approach [J. Chem. Theory Comput. 2021, 17, 4846], which maintains the same CAS and meanwhile provides localized orbitals along the whole reaction. Taking dithiobenzoate as a representative of RAFT agents, it is found here that electron-donating substitution (by methoxy) clearly outperforms both electron-standing (by methyl) and electron-withdrawing (by cyano) substitutions in facilitating photo-RAFT polymerization, by narrowing the gap between the π* and σ* orbitals, so as to facilitate the π* → σ* charge transfer dominating both the photolysis and chain transfer processes. Such findings are of general values.

iCISCF: An Iterative Configuration Interaction-Based Multiconfigurational Self-Consistent Field Theory for Large Active Spaces

An iterative configuration interaction (iCI)-based multiconfigurational self-consistent field (SCF) theory, iCISCF, is proposed to handle systems that require large active spaces. The success of iCISCF stems from three ingredients: (1) efficient selection of individual configuration state functions spanning the active space while maintaining full spin symmetry; (2) the use of Jacobi rotation for optimization of the active orbitals in conjunction with a quasi-Newton algorithm for the core/active-virtual and core-active orbital rotations; (3) a second-order perturbative treatment of the residual space left over by the selection procedure (i.e., iCISCF(2)). Several examples that go beyond the capability of CASSCF are taken as showcases to reveal the efficacy of iCISCF and iCISCF(2), facilitated by iCAS for imposed automatic selection and localization of active orbitals.

iCAS: Imposed Automatic Selection and Localization of Complete Active Spaces

It is shown that in the spirit of "from fragments to molecule" for localizing molecular orbitals [J. Chem. Theory Comput. 2011, 7, 3643], a prechosen set of occupied/virtual valence/core atomic/fragmental orbitals can be transformed to an equivalent set of localized occupied/virtual pre-localized molecular orbitals (pre-LMO), which can then be taken as probes to select the same number of maximally matching localized occupied/virtual Hartree-Fock (HF) or restricted open-shell HF (ROHF) molecular orbitals as the initial local orbitals spanning the desired complete active space (CAS). In each cycle of the self-consistent field (SCF) calculation, the CASSCF orbitals can be localized by means of the noniterative "top-down least-change" algorithm for localizing ROHF orbitals [J. Chem. Phys. 2017, 146, 104104] such that the maximum matching between the orbitals of two adjacent iterations can readily be monitored, leading finally to converged localized CASSCF orbitals that overlap most the guess orbitals. Such an approach is to be dubbed as "imposed CASSCF" (iCASSCF or simply iCAS in short) for good reasons: (1) it has been assumed that only those electronic states that have largest projections onto the active space defined by the prechosen atomic/fragmental orbitals are to be targeted. This is certainly an imposed constraint but has wide applications in organic and transition metal chemistry where valence (or core) atomic/fragmental orbitals can readily be identified. (2) The selection of both initial and optimized local active orbitals is imposed from the very beginning by the pre-LMOs (which span the same space as the prechosen atomic/fragmental orbitals). Apart from the (imposed) automation and localization, iCAS has two additional merits: (a) the guess orbitals are guaranteed to be the same for all geometries, for the pre-LMOs do not change in character with geometry and (b) the use of localized orbitals facilitates the SCF convergence, particularly for large active spaces. Both organic molecules and transition-metal complexes are taken as showcases to reveal the efficacy of iCAS.