Matthieu Saubanère

The one-particle reduced density-matrix (1-RDM) functional theory is a promising alternative to density-functional theory (DFT) that uses the 1-RDM rather than the electronic density as a basic variable. However, long-standing challenges such as the lack of the Kohn–Sham scheme and the complexity of the pure N-representability conditions are still impeding its wild utilization. Fortunately, ensemble N-representability conditions derived in the natural orbital basis are known and trivial such that almost every functional of the 1-RDM is actually natural orbital functional, which does not perform well for all the correlation regimes. In this work, we propose a variational minimization scheme in the ensemble N-representable domain that is not restricted to the natural orbital representation of the 1-RDM. We show that splitting the minimization into the diagonal and off-diagonal parts of the 1-RDM can open the way toward the development of functionals of the orbital occupations, which remains a challenge for the generalization of site-occupation functional theory in chemistry. Our approach is tested on the uniform Hubbard model using the Müller and the Töws–Pastor functionals, as well as on the dihydrogen molecule using the Müller functional.

Mechanical degradation in electrode materials during successive electrochemical cycling is critical for battery lifetime and aging properties. A common strategy to mitigate electrode mechanical degradation is to suppress the volume variation induced by Li/Na intercalation/deintercalation, thereby designing strain-less electrodes. In this study, we investigate the electrochemically-induced volume variation in layered and spinel compounds used in Li-ion and Na-ion battery electrode materials through density functional theory computations. Specifically, we propose to decompose the volume variation into electronic, ionic, and structural contributions. Based on this analysis, we suggest methods to separately influence each contribution through strategies such as chemical substitution, doping, and polymorphism. Altogether, we conclude that volume variations can be controlled by designing either mechanically hard or compact electrode materials.

The performance of embedding methods is directly tied to the quality of the bath orbital construction. In this paper, we develop a versatile framework, enabling the investigation of the optimal construction of the orbitals of the bath. As of today, in state-of-the-art embedding methods, the orbitals of the bath are constructed by performing a Singular Value Decomposition (SVD) on the impurity-environment part of the one-body reduced density matrix, as originally presented in Density Matrix Embedding Theory. Recently, the equivalence between the SVD protocol and the use of unitary transformation, the so-called Block-Householder transformation, has been established. We present a generalization of the Block-Householder transformation by introducing additional flexible parameters. The additional parameters are optimized such that the bath-orbitals fulfill physically motivated constraints. The efficiency of the approach is discussed and exemplified in the context of the half-filled Hubbard model in one-dimension.

We derive recursive relations for the Schrieffer-Wolff (SW) transformation applied to the half-filled Hubbard dimer. While the standard SW transformation is set to block-diagonalize the transformed Hamiltonian solely at the first order of perturbation, we infer from recursive relations two types of modifications, variational or iterative, that approximate, or even enforce for the homogeneous case, the desired block-diagonalization at infinite order of perturbation. The modified SW unitary transformations are then used to design an test quantum algorithms adapted to the noisy and fault-tolerant era. This work paves the way toward the design of alternative quantum algorithms for the general Hubbard Hamiltonian.

Quantum Chemistry and Physics have been pinpointed as killer applications for quantum computers, and quantum algorithms have been designed to solve the Schrödinger equation with the wavefunction formalism. It is yet limited to small systems, as their size is dictated by the number of qubits available. Computations on large systems rely mainly on mean-field-type approaches such as density functional theory, for which no quantum advantage has been envisioned so far. In this work, we question this a priori by proposing a counter-intuitive mapping from the non-interacting to an auxiliary interacting Hamiltonian that may provide the desired advantage.

We derive recursive relations for the Schrieffer--Wolff (SW) transformation applied to the half-filled Hubbard dimer. While the standard SW transformation is set to block-diagonalize the transformed Hamiltonian solely at the first order of perturbation, we infer from recursive relations two types of modifications, variational or iterative, that approach, or even enforce for the homogeneous case, the desired block-diagonalization at infinite order of perturbation. The modified SW unitary transformations are then used to design an test quantum algorithms adapted to the noisy and fault-tolerant era. This work paves the way toward the design of alternative quantum algorithms for the general Hubbard Hamiltonian.

Recently, some of the authors introduced the use of the Householder transformation as a simple and intuitive method for embedding local molecular fragments [see Sekaran et al., Phys. Rev. B 104, 035121 (2021) and Sekaran et al., Computation 10, 45 (2022)]. In this work, we present an extension of this approach to the more general case of multi-orbital fragments using the block version of the Householder transformation applied to the one-body reduced density matrix, unlocking the applicability to general quantum chemistry/condensed matter physics Hamiltonians. A step-by-step construction of the block Householder transformation is presented. Both physical and numerical areas of interest of the approach are highlighted. The specific mean-field (noninteracting) case is thoroughly detailed as it is shown that the embedding of a given N spin–orbital fragment leads to the generation of two separated sub-systems: (1) a 2N spin–orbitals “fragment+bath” cluster that exactly contains N electrons and (2) a remaining cluster’s “environment” described by so-called core electrons. We illustrate the use of this transformation in different cases of embedding scheme for practical applications. We particularly focus on the extension of the previously introduced Local Potential Functional Embedding Theory and Householder-transformed Density Matrix Functional Embedding Theory to the case of multi-orbital fragments. These calculations are realized on different types of systems, such as model Hamiltonians (Hubbard rings) and ab initio molecular systems (hydrogen rings).