Anisotropic Kondo Model Research Articles

Variational principle for quantum impurity systems in and out of equilibrium: Application to Kondo problems

We provide a detailed formulation of the recently proposed variational approach [Y. Ashida et al., Phys. Rev. Lett. 121, 026805 (2018)] to study ground-state properties and out-of-equilibrium dynamics for generic quantum spin-impurity systems. Motivated by the original ideas by Tomonaga, Lee, Low, and Pines, we construct a canonical transformation that completely decouples the impurity from the bath degrees of freedom. By combining this transformation with a Gaussian ansatz for the fermionic bath, we obtain a family of variational many-body states that can efficiently encode the strong entanglement between the impurity and fermions of the bath. We give a detailed derivation of equations of motions in the imaginary- and real-time evolutions on the variational manifold. We benchmark our approach by applying it to investigate ground-state and dynamical properties of the anisotropic Kondo model and compare results with those obtained using matrix-product state (MPS) ansatz. We show that our approach can achieve an accuracy comparable to MPS-based methods with several orders of magnitude fewer variational parameters than the corresponding MPS ansatz. Comparisons to the Yosida ansatz and the exact solution from the Bethe ansatz are also discussed. We use our approach to investigate the two-lead Kondo model and analyze its long-time spatiotemporal behavior and the conductance behavior at finite bias and magnetic fields. The obtained results are consistent with the previous findings in the Anderson model and the exact solutions at the Toulouse point.

Exploring the anisotropic Kondo model in and out of equilibrium with alkaline-earth atoms

We propose a scheme to realize the Kondo model with tunable anisotropy using alkaline-earth atoms in an optical lattice. The new feature of our setup is Floquet engineering of interactions using time-dependent Zeeman shifts, that can be realized either using state-dependent optical Stark shifts or magnetic fields. The properties of the resulting Kondo model strongly depend on the anisotropy of the ferromagnetic interactions. In particular, easy-plane couplings give rise to Kondo singlet formation even though microscopic interactions are all ferromagnetic. We discuss both equilibrium and dynamical properties of the system that can be measured with ultracold atoms, including the impurity spin susceptibility, the impurity spin relaxation rate, as well as the equilibrium and dynamical spin correlations between the impurity and the ferromagnetic bath atoms. We analyze the non-equilibrium time evolution of the system using a variational non-Gaussian approach, which allows us to explore coherent dynamics over both short and long timescales, as set by the bandwidth and the Kondo singlet formation, respectively. In the quench-type experiments, when the Kondo interaction is suddenly switched on, we find that real-time dynamics shows crossovers reminiscent of poor man's renormalization group flow used to describe equilibrium systems. For bare easy-plane ferromagnetic couplings, this allows us to follow the formation of the Kondo screening cloud as the dynamics crosses over from ferromagnetic to antiferromagnetic behavior. On the other side of the phase diagram, our scheme makes it possible to measure quantum corrections to the well-known Korringa law describing the temperature dependence of the impurity spin relaxation rate. Theoretical results discussed in our paper can be measured using currently available experimental techniques.

Universal spatial correlations in the anisotropic Kondo screening cloud: Analytical insights and numerically exact results from a coherent state expansion

We analyze the spatial correlations in the spin density of an electron gas in the vicinity of a Kondo impurity. Our analysis extends to the spin-anisotropic regime, which was not investigated in the literature. We use an original and numerically exact method, based on a systematic coherent-state expansion of the ground state of the underlying spin-boson Hamiltonian, which we apply to the computation of observables that are specific to the fermionic Kondo model. We also present an important technical improvement to the method, that obviates the need to discretize modes of the Fermi sea, and allows one to tackle the problem in the thermodynamic limit. One can thus obtain excellent spatial resolution over arbitrary length scales, for a relatively low computational cost, a feature that gives the method an advantage over popular techniques such as NRG and DMRG. We find that the anisotropic Kondo model shows rich universal scaling behavior in the spatial structure of the entanglement cloud. First, SU(2) spin-symmetry is dynamically restored in a finite domain in parameter space in vicinity of the isotropic line, as expected from poor man's scaling. We are also able to obtain in closed analytical form a set of different, yet universal, scaling curves for strong exchange asymmetry, which are parametrized by the longitudinal exchange coupling. Deep inside the cloud, i.e. for distances smaller than the Kondo length, the correlation between the electron spin density and the impurity spin oscillates between ferromagnetic and antiferromagnetic values at the scale of the Fermi wavelength, an effect that is drastically enhanced at strongly anisotropic couplings. Our results also provide further numerical checks and alternative analytical approximations for the recently computed Kondo overlaps [PRL 114, 080601 (2015)].

