Quantum phase transitions in a resonant-level model with dissipation: Renormalization-group studies

Quantum phase transitions in a resonant-level model with dissipation:
Renormalization-group studies

Chung-Hou Chung Electrophysics Department, National Chiao-Tung University, Hsinchu, Taiwan, R.O.C.    Matthew T. Glossop Department of Physics, University of Florida, Gainesville, FL 32611-8440, USA    Lars Fritz Institut für Theoretische Physik, Universität Köln, Zülpicher Straße 77, 50937 Köln, Germany Department of Physics, Harvard University, Cambridge MA 02138, USA    Marijana Kirćan Max-Planck-Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany    Kevin Ingersent Department of Physics, University of Florida, Gainesville, FL 32611-8440, USA    Matthias Vojta Institut für Theoretische Physik, Universität Köln, Zülpicher Straße 77, 50937 Köln, Germany
July 9, 2019
Abstract

We study a spinless level that hybridizes with a fermionic band and is also coupled via its charge to a dissipative bosonic bath. We consider the general case of a power-law hybridization function with , and a bosonic bath spectral function with . For and , this Bose-Fermi quantum impurity model features a continuous zero-temperature transition between a delocalized phase, with tunneling between the impurity level and the band, and a localized phase, in which dissipation suppresses tunneling in the low-energy limit. The phase diagram and the critical behavior of the model are elucidated using perturbative and numerical renormalization-group techniques, between which there is excellent agreement in the appropriate regimes. For this model’s critical properties coincide with those of the spin-boson and Ising Bose-Fermi Kondo models, as expected from bosonization.

pacs:
75.20.Hr,74.72.-h

I Introduction

Quantum phase transitionssubirbook () in mesoscopic systems form a growing area of condensed matter research. From a theoretical perspective, it is known that models of a finite system (the “impurity”) coupled to infinite baths may exhibit boundary quantum phase transitions (QPTs), at which only a subset of the degrees of freedom becomes critical.mvreview () Such models help to advance our understanding of quantum criticality in strongly correlated systems: Concepts and solution techniques developed in the impurity context may be applied to lattice models, e.g., within the framework of dynamical mean-field theory (DMFT)dmft-rmp () and its extensions. This approach has been followed in connection with the “local criticality” proposed to underlie the anomalous non-Fermi-liquid behavior of several heavy-fermion systems.lcqpt () On the experimental side, QPTs in mesoscopic few-level systems are of great interest, both for the unprecedented opportunity to probe quantum criticality in a direct and highly controlled fashion,potok (); dias () and for their numerous potential technological applications, e.g., in nanoelectronics and quantum information processing.kondo (); device (); noise ()

In recent years, QPTs have been identified and studied in a number of quantum impurity models.mvreview () Such models can contain both fermionic bands (e.g., conduction-electron quasiparticles) and bosonic baths (e.g., phonons, spin fluctuations, or electromagnetic noise). Analytical and numerical techniques have been refined to analyze the critical behavior of these models. Analytical approaches based on bosonization or conformal field theory have been used extensively, although their applicability is limited, e.g., to certain forms of the bath spectrum. For other situations, powerful epsilon-expansion techniques have been developed. As such expansions are asymptotic in character, a comparison with numerical results is mandatory to assess their reliability.

An example with especially rich behavior is the fermionic pseudogap Kondo model,withoff () which features QPTs between Kondo-screened and local-moment ground states.withoff (); cassa (); GBI (); insi (); lars (); larslong () Essentially perfect agreement between the results of various epsilon expansions (around different critical dimensions) and numerical renormalization-group (NRG) calculations has been found in critical exponents as well as universal amplitudes such as the residual impurity entropy.lars (); larslong ()

Impurity models that include bosons are harder to tackle numerically than pure-fermionic problems due to the large Hilbert space, and fewer results are available. The development of a bosonic versionBTV (); BLTV () of Wilson’s NRG approachBulla:07 () has made possible a detailed nonperturbative study of the spin-boson model, where tunneling in a two-state system competes with dissipation.leggett () For the case of Ohmic dissipation, the spin-boson model has long been known to display a QPT of the Kosterlitz-Thouless type. In the sub-Ohmic case, the model instead exhibits a line of continuous QPTs governed by interacting quantum critical points (QCPs).BTV (); BLTV (); VTB () (The latter lie in a different universality class than the QCP of the pseudogap Kondo model.)

