Distinguishing homolytic vs heterolytic bond dissociation of phenylsulfonium cations with localized active space methods.

Modeling chemical reactions with quantum chemical methods is challenging when the electronic structure varies significantly throughout the reaction and when electronic excited states are involved. Multireference methods, such as complete active space self-consistent field (CASSCF), can handle these multiconfigurational situations. However, even if the size of the needed active space is affordable, in many cases, the active space does not change consistently from reactant to product, causing discontinuities in the potential energy surface. The localized active space SCF (LASSCF) is a cheaper alternative to CASSCF for strongly correlated systems with weakly correlated fragments. The method is used for the first time to study a chemical reaction, namely the bond dissociation of a mono-, di-, and triphenylsulfonium cation. LASSCF calculations generate smooth potential energy scans more easily than the corresponding, more computationally expensive CASSCF calculations while predicting similar bond dissociation energies. Our calculations suggest a homolytic bond cleavage for di- and triphenylsulfonium and a heterolytic pathway for monophenylsulfonium.

A Synergistic Approach towards Optimization of Coupled Cluster Amplitudes by Exploiting Dynamical Hierarchy.

The coupled cluster iteration scheme for determining the cluster amplitudes involves a set of nonlinearly coupled difference equations. In the space spanned by the amplitudes, the set of equations are analyzed as a multivariate time-discrete map where the concept of time appears in an implicit manner. With the observation that the cluster amplitudes have difference in their relaxation timescales with respect to the distributions of their magnitudes, the coupled cluster iteration dynamics are considered as a synergistic motion of coexisting slow and fast relaxing modes, manifesting a dynamical hierarchical structure. With the identification of the highly damped auxiliary amplitudes, their time variation can be neglected compared to the principal amplitudes which take much longer time to reach the fixed points. We analytically establish the adiabatic approximation where each of these auxiliary amplitudes are expressed as unique parametric functions of the collective principal amplitudes, allowing us to study the optimization with the latter taken as the independent degrees of freedom. Such decoupling of the amplitudes significantly reduces the computational scaling without sacrificing the accuracy in the ground state energy as demonstrated by a number of challenging molecular applications. A road-map to treat higher order post-adiabatic effects is also discussed.

A hybrid coupled cluster-machine learning algorithm: Development of various regression models and benchmark applications.

The iterative solution of the coupled cluster equations exhibits a synergistic relationship among the various cluster amplitudes. The iteration scheme is analyzed as a multivariate discrete time propagation of nonlinearly coupled equations, which is dictated by only a few principal cluster amplitudes. These principal amplitudes usually correspond to only a few valence excitations, whereas all other cluster amplitudes are enslaved and behave as auxiliary variables [Agarawal et al., J. Chem. Phys. 154, 044110 (2021)]. We develop a coupled cluster-machine learning hybrid scheme where various supervised machine learning strategies are introduced to establish the interdependence between the principal and auxiliary amplitudes on-the-fly. While the coupled cluster equations are solved only to determine the principal amplitudes, the auxiliary amplitudes, on the other hand, are determined via regression as unique functionals of the principal amplitudes. This leads to significant reduction in the number of independent degrees of freedom during the iterative optimization, which saves significant computation time. A few different regression techniques have been developed, which have their own advantages and disadvantages. The scheme has been applied to several molecules in their equilibrium and stretched geometries, and our scheme, with all the regression models, shows a significant reduction in computation time over the canonical coupled cluster calculations without unduly sacrificing the accuracy.

An approximate coupled cluster theory via nonlinear dynamics and synergetics: The adiabatic decoupling conditions.

The coupled cluster iteration scheme is analyzed as a multivariate discrete time map using nonlinear dynamics and synergetics. The nonlinearly coupled set of equations to determine the cluster amplitudes are driven by a fraction of the entire set of cluster amplitudes. These driver amplitudes enslave all other amplitudes through a synergistic inter-relationship, where the latter class of amplitudes behave as the auxiliary variables. The driver and the auxiliary variables exhibit vastly different time scales of relaxation during the iteration process to reach the fixed points. The fast varying auxiliary amplitudes are small in magnitude, while the driver amplitudes are large, and they have a much longer time scale of relaxation. Exploiting their difference in relaxation time scale, we employ an adiabatic decoupling approximation, where each of the fast relaxing auxiliary modes is expressed as a unique function of the principal amplitudes. This results in a tremendous reduction in the independent degrees of freedom. On the other hand, only the driver amplitudes are determined accurately via exact coupled cluster equations. We will demonstrate that the iteration scheme has an order of magnitude reduction in computational scaling than the conventional scheme. With a few pilot numerical examples, we would demonstrate that this scheme can achieve very high accuracy with significant savings in computational time.

