Molecular Physics

Excitation spectra of fully correlated donor-acceptor complexes by density matrix renormalisation group

Gergely Barcza, Anton Pershin, Ádám Gali, Örs Legeza

Published Oct 7, 2022• 2 citations

Abstract

In this paper, we investigate the convergence properties of density matrix renormalisation group (DMRG) calculations for vertical electronic excitations performed on three bimolecular complexes motivated by Szalay's recent benchmark results [JCTC. 2020;16:7]. Besides the high-level coupled-cluster reference spectrum, the extrapolated truncation-free DMRG results are compared against several alternative solutions based on multireference configuration interaction and coupled-cluster theories. By taking advantage of the error cancellation effects, which are inherent in the state-averaged computational scheme, we demonstrate that the extrapolated solution based on only low bond-dimensional DMRG data could already provide quantitative predictions even for large active spaces. We also show that orbital optimisation improves further the accuracy of the DMRG results which systematically approach the excitation energies of coupled cluster single-double-triple at similar computational costs.

DOI: 10.1080/00268976.2022.2130834

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

Computing Exchange Coupling Constants in Transition Metal Complexes with Tensor Product Selected Configuration Interaction

Arnab Bachhar, Nicholas J. Mayhall

Transition metal complexes present significant challenges for electronic structure theory due to strong electron correlation arising from partially filled d-orbitals. We compare our recently developed Tensor Product Selected Configuration Interaction (TPSCI) with Density Matrix Renormalization Group (DMRG) for computing exchange coupling constants in six transition metal systems, including dinuclear Cr, Fe, and Mn complexes and a tetranuclear Ni-cubane. TPSCI uses a locally correlated tensor product state basis to capture electronic structure efficiently while maintaining interpretability. From calculations on active spaces ranging from (22e,29o) to (42e,49o), we find that TPSCI consistently yields higher variational energies than DMRG due to truncation of local cluster states, but provides magnetic exchange coupling constants (J) generally within 10–30 cm –1 of DMRG results. Key advantages include natural multistate capability enabling direct J extrapolation with smaller statistical errors, and computational efficiency for challenging systems. However, cluster state truncation represents a fundamental limitation requiring careful convergence testing, particularly for large local cluster dimensions. We identify specific failure cases where current truncation schemes break down, highlighting the need for improved cluster state selection methods and distributed memory implementations to realize TPSCI’s full potential for strongly correlated systems.

arXiv:2508.13002Journal of Chemical Theory and ComputationFeb 4, 20260 citations

Multireference electron correlation methods with density matrix renormalisation group reference functions

Yuki Kurashige

Recent advances in quantum chemical density matrix renormalisation group (DMRG) theory are presented. The DMRG, originally devised as an alternative to the exact diagonalisation in condensed matter physics, has become a powerful quantum chemical method for molecular systems that exhibit a multireference character, e.g., excited states, π-conjugated systems, transition metal complexes, and in particular for large systems by combining it with conventional multireference electron correlation methods. The capability of the current quantum chemical DMRG is demonstrated for an application involving the potential energy curve of the chromium dimer, which is one of the most demanding multireference systems and thus requires the best electronic structure treatment for non-dynamical and dynamical correlation as well as large basis sets.

Molecular PhysicsSep 16, 201392 citations

DMRG control using an automated Richardson-type error protocol

Konrad Heinrich Marti, Markus Reiher

The density matrix renormalisation group (DMRG) algorithm solves the many-particle electronic Schrödinger equation as accurately as possible in a full-configuration-interaction (FCI) sense for a given one-electron basis. The FCI limit is approached by increasing the number of many-particle DMRG states which span the Hilbert space of interest. In this work we investigate an automated DMRG error protocol which extrapolates the electronic energy using Richardson's deferred approach to the limit. The key idea is to consider the result of a DMRG algorithm as an analytical function of an adjustable parameter like the number of DMRG states. We can then probe this analytical function by performing calculations for different set sizes of renormalised DMRG basis states. None of these calculations has to actually provide the desired accuracy but after we have collected enough information about the function's behaviour, we can represent it by an analytic rational function that may then be used to extrapolate to the converged energy. The advantage of this approach is that it delivers an error estimate for the electronic energy which can be used for accuracy control and as a convergence criterion. In addition, one may detect convergence to local energy minima by automatically increasing the analytic parameter (i.e. the number of renormalised DMRG states). The error estimates also allows us to aim at relative energies from extrapolated total energies of similar accuracy that is then independent of the number of renormalised DMRG states or of the truncation error of the reduced density matrix.

Molecular PhysicsFeb 10, 201025 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

State-of-the-art density matrix renormalization group and coupled cluster theory studies of the nitrogen binding curve

Garnet Kin‐Lic Chan, Mihály Kállay, Jürgen Gauß

We study the nitrogen binding curve with the density matrix renormalization group (DMRG) and single-reference and multireference coupled cluster (CC) theory. Our DMRG calculations use up to 4000 states and our single-reference CC calculations include up to full connected hextuple excitations. Using the DMRG, we compute an all-electron benchmark nitrogen binding curve, at the polarized, valence double-zeta level (28 basis functions), with an estimated accuracy of 0.03 mEh. We also assess the performance of more approximate DMRG and CC theories across the nitrogen curve. We provide an analysis of the relative strengths and merits of the DMRG and CC theory under different correlation conditions.

The Journal of Chemical PhysicsSep 27, 2004203 citations