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:2610.07512

X2C-inspired approximation of the Dirac coupling operator in a multiwavelet basis

Quentin Pitteloud, Jacopo Masotti, Luca Frediani

quant-phphysics.chem-phSubmitted Oct 5, 2026• Announced Oct 7, 2026
arXivHTMLPDF

Abstract

We present a method to obtain the proven exact 2-component (X2C) method, normally applied to Gaussian-type fixed Hilbert bases, in an adaptative non-uniform grid basis of the Hilbert space, namely Multiwavelets. The X2C method manages to block-diagonalise the Dirac Hamiltonian restricted to the subspace of the basis, in an energy-independent manner. The method takes inspiration from the success of the atomic mean-field X2C method on the Gaussian bases to represent the coupling operator as a projector onto the small components of the eigenstates of the constituent atoms of the system.

Comments
REHE2026

Related papers

Numerical method for solving Schrödinger equations using wavelet basis functions

Myoung Won Cho

Abstract This paper presents a numerical method for solving Schrödinger equations with arbitrary potentials with the adoption of the qunatum diagonalization scheme and basis functions in the discrete wavelet transform. The quantum diagonalization scheme, developed to solve the Hubbard model and so on, is known to allow the diagonalization of a Hamiltonian matrix of a much large size because it does not require to maintain the entire matrix in memory. And, the adoption of wavelet basis functions enables representing eigenstates efficiently using a number of bases much smaller than the number of grid points to integral and reducing the computation time to perform the diagonalization process. With the help of another additional techniques, the method allows to find solutions of a Schrödinger equation much exactly and efficiently. The validity of the method and its efficiency depending on several conditions are proven and surveyed by applying it to several problems for which the exact solutions are known.

•Jun 30, 2025•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

Quantum chemistry in Fock space. II. Effective Hamiltonians in Fock space

Werner Kutzelnigg, Sigurd Koch

The concept of an effective Hamiltonian in Fock space is introduced. It is based on the division of the entire one-particle space into subspaces of ‘‘active’’ and ‘‘inactive’’ orbitals. The effective Fock space Hamiltonian has—for active model states—the same eigenvalues as the full Hamiltonian. The theory outlined in this context differs from that of paper I mainly in a different definition of the ‘‘diagonal part’’ of an operator, and in the fact that the ‘‘quasidegenerate case’’ applies throughout. The separation theorem, and as a consequence the connected diagram theorem, is shown to hold, in a more limited sense though, even for those normalizations where it did not in the context of universal wave and energy operators. Unlike in the theory of the ‘‘universal’’ operators of paper I the Fock space and n-particle Hilbert space approaches with analogous normalizations are no longer equivalent. In particular, the Primas normalization with a fully Lie-algebraic structure does not lead to a connected diagram expansion if it is formulated in n-particle Hilbert space, only so in a Fock space formulation. In n-particle Hilbert space with the present definition of the diagonal part of an operator the normalizations b (‘‘canonical’’) and c (‘‘Primas’’) happen to agree. As an alternative to the construction of the wave and energy operator W and L by perturbation theory the nonperturbative approach is presented as a generalization of the coupled-cluster method, in detail both in the intermediate and in the unitary normalization. In the unitary variant only a linear system for σ (the logarithm of the wave operator) has to be solved in order to get L correct through fifth order in perturbation theory with important contributions of higher orders included. A generalization of the Hartree–Fock method to Fock space theory is outlined, which guarantees stationarity of all (active) eigenstates with respect to one-particle transformations. A generalized electron pair theory is also defined. An analysis of the necessary computational steps shows that the nonperturbative approaches do not require significantly more computational effort than perturbation theory to the corresponding order. As a numerical example the H2 molecule in a small basis is discussed.

The Journal of Chemical Physics•Nov 1, 1983•171 citations

Quantum system partitioning at the single-particle level

Adrian H. Mühlbach, Markus Reiher

