Fermionic Bath Research Articles

Self-oscillating open quantum systems

Evolution of occupation number is studied for a bosonic oscillator (with one and two degrees of freedom) linearly fully coupled to fermionic and bosonic heat baths. The lack of equilibrium in this oscillator is discussed in the light of the creation of an energy source. The connection of such a system with known nonlinear self-oscillating systems is shown.

Relaxation to persistent currents in a Hubbard trimer coupled to fermionic baths

Relaxation to persistent currents in a Hubbard trimer coupled to fermionic baths

Floquet Fermi Liquid.

We demonstrate the existence of a nonequilibrium "Floquet Fermi liquid" state arising in partially filled Floquet Bloch bands weakly coupled to ideal fermionic baths, which possess a collection of "Floquet Fermi surfaces" enclosed inside each other, resembling matryoshka dolls. We elucidate several properties of these states, including their quantum oscillations under magnetic fields which feature slow beating patterns of their amplitude reflecting the different areas of the Floquet Fermi surfaces, consistent with those observed in microwave induced resistance oscillation experiments. We also investigate their specific heat and thermodynamic density of states and demonstrate how by controlling properties of the drive, such as its frequency, one can tune some of the Floquet Fermi surfaces toward nonequilibrium Van Hove singularities without changing the electron density.

Perturbative solution of fermionic sign problem in quantum Monte Carlo computations

We have developed a strong-coupling perturbation scheme for a general doped Hubbard model around a particle-hole-symmetric reference system, which is free from the fermionic sign problem. Our approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in both continuous and discrete time versions for large periodic clusters in a fermionic bath. By considering the first-order perturbation in the shift of the chemical potential and the second-neighbor hopping, we are able to obtain an accurate electronic spectral function for a range of parameters that correspond to optimally doped cuprate systems at temperatures of up to T = 0.1t, which are challenging to access using straightforward lattice QMC calculations. We also discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system with the interaction parameter U equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter t′\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${t}^{{\\prime} }$$\\end{document} for high-temperature superconducting cuprates.

Stochastic Schrödinger Equation Approach to Real-Time Dynamics of Anderson-Holstein Impurities: An Open Quantum System Perspective.

We develop a stochastic Schrödinger equation (SSE) framework to simulate the real-time dynamics of Anderson-Holstein (AH) impurities coupled to a continuous Fermionic bath. The bath degrees of freedom are incorporated through fluctuating terms determined by exact system-bath correlations, which is derived in an ab initio manner. We show that such an SSE treatment provides a middle ground between numerically expansive microscopic simulations and macroscopic master equations. Computationally, the SSE model enables efficient numerical methods for propagating stochastic trajectories. We demonstrate that this approach not only naturally provides microscopically detailed information unavailable from reduced models but also captures effects beyond master equations, thus serving as a promising tool to study open quantum dynamics emerging in physics and chemistry.

Signature of topology via heat transfer analysis in the Su–Schrieffer–Heeger (SSH) model

In this work, we explore how thermodynamics can be a potential tool for identifying the topological phase transition. Specifically, we focus on a one-dimensional Su–Schrieffer–Heeger (SSH) chain sandwiched between two fermionic baths. To investigate distinctive thermodynamic signatures associated with the topological phase, we employ heat flow analysis. Our results, derived using a global master equation, unveil a significant suppression of heat flow as we transition from the trivial to the topological phase. This decline in heat flow can be attributed to the reduction in transmission coefficients of non-zero energy modes within the topological phase. It may serve as an indicator of a phase transition. Furthermore, we investigate the heat flow asymmetry to search for phase transition indicators. Interestingly, no asymmetry is observed when employing fermionic baths. However, upon substituting fermionic baths with bosonic ones, we report a non-zero heat flow asymmetry. For SSH model with few fermionic sites, this asymmetry is more pronounced in the topological phase compared to the trivial phase. Therefore, the observed behavior of the heat diode provides an additional means of distinguishing between the topological and trivial phases. Finally, we delve into the contributions from both bulk and edge effects in heat flow and rectification to explore the impact of small system sizes on our findings.

Strongly coupled fermionic probe for nonequilibrium thermometry