Of particular interest, both for mesoscopics and in the context of extended DMFT for correlated lattice-systems,edmft (); chitra () are impurity models with fermionic and bosonic baths. The best-studied member of this class is the Bose-Fermi Kondo model,bfk (); sengupta (); bfknew (); kircan2 (); kirchner () with a spin- local moment coupled to fermionic quasiparticles (the regular Kondo model) as well as to a bosonic bath. The latter may describe spin or charge fluctuations of the bulk system in which the impurity is embedded. The scope of NRG applications has recently been widened to provide a comprehensive treatment of an Ising-symmetric version of the Bose-Fermi Kondo model.Glossop:05 (); Glossop:07 ()

The purpose of this paper is to investigate a somewhat simpler quantum impurity model containing both fermionic and bosonic baths, namely a resonant-level model of spinless electrons, with the impurity charge coupled to a dissipative reservoir. In standard notation, its Hamiltonian is

(1)

with characterizing the hybridization between conduction electrons of energy and the impurity level at energy , and coupling bosons of energy to the impurity occupancy. Without loss of generality, and are taken to be real and non-negative. Equation (I) represents perhaps the simplest nontrivial Bose-Fermi quantum impurity model, making it a paradigm for this class and an ideal problem for detailed comparison between analytical and numerical results.

The model is completely specified by the impurity level energy , the hybridization function

(2)
and the bosonic bath spectral function
(3)

with and acting as fermionic and bosonic cutoffs, respectively. Thus, in addition to a power-law spectrum for the bosonic bath density of states (DOS) characterized by an exponent , we consider a nonconstant particle-hole (p-h) symmetric hybridization function characterized by an exponent . Increasing (and hence depleting the hybridization function around the Fermi level ) and increasing both act to suppress tunneling between the local level and the conduction band. For most of the numerical work presented in Sec. III, we fix , , and the hybridization strength , then tune the dissipation strength to the vicinity of a QPT.

Although the bath densities of states and , do not require separate specification, it will facilitate comparison between numerical and perturbative results to assume that , for all , . In this case, and , with the fermionic and bosonic DOS given, respectively, by

(4)
(5)

where and are normalization factors. Thus, and . The metallic case is recovered for , and Ohmic dissipation corresponds to taking .

It is convenient to identify a pseudospin—making clear the close relationship between model (I) and the spin-boson model and its variants—by writing

(6)

In the model described by Eq. (I), the friction caused by the bosonic bath competes with the resonant tunneling of electrons. In contrast to the simpler spin-boson model,leggett () the tunneling properties are determined by the hybridization function .

Figure 1: Schematic phase diagram of the dissipative resonant-level model (I), in the parameter space spanned by exponents and characterizing the low-energy behavior of fermionic and bosonic baths, respectively. (A finite coupling to both baths is assumed.) For , the model shows a boundary quantum phase transition (as the couplings and are varied) between a delocalized phase and a localized phase. (The physics along the line is identical to that of the spin-boson model.) In contrast, for and the system is generically delocalized, whereas it is localized in the rest of the parameter regime (shaded). The impurity entropy is discussed in the text. Two perturbative RG expansions are employed: around the free-impurity fixed point, where the expansion is controlled about (Sec. II.2) and around the resonant-level fixed point, where the expansion is controlled in (Sec. II.3).

For the model features a Z symmetry of particle-hole type [assuming as noted above], namely , , and . Then, we expect that the competition between resonant tunneling and dissipation yields a QPT between a “delocalized” phase (), in which the principal effect of dissipation is to renormalize the tunneling amplitude, and a “localized” phase () with a doubly degenerate ground state, where the tunneling amplitude renormalizes to zero in the low-energy limit. We note that for the case of a metallic fermionic bath [ in Eq. (4)], bosonization techniques can be used to map the model (I) to the spin-boson model.hur1 () (The same applies to the Ising-symmetric Bose-Fermi Kondo model with , and this equivalence has been verified using NRG.Glossop:05 (); Glossop:07 ())

In this paper, we employ renormalization-group (RG) techniques to map out the phase diagram of the Hamiltonian (I) and to establish over what range of bath exponents and the model can be tuned to a delocalized-to-localized QPT, akin to that of the spin-boson model. We do so using both perturbative RG methods, based on epsilon-expansion techniques developed in the context of the pseudogap Kondo and Anderson models,larslong () and the Bose-Fermi extensionGlossop:05 (); Glossop:07 () of the NRG approach, which allows us to access the entire parameter range of the model.

Our main result is summarized in Fig. 1, which illustrates the qualitative behavior of the model in the plane spanned by the bath exponents and . A delocalized-to-localized transition—which for is identically that of the spin-boson model—is present at as well. A more detailed discussion is given in Sec. II.3.3.