We discuss the partitioning of a quantum system through subsystem separation by unitary block-diagonalization (SSUB) applied to a Fock operator. For a one-particle Hilbert space, this separation can be formulated in a very general way. Therefore, it can be applied to very different partitionings ranging from those driven by features in the molecular structure (such as a solute surrounded by solvent molecules or an active site in an enzyme) to those that aim at an orbital separation (such as core-valence separation). Our framework embraces recent developments of Manby and Miller as well as the older ones of Huzinaga and Cantu. Projector-based embedding is simplified and accelerated by SSUB. Moreover, it directly relates to decoupling approaches for relativistic four-component many-electron theory. For a Fock operator based on the Dirac one-electron Hamiltonian, one would like to separate the so-called positronic (negative-energy) states from the electronic bound and continuum states. The exact two-component (X2C) approach developed for this purpose becomes a special case of the general SSUB framework and may therefore be viewed as a system-environment decoupling approach. Moreover, for SSUB, there exists no restriction with respect to the number of subsystems that are generated-in the limit, decoupling of all single-particle states is recovered, which represents exact diagonalization of the problem. The fact that a Fock operator depends on its eigenvectors poses challenges to all system-environment decoupling approaches and is discussed in terms of the SSUB framework. Apart from improved conceptual understanding, these relations bring about technical advances as developments in different fields can immediately cross-fertilize one another. As an important example, we discuss the atomic decomposition of the unitary block-diagonalization matrix in X2C-type approaches that can inspire approaches for the efficient partitioning of large total systems based on SSUB.

arXiv:1809.04476•The Journal of Chemical Physics•Nov 12, 2018•19 citations

Two-center Dirac equation in a Gaussian basis set

Dávid Ferenc

A method for solving the two-center Dirac equation using a Gaussian-product basis set tailored to the axial symmetry of diatomic systems is presented. For the ground state of H$_2^+$ at an internuclear distance $R=2.0~a_0$, we obtain a relativistic energy with relative precision at the $10^{-18}$ level. The properties of the basis set allow for an analytic Fourier transform of the wave function. As a first application we evaluate, for the first time, the zero-potential contribution to the one-loop self-energy for H$_2^+$ in Feynman and Coulomb gauges. Although this contribution is not separately gauge invariant, its calculation provides a proof of principle for the use of the present framework in prospective QED calculations.

quant-pharXiv:2610.03217•Oct 5, 2026

Two-dimensional quantum propagation using wavelets in space and time

Douglas K. Sparks, Bruce R. Johnson

A recent method for solving the time-dependent Schrodinger equation has been developed using expansions in compact-support wavelet bases in both space and time [H. Wang et al., J. Chem. Phys. 121, 7647 (2004)]. This method represents an exact quantum mixed time-frequency approach, with special initial temporal wavelets used to solve the initial value problem. The present work is a first extension of the method to multiple spatial dimensions applied to a simple two-dimensional (2D) coupled anharmonic oscillator problem. A wavelet-discretized version of norm preservation for time-independent Hamiltonians discovered in the earlier one-dimensional investigation is verified to hold as well in 2D and, by implication, in higher numbers of spatial dimensions. The wavelet bases are not restricted to rectangular domains, a fact which is exploited here in a 2D adaptive version of the algorithm.