We characterise the measurement sensitivity, quantified by the quantum Fisher information (QFI), of a single-fermionic thermometric probe strongly coupled to the sample of interest, a fermionic bath, at temperature T. For nonequilibrium protocols, in which the probe is measured before reaching equilibrium with the sample, we find new behaviour of the measurement sensitivity arising due to non-Markovian dynamics. First, we show that the QFI displays a highly non-monotonic behaviour in time, in contrast to the Markovian case where it grows monotonically until equilibrium, so that non-Markovian revivals can be exploited to reach a higher QFI. Second, the QFI rate is maximised at a finite interrogation time t∗ , which we characterize, in contrast to the solution t∗→0 known in the Markovian limit (Pavel Sekatski and Martí Perarnau-Llobet 2022 Quantum 6 869). Finally, we consider probes make up of few fermions and discuss different collective enhancements in the measurement precision.

Theoretical study of nonadiabatic hydrogen atom scattering dynamics on metal surfaces using the hierarchical equations of motion method.

Hydrogen atom scattering on metal surfaces is investigated based on a simplified Newns-Anderson model. Both the nuclear and electronic degrees of freedom are treated quantum mechanically. By partitioning all the surface electronic states as the bath, the hierarchical equations of motion method for the fermionic bath is employed to simulate the scattering dynamics. It is found that, with a reasonable set of parameters, the main features of the recent experimental studies of hydrogen atom scattering on metal surfaces can be reproduced. Vibrational states on the chemisorption state whose energies are close to the incident energy are found to play an important role, and the scattering process is dominated by a single-pass electronic transition forth and back between the diabatic physisorption and chemisorption states. Further study on the effects of the atom-surface coupling strength reveals that, upon increasing the atom-surface coupling strength, the scattering mechanism changes from typical nonadiabatic transitions to dynamics in the electronic friction regime.

Fermi-Dirac staircase occupation of Floquet bands and current rectification inside the optical gap of metals: An exact approach

We consider a model of a Bloch band subjected to an oscillating electric field and coupled to a featureless fermionic heat bath, which can be solved exactly. We demonstrate rigorously that in the limit of vanishing coupling to this bath (so that it acts as an ideal thermodynamic bath) the occupation of the Floquet band is not a simple Fermi-Dirac distribution function of the Floquet energy, but instead it becomes a ``staircase'' version of this distribution. We show that this distribution generically leads to a finite rectified electric current within the optical gap of a metal even in the limit of vanishing carrier relaxation rates, providing a rigorous demonstration that such rectification is generically possible and clarifying previous statements in the optoelectronics literature. We show that this current remains non-zero even up to the leading perturbative second order in the amplitude of electric fields, and that it approaches the standard perturbative expression of the Jerk current obtained from a simpler Boltzmann description within a relaxation time approximation when the frequencies are small compared to the bandwidth.

Quantum thermodynamics with fast driving and strong coupling via the mesoscopic leads approach

Understanding the thermodynamics of driven quantum systems strongly coupled to thermal baths is a central focus of quantum thermodynamics and mesoscopic physics. A variety of different methodological approaches exist in the literature, all with their own advantages and disadvantages. The mesoscopic leads approach was recently generalised to steady state thermal machines and has the ability to replicate Landauer B\"uttiker theory in the non-interacting limit. In this approach a set of discretised lead modes, each locally damped, provide a markovian embedding for the baths. In this work we further generalise this approach to incorporate an arbitrary time dependence in the system Hamiltonian. Following a careful discussion of the calculation of thermodynamic quantities we illustrate the power of our approach by studying several driven mesoscopic examples coupled to finite temperature fermionic baths, replicating known results in various limits. In the case of a driven non interacting quantum dot we show how fast driving can be used to induce heat rectification.

Superconductivity and fermionic dissipation in quantum Hall edges

