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

arXiv:2609.08838

Orbital choice in constructing model Hamiltonians

David Wilian Oliveira de Sousa, Elvira R. Sayfutyarova

Published Sep 8, 2026• 0 citations
arXivHTMLPDF

Abstract

Model Hamiltonians (e.g., the Hubbard or Heisenberg models) provide a simple yet physically meaningful description of complex electronic phenomena through a small number of parameters, such as electron hopping integrals or magnetic exchange couplings. Their construction typically exploits the local nature of electron correlation effects and therefore relies on localized orbitals. Although many localization schemes are available, the resulting localized orbitals depend not only on the chosen localization functional but also on the orbital space to which the localization procedure is applied. These choices are often treated as technical details and are rarely discussed explicitly. Here, we demonstrate that different, chemically reasonable choices of localized orbitals can lead to substantially different model Hamiltonian parameters, even when they produce nearly identical electronic energies. Using Density Matrix Downfolding, we systematically investigate how the localization functional and, more importantly, the choice of model orbital space affect effective Hamiltonians derived from ab initio calculations for representative $π$-conjugated systems and transition metal complexes.

DOI
10.48550/arxiv.2609.08838

Related papers

New electron correlation theories for transition metal chemistry

Konrad Heinrich Marti, Markus Reiher

Electronic structure theory faces many computational challenges in transition metal chemistry. Usually, density functional theory is the method of choice for theoretical studies on transition metal complexes and clusters mostly because it is the only feasible one, although its results are not systematically improvable. By contrast, multireference ab initio methods could provide a correct description of the electronic structure, but are limited to small molecules because of the tremendous computational resources required. In recent years, conceptually new ab initio methods emerged that turned out to be promising for theoretical coordination chemistry. We review and discuss two efficient parametrization schemes for the electronic wave function, the matrix product states and the complete-graph tensor network states. Their advantages are demonstrated at example transition metal complexes. Especially, tensor network states might provide the key to accurately describe strongly correlated and magnetic molecular systems in transition metal chemistry.

Physical Chemistry Chemical Physics•Jan 1, 2011•130 citations

Orbital Entanglement Analysis of Exchange-Coupled Systems

Christopher Johannes Stein, Dimitrios A. Pantazis, Vera Krewald

A new tool for the interpretation of multiconfigurational wave functions representing the spin states of exchange-coupled transition metal complexes is introduced. Based on orbital entanglement measures, herein derived from multiconfigurational density matrix renormalization group calculations, the complexity of the wave function is reduced, thus facilitating a connection with established concepts for the interpretation of magnetically coupled systems. We show that the entanglement of localized orbitals with a small basis set is a good representation of the magnetic coupling topology and that it is sensitive to chemical changes in homologous complexes. Furthermore, we introduce a measure for the magnetic relevance of orbitals in the active subspace and a concept for the quantitative comparison of different chemical species. The approach presented here will be easily applicable to higher nuclearity clusters, providing a direct insight into all states of the Heisenberg spin ladder for systems previously accessible only by single-configurational methods.

The Journal of Physical Chemistry Letters•Oct 15, 2019•23 citations

Noneigenstate Ab Initio Density Matrix Downfolding for Constructing Model Hamiltonians in Quantum Chemistry

David Wilian Oliveira de Sousa, Elvira R. Sayfutyarova

