The Phase Diagram in Electron-Doped LaCeCuO
Superconductors are a striking example of a quantum phenomenon in which electrons move coherently over macroscopic distances without scattering. The high-temperature superconducting oxides (cuprates) are the most studied class of superconductors, composed of two-dimensional CuO planes separated by other layers which control the electron concentration in the planes. A key unresolved issue in cuprates is the relationship between superconductivity and magnetism. In this paper, we report a sharp phase boundary of static three-dimensional magnetic order in the electron-doped superconductor LaCeCuO where small changes in doping or depth from the surface switch the material from superconducting to magnetic. Using low-energy spin polarized muons, we find static magnetism disappears close to where superconductivity begins and well below the doping where dramatic changes in the transport properties are reported. These results indicate a higher degree of symmetry between the electron and hole-doped cupratets than previously thought.
The electron-doped high- superconductors are less well-understood than their hole-doped counterparts (1). Nevertheless experiments on LaCeCuO (LCCO), PrCeCuO (PCCO) and NdCeCuO (NCCO) have great significance in the field of condensed matter physics because they provide a way to investigate particle-hole symmetry in the phase diagram where superconductivity emerges from doping a highly correlated antiferromagnetic Mott insulator or charge transfer insulator (2). In both and type materials, there appears to be a quantum critical transition as a function of doping which is characterized by sharp maxima in a variety of properties near absolute zero (3). Remarkably, the superconducting phases of the hole- and electron-doped materials have the same -wave pairing symmetry (4). However, there is uncertainty about the position of the critical doping and even the number of critical points. This results from the difficulty to monitor how the antiferromagnetic phase of the parent compound evolves and changes as a function of doping, especially on the electron-doped side (5); (6); (7). This issue remains of central importance since fluctuations associated with the quantum critical point (QCP) may be the origin of strong superconducting pairing and unusual properties in the normal state (8); (9). Although there are many similarities in the generic phase diagrams of hole and electron-doped cuprates, there are also considerable differences (10). In particular, the magnetism is much more prominent on the electron-doped side and appears to overlap significantly with superconductivity. Also the superconducting dome is much smaller and narrower in electron-doped systems (2). Both systems show a pseudo-gap opening below a temperature in the underdoped region which indicates some kind of fluctuating charge or spin order (6); (11); (12); (13); (14). However it is still unclear if the character of the pseudo-gap phase is the same in both the electron- and hole-doped cuprates (2).
As mentioned previously, determining the phase diagram in electron-doped cuprates is complicated by the difficulty in determining how the magnetism evolves and eventually disappears as a function of doping. For example in NCCO, inelastic magnetic neutron-scattering results show that the long-range antiferromagnetic order disappears close to where superconductivity first appears (6). However, both Shubnikov-de Haas Oscillations (15) and ARPES measurements (16) find a transition well inside the superconducting dome where the Fermi surface reconstructs. It is unclear what happens in between these two dopings. One possibility is that neutrons detect the disappearance of the long range 3D antiferromagnetic order where superconductivity begins. However, another less well ordered phase persists up to a larger value of where the Fermi surface reconstructs, leading to sudden changes in transport properties (17); (18). LCCO has very similar properties to PCCO and shows the characteristic -linear resistivity below a critical doping of 0.17 and Fermi liquid form () above (5). Also recent angular magneto-resistance (AMR) measurements reveal evidence for the disappearance of static magnetism above an optimal doping close to where many transport properties change abruptly (19).
A major difficulty in determining how the magnetism evolves in LaCeCuO is that it can only be made in a thin film form. Therefore, the traditional experimental techniques for studying magnetism such as neutron scattering, NMR, and bulk SR are not applicable. Recently, the technique of low-energy muon spin rotation (LE-SR) has been developed (20). In LE-SR, the mean implantation depth of the muons can be controlled from a few nm to a few hundred nm (see Methods). It is now well established that LE-SR is a powerful way to investigate both the magnetic and electronic properties of quantum materials (21); (22); (23); (24).
In this paper, we report zero-field (ZF) LE-SR measurements on thin films of LaCeCuO with Ce concentrations near the antiferromagnetic-superconducting boundary. The measurements were performed at the E4 beam-line of the Swiss Muon Source (25), at the Paul Scherrer Institute, in Switzerland. Samples with Ce concentration of , 0.08, and 0.10 were studied and had thicknesses of 200 nm (see Methods). For comparison, a fourth film with a Ce concentration of 0.08 and thickness 100 nm was also studied. Monte Carlo simulations using TRIM.SP (26); (27) show that the average muon implantation depth in LaCeCuO ranges from 7 to 110 nm with corresponding straggling of 3 and 23 nm, for implantation energies from 1 to 25 keV, respectively. Both the mean depth and range straggling are close to being linear functions of the implantation energy.
ii.1 Energy dependence of static magnetism
Typical ZF-SR spectra at 5 K and energies 5 and 21 keV in all the films are shown in Fig. 1. Note that at 21 keV, corresponding to an average implantation depth of nm and range straggling nm, one can clearly see a fast relaxing signal in the sample (red points in the Fig. 1(a)). Whereas no such fast component is observable in the two higher dopings at this energy (red points in the Figs. 1(b), Figs. 1(c) and 1(d)). This demonstrates that there is bulk static magnetism in the film which is absent at the two higher dopings. The behaviour is markedly different at the lower energy of 5 keV, corresponding to nm with nm (see blue points in Fig. 1). In particular, the fast relaxing signal is present in both the and both samples, but absent in the sample with the . The magnetic fraction is slightly larger in the 200 nm thick film compared to the 100 nm film but the transition temperature for the magnetism to disappear is the same (discussed later). It is more subtle but still clear that the signal is less damped at 5 keV compared to 21 keV in the sample (red and blue points in Fig. 1(a)). In general, static magnetism is enhanced close to the surface in all four samples, with the contrast being largest in the two samples.
The observed relaxation reported here is due to quasi-static magnetic fields. This was established with longitudinal field measurements where the applied magnetic field is along the initial spin polarization direction. In particular, a longitudinal field of 100 G was enough to quench the observed relaxation in Figs. 1. This implies the internal magnetic field must be static on the time scale of muon Larmor frequency s which is typical of 3D ordering.
The energy dependence of the magnetic volume fraction and average internal field are shown in Fig. 2. The sample is fully magnetic at all energies. The internal field is about 60 G at low energy (near the surface) and decreases only slightly at the higher energy. The sample appears magnetic below 10 keV, where the mean implantation depth is 40 nm. However, the internal field drops abruptly to zero above 10 keV. The results on the 100 nm thick film with (not shown here) are very similar to the 200 nm thick film indicating these affects are not dependent on the film thickness. In the sample there is only a very weak magnetic relaxation which is enhanced slightly close to the surface (see Fig. 1(d)).
ii.2 Temperature dependence of static magnetism
The maximum magnetic ordering temperature in the region probed by the muons is obtained from the temperature dependence of ZF-SR spectra, examples of which are shown in Fig. 3. This shows how the spectrum evolves in the sample going through the magnetic ordering transition. Below the transition there is a clear evidence for an over-damped oscillation due to quasi-static electronic moments which give rise to a broad distribution of the local internal fields. As one approaches the transition, the magnitude of this average field drops to zero. Above the transition, in the paramagnetic state, the electronic moments are rapidly fluctuating, so the observed weak time dependence of the muon polarization is due mostly to weak quasi-static nuclear dipolar fields. In the sample at 20 keV (not shown in Fig. 4), the fitted magnetic fraction and the magnitude of the internal field approach zero near 65 K indicating that the entire sample is paramagnetic above this temperature. At 5 keV, much closer to the surface, the transition temperature is much higher (see purple points in Fig. 4). In the sample, the surface has a magnetic ordering temperature close to K (see green points in Fig. 4), but the bulk is non magnetic with no sign of static magnetism. This is evident from the fast relaxing signal at low energies (below 10 keV) as seen earlier in Figs. 1 and 2, compared to the slowly relaxing signal at higher energies. At the highest Ce concentration () there is no magnetic ordering at any temperature. It is important to note that the muon is sensitive to any static order both short and long range. The absence of any static magnetic field in a SR experiment (e.g. in the sample) implies there is no static 3D magnetism of any type (long range antiferromagnetism, spin density wave, or spin glass). On the other hand it does not exclude quasi 2D magnetic order in the planes as long as the internal fields are still fluctuating fast compared to s.
ii.3 Phase diagram
The results are summarized in the phase diagram in Fig. 5. The green circles (5 keV), yellow squares (11 keV), and black diamonds (21 keV) are LE-SR results for 200 nm thick films. The red and blue triangles are for the 100 nm thick film. The brown band defines the magnetic phase boundary from the LE-SR measurements, with static magnetism to the left. The width of the band originates from the depth dependence of the magnetism in the film (i.e. the phase boundary for the bulk of the film is the left side of the brown band). The dashed blue curve is the superconducting transition temperature obtained from resistivity, while the cyan dashed curve is the magnetic transition temperature obtained from angular magneto-resistance (AMR) measurements (19).
There is a considerable difference in the magnetic phase boundary seen with ZF-SR and AMR as reported in Fig. 5. In particular LE-SR shows a much narrower region of overlap of the static magnetism with superconductivity. Comparing the results from these two techniques provides an important insight into the nature of the magnetism. Angular magneto resistance (AMR) and LE-SR are sensitive to magnetic fluctuations on a much different time scale. As mentioned above LE-SR measures the temperature where the internal fields become static on the scale of the period for the Larmor frequency associated with the internal fields, which in this case is s. On the other hand the AMR detects magnetism when the fluctuations become much slower than the relaxation time ( s). Normally this would not shift the ordering temperature very much. However LCCO is highly anisotropic so that the coupling between moments in the plane is considerably larger than the coupling between planes. Thus the ordering temperature detected by AMR is determined by the in-plane coupling between spins . Just below this temperature the moments are highly correlated in the CuO plane but are still dynamic on the time scale detectable with LE-SR. The ordering temperature in LE-SR defines the boundary for static 3D ordering, which depends and varies logarithmically with the much smaller magnetic coupling between planes ’ (29). This gives new insight into the relationship between the magnetism in the region of the SC dome before the Fermi surface reconstructs. Apparently the moments are still fully developed above but are very dynamic on the muon time scale. Finally we note that the LE-SR measurements are performed in zero applied magnetic field whereas the AMR measurements are made in a magnetic field of 14 T applied in plane. It is possible that the large magnetic field used in AMR would frustrate the 2D antiferromagnetic order and thereby suppress ordering temperature. This could explain why the AMR 2D phase boundary seems to cross the LE-SR 3D phase boundary at .
From the LE-SR results in the bulk of the film, we conclude that the static 3D static magnetism disappears close to where superconductivity begins and thus there is little overlap region between static 3D magnetism and superconductivity. This is similar to what is found in LaSrCuO except of the evidence of phase separation in bulk, such that a weak spin glass phase with a small volume fraction extends into the superconducting dome in LaSrCuO (30). It is possible that a similar spin glass phase exists in LCCO in bulk but the glass temperature would have to be lower than our base temperature of 4 K. Thus the present results are consistent with a sharp boundary and competition between the superconductivity and magnetism without microscopic coexistence. However comparing these results with AMR leads us to the conclusion that above the SC dome the moments are fully developed and highly correlated within the CuO plane but are still fluctuating without any 3D order. Striking similarities between hole and electron-doped cuprates are also reported by transport studies (5). It is now clear there are at least two quantum critical points in LaCeCuO. The lower one is near where the antiferromagnetic disappears, whereas, the Fermi surface reconstructs at a much higher doping level above (7); (5). It is worth noting that recent studies of certain heavy fermions show a similar phase diagram with two critical points or possibly a line of quantum critical points called a quantum critical phase in between (31).
The difference in the magnetic behaviour of the surface compared to the bulk of the film is also very interesting. It is possible that small structural changes, expected near a free surface, are enough to tip the balance towards magnetism. This would be consistent with a small difference in (free) energy between competing phases near a critical point. As mentioned previously LE-SR measurements were also taken on a thinner (100 nm) sample with . The magnetic properties of the free surface were the same as in the 200 nm thick film with the same value of . However in the thinner film it was possible, using 21 keV muons, to probe the interface between the film and the substrate. These measurements showed this interface was non-magnetic, implying the magnetism is not simply a consequence of the broken translational symmetry and resulting boundary conditions.
In conclusion, we have conducted a depth-resolved LE-SR study of magnetism in LaCeCuO films close to the magnetic-superconducting transition region. We find that the near surface region tends to be more magnetic than the bulk of the films. The enhanced magnetism in the near surface region of is a property of the free surface and absent near the substrate interface. The disappearance of strong static magnetic order in the bulk occurs just below , which is close to where superconductivity appears. Above this critical value of the moments may exist but are fluctuating rapidly on the muon time scale. The AMR results there is a high degree of 2D order above the superconducting dome. Thus above the superconducting transition the system is in a dynamic magnetic state with strong antiferromagnetic correlations. This is similar to the hole-doped counterpart LaSrCuO and consistent with a competition between order parameters. These results suggest there may be a higher degree of particle-hole symmetry in this critical region of the phase diagram than previously thought.
LE-SR experiment: The measurements are performed using low-energy SR technique, where an intense high-energy beam of muons is moderated in a solid Ar film. A small percentage of the incident muons emerge from the argon surface at low-energy (10 eV). These are subsequently accelerated to 15 keV and transported to the sample chamber and SR spectrometer. The samples are mounted onto a Ag coated metal plate which is electrically isolated from the cold finger of the cryostat and biased to a high voltage ranging from to kV. This allows the implantation energy of the muons to be adjusted between about and keV. The measurements reported here were performed in zero external magnetic field (ZF) such that any stray magnetic field at the sample was less than 0.01 mT. Measurements in a longitudinal field and transverse field were also performed but are not shown here. The time evolution of the muon spin polarization , which is measured through the properties of the muon decay, depends on the static and fluctuating components of the internal magnetic field at the site of the muon (32).
Samples: The -axis-oriented LaCeCuO films were deposited directly on insulating (100) SrTiO substrates by a pulsed laser deposition technique utilizing a KrF excimer laser as the exciting light source. The three films had Ce concentrations of , and . Since the oxygen content has an influence on both the superconducting and normal state properties of the material, the annealing process is optimized for each . We studied three Ce concentration: 0.07 ( K), ( K), and 0.10 ( K). The sample with is metallic (with 0) above with a narrow transition width. The and samples showed an upturn (either in a field or ZF) at low temperature. Each LaCeCuO sample consisted of four identical pieces with a total area of 22 cm and was attached to the sample holder with a conductive Ag paint.
Analysis: The ZF LE-SR data is fit (33) to a sum of two functions; , corresponding to the signal from the fraction of muons stopping in LCCO, and the remaining fraction of muons ( ) stopping in the Ag mounting plate. The latter is accounted for with an exponential with a small relaxation rate of s. The signal from the LCCO may be decomposed into a magnetic fraction and a non-magnetic or paramagnetic fraction ;
The magnetic fraction is fit to a phenomenological function consisting of three relaxation functions: (1) a weakly relaxing Kubo-Toyabe function due to nuclear moments ( s). This term dominates well above the ordering temperature when the electronic moments are rapidly fluctuating but has little effect at low temperatures where the electronic moments dominate. (2) A slowly relaxing exponential, , attributed mostly to slow fluctuations in the nuclear dipolar field. (3) A term which takes into account a broad distribution of large quasi-static internal magnetic fields from the magnetism; . For simplicity, we assume the direction of the internal magnetic field is random in orientation as in a spin glass, giving rise to the geometric “” non precessing component appearing in . The parameter represents the average internal magnetic field from the static ordering, whereas characterizes the width of this field distribution. This term dominates at low temperatures in the ordered state. Since there is an interplay between and , the latter is parametrized as , where is a fitted parameter. This assumes that the static width of the magnetic field distribution is proportional to the average field. In the non-magnetic fraction of the sample only the first two terms in the relaxation function contribute. It should be noted that other fitting functions have been tested, leading to very similar results.
Acknowledgments: The authors H.S. and E.M. acknowledge the financial support of the MANEP program, and R.F.K. the support of NSERC. The authors J.Y. and R.L.G. acknowledge the support of Maryland Center for Nanophysics and Advanced Materials and the NSF grant under DMR-1104256. We would like to thank Johnpierre Paglione for important discussions at the early stages of this work, and Mohamed Azzouz for commenting and critical reading of this manuscript.
Contributions: Project planning: R.F.K., H.S., R.L.G., and Z.S.; Sample growth: J.P., K.J., and R.L.G.; LE-SR experiment: H.S., Z.S., R.F.K., E.M., T.P., H.L., A.S., and W.A.M.; Data analysis and interpretation: H.S., and R.F.K.; Draft writing: H.S. and R.F.K., with contributions and/or comments from all authors.
Competing financial interests: The authors declare no competing financial interests.
- Bennemann, K. H., Ketterson, J. B. (Editors). Superconductivity. (Springer, 2008).
- Armitage, N. P., Fournier, P., and Greene, R. L. Progress and perspectives on the electron-doped cuprates. Rev. Mod. Phys. 82, 2421 (2012), and references therein.
- Sachdev, S., and Keimer, B. Quantum criticality. Phys. Today 64, 29 (2011).
- Tsuei, C. C., and Kirtley, J. R. Phase-sensitive evidence for -wave pairing symmetry in electron-doped cuprate superconductors. Phys. Rev. Lett. 85, 182 (2000).
- Jin, K. et al. Link between spin fluctuations and electron pairing in copper oxide superconductors. Nature 476, 73 (2011).
- Motoyama, E. M. et al. Spin correlations in the electron-doped high-transition-temperature superconductor NdCeCuO. Nature 445, 186 (2007).
- Helm, T. et al. Magnetic breakdown in the electron-doped cuprate superconductor NdCeCuO: the reconstructed Fermi surface survives in the strongly overdoped regime. Phys. Rev. Lett. 105, 247002 (2010).
- Sachdev, S. Quantum phase transitions (Cambridge University Press, Cambridge, UK, 2011).
- Tremblay, A.-M. S. In theoretical methods for strongly correlated systems, edited by A. Avella and F. Mancini (Springer, Berlin, 2011).
- Weber, C., Haule, K., and Kotliar, G. Strength of correlations in electron- and hole-doped cuprates. Nature Phys. 6, 574 (2010).
- Onose, Y et al. Charge dynamics in underdoped NdCeCuO: Pseudogap and related phenomena. Phys. Rev. B 69, 024504 (2004).
- Armitage, N. P. et al. Doping dependence of an n-Type cuprate superconductor investigated by angle-resolved photoemission Spectroscopy. Phys. Rev. Lett. 88, 257001 (2002).
- Zimmers, A. et al. Infrared properties of electron-doped cuprates: Tracking normal-state gaps and quantum critical behavior in PrCeCuO. Europhys. Lett. 70, 225 (2005).
- Kyung, B et al. Pseudogap and spin fluctuations in the normal state of the electron-doped cuprates. Phys. Rev. Lett. 93, 147004 (2004).
- Helm, T. et al. Evolution of the Fermi surface of the electron-doped high-temperature superconductor NdCeCuO revealed by Shubnikov-de Haas oscillations. Phys. Rev. Lett. 103, 157002 (2009).
- Matsui, H. et al. Evolution of the pseudogap across the magnet-superconductor phase boundary of NdCeCuO. Phys. Rev. B 75, 224514 (2007).
- Dagan, Y. et al. Evidence for a quantum phase transition in PrCeCuO from transport measurements. Phys. Rev. Lett. 92, 167001 (2004).
- Dagan, Y., Qazilbash, M. M., and Greene, R. L. Tunneling into the normal state of PrCeCuO. Phys. Rev. Lett. 94, 187003 (2005).
- Jin, K. et al. Evidence for antiferromagnetic order in LaCeCuO from angular magnetoresistance measurements. Phys. Rev. B 80, 012501 (2009).
- Morenzoni, E. et al. Low-energy SR at PSI: present and future. Physica B 289, 653 (2000).
- Boris, A. V. et al. Dimensionality control of electronic phase transitions in Nickel-Oxide superlattices. Science 332, 937 (2011).
- Suter, A. et al. Magnetism in the 2D limit and interface superconductivity in metal-insulator La(2-x)Sr(x)CuO4 superlattices. Phys. Rev. Lett. 106, 237003 (2011).
- Dunsiger, S. et al. Spatially homogeneous ferromagnetism of (Ga,Mn)As. Nature Materials 9, 299 (2010).
- Saadaoui, H. et al. Absence of spontaneous magnetism associated with a possible time-reversal symmetry breaking state beneath the surface of (110)-oriented YBaCuO superconducting films. Phys. Rev. B 88, 180501(R) (2013).
- Prokscha, T. et al. The new E4 beam at PSI: a hybrid-type large acceptance channel for the generation of a high intensity surface-muon beam. Nuclear Instruments and Methods A 595, 317 (2008).
- Eckstein, W. Computer Simulation of Ion-Solid Interactions (Springer, Berlin, 1991).
- Morenzoni, E. et al. Implantation studies of keV positive muons in thin metallic films. Nuclear Instruments and Methods B 192, 254 (2002).
- Matsuura, M. et al. Effect of a magnetic field on the long-range magnetic order in insulating NdCuO and non-superconducting and superconducting NdCeCuO. Phys. Rev. B 68, 144503 (2003).
- Yasuda, C. et al. Nel temperature of quasi-low-dimensional Heisenberg antiferromagnets. Phys. Rev. Lett. 94, 217201 (2005).
- Julien, M.-H. Magnetic order and superconductivity in LaSrCuO: a review. Physica B 329, 693 (2003).
- Coleman, P., and Nevidomskyy, A. H. Frustration and the Kondo effect in heavy fermion materials. J. Low Temp. Phys. 161, 182 (2010).
- Yaouanc, A., and Dalmas de Rotier, P. Muon spin rotation, relaxation, and resonance, (Oxford University Press, 2011).
- Suter, A, and Wojek, B. M. Musrfit: a free platform-independent framework for SR data analysis. Physics Procedia 30, 69 (2012).