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.
Simulating dissipative quantum many-body systems presents a significant theoretical challenge. Here, the authors develop an efficient nonperturbative framework that is tailored to dissipative impurity problems, combining a superposition of Gaussian states (SGS) variational ansatz with the quantum trajectory approach. They showcase the method with a model of a spinful impurity subject to two-body losses embedded in a fermionic bath, for which they explore the full crossover from weak to strong dissipation. Their approach captures the emergence of a loss-induced Kondo effect under strong dissipation, and reveals an exotic negative conductance phenomenon at intermediate loss rates.
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), 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.