Model Hamiltonians represent a convenient way of reducing complex problems of many-electron quantum mechanics to much simpler problems: they can fully reproduce the core behaviors of a system of interest by encoding only the dominant physical interactions and using only a small number of associated parameters. Model Hamiltonians have been successfully applied to describe many chemical and physical phenomena. Density matrix downfolding (DMD) [ J. Chem. Phys. 2015, 143 (10), 102814 ] allows the derivation of model Hamiltonians of any form in a systematically improvable fashion by matching the energy spectrum of ab initio Hamiltonians with those of the model Hamiltonians. This method allows not only the improvement of existing models but also the construction of accurate and efficient physical models for various systems. While DMD looks like a promising approach, it has rarely been applied within chemistry, and neither its limits nor its practical performance is well-understood. In this work, we evaluated the performance of DMD, based on noneigenstates of ab initio Hamiltonians, for several realistic chemical systems: benzene, naphthalene, FeSe, and a prototypical Fe(IV)═O complex found in the active sites of 2-oxoglutarate-dependent oxygenases. Our results show that DMD is a reliable and computationally efficient tool for obtaining optimized model Hamiltonians in quantum chemistry. This not only opens the door to studying complex systems at reduced computational cost but also to isolating and understanding the physical core principles that dominate their behavior─this might offer new insights for tuning or even designing chemical systems for applications ranging from biochemistry to catalysis.

Journal of Chemical Theory and Computation•Aug 28, 2025•0 citations

A closed local-orbital unified description of DFT and many-body effects

Fernando Flores, Diego Soler‐Polo, José Ortega

Abstract Density functional theory (DFT) is usually formulated in terms of the electron density as a function of position n(r). Here we discuss an alternative formulation of DFT in terms of the orbital occupation numbers {n α } associated with a local-orbital orthonormal basis set {ϕ α }. First, we discuss how the building blocks of DFT, namely the Hohenberg–Kohn theorems, the Levy–Lieb approach and the Kohn–Sham method, can be adapted for a description in terms of {n α }. In particular, the total energy is now a function of {n α }, E[{n α }], and a Kohn–Sham-like Hamiltonian is derived introducing the effects of the electron–electron interactions via effective potentials, { V α eff = ∂ E e e [ { n β } ] / ∂ n α } . In a second step we consider the Hartree and exchange energies and discuss how to describe them, in the spirit of a DFT approach, in terms of the orbital occupation numbers. In this contribution special attention is paid to the description of the (intra-atomic) correlation energy and corresponding correlation potentials {V corr,α }. For this purpose, a model system is analyzed in detail, whereby an atomic Hamiltonian interacts with the environment via a simplified model; the use of this model allows us to obtain the correlation energy and potentials (in terms of {n α }) for different cases corresponding to low, intermediate and high electron correlations.

Journal of Physics Condensed Matter•May 11, 2022•3 citations

Non-eigenstate Ab Initio Density Matrix Downfolding to produce Model Hamiltonians in Quantum Chemistry

David Wilian Oliveira de Sousa, Elvira R. Sayfutyarova

Model Hamiltonians represent a convenient way of reducing complex problems of many-electron quantum mechanics to much simpler problems: they can fully reproduce the core behaviors of a system of interest by encoding only the dominant physical interactions and using only a small number of associated parameters. Model Hamiltonians have been successfully applied to describe many chemical and physical phenomena. Density Matrix Downfolding (DMD) [J. Chem. Phys. 2015, 143 (10), 102814] allows the derivation of model Hamiltonians, of any form, in a systematically improvable fashion, by matching the energy spectrum of ab initio Hamiltonians with those of the model Hamiltonians. This method allows not only the improvement of existing models but also the construction of accurate and efficient physical models for various systems. While DMD looks like a promising approach, it has rarely been applied within chemistry, and neither its limits nor practical performance are well understood. In this work we evaluated the performance of DMD, based on non-eigenstates of ab initio Hamiltonians, for several realistic chemical systems: benzene, naphthalene, FeSe, and a prototypical Fe(IV)=O complex found in the active sites of 2-oxoglutarate-dependent oxygenases. Our results show that DMD is a reliable and computationally efficient tool for obtaining optimized model Hamiltonians in quantum chemistry. This not only opens the door to studying complex systems at reduced computational cost but also to isolating and understanding the physical core principles that dominate their behavior — this might offer new insights for tuning, or even designing, chemical systems for applications ranging from biochemistry to catalysis.

•Apr 21, 2025•0 citations

Effective Hamiltonians derived from equation-of-motion coupled-cluster wave functions: Theory and application to the Hubbard and Heisenberg Hamiltonians

