The Journal of Chemical Physics

Ab initio density matrix renormalization group study of magnetic coupling in dinuclear iron and chromium complexes

Travis V. Harris, Yuki Kurashige, Takeshi Yanai, Keiji Morokuma

Published Feb 3, 2014• 45 citations

Abstract

The applicability of ab initio multireference wavefunction-based methods to the study of magnetic complexes has been restricted by the quickly rising active-space requirements of oligonuclear systems and dinuclear complexes with S > 1 spin centers. Ab initio density matrix renormalization group (DMRG) methods built upon an efficient parameterization of the correlation network enable the use of much larger active spaces, and therefore may offer a way forward. Here, we apply DMRG-CASSCF to the dinuclear complexes [Fe2OCl6](2-) and [Cr2O(NH3)10](4+). After developing the methodology through systematic basis set and DMRG M testing, we explore the effects of extended active spaces that are beyond the limit of conventional methods. We find that DMRG-CASSCF with active spaces including the metal d orbitals, occupied bridging-ligand orbitals, and their virtual double shells already capture a major portion of the dynamic correlation effects, accurately reproducing the experimental magnetic coupling constant (J) of [Fe2OCl6](2-) with (16e,26o), and considerably improving the smaller active space results for [Cr2O(NH3)10](4+) with (12e,32o). For comparison, we perform conventional MRCI+Q calculations and find the J values to be consistent with those from DMRG-CASSCF. In contrast to previous studies, the higher spin states of the two systems show similar deviations from the Heisenberg spectrum, regardless of the computational method.

DOI: 10.1063/1.4863345

Related papers

High-performance ab initio density matrix renormalization group method: Applicability to large-scale multireference problems for metal compounds

Yuki Kurashige, Takeshi Yanai

This article presents an efficient and parallelized implementation of the density matrix renormalization group (DMRG) algorithm for quantum chemistry calculations. The DMRG method as a large-scale multireference electronic structure model is by nature particularly efficient for one-dimensionally correlated systems, while the present development is oriented toward applications for polynuclear transition metal compounds, in which the macroscopic one-dimensional structure of electron correlation is absent. A straightforward extension of the DMRG algorithm is proposed with further improvements and aggressive optimizations to allow its application with large multireference active space, which is often demanded for metal compound calculations. Special efficiency is achieved by making better use of sparsity and symmetry in the operator and wave function representations. By accomplishing computationally intensive DMRG calculations, the authors have found that a large number of renormalized basis states are required to represent high entanglement of the electron correlation for metal compound applications, and it is crucial to adopt auxiliary perturbative correction to the projected density matrix during the DMRG sweep optimization in order to attain proper convergence to the solution. Potential energy curve calculations for the Cr(2) molecule near the known equilibrium precisely predicted the full configuration interaction energies with a correlation space of 24 electrons in 30 orbitals [denoted by (24e,30o)]. The energies are demonstrated to be accurate to 0.6mE(h) (the error from the extrapolated best value) when as many as 10,000 renormalized basis states are employed for the left and right DMRG block representations. The relative energy curves for [Cu(2)O(2)](2+) along the isomerization coordinate were obtained from DMRG and other correlated calculations, for which a fairly large orbital space (32e,62o) is modeled as a full correlation space. The DMRG prediction nearly overlaps with the energy curve from the coupled cluster with singles, doubles, and perturbative triple [CCSD(T)] calculations, while the multireference complete active space self-consistent field (CASSCF) calculations with the small reference configuration (8e,8o) are found to overestimate the biradical character of the electronic state of [Cu(2)O(2)](2+) according to the one-electron density matrix analysis.

The Journal of Chemical PhysicsJun 19, 2009193 citations

Calculation of NMR Shieldings for Strongly Correlated Molecules Using the Density Matrix Renormalization Group Self-Consistent-Field Method with Large Active Spaces

Jiaao Su, Wei Huang, Yinxuan Song +1 more

Accurate prediction of nuclear magnetic resonance (NMR) shielding constants for strongly correlated systems remains challenging. In this context, we present a density matrix renormalization group self-consistent-field (DMRG-SCF) implementation for NMR shielding calculations with large active spaces. Within a gauge-including atomic orbital DMRG-SCF (GIAO-DMRG-SCF) framework, we formulate and solve the coupled-perturbed DMRG-SCF (CP-DMRG-SCF) equations using local-site effective Hamiltonians. This approach enables the direct evaluation of the wave function's response to an external magnetic field. In the limit of small active spaces, the computed NMR shieldings are in agreement with conventional complete active space self-consistent-field (CASSCF) results, while the present approach naturally extends to larger active spaces inaccessible to standard multiconfigurational methods. Applications to a series of molecular systems, particularly strongly correlated transition metal complexes, demonstrate that the proposed implementation can be both reliable and accurate.

