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:2105.01088 · Physical Review Research

Quantum impurity models using superpositions of fermionic Gaussian states: Practical methods and applications

Samuel Boutin, Bela Bauer

Published Aug 24, 2021• 9 citations
arXivHTMLPDF

Abstract

The coherent superposition of nonorthogonal fermionic Gaussian states has been shown to be an efficient approximation to the ground states of quantum impurity problems [Bravyi and Gosset, Commun. Math. Phys. 356, 451 (2017)]. We propose a practical approach for performing a variational calculation based on such states. Our method is based on approximate imaginary-time equations of motion that decouple the dynamics of each Gaussian state forming the ansatz. It is independent of the lattice connectivity of the model and the implementation is highly parallelizable. To benchmark our variational method, we calculate the spin-spin correlation function and R\'enyi entanglement entropy of an Anderson impurity, allowing us to identify the screening cloud and compare to density matrix renormalization group calculations. Secondly, we study the screening cloud of the two-channel Kondo model, a problem difficult to tackle using existing numerical tools.

DOI
10.1103/physrevresearch.3.033188

Related papers

Efficient impurity-bath trial states from superposed Slater determinants

Izak Snyman, Serge Florens

The representation of ground states of fermionic quantum impurity problems as superpositions of Gaussian states has recently been given a rigorous mathematical foundation [S. Bravyi and D. Gosset, Commun. Math. Phys. 356, 451 (2017)]. It is natural to ask how many parameters are required for an efficient variational scheme based on this representation. An upper bound is $O({N}^{2})$, where $N$ is the system size, which corresponds to the number parameters needed to specify an arbitrary Gaussian state. We provide an alternative representation, with more favorable scaling, only requiring $O(N)$ parameters, that we illustrate for the interacting resonant-level model. We achieve the reduction by associating mean-field-like parent Hamiltonians with the individual terms in the superposition, using physical insight to retain only the most relevant channels in each parent Hamiltonian. We benchmark our variational ansatz against the numerical renormalization group, and compare our results to other variational schemes of a similar nature to ours. Apart from the ground-state energy, we also study the spectrum of the covariance matrix---a very stringent measure of accuracy. Our approach outperforms some existing variational schemes and remains quantitatively accurate in the numerically challenging near-critical regime.

arXiv:2106.14027•Physical review. B./Physical review. B•Nov 19, 2021•5 citations

Variational principle for quantum impurity systems in and out of equilibrium: Application to Kondo problems

Yuto Ashida, Tao Shi, Mari Carmen Bañuls +2 more

We provide a detailed formulation of the recently proposed variational approach [Y. Ashida, T. Shi, M.-C. Ba\~nuls, J. I. Cirac, and E. Demler, Phys. Rev. Lett. 121, 026805 (2018)] to study ground-state properties and out-of-equilibrium dynamics for generic quantum spin-impurity systems. Motivated by the original ideas of Tomonaga, Lee, Low, and Pines, we construct a canonical transformation that completely decouples the impurity from the bath degrees of freedom. By combining this transformation with a Gaussian ansatz for the fermionic bath, we obtain a family of variational many-body states that can efficiently encode the strong entanglement between the impurity and fermions of the bath. We give a detailed derivation of equations of motions in the imaginary- and real-time evolutions on the variational manifold. We benchmark our approach by applying it to investigate ground-state and dynamical properties of the anisotropic Kondo model and compare results with those obtained using the matrix-product state (MPS) ansatz. We show that our approach can achieve an accuracy comparable to MPS-based methods with several orders of magnitude fewer variational parameters than the corresponding MPS ansatz. Comparisons to the Yosida ansatz and the exact solution from the Bethe ansatz are also discussed. We use our approach to investigate the two-lead Kondo model and analyze its long-time spatiotemporal behavior and the conductance behavior at finite bias and magnetic fields. The obtained results are consistent with the previous findings in the Anderson model and the exact solutions at the Toulouse point.

arXiv:1802.03861•Physical review. B./Physical review. B•Jul 11, 2018•37 citations

Disentangling interacting systems with fermionic Gaussian circuits: Application to quantum impurity models