The remainder of the paper is organized as follows. The perturbative RG analysis is outlined in Sec. II, where results for various critical exponents are obtained by expansion around two distinct fixed points. In Sec. III, we provide nonperturbative NRG results for the model, including discussion of the phase diagram, the response to a local field, and the single-particle spectral function. We find excellent quantitative agreement between analytical and numerical results in the appropriate limits. Although the critical properties of the model (I) for are established via the mapping to the spin-boson model, we confirm the equivalence by direct calculation.

Ii Perturbative renormalization group

ii.1 Zero-temperature phases

We begin by discussing the trivial fixed points of the model (I) in the presence of p-h symmetry, . As a characterization, we will refer to the residual impurity entropy , which is defined as the impurity contribution to the total entropy in the limit temperature .mvreview ()

For , the impurity is decoupled from both baths. We denote this free-impurity fixed point by FImp. The ground state is doubly degenerate: .

For and one has a resonant-level model with a power-law conduction-band DOS given by Eq. (4). The hybridization is relevant in the RG sense (w.r.t. FImp) for , and hence the impurity charge strongly fluctuates.GBI (); larslong () We refer to this as the delocalized fixed point (Deloc), which, as discussed in Ref. larslong, , is located at intermediate RG coupling, . Somewhat surprisingly, the impurity entropy is , and vanishes only in the metallic case . For , by contrast, the hybridization is RG-irrelevant, and the delocalized fixed point merges with FImp.spinfoot ()

The dissipative coupling turns out to be RG-relevant at the FImp fixed point for (see, e.g., Refs. leggett, and BTV, ). It tends to suppress tunneling in the low-energy limit. By analogy with the spin-boson model, this can be expected to result in a doubly degenerate ground state, , i.e., a phase with broken symmetry. This localized fixed point (Loc) corresponds to coupling values . (Note that for the effect of the bosonic bath is weak, not causing localization.)

The preceding discussion suggests that, for and , a QPT separates a delocalized (small-dissipation) phase from a localized (large-dissipation) phase. Clearly, this applies only to the case of p-h symmetry, . Otherwise the symmetry of the Hamiltonian is broken from the outset, and the phase transition upon variation of the dissipation strength will be smeared into a crossover; this is analogous to the behavior of the spin-boson model in the presence of a finite bias. Furthermore, in situations where the system is localized at , there will be a first-order transition upon tuning from positive to negative values (as in an ordered magnet subject to a field).

We now proceed with an RG treatment of the model (I), carried out without recourse to bosonization. We can access quantum-critical properties via two distinct expansions: (i) an expansion around the free-impurity fixed point (Sec. II.2), which is formally valid provided that the couplings to both baths are small, and (ii) an expansion around the resonant-level fixed point (Sec. II.3), performed after exactly integrating out the fermions. The second approach proves to have the wider range of applicability.

ii.2 RG expansion around the free-impurity limit

In this subsection, we apply an RG epsilon expansion for weak couplings near the free-impurity fixed point where .

ii.2.1 RG equations

We model the bosonic bath by a relativistic scalar field, , in dimensions, with the action

(7)

being a momentum-space cutoff (related to the energy cutoff via with being a velocity). This produces a DOS of the form

(8)

for , with . [Note that is just a symmetrized version of defined in Eq. (5).] Similarly, we represent the fermionic bath by Dirac fermions in dimensions:

(9)

with and being the (Fermi) velocity, which reproduces the DOS defined in Eq. (4). A path-integral representation of Eq. (I) reads

(10)

Power counting yields the bare scaling dimensions of fields and couplings with respect to : , , , , and . Thus, we can carry out an RG expansion around and , where both and become marginal, defining

(11)

In order to proceed with the RG analysis, we define a renormalized field and couplings and according to

(12)

where is an arbitrary renormalization energy scale and , , and are renormalization factors. As is usual for impurity problems, there is no renormalization of the bosonic and fermionic bulk propagators, since the impurity only provides a one-over-volume correction to the bulk properties. The relevant diagrams for obtaining the one-loop RG beta functions are shown in Fig. 2.

Figure 2: Diagrams appearing in the perturbative expansion for the dissipative resonant-level model. Dashed, solid, and wiggly lines denote respectively , , and propagators. The gray (black) circles are the interaction vertices (). (a) and (b): fermion self-energy diagrams to one-loop order. (c) One-loop vertex renormalization of .

Following standard procedures,bgz () the one-loop RG beta functions of the dissipative resonant-level model are given by

(13)

where the calculation parallels that of Ref. larslong, . The corresponding factors, to one-loop accuracy, are , , and .