The Journal of Chemical Physics•Sep 20, 2006•8 citations
We present a method to obtain the proven exact 2-component (X2C) method, normally applied to Gaussian-type fixed Hilbert bases, in an adaptative non-uniform grid basis of the Hilbert space, namely Multiwavelets. The X2C method manages to block-diagonalise the Dirac Hamiltonian restricted to the subspace of the basis, in an energy-independent manner. The method takes inspiration from the success of the atomic mean-field X2C method on the Gaussian bases to represent the coupling operator as a projector onto the small components of the eigenstates of the constituent atoms of the system.
Abstract This paper presents a numerical method for solving Schrödinger equations with arbitrary potentials with the adoption of the qunatum diagonalization scheme and basis functions in the discrete wavelet transform. The quantum diagonalization scheme, developed to solve the Hubbard model and so on, is known to allow the diagonalization of a Hamiltonian matrix of a much large size because it does not require to maintain the entire matrix in memory. And, the adoption of wavelet basis functions enables representing eigenstates efficiently using a number of bases much smaller than the number of grid points to integral and reducing the computation time to perform the diagonalization process. With the help of another additional techniques, the method allows to find solutions of a Schrödinger equation much exactly and efficiently. The validity of the method and its efficiency depending on several conditions are proven and surveyed by applying it to several problems for which the exact solutions are known.
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.
The concept of an effective Hamiltonian in Fock space is introduced. It is based on the division of the entire one-particle space into subspaces of ‘‘active’’ and ‘‘inactive’’ orbitals. The effective Fock space Hamiltonian has—for active model states—the same eigenvalues as the full Hamiltonian. The theory outlined in this context differs from that of paper I mainly in a different definition of the ‘‘diagonal part’’ of an operator, and in the fact that the ‘‘quasidegenerate case’’ applies throughout. The separation theorem, and as a consequence the connected diagram theorem, is shown to hold, in a more limited sense though, even for those normalizations where it did not in the context of universal wave and energy operators. Unlike in the theory of the ‘‘universal’’ operators of paper I the Fock space and n-particle Hilbert space approaches with analogous normalizations are no longer equivalent. In particular, the Primas normalization with a fully Lie-algebraic structure does not lead to a connected diagram expansion if it is formulated in n-particle Hilbert space, only so in a Fock space formulation. In n-particle Hilbert space with the present definition of the diagonal part of an operator the normalizations b (‘‘canonical’’) and c (‘‘Primas’’) happen to agree. As an alternative to the construction of the wave and energy operator W and L by perturbation theory the nonperturbative approach is presented as a generalization of the coupled-cluster method, in detail both in the intermediate and in the unitary normalization. In the unitary variant only a linear system for σ (the logarithm of the wave operator) has to be solved in order to get L correct through fifth order in perturbation theory with important contributions of higher orders included. A generalization of the Hartree–Fock method to Fock space theory is outlined, which guarantees stationarity of all (active) eigenstates with respect to one-particle transformations. A generalized electron pair theory is also defined. An analysis of the necessary computational steps shows that the nonperturbative approaches do not require significantly more computational effort than perturbation theory to the corresponding order. As a numerical example the H2 molecule in a small basis is discussed.
We discuss the partitioning of a quantum system through subsystem separation by unitary block-diagonalization (SSUB) applied to a Fock operator. For a one-particle Hilbert space, this separation can be formulated in a very general way. Therefore, it can be applied to very different partitionings ranging from those driven by features in the molecular structure (such as a solute surrounded by solvent molecules or an active site in an enzyme) to those that aim at an orbital separation (such as core-valence separation). Our framework embraces recent developments of Manby and Miller as well as the older ones of Huzinaga and Cantu. Projector-based embedding is simplified and accelerated by SSUB. Moreover, it directly relates to decoupling approaches for relativistic four-component many-electron theory. For a Fock operator based on the Dirac one-electron Hamiltonian, one would like to separate the so-called positronic (negative-energy) states from the electronic bound and continuum states. The exact two-component (X2C) approach developed for this purpose becomes a special case of the general SSUB framework and may therefore be viewed as a system-environment decoupling approach. Moreover, for SSUB, there exists no restriction with respect to the number of subsystems that are generated-in the limit, decoupling of all single-particle states is recovered, which represents exact diagonalization of the problem. The fact that a Fock operator depends on its eigenvectors poses challenges to all system-environment decoupling approaches and is discussed in terms of the SSUB framework. Apart from improved conceptual understanding, these relations bring about technical advances as developments in different fields can immediately cross-fertilize one another. As an important example, we discuss the atomic decomposition of the unitary block-diagonalization matrix in X2C-type approaches that can inspire approaches for the efficient partitioning of large total systems based on SSUB.
A method for solving the two-center Dirac equation using a Gaussian-product basis set tailored to the axial symmetry of diatomic systems is presented. For the ground state of H2+_2^+2+​ at an internuclear distance R=2.0 a0R=2.0~a_0R=2.0 a0​, we obtain a relativistic energy with relative precision at the 10−1810^{-18}10−18 level. The properties of the basis set allow for an analytic Fourier transform of the wave function. As a first application we evaluate, for the first time, the zero-potential contribution to the one-loop self-energy for H2+_2^+2+​ in Feynman and Coulomb gauges. Although this contribution is not separately gauge invariant, its calculation provides a proof of principle for the use of the present framework in prospective QED calculations.
A recent method for solving the time-dependent Schrodinger equation has been developed using expansions in compact-support wavelet bases in both space and time [H. Wang et al., J. Chem. Phys. 121, 7647 (2004)]. This method represents an exact quantum mixed time-frequency approach, with special initial temporal wavelets used to solve the initial value problem. The present work is a first extension of the method to multiple spatial dimensions applied to a simple two-dimensional (2D) coupled anharmonic oscillator problem. A wavelet-discretized version of norm preservation for time-independent Hamiltonians discovered in the earlier one-dimensional investigation is verified to hold as well in 2D and, by implication, in higher numbers of spatial dimensions. The wavelet bases are not restricted to rectangular domains, a fact which is exploited here in a 2D adaptive version of the algorithm.