Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Chemical Theory and Computation

DMRG for Multi-Component Quantum Systems: A Symmetry-Aware and Unified Framework

Xinyu Sun, Haibo Ma

Published Sep 23, 2026• 0 citations
Publisher

Abstract

Recent advances in high-resolution spectroscopy and cavity quantum technologies have driven growing interest in multicomponent quantum chemistry, where electronic degrees of freedom are treated on an equal footing with other quantum particles. Here, we present a unified symmetry-aware density matrix renormalization group (DMRG) framework for multicomponent quantum chemistry. By decoupling physical model specification from tensor-network representation and numerical optimization, new multicomponent models are introduced through symmetry quantum numbers, local basis, and operators, without modifying the underlying tensor algorithms. Built upon a mature tensor infrastructure, the framework provides a reusable symmetry-aware numerical engine, allowing new models to be implemented through concise Python-level specifications while fully reusing the optimized computational backend. The framework is demonstrated through symmetry-adapted nuclear-electronic orbital (NEO) and cavity quantum electrodynamics (CQED) DMRG calculations, yielding accurate electronic structures, excited-state energies, and polaritonic properties. This work establishes a reusable and extensible computational framework for symmetry-aware multicomponent quantum chemistry.

DOI
10.1021/acs.jctc.6c01338

Related papers

DMRG for multicomponent quantum chemistry

Mikuláš Matoušek

One of the more modern directions in quantum chemistry, which is mostly unexplored, is using Hamiltonians which include not only the electrons of the molecule, but also other quantum particles. These are also termed multicomponent Hamiltonians. The Density Matrix Renormalization Group method has become one of the mainline methods in multireference quantum chemistry. In this work we describe a way to extend the DMRG algorithm to use multicomponent Hamiltonians. We start by a description of the DMRG algorithm, showcasing the usage of quantum chemical DMRG on multiple challenging systems. We include a brief description how to extend the DMRG algorithm to also treat dynamical correlation with the Adiabatic Connection method. We then describe the extension of DMRG to multicomponent Hamiltonians which has been implemented by the author of this work, and finally show how DMRG can be useful for polaritonic chemistry, where the molecule interacts with the quantized states of light in an optical cavity.

•Jan 1, 2025•0 citations

The software landscape for the density matrix renormalization group

Per Sehlstedt, Jan Brandejs, Paolo Bientinesi +1 more

The density matrix renormalization group (DMRG) algorithm is a cornerstone computational method for studying quantum many-body systems, renowned for its accuracy and adaptability. Because DMRG provides a general framework applicable across various fields such as materials science, quantum chemistry, and quantum computing, one might expect a shared, flexible library to serve most users. Nevertheless, numerous independent implementations continue to appear, resulting in significant duplication of effort. To identify collaboration opportunities that can promote a more unified approach, we map the rapidly expanding DMRG software landscape and provide a comprehensive comparison of features across 37 existing packages. When comparing key features, such as parallelism strategies for high-performance computing and symmetry-adapted formulations that enhance efficiency, we found significant overlap among the packages. This overlap suggests opportunities for collaboration to modularize common functionality—e.g., tensor operations, symmetry representations, and eigensolvers—as the packages are mostly independent and share few third-party library dependencies. More collaboration on modularization could reduce duplication of effort, improve interoperability, and enable prioritization and quicker spread of new advances. We believe the current lack of modularity is more socially driven than a technical issue; hence, we see raising awareness about the existing implementations as a first step in the right direction. Ultimately, this work emphasizes the value of greater cohesion through modularity, which would benefit DMRG software and related tensor-network-centered software, enabling the solution of more complex and ambitious problems.

arXiv:2506.12629•Computer Physics Communications•Mar 20, 2026•1 citation

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 Physics•Jun 19, 2009•193 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 Science•Jan 1, 2014•0 citations

Nuclear-electronic all-particle density matrix renormalization group

Andrea Muolo, Alberto Baiardi, Robin Feldmann +1 more

We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians’ positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H2 and H3+, and one larger system, namely, BH3.

arXiv:2003.04446•The Journal of Chemical Physics•May 26, 2020•32 citations

QCMaquis 4.0: Multipurpose Electronic, Vibrational, and Vibronic Structure and Dynamics Calculations with the Density Matrix Renormalization Group

Kalman Szenes, Nina Glaser, Mihael Eraković +6 more