The RG flows arising from Eqs. (13) are plotted in Fig. 3. In this subsection, we consider the case ; the regime is discussed in Sec. II.3.2. Fixed points at and describe the delocalized (Deloc) and localized (Loc) phases, respectively. For , where

(14)

both these fixed points are stable: For small and large , the ground state is delocalized, characterized by strong local charge fluctuations due to resonant tunneling between the impurity and the conduction electron bath (). In the opposite limit of small and large , we find a localized ground state where charge tunneling renormalizes to zero in the low-energy limit (). An unstable critical fixed point [Cr], located at , controls the QPT between these two phases. This critical fixed point lies on the separatrix specifying the phase boundary in the - plane between the delocalized and localized phases.

As approaches from below, the critical fixed point merges with the delocalized fixed point (which itself merges with FImp as from below). Hence, no transition occurs for : Deloc and FImp are unstable w.r.t. infinitesimal bosonic coupling, such that the ground state is always localized for .

Figure 3: Schematic RG flow diagrams for the dissipative resonant-level model with p-h symmetry. Although these diagrams are obtained by expansion about the free-impurity fixed point and hence are formally valid as , they are confirmed by expansion about the delocalized fixed point (Sec. II.3) and NRG calculations (Sec. III) to capture the correct physics for all and . The horizontal axis denotes the renormalized bosonic coupling ; the vertical axis denotes the renormalized hybridization . (a) : Stable fixed points at Deloc and Loc describe the delocalized and localized phases, respectively. The continuous impurity QPT is controlled by the critical fixed point Cr. (b) : The delocalized (Deloc) fixed point is unstable against finite . As , the Deloc fixed point merges with the free-impurity fixed point (FImp). (c) : is irrelevant. In both (b) and (c), the flow is toward Loc for any finite .

ii.2.2 Correlation-length exponent

In the following, we discuss the properties of the boundary QPT, controlled by the critical fixed point Cr. We start with the correlation-length exponent , describing the flow away from criticality: The characteristic energy scale above which quantum-critical behavior is observed vanishes assubirbook ()

(15)

where is a dimensionless measure of the distance to criticality, defined such that () corresponds to the localized (delocalized) phase. Upon linearization of the RG beta functions around the Cr fixed point, we obtain

(16)

Clearly, diverges as and together. By expanding the square-root in Eq. (16), the inverse correlation length exponent can be approximated as . The same result, valid for small , is also obtained in Sec. II.3 following an RG expansion valid near the strong-coupling fixed point. The divergence of as is demonstrated numerically in Sec. III.2.1 and the form compared to Eq. (16).

ii.2.3 Response to a local field

The local impurity susceptibility is the impurity response to a field applied only to the impurity.mvreview () Here, for the spinless resonant-level model under consideration, the level energy plays the role of a local electric field. Defining the impurity “magnetization” , with the pseudospin as specified in Eq. (6), it follows that

(17)

is nothing other than the impurity capacitance.

Near criticality, is expected to follow a power-law form

(18)

up to a nonuniversal cutoff scale . This relation defines the anomalous exponent , which governs the anomalous decay of the impurity “spin-spin” correlation function and is calculated via

(19)

The renormalization factor obeys the exact relationmvreview (); bfknew ()

(20)

which is graphically represented in Fig. 4(a). This allows us to derive the exact result

(21)

at the Cr fixed point, a relation that is borne out by the numerical results presented in Sec. III.

ii.2.4 Conduction electron -matrix

The conduction electron -matrix, describing the scattering of the electrons off the impurity, is another important observable, being central to the calculation of transport properties. For a resonant-level model, the -matrix is given by where is the full impurity (-electron) Green’s function, graphically represented in Fig. 4(b). As with the local susceptibility, we expect a power-law behavior of the -matrix spectral density near criticality:

(22)

It has been shownlarslong () that all critical fixed points for in the pseudogap Anderson and Kondo models display as , which behavior has been observed in a number of separate studies.Bulla:97 (); Logan:00 (); Bulla:00 ()

Using the exact relation , we can derive an exact result for the critical point of the dissipative resonant level model:

(23)

Thus, even though the multiplicative prefactor of the behavior (22) is expected to exhibit both and dependence, the power law followed at criticality is identical to that of the pseudogap Kondo and Anderson models.

ii.2.5 Hyperscaling and other critical exponents

The QCP is expected to satisfy hyperscaling relations characteristic of an interacting fixed point, including scaling in dynamical quantities.mvreview () It follows that the correlation-length exponent and the anomalous exponent are sufficient to determine all critical exponents associated with the application of a local field.insi (); mvreview () For example, one can define exponents and through the limit of the local susceptibility near criticality:

