Quantum-Zeno Fermi polaron in the strong dissipation limit
Tomasz Wasak,1 Richard Schmidt,2 and Francesco Piazza1 1Max-Planck-Institut f¨urPhysik komplexer Systeme, N¨othnitzer Str. 38, 01187 Dresden, Germany 2Max-Planck-Institut f¨urQuantenoptik, Hans-Kopfermann-Str. 1, 85748 Garching, Germany The interplay between measurement and quantum correlations in many-body systems can lead to novel types of collective phenomena which are not accessible in isolated systems. In this work, we merge the Zeno-paradigm of quantum measurement theory with the concept of polarons in condensed-matter physics. The resulting quantum-Zeno Fermi-polaron is a quasi-particle which emerges for lossy impurities interacting with a quantum-degenerate bath of fermions. For loss rates of the order of the impurity-fermion binding energy the quasi-particle is short lived. However, we show that in the strongly dissipative regime of large loss rates a long-lived polaron branch re- emerges. This quantum-Zeno Fermi-polaron originates from the nontrivial interplay between the Fermi-surface and the surface of the momentum region forbidden by the quantum Zeno projection. The situation we consider here is realized naturally for polaritonic impurities in charge-tuneable semiconductors and can be also implemented using dressed atomic states in ultracold gases.
I. INTRODUCTION tion in mind: how are the many-body states affected by the emergence of a dissipation-induced constraint for the Hilbert space as the quantum-Zeno limit is approached? The effect of measurement on the time-evolution of a We focus on the paradigm of mobile quantum impurity system is one of the most puzzling aspects of quantum problems at the border between few- and many-body dynamics [1,2], making it drastically distinct from its physics and extend it to the regime of gain and loss classical counterpart. One paradigmatic example is the of impurity particles. We demonstrate how the intrin- quantum Zeno effect [3–5], whereby a frequent-enough sic correlation scales, governing the impurity problem, measurement of a quantum system results in a locali- have a nontrivial interplay with the scale brought in by sation of its state in Hilbert space [6]. Increasing the the dissipation-induced constraint, leading to qualitative measurement rate from zero, the system dynamics first modification of the many-body quantum states. gains a non-unitary character due to the loss of infor- The concept of polarons plays a central role in con- mation to the measurement device. However, above a densed matter physics, developed to understand the mo- threshold-rate the trend is reversed: the dynamics tends tion of an electron moving in a dielectric crystal [25–28]. asymptotically to a unitary limit within a Hilbert space The Fermi polaron, describing a mobile impurity inter- reduced by the region that is forbidden due the measure- acting with a bath of Fermions [29], is particularly im- ment process. portant for the quantitative description of the physics of Recently, growing attention is being devoted to the ef- imbalanced Fermi mixtures [30, 31] and sets the basis fect of repeated measurements on the dynamics of quan- for our understanding of other strongly correlated sys- tum many-body systems [7–13]. This has been motivated tems [32–35]. The physics of Fermi polarons has become by experimental developments in controlling and selec- the center of significant interest due to recent experi- tively coupling complex quantum systems to their envi- mental advances that allowed for its observation both ronment, as demonstrated with superconducting qubits with ultracold atoms [36–38] and more recently in charge- [14–16], ultracold atoms [17–21] and molecules [22, 23]. tunable monolayer semiconductors [39, 40]. The interplay between measurement and intrinsic quan- So far, theoretical studies of the polaron paradigm are tum correlations in many-body systems can lead to novel assuming a conserved number of impurities in a closed collective phenomena. Recently, two instances of such system. However, in addition to the fundamental interest collective phenomena appearing in the quantum-Zeno in exploring the effect of strong measurement on quantum limit have been predicted. First, analyzing the entangle- many-body systems in general, the recent experiments ment entropy, it was found that by increasing the mea- in charge-tunable monolayer semiconductors open a new surement rate of a closed quantum many-body system, frontier to study polarons in the presence of impurity loss a phase transition from a volume-law entangled phase to and gain, as the corresponding rates can become compa- arXiv:1912.06618v2 [cond-mat.quant-gas] 13 Jan 2021 a quantum Zeno phase with area-law entanglement can rable to the timescales of unitary dynamics [39]. In this take place [11, 13]. Second, a study of one-dimensional work we develop a theoretical approach to describe such fermionic models with a localized dissipation (of Dirac quantum driven-dissipative impurity dynamics. Using a delta type) revealed that the quantum Zeno projection is diagrammatic, real-time Keldysh approach we show that, enhanced in the vicinity of the Fermi surface due to the while for intermediate loss rates the polaron becomes ex- existence of gapless modes [24]. pectedly short-lived, in the Quantum Zeno limit of large In this work, we consider a situation where a strong loss rates a long-lived polaron branch re-emerges. This dissipation acts on a subspace of the Hilbert space of Quantum-Zeno Fermi polaron is not simply a copy of an interacting quantum system, with the following ques- the original polaron in absence of losses, but originates 2 from the nontrivial interplay between the Fermi-surface repumped into the system. To model the action of dissi- and the surface of the region forbidden by the quantum pation and pump, we employ the master equation Zeno-projection. This indeed leads to qualitative modifi- ∂ %ˆ(t) = i[H,ˆ %ˆ(t)] + %ˆ(t). (3) cations of the dispersion of the impurity-fermion bound t − Ld states as well as of the scattering continuum. The dissipative Lindblad operator d, acting on the den- This complex interplay underlying the formation of the sity matrix of the whole system% ˆ(t),L splits into two parts, Quantum-Zeno Fermi-polaron requires a loss profile that P i.e., d%ˆ = k γ(k)D[ˆck]+Ω(k)P [ˆck] %ˆ, which describe can single-out a non-trivial region in momentum space. lossL of particles{ with rate γ(k) and their} reinjection with This can be achieved by mixing the impurity with an rate Ω(k). The Markov form of the master equation (3) additional lossy degree of freedom of much smaller mass. is valid as long as the reservoir of impurity-particles has This is naturally the case for exciton-polariton impurities a correlation time much shorter than all other timescales in charge-tunable monolayer semiconductors. It can also characterizing the system [1]. We model here the dis- be engineered for ultracold atoms by coupling a stable sipation strength by a momentum dependent profile, as electronic level to a short-lived level, and additionally it may be generated by a coupling of an impurity to a using level-depending trapping potentials. lossy degree of freedom, for example excitons to micro- This manuscript is structured as follows. In SectionII, cavity photons or ground state atoms to short-lived ex- we introduce a model describing the driven-dissipative cited states. Fermi-polaron. In SectionIII, we discuss the results for We consider loss and injection of single impurities. the spectral response of the system. We show the emer- These two incoherent processes are expressed by the gence of the quantum-Zeno Fermi-polaron for large loss 1 super-operators D[ˆc]ˆ% cˆ†%ˆcˆ cˆ†c,ˆ %ˆ and P [ˆc] rates and illustrate the interplay between the surface of ≡ − 2 ≡ D[ˆc] + D[ˆc ]. The action of each super-operator gener- the Zeno-forbidden region and the Fermi surface. In Sec- † ates two terms. The anti-commutator term corresponds tionIV, we discuss in detail the concrete implementation to a non-Hermitian Hamiltonian evolution describing the of the Fermi-polaron subject to strong dissipation. Fi- damping/growth of the number of impurities, while the nally, in SectionV, we offer some concluding remarks other term describes the noise which necessarily accom- and an outlook. Technical details are presented in Ap- panies dissipation. In a stochastic quantum trajectory pendices, where, in particular, we introduce a Keldysh picture [1], the noise-term generates projections of the diagrammatic approach to non-equilibrium Green’s func- wave function at random times induced by the measure- tions. ment, which, in our case, corresponds to the detection of an impurity particle exiting or entering the system. In this picture, both the measurement and the damping (or II. DRIVEN-DISSIPATIVE IMPURITY IN A growth) rate are set by γ(k) (or Ω(k)). For the present FERMI BATH study, the specific form of the pump term P is irrelevant since the excitation spectrum of the system does not de- We consider a system composed of a bath of fermions pend on it as long as the average density of impurities is interacting with a vanishingly small density of impurity small. particles. The Hamiltonian of the system consists of In the limit of vanishing dissipation and drive, our three parts: model corresponds to the Fermi-polaron problem treated as a closed system. One approach to describe the pola- Hˆ = Hˆc + Hˆf + Hˆint, (1) ronic states of an impurity coupled to a bath of fermions in this closed system, is to employ a variational approach ˆ P ˆ P ˆ ˆ where Hc = k εc(k)ˆck† cˆk and Hf = k εf (k)fk†fk are based on the so-called Chevy ansatz for the wave func- free-particle Hamiltonians with kinetic energies εc(k) = tion [30, 35, 44–47]. In this ansatz the impurity motional 2 2 k /2mc and εf (k) = k /2mf of impurities and fermions, degree of freedom is entangled with a single particle- respectively. We assume that the interaction is a contact hole excitation of the Fermi sea. The variation with potential characterized by a coupling strength U, i.e., respect to the expansion coefficients yields an equation that describes the properties of the lowest lying polaron Z state, the so-called attractive polaron. Further varia- Hˆ = U dr cˆ (r)ˆc(r)fˆ (r)fˆ(r). (2) int † † tional methods, including higher-order particle-hole exci- tations, have also been developed to improve the Chevy In this manuscript, all concrete calculations will be given expansion [31]. The variational approach [48] and a for the two-dimensional case which has been realized in related non-self-consistent T -matrix approach [49] were ultracold atoms [38, 41] and two-dimensional semicon- used to capture the properties of an excited state of the ductors [39, 40]. polaron, known as the repulsive polaron. All these theo- Driven-dissipative impurities which are localized and retical studies rely on the Hamiltonian description of the immobile in a one-dimensional system have been recently system and, therefore, are not directly applicable to the studied in Refs. [42, 43]. Here we consider mobile impu- case of driven-dissipative impurities. rities which are subject to loss and thus require to be In this work, we develop a diagrammatic approach 3 based on non-equilibrium Green’s functions (GFs) which given by the following formula: extends the non-self-consistent T -matrix approach to (k, ω) = 2Im[GR(k, ω)]. (5) driven-dissipative systems. We employ a path-integral Ac − c formulation of our problem on the Keldysh real-time con- tour which allows for a straightforward inclusion of the The spectral function satisfies the exact sum rule R (p)dω/(2π) = 1. It exhibits resonant peaks in cor- Lindblad terms of the master equation [50, 51]. Within c respondenceA to quasi-particle excitations with a width this formalism, diagrammatic approaches are applicable related to their inverse lifetime. In the absence of in- in a way analogous to the equilibrium Matsubara for- teractions with the fermions, the impurity spectral func- malism. In the following, we analyze the properties of tion = 2Im[GR ] has resonance at its bare dis- the system in its steady state, where the GFs are time- 0,c 0,c persionA ω = −ε (k), with a width determined by γ(k). translation invariant. The details of the calculations can c In presence of interactions, resonances appear at the be found in AppendixA. renormalized frequencies ω = ε (k) Re[ΣR(k, ω)]. In what follows, we will focus on the impurity limit c c In this case, the inverse lifetime is set− by the imagi- where the number of impurities in the system is van- nary part of the pole position on the complex plane: ishingly small and the excitation spectrum of the system Im[ω(k)] satisfying [GR(k, ω(k))] 1 = 0. For sufficiently will not depend on the pump Ω(k) or the quantum statis- c − large inverse-lifetime the latter is approximately given by tics of impurities. R 1 R Im [Gc (k, ω)]− = Im[ Σc (k, ω) + iγ(k)/2]. An{ example of} the impurity− spectral function for grad- ually increasing γ is presented in Fig.1 for the mass bal- III. FERMI-POLARON WITH A LOSSY 0 REGION IN IMPURITY-MOMENTUM SPACE anced case. For small γ0, two quasi-particle resonances are visible. The lower one in energy is the so-called at- tractive polaron. It corresponds to a state where the A. Impurity response function: quantum-Zeno impurity attracts a cloud of surrounding fermions. The Fermi-polaron attractive polaron is long-lived (the width of the peak vanishing for k 0) and ceases to exist for large mo- As we shall see, the emergence of a quantum-Zeno menta, when the| | resonance→ enters the molecule-hole con- Fermi-polaron relies on impurity losses which single out tinuum. The latter appears as blurred region between the at least one subset of momentum space. Such a subset two polaron resonances. The second resonance is higher we call the dissipative subspace (DS). In this section, we in energy and corresponds to a repulsive polaron state consider the simplest case of a loss function γ(k) that is where the surrounding cloud is repelled by the impurity. characterized by a single scale given by a cutoff momen- For low-enough losses and small enough momenta, the re- tum k , and we assume that γ(k) = γ if k < k and is γ 0 | | γ pulsive polaron has a shorter lifetime than the attractive vanishing otherwise. Our approach is valid for any value polaron. of the loss rate γ0. However, we will be interested in the It is important to note that at any finite value of γ0, case of large loss rates: γ E with E the impurity- 0 B B while the losses affect bare impurity only within the DS fermion binding energy. This leads to effective blocking defined by k < kγ , they do affect the polaron also out- of transitions between states involving the DS, which can side the DS.| The| reason for this is that the polarons are be interpreted in terms of the quantum Zeno effect. The impurities dressed by the interactions with the surround- exact form of the loss profile is not crucial for the physics ing fermions and therefore involve a superposition of all described in what follows, and it is sufficient that the the allowed momentum states. This statement is partic- loss rate decays fast enough over the momentum scale p ularly explicit in the simple variational ansatz introduced qB = mf EB. We will demonstrate this by considering at the end of AppendixA2. Impurity losses thus affect experimentally realistic cases in Sec.IV. Finally, we re- the polarons in a non-trivial manner via the coupling be- mark that the coupling strength U is eliminated from the tween the impurity and the molecular state. The latter is equations by finding the impurity-fermion binding energy very sensitive to dissipation. We note that for the value EB (based on the poles of the vacuum T -matrix) in the of kγ used in Fig.1 and for momenta outside the DS, the limit of vanishing fermion density and with dissipation repulsive polaron remains largely unaffected by the losses strength set to zero. (see also Fig.2). This is due to small coupling of the im- The spectral response of the system is encoded in the purity with the bound state at high energies, where the retarded impurity GF, which is given by repulsive polaron recovers the impurity dispersion rela- tion. The effect of losses outside the DS is most clearly R 1 visible in panels (c) and (d) in Fig.1, where, for γ0 EB, Gc (k, ω) = k , (4) ω ε (k) + i γ( ) ΣR(p) + i0 the attractive-polaron peak is almost fully washed∼ out. − c 2 − c Quite remarkably, for even larger losses γ0 EB, a R where the retarded self-energy Σc (p) of the impurity resonance reappears below the repulsive polaron just out- contains information about the interaction with the sur- side the DS: k > kγ . This resonance has a width which rounding bath of fermions; here and below the four- does not depend anymore on the loss rate γ0 and is rather momentum p (k, ω). The spectral function is then given only by the coupling to the molecule-hole contin- ≡ 4
c 6 (a) (b) (c) (d) (e) A 8 4 B 6
/E 2 4 ! 0 quantum-Zeno Fermi-polaron 2 2 0 =0 0 =0.5EB 0 =1.0EB 0 =4.0EB region forbidden 0 = 20EB 0 1 2 3 0 1 2 3 0 1 2 3 0 1 2 3 0 1 2 3 k/qB k/qB k/qB k/qB k/qB
FIG. 1. The impurity spectral function Ac(k, ω) (in the units of EB ) for a 2D system with the mass balance mc = mf . In the panels from (a)–(e) the strength of the loss is increased: γ0/EB = 0, 0.5, 1.0, 4.0, 20. The energy is expressed in the units of p 2 the impurity-fermion binding energy EB , and momenta are given in qB = mf EB . The Fermi energy F = kF /2mf is set by the Fermi wavevector kF = 0.7qB . The color coding of the values for Ac is on the right. The cutoff momentum of the loss 2 profile is at kγ = qB . The dashed red curve is the free dispersion relation ω = εc(k) ≡ k /2mc.
c 6 (a) (b) (c) (d) (e) A 8 4 B 6 2 /E 4 ! 0 2 2 k =0 k =0.5qB k =1.0qB k =1.5qB k =2.0qB 0 1 2 3 0 1 2 3 0 1 2 3 0 1 2 3 0 1 2 3 0 k/qB k/qB k/qB k/qB k/qB
−1 FIG. 2. Impurity spectral function Ac(k, ω) (in units of EB ) as a function of |k| and ω (in units of qB and EB , respectively) 2 for different sizes of the dissipative subspace: kγ /qB = 0, 0.5, 1.0, 1.5, 2.0 [from (a) to (e)]. The Fermi energy F = kF /2mf is set by kF = qB , and γ0 = 20 EB . The dotted horizontal line depicts ω = 0. The dashed vertical line [in the panels (b)–(e)] 2 displays kγ . The dashed red parabolic curve shows the dispersion relation ω = k /2mc of the free impurity. uum, as would be the case without losses. At the same of k compared to the lossless case and is also blueshifted. time, the DS does not contain any spectral weight and is This is evident from Fig.2(d) where the resonance ap- fully forbidden for the excitations. This phenomenology pears also for ω > 0. In addition, the spectral weight can be understood within the paradigm of the quantum- of the repulsive polaron is transferred to the attractive Zeno effect. Indeed, when the impurity occupies a state branch. By increasing kγ the attractive polaron peak be- with k DS it is lost immediately, which is equivalent comes broader due to the entrance into the molecule-hole to a frequent∈ measurement performed within the DS. As continuum. a result, the dynamics of the particle is confined to the orthogonal subspace, i.e. k / DS. We note that for in- finite dissipation strength, the∈ equation for the pole of B. Molecular response function: quantum-Zeno the impurity Green’s function can be reproduced by a scattering modification of the Chevy polaron ansatz for the wave function in which the momenta from DS are excluded. In the interaction between an impurity and a fermion, the formation of paired (or molecular) states plays a cru- The loss rate γ0 does not appear explicitly anymore as the dynamics becomes unitary again in the quantum cial role. Our approach based on an auxiliary molecular Zeno limit γ E (see also Fig.4 in the next section). field allows to rewrite the impurity self-energy in the fol- 0 B lowing way, see AppendixA: Still, the size of the DS kγ remains as further scale which plays a central role in the problem, as does the Fermi R 1 X R Σc (k, ω) = G∆ k + k0, ω + εf (k0) , (6) wavevector kF . One effect of the non-trivial interplay V k0 FS between these scales is demonstrated in Fig.2, where we ∈ show how the polaron dispersion depends on the size of from which we see that the renormalization of the impu- the DS. By increasing kγ the attractive-polaron disper- rity is fully determined by the molecular GF. The molec- sion gradually changes from quadratic to linear. More- ular GF is at the same time directly related to the T - over, the branch remains well defined up to larger values matrix describing the scattering between the impurity 5
(a) k
k
k
FS
γ0 & EB. The increase of the DS leads to a fast increase black lines in Fig.5a,b). of the threshold value of γ0,min. Comparing this to the For finite loss rates (blue-dashed lines), the jump in Fig.4a, we see that to reach the Zeno limit when s 1 the T -matrix-phase has a local maximum at the bound | | ≈ requires an order of magnitude of increase of γ0 as com- state but does not reach the maximal value π. By fur- pared to γ0,min. We note that in this figure, kF = 0 and ther increasing the energy the T -matrix-phase then grows the value of ω is set to its limiting case ωth. For larger ω again to reach π at the threshold value and then decreases the value of s increases, and similarly s grows with kF . monotonically and faster than in the lossless case above | | | | threshold. The T -matrix-phase behaves nontrivially upon increas- 1.0 ing the loss rate. The double-peak structure present be- 0.9 low threshold develops, in the quantum-Zeno limit of in- finite loss rate, into a single non-analytical plateau iden- 0.8
| tical to the one present in the lossless case, apart from s | 0.7 a blue shift of the threshold and the bound state. As we argued above, the vacuum scattering in the quantum- 0.6 Zeno limit is the same as the lossless scattering in pres- 0.5 (a) ence of a Fermi surface, where the role of kF is played by k . Consistently with this, the phase of the T -matrix 0.4 γ 0.001 0.010 0.100 1 10 100 1000 in the quantum-Zeno limit and in particular the scatter- 0/EB ing phase shift indeed behaves as its lossless counterpart γ/EB in the presence of a medium, which has been studied 1.0 20 0 in [53, 54, 58]. B
E In particular, one can quantify the distance from the 0.9 / 15 0.8 10 quantum-Zeno limit by the size of the jump of the T - | 0,min 5 matrix-phase across the bound state. This is shown in 0.7 γ min Fig. (5)c, where we plot Arg[ T (k, ω)] maximized over
s 0 | − 0.6 0 2 4 6 8 δω. We see that for increasing γ0, the maximum jump k 0/E/qBB of the phase of T drops, and for γ E it reaches a 0.5 / 0 ∼ B kγ EB minimum, after which the further increase of γ0 leads to 0.4 (b) recovery of the maximal jump of the phase. 0.3 The recovery of unitary scattering dynamics explains 0 2 4 6 8 the re-appearance of an attractive-polaron branch for k /qB kγ/EB large loss rates shown in Fig.1. As already discussed FIG. 4. (a) The S-matrix: |s(k, ω)| = |eiδ(k,ω)| as a function in sectionIIIA, this quantum-Zeno polaron features a clearly different dispersion affected by the presence of of the loss rate γ0 (in units of EB , computed at threshold ω = ωth for k/qB = 0.05 (solid black), 0.5 (dashed red), 1 the additional scale kγ , as demonstrated by Fig.2. (orange dotted), 2 (dotted-dashed green), 5 (double-dotted- We now turn to the discussion of the interplay between dashed blue); here kγ = 2.5qB . The dotted gray horizontal kγ and the Fermi scale kF , which is best understood by line is for infinitely fast losses, see Eq. (11). (b) The minimum analyzing the molecular spectral function. In Fig.6 we of |s | ≡ min |s(γ )| as a function of k calculated for min γ0 0 γ show the spectral function ∆ for various values of kγ k = 0.05qB . The inset of panel (b) shows the value of γ0 A at a fixed density of fermions kF = qB. Moreover, we for which |s| takes the minimum value, denoted by γ0,min. In choose a finite but large loss rate γ = 20E , where the both panels, m = m . 0 B c f scattering is already close to the quantum-Zeno regime. For kγ = 0, i.e. without losses (see. Fig.6a), the spec- A further characterization of the scattering can be ob- trum consists of a molecular bound state (in this dis- tained by analyzing the phase of the T -matrix. We com- sipationless case, we artificially set a finite lifetime for pute Arg[ T (k, ω)] for fixed k as a function of δω, where the bound state so that it is visible in the plot) and − we define ω = k2/2M + δω. In Fig.5 we consider the a continuum. In the non-self-consistent T -matrix ap- vacuum case at k = 0. In absence of losses (dotted-red proximation employed here, the boundaries of the con- lines in Fig.5a,b), the phase experiences a π-jump when tinuum can simply be understood using energy and mo- crossing the bound state. Here, the phase is zero for mentum conservation in the processes of Fig.3. This δω < EB, and is π for EB < δω < 0 , and for δω > 0 has been discussed already in previous works [49, 54], − − it decays monotonically. but we reproduce the arguments here since we will ex- In presence of losses, as discussed above, the scattering tend them below to the quantum-Zeno limit. A molecule 2 threshold is shifted by kγ , since impurities can escape with energy ω and momentum k = k can decay into only for momenta outside the DS. Therefore, the phase the impurity and a fermion. The| | conservation laws re- of the T -matrix is analytical and the scattering phase quire ω = εf (k0) + εc(k00) and k = k0 + k00. The Pauli 2 shift defined only for ω > k (see blue-dashed and solid blocking sets an additional constraint k0 > k . For γ | | F 8
40 k < mckF /mr, the continuum boundary has the form: 2 2 (kF mrk/mc) /2mr + k /2M, and for k > mckF /mr (a) it is:−k2/2M, where M = m + m is the total mass. 30 c f B Consequently, for kγ = 0, the continuum has a mini- B E mum at k = kF and ω = F . For smaller k < kF , it is E
Δ 20 blueshifted due to Pauli blocking, and for k = 0 it starts A
2 at ω = kF /2mr. 10 For a finite size of the DS and in the quantum-Zeno th limit, several new features emerge in the molecular spec- Erel trum with respect to the lossless case. For small k , the 0 γ -4 -2 0 2 4 6 8 10 spectrum is affected mostly in the vicinity of the Fermi momentum, namely in the region kF kγ < k < kF +kγ , !/EB − δω/EB as can be seen from panels (b) and (c) in Fig.6. By in- creasing kγ , the spectral weight of the continuum trans- fers to ω k2/2m + k2 /2m in the region k k . Due 3.0 (b) > γ c F f F to the presence of the DS, for smaller k than≈ k k 2.5 | F − γ | ] ] the spectrum is blueshifted more rapidly, and at k = 0 it T T 2 2.0 reaches Q /2mr, where Q = max(kF , kγ ). Notice that, [- for k lying in the interval k k < k < k + k , the 1.5 F − γ F γ
Arg bound state dispersion is linear in k. For even larger kγ , Arg[ 1.0 the effect of the Pauli blocking is enhanced by the pres- th 0.5 Erel ence of the DS, and the bound state energy is blueshifted 0.0 even more. With increasing k the bound state energy drops, and enters the continuum for larger k. Due to -4 -2 0 2 4 6 8 10 the forbidden regions in the formation of a molecule, !/EB δω/EB the impurity and the fermion must have momenta not 1.0 smaller than kγ and kF , respectively. Therefore, the molecule with the smallest energy and momentum will
⇡ (c) have k = kF kγ . Since the continuum starts now at / 0.9 2 | 2 − | kγ /2mc +kF /2mf , the molecular state can appear at this ]/
max energy also for k = kγ +kF . This gives rise to the appear- ] T 0.8 T
[- ance of two minima in the bound-state dispersion visible
in Fig.6e. This analysis makes apparent the crucial role
Arg 0.7 played by the size of the DS in the quantum-Zeno regime,
Arg[ where the molecule, and in turn, the polaron dispersion 0.6 can be modified qualitatively. The results presented in this section combined with 0.001 0.010 0.100 1 10 100 1000 Figs.1 and2 highlight an important ingredient for the 0/EB appearance of strong qualitative effects of dissipation on 0/EB molecular and polaron states. It is namely the molecu- FIG. 5. (a) Spectral function of the molecule A∆ = lar bound state which is strongly modified by the losses, −2Im[T (k, ω)] calculated for k = 0 and ω = k2/2M + δω; while the continuum experiences weaker modifications. here kF = 0, kγ = 1.5qB . The dotted red line – the vacuum, This is reflected by the behavior of the polarons: only the lossless case, cf. Eq (9). The dashed blue line – A in the attractive polaron is strongly affected by loss. This leads presence of moderate loss (γ = E ). The black line – infi- 0 B us to the expectation that the phenomenology we pro- nite loss case, cf. Eq. (11). The arrow indicates the threshold th pose here will crucially depend on the presence/absence Erel = 2.25EB . (b) The phase of the T -matrix, i.e., Arg[−T ]. The color-code and parameters are as in the panel (a). (c) of the molecular bound state and therefore on the dimen- The maximum value of Arg[−T ], calculated with respect to sionality. parameter δω, as a function of γ0. The color-coding is the Finally, we end this discussion by noting that all the same as in Fig.4a. In all the panels, mc = mf . results are obtained within the non-self-consistent T- matrix approach, whose limitations are known [49]. The full solution would require to employ self-consistency of example, setting k = 0, we have k0 = k00 and the en- the Green’s functions and would lead to modification of 2 − 2 ergy conservation leads to ω = k00 /2mr 6 kF /2mr. the molecular spectral function shifting the position of 2 | | Therefore, for ω > kF /2mr the molecule can decay into the molecular and polaron resonances and affecting their an impurity-fermion pair. Moreover, one can show that lifetimes. These corrections, however, in the configu- for k kF , the threshold for the molecular decay be- rations and regimes considered affect the quantities of haves≈ as: k2 /2m + (k k )2/2m . In general, for present interest only quantitatively and do not change F f − F c 9
4 3 A 50
B 2 40
/E 1 30 ! 0 k =0 k =0.5qB k =1.0qB k =1.5qB k =2.0qB 20 1 (a) (b) (c) (d) (e) 10 2 0 1 2 3 4 50 1 2 3 4 50 1 2 3 4 50 1 2 3 4 50 1 2 3 4 5 0 k/qB k/qB k/qB k/qB k/qB
2 FIG. 6. Molecular spectral function A∆(k, ω) (in units of qB EB ) for mc = mf as a function of |k| (in units of qB ) and ω (in units of EB ) for different sizes of the dissipative subspace: kγ /qB = 0, 0.5, 1.0, 1.5, 2.0 [from (a) to (e)]. The Fermi energy 2 F = kF /2mf is set to kF = qB , and γ0 = 20 EB . The dotted/dotted-dashed vertical lines indicate |kγ ∓ kF |. The dashed red 2 curve is the molecular free dispersion given by ω = −EB + k /2M. the overall physical picture. corresponding to the Lindblad operator X = γ D[dˆk]. (18) Ldiss d IV. IMPLEMENTATION OF THE MODEL k
The model just introduced also describes atoms featur- In this section, we discuss how the physics discussed so ing a (meta)stable internal electronic state coupled by a far could be realized in experiment. Starting from a setup laser to a short-lived excited state, so that δ is the de- featuring a density-density interaction between fermions tuning between the laser frequency and the atomic tran- and impurity particles, the goal is to realize momentum- sition frequency. A situation where the short-lived state dependent losses for the impurities. The loss profile needs has also a much smaller mass, as is naturally the case for to single out a given region of momentum space by being excitons and photons, can be implemented with atoms much larger there than elsewhere. by trapping the stable state in an optical lattice while Recently it was demonstrated that it is possible to leaving the excited state untrapped. explore the two-dimensional Fermi polaron problem in The Lindblad master equation resulting from H and transition metal dichalcogenide monolayer semiconduc- eff translates in the Keldysh approach to the inverse tor embedded into microcavities [39, 40]. In this type of diss retardedL GF: setup, excitons are deeply bound and hybridize with the cavity photons, forming polaritons. By optically exciting ω (k)+i0 Ω [GR ] 1= c R . a small number of polaritons it is possible to create impu- cd − − Ω ω [ (k) −δ]+iγ /2+i0 − R − d − d rities that feature a mass that is much smaller than that (19) of the surrounding electrons, since photons propagate Diagonalization of this matrix leads to the dispersion re- much faster. More important for the present discussion lation and loss rate of the polaritons. One of the complex is the fact that the photons feature roughly momentum- eigenvalues, which we denote by εc(k), converges to c(k) independent losses due to mirror leakage. These are for large k ; the other converges to (k)+iγ /2 for large inherited by the polaritons in a momentum-dependent | | d d k and we denote it with εd(k). The real part of the fashion, as the photonic admixture indeed depends on eigenvalue| | ε(k) describes the position of the resonance, momentum. This means that for one of the two polari- whereas its imaginary part characterizes its width, and tons the loss rate will decrease by increasing momentum, hence the loss rate. For simplicity, we assume that the while the opposite is true for the other polariton. low momenta are sufficient to describe the dynamics so The simplest model for the above scenario is given by 2 that we may expand c,d(k) k /2mc,d, where mc and the following Hamiltonian ≈ md are the effective masses of the uncoupled particles. X Now, by changing the ratio mc/md, the detuning δ and H = (k)ˆc† cˆk + [ δ + (k)]dˆ† dˆk + (16) Rabi frequency Ω , one can adjust the regions in mo- eff c k − d k R k mentum space where the lossy state d interacts with the long-lived state c, thereby tuning the loss-profile of the ˆ ˆ +ΩR(ˆck† dk + dk† cˆk) , (17) polaritons. In the actual numerical calculations, discussed below, where c,d(k) is the dispersion of the exciton and the we employ the same set of Dyson’s equations introduced photon, respectively, and δ the relative energy offset. ΩR in sectionII where the impurity dispersion and loss rate is the Rabi frequency of interconversion between exciton are the one of the c-polariton (the one showing a loss rate and photon. The latter is subject to losses at a rate γd, which is larger at small momenta), which corresponds to 10
3 c 2 (a) (b) (c) (d) (e) A 8 B 1 6
/E 0 ! 4 1 2 2 3 0 0 0.5 1.0 1.5 2 0 0.5 1.0 1.5 2 0 0.5 1.0 1.5 2 0 0.5 1.0 1.5 2 0 0.5 1.0 1.5 2 k/qB k/qB k/qB k/qB k/qB B 44 4 4 4 4 B B B B B /E
E 3 E 3 E 3 E 3 E 3
)] 3 (f) (g) (h) (i) (j) )]/ )]/ )]/ )]/ )]/ k k k k k k ( ( ( ( ( c c c c c ( 2 2 2 2 2 [ [ [ [ [ c 2 "
2Im 1 2Im 1 2Im 1 2Im 1 2Im 1
- 1 - - - -
2Im[ 0 0 0 0 0 00.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0 0.5 1.0 1.5 2 0 0.5 1.0 1.5 2 0 0.5 1.0 1.5 2 0 0.5 1.0 1.5 2 0 0.5 1.0 1.5 2 k/qB k/qB k/qB k/qB k/qB k/qB k/qB k/qB k/qB k/qB
FIG. 7. (Upper row) The spectral function Ac(k, ω) of the impurity polariton and (lower row) its decay rate (without interactions with the fermions) for different Rabi frequencies ΩR:ΩR/EB = 5.59, 11.18, 13.69, 17.68, 25.00 [from (a) to (e) and from (f) to 2 (j)]. The Fermi energy F = kF /2mf is set by kF = 0.22qB . The polariton mixes an excitonic state c (or stable internal state in the atomic context) with a photonic mode d (or short-lived internal state in the atomic context). The short-lived state is −4 characterized by the decay rate γd = 333 EB and mass md/mc = 10 ; the fermion mass is mf = mc. The offset is set to zero, 2 δ = 0. In the panels (a)–(e), the dashed red line is the dispersion relation c(k) = k /2mc of the exciton, the dotted-dashed red line shows the lower polariton branch Re[εc(k)], which for large k approaches c(k). In the panels (f)–(j), the dotted vertical 2 red line is kF , and the dotted horizontal gray line is 4ΩR/γd.
R replacing c(k) + iγ(k)/2 with εc(k) in G0 of equation same time induces polariton losses concentrated around (A3). Here we are assuming that only one of the two zero momentum. At sufficiently large ΩR the quantum polariton branches is excited and neglecting the mixing Zeno effect creates a forbidden region and the dispersion between the branches induced by the interaction with of the polarons is significantly altered. Upon entering the electrons. We are also neglecting the fact that the the quantum-Zeno regime, both polaron-branches show polariton-electron coupling strength depends on momen- an inversion of their curvature i.e. of the polaron effec- tum. The inclusion of these features, which in general tive mass from positive to negative. In particular, for the leads to appreciable quantitative differences, is however attractive-polaron the negative curvature becomes very beyond the scope of the present work. large in the closest vicinity of the forbidden region, which In Fig.7, we show the impurity spectral function seems to result from a mutual repulsion between the lat- (k, ω) for different Rabi frequencies [see the panels (a)– ter and the polaron branch. Ac (e)] ΩR/EB = 5.59, 11.18, 13.69, 17.68, 25.00. The pho- ton (or short-lived state in the atomic context) has a large Since the quantum-Zeno Fermi-polaron branch has 4 decay rate γd = 333 EB and a small mass md/mc = 10− ; weight only for finite momenta lying outside the forbid- the detuning is set to δ = 0. In this parameter regime, the den region, the excitation of finite-momentum states is 2 maximum decay rate is approximately given by 2ΩR/γd. necessary for experimental probing. In transition metal The parameters are chosen in such a way to have a steep dichalcogenide monolayer semiconductors, optical exci- dispersion relation d(k) compared to c(k), so that the tation cannot transfer momenta which are comparable tails of the polariton loss-profile drops on a scale smaller with the Fermi momentum [39], so that a direct probe of than EB. This feature is important to have a well de- the attractive-polaron branch at large-enough momenta fined DS and thus an almost forbidden region at large loss seems unfeasible. This portion of the attractive-polaron rates. The impurity and bath fermions are of the same branch is however accessible indirectly via the interme- mass. The Fermi wavevector is kF = 0.22qB. In the pan- diate excitation of the repulsive-polaron branch at zero els (f)–(j), we show the polariton loss-profiles correspond- momentum and subsequent electron-mediated decay into ing to the panels (a)–(e). At small Rabi frequency, the the attractive branch at momenta larger than kF . For polariton dispersion shows minimal deviations from the ultracold atoms, the radiofrequency-spectroscopy [36– original exciton dispersion and has also negligible losses, 38, 41] usually employed for probing Fermi-polarons can- so that the two polaron-polariton peaks are essentially not be used here if one wants to excite finite-momentum the same as in Fig.1a. By increasing the Rabi frequency states. Instead, two-photon Raman transitions [59] can the impurity-polariton dispersion develops a pronounced be exploited for a direct excitation of the attractive- minimum due to the photon admixture. The latter at the polaron branch at a given finite momentum. 11
V. CONCLUSIONS TABLE I. Diagrammatic representation of the Green’s func- tions for the impurity (left column), fermions (middle col- In this work, we have studied the 2D Fermi-polaron umn) and molecules (right column). To each line energy and problem in presence of impurity loss and drive. Our momentum variables are attached, which are taken as the ar- approach, employing diagrammatic methods for non- gument of the corresponding Green’s function. The double equilibrium Green’s function on the Keldysh time- arrows on the molecular lines indicate their composite na- contour [50, 51], allows to consistently take into account ture. Keldysh GFs are drawn with solid lines, whereas solid- dissipation and pump and study the strongly dissipa- dashed (dashed-solid) lines represent retarded (advanced) tive limit. It provides a framework to determine non- Green functions. equilibrium distributions of fermions, molecules and im- Impurity Fermion Molecule R R R purities. We emphasize they cannot be determined on Gc Gf G phenomenological grounds basing on equilibrium quan- tum field theory and require self-consistent treatment A A A G Gf G involving the Keldysh component of the Green func- c tions. We leave the investigation of these intrinsically K K K non-equilibrium features to future work. Gc Gf G In the strongly dissipative regime, we observed the reemergence of unitary Fermi-polaron-dynamics, but with an additional length scale set by the size of the the space-time variable by x = (r, t). The correlator region forbidden by the quantum-Zeno effect. The non- in (A1a) is the retarded GF, which is non-zero only R trivial interplay between this scale and the Fermi scale for t > t0. Related to G is the advanced GF denoted induces large qualitative modifications of the impurity- A R by Gα (x, x0) = [Gα (x0, x)]∗. The correlator in (A1b) fermion scattering and consequently of the polaron dis- is the so-called Keldysh GF or statistical propagator. persion. We expect this phenomenology to be generic While the retarded GF GR contains information only for situations where the losses single-out a region of mo- about the excitation spectrum of the system (being mentum space. In particular, we have shown this to be the expectation value of a (anti)commutator for bosons naturally the case for exciton-polariton impurities and (fermions) it is indeed normalized to 1), the Keldysh also possible to engineer for atomic gases. GF also depends on the occupation of excitation modes. In the same way as the polaron concept has been in- The Keldysh GF is in general independent of GR/A strumental in developing a systematic understanding of and has to satisfy its own evolution equation. Only in fundamental phenomena in quantum many-body physics, thermal equilibrium the fluctuation-dissipation relation polarons undergoing a Quantum-Zeno projection pro- (FDR) allows to express the value of GK in terms of vide a step towards the understanding of the many-body GR and GA. It is only in this case that the latter are physics in presence of strong measurements. In partic- sufficient for the description of the problem. Away from ular, our findings identify the emergence of a new scale thermal equilibrium, the introduction of the Keldysh associated with the Zeno-forbidden region in an interact- GF is necessary and in particular in driven-dissipative ing system as a fundamental and generic ingredient able systems even in stationary states [51]. to give rise to novel collective phenomena. In our system only the impurities are subject to loss and gain, while the fermions are stable particles. Thus, since the number of impurities is very small compared to Appendix A: Non-equilibrium Green’s functions the density of fermions, we can assume that the fermions remain in a thermal equilibrium state. The Keldysh GF 1. Impurity, fermions and pairing can then be related by FDR to the remaining two GFs as follows: The non-equilibrium Green’s functions (GFs) are the eq GK (p) = F (ω)[GR(p) GA(p)], (A2) building blocks of our theory. For the details of the the- f f f − f ory in the path-integral formulation, and some calcula- eq with the thermal distribution Ff (ω) = 1 2nf (ω) where tions, we refer the reader to AppendixB. The GFs are − nf (ω) = θ(F ω) is the Fermi-Dirac distribution at defined by the relations [50]: − T = 0, θ the unit step function, and F the Fermi energy. R D E Using the time- and space-translation invariance in the iGα (x, x0) θ(t t0) αˆ(x), αˆ†(x0) , (A1a) steady state, all GFs are expressed by their Fourier trans- ≡ − ∓ R K D E forms in the relative coordinate, for example, Gf (x, 0) iGα (x, x0) αˆ(x), αˆ†(x0) , (A1b) R R P ip x R ≡ ≡ G (r, t) G (x) = e · G (p), and hereafter, fol- ± f ≡ f p f where the subscript +/ indicates the anti- lowing Ref. [50], we denote the momentum-energy vector − p = (k, ω), and the scalar product p x = k r ωt; commutator/commutator, andα ˆ† creates a particle P P R · · − in one of the degrees of freedom, that is, for a fermion the sum here is p = k dω/2π. Finally, the retarded R/A 1 αˆ = fˆ or for the impurityα ˆ =c ˆ ; hereafter, we denote (advanced) GF is given by G = (ω ε (k) i0)− , † † † † f − f ± 12
(a) f f erly introduced within the Keldysh functional-integral formalism by performing a Hubbard-Stratonovich trans- R c c c c ⌃c = + formation. This leads to the following interaction part of the functional-integral action: Z Z S = dt dr (¯cf¯∆ + ∆¯ fc), (A4) (b) ∆,f,c − f f C R where is the Keldysh time-contour [50]. The molecular ⌃ = + C c c field ∆ introduced via the Hubbard-Stratonovich trans- formation has actually no own dynamics, the latter be- FIG. 8. The diagrams showing the scattering processes that ing generated only by the interaction with impurities and give rise to the impurity (a) and the molecule (b) self-energies fermions. Formally, this translates in the bare molecular R/A 1 (after truncation of external lines). The Keldysh structure, propagator being simply given by G0,∆ = (1/U i0)− . i.e., the types of Green functions, is explicitly presented, cf. The non-interacting Keldysh component of the± molecu- Tab.I. In (a), during propagation of the impurity, the particle lar action is vanishingly small and is overwritten by the scatters with a fermion f to form a molecule ∆, and then interaction corrections as expressed by the self-energies, the molecule decays back into the fermion and the impurity. R which we introduce next. This process yields Σc . In (b), a propagating molecule decays into a fermion and the impurity, which subsequently collide to form the molecule. This process gives rise to ΣR . To highlight ∆ 2. Dyson equation and diagrammatic expansion the compound molecular structure, we denote its propagator with the line with two arrows. The interactions in the system lead to the modifica- tion (so-called dressing) of the bare GFs. We compute where the infinitesimal imaginary number, denoted with these corrections from Dyson’s equation, for details see i0, ensures the proper support in real time. ± AppendixB3. Within the Keldysh formalism, one has Neglecting the impurity-fermion interaction, it is to solve two coupled Dyson’s equations for the retarded straightforward to evaluate the impurity GFs. These and the Keldysh GFs, respectively, since, as explained GFs, which are the building blocks for our diagrammatic above, the two are in general independent when the out- approach, are referred to as “bare” (denoted by an ad- of-equilibrium case is considered. In what follows, we ditional subscript 0). Since we assume that the fermions are interested in the spectral signatures of the polarons, are not renormalized by the presence of a small number which are contained in the retarded, dressed GF of the of impurities, we drop the subscript “0” from fermionic impurity GFs. The impurity GFs read R 1 1 Gc (k, ω) = . (A5) R/A γ(k) R G0,c (k, ω) = , (A3a) ω εc(k) + i Σ (p) + i0 ω ε (k) i γ(k) − 2 − c − c ± 2 K R A Here the retarded self-energy ΣR results from the inter- G0,c(k, ω) = F0,c(k)[G0,c(p) G0,c(p)], (A3b) c − actions, the relevant process being the one where the with F0,c(k) = 1 + 2n0,c(k) and n0,c(k) = Ω(k)/γ(k). impurity scatters with a fermion, temporarily forms a The distribution function F0,c has the same structure molecule, which then decays back into an impurity and a eq as Ff (the opposite relative sign is due to the bosonic fermion. This process is displayed in Fig.8a in terms of statistics for the impurities). However, in place of the Feynman diagrams, where, as summarized in TabI, each thermal occupation the ratio of pump to loss rate ap- line corresponds to a GF and the intersection points to pears in F0,c. This ratio gives the maximal number of interaction vertices. In particular, dressed and bare GFs impurities n0,c(k) allowed to occupy a given momentum are depicted as thick and thin lines, respectively. The state k. choice of the diagrams in Fig.8 is a direct generalization In the Fermi-polaron problem, bound-state formation of a non-self-consistent T -matrix approach to a driven- between the impurity and a bath-fermion plays a promi- dissipative setting. In particular, if we assume thermal nent role [60]. Thus, featuring a pole already in the equilibrium, i.e., set the impurity losses to zero, it re- two-body problem, the particle-particle scattering vertex produces the result of the non-self-consistent T -matrix dominates the interactions in the Fermi polaron problem. approach in the impurity limit. To capture this pairing channel we introduce an auxil- Every diagram corresponds to a mathematical expres- iary molecular field, which provides access to molecular sion for the self-energy that can be obtained as follows. degrees of freedom having the bound state as a long lived To each line an energy-momentum vector is assigned. For excitation. As already known from thermal equilibrium, internal lines, the latter is taken as the argument of the this approach to the interactions between impurity and GF, while for the external lines it sets the argument of bath is fully equivalent to the T -matrix approach. As we the self-energy we are computing. At each vertex energy show in AppendixB, the molecular field can be prop- and momentum conservation is imposed. Finally, a sum 13 is performed over all energy-momentum vectors of the propagating molecule can decay into an impurity-fermion internal lines. pair, but due to Pauli blocking, the fermion cannot oc- With this procedure, we obtain from Fig.8a the fol- cupy any state from the Fermi Sea. Due to the composite lowing self-energy: nature of the molecule, its GF is affected by impurity loss. In turn, the impurity self-energy in Eq. (A8) is affected R i X K R K A by loss through ΣR . Since GK (p ) does not depend on Σ (p) = G (p0)G (p+p0)+G (p0)G (p0 p) , ∆ 0,c 0 c −2V f ∆ ∆ f − Ω(k) in the impurity limit, there is no pump-dependence p0 (A6) left in any of the retarded GFs, which are discussed in where V is the area of the two-dimensional system. The the main text. advanced/retarded GF of each particle is always multi- Finally, let us make more explicit the connection be- plied by a Keldysh GF of the other particle, which carries tween the above results and the known case of the loss- information about occupations of the states. less Fermi polaron. For this purpose, let us consider In order to compute the impurity self-energy from the quasi-particle excitations of the system, whose en- Eq. (A6), we need the molecular GF. The retarded com- ergy and damping is given by the complex poles of the ponent is given by retarded impurity GF. By setting γ(k) to zero we ob- tain the following equation determining the quasi-particle 1 dispersion: ω(k) = ε (k) ΣR(k, ω(k)). Exactly the GR (p) = , (A7) c − c ∆ 2 R same equation [30, 35, 44, 47] is derived when the vari- ν Σ∆(p) + i0 − − ational approach is employed with the Chevy polaron 2 R P ˆ ˆ with ν = 1/U and Σ∆(p) the molecular self-energy. ansatz Ψ(k) = χcˆk† FS + k0,q χk0qcˆk† k0+qfk†0 fq FS , − | i | i − | i The molecular Keldysh GF can be parametrized again where k0 / FS and q FS. This trial wave function, K R A ∈ ∈ as G = F∆(G G ). If the number of impurities is with variational parameters χ and χ 0 , includes single ∆ ∆ − ∆ k q very small, the number of molecules is also very small, particle-hole excitations and is able to capture the influ- resulting in F∆ 1. In the impurity limit and assuming ence of the molecular state on the impurity propagation. the Fermi-bath to≈ be at T = 0, we thus have 1 X 1 ΣR(k, ω)= , Appendix B: Keldysh path-integral formulation of c V ν2 ΣR (k+k , ω+ε (k )) + i0 the theory k0 FS − − ∆ 0 f 0 ∈ (A8) where the summation is restricted to momenta from in- In our work, we follow the approach presented in side the Fermi Sea (FS) for which k0 6 F . This con- Refs. [50] and [51]. The quantum average of any opera- straint on the fermion momenta is| a direct| consequence tor can be calculated using the path-integral formulation. K of the presence of Gf in Eq. (A6). Physically, during the Here, the quantum averages over a state, given by a den- propagation of the impurity a transient molecular state sity matrix, are calculated by summing over all possible can be formed if the impurity collides with a fermion, realizations of fields weighted by a complex exponent that and this can be taken only from the Fermi Sea. can be interpreted as an action. Below, we present the R The molecular self-energy Σ∆(p) is needed to evaluate structure of the theory and derive the formulas presented R R Σc (p) and in turn Gc (p). The dominant process that in the main text. leads to the modification of the propagation of a free We first consider the impurities. Although the statis- molecule is depicted in Fig.8b. A propagating molecule tics of the impurities is irrelevant, we choose them to be decays virtually into an impurity and a fermion, which bosons, since this simplifies the functional-integral for- subsequently collide to form the molecule again. The mulation. We then treat the fermionic bath using Grass- expression for the self-energy resulting from the Feynman mann fields. Finally, we include the interactions between diagram reads impurities and fermions, which can be reformulated by introducing the molecular fields and their interaction ver- R i X K R K R tex with fermions and impurities. Σ (p)= G (p0)G (p p0)+G (p0)G (p p0) , ∆ 2V 0,c f − f 0,c − p0 (A9) 1. The partition function where we use the bare impurity GFs. Since all functions R entering the equation for Σ∆ are known analytically, the self-energy can be computed explicitly in the impurity Our starting point is the partition function of the sys- limit: tem, written as an integral over field configurations:
Z ¯ ¯ 1 X 1 iS[¯c+,c+,c¯−,c−,f+,f+,f−,f−] ΣR (k, ω)= . Z = [c, f]e , (B1) ∆ γ(k k0) D V 0 ω ε (k k ) ε (k )+i − +i0 k /FS − c − 0 − f 0 2 ∈ (A10) where for simplicity in the integration measure [c, f] D ≡ Notice that the sum over momenta extends outside the [¯c+, c+, c¯ , c , f¯+, f+, f¯ , f ] we suppressed depen- Fermi Sea. This is a consequence of the fact that the denceD on position− − (or momentum)− − and time of the fields. 14
The action S = Sc + Sf + Sint can be decomposed into We emphasize that in the Keldysh formulation of the 1 three parts that describe the bosonic impurities, fermions theory the cl cl component of G0−,c is zero. This prop- and the interaction, respectively. Below, we describe erty applies also− to the full dressed GF when the inter- them separately. actions are included and the average is taken over the The bosonic impurities are governed by total action S. This is referred to as the “causal struc- ture” and originates from the detailed derivation of the Z X path-integral formalism [50]. The causal structure of the Sc= dt c¯+[i∂t εc(k)]c+ c¯ [i∂t εc(k)]c (B2) − − − − − k theory is reflected in the vanishing of q q component of the GF. − 1 i[γ(k) + Ω(k)] c+c¯ (¯c+c+ +c ¯ c ) (B3) Now, we proceed to the fermionic action given by − − − 2 − − Z 1 X Sc = dt f¯(k, t)[i∂t εf (k)]f(k, t). (B12) iΩ(k) c¯+c (¯c+c+ +c ¯ c ) , (B4) − − − − − − 2 C k where the complex-valued field c c (k, t) resides on Here, the fields f and f¯ are independent Grassmann non- the forward (backward) branch of± the≡ Keldysh± contour , commuting numbers. Next, we write explicitly the fields which goes from t = to t = (forward branch) andC on forward and backward branches of the Keldysh con- then from t = + to−∞t = (backward∞ branch). Next, tour, and perform the Keldysh rotation according to: we perform the Keldysh∞ rotation−∞ according to f f f f f = 1 ± 2 and f¯ = 2 ∓ 1 . (B13) c+ c c¯+ c¯ ± √2 ± √2 ccl/q = ± − andc ¯cl/q = ± − , (B5) √2 √2 As a result, after taking the Fourier transform, we obtain the action: where the indices cl and q stand for “classical” and X 1 “quantum”, respectively. Since we are interested in the ¯ 0 Sf = fa(p)Gf,aa− 0 (p)fa (p), (B14) stationary states, we perform the Fourier transform and p express the fields in terms of the frequency variable. In the new basis, the action takes the form: where the Keldysh indices a and a0 can be either 1 or 2. The inverse of the fermionic GF is given by X 1 S = c¯ (p)G− 0 (p)c 0 (p), (B6) c α 0,c,αα α R 1 1 K ! p 1 [Gf ]− [Gf− ] Gf− = A 1 , (B15) 0 [Gf ]− where the sum over repeated indices is implied here and 1 in the following. The object G0−,c,αα0 (p) is the bare inverse which, after inversion, leads to the following GF: of the impurity GF, which reads: GR GK G = f f , (B16) A 1 f 0 GA 1 0 [G0,c]− f G0−,c = R 1 1 K . (B7) [G ]− [G− ] 0,c 0,c and, assuming fermions of infinite lifetime, the entries R/A 1 K are given by G (p) = [ω ε (k) i0]− and G is We note that the q q term is the Keldysh component f − f ± f of the inverse of the− GF, and its inverse is not equal in Eq. (A2). The causal structure of the theory leads to [G ] 0. to the Keldysh component of the GF (see below). The f 21 ≡ components of the matrix are given by 2. Interaction and molecules R/A 1 γ(k) [G ]− = ω ε (k) i , (B8a) 0,c − c ± 2 1 K The action corresponding to the interaction in Eq. (2) [G0−,c] = i[γ(k) + Ω(k)]. (B8b) is given by The inversion of Eq. (B7) results in Eq. (A3). In this Z Z way, we have S = U dt dr c¯(x)c(x)f¯(x)f(x). (B17) int − C K R G0,c G0,c We proceed with the Hubbard-Strantonovich transforma- G0,c = A , (B9) G0,c 0 tion that introduces the molecular field which mediates the contact interaction in Eq. (B17). and the entries are given by To this end, we introduce an auxiliary fermionic field ∆(x) that decouples the interaction potential: G0,c,αα0 (p) = cα(p)¯cα0 (p) (B10) Zh i Z iSc[c,c¯] iSint iS∆+iS∆,f,c = [¯c, c] c (p)¯c 0 (p)e . (B11) e = [∆¯ , ∆]e , (B18) D α α D 15 with p1,u p3, Z Z 1 S∆ = dt dr ∆(¯ x)U − ∆(x), (B19a) CZ Z a q b S = dt dr (¯cf¯∆ + ∆¯ fc). (B19b) ∆,f,c − C Now, the action is linear in the fermionic fields, and thus the fermionic degrees of freedom can be integrated out. Consequently, we write p4, ↵ p2,v Z [f,¯ f]eiSf +iS∆,f,c = eiS∆c , (B20) D FIG. 9. The vertex generated by the action S∆c, see Eq. (B24), describing interaction between a molecule (denoted where the effective resulting action is given by with a solid line with a double arrow) and impurities (denoted Z Z with a solid line with a single arrow). Each line carries the 1 α S∆c = dx0 dx0c¯β(x0)∆u(x)γa,u (B21) momentum variable and the Keldysh index. The interaction −2 × mediated by fermions is denoted by the wavy line, which is β labeled by a momentum variable q and two Keldysh indices Gf,ab(x, x0)γb,v∆v(x0)cα(x). (B22) × (together they determine Gf,ab(q)). To the diagram we assign The sum over indices a, b, u, v runs over 1 and 2, whereas the expression that is shown in Eq. (B25). the sum over α, β runs over cl and q. The matrix γcl is the cl identity matrix in the Keldysh indices, i.e, γa,b = δa,b and q q γ has unit entries on the off-diagonal, i.e, γ = 1 δ . where q = p1 p4 = p2 p3. ab − a,b − − Finally, the molecular action S∆ leads to the inverse bare From the vertex shown in Fig.9 we can evaluate the molecular GF: self-energies. To calculate the molecular self-energy, we connect the external impurity lines, i.e., we set p3 = p4, R 1 1 K [G ] [G− ] 1 0,∆ − 0,∆ and introduce G0,c,αβ(p3) for that line. In this way, we G0−,∆ = A 1 , (B23) 0 [G0,∆]− arrive at the self-energy shown in Fig.8b. When the sum over Keldysh indices is evaluated, we obtain a sum of dia- R/A 1 where [G ]− = 1/U i0. The action in the Keldysh grams in which retarded/advanced and Keldysh GFs ap- 0,∆ ± indices and in momentum space can be written as S∆ = pear, see Fig.8b. Analogously, by joining the molecular P ¯ 1 p ∆a(p)G0−,∆,ab(p)∆b(p); here a, b takes values 1 or 2. line, i.e., by setting p1 = p2 and introducing G∆,uv(p1), The off-diagonal Keldysh component is vanishingly small we obtain the diagram from Fig.8a. Taking the sum over for the bare propagator, but it is renormalized when the the Keldysh indices we arrive at the sum of diagrams, see interactions are taken into account. Fig.8a. However, in the resummation scheme employed Going back to the momentum representation the ac- in this work, in the calculation of the impurity self-energy tion reads the dressed molecular GF is used instead of the bare one.
1 X α β S = 0 G (p p )γ γ ∆c −2V f,ab 1 − 4 u,a b,v × p1,p2,p3,p4 3. Dyson equation ∆¯ (p )∆ (p )¯c (p )c (p ), (B24) × u 1 v 2 β 3 α 4 The action of the system, after tracing out the fermions where the prime on the sum sign indicates the constraint and introducing the molecular channel consists of three p + p = p + p ; the sum over repeated discrete indices 1 2 3 4 parts a, b, u, v, α, β is implicit. In Fig.9, we display the vertex generated by the action S = Sc + S∆ + S∆c. (B26) S∆c defined by Eq. (B24). Each line carries a momen- tum and a Keldysh index. A solid line with an arrow The first two terms on the right-hand side generate the denotes the impurity field, while double arrows represent impurity and molecular bare Green functions, respec- the molecular fields. When forming diagrams, two lines tively. The third term is responsible for the conversion are connected to form a GF that, depending on the car- of a fermion and an impurity in an impurity-fermion ried Keldysh indices, is taken as the retarded/advanced molecule. or Keldysh GF, as shown in Tab.I. In the figure, the The evaluation of the dressed GFs, which are given by wavy line stands for the fermion GF. It carries an index average over the field weighted by eiS, can be understood at both ends, denoted by a and b, and a momentum q. as averaging over eiSc+iS∆ and an expansion of eiS∆c in The diagram is thus translated into the expression powers of the interaction. The usual diagrammatic for- 1 1 mulation of the theory is then recovered when we assign γα G (q)γβ , (B25) − 2 V u,a f,ab b,v lines to the GFs as shown in Tab.I. 16
As a result of the summation of the diagrams, we ob- The self-energy of the molecule is generated by the tain Dyson equations for the GFs: diagram shown in Fig.8b. Employing the vertex from the previous appendix, we arrive at G = G + G Σ G , (B27a) c 0,c 0,c · c · c i X α β G∆ = G0,∆ + G0,∆ Σ∆ G∆, (B27b) Σ (p) = γ 0 γ 0 · · ∆,ab 2V aa bb × (B27c) p0
G 0 0 (p p0)G (p0). (B30) where by we denote the multiplication of the Keldysh f,a b − 0,c,αβ · matrices. Here, the self-energies Σc and Σ∆ are matrices R The retarded self-energy is thus given by Σ∆ = [Σ∆]11, generated by irreducible diagrams for the impurity and which leads to Eq. (A9). molecular GFs, respectively. They have a causal struc- Analogously, the diagram from Fig.8 generates the ture similar to the inverse GFs [50]: self-energy for the impurity. Here, we take the dressed
R K GFs for the molecule with self-energy given by (B30). As Σ∆ Σ∆ Σ∆ = A (B28) a result, we obtain: 0 Σ∆ αβ i X β α Σ (p) = − Tr γ G (p0) γ G (p + p0) , for the molecule, and x 2V · f · · ∆ p0 A 0 Σc (B31) Σc = R K (B29) Σc Σc where the trace acts in the Keldysh indices. Finally, us- R ing Eq. (B29) we find Σx = [Σx]q,cl, from which follows for the impurity. Eq. (A6).
[1] H.-P. Breuer and F. Petruccione, B 99, 224307 (2019). The theory of open quantum systems (Oxford Uni- [13] B. Skinner, J. Ruhman, and A. Nahum, “Measurement- versity Press, Oxford, 2002). induced phase transitions in the dynamics of entangle- [2] H. M. Wiseman and G. J. Milburn, ment,” Phys. Rev. X 9, 031009 (2019). Quantum measurement and control (Cambridge Univer- [14] N. Katz, M. Ansmann, R. C. Bialczak, E. Lucero, R. Mc- sity Press, 2009). Dermott, M. Neeley, M. Steffen, E. M. Weig, A. N. Cle- [3] B. Misra and E. C. G. Sudarshan, “The Zeno’s paradox land, J. M. Martinis, and A. N. Korotkov, “Coherent in quantum theory,” J. Math. Phys. 18, 756–763 (1977). state evolution in a superconducting qubit from partial- [4] W. M. Itano, D. J. Heinzen, J. J. Bollinger, and D. J. collapse measurement,” Science 312, 1498 (2006). Wineland, “Quantum Zeno effect,” Phys. Rev. A 41, [15] R. Vijay, C. Macklin, D. H. Slichter, S. J. Weber, K. W. 2295 (1990). Murch, R. Naik, A. N. Korotkov, and I. Siddiqi, “Stabi- [5] C. Presilla, R. Onofrio, and U. Tambini, “Measurement lizing Rabi oscillations in a superconducting qubit using quantum mechanics and experiments on quantum Zeno quantum feedback,” Nature 490, 77 (2012). effect,” Ann. Phys. 248, 95 (1996). [16] P. Campagne-Ibarcq, P. Six, L. Bretheau, A. Sarlette, [6] P. Facchi and S. Pascazio, “Quantum Zeno subspaces,” M. Mirrahimi, P. Rouchon, and B. Huard, “Observing Phys. Rev. Lett. 89, 080401 (2002). quantum state diffusion by heterodyne detection of fluo- [7] T. J. Elliott, W. Kozlowski, S. F. Caballero-Benitez, and rescence,” Phys. Rev. X 6, 011002 (2016). I. B. Mekhov, “Multipartite entangled spatial modes of [17] N. Syassen, D. M. Bauer, M. Lettner, T. Volz, D. Di- ultracold atoms generated and controlled by quantum etze, J. J. Garcia-Ripoll, J. I. Cirac, G. Rempe, and measurement,” Phys. Rev. Lett. 114, 113604 (2015). S. D¨urr,“Strong dissipation inhibits losses and induces [8] G. Mazzucchi, W. Kozlowski, S. F. Caballero-Benitez, correlations in cold molecular gases,” Science 320, 1329 T. J. Elliott, and I. B. Mekhov, “Quantum measurement- (2008). induced dynamics of many-body ultracold bosonic and [18] S. D¨urr,J. J. Garc´ıa-Ripoll, N. Syassen, D. M. Bauer, fermionic systems in optical lattices,” Phys. Rev. A 93, M. Lettner, J. I. Cirac, and G. Rempe, “Lieb-Liniger 023632 (2016). model of a dissipation-induced Tonks-Girardeau gas,” [9] S. Dhar and S. Dasgupta, “Measurement-induced phase Phys. Rev. A 79, 023614 (2009). transition in a quantum spin system,” Phys. Rev. A 93, [19] J. J. Garc´ıa-Ripoll, S. D¨urr,N. Syassen, D. M. Bauer, 050103 (2016). M. Lettner, G. Rempe, and J. I. Cirac, “Dissipation- [10] X. Cao, A. Tilloy, and A. De Luca, “Entanglement in induced hard-core boson gas in an optical lattice,” New a fermion chain under continuous monitoring,” SciPost J. Phys. 11, 013053 (2009). Phys. 7, 24 (2019). [20] G. Barontini, R. Labouvie, F. Stubenrauch, A. Vogler, [11] Y. Li, X. Chen, and M. P. A. Fisher, “Quantum Zeno ef- V. Guarrera, and H. Ott, “Controlling the dynamics fect and the many-body entanglement transition,” Phys. of an open many-body quantum system with localized Rev. B 98, 205136 (2018). dissipation,” Phys. Rev. Lett. 110, 035302 (2013). [12] A. Chan, R. M. Nandkishore, M. Pretko, and G. Smith, [21] R. Labouvie, B. Santra, S. Heun, and H. Ott, “Bistabil- “Unitary-projective entanglement dynamics,” Phys. Rev. ity in a driven-dissipative superfluid,” Phys. Rev. Lett. 17
116, 235302 (2016). 021011 (2020). [22] B. Yan, S. A. Moses, B. Gadway, J. P. Covey, K. R. A. [41] Y. Zhang, W. Ong, I. Arakelyan, and J. E. Thomas, Hazzard, A. M. Rey, D. S. Jin, and J. Ye, “Observation of “Polaron-to-polaron transitions in the radio-frequency dipolar spin-exchange interactions with lattice-confined spectrum of a quasi-two-dimensional Fermi gas,” Phys. polar molecules,” Nature 501, 521 (2013). Rev. Lett. 108, 235302 (2012). [23] B. Zhu, B. Gadway, M. Foss-Feig, J. Schachenmayer, [42] F. Tonielli, R. Fazio, S. Diehl, and J. Marino, “Orthog- M. L. Wall, K. R. A. Hazzard, B. Yan, S. A. Moses, onality catastrophe in dissipative quantum many-body J. P. Covey, D. S. Jin, J. Ye, M. Holland, and A. M. systems,” Phys. Rev. Lett. 122, 040604 (2019). Rey, “Suppressing the loss of ultracold molecules via the [43] W. Berdanier, J. Marino, and E. Altman, “Universal continuous quantum Zeno effect,” Phys. Rev. Lett. 112, dynamics of stochastically driven quantum impurities,” 070404 (2014). Phys. Rev. Lett. 123, 230604 (2019). [24] H. Fr¨oml, A. Chiocchetta, C. Kollath, and S. Diehl, [44] R. Combescot, A. Recati, C. Lobo, and F. Chevy, “Nor- “Fluctuation-induced quantum Zeno effect,” Phys. Rev. mal state of highly polarized Fermi gases: simple many- Lett. 122, 040402 (2019). body approaches,” Phys. Rev. Lett. 98, 180402 (2007). [25] L. D. Landau, “Uber¨ die Bewegung der Elektronen in [45] S. Z¨ollner,G. M. Bruun, and C. J. Pethick, “Polarons Kristallgitter,” Phys. Z. Sowjetunion 3, 644 (1933). and molecules in a two-dimensional Fermi gas,” Phys. [26] L. D. Landau and S. I. Pekar, “Effective mass of a po- Rev. A 83, 021603 (2011). laron,” Zh. Eksp. Teor. Fiz 18, 419 (1948). [46] M. M. Parish, “Polaron-molecule transitions in a two- [27] H. Fr¨ohlich, H. Pelzer, and S. Zienau, “Properties of slow dimensional Fermi gas,” Phys. Rev. A 83, 051603(R) electrons in polar materials,” The London, Edinburgh, (2011). and Dublin Philos. Mag. J. Sci. 41, 221 (1950). [47] P. Kroiss and L. Pollet, “Diagrammatic Monte Carlo [28] R. P. Feynman, “Slow electrons in a polar crystal,” Phys. study of quasi-two-dimensional Fermi polarons,” Phys. Rev. 97, 660 (1955). Rev. B 90, 104510 (2014). [29] A. Rosch and T. Kopp, “Heavy particle in a d- [48] V. Ngampruetikorn, J. Levinsen, and M. M. Parish, “Re- dimensional fermionic bath: A strong coupling ap- pulsive polarons in two-dimensional Fermi gases,” EPL proach,” Phys. Rev. Lett. 75, 1988 (1995). 98, 30005 (2012). [30] F. Chevy, “Universal phase diagram of a strongly in- [49] R. Schmidt, T. Enss, V. Pietil¨a, and E. Demler, “Fermi teracting Fermi gas with unbalanced spin populations,” polarons in two dimensions,” Phys. Rev. A 85, 021602 Phys. Rev. A 74, 063628 (2006). (2012). [31] R. Combescot and S. Giraud, “Normal state of highly [50] A. Kamenev, Field theory of non-equilibrium systems polarized Fermi gases: full many-body treatment,” Phys. (Cambridge University Press, 2011). Rev. Lett. 101, 050404 (2008). [51] L. M. Sieberer, M. Buchhold, and S. Diehl, “Keldysh [32] R. A. Duine and A. H. MacDonald, “Itinerant ferromag- field theory for driven open quantum systems,” Rep. netism in an ultracold atom Fermi gas,” Phys. Rev. Lett. Prog. Phys. 79, 096001 (2016). 95, 230403 (2005). [52] M. Randeria, J.-M. Duan, and L.-Y. Shieh, “Bound [33] X. Cui and H. Zhai, “Stability of a fully magnetized ferro- states, Cooper pairing, and Bose condensation in two di- magnetic state in repulsively interacting ultracold Fermi mensions,” Phys. Rev. Lett. 62, 981 (1989). gases,” Phys. Rev. A 81, 041602 (2010). [53] J. R. Engelbrecht and M. Randeria, “New collective mode [34] V. Pietil¨a,D. Pekker, Y. Nishida, and E. Demler, “Pair- and corrections to Fermi-liquid theory in two dimen- ing instabilities in quasi-two-dimensional Fermi gases,” sions,” Phys. Rev. Lett. 65, 1032 (1990). Phys. Rev. A 85, 023621 (2012). [54] J. R. Engelbrecht and M. Randeria, “Low-density repul- [35] P. Massignan, M. Zaccanti, and G. M. Bruun, “Polarons, sive Fermi gas in two dimensions: Bound-pair excita- dressed molecules and itinerant ferromagnetism in ultra- tions and Fermi-liquid behavior,” Phys. Rev. B 45, 12419 cold Fermi gases,” Rep. Prog. Phys. 77, 034401 (2014). (1992). [36] A. Schirotzek, C.-H. Wu, A. Sommer, and M. W. Zwier- [55] D. K. Efimkin and A. H. MacDonald, “Many-body theory lein, “Observation of Fermi polarons in a tunable Fermi of trion absorption features in two-dimensional semicon- liquid of ultracold atoms,” Phys. Rev. Lett. 102, 230402 ductors,” Phys. Rev. B 95, 035417 (2017). (2009). [56] M. Randeria, J.-M. Duan, and L.-Y. Shieh, “Supercon- [37] C. Kohstall, M. Zaccanti, M. Jag, A. Trenkwalder, ductivity in a two-dimensional Fermi gas: Evolution from P. Massignan, G. M. Bruun, F. Schreck, and R. Grimm, Cooper pairing to Bose condensation,” Phys. Rev. B 41, “Metastability and coherence of repulsive polarons in a 327 (1990). strongly interacting Fermi mixture,” Nature 485, 615 [57] E. Braaten, H.-W. Hammer, and G. P. Lepage, “Open (2012). effective field theories from highly inelastic reactions,” [38] M. Koschorreck, D. Pertot, E. Vogt, B. Fr¨ohlich, M. Feld, (2016), arXiv:1612.08047 [hep-ph]. and M. K¨ohl,“Attractive and repulsive Fermi polarons [58] P. Nozieres and S. Schmitt-Rink, “Bose condensation in in two dimensions,” Nature 485, 619 (2012). an attractive fermion gas: From weak to strong coupling [39] M. Sidler, P. Back, O. Cotlet, A. Srivastava, T. Fink, superconductivity,” J. Low Temp. Phys. 59, 195 (1985). M. Kroner, E. Demler, and A. Imamoglu, “Fermi [59] C. Shkedrov, G. Ness, Y. Florshaim, and Y. Sagi, “In polaron-polaritons in charge-tunable atomically thin situ momentum-distribution measurement of a quantum semiconductors,” Nat. Phys. 13, 255 (2017). degenerate Fermi gas using Raman spectroscopy,” Phys. [40] L. B. Tan, O. Cotlet, A. Bergschneider, R. Schmidt, Rev. A 101, 013609 (2020). P. Back, Y. Shimazaki, M. Kroner, and A. Imamoglu, [60] N. Prokof’ev and B. Svistunov, “Fermi-polaron problem: “Interacting polaron-polaritons,” Phys. Rev. X 10, Diagrammatic Monte Carlo method for divergent sign- alternating series,” Phys. Rev. B 77, 020408 (2008).