Easy-axis ferromagnetic chain on a metallic surface
The phases and excitation spectrum of an easy-axis ferromagnetic chain of magnetic impurities built on the top of a clean metallic surface are studied. As a function of the (Kondo) coupling to the metallic surface and at low temperatures, the spin chain exhibits a quantum phase transition from an Ising ferromagnetic phase with long-range order to a paramagnetic phase where quantum fluctuations destroy the magnetic order. In the paramagnetic phase, the system consists of a chain of Kondo-singlets where the impurities are completely screened by the metallic host. In the ferromagnetic phase, the excitations above the Ising gap are damped magnons, with a finite lifetime arising due to the coupling to the substrate. We discuss the experimental consequences of our results to spin-polarized electron energy loss spectroscopy (SPEELS), and we finally analyze possible extensions to spin chains with .
pacs:75.10.Pq, 75.20.Hr, 75.40.-s
Recent progress in scanning tunneling microscopy (STM) and other surface-manipulation techniques at the atomic scale has made it possible to build low-dimensional magnetic structures with unprecedented control [1, 2]. Atomic-scale magnetic structures possess an enormous potential for applications in spintronics as well as in classical and quantum computation due to their ability to store and process (quantum) information in the atomic spins. However, in order to build reliable atomic-size memories or qubits out of such systems, it is crucial to understand how the coupling to the environment affects their properties. For instance, it is known that quantum nanomagnets coupled to electronic degrees of freedom are prone to the effects of decoherence induced by Ohmic quantum dissipation [3, 4].
In addition to their interest for information technologies, low-dimensional magnetic systems coupled to dissipative environments are also very important from the point of view of fundamental physics. This is because they allow to study the crossover from quantum to classical behavior as the coupling to the environment is increased [5, 4]. Current theoretical models that describe local quantum degrees of freedom coupled to sources of dissipation predict rich and interesting physical phenomena at low temperatures, including quantum phase transitions and quantum criticality [6, 7, 8, 9], as well as broken continuous symmetries and long-range order (LRO) in low-dimensional systems [10, 11, 12, 13, 14, 15]. The interplay between quantum fluctuations, which are induced by quantum confinement and reduced dimensionality, and quantum dissipation, caused by the interaction with the environment, is at the core of this exotic behavior. From this perspective, magnetic nanostructures built on the top of clean metals offer a unique platform to test our understanding on dissipative quantum systems. For instance, the simplest zero-dimensional (0D) case of a single magnetic impurity coupled to a metallic environment [such like a Co or Mn adatom sitting on a Au(111) surface] is a practical realization of the celebrated Kondo effect, one of the most paradigmatic phenomena in the physics of strongly correlated electrons . This phenomenon is the spin-compensation by the conduction electrons of the magnetic moment of an impurity embedded in a metallic host at temperatures , where is the so-called Kondo temperature. Eventually, at , the ground state of the system is a many-body singlet (termed a ‘Kondo singlet’), formed by a linear superposition of and the spin of a cloud of conduction electrons that is anti-ferromagnetically correlated with the impurity local moment. The original magnetic impurity is thus completely screened by the metallic environment. In the context of magnetic adatoms deposited on clean metallic surfaces, this is a well-established phenomenon, which reveals itself as a narrow Fano resonance at the Fermi energy in the STS spectra [17, 18, 19].
Unfortunately, the extension of the single Kondo impurity to the case of many interacting magnetic impurities embedded in a metallic host represents a formidable theoretical challenge. Due to the many-body nature of the problem, exact methods are only limited to simple cases (i.e., consisting of few magnetic impurities), and at the cost expensive numerical calculations . On the other hand, dynamical mean-field theory (DMFT) methods  provide reliable results in the case of bulk 3D compounds, but are much less successful when applied to low-dimensional systems.
One particular class of such low-dimensional systems, namely, one-dimensional (1D) chains of magnetic atoms or organic molecules with active magnetic atoms, are currently under intensive experimental investigation [22, 23, 24, 25, 26]. Thus, by depositing monolayer (ML) of Co on Pt(997), which is a vicinal surface containing steps, chains of magnetic atoms have been recently created . The magnetic Co atoms decorate the steps and form ferromagnetic chains, which, when probed by an external magnetic field, typically exhibit superparamagnetism. However, in the case of Co chains at temperatures K, the magneto-crystalline anisotropy results in relaxation times which are longer than the experimental typical timescale, and the response of the system becomes effectively ferromagnetic . Another recently reported system, consisting of chains of Co atoms on Cu(775), has been shown to undergo a Peierls distortion leading to dimerization . This dimerized phase has been explained by the strong correlations existing in the partially filled shells of the Co atoms, which are in the high spin configuration. Maintaining the high spin in a geometry where the overlap between the d orbitals is small and correlations are strong, favors a ferromagnetic ground state, which in turn makes the charge-density wave instability possible.
In the above-mentioned systems, the coupling to the metallic substrate was assumed to be weak, and therefore was neglected in their analysis. However, as we show below, the metallic environment can have dramatic effects for the phase diagram at low temperatures, and it cannot be neglected when studying the magnetic excitations of the chain. With these systems in mind, in this contribution we focus on the analysis of an easy-axis ferromagnetic (FM) 1D spin chain deposited on the top of a metallic surface. Interestingly, as a function of the coupling with the substrate, we predict a quantum phase transition from an Ising ferromagnetic phase with long range order to a paramagnetic phase where quantum fluctuations proliferate and destroy the ordered phase. Although modifying the coupling to the substrate in the experimental systems created so far may be hard, in the future other systems can be created where the coupling to the substrate could be tuned, uncovering even more interesting phenomena.
The article is organized as follows: In the next Section, we introduce a simple model to describe an easy-axis ferromagnetic chain on a metallic substrate. This model is studied in the framework of mean-field theory in Section 3, and the quantum phase diagram is thus obtained. The phase diagram is displayed in Fig. 4, and it is one of the most important results of this work. In Section 4 we consider the excitation spectra of the ferromagnetic phase. We show that a weak coupling with the metallic substrate can have important consequences for the magnetic excitations of the ferromagnetic chain. This properties of the excitations, and in particular, the substrate-induced damping, can be accessible by experimental probes such as spin-polarized electron energy loss spectroscopy (SPEELS). In Fig. 5, we plot the spin-response of the 1D spin chain due the substrate-modified magnon excitation spectrum, which is also an important result of this work. Finally, in Section 6, we offer our conclusions and outlook for future research directions. In the A, we give more details about the technical aspects of our calculations.
Let us consider a chain of magnetic atoms deposited on top of metallic surface (see Fig. 1). We model the magnetic interactions between the atoms are using the spin- XXZ Hamiltonian:
In the above expression, we have assumed anisotropic exchange interactions, which originate from the reduced symmetry of the surface crystal environment, as well as from spin-orbit interactions of the electrons in the metallic substrate . Assuming that the impurity spin is implies that the single-ion anisotropy reduces to a constant term. Furthermore, for , exchange terms like can be expressed in terms of . However, for impurities with such terms must be taken into account, but we shall not pursue their analysis here.
The nature of the ground state of Eq. (1) depends crucially on the ratio , where we assume and . For , Eq. (1) describes an easy-axis (i.e. Ising-like) ferromagnet (anti-ferromagnet). In this system, the spin waves (magnons) and their bound states are separated from the ground state by a gap [28, 29]. On the other hand, for , the spin chain Hamiltonian (1) exhibits a XY phase with a gapless spin-wave spectrum [30, 31]. However, as anticipated in the introduction, in what follows we will focus only on the FM Ising limit of (1), and we refer the reader to Ref.  for a complementary study in the case of easy-plane anisotropy, where the spin chain exhibits an XY phase. Nevertheless, we note that the results reported in Section 3 are also applicable to the antiferromagnetic (AFM) regime corresponding to (see Section 6). Finally, for the sake of simplicity, we have neglected longer-range magnetic exchange interactions in Eq. (1), may be present due to the Rudermann-Kittel-Kasuya-Yosida (RKKY) mechanism mediated by the metallic substrate . We note that these interactions can be easily included, as long as they do not modify the nature of the Ising FM ground state.
The coupling between the spin chain and the metallic substrate is described by an anisotropic Kondo-exchange Hamiltonian :
where () is the spin-density of conduction electrons at the position of the -th spin impurity in the surface (with the lattice parameter of the chain). Note that this model bears some similarity with the Kondo-lattice model, which is used to describe the properties of heavy-fermion systems [16, 33, 34]. In the case of impurities with , the Kondo exchange involves more channels (with channel-dependent couplings) and results in a dynamical screening of the impurity spins occurring in several stages, involving different Kondo temperatures . We shall not pursue such an analysis here.
Despite its simplicity, the model introduced above still represents a formidable theoretical challenge, and in order to make analytical progress we shall introduce some simplifying assumptions. We approximate the metallic substrate by a collection of one-dimensional metals (cf. Fig. 2) described by the Hamiltonian
where creates an electron with linear momentum and spin projection . The index indicates that the electron reservoir described by the operators is only coupled to the impurity located at and not to others with . For this reason, we dubbed this approximation the “local bath approximation” (LBA) . Within the LBA, the spin-density operator reads:
where are the Pauli matrices.
The dimensionality of the substrate is crucial for the justification of the LBA. When the magnetic impurities embedded in a high-dimensional (2D or 3D) metallic substrate are separated by a distance of the order of a few , with the Fermi wavevector of conduction electrons, the interference of Kondo screening clouds that dynamically quench the magnetic moments, becomes negligible [36, 37, 38]. Consequently, at distances the impurities can be considered as independently screened . This is an important difference with respect to the case of a strictly 1D metallic substrate, as in the case of the 1D Kondo-lattice model . In that case, the single Kondo-impurity limit is reached only at distances . Since the Kondo temperature is an exponentially small scale, in the purely 1D geometry the single-impurity regime is only reached at extremely dilute impurity-spin concentrations [36, 37, 38].
In actual experiments, although the magnetic nanostructures are built on the top of a metallic surface, is determined by the 3D bulk electron density due to a non-vanishing overlap between the bulk and surface conduction states . For instance, the Fermi wave number for Au bulk free conduction electrons is , while for Au(111) surface states is . This fact dramatically increases the range of applicability of the LBA. The physical picture provided by the LBA is also supported by the behavior of the STS Fano line shapes in experiments on magnetic Co atoms deposited on Cu and separated by distances , which are identical to the single-impurity STS line shapes .
In addition, although in the actual experimental realizations the limit might not be strictly realized, the model resulting from applying the LBA to Eqs. (1,2) can be considered as the simplest “toy model” that captures the competition between the Kondo physics, which is a local quantum critical phenomenon, and the magnetic interactions along the chain . Concerning the latter, it is worth noting that the LBA cannot reproduce the non-local RKKY exchange coupling mediated by the metal and arising at order . As mentioned before, this is not a crucial drawback, since one can always redefine the couplings in Eq. (1) to take the RKKY contribution explicitly into account. More importantly, what eventually justifies the separate treatment of RKKY and Kondo interactions is that, although they both originate in the same Hamiltonian (2), the RKKY interaction results from electronic states deep inside the Fermi sea and is a static coupling, while the Kondo effect is purely dynamical effect originated predominantly in the Fermi surface [42, 34].
3 Mean-field analysis of the phase diagram
In order to gain insight into the low-temperature phase diagram of the system, we now focus on the limit of Eq. (1), describing an isolated FM Ising spin chain. In that case, the ground state of the isolated chain is separated from the single-magnon excitation spectrum by a gap and by a gap from the two-magnon bound states [28, 29]. This means that the long-range FM order in the isolated chain is stable at low temperatures and, to a first approximation, it is safe to take in Eq. (1).
We now introduce the coupling to the substrate Eq. (2), which has important consequences for the low-temperature behavior of the the spin chain. For instance, the nature of the ferromagnetic to paramagnetic phase transition is profoundly modified. As we show in the A, in the case the model for the Ising chain coupled to the substrate can be mapped onto the 1D dissipative quantum Ising model, which is in the universality class of the classical 3D (instead of 1D) Ising model. In the following, we exploit this fact to introduce a mean-field (MF) approach, which therefore provides a good description of the phases of the spin chain in the (easy-axis) Ising limit.
We introduce the MF approximation to the Ising term in Eq. (1) by making the replacement , which leads to the MF Hamiltonian of the spin chain:
is the local Weiss field. The MF Hamiltonian of the system is therefore written as
with describing the -th conduction electron bath. In the following we drop the lattice sub-index in and in order to lighten the notation.
We note that the MF approximation together with the LBA leads to a set of independent Kondo-impurity problems in an effective magnetic field , which is a quantity that must be obtained self-consistently. The MF thermodynamical properties of the spin-chain can therefore be obtained from the free-energy for a single Kondo impurity in a magnetic field . Note that the Weiss field does not act upon the conduction electrons, it only acts on the impurity. The mean-field self-consistency condition is obtained from the relation between the impurity magnetization,
and the definition of the Weiss field Eq. (6):
To make further progress we now need to know the function for a single Kondo impurity. Introducing the Abelian bosonization method  to describe the metallic 1D chains in Fig. 2, and the unitary Emery-Kivelson transformation  [cf. Eq. (48)], the anisotropic Kondo model can be mapped onto the resonant-level model [44, 45, 46]. The resonant-level model is exactly solvable when the parameter is chosen at the so-called “Toulouse line” , defined by , where and are the density of states and lattice parameter of the bosonized 1D chains, respectively (cf. A). In combination with the renormalization group (RG) approach , which describes the flow of parameters of the anisotropic Kondo model as a function of the conduction-band cutoff parametrized by , the mapping to the resonant-level model allows to describe the physical properties of the strong coupling limit of the Kondo model. One starts with the “bare” couplings , and renormalizes up to the point where , and at that point we use the exact solution of the resonant-level model . However, we stress that can be obtained in a more general case (and with more sophistication) using any “impurity-solver” method, such as the numerical renormalization group (NRG) . The mapping to the resonant-level model yields the result [44, 45, 46]
where is the digamma function  and the broadening of the resonant level can be interpreted as .
The phase boundary between the FM phase where , and the paramagnetic phase with is obtained from (10) in the limit (cf. Fig 3). For , the impurity exhibits a characteristic (Pauli) paramagnetic response , where
where is the magnetic susceptibility of the impurity at . On the other hand, for , approaches the saturation limit of as
The above result can be also obtained using the scaling equations, and it is modified by anisotropy at the Toulouse point. Thus, a non-vanishing solution for is found from Eq. (10) in the limit
From the above considerations, the phase boundary is therefore given by the critical line
and from here we can obtain the MF phase diagram of the spin chain (See Fig. 4).
It is interesting to note that, in the FM phase, as increases from the phase boundary, the value of the Weiss field crosses over from a regime where the Kondo singlet is still robust, to the regime of large where the Kondo singlet is destroyed (cf. Fig. 3). Therefore, the above mean-field theory yields a regime where the Kondo effect and a FM LRO coexist. In the following Section, we shall return to this point again.
It is also worth mentioning that the result in Eq. (16) clearly contrasts with the case of the isolated Ising chain, for which the analytical form of Eq. (16), is very different from the expression for encountered in the study of the classical Ising model, where . Physically, in the present case, the spin-flip processes induced by the Kondo coupling to the metal Eq. (2) can be considered as quantum fluctuations that tend to destroy the LRO in the Ising chain as for (See Fig. 4).
In the following Section we characterize the excitations of the ordered FM phase.
4 Magnetic excitations in the FM phase
We next focus on the FM phase, in the regime (right-bottom region in Fig. 4), where the spectrum of the system is gapped and the effective Weiss field is large enough so as to inhibit the Kondo screening of the impurities. We shall focus in particular on single magnon excitations that can be created by probes (such like spin-polarized electron energy loss spectroscopy, SPEELS, see Section 5) that couple weakly to the ferromagnetic chains. Such probes are assumed to produce single spin-flips, and thus the calculation of the system response can be carried out within linear response theory. This approximation allows us to neglect the effect of magnon bound states in the spectrum of the FM XXZ chain. Thus, the low-lying excitations involving a single spin-flip can be described using the Holstein-Primakoff (HP) representation of the spin operators :
where are bosonic annihilation and destruction operators, obeying , which physically represent a single spin-flip. For the XXZ chain introduced in Sec. 2, we should take in the above expressions. However, below we shall work assuming that is large, which is known to provide an accurate description of the lowest energy excitations . Assuming large means that our results also apply to spin chains with , provided the Hamiltonian is of the XXZ form and we assume that the Kondo-exchange coupling with the substrate describes the strongest channel of substrate electrons that couple to each impurity along the chain. Thus, working to leading order in (i.e. in the large- approximation), we find:
which can be easily diagonalized in a running-wave basis
and leads to
is the dispersion relation for magnons, and is the lattice parameter. In addition, in Eq. (22) we have neglected terms of , which involve interactions between the magnons. From Eq. (23), we see that for (i.e. ) the excitation spectrum of the spin chain is gapped, and the spin waves exhibit a quadratic dispersion for momentum .
Let us next consider on the Kondo-exchange coupling to the metal, Eq. (2). We shall treat this coupling using the HP representation in the large- limit:
As we will see, although this Hamiltonian no longer describe the Kondo effect (we recall that we are working in the limit of large ), the coupling to the gapless degrees of freedom in the metal still can have important consequences to the spin chain. To make further progress, we now need to integrate out the electronic degrees of freedom in the conduction band. To this end, we shall rely upon the coherent-state functional integral representation of the partition function :
where is the effective Euclidean action of the spin chain
where the first term is the Euclidean action corresponding to (22)
and is the contribution arising from Hamiltonian (24), where we integrate over the fermions . Note that this can be done exactly, as is quadratic (we assume that Fermi liquid theory applies) and leads to
where we have introduced the spin correlation functions in the conduction band:
Note that these correlation functions do not depend on because all the local fermion baths are identical. Now, since these baths are assumed to be Fermi liquids, we have that at and for (where is the Fermi energy of the metallic substrate):
which follows from the existence of a (where is the excitation energy) spectrum of particle-hole excitations in the Fermi sea at low energies. Note that since and have units of (where and denote energy and time units, respectively), which stem from the prefactor of in Eq. 30 and 31; then have units . In Fourier representation, we express (32) as:
where are the Matsubara frequencies, and therefore Eq. (29) becomes
where we have dropped the last term in Eq. (29) because it amounts to a (retarded) magnon-magnon interaction. This is consistent with neglecting magnon interactions in (23). The above contribution is quadratic, and replacing it into (26) leads to the effective action
where we have defined the magnon propagator 
with a dimensionless parameter quantifying the dissipation caused by the conduction electrons on the spin chain. We next perform the analytic continuation to real frequencies, in order to obtain the retarded magnon propagator. Thus, replacing and , yields
Thus, the pole of the propagator becomes , for small , which implies that the coupling to the metal Eq. (24) induces a finite lifetime and damping on the magnon excitations. Physically, this happens due to the presence of gapless particle-hole excitations in the excitation spectrum of the metallic substrate, and it can be interpreted as Landau damping [53, 54, 55]. Interestingly, this damping does not vanish in the limit . At first sight, this may seem surprising, since in the limit, magnon excitations correspond to a uniform magnetization along the direction. Thus, the total projection of the spin along , i.e., , is a conserved quantity in the absence of coupling to the substrate. In other words, the number of magnons is a conserved quantity of Hamiltonian in (20). However, in the presence of [Eq. (24)], is no longer conserved (i.e. the magnons can “leak out” from the chain), and therefore total magnetization fluctuations also become damped out by the coupling to the substrate.
The above treatment is valid provided the coupling to the substrate is weak and can be treated perturbatively. At stronger coupling, we have to worry about the possibility that the magnetic moments of the impurities are Kondo-screened by the substrate electrons. To deal with this situation, we shall rely on the non-perturbative treatment discussed in the A. As shown there, in the strong-coupling limit where Kondo correlations are important, it is possible to map the model introduced in Section 2 to a dissipative quantum Ising model. When coarse-grained over distances , this model can be described by a dissipative field theory, whose partition function is , where the action is given by the expression :
being () is the coarse-grained magnetization along the easy axis, i.e. , where . This result suggests that the quantum critical point (QPC) separating the paramagnetic and ferromagnetic ground states (cf. Fig 4) can be described by this field theory, and therefore the critical exponents around the QCP will correspond to those of the Wilson-Fisher fixed-point of the renormalization group (see e.g. Ref. , for a survey). However, in deriving the map to the dissipative quantum Ising model, we have neglected the Heisenberg term from Eq. (1). If this term is included, we cannot rule out that the actual values of the critical exponents near the QCP turn out to be different from those of the Wilson-Fisher universality class. Addressing this question in detail require a numerical investigation of the full model, which is beyond the scope of the present work.
Nevertheless, the field theory of Eq. (39), is capable of describing the phases of the system. Focusing on the ferromagnetic phase where , which corresponds to , it can be seen that the magnon excitations are indeed described by the propagator:
where is the energy gap. This magnon propagator indicates that the magnetic excitations of the chain in the strong coupling regime are overdamped by their coupling to the metallic substrate. Physically, this makes sense as the impurity spins are screened by the metal electrons and thus their magnetic moments undergo strong damping by the collective nature of their screening clouds. This regime must be contrasted with the weak-coupling regime, where we found the damping to be weak. Thus, a physical picture emerges that allows us to distinguish the two mean-field solutions that we briefly described at the end of Section 3 (cf. Fig. 3): In the weak coupling regime, the strong Weiss field (i.e. ) due to the FM coupling of the impurities along the chain does not allow for the Kondo screening of the local moments of the atoms in the chain. As a consequence, the magnetic excitations of the chain are only weakly affected by their coupling to the substrate, and thus acquire a small line-width . On the other hand, in the strong coupling regime, the Weiss field is much weaker (i.e. ) and the magnetic moments are effectively screened by the metal electrons, leading to overdamped magnetic excitations. In the next Section, we shall study some experimental consequences of our findings.
5 Relation to SPEELS experiments
In this Section, we focus on the effects of the dissipative environment on observable quantities. The magnetic properties of low dimensional spin systems deposited on metals have been studied with a variety of experimental techniques, like X-ray magnetic circular dichroism (XMCD), used to measure directly the magnetization of FM 1D spin chains , or local probes, like spin polarized STM, which can provide information on the magnon dispersion relation . Here we focus on the spin-polarized electron energy-loss spectroscopy (SPEELS) experiment [58, 59, 60], which provides direct information of the magnon dispersion relation.
The SPEELS cross section corresponds to the fractional number of electrons which emerge onto the solid angle after being scattered by a magnetic excitation in the energy range . According to the theory of the SPEELS experiment , is related to the spin response function by
where is the retarded spin correlation function obtained by analytic continuation from with
In Fig. (5) we show the spin response function as a function of and . The curves show resonances centered at the magnon frequencies (red line in the bottom plane), which are broadened by the effect of , the dimensionless coupling to the metal.
6 Conclusions and Outlook
We have studied the phase diagram and the excitation spectrum of a magnetic chain of atoms or molecules with easy-axis ferromagnetic interactions. We focused on a simple model of magnetic impurities displaying a Kondo-exchange interaction with the substrate. In the Ising limit, which provides a good approximation to the ground state of an easy-axis ferromagnet, we obtained the phase diagram using a mean-field theory approach. We find that this system exhibits two possible phases at zero temperature: a paramagnetic phase where the impurity spins are screened by the substrate, and a gapped ferromagnetic phase whose excitations are damped by the magnetic (Kondo) exchange with the metallic substrate. Using a bosonic representation for the spins as well as various types of mathematical mappings between models, we have also investigated the excitation spectrum in the ferromagnetic phase, which may be accessible through probes like spin-polarized electron energy loss (SPEELS). Although we have focused mainly on an easy-axis ferromagnetic (FM) chain, it is worth describing how our results can be modified in the case of the anti-ferromagnetic (AFM) chain corresponding to the parameter regime of the XXZ chain described by Eq. (1). An experimental motivation to extend our results to this regime can be found, for instance, in Ref. . Relying on the local bath approximation (LBA) and in the Ising limit, it is possible to apply the transformation , which maps the AFM to onto FM Ising chain. In presence of the Kondo exchange, this transformation must be supplemented by another rotation which takes and provided that . For , the excitations are separated from the ground state by a gap . Using the Holstein-Primakoff representation , we can compute the spectrum of an antiferromagnetic XXZ chain, which yields . The results of Sections 4 and 5 carry over to the AFM () case provided we replace by . A feature of the AFM regime that is also worth noticing is that bound states of magnons are absent from the spectrum [28, 29].
Finally, let us discuss our outlook for future research directions. One such direction is the study of the phases of chains of impurities with higher (i.e. ) spin. In this case, it is also possible to apply the same mean-field theory employed in Section 3. However, in order to obtain the impurity magnetization we must in general resort to an impurity solver such as Wilson’s numerical renormalization group [20, 16]. This is specially true for , which will will undergo a several stage Kondo effect [16, 35] and possibly underscreening . Another interesting direction will be to asses the accuracy of the local bath approximation employed throughout. This remains quite a challenge, which will most likely require to rely on numerical methods, such as quantum Monte Carlo. Indeed, it is possible that this approximation only provides an rough picture of the paramagnetic phase, which must be corrected by the inclusion of the coupling between the different baths that screen the magnetic impurities individually. However, we believe that such effects will be less important in the ferromagnetic phase (specially for large ), due to the existence of a gap that separates the ferromagnetic ground state from the excitations and therefore protects the ground state from perturbations.
Appendix A Bosonization of the fermionic chains and mapping to the dissipative quantum Ising model
In this Appendix, we implement the Abelian bosonization approach to the semi-infinite fermionic 1D chains along the direction in Fig. (2) [31, 62]. This procedure is standard and has been successfully applied to describe the low-energy properties of the single Kondo-impurity. We refer the reader to the standard bibliography [31, 62, 44, 43, 63, 64]. In the present case, the LBA greatly simplifies the complexity of the problem and enables a straightforward generalization to the case of many-impurities, each one coupled to an independent fermionic bath. At low temperature, in the bosonic representation  the Hamiltonians (3) and (2) become
where are bosonic chiral fields which obey the commutation relations , and which are related to charge and spin density-fluctuations in the 1D fermionic chains through the relations and , respectively [31, 64]. In Eq. (46) is the Fermi velocity, and in Eq. (47) is the scattering phase-shift associated with the potential , the conduction electron density of states at the Fermi energy, and the lattice parameter in the fermionic chain. For simplicity we assume these parameters to be identical for all chains. We then introduce the (Emery-Kivelson) unitary transformation 
under which the bosonic field and the spin operator transform as
Upon this transformation, the total Hamiltonian of the spin-chain coupled to the metallic bath, , transforms as , with
where we have defined . Note that the local bath approximation Eq. (3) is crucial to implement bosonization along the chains, and to put these ideas on a clear mathematical framework. It is also interesting to note that in the transformed representation, the quantum dynamics of the bath [represented by the chiral field ] appears explicitly in the Heisenberg term . Physically, this means that the Heisenberg interaction is now “dressed” by the spin-density fluctuations of the electron gas. In the case of an easy-axis spin chain in Ising limit , the effect of this term is negligible and can be ignored in a first approximation.
We now exploit the freedom to choose and set . In this case, the Hamiltonian reads
This Hamiltonian corresponds to the 1D dissipative quantum Ising model (DQIM), where now the Kondo Hamiltonian is equivalent to the spin-boson model with Ohmic dissipation [63, 5, 4, 14] with related to the dissipative parameter in the context of macroscopic quantum coherence through , and with the in-plane Kondo interaction playing the role of a magnetic field along the axis . At , the 1D DQIM is known to display a paramagnetic to ferromagnetic quantum phase transition which is in the universality class of quantum dissipative systems and whose dynamical critical exponent is . Note that this is very different with the case of the 1D quantum Ising chain, which is in the universality class of the 2D classical Ising model and where . The critical properties of this model near the quantum phase transition have been studied in the context of antiferromagnetic instabilities of Fermi liquids [56, 11] using the framework of the Hertz-Moriya-Millis theory [65, 66, 67] This theory describes critical quantum fluctuations of the order parameter around the Gaussian fixed point of the theory. The predicted value of the critical dynamical exponent and has been confirmed numerically with Monte Carlo simulations [12, 13]. Physically speaking, such a dynamical exponent implies that the effective dimensionality of the 1D spin chain coupled to the metallic bath is