QCMaquis is a quantum chemistry software package for general molecular structure calculations in a matrix product state/matrix product operator formalism of the density matrix renormalization group (DMRG). It supports a wide range of features for electronic structure, multicomponent (pre-Born-Oppenheimer), anharmonic vibrational structure, and vibronic calculations. In addition to the ground and excited state solvers, QCMaquis allows for time propagation of matrix product states based on the tangent-space formulation of time-dependent DMRG. The latest developments include transcorrelated electronic structure calculations, very recent vibrational and vibronic models, and a convenient Python wrapper, facilitating the interface with external libraries. This paper reviews all the new features of QCMaquis and demonstrates them with new results.

arXiv:2505.01405•The Journal of Physical Chemistry A•Aug 1, 2025•6 citations
Recent advances in high-resolution spectroscopy and cavity quantum technologies have driven growing interest in multicomponent quantum chemistry, where electronic degrees of freedom are treated on an equal footing with other quantum particles. Here, we present a unified symmetry-aware density matrix renormalization group (DMRG) framework for multicomponent quantum chemistry. By decoupling physical model specification from tensor-network representation and numerical optimization, new multicomponent models are introduced through symmetry quantum numbers, local basis, and operators, without modifying the underlying tensor algorithms. Built upon a mature tensor infrastructure, the framework provides a reusable symmetry-aware numerical engine, allowing new models to be implemented through concise Python-level specifications while fully reusing the optimized computational backend. The framework is demonstrated through symmetry-adapted nuclear-electronic orbital (NEO) and cavity quantum electrodynamics (CQED) DMRG calculations, yielding accurate electronic structures, excited-state energies, and polaritonic properties. This work establishes a reusable and extensible computational framework for symmetry-aware multicomponent quantum chemistry.
One of the more modern directions in quantum chemistry, which is mostly unexplored, is using Hamiltonians which include not only the electrons of the molecule, but also other quantum particles. These are also termed multicomponent Hamiltonians. The Density Matrix Renormalization Group method has become one of the mainline methods in multireference quantum chemistry. In this work we describe a way to extend the DMRG algorithm to use multicomponent Hamiltonians. We start by a description of the DMRG algorithm, showcasing the usage of quantum chemical DMRG on multiple challenging systems. We include a brief description how to extend the DMRG algorithm to also treat dynamical correlation with the Adiabatic Connection method. We then describe the extension of DMRG to multicomponent Hamiltonians which has been implemented by the author of this work, and finally show how DMRG can be useful for polaritonic chemistry, where the molecule interacts with the quantized states of light in an optical cavity.
The density matrix renormalization group (DMRG) algorithm is a cornerstone computational method for studying quantum many-body systems, renowned for its accuracy and adaptability. Because DMRG provides a general framework applicable across various fields such as materials science, quantum chemistry, and quantum computing, one might expect a shared, flexible library to serve most users. Nevertheless, numerous independent implementations continue to appear, resulting in significant duplication of effort. To identify collaboration opportunities that can promote a more unified approach, we map the rapidly expanding DMRG software landscape and provide a comprehensive comparison of features across 37 existing packages. When comparing key features, such as parallelism strategies for high-performance computing and symmetry-adapted formulations that enhance efficiency, we found significant overlap among the packages. This overlap suggests opportunities for collaboration to modularize common functionality—e.g., tensor operations, symmetry representations, and eigensolvers—as the packages are mostly independent and share few third-party library dependencies. More collaboration on modularization could reduce duplication of effort, improve interoperability, and enable prioritization and quicker spread of new advances. We believe the current lack of modularity is more socially driven than a technical issue; hence, we see raising awareness about the existing implementations as a first step in the right direction. Ultimately, this work emphasizes the value of greater cohesion through modularity, which would benefit DMRG software and related tensor-network-centered software, enabling the solution of more complex and ambitious problems.
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.
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.
We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians’ positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H2 and H3+, and one larger system, namely, BH3.
QCMaquis is a quantum chemistry software package for general molecular structure calculations in a matrix product state/matrix product operator formalism of the density matrix renormalization group (DMRG). It supports a wide range of features for electronic structure, multicomponent (pre-Born-Oppenheimer), anharmonic vibrational structure, and vibronic calculations. In addition to the ground and excited state solvers, QCMaquis allows for time propagation of matrix product states based on the tangent-space formulation of time-dependent DMRG. The latest developments include transcorrelated electronic structure calculations, very recent vibrational and vibronic models, and a convenient Python wrapper, facilitating the interface with external libraries. This paper reviews all the new features of QCMaquis and demonstrates them with new results.