(24)

One can also determine critical exponents and associated with the local magnetization :

(25)

Thus, near criticality

(26)

and

(27)

where, in contrast to Eqs. (21) and (23), the higher-order corrections do not cancel. Section III reports NRG results for several of these critical exponents that demonstrably obey the hyperscaling relations.

Figure 4: (a) Exact relation for the local susceptibility. The black triangle denotes the full vertex function and the dashed double line denotes the full impurity level propagator. (b) The large dot denotes the full hybridization vertex.

ii.3 RG expansion around the delocalized fixed point

In addition to the RG expansion for and , as described in Sec. II.2, a second epsilon expansion can be performed around the Deloc fixed point.

ii.3.1 RG equations

To begin, we integrate out the conduction electrons, which is an exact operation for the present model. The resulting action islarslong ()

(28)

where the local fermions are now “dressed” by the conduction lines,

(29)

is a nonuniversal energy scale, and . For , the term dominates the propagator at low energies. Then, dimensional analysis of the bosonic coupling (here w.r.t. the Deloc fixed point) yields

(30)

which implies that an RG expansion can be controlled in the smallness of

(31)

We introduce a dimensionless coupling according

(32)

and, following the procedure described in Sec. II.2, we find that the only contribution to is that shown in Fig. 2(c), which reads (note that to this order)

(33)

The RG beta function for is

(34)

It is clear from Eq. (34) that for and , there exists a critical fixed point at

(35)

which controls the delocalized-to-localized transition. The RG flow diagram is sketched in Fig. 5.

Note that the critical coupling approaches zero as and/or as , suggesting that beyond these limiting cases the delocalized fixed point is unstable towards the localized fixed point. The same instability has already been deduced for [i.e., for ], based on expansion about the free-impurity fixed point (see Sec. II.2). The behavior for is analyzed in the next section.

Figure 5: RG flow diagram of the dissipative resonant-level model near the delocalized (Deloc) fixed point for , . The two stable phases are governed by the delocalized () and localized () fixed points, separated by the critical fixed point [ specified in Eq. (35)].

ii.3.2 The regime

For , the perturbation theory described in Sec. II.3 is singular due to the divergent DOS in the bosonic propagator. In this range of , the delocalized fixed point is always unstable against any infinitesimal bosonic coupling , which favors the localized fixed point.

We can gain a better understanding of this instability by considering the local bosonic propagator in the presence of the impurity. Including impurity effects via the boson self-energy, the local boson propagator is given by

(36)

where is a momentum cutoff energy scale. Let us discuss first. For , the local boson propagator is massive, meaning that the ground state for the bulk is just the empty state. For , by contrast, the local boson propagator has “negative mass”, as a consequence of which the local boson condenses at zero temperature with an expectation value . This drives the system to the localized phase where the pseudospin operator also assumes a nonzero expectation value. This reasoning supports the existence of a QPT for , with criticality reached at . For , the local boson propagator always has a negative mass, i.e., the impurity is localized. (Technically, the impurity induces a bound state in .) The observation that the ground state is always localized for is consistent with previous studies of the spin-boson modelBTV (); VTB () and the Bose-Fermi Kondo model,Glossop:05 (); Glossop:07 () which belong to the same universality class as the dissipative resonant-level model in the metallic limit .

ii.3.3 Phase diagram

The RG flow allows us to deduce that the qualitative phase diagram of the dissipative resonant-level model in the parameter space specified by and is as shown in Fig. 1. The solid line denotes the locus of points satisfying . In the unshaded region to the left of the line [i.e., for , or equivalently with defined in Eq. (14)], the RG expansion predicts a continuous QPT as and are varied. For (shaded area), the ground state of the model is always localized for any finite bosonic coupling . This is consistent with the RG flow diagrams presented in Fig. 3, where the RG expansion is carried out for . The phase diagram is confirmed by NRG results in Sec. III.

ii.3.4 Critical exponents

By linearizing the RG equation around the fixed point, the correlation-length exponent at the critical point is found to satisfy

(37)

For the anomalous exponent associated with the local susceptibility [Eq. (18)], we again have the exact property Eq. (20) [see also Fig. 4(a)], from which it follows that

(38)

The exponents and can be obtained from the hyperscaling relations (25):

(39)

and

(40)

The exponent , associated with conduction-electron -matrix, is also found to obey [see Eq. (23)]. Of course, all critical exponents for the two RG expansions (one for and one for ) are expected to be compatible since the expansions describe the same QPT. In the limit , the square root of Eq. (16) may be expanded to yield Eq. (37). The equivalences of Eqs. (26) and (39) for and of Eqs. (27) and (40) for are also readily verified.