Journal of Chemical Theory and ComputationSep 2, 20260 citations

Second-order perturbation theory with a density matrix renormalization group self-consistent field reference function: Theory and application to the study of chromium dimer

Yuki Kurashige, Takeshi Yanai

We present a second-order perturbation theory based on a density matrix renormalization group self-consistent field (DMRG-SCF) reference function. The method reproduces the solution of the complete active space with second-order perturbation theory (CASPT2) when the DMRG reference function is represented by a sufficiently large number of renormalized many-body basis, thereby being named DMRG-CASPT2 method. The DMRG-SCF is able to describe non-dynamical correlation with large active space that is insurmountable to the conventional CASSCF method, while the second-order perturbation theory provides an efficient description of dynamical correlation effects. The capability of our implementation is demonstrated for an application to the potential energy curve of the chromium dimer, which is one of the most demanding multireference systems that require best electronic structure treatment for non-dynamical and dynamical correlation as well as large basis sets. The DMRG-CASPT2/cc-pwCV5Z calculations were performed with a large (3d double-shell) active space consisting of 28 orbitals. Our approach using large-size DMRG reference addressed the problems of why the dissociation energy is largely overestimated by CASPT2 with the small active space consisting of 12 orbitals (3d4s), and also is oversensitive to the choice of the zeroth-order Hamiltonian.

The Journal of Chemical PhysicsSep 2, 2011312 citations

Density matrix renormalization group for ab initio Calculations and associated dynamic correlation methods: A review of theory and applications

Takeshi Yanai, Yuki Kurashige, Wataru Mizukami +3 more

The recent advent of the density matrix renormalization group (DMRG) theory has delivered a new capability to compute multireference (MR) wave function with large configuration space, which far exceeds the limitation of conventional approaches. Here, we provide an overview of our recent work on the developments of ab initio DMRG methods in the context of the active space approaches and their applications to MR chemical systems. © 2014 Wiley Periodicals, Inc.

International Journal of Quantum ChemistryOct 23, 2014172 citations

Density matrix renormalization group (DMRG) method as a common tool for large active-space CASSCF/CASPT2 calculations

Naoki Nakatani, Sheng Guo

This paper describes an interface between the density matrix renormalization group (DMRG) method and the complete active-space self-consistent field (CASSCF) method and its analytical gradient, as well as an extension to the second-order perturbation theory (CASPT2) method. This interfacing allows large active-space multi-reference computations to be easily performed. The interface and its extension are both implemented in terms of reduced density matrices (RDMs) which can be efficiently computed via the DMRG sweep algorithm. We also present benchmark results showing that, in practice, the DMRG-CASSCF calculations scale with active-space size in a polynomial manner in the case of quasi-1D systems. Geometry optimization of a binuclear iron-sulfur cluster using the DMRG-CASSCF analytical gradient is demonstrated, indicating that the inclusion of the valence p-orbitals of sulfur and double-shell d-orbitals of iron lead to non-negligible changes in the geometry compared to the results of small active-space calculations. With the exception of the selection of M values, many computational settings in these practical DMRG calculations have been tuned and black-boxed in our interface, and so the resulting DMRG-CASSCF and DMRG-CASPT2 calculations are now available to novice users as a common tool to compute strongly correlated electronic wavefunctions.

The Journal of Chemical PhysicsMar 1, 201791 citations

Ab initio Density Matrix Renormalization Group Theory for Multireference Quantum Chemistry

Takeshi Yanai

A major challenge in molecular electronic structure theory arises from chemical systems where the mean-field picture based upon a single electronic configuration fails to properly describe mechanics of many electrons. Such complex electronic states should be accounted for using a correlated superposition of multiple determinants and are referred to as multireference in quantum chemistry. Application domain of traditional multireference methods is severely limited to systems containing a small number of correlated electrons because of their exponential complexity. In recent years, much attention has been drawn to the density matrix renormalization group (DMRG) theory as a promising alternative. It has been introduced in ab initio quantum chemistry calculations and shown to be a highly-scalable multireference method that can handle full quantum degrees of freedom for the strongly-correlated wavefunction with unprecedented large configuration space. In this article, we provide an explicative introduction to the underlying theory and formalism built into the DMRG method and give a brief overview of our methodological extensions that combine DMRG and active-space correlation models.

Molecular ScienceJan 1, 20140 citations