Accelerating coupled cluster calculations with nonlinear dynamics and supervised machine learning.

In this paper, the iteration scheme associated with single reference coupled cluster theory has been analyzed using nonlinear dynamics. The phase space analysis indicates the presence of a few significant cluster amplitudes, mostly involving valence excitations, that dictate the dynamics, while all other amplitudes are enslaved. Starting with a few initial iterations to establish the inter-relationship among the cluster amplitudes, a supervised machine learning scheme with a polynomial kernel ridge regression model has been employed to express each of the enslaved amplitudes uniquely in terms of the former set of amplitudes. The subsequent coupled cluster iterations are restricted solely to determine those significant excitations, and the enslaved amplitudes are determined through the already established functional mapping. We will show that our hybrid scheme leads to a significant reduction in the computational time without sacrificing the accuracy.

Total Synthesis of the Meroterpenoid Manginoid†A as Fueled by a Challenging Pinacol Coupling and Bicycle-forming Etherification.

The manginoids are a unique collection of bioactive natural products whose structures fuse an oxa-bridged spirocyclohexanedione with a heavily substituted trans-hydrindane framework. Herein, we show that such architectures can be accessed through a strategy combining a challenging pinacol coupling and bicycle-forming etherification with several additional chemo- and regioselective reactions. The success of these key events proved to be highly substrate and condition specific, affording insights for their application to other targets. As a result, not only has a 19-step total synthesis of manginoid A been achieved, but a potential roadmap to access other members of the family and related natural products has also been identified.

Stability analysis of a double similarity transformed coupled cluster theory.

In this paper, we have analyzed the time series associated with the iterative scheme of a double similarity transformed coupled cluster theory. The coupled iterative scheme to solve the ground state Schrödinger equation is cast as a multivariate time-discrete map, and the solutions show the universal Feigenbaum dynamics. Using recurrence analysis, it is shown that the dynamics of the iterative process is dictated by a small subgroup of cluster operators, mostly those involving chemically active orbitals, whereas all other cluster operators with smaller amplitudes are enslaved. Using synergetics, we will indicate how the master-slave dynamics can suitably be exploited to develop a novel coupled-cluster algorithm in a much reduced dimension.

Automatic State Interaction with Large Localized Active Spaces for Multimetallic Systems.

The localized active space self-consistent field method factorizes a complete active space wave function into an antisymmetrized product of localized active space wave function fragments. Correlation between fragments is then reintroduced through localized active space state interaction (LASSI), in which the Hamiltonian is diagonalized in a model space of LAS states. However, the optimal procedure for defining the LAS fragments and LASSI model space is unknown. We here present an automated framework to explore systematically convergent sets of model spaces, which we call LASSI[r, q]. This method requires the user to select only r, the number of electron hops from one fragment to another, and q, the number of fragment basis functions per Hilbert space, which converges to CASCI in the limit of r, q → ∞. Numerical tests of this method on the trimetal oxo-centered complexes [Fe(III)Al(III)Fe(II)(µ3-O)(HCOO)6] and [Fe(III)2Fe(II)(µ3-O)(HCOO)6] show efficient convergence to the CASCI limit with 4-10 orders of magnitude fewer states than CASCI.

Periodic Density Matrix Embedding for CO Adsorption on the MgO(001) Surface.

The adsorption of simple gas molecules to metal oxide surfaces is a primary step in many heterogeneous catalysis applications. Quantum chemical modeling of these reactions is a challenge in terms of both cost and accuracy, and quantum-embedding methods are promising, especially for localized chemical phenomena. In this work, we employ density matrix embedding theory (DMET) for periodic systems to calculate the adsorption energy of CO to the MgO(001) surface. Using coupled-cluster theory with single and double excitations and second-order Møller-Plesset perturbation theory as quantum chemical solvers, we perform calculations with embedding clusters up to 266 electrons in 306 orbitals, with the largest embedding models agreeing to within 1.2 kcal/mol of the non-embedding references. Moreover, we present a memory-efficient procedure of storing and manipulating electron repulsion integrals in the embedding space within the framework of periodic DMET.