Iii Numerical renormalization group

The NRG methodBulla:07 () has recently been extended to provide nonperturbative results for the Bose-Fermi Kondo model.Glossop:05 (); Glossop:07 () In the following, we implement the same approach for the spinless resonant-level model (I), which also involves both fermionic and bosonic baths.

Figure 6: (Color online) (a) The lowest NRG eigenstates vs even iteration number for , hybridization strength , and a range of dissipation strengths . The flows are typical of those for with . The levels at the critical coupling are shown as bold dotted lines while those nearby in the delocalized () [localized ()] phase are shown as solid [dashed] lines. As approaches in either phase, the levels follow those of the unstable critical fixed point down to progressively lower temperatures, before crossing over to the levels characteristic of the delocalized or localized stable fixed point. (b) NRG level flows for . In this case, and more generally for with , the flow is towards the localized fixed point for any , but follows the delocalized fixed point down to progressively lower temperatures as is reduced towards zero. The solid lines show the flow for .

There are three essential features of the NRG: (i) The energy axis is logarithmically discretized, introducing a discretization parameter . (ii) The Hamiltonian is then mapped to a chain form, with the impurity degrees of freedom coupled to the first site only of one or more tight-binding chains. (iii) Owing to the discretization, the tight-binding coefficients decay exponentially with increasing chain length. This allows the problem to be solved in an iterative fashion, diagonalizing progressively longer finite-length chains and thereby including exponentially smaller energy scales, , at each iterative step , , , . The RG transformation relating the effective Hamiltonians at consecutive iterations eventually reaches a scale-invariant fixed point that determines the low-temperature properties of the system.

In all applications of the NRG, the maximum number of many-body eigenstates retained from iteration to form basis states for iteration must be truncated for sufficiently large due to the limitations of finite computational power. The presence of one or more bosonic chains introduces additional considerations. First, the bosonic Hilbert space must be truncated even at iteration , allowing a maximum of bosons per site of a bosonic chain. Second, for problems involving both fermionic and bosonic chains, the fact that the bosonic tight-binding coefficients decay as the square of those for fermionic chains must be reflected in the specific iterative scheme employed. That is, only (bosonic and fermionic) excitations of the same energy scale should be considered at the same iterative step. Thus, while the fermionic chain is extended at each iteration, the bosonic chain is extended only at every second iteration. These issues, together with further details of the implementation of the Bose-Fermi NRG, are discussed in detail in Ref. Glossop:07, .

The NRG method has provided a comprehensive numerical account of the quantum-critical properties of a number of impurity problems, e.g., the fermionic pseudogap Kondo and Anderson models, the spin-boson model, and the Bose-Fermi Kondo model. In all cases it is found that the critical properties (such as exponents) are insensitive to the discretization parameter and converge rapidly with the number of retained states . For models involving bosonic baths, critical exponents also rapidly converge with increasing bosonic truncation parameter . In the following we take , with all data suitably converged for the choice and . For convenience we set .

iii.1 Phase diagram

Figure 6 shows the flow of the lowest NRG eigenstates of the effective Hamiltonian at even iteration numbers for two representative cases for : (a) and (b) . Figure 6(a) shows data obtained for and . Here, and for any , the flow is schematized by Fig. 3(a), which follows from the perturbative analysis. For , the NRG flow is towards the delocalized fixed point, where the spectrum coincides with that for coupling to the bosonic bath. For the NRG flow is towards the localized fixed point, where the spectrum coincides with that for coupling to the fermionic band. For close to , as considered in Fig. 6(a), the flow in either case is first towards the critical spectrum. The departure from the critical flow, at a crossover scale that vanishes at , is governed by the correlation-length exponent discussed in Sec. III.2.1.

Figure 6(b) shows NRG level flows for and . These flows are typical of those for any and correspond to the perturbative RG flows of Fig. 3(b). The localized ground state obtains for any . As is reduced towards zero, the levels follow those of the delocalized fixed point (obtained for ) down to progressively lower energy scales.

Figure 7: (Color online) Phase diagram in the - plane, obtained using NRG for the fixed bosonic bath exponent and the hybridization strength shown in the legend. For , we find a continuous QPT between delocalized (Deloc) and localized (Loc) phases. The critical dissipation strength is found to vanish continuously at . For only the localized phase can be accessed for . The inset shows the vanishing of with decreasing in each case, compared to the results obtained from the perturbative analysis.