Proximity-induced superconductivity in fractional quantum Hall edges is a prerequisite to proposed realizations of parafermion zero modes. A recent experimental work [G\"ul et al., Phys. Rev. X 12, 021057 (2022)] provided evidence for such coupling, in the form of a crossed Andreev reflection signal, in which electrons enter a superconductor from one chiral mode and are reflected as holes to another, counterpropagating chiral mode. Remarkably, while the probability for crossed Andreev reflection was small, it was stronger for $\ensuremath{\nu}=1/3$ fractional quantum Hall edges than for integer ones. We theoretically explain these findings, including the relative strengths of the signals in the two cases and their qualitatively different temperature dependencies. An essential part of our model is the coupling of the edge modes to normal states in the cores of Abrikosov vortices induced by the magnetic field, which provide a fermionic bath. We find that the stronger crossed Andreev reflection in the fractional case originates from the suppression of electronic tunneling between the fermionic bath and the fractional quantum Hall edges. Our theory shows that the mere observation of crossed Andreev reflection signal does not necessarily imply the presence of localized parafermion zero modes, and suggests ways to identify their presence from the behavior of this signal in the low-temperature regime.

Correlated Mott insulators in a strong electric field: The effects of phonon renormalization

We characterize the response of a Mott insulating system to a static electric field in terms of its conducting and spectral properties. Dissipation is included by a coupling to fermionic baths and to either optical or acoustic phonons. This paper extends and completes the analysis made in a previous work by the authors [Mazzocchi et al., Phys. Rev. B 106, 125123 (2022)]. In the present work, phonons are included diagrammatically within the Migdal approximation by also including self-consistency from the electronic feedback. The nonequilibrium steady state is addressed by means of the dynamical mean-field theory based on the nonequilibrium Green's function approach, while the so-called auxiliary master equation approach is employed as impurity solver. With optical phonons the self-consistency suppresses the steady-state current for field strengths comparable to the band gap with respect to the non-self-consistent case. This is due to the interaction of phonons with the hot electrons of the lattice which increases their temperature, thus providing a less effective relaxation channel for the current-induced Joule heat. In addition, in the case of optical phonons, the results are essentially independent of the temperature of the fermionic baths, as the latter is sensibly smaller than their characteristic frequency. On the other hand, with acoustic phonons the steady-state current is slightly suppressed by the self-consistent treatment only at field strengths close to half of the gap and especially at very small phonon frequency. Also, in this case, the results seem to slightly depend on the temperature of the fermionic baths.

Cooling with fermionic thermal reservoirs.

The quantum reservoirs commonly considered in open-quantum systems theory are those modeled by quantum harmonic oscillators, which are called bosonic reservoirs. Recently, quantum reservoirs modeled by two-level systems, the so-called fermionic reservoirs, have received attention due to their features. Given that the components of these reservoirs have a finite number of energy levels, unlike bosonic reservoirs, some studies are being carried out to explore the advantages of using this type of reservoir, especially in the operation of heat machines. In this paper, we carry out a case study of a quantum refrigerator operating in the presence of bosonic or fermionic thermal reservoirs, and we show that fermionic baths have advantages over bosonic ones.

Non-Fermi liquid behavior and quantum criticality in cubic heavy fermion systems with non-Kramers multipolar local moments

Notable non-Fermi liquid and quantum critical behaviors are observed in rare-earth metallic systems with non-Kramers local moments supporting a number of different multipolar moments. A prominent example is $\text{Pr(Ti,V)}_{2}\text{Al}_{20}$, where the non-Kramers doublet of the $\text{Pr}^{3+}$ ion allows quadrupolar and octupolar moments, but lacks a dipolar moment. Previous theoretical studies show that a single impurity Kondo problem with such an unusual local moment leads to novel non-Fermi liquid states. In this work, we investigate possible quantum critical behaviors arising from the competition between non-Fermi liquid states and multipolar-ordered phases induced by the RKKY interaction. We consider a local version of the corresponding Kondo lattice model, namely the Bose-Fermi Kondo model. Here, the multipolar local moments are coupled to fermionic and bosonic bath degrees of freedom representing the multipolar Kondo effect and RKKY interactions. Using a perturbative renormalization group (RG) study up to two loop order, we find critical points between non-Fermi liquid Kondo fixed points and a quadrupolar ordered fixed point. The critical points describe quantum critical behaviors at the corresponding phase transitions and can be distinguished by higher order corrections in the octupolar susceptibility that can be measured by ultrasound experiments. Our results imply the existence of a rich expansion of the phases and quantum critical behaviors in multipolar heavy fermion systems.

Correlated Mott insulators in strong electric fields: Role of phonons in heat dissipation

We study the spectral and transport properties of a Mott insulator driven by a static electric field into a non-equilibrium steady state. For the dissipation, we consider two mechanisms: Wide-band fermion reservoirs and phonons included within the Migdal approximation. The electron correlations are treated via non-equilibrium dynamical mean field theory with an impurity solver suitable for strong correlations. We find that dissipation via phonons is limited to restricted ranges of field values around Wannier-Stark resonances. To cover the full range of field strengths, we allow for a small coupling with fermionic baths, which stabilizes the solution. When considering both dissipation mechanisms, we find that phonons enhance the current for field strengths close to half of the gap while lowering it at the gap resonance as compared to the purely electronic dissipation used by Murakami and Werner [arXiv:1804.08257]. Once phonons are the only dissipation mechanism, the current in the metallic phase is almost one order of magnitude smaller than the typical values obtained by coupling to a fermionic bath. In this case, the transport regime is characterized by an accumulation of charge in the upper Hubbard band described by two effective chemical potentials.

Quantum Collision Models: A Beginner Guide.

In recent years, quantum collision models, sometimes dubbed repeated interaction models, have gained much attention due to their simplicity and their capacity to convey ideas without resorting to technical complications typical of many approaches and techniques used in the field of open quantum systems. In this tutorial, we show how to use these models, highlighting their strengths and some technical subtleties often overlooked in the literature. We do this by deriving the Markovian master equation and comparing the standard collisional derivation with the standard microscopic one. We then use the collision model to derive the master equation of a two-level system interacting with either a bosonic or fermionic bath to give the reader a flavour of the real use of the model.

Nonequilibrium scalar field dynamics starting from Fock states: Absence of thermalization in one-dimensional phonons coupled to fermions

We propose a new method to study non-equilibrium dynamics of scalar fields starting from non-Gaussian initial conditions using Keldysh field theory. We use it to study dynamics of phonons coupled to non-interacting bosonic and fermionic baths, starting from initial Fock states. We find that in one dimension long wavelength phonons coupled to fermionic baths do not thermalize both at low and high bath-temperatures. At low temperature, constraints from energy-momentum conservation lead to a narrow bandwidth of particle-hole excitations and the phonons effectively do not see this bath. On the other hand, the strong band-edge divergence of the particle-hole density of states leads to an undamped polariton-like mode of the dressed phonons above the band edge of the particle-hole excitations. These undamped modes contribute to the lack of thermalization of long wavelength phonons at high temperatures. In higher dimensions, these constraints and the divergence of density of states are weakened and lead to thermalization at all wavelengths.

Non-equilibrium dynamics of a dissipative two-site Hubbard model simulated on IBM quantum computers

Many-body physics is one very well suited field for testing quantum algorithms and for finding working heuristics on present quantum computers. We have investigated the non-equilibrium dynamics of one- and two-electron systems, which are coupled to an environment that introduces decoherence and dissipation. In our approach, the electronic system is represented in the framework of a two-site Hubbard model while the environment is modelled by a spin bath. To simulate the non-equilibrium population probabilities of the different states on a quantum computer we have encoded the electronic states and environmental degrees of freedom into qubits and ancilla qubits (bath), respectively. The total evolution time was divided into short time intervals, during which the system evolves. After each of these time steps, the system interacts with ancilla qubits representing the bath in thermal equilibrium. We have specifically studied spin baths leading to both, unital and non-unital dynamics of the electronic system and have found that electron correlations clearly enhance the electron transfer rates in the latter case. For short time periods, the simulation on the quantum computer is found to be in very good agreement with the exact results if error mitigation methods are applied. Our method to simulate also non-unitary time-evolution on a quantum computer can be well extended to simulate electronic systems in correlated spin baths as well as in bosonic and fermionic baths.

Stochastic dynamics and bound states of heavy impurities in a Fermi bath

We investigate the dynamics of heavy impurities embedded in an ultracold Fermi gas by using a generalized Langevin equation. The latter---derived by means of influence functional theory---describes how the stochastic classical dynamics of the impurities and the quantum nature of the fermionic bath manifests in the emergent interaction between the impurities and in the viscosity tensor. By focusing on the two-impurity case, we predict the existence of bound states, in different conditions of coupling and temperature, whose lifetime can be analytically estimated. Our predictions should be testable using cold-gases platforms within current technology.

Environment-induced decay dynamics of antiferromagnetic order in Mott-Hubbard systems

We study the dissipative Fermi-Hubbard model in the limit of weak tunneling and strong repulsive interactions, where each lattice site is tunnel-coupled to a Markovian fermionic bath. For cold baths at intermediate chemical potentials, the Mott insulator property remains stable and we find a fast relaxation of the particle number towards half filling. On longer time scales, we find that the anti-ferromagnetic order of the Mott-N\'eel ground state on bi-partite lattices decays, even at zero temperature. For zero and non-zero temperatures, we quantify the different relaxation time scales by means of waiting time distributions which can be derived from an effective (non-Hermitian) Hamiltonian and obtain fully analytic expressions for the Fermi-Hubbard model on a tetramer ring.