Nonequilibrium Dynamics of a Periodically-driven Impurity Spin Coupled to the Bath of Ultracold Fermions by Kondo Coupling

The dynamics of the anisotropic Kondo model under the periodic driving is studied. The periodic driving we consider is the local magnetic field which couples with the impurity spin and the in-plane exchange interaction. At the special point in the parameter space (Toulouse limit), the anisotropic Kondo model can be transformed to the non-interacting resonant level model and thereby the exact results for arbitrary external fields can be obtained. It is shown that characteristic behavior of the dynamics changes with the field intensity, the field frequency, and the Kondo temperature. When the local magnetic field intensity is very large compared with the field frequency and the Kondo temperature, the time dependence of the impurity spin and the spin density of the fermionic bath shows an oscillation whose period is inversely proportional to the field intensity. We also find that another time scale, which is given by the inverse of the Kondo temperature, governs the response of the system to the external field.

Quantum criticality in a Kondo quantum dot coupled to helical edge states of interacting 2D topological insulators

We investigate theoretically the quantum phase transition (QPT) between the one-channel Kondo (1CK) and two-channel Kondo (2CK) fixed points in a quantum dot coupled to helical edge states of interacting 2D topological insulators (2DTI) with Luttinger parameter . The model was studied by Law et al (2010 Phys. Rev. B 81 041305(R)), and was mapped onto an anisotropic two-channel Kondo model via bosonization. For , the strong coupling 2CK fixed point was argued to be stable for infinitesimally weak tunnelings between the dot and the 2DTI based on a simple scaling dimensional analysis (Law et al 2010 Phys. Rev. B 81 041305(R). We re-examine this model beyond the bare scaling dimension analysis via a one-loop renormalization group (RG) approach combined with bosonization and re-fermionization techniques near weak-coupling and strong-coupling (2CK) fixed points. We find for a fixed value of that the 2CK fixed point can be unstable towards the 1CK fixed point and the system is expected to undergo a quantum phase transition between 1CK and 2CK fixed points with changing Kondo couplings. Our RG approach is controlled near K = 1. In general, this QPT can also occur upon tuning the Luttinger parameter K to a critical value smaller than unity () for fixed Konodo couplings. The QPT in our model comes as a result of the combined Kondo and the helical Luttinger physics in 2DTI, and it serves as the first example of the 1CK-2CK QPT that is accessible by the controlled RG approach. We extract quantum critical and crossover behaviors from various thermodynamical quantities near the transition. Our results are robust against particle-hole asymmetry for .

Anisotropic charge Kondo effect in a triple quantum dot.

We predict that an anisotropic charge Kondo effect appears in a triple quantum dot, when the system has twofold degenerate ground states of (1,1,0) and (0,0,1) charge configurations. Using bosonization and refermionization methods, we find that at low temperature the system has the two different phases of massive charge fluctuations between the two charge configurations and vanishing fluctuations, which are equivalent with the Kondo-screened and ferromagnetic phases of the anisotropic Kondo model, respectively. The phase transition is identifiable by electron conductance measurement, offering the possibility of experimentally exploring the anisotropic Kondo model. Our charge Kondo effect has a similar origin to that in a negative-U Anderson impurity.

Scaling Theory vs Exact Numerical Results for Spinless Resonant Level Model

The continuous-time quantum Monte Carlo method is applied to the interacting resonant level model (IRLM) using double expansion with respect to Coulomb interaction Ufc and hybridization V. Thermodynamics of the IRLM without spin is equivalent to the anisotropic Kondo model in the low-energy limit. Exact dynamics and thermodynamics of the IRLM are derived numerically for a wide range of Ufc with a given value of V. For negative Ufc, excellent agreement including a quantum critical point is found with a simple scaling formula that deals with V in the lowest-order, and Ufc up to infinite order. As Ufc becomes positive and large, lower order scaling results deviate from exact numerical results. Possible relevance of the results is discussed to certain Samarium compounds with unusual heavy-fermion behavior.

Kondo effect in monolayer and bilayer graphene: Physical realizations of the multichannel Kondo models

We perform a general group-theoretical study of the Kondo problem in monolayer and bilayer graphene around the charge neutrality point. Utilizing the group representation theory, we derive from symmetry considerations a family of the Kondo models for all symmetric placements with either 3- or 6-fold rotational axis of an impurity atom in an arbitrary orbital state. We find six possible classes of the partially anisotropic four-channel Kondo model. As the key result, we argue several possibilities to realize the regime of the dominant channel-symmetric two-channel Kondo effect, protected by the local symmetry and specifics of the graphene band structure. Our findings open prospects for the observation of the rich multi-channel Kondo physics in graphene and the associated non-Fermi-liquid behavior.

Entanglement in quantum impurity problems is nonperturbative

We study the entanglement entropy of a region of length 2L with the remainder of an infinite one dimensional gapless quantum system in the case where the region is centered on a quantum impurity. The coupling to this impurity is not scale invariant, and the physics involves a crossover between weak and strong coupling regimes. While the impurity contribution to the entanglement has been computed numerically in the past, little is known analytically about it, since in particular the methods of conformal invariance cannot be applied because of the presence of a crossover length. We show in this paper that the small coupling expansion of the entanglement entropy in this problem is quite generally plagued by strong infrared divergences, implying a non-perturbative dependence on the coupling. The large coupling expansion turns out to be better behaved, thanks to powerful results from the boundary CFT formulation and, in some cases, the underlying integrability of the problem. However, it is clear that this expansion does not capture well the crossover physics. In the integrable case -- which includes problems such as an XXZ chain with a modified link, the interacting resonant level model or the anisotropic Kondo model -- a non perturbative approach is in principle possible using form-factors. We adapt in this paper the ideas of [1,2] to the gapless case and show that, in the rather simple case of the resonant level model, and after some additional renormalizations, the form factors approach yields remarkably accurate results for the entanglement all the way from short to large distances. This is confirmed by detailed comparison with numerical simulations. Both our form factor and numerical results are compatible with a non-perturbative form at short distance.

Numerical renormalization group calculation of impurity internal energy and specific heat of quantum impurity models

We introduce a method to obtain the specific heat of quantum impurity models via a direct calculation of the impurity internal energy requiring only the evaluation of local quantities within a single numerical renormalization group (NRG) calculation for the total system. For the Anderson impurity model, we show that the impurity internal energy can be expressed as a sum of purely local static correlation functions and a term that involves also the impurity Green function. The temperature dependence of the latter can be neglected in many cases, thereby allowing the impurity specific heat, $C_{\rm imp}$, to be calculated accurately from local static correlation functions; specifically via $C_{\rm imp}=\frac{\partial E_{\rm ionic}}{\partial T} + 1/2\frac{\partial E_{\rm hyb}}{\partial T}$, where $E_{\rm ionic}$ and $E_{\rm hyb}$ are the energies of the (embedded) impurity and the hybridization energy, respectively. The term involving the Green function can also be evaluated in cases where its temperature dependence is non-negligible, adding an extra term to $C_{\rm imp}$. For the non-degenerate Anderson impurity model, we show by comparison with exact Bethe ansatz calculations that the results recover accurately both the Kondo induced peak in the specific heat at low temperatures as well as the high temperature peak due to the resonant level. The approach applies to multiorbital and multichannel Anderson impurity models with arbitrary local Coulomb interactions. An application to the Ohmic two state system and the anisotropic Kondo model is also given, with comparisons to Bethe ansatz calculations. The new approach could also be of interest within other impurity solvers, e.g., within quantum Monte Carlo techniques.

Kondo resonance of a microwave photon

We emulate renormalization group models, such as the Spin-Boson Hamiltonian or the anisotropic Kondo model, from a quantum optics perspective by considering a superconducting device. The infra-red confinement involves photon excitations of two tunable transmission lines entangled to an artificial spin-1/2 particle or double-island charge qubit. Focusing on the propagation of microwave light, in the underdamped regime of the Spin-Boson model, we identify a many-body resonance where a photon is absorbed at the renormalized qubit frequency and reemitted forward in an elastic manner. We also show that asymptotic freedom of microwave light is reached by increasing the input signal amplitude at low temperatures which allows the disappearance of the transmission peak.

Magnetic and transport properties of CePt3Ge Kondo lattice in crystalline and sub-micron state

We have studied a ternary phase of the CePt3±xGe composition, with the maximum values of x=0.2. A series of samples was prepared by arc melting, Czochralski method and splat cooling, and characterized by scanning electron microscopy and X-ray diffraction. The arc-melted samples are well-crystalline and single-phase with slightly varying stoichiometry close to the nominal CePt3Ge, while the splat-cooled samples form an amorphous or a nanocrystalline constitution of the ∼CePt3Ge composition. Detailed measurements of magnetization, specific heat and electrical resistivity were performed in the temperature range 0.40–300K, and magnetic fields up to 9T. In addition, a delicate neutron diffraction experiment down to 70mK was done, which excluded long-range magnetic order in the CePt3Ge compound. The specific heat data were tested to anisotropic Kondo model, and the resistivity behavior has been treated by means of electron localization effects. The experiments evidenced that the CePt3Ge compound can be considered as a Kondo lattice with low-temperature magnetic correlations, and the size effect seems to have minor role in control of the Kondo scattering phenomena in this compound.

Dynamical properties of a nonequilibrium quantum dot close to a dissipative quantum phase transition

We calculate the dynamical decoherence rate and susceptibility of a nonequilibrium quantum dot close to the delocalized-to-localized quantum phase transitions. The setup concerns a resonance-level coupled to two spinless fermionic baths with a finite bias voltage and an Ohmic bosonic bath representing the dissipative environment. The system is equivalent to an anisotropic Kondo model. As the dissipation strength increases, the system at zero temperature and zero bias show quantum phase transition between a conducting delocalized phase to an insulating localized phase. Within the nonequilibrium functional Renormalization Group (FRG) approach, we address the finite bias crossover in dynamical decoherence rate and charge susceptibility close to the phase transition. We find the dynamical decoherence rate increases with increasing frequency. In the delocalized phase, it shows a singularity at frequencies equal to positive or negative bias voltage. As the system crossovers to the localized phase, the decoherence rate at low frequencies get progressively smaller and this sharp feature is gradually smeared out, leading to a single linear frequency dependence. The dynamical charge susceptibility shows a dip-to-peak crossover across the delocalized-to-localized transition. Relevance of our results to the experiments is discussed.

Quantum impurity spin in Majorana edge fermions

We show that Majorana edge modes of two-dimensional spin-triplet topological superconductors have Ising-like spin density whose direction is determined by the d-vector characterizing the spin-triplet pairing symmetry. Exchange coupling between an impurity spin (S=1/2) and Majorana edge modes is thus Ising-type. Under external magnetic field perpendicular to the Ising axis, the system can be mapped to a two-level system with Ohmic dissipation, which is equivalent to the anisotropic Kondo model. The magnetic response of the impurity spin can serve as a local experimental probe for the order parameter.

Nonequilibrium occupation number and charge susceptibility of a resonance level close to a dissipative quantum phase transition

Based on the recent paper [Phys. Rev. Lett. 102, 216803 (2009)], we study the nonequilibrium occupation number ${n}_{d}$ and charge susceptibility $\ensuremath{\chi}$ of a resonance level close to dissipative quantum phase transition of the Kosterlitz-Thouless (KT) type between a delocalized phase for weak dissipation and a localized phase for strong dissipation. The resonance level is coupled to two spinless fermionic baths with a finite bias voltage and an Ohmic bosonic bath representing the dissipative environment. The system is equivalent to an effective anisotropic Kondo model out of equilibrium. Within the nonequilibrium renormalization-group approach, we calculate nonequilibrium magnetization $M$ and spin susceptibility $\ensuremath{\chi}$ in the effective Kondo model, corresponding to $2{n}_{d}\ensuremath{-}1$ and $\ensuremath{\chi}$ of a resonance level, respectively. We demonstrate the smearing of the KT transition in the nonequilibrium magnetization $M$ as a function of the effective anisotropic Kondo couplings, in contrast to a perfect jump in $M$ at the transition in equilibrium. In the limit of large bias voltages, we find $M$ and $\ensuremath{\chi}$ at the KT transition and in the localized phase show deviations from the equilibrium Curie-law behavior. As the system gets deeper in the localized phase, both ${n}_{d}\ensuremath{-}1/2$ and $\ensuremath{\chi}$ decrease more rapidly to zero with increasing bias voltages.

Many-particle effects in adsorbed magnetic atoms with easy-axis anisotropy: the case of Fe on the CuN/Cu(100) surface

We study the effects of the exchange interaction between an adsorbed magnetic atom with easy-axis magnetic anisotropy and the conduction-band electrons from the substrate. We model the system using an anisotropic Kondo model and we compute the impurity spectral function, which is related to the differential conductance (dI/dV) spectra measured using a scanning tunneling microscope. To make contact with the known experimental results for iron atoms on the CuN/Cu(100) surface (Hirjibehedin et al 2007 Science 317 1199), we calculated the spectral functions in the presence of an external magnetic field of varying strength applied along all three spatial directions. It is possible to establish an upper bound on the coupling constant J: in the range of the magnetic fields for which the experimental results are currently known (up to 7 T), the low-energy features in the calculated spectra agree well with the measured dI/dV spectra if the exchange coupling constant J is at most half as large as that for cobalt atoms on the same surface. We show that for an even higher magnetic field (between 8 and 9 T) applied along the ‘hollow direction’, the impurity energy states cross, giving rise to a Kondo effect which takes the form of a zero-bias resonance. The coupling strength J could be determined experimentally by performing tunneling spectroscopy in this range of magnetic fields. On the technical side, the paper introduces an approach for calculating the expectation values of global spin operators and all the components of the impurity magnetic susceptibility tensor (including the out-of-diagonal ones) in numerical renormalization group (NRG) calculations with no spin symmetry. An appendix contains a density functional theory (DFT) study of the Co and Fe adsorbates on the CuN/Cu(100) surface: we compare magnetic moments, as well as orbital energies, occupancies, centers and spreads, by calculating the maximally localized Wannier orbitals of the adsorbates.

Relaxation versus Decoherence: Spin and Current Dynamics in the Anisotropic Kondo Model at Finite Bias and Magnetic Field

Using a nonequilibrium renormalization group method we study the real-time evolution of spin and current in the anisotropic Kondo model (both antiferromagnetic and ferromagnetic) at a finite magnetic field h(0) and bias voltage V. We derive analytic expressions for all times in the weak-coupling regime max{V, h(0),1/t} >> T(c) (T(c) is the strong-coupling scale). We find that all observables decay both with the spin relaxation and decoherence rates Gamma(1/2). Various V-dependent logarithmic, oscillatory, and power-law contributions are predicted. The low-energy cutoff of logarithmic terms is generically identified by the difference of transport decay rates. For small times t << max{V, h(0)}(-1), we obtain universal dynamics for spin and current.

Erratum: Real-time renormalization group in frequency space: A two-loop analysis of the nonequilibrium anisotropic Kondo model at finite magnetic field [Phys. Rev. B80, 045117 (2009)

Erratum: Real-time renormalization group in frequency space: A two-loop analysis of the nonequilibrium anisotropic Kondo model at finite magnetic field [Phys. Rev. B<b>80</b>, 045117 (2009)

Real-time renormalization group in frequency space: A two-loop analysis of the nonequilibrium anisotropic Kondo model at finite magnetic field

We apply a recently developed nonequilibrium real-time renormalization group (RG) method in frequency space to describe nonlinear quantum transport through a small fermionic quantum system coupled weakly to several reservoirs via spin and/or orbital fluctuations. Within a weak-coupling two-loop analysis, we derive analytic formulas for the nonlinear conductance and the kernel determining the time evolution of the reduced density matrix. A consistent formalism is presented how the RG flow is cut off by relaxation and dephasing rates. We apply the general formalism to the nonequilibrium anisotropic Kondo model at finite magnetic field. We consider the weak-coupling regime, where the maximum of voltage and bare magnetic field is larger than the Kondo temperature. In this regime, we calculate the nonlinear conductance, the magnetic susceptibility, the renormalized spin relaxation and dephasing rates, and the renormalized $g$ factor. All quantities are considered up to the first logarithmic correction beyond leading order at resonance. Up to a redefinition of the Kondo temperature, we confirm previous results for the conductance and the magnetic susceptibility in the isotropic case. In addition, we present a consistent calculation of the resonant line shapes, including the determination whether the spin relaxation or dephasing rate cuts off the logarithmic divergence. Furthermore, we calculate quantities characterizing the exponential decay of the time evolution of the magnetization. In contrast to the conductance, we find that the derivative of the spin relaxation (dephasing) rate with respect to the magnetic field is logarithmically enhanced (suppressed) for voltages smaller (larger) than the renormalized magnetic field, and that the logarithmic divergence is cut off by the opposite rate. The renormalized $g$ factor is predicted to show a symmetric logarithmic suppression at resonance, which is cut off by the spin relaxation rate. We propose a three-terminal setup to measure the suppression at resonance. For all quantities, we analyze also the anisotropic case and find additional nonequilibrium effects at resonance.

Quantum Criticality Perspective on the Charging of Narrow Quantum-Dot Levels

Understanding the charging of exceptionally narrow levels in quantum dots in the presence of interactions remains a challenge within mesoscopic physics. We address this fundamental question in the generic model of a narrow level capacitively coupled to a broad one. Using bosonization we show that for arbitrary capacitive coupling charging can be described by an analogy to the magnetization in the anisotropic Kondo model, featuring a low-energy crossover scale that depends in a power-law fashion on the tunneling amplitude to the level. Explicit analytical expressions for the exponent are derived and confirmed by detailed numerical and functional renormalization-group calculations.