Figure 7 shows the phase diagram of the model on the - plane for three different combinations of the bosonic bath exponent and the hybridization strength . For all and pairs considered, the phase-boundary value of decreases monotonically with increasing from that found for a metallic conduction band (). This is particularly clear from the data set obtained for and (circles in Fig. 7), where the metallic system undergoes a continuous QPT at a critical . With increasing , and hence growing depletion of the conduction electron density of states around the Fermi level, the critical dissipation strength required to localize the system is reduced, as expected on physical grounds. is found to vanish continuously at , with as defined in Eq. (14). This vanishing is illustrated in the inset to Fig. 7, which shows vs on a logarithmic scale.

For , localized solutions are found for arbitrarily small dissipation strength . The symbols at the largest () in each case, which lie at , mark the point at and above which no delocalized solutions can be found with . Thus, we find that we can tune the system to a QPT if, and only if, and , in complete agreement with the scenario deduced via the perturbative analyses and illustrated in Fig. 1.

For and , we find a line of Kosterlitz-Thouless-like transitions between delocalized and localized ground states, and for only the delocalized phase is accessed (provided ). For and , the essential physics is controlled by the free-impurity fixed point, regardless of the couplings and .

Figure 8: Critical dissipation strength vs hybridization strength for the pairs specified in the legend. We find that , with .

For a given pair that exhibits a continuous QPT, the critical dissipation strength varies with the hybridization strength as

(41)

provided that all scales are small compared to the cutoffs. This result, which follows from dimensional arguments [Eq. (41) can readily be obtained using Eq. (13)] and is confirmed numerically in Fig. 8, identifies as the tunneling amplitude analogous to of the spin-boson model, whereBTV () the critical dissipation strength is . A similar result for the Bose-Fermi Kondo model finds , with the bare Kondo temperature serving as a tunneling amplitude between impurity spin states.Glossop:05 (); Glossop:07 ()

It is interesting to compare the location of the phase boundary obtained using NRG with that inferred from analytical expansion. We have in mind fixing the hybridization strength and the bosonic-bath exponent (as in Fig. 7), and finding the critical coupling as a function of the conduction-band exponent . However, an analysis of the expansion around the free-impurity fixed point (Sec. II.2) reveals no simple analytical expression for the phase boundary, due to the fact that the problem is described by a two-parameter flow, which cannot be linearized in general. We have therefore analyzed the coupled differential flow equations numerically. The phase boundary can be obtained by determining the eigenvalues and eigenvectors of the linearized RG equations near the critical point and then following the RG flow backwards along the separatrix.

The inset of Fig. 7 compares phase boundaries determined via NRG (symbols) with those obtained via the perturbative RG equations (13) (dashed lines). For the range of considered by NRG, appears to vanish as a power law, with an exponent that depends on both the bosonic bath exponent and the hybridization . This apparent power law does not reflect the asymptotic behavior, revealed by the perturbative calculations to be as . (This regime is inaccessible to NRG because the merging of the critical and delocalized fixed points with decreasing make it impossible to reliably determine the critical coupling .) Nevertheless, we find the level of agreement remarkable and stress that there is no fitting procedure involved in making this comparison.

From the expansion around the delocalized fixed point (Sec. II.3), where we have a one-parameter flow, it seems possible to obtain an analytical expression for the phase boundary. However, we have to keep in mind that the dressed propagator in Eq. (II.3.1) contains terms with different frequency dependencies, and is dominated by in the low-energy limit only. (The coefficient is nonzero in general, except right at the Deloc fixed point.) The interplay of the and terms introduces a nonuniversal crossover scale into the problem, and a proper treatment including elevated energies would require a multistage RG scheme, which is beyond the scope of this paper.

iii.2 Critical exponents

iii.2.1 Correlation-length exponent

The correlation-length exponent defined in Eq. (15) is readily extracted from the crossover scale in the NRG level flows between the unstable and either of the stable fixed points. Here, denotes the NRG iteration number at which crossover is observed in a chosen NRG eigenvalue . (See Refs. Glossop:05, and Glossop:07, for further details.) Figure 9 shows vs for the pairs specfied in the legend. The dashed lines are linear fits to the log-log data, which yield the correlation length exponent , independent of the hybridization strength and the phase (Deloc or Loc) from which the QCP is accessed.

Figure 9: Crossover scale vs for the pairs specified in the legend. In the vicinity of the transition (), . The correlation-length exponent is independent both of the hybridization and of the phase from which the QCP is approached.

