The Effect of Geometry, Spin and Orbital Optimization in Achieving Accurate, Fully-Correlated Results for Iron-Sulfur Cubanes
arXiv:2105.01754 · doi:10.1021/acs.jctc.1c00830
Abstract
Iron-sulfur clusters comprise an important functional motif of the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from fully-correlated calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. To provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of cubane-based clusters using two complementary, fully-correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (< 1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the lowest singlet and the highest spin states is strongly geometry dependent, changing by an order of magnitude upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap by a factor of two compared to typical mean-field orbital approaches. Our results clearly demonstrate the need for performing highly correlated calculations to unveil the challenging electronic structure of these complex catalytic centers.
arXiv comment16 pages, 7 figures, 6 tables plus SI (12 pages, 3 figures, 9 tables)
References in corpus (25)
- The density-matrix renormalization group in the age of matrix product states
- NWChem: Past, Present, and Future
- Heat-bath Configuration Interaction: An efficient selected CI algorithm inspired by heat-bath sampling
- Orbital Optimization in the Density Matrix Renormalization Group, with applications to polyenes and β-carotene
- The Density Matrix Renormalization Group in Chemistry and Molecular Physics: Recent Developments and New Challenges
- A deterministic alternative to the full configuration interaction quantum Monte Carlo method
- Low-energy spectrum of iron-sulfur clusters directly from many-particle quantum mechanics
- The density matrix renormalization group for ab initio quantum chemistry
- Tensor product methods and entanglement optimization for ab initio quantum chemistry
- A spin-adapted Density Matrix Renormalization Group algorithm for quantum chemistry
- Controlling the accuracy of the density matrix renormalization group method: The Dynamical Block State Selection approach
- Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method
- CheMPS2: a free open-source spin-adapted implementation of the density matrix renormalization group for ab initio quantum chemistry
- CASSCF with Extremely Large Active Spaces using the Adaptive Sampling Configuration Interaction Method
- The Ground State Electronic Energy of Benzene
- Fast Semistochastic Heat-Bath Configuration Interaction
- The electronic complexity of the ground-state of the FeMo cofactor of nitrogenase as relevant to quantum simulations
- Orbital entanglement in quantum chemistry
- Spin-adapted Matrix Product States and Operators
- Electronic landscape of the P-cluster of nitrogenase as revealed through many-electron quantum wavefunctions
- DMRG-SCF study of the singlet, triplet, and quintet states of oxo-Mn(Salen)
- Fermionic orbital optimisation in tensor network states
- Investigation of metal-insulator like transition through the ab initio density matrix renormalization group approach
- Are multi-quasiparticle interactions important in molecular ionization?
- Dynamical Mean-Field Theory Simulations with the Adaptive Sampling Configuration Interaction Method