Ang-Kun Wu, Benedikt Kloss, Wladislaw Krinitsin +3 more

Tensor network quantum states are powerful tools for strongly correlated systems, tailored to capture local correlations such as in ground states with entanglement area laws. When applying tensor network states to interacting fermionic systems, a proper choice of the basis or orbitals can reduce the bond dimension of tensors and provide physically relevant orbitals. We introduce such a change of basis with unitary gates obtained from compressing fermionic Gaussian states into quantum circuits corresponding to various tensor networks. These circuits can reduce the ground-state entanglement entropy and improve the performance of algorithms such as the density matrix renormalization group. We study the Anderson impurity model with one and two impurities to show the potential of the method for improving computational efficiency and interpreting impurity physics. Furthermore, fermionic Gaussian circuits can also suppress entanglement during the time evolution out of low-energy state. Last, we consider Gaussian multiscale entanglement renormalization ansatz (GMERA) circuits which compress fermionic Gaussian states hierarchically. The emergent coarse-grained physical models from these GMERA circuits are studied in terms of their entanglement properties and suitability for performing time evolution.

arXiv:2212.09798•Physical review. B./Physical review. B•Jan 9, 2025•11 citations

Block Lanczos density-matrix renormalization group method for general Anderson impurity models: Application to magnetic impurity problems in graphene

Tomonori Shirakawa, Seiji Yunoki

We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multiorbital Anderson impurity model where not only inter- and intraorbital Coulomb interactions but also Hund's coupling and pair hopping interactions are included. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene with a structural defect and with a single hydrogen or fluorine absorbed, where a single Anderson impurity is coupled to conduction electrons in the honeycomb lattice. These models include (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis of the local dynamical magnetic susceptibility and the local density of states at the impurity site reveals that, for the particle-hole symmetric case at half-filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening, while the ground state of the other two models forms a spin-singlet state where the impurity moment is screened by the conduction electrons. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence or absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}3}$, differently from the noninteracting limit ($\ensuremath{\propto}$ ${r}^{\ensuremath{-}2}$), for model (i) and as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}4}$, exactly the same as the noninteracting limit, for models (ii) and (iii) in the asymptotic $r$, where $r$ is the distance between the impurity site and the conduction site. Finally, based on our results, we shed light on recent experiments on graphene where the formation of local magnetic moments as well as the Kondo-like behavior have been observed.

arXiv:1409.7995•Physical Review B•Nov 5, 2014•33 citations

Solving Quantum Impurity Problems in and out of Equilibrium with the Variational Approach

Yuto Ashida, Tao Shi, Mari Carmen Bañuls +2 more

A versatile and efficient variational approach is developed to solve in- and out-of-equilibrium problems of generic quantum spin-impurity systems. Employing the discrete symmetry hidden in spin-impurity models, we present a new canonical transformation that completely decouples the impurity and bath degrees of freedom. Combining it with Gaussian states, we present a family of many-body states to efficiently encode nontrivial impurity-bath correlations. We demonstrate its successful application to the anisotropic and two-lead Kondo models by studying their spatiotemporal dynamics and universal behavior in the correlations, relaxation times, and the differential conductance. We compare them to previous analytical and numerical results. In particular, we apply our method to study new types of nonequilibrium phenomena that have not been studied by other methods, such as long-time crossover in the ferromagnetic easy-plane Kondo model. The present approach will be applicable to a variety of unsolved problems in solid-state and ultracold-atomic systems.

arXiv:1801.05825•Physical Review Letters•Jul 11, 2018•53 citations

Variational study of fermionic and bosonic systems with non-Gaussian states: Theory and applications

Shi, T., Eugene Demler, J. I. Cirac