Figure 10: (Color online) (a) Correlation-length exponent vs conduction-band exponent for two values of the bosonic-bath exponent , as shown in the legend. The symbols show NRG data, while the dashed lines are the corresponding perturbative results [Eq. (16)], expanding about the free-impurity fixed point. We find that vanishes at , in keeping with the qualitatively distinct behavior for . (b) The same data plotted vs . For small , , a result [Eq. (37)] (shown as a dotted line) obtained by a perturbative expansion about the delocalized fixed point.

The dependence of the correlation-length exponent is demonstrated in Fig. 10(a) for two values of the bosonic bath exponent . As anticipated, for we find that within our estimated numerical error of about 1%, is in essentially exact agreement with for the spin-boson modelBTV (); VTB () (and the Ising-symmetry Bose-Fermi Kondo model, demonstrated in Ref. Glossop:05, to share the same universality class). By increasing we find that diverges as from below, i.e., as . The dashed lines are the corresponding perturbative results [Eq. (16)], with which there is excellent agreement for approaching . Figure 10(b) shows the same data plotted vs . With decreasing , the curves approach the result (shown as a dotted line), as obtained in Sec. II.C.3 by an expansion about the delocalized fixed point.

iii.2.2 Response to a local field

As discussed in Sec. II.2.3, the response to a field applied only at the impurity provides a useful probe of the locally critical properties of the model. The inset to Fig. 11(a) shows vs for and hybridization strength . Behaving as a suitable order parameter for the problem, is finite in the localized phase (), saturating to for and vanishing continuously as . In the delocalized phase (), . The main part of Fig. 11 shows vs on a logarithmic scale, from which the power-law behavior Eq. (25) is clearly apparent. The exponent is found to be . At the QCP (), the dependence of on the field defines the exponent according to Eq. (25). We typically observe such power-law behavior over several orders of magnitude of , as shown in Fig. 11(b). For , .

We note that for and , undergoes a jump at the critical point . Here, the essential behavior has been discussed in Refs. hur1, ; Glossop:05, ; Glossop:07, and Borda:05, for the case relevant to charge fluctuations on a metallic island subject to electromagnetic noise.

Figure 11: (Color online) Critical exponents and , defined in Eq. (25), for and , where . (a) Continuous vanishing of order parameter vs as with characteristic exponent (extracted as the limiting slope of the data on a logarithmic scale). The inset shows the data on an absolute scale. (b) Variation of with local level energy at the critical point . The data clearly follow a power law for small , defining the exponent .

Figure 12: Static local susceptibilty vs for , , and (circles), (stars), and (squares). The anomalous exponent in the quantum-critical regime is found to be , independent of . See text for further discussion.

We calculate the static local susceptibility via

(42)

In the delocalized phase , vanishes linearly with and thus for . In the localized phase , is nonzero as with . In the quantum-critical regime , diverges as a power law with an anomalous exponent defined in Eq. (18). For all pairs considered (such that and a critical fixed point exists), we find that

(43)

independent of . The behavior described above is clearly illustrated in Fig. 12, which shows three data sets for : one at the critical coupling and one close to it in either phase. In this example, we extract .

iii.2.3 Hyperscaling

As discussed in Sec. II.2.5, critical exponents for the present model are expected to obey hyperscaling relations derived via a scaling ansatz for the critical part of the free energy that assumes the critical fixed point is interacting.insi () This expectation is borne out by the numerical analysis: we find hyperscaling relations to be obeyed to within the estimated error (typically less than 1%) across the range of displaying critical behavior. For example, for the case , and . Thus, the values and extracted from the data presented in Fig. 11 obey Eqs. (25) to within numerical uncertainty.

iii.3 Spectral function

We now turn to the single-particle spectral function , calculated via

(44)

where is a many-body eigenstate of NRG iteration , and is the partition function; for the p-h symmetric parameters studied. The discrete delta-functions are Gaussian broadened on a logarithmic scale: a standard NRG procedure discussed, e.g., in Ref. Bulla:07, . We set the broadening parameter such that for the simplest resonant-level model (with , , and ) is in optimal agreement with the exact result .

Figure 13: (Color online) (a) Spectral function vs for , , , and three values of the dissipation strength: (dotted line), (solid line), and (thick dashed line). At the quantum critical point , , which behavior is also followed for close to and . In the delocalized phase , there is a crossover in to the behavior Eq. (45) for . For the data shown, . (b) The crossover behavior is more readily seen in the modified spectral function , which shows the ultimate low- behavior throughout the delocalized phase , for , and throughout the localized phase . In order of decreasing crossover scale, delocalized-phase spectra are shown for (dotted line) and for , , , and (solid lines); localized-phase spectra (dashed lines) are shown for and , , , and . The critical spectrum is shown as a thick dashed line.

Figure 13(a) shows vs