We present a method for the modeling of fermionic reservoirs using a new class of ancillary damped fermions, dubbed purified pseudofermions, which exhibit unusual free correlations. We show that this key feature, when combined with existing efficient decomposition algorithms for the reservoir correlation functions, enables the development of an easily implementable and accurate scheme for constructing effective models of fermionic reservoirs. We numerically demonstrate the validity, accuracy, efficiency and potential use of our method by studying the particle transport of spinless fermions in a one-dimensional chain. Beyond its utility as a quantum impurity solver, our method holds promise for addressing a wide range of problems involving extended systems in fields like quantum transport, quantum thermodynamics, thermal engines and nonequilibrium phase transitions.
We present a method for the modeling of fermionic reservoirs using a new class of ancillary damped fermions, dubbed purified pseudofermions, which exhibit unusual free correlations. We show that this key feature, when combined with existing efficient decomposition algorithms for the reservoir correlation functions, enables the development of an easily implementable and accurate scheme for constructing effective models of fermionic reservoirs. We numerically demonstrate the validity, accuracy, efficiency and potential use of our method by studying the particle transport of spinless fermions in a one-dimensional chain. Beyond its utility as a quantum impurity solver, our method holds promise for addressing a wide range of problems involving extended systems in fields like quantum transport, quantum thermodynamics, thermal engines and nonequilibrium phase transitions.
We develop a discrete fermion approach for modeling the strong interaction of an arbitrary system interacting with continuum electronic reservoirs. The approach is based on a pseudofermion decomposition of the continuum bath correlation functions and is only limited by the accuracy of this decomposition. We show that to obtain this decomposition, one can allow for imaginary pseudofermion parameters, and strong damping in individual pseudofermions, without introducing unwanted approximations. For a noninteracting single-resonant level, we benchmark our approach against an analytical solution and an exact hierarchical-equations-of-motion approach. We also show that, for the interacting case, this simple method can capture the strongly correlated low-temperature physics of Kondo resonance, even in the difficult scaling limit, by employing matrix product state techniques.
The hierarchical equations of motion (HEOM) approach is an accurate method to simulate open system quantum dynamics, which allows for systematic convergence to numerically exact results. To represent effects of the bath, the reservoir correlation functions are usually decomposed into summation of multiple exponential terms in the HEOM method. Since the reservoir correlation functions become highly non-Markovian at low temperatures or when the bath has complex band structures, a present challenge is to obtain accurate exponential decompositions that allow efficient simulation with the HEOM. In this work, we employ the barycentric representation to approximate the Fermi function and hybridization functions in the frequency domain. This method, by approximating these functions with optimized rational decomposition, greatly reduces the number of basis functions in decomposing the reservoir correlation functions, which further allows the HEOM method to be applied to ultralow temperature and general band structures. We demonstrate the efficiency, accuracy, and long-time stability of this decomposition scheme by applying it to the Anderson impurity model in the low-temperature regime with the Lorentzian and tight-binding hybridization functions.
We discuss a new method to perform numerical simulations of one-dimensional systems with fermion and boson degrees of freedom. The method is based on a direct-space, imaginary-time representation of the fermion field. It is fast so that systems having up to 100 sites can easily be simulated. In addition, the method provides an intuitive physical "picture" of the ground state of a one-dimensional many-body system. We discuss in detail how to implement the method and how to compute various physical quantities. In particular, we show how to extend the method to study averages of off-diagonal quantities in an occupation-number representation. To assess the accuracy of our procedure, we apply it to free fermions in one dimension and compare with exact results. We then study a model of spinless interacting fermions and obtain the expected phase structure and behavior of correlation functions. We also consider the extended Hubbard model at various points in its phase diagram and study the behavior of spin-density, charge-density, and pairing correlation functions. We then study the Gross-Neveu model and show how the behavior depends on the number of fermion flavors. Finally, we consider an electron-phonon model and study its behavior both in the one-particle polaron sector and in the half-filled-band case. Along the way we show pictures of the ground-state configurations that give physical insight into the properties of the systems, like charge-density-wave, spin-density-wave, and superconducting states, "fractional charges," and solitons. We conclude by comparing our method with other methods and discuss the possibility of extending it to higher dimensions.
We present a methodology to simulate the quantum thermodynamics of thermal machines which are built from an interacting working medium in contact with fermionic reservoirs at a fixed temperature and chemical potential.Our method works at a finite temperature, beyond linear response and weak systemreservoir coupling, and allows for nonquadratic interactions in the working medium.The method uses mesoscopic reservoirs, continuously damped toward thermal equilibrium, in order to represent continuum baths and a novel tensor-network algorithm to simulate the steady-state thermodynamics.Using the example of a quantum-dot heat engine, we demonstrate that our technique replicates the well-known Landauer-Büttiker theory for efficiency and power.We then go beyond the quadratic limit to demonstrate the capability of our method by simulating a three-site machine with nonquadratic interactions.Remarkably, we find that such interactions lead to power enhancement, without being detrimental to the efficiency.Furthermore, we demonstrate the capability of our method to tackle complex many-body systems by extracting the superdiffusive exponent for high-temperature transport in the isotropic Heisenberg model.Finally, we discuss transport in the gapless phase of the anisotropic Heisenberg model at a finite temperature and its connection to charge conjugation parity, going beyond the predictions of single-site boundary driving configurations.
Strongly interacting fermions underpin some of the most challenging problems in condensed matter physics, such as high-temperature superconductivity. The low-energy states of these systems encode their essential microscopic properties, yet remain largely inaccessible to classical methods. Quantum simulation offers a promising path forward, and among state-preparation strategies, engineered dissipation has emerged as a particularly compelling approach. Existing cooling protocols, however, typically rely on knowledge of the quasiparticle spectrum or mappings to free-fermion limits. In this letter, we introduce a randomized, symmetry-preserving cooling algorithm that requires no spectral information, using only local coupling operators to ancilla degrees of freedom with randomly sampled energy splittings to drive generic fermionic systems toward their low-energy manifold. We benchmark the protocol on canonical correlated fermionic models relevant to high-temperature superconductors, spanning metallic, density-wave, paired, superconducting, and phase-separated phases. Across all models, we observe universal cooling behavior: monotonic energy relaxation, concentration of spectral weight at low energies, and stabilization of correlated ground-state order. Our results establish randomized dissipative cooling as a general strategy for preparing strongly correlated fermionic states on programmable quantum devices.