Pavel Pokhilko, Anna I. Krylov

Effective Hamiltonians, which are commonly used for fitting experimental observables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not trivial. In this contribution, we apply Bloch's formalism to equation-of-motion coupled-cluster wave functions to rigorously derive effective Hamiltonians in Bloch's and des Cloizeaux's forms. We report the key equations and illustrate the theory by application to systems with two or three unpaired electrons, which give rise to electronic states of covalent and ionic characters. We show that Hubbard's and Heisenberg's Hamiltonians can be extracted directly from the so-obtained effective Hamiltonians. By establishing a quantitative connection between many-body states and simple models, the approach facilitates the analysis of the correlated wave functions. We propose a simple diagnostic for assessing the validity of the model space choice based on the overlaps between the target- and model-space states. Artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed.

The Journal of Chemical Physics•Mar 3, 2020•27 citations
Model Hamiltonians (e.g., the Hubbard or Heisenberg models) provide a simple yet physically meaningful description of complex electronic phenomena through a small number of parameters, such as electron hopping integrals or magnetic exchange couplings. Their construction typically exploits the local nature of electron correlation effects and therefore relies on localized orbitals. Although many localization schemes are available, the resulting localized orbitals depend not only on the chosen localization functional but also on the orbital space to which the localization procedure is applied. These choices are often treated as technical details and are rarely discussed explicitly. Here, we demonstrate that different, chemically reasonable choices of localized orbitals can lead to substantially different model Hamiltonian parameters, even when they produce nearly identical electronic energies. Using Density Matrix Downfolding, we systematically investigate how the localization functional and, more importantly, the choice of model orbital space affect effective Hamiltonians derived from ab initio calculations for representative πππ-conjugated systems and transition metal complexes.
Electronic structure theory faces many computational challenges in transition metal chemistry. Usually, density functional theory is the method of choice for theoretical studies on transition metal complexes and clusters mostly because it is the only feasible one, although its results are not systematically improvable. By contrast, multireference ab initio methods could provide a correct description of the electronic structure, but are limited to small molecules because of the tremendous computational resources required. In recent years, conceptually new ab initio methods emerged that turned out to be promising for theoretical coordination chemistry. We review and discuss two efficient parametrization schemes for the electronic wave function, the matrix product states and the complete-graph tensor network states. Their advantages are demonstrated at example transition metal complexes. Especially, tensor network states might provide the key to accurately describe strongly correlated and magnetic molecular systems in transition metal chemistry.
A new tool for the interpretation of multiconfigurational wave functions representing the spin states of exchange-coupled transition metal complexes is introduced. Based on orbital entanglement measures, herein derived from multiconfigurational density matrix renormalization group calculations, the complexity of the wave function is reduced, thus facilitating a connection with established concepts for the interpretation of magnetically coupled systems. We show that the entanglement of localized orbitals with a small basis set is a good representation of the magnetic coupling topology and that it is sensitive to chemical changes in homologous complexes. Furthermore, we introduce a measure for the magnetic relevance of orbitals in the active subspace and a concept for the quantitative comparison of different chemical species. The approach presented here will be easily applicable to higher nuclearity clusters, providing a direct insight into all states of the Heisenberg spin ladder for systems previously accessible only by single-configurational methods.
Model Hamiltonians represent a convenient way of reducing complex problems of many-electron quantum mechanics to much simpler problems: they can fully reproduce the core behaviors of a system of interest by encoding only the dominant physical interactions and using only a small number of associated parameters. Model Hamiltonians have been successfully applied to describe many chemical and physical phenomena. Density matrix downfolding (DMD) [ J. Chem. Phys. 2015, 143 (10), 102814 ] allows the derivation of model Hamiltonians of any form in a systematically improvable fashion by matching the energy spectrum of ab initio Hamiltonians with those of the model Hamiltonians. This method allows not only the improvement of existing models but also the construction of accurate and efficient physical models for various systems. While DMD looks like a promising approach, it has rarely been applied within chemistry, and neither its limits nor its practical performance is well-understood. In this work, we evaluated the performance of DMD, based on noneigenstates of ab initio Hamiltonians, for several realistic chemical systems: benzene, naphthalene, FeSe, and a prototypical Fe(IV)═O complex found in the active sites of 2-oxoglutarate-dependent oxygenases. Our results show that DMD is a reliable and computationally efficient tool for obtaining optimized model Hamiltonians in quantum chemistry. This not only opens the door to studying complex systems at reduced computational cost but also to isolating and understanding the physical core principles that dominate their behavior─this might offer new insights for tuning or even designing chemical systems for applications ranging from biochemistry to catalysis.
Abstract Density functional theory (DFT) is usually formulated in terms of the electron density as a function of position n(r). Here we discuss an alternative formulation of DFT in terms of the orbital occupation numbers {n α } associated with a local-orbital orthonormal basis set {ϕ α }. First, we discuss how the building blocks of DFT, namely the Hohenberg–Kohn theorems, the Levy–Lieb approach and the Kohn–Sham method, can be adapted for a description in terms of {n α }. In particular, the total energy is now a function of {n α }, E[{n α }], and a Kohn–Sham-like Hamiltonian is derived introducing the effects of the electron–electron interactions via effective potentials, { V α eff = ∂ E e e [ { n β } ] / ∂ n α } . In a second step we consider the Hartree and exchange energies and discuss how to describe them, in the spirit of a DFT approach, in terms of the orbital occupation numbers. In this contribution special attention is paid to the description of the (intra-atomic) correlation energy and corresponding correlation potentials {V corr,α }. For this purpose, a model system is analyzed in detail, whereby an atomic Hamiltonian interacts with the environment via a simplified model; the use of this model allows us to obtain the correlation energy and potentials (in terms of {n α }) for different cases corresponding to low, intermediate and high electron correlations.
Model Hamiltonians represent a convenient way of reducing complex problems of many-electron quantum mechanics to much simpler problems: they can fully reproduce the core behaviors of a system of interest by encoding only the dominant physical interactions and using only a small number of associated parameters. Model Hamiltonians have been successfully applied to describe many chemical and physical phenomena. Density Matrix Downfolding (DMD) [J. Chem. Phys. 2015, 143 (10), 102814] allows the derivation of model Hamiltonians, of any form, in a systematically improvable fashion, by matching the energy spectrum of ab initio Hamiltonians with those of the model Hamiltonians. This method allows not only the improvement of existing models but also the construction of accurate and efficient physical models for various systems. While DMD looks like a promising approach, it has rarely been applied within chemistry, and neither its limits nor practical performance are well understood. In this work we evaluated the performance of DMD, based on non-eigenstates of ab initio Hamiltonians, for several realistic chemical systems: benzene, naphthalene, FeSe, and a prototypical Fe(IV)=O complex found in the active sites of 2-oxoglutarate-dependent oxygenases. Our results show that DMD is a reliable and computationally efficient tool for obtaining optimized model Hamiltonians in quantum chemistry. This not only opens the door to studying complex systems at reduced computational cost but also to isolating and understanding the physical core principles that dominate their behavior — this might offer new insights for tuning, or even designing, chemical systems for applications ranging from biochemistry to catalysis.
Effective Hamiltonians, which are commonly used for fitting experimental observables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not trivial. In this contribution, we apply Bloch's formalism to equation-of-motion coupled-cluster wave functions to rigorously derive effective Hamiltonians in Bloch's and des Cloizeaux's forms. We report the key equations and illustrate the theory by application to systems with two or three unpaired electrons, which give rise to electronic states of covalent and ionic characters. We show that Hubbard's and Heisenberg's Hamiltonians can be extracted directly from the so-obtained effective Hamiltonians. By establishing a quantitative connection between many-body states and simple models, the approach facilitates the analysis of the correlated wave functions. We propose a simple diagnostic for assessing the validity of the model space choice based on the overlaps between the target- and model-space states. Artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed.