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

The Journal of Chemical Physics

Local expansion of N-representable one-particle density matrices yielding a prescribed electron density

Jānos Pipek, Szilvia Nagy

Published Oct 16, 2003• 2 citations
Publisher

Abstract

Multiresolution (or wavelet) analysis offers a strictly local basis set for a systematic introduction of new details into Hilbert space operators. Using this tool we have previously developed an expansion method for density matrices. The set of density operators providing a given electron density plays an essential role in density functional theory, in the minimization of energy expectation values with the constraint that the electron density is fixed. In this contribution, using multiresolution analysis, we present an excellent quality density matrix expansion yielding a prescribed electron density, and compare it to other known methods. Due to the strictly local nature of the applied basis functions, our construction has the specific advantage that the resulting density matrix is correlated and N-representable in the infinite resolution limit. As a further consequence of this scheme we can conclude that the deviation of the exact kinetic energy functional from the Weizsäcker term is not a necessary consequence of the particle statistics.

DOI
10.1063/1.1611176

Related papers

Quantum Mechanical Operators in Multiresolution Hilbert Spaces

Jānos Pipek, Theodore E. Simos, George Maroulis

Wavelet analysis, which is a shorthand notation for the concept of multiresolution analysis (MRA), becomes increasingly popular in high efficiency storage algorithms of complex spatial distributions. This approach is applied for describing wave functions of quantum systems. At any resolution level of MRA expansions a physical observable is represented by an infinite matrix which is “canonically” chosen as the projection of its operator in the Schrödinger picture onto the subspace of the given resolution. It is shown that this canonical choice is only a particular member of possible operator representations. Among these, there exits an optimal choice, usually different from the canonical one, which gives the best numerical values in eigenvalue problems. This construction works even in those cases, where the canonical definition is unusable. The commutation relation of physical operators is also studied in MRA subspaces. It is shown that the required commutation rules are satisfied in the fine resolution limit, whereas in coarse grained spaces a correction appears depending only on the representation of the momentum operator.

AIP conference proceedings•Jan 1, 2007•0 citations

Wavefunction Optimization at the Complete Basis Set Limit with Multiwavelets and DMRG

Martina Nibbi, Luca Frediani, Evgueni Dinvay +1 more

High Resolution Image Download MS PowerPoint Slide The density matrix renormalization group (DMRG) is a powerful numerical technique to solve strongly correlated quantum systems: it deals well with systems which are not dominated by a single configuration (unlike Coupled Cluster) and it converges rapidly to the Full Configuration Interaction (FCI) limit (unlike truncated Configuration Interaction (CI) expansions). In this work, we develop an algorithm integrating DMRG within the multiwavelet-based multiresolution analysis (MRA). Unlike fixed basis sets, multiwavelets offer an adaptive and hierarchical representation of functions, approaching the complete basis set limit to a specified precision. As a result, this combined technique leverages the multireference capability of DMRG and the complete basis set limit of MRA and multiwavelets. More specifically, we adopt a pre-existing Lagrangian optimization algorithm for orbitals represented in the MRA domain and improve its computational efficiency by replacing the original CI calculations with DMRG. Additionally, we substitute the reduced density matrices computation with the direct extraction of energy gradients from the DMRG tensors. We apply our method to small systems such H 2, He, HeH 2, BeH 2 and N 2 . The results demonstrate that our approach reduces the final energy while keeping the number of orbitals low compared to FCI calculations on an atomic orbital basis set.

arXiv:2503.10808•The Journal of Physical Chemistry A•Nov 13, 2025•1 citation

Diagrammatic multiresolution analysis for electron correlations

Heinz‐Jürgen Flad, Wolfgang Hackbusch, Hongjun Luo +1 more

We present a multiscale treatment of electron correlations based on hyperbolic wavelet expansions of Jastrow-type correlation functions. Wavelets provide hierarchical basis sets that can be locally adapted to the length- and energy-scales of physical phenomena. Combined with hyperbolic tensor products and local adaptive refinement near the interelectron cusp, these wavelet bases enable sparse representations of Jastrow factors. The computational efficiency of wavelets in electronic structure calculations is demonstrated within the coupled electron-pair approximation (local ansatz). Based on a diagrammatic multiresolution analysis, we discuss various kinds of sparsity features for matrix elements required by the local ansatz. Sparsity originates from the hierarchical structure and vanishing moments property of wavelet bases. This led us to a recurrence scheme for the evaluation of matrix elements with almost linear computational complexity with respect to the size of the underlying isotropic 3d-wavelet basis. Numerical studies for selected diagrams are presented for a homogeneous electron gas.

Physical Review B•Mar 18, 2005•15 citations

Scale‐adaptive tensor algebra for local many‐body methods of electronic structure theory

Dmitry I. Lyakh

While the formalism of multiresolution analysis, based on wavelets and adaptive integral representations of operators, is actively progressing in electronic structure theory (mostly on the independent‐particle level and, recently, second‐order perturbation theory), the concepts of multiresolution and adaptivity can also be utilized within the traditional formulation of correlated (many‐particle) theory based on second quantization and the corresponding (generally nonorthogonal) tensor algebra. In this article, we present a formalism called scale‐adaptive tensor algebra, which introduces an adaptive representation of tensors of many‐body operators via the local adjustment of the basis set quality. Given a series of locally supported fragment bases of a progressively lower quality, we formulate the explicit rules for tensor algebra operations dealing with adaptively resolved tensor operands. The formalism suggested is expected to enhance the applicability of certain local correlated many‐body methods of electronic structure theory, for example, those directly based on atomic orbitals (or any other localized basis functions in general). © 2014 Wiley Periodicals, Inc.

International Journal of Quantum Chemistry•Jul 7, 2014•4 citations

Representation of the Kato electron-electron cusp condition by wavelet-based density-operator expansions

Jānos Pipek, Szilvia Nagy

Since Kato proved his singularity condition for Coulomb potentials in 1957, there has been interest in the creation of wave functions that meet the prescriptions of the cusp conditions, necessary for high-precision quantum-mechanical calculations. It is well known, that wave-function expansions based on Slater determinants of one-electron functions are poorly convergent with respect to satisfying the electron-electron cusp condition. In this contribution we show that with the wavelet expansion of density operators even the local form of the electron-electron cusp condition is easily representable by Slater determinants of one-electron wavelet functions with a proper asymptotics of the expansion coefficients, which is explicitly calculated for Haar wavelets.

Physical Review A•Oct 12, 2001•10 citations

Generalized density functional theories using the k-electron densities: Development of kinetic energy functionals

Paul W. Ayers

Several explicit formulas for the kinetic energy of a many-electron system as a functional of the k-electron density are derived, with emphasis on the electron pair density. The emphasis is on general techniques for deriving approximate kinetic energy functionals and features generalized Weisacker bounds and methods using density-matrix reconstruction. Adapting results from statistical mechanics, a hierarchy of equations is derived that links electron pairs, triplets, quadruplets, etc.; this may be used to derive more accurate approximations. Several methods for defining the exact kinetic energy functional are presented, including the generalizations of the Levy and Lieb formulations of density-functional theory. Together with N-representability constraints on the k-density, this paper provides the basis for “generalized density functional theories” based on the electron pair density. There are also implications for conventional density-functional theory, notably regarding the development of more accurate density functionals for the kinetic energy.

Journal of Mathematical Physics•May 24, 2005•82 citations
Multiresolution (or wavelet) analysis offers a strictly local basis set for a systematic introduction of new details into Hilbert space operators. Using this tool we have previously developed an expansion method for density matrices. The set of density operators providing a given electron density plays an essential role in density functional theory, in the minimization of energy expectation values with the constraint that the electron density is fixed. In this contribution, using multiresolution analysis, we present an excellent quality density matrix expansion yielding a prescribed electron density, and compare it to other known methods. Due to the strictly local nature of the applied basis functions, our construction has the specific advantage that the resulting density matrix is correlated and N-representable in the infinite resolution limit. As a further consequence of this scheme we can conclude that the deviation of the exact kinetic energy functional from the Weizsäcker term is not a necessary consequence of the particle statistics.
Wavelet analysis, which is a shorthand notation for the concept of multiresolution analysis (MRA), becomes increasingly popular in high efficiency storage algorithms of complex spatial distributions. This approach is applied for describing wave functions of quantum systems. At any resolution level of MRA expansions a physical observable is represented by an infinite matrix which is “canonically” chosen as the projection of its operator in the Schrödinger picture onto the subspace of the given resolution. It is shown that this canonical choice is only a particular member of possible operator representations. Among these, there exits an optimal choice, usually different from the canonical one, which gives the best numerical values in eigenvalue problems. This construction works even in those cases, where the canonical definition is unusable. The commutation relation of physical operators is also studied in MRA subspaces. It is shown that the required commutation rules are satisfied in the fine resolution limit, whereas in coarse grained spaces a correction appears depending only on the representation of the momentum operator.
High Resolution Image Download MS PowerPoint Slide The density matrix renormalization group (DMRG) is a powerful numerical technique to solve strongly correlated quantum systems: it deals well with systems which are not dominated by a single configuration (unlike Coupled Cluster) and it converges rapidly to the Full Configuration Interaction (FCI) limit (unlike truncated Configuration Interaction (CI) expansions). In this work, we develop an algorithm integrating DMRG within the multiwavelet-based multiresolution analysis (MRA). Unlike fixed basis sets, multiwavelets offer an adaptive and hierarchical representation of functions, approaching the complete basis set limit to a specified precision. As a result, this combined technique leverages the multireference capability of DMRG and the complete basis set limit of MRA and multiwavelets. More specifically, we adopt a pre-existing Lagrangian optimization algorithm for orbitals represented in the MRA domain and improve its computational efficiency by replacing the original CI calculations with DMRG. Additionally, we substitute the reduced density matrices computation with the direct extraction of energy gradients from the DMRG tensors. We apply our method to small systems such H 2, He, HeH 2, BeH 2 and N 2 . The results demonstrate that our approach reduces the final energy while keeping the number of orbitals low compared to FCI calculations on an atomic orbital basis set.
We present a multiscale treatment of electron correlations based on hyperbolic wavelet expansions of Jastrow-type correlation functions. Wavelets provide hierarchical basis sets that can be locally adapted to the length- and energy-scales of physical phenomena. Combined with hyperbolic tensor products and local adaptive refinement near the interelectron cusp, these wavelet bases enable sparse representations of Jastrow factors. The computational efficiency of wavelets in electronic structure calculations is demonstrated within the coupled electron-pair approximation (local ansatz). Based on a diagrammatic multiresolution analysis, we discuss various kinds of sparsity features for matrix elements required by the local ansatz. Sparsity originates from the hierarchical structure and vanishing moments property of wavelet bases. This led us to a recurrence scheme for the evaluation of matrix elements with almost linear computational complexity with respect to the size of the underlying isotropic 3d-wavelet basis. Numerical studies for selected diagrams are presented for a homogeneous electron gas.
While the formalism of multiresolution analysis, based on wavelets and adaptive integral representations of operators, is actively progressing in electronic structure theory (mostly on the independent‐particle level and, recently, second‐order perturbation theory), the concepts of multiresolution and adaptivity can also be utilized within the traditional formulation of correlated (many‐particle) theory based on second quantization and the corresponding (generally nonorthogonal) tensor algebra. In this article, we present a formalism called scale‐adaptive tensor algebra, which introduces an adaptive representation of tensors of many‐body operators via the local adjustment of the basis set quality. Given a series of locally supported fragment bases of a progressively lower quality, we formulate the explicit rules for tensor algebra operations dealing with adaptively resolved tensor operands. The formalism suggested is expected to enhance the applicability of certain local correlated many‐body methods of electronic structure theory, for example, those directly based on atomic orbitals (or any other localized basis functions in general). © 2014 Wiley Periodicals, Inc.
Since Kato proved his singularity condition for Coulomb potentials in 1957, there has been interest in the creation of wave functions that meet the prescriptions of the cusp conditions, necessary for high-precision quantum-mechanical calculations. It is well known, that wave-function expansions based on Slater determinants of one-electron functions are poorly convergent with respect to satisfying the electron-electron cusp condition. In this contribution we show that with the wavelet expansion of density operators even the local form of the electron-electron cusp condition is easily representable by Slater determinants of one-electron wavelet functions with a proper asymptotics of the expansion coefficients, which is explicitly calculated for Haar wavelets.
Several explicit formulas for the kinetic energy of a many-electron system as a functional of the k-electron density are derived, with emphasis on the electron pair density. The emphasis is on general techniques for deriving approximate kinetic energy functionals and features generalized Weisacker bounds and methods using density-matrix reconstruction. Adapting results from statistical mechanics, a hierarchy of equations is derived that links electron pairs, triplets, quadruplets, etc.; this may be used to derive more accurate approximations. Several methods for defining the exact kinetic energy functional are presented, including the generalizations of the Levy and Lieb formulations of density-functional theory. Together with N-representability constraints on the k-density, this paper provides the basis for “generalized density functional theories” based on the electron pair density. There are also implications for conventional density-functional theory, notably regarding the development of more accurate density functionals for the kinetic energy.