We present a new variational method for investigating the ground state and out of equilibrium dynamics of quantum many-body bosonic and fermionic systems. Our approach is based on constructing variational wavefunctions which extend Gaussian states by including generalized canonical transformations between the fields. The key advantage of such states compared to simple Gaussian states is presence of non-factorizable correlations and the possibility of describing states with strong entanglement between particles. In contrast to the commonly used canonical transformations, such as the polaron or Lang-Firsov transformations, we allow parameters of the transformations to be time dependent, which extends their regions of applicability. We derive equations of motion for the parameters characterizing the states both in real and imaginary time using the differential structure of the variational manifold. The ground state can be found by following the imaginary time evolution until it converges to a steady state. Collective excitations in the system can be obtained by linearizing the real-time equations of motion in the vicinity of the imaginary time steady-state solution. Our formalism allows us not only to determine the energy spectrum of quasiparticles and their lifetime, but to obtain the complete spectral functions and to explore far out of equilibrium dynamics such as coherent evolution following a quantum quench. We illustrate and benchmark this framework with several examples: a single polaron in the Holstein and Su-Schrieer-Heeger models, non-equilibrium dynamics in the spin-boson and Kondo models, the superconducting to charge density wave phase transitions in the Holstein model.

•Mar 1, 2018•106 citations
The coherent superposition of nonorthogonal fermionic Gaussian states has been shown to be an efficient approximation to the ground states of quantum impurity problems [Bravyi and Gosset, Commun. Math. Phys. 356, 451 (2017)]. We propose a practical approach for performing a variational calculation based on such states. Our method is based on approximate imaginary-time equations of motion that decouple the dynamics of each Gaussian state forming the ansatz. It is independent of the lattice connectivity of the model and the implementation is highly parallelizable. To benchmark our variational method, we calculate the spin-spin correlation function and R'enyi entanglement entropy of an Anderson impurity, allowing us to identify the screening cloud and compare to density matrix renormalization group calculations. Secondly, we study the screening cloud of the two-channel Kondo model, a problem difficult to tackle using existing numerical tools.
The representation of ground states of fermionic quantum impurity problems as superpositions of Gaussian states has recently been given a rigorous mathematical foundation [S. Bravyi and D. Gosset, Commun. Math. Phys. 356, 451 (2017)]. It is natural to ask how many parameters are required for an efficient variational scheme based on this representation. An upper bound is O(N2)O({N}^{2})O(N2), where NNN is the system size, which corresponds to the number parameters needed to specify an arbitrary Gaussian state. We provide an alternative representation, with more favorable scaling, only requiring O(N)O(N)O(N) parameters, that we illustrate for the interacting resonant-level model. We achieve the reduction by associating mean-field-like parent Hamiltonians with the individual terms in the superposition, using physical insight to retain only the most relevant channels in each parent Hamiltonian. We benchmark our variational ansatz against the numerical renormalization group, and compare our results to other variational schemes of a similar nature to ours. Apart from the ground-state energy, we also study the spectrum of the covariance matrix---a very stringent measure of accuracy. Our approach outperforms some existing variational schemes and remains quantitatively accurate in the numerically challenging near-critical regime.
We provide a detailed formulation of the recently proposed variational approach [Y. Ashida, T. Shi, M.-C. Ba\~nuls, J. I. Cirac, and E. Demler, Phys. Rev. Lett. 121, 026805 (2018)] to study ground-state properties and out-of-equilibrium dynamics for generic quantum spin-impurity systems. Motivated by the original ideas of Tomonaga, Lee, Low, and Pines, we construct a canonical transformation that completely decouples the impurity from the bath degrees of freedom. By combining this transformation with a Gaussian ansatz for the fermionic bath, we obtain a family of variational many-body states that can efficiently encode the strong entanglement between the impurity and fermions of the bath. We give a detailed derivation of equations of motions in the imaginary- and real-time evolutions on the variational manifold. We benchmark our approach by applying it to investigate ground-state and dynamical properties of the anisotropic Kondo model and compare results with those obtained using the matrix-product state (MPS) ansatz. We show that our approach can achieve an accuracy comparable to MPS-based methods with several orders of magnitude fewer variational parameters than the corresponding MPS ansatz. Comparisons to the Yosida ansatz and the exact solution from the Bethe ansatz are also discussed. We use our approach to investigate the two-lead Kondo model and analyze its long-time spatiotemporal behavior and the conductance behavior at finite bias and magnetic fields. The obtained results are consistent with the previous findings in the Anderson model and the exact solutions at the Toulouse point.
Tensor network quantum states are powerful tools for strongly correlated systems, tailored to capture local correlations such as in ground states with entanglement area laws. When applying tensor network states to interacting fermionic systems, a proper choice of the basis or orbitals can reduce the bond dimension of tensors and provide physically relevant orbitals. We introduce such a change of basis with unitary gates obtained from compressing fermionic Gaussian states into quantum circuits corresponding to various tensor networks. These circuits can reduce the ground-state entanglement entropy and improve the performance of algorithms such as the density matrix renormalization group. We study the Anderson impurity model with one and two impurities to show the potential of the method for improving computational efficiency and interpreting impurity physics. Furthermore, fermionic Gaussian circuits can also suppress entanglement during the time evolution out of low-energy state. Last, we consider Gaussian multiscale entanglement renormalization ansatz (GMERA) circuits which compress fermionic Gaussian states hierarchically. The emergent coarse-grained physical models from these GMERA circuits are studied in terms of their entanglement properties and suitability for performing time evolution.
We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multiorbital Anderson impurity model where not only inter- and intraorbital Coulomb interactions but also Hund's coupling and pair hopping interactions are included. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene with a structural defect and with a single hydrogen or fluorine absorbed, where a single Anderson impurity is coupled to conduction electrons in the honeycomb lattice. These models include (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis of the local dynamical magnetic susceptibility and the local density of states at the impurity site reveals that, for the particle-hole symmetric case at half-filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening, while the ground state of the other two models forms a spin-singlet state where the impurity moment is screened by the conduction electrons. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence or absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}3}$, differently from the noninteracting limit ($\ensuremath{\propto}$ ${r}^{\ensuremath{-}2}$), for model (i) and as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}4}$, exactly the same as the noninteracting limit, for models (ii) and (iii) in the asymptotic rrr, where rrr is the distance between the impurity site and the conduction site. Finally, based on our results, we shed light on recent experiments on graphene where the formation of local magnetic moments as well as the Kondo-like behavior have been observed.
A versatile and efficient variational approach is developed to solve in- and out-of-equilibrium problems of generic quantum spin-impurity systems. Employing the discrete symmetry hidden in spin-impurity models, we present a new canonical transformation that completely decouples the impurity and bath degrees of freedom. Combining it with Gaussian states, we present a family of many-body states to efficiently encode nontrivial impurity-bath correlations. We demonstrate its successful application to the anisotropic and two-lead Kondo models by studying their spatiotemporal dynamics and universal behavior in the correlations, relaxation times, and the differential conductance. We compare them to previous analytical and numerical results. In particular, we apply our method to study new types of nonequilibrium phenomena that have not been studied by other methods, such as long-time crossover in the ferromagnetic easy-plane Kondo model. The present approach will be applicable to a variety of unsolved problems in solid-state and ultracold-atomic systems.
We present a new variational method for investigating the ground state and out of equilibrium dynamics of quantum many-body bosonic and fermionic systems. Our approach is based on constructing variational wavefunctions which extend Gaussian states by including generalized canonical transformations between the fields. The key advantage of such states compared to simple Gaussian states is presence of non-factorizable correlations and the possibility of describing states with strong entanglement between particles. In contrast to the commonly used canonical transformations, such as the polaron or Lang-Firsov transformations, we allow parameters of the transformations to be time dependent, which extends their regions of applicability. We derive equations of motion for the parameters characterizing the states both in real and imaginary time using the differential structure of the variational manifold. The ground state can be found by following the imaginary time evolution until it converges to a steady state. Collective excitations in the system can be obtained by linearizing the real-time equations of motion in the vicinity of the imaginary time steady-state solution. Our formalism allows us not only to determine the energy spectrum of quasiparticles and their lifetime, but to obtain the complete spectral functions and to explore far out of equilibrium dynamics such as coherent evolution following a quantum quench. We illustrate and benchmark this framework with several examples: a single polaron in the Holstein and Su-Schrieer-Heeger models, non-equilibrium dynamics in the spin-boson and Kondo models, the superconducting to charge density wave phase transitions in the Holstein model.