## Magnetic and superconducting anisotropy in Ni-doped RbEuFe As 4 4 single crystals

K. Willa, 1, 2 M.P. Smylie, 1, 3 Y. Simsek, 1, 4 J.-K. Bao, 1 D. Y. Chung, 1 M. G. Kanatzidis, 1, 5 W.-K. Kwok, 1 and U. Welp 1

1

Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA Institute for Solid State Physics, Karlsruhe Institute of Technology, Karlsruhe D-76021, Germany 3 Department of Physics and Astronomy, Hofstra University, Hempstead, New York 11549 4 SUNUM-Sabanci University Nanotechnology Research and Application Center, Istanbul, Turkey 5 Department of Chemistry, Northwestern University, Evanston, Illinois, 60208, USA

2

(Dated: August 5, 2019)

We investigate the effect of Ni doping on the Fe-site in single crystals of the magnetic superconductor RbEuFe As 4 4 for doping concentrations of up to 4%. A clear suppression in the superconducting transition temperature is observed in specific heat, resistivity and magnetization measurements. Upon Ni-doping, the resistivity curves shift up in a parallel fashion indicating a strong increase of the residual resistivity due to scattering by charged dopand atoms while the shape of the curve and thus the electronic structure appears largely unchanged. The observed step ∆ C/T c at the superconducting transition decreases strongly for increasing Ni doping in agreement with expectations based on a model of multi-band superconductivity and strong inter-band pairing. The upper critical field slopes are reduced upon Ni doping for in- as well as out-of-plane fields leading to a small reduction in the superconducting anisotropy. The specific heat measurements of the magnetic transition reveal the same BKT behavior close to the transition temperature T m for all doping levels. The transition temperature is essentially unchanged upon doping. The in to out-of-plane anisotropy of Eu-magnetism observed at small magnetic fields is unaltered as compared to the undoped compound. All of these observations indicate a decoupling of the Eu magnetism from superconductivity and essentially no influence of Ni doping on the Eu magnetism in this compound.

## I. INTRODUCTION

The family of iron-based superconductors has been constantly growing, now comprising many members of the so-called 11, 111, 1111 and 122 materials 1-9 . In the majority of these materials, superconductivity arises from an antiferromagnetic parent compound upon electron, hole or isovalent doping or mechanical pressure 1,2,4-14 . The so-called 1144 compounds with composition AAe Fe As 4 4 ( A = K, Rb, Cs; Ae = Ca, Sr, Eu) 15-19 form the most recently discovered members of this family. In these materials, the large difference in the ionic radii of the A and Ae atoms leads to the formation of alternating A and Ae-layers between the Fe 2 As -planes producing an asymmetric environment 2 for the Fe 2 As -layers. 2 Contrary to the doped 122 compounds, the 1144 compounds are superconducting in their stoichiometric state with T c reaching ≈ 37 K. A formal charge count yields that stoichiometric 1144 are intrinsically hole doped to ≈ 0.25 holes/Fe, close to the doping level that yields optimum T c in the 122 compounds. Recently, it has been realized 20 that the high purity and symmetry properties of 1144 materials such as CaKFe As 4 4 may open a new platform for the observation of topological bandstructures and Majorana zero modes. The A EuFe As 4 4 compounds represent a peculiar subgroup of this new family as the planes of Eu 2+ ions can order magnetically as has been shown previously in EuFe As 2 2 21-25 . Initial experiments on polycrystalline RbEuFe As samples indicate ferromagnetic in-plane or4 4 dering at 15K deep within the superconducting phase

( T c = 37 K) 17,18 . Recent studies on single crystals 26 affirmed these findings and allowed for the study of the magnetic and superconducting anisotropy 27,28 , demonstrating a low superconducting anisotropy of 1.8, highly anisotropic quasi-2D Eu-magnetism, and an associated Berezinskii-Kosterlitz-Thouless (BKT) transition of the Eu 2+ moments. Large upper critical fields were observed in pulsed field measurements 29 . Recent neutron scattering data indicate a modulated, possibly helical stacking of ferromagnetic Eu-layers along the c axis with a period of four c -axis lattice constants 30 . Furthermore, a study on polycrystalline samples 31 revealed that T c is suppressed upon Ni-substitution on the Fe-site, while the magnetic ordering temperature is unchanged. Meanwhile, Ca-substitution on the Eu-site 32 suppressed the magnetic ordering temperature in polycrystalline samples without changing T c . These results and recent high pressure studies 33 suggest an almost complete decoupling of superconductivity hosted by the Fe As 2 2 planes and magnetism within the Eu layers. This model system therefore provides an ideal playground to selectively tune the relative importance of magnetism and superconductivity, rendering these compounds promising for further investigations. In this work we investigate single crystals of Ni-doped RbEuFe As 4 4 with a particular focus on the anisotropy of the superconducting and magnetic transition and its doping dependence.

FIG. 1. XRD measurements of (00l) peaks on single-crystal RbEu(Fe 1 -x Ni ) As x 4 4 for x = 0 00 (black), . x = 0 02 (blue), . x = 0 03 (orange), . x = 0 04 (green). . As doping increases, the lattice peaks systematically shift to higher angles as the lattice shrinks; compare table I. The inset shows the (009) peak. The asterisk marks a small amount of 122-phase compound that grows at the surface of the 1144-phase. No other phases have been observed.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-b66e3d17-d5a9-438c-96dc-af8d51f85ef9.jpg)

## II. EXPERIMENTAL METHODS

High quality single crystals RbEu(Fe 1 -x Ni ) As x 4 4 were grown out of RbAs flux with a similar method to the one used for the undoped RbEuFe As 4 4 in Ref. 26 . In order to provide the Ni doping, an extra NiAs precursor was synthesized by reacting Ni and As powder under a 1:1 ratio in an evacuated quartz tube at 823K for 12 hours. The sintered NiAs precursor was then ground into a fine powder in an agate mortar. To achieve a nominal Ni doping x ( x = 0 03, 0.05, 0.07), EuAs:Fe . 2 As:NiAs:RbAs = 1:1x :2 x :15 were used for the RbAs flux growth. After the growth the RbEu(Fe 1 -x Ni ) As x 4 4 crystals are obtained by removing the RbAs flux with reagent alcohol. This procedure yields thin flat plates with typical sizes of 300 × 300 × 50 µ m , with the tetragonal [110] and [110] 3 orientations parallel to the long edges, and the [001] orientation perpendicular to the plate.

We performed x-ray diffraction (XRD) measurements on single crystals using a Powder X-ray PANalytical diffractometer X'pert with a Cu-K α source. Plate-like crystals were well oriented in a flat zero-background single crystalline silicon sample holder to make sure that peaks with the index (00l) were observed, see Figure 1. Magnetization measurements were performed in a Quantum Design MPMS-7 system with samples mounted in a notched quartz rod for H ‖ [001] and on a smooth quartz rod for H ‖ [110]. For low-field measurements, careful background field calibration was performed prior to each run with a Pd reference sample. For magnetotransport measurements, four-bar gold patterns were deposited onto thin rectangular-shaped crystals with the help of a shadow-mask technique then thin gold wires were attached onto the gold bars with silver epoxy. The measurements were performed inside a 9 -1 -1T 3-axis

FIG. 2. Resistivity as a function of temperature of single crystals of RbEu(Fe 1 -x Ni ) As x 4 4 with current in-plane. As the doping increases, the normal-state resistivity increases (see inset) and the transition temperature decreases. Qualitatively the temperature dependence of the resistivity remains unchanged.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-ffbe5c3b-f1f6-429b-91fc-702f4cc40a4e.jpg)

AMI superconducting vector magnet in a 4-point geometry and applying DC currents of 1 mA along the in-plane direction. The ac specific heat is measured on a SiN calorimetric membrane 34,35 . For this purpose, a small ( ≈ 100 × 100 × 25 µ m ) platelet-shaped single crystal was 3 mounted on the nanocalorimeter platform with apiezon grease and inserted into the same three axis superconducting vector magnet as used for resistivity measurements. The sample was then subjected to a small oscillatory heating current with frequencies usually in the range of 1Hz and the response of the sample was measured using a SynkTek MCL1-540 multichannel lock-in system.

## III. EXPERIMENT

## A. Diffraction

We used Energy-dispersive X-ray spectroscopy (EDS) to determine the actual Ni content of the RbEu(Fe 1 -x Ni ) As x 4 4 single crystals. For this, about 10 crystals from each doping series were investigated. We found that the nominal dopings of x = 0 03, . x = 0 05, and . x = 0 07 correspond to an actual doping . concentration of x = 0 02, . x = 0 03, . and x = 0 04 . respectively; from here on, samples will be referred to by their average EDS doping values. The spread of doping concentrations in any batch is around ± 0.1 %. In XRD measurements using Cu-K α 1 radiation we observe (00l) peaks of the of the single-crystal RbEu(Fe 1 -x Ni ) As x 4 4 with doping levels of x = 0 00 (black), . x = 0 02 (blue), . x = 0 03 (orange), . x = 0 04 (green), as shown in Fig. 1. . With increasing doping, the lattice peaks systematically shift to higher angles, see inset of Fig. 1 where the (009) peak is shown. This indicates a shrinking of the c -axis lattice constant.

FIG. 3. Temperature dependence of the resistivity of RbEu(Fe 1 -x Ni ) As x 4 4 in in fields of up to 9T. Upper panels are for fields in-plane, lower panels for fields along the c axis. The color code is the same for all panels. For all doping level, the superconducting transition is suppressed and broadens with increasing fields, more for fields along the c axis then for in-plane fields. The suppression with fields gets larger for increasing Ni doping. The temperature axis of all panels spans 5 K whereas the vertical axis is adjusted to accommodate the large increase in resistivity with doping.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-0755b17a-2717-41b8-8030-6d03bd221cfe.jpg)

We determine the c -axis lattice parameter from leastsquares fits with a zero-shift parameter of all observed peaks as shown in Fig. 1. We estimate the uncertainty to be less than 7 × 10 -4 ˚. A A few samples show a small amount of nonsuperconducting Eu122-impurity phase that grows on the surface of the RbEuFe As crystals. 4 4 No other impurity phases were observed. The c axis decreases almost linearly with increasing Ni-doping at a rate of ≈ -8 × 10 -3 ˚ %Ni which is in good agreement A / with low-doping data on polycrystalline samples 31 and is consistent with Vegard's law implying the uniform incorporation of Ni.

## B. Resistivity

Figure 2 shows the temperature dependence of the in-plane electrical resistivity ρ ab of single crystals of RbEu(Fe 1 -x Ni ) As x 4 4 at different doping levels. The superconducting transition temperature T c is clearly suppressed upon increasing Ni doping at a rate of ≈ -2.4 K/%Ni, consistent with the polycrystalline data at low Ni-doping 31 . In all samples the resistive transitions remain sharp with a transition width below 0.5 K, indicating uniform doping and single-phase material. As seen in the inset of Fig. 2, the salient effect of Ni-doping is the remarkable parallel upwards shift of the resistivity curve. The resistivity at T c vs. T c can be seen in the inset of Fig. 4. While the residual resistivity increases strongly, the temperature dependence of the resistivity is not altered. In particular, an upturn (on decreasing temperature) of the resistivity which was observed around 30 K on polycrystalline samples 31 at doping levels of 3% and above and interpreted as the reemergence of the spin density wave on Fe, is not apparent in our data. As electron scattering by phonons as well as spin fluctuations depend on the details of the electronic band structure, our results imply that the effects of doping-induced changes of the band structure on electron transport are small. Recent first principals calculations 36 of the electronic properties of Ni-doped RbEuFe As 4 4 revealed that Ni-doping induces electron doping accompanied by an upwards shift of the Fermi energy. For 6.25% Ni, the partial density of states derived from the Fe-d x y z ( ) , d z 2 and d x 2 -y 2 orbitals were found to decrease by approximately 10 to 20% while for the d xy states it increases by 4%. Since T c depends exponentially on the density of states, these small changes result in a clear suppression of T c . In contrast, the resistivity depends much weaker on the density of states, and the large observed increases of the resistivity are attributable to enhanced disorder scattering due to the charged Ni-dopants. For all doping levels, we do not observe a reentrant resistive state associated with the onset of Eu sublattice magnetic order;

FIG. 4. Temperature dependence of C/T of RbEu(Fe 1 -x Ni ) As x 4 4 in zero field. A clear step marks the onset of superconductivity, while the magnetic transition of the Eu moments is seen as a kink around T m = 15K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-488d0e36-be22-43a5-89ba-9c8e28b0698a.jpg)

this is in agreement with earlier reports on the polycrystalline material 31 but in contrast to other Eu-containing iron-based superconductors e.g. Eu(Fe 1 -x Ir x ) 2 As 2 37 .

Fig. 3 shows the the superconducting transitions of the pristine and doped samples in in-plane (top row) and out-of-plane fields (bottom row). For all samples, the field-induced suppression and broadening of the superconducting transition is stronger for fields along the c -axis then for in-plane fields and also increases for increasing Ni content yielding a modest superconducting anisotropy of 1.8 in the pristine sample which slightly decreases upon Ni-doping (see below, Fig. 7). The superconducting transition temperature as a function of field was extracted as the mid-point of the transition and is used to construct the phase boundaries, which can be seen in figure 7. The high-field transitions of the pristine sample are characterized by an unusual sharpening on decreasing temperature, clearly seen in the 9T ‖ c -data. This feature is associated with the vortex lattice melting transition 38 and is suppressed due to the increased disorder in the doped samples. Furthermore, the negative normal state magnetoresistance that is clearly seen in the pristine sample (see also Ref. 27 ) is strongly reduced with increasing Ni-doping, consistent with the increase in impurity scattering.

## C. Specific Heat

Figure 4 shows the zero field specific heat of the pure RbEuFe As 4 4 and for the three different Ni doping levels. Due to the small sample size, the determination of the molar volume for each sample introduces uncertainties. Therefore, in order to facilitate the comparison of the data from different samples, we scale the normal state molar specific heat of the Ni-doped samples to the normal state specific heat of the pure RbEuFe As compound. 4 4 This step is motivated by the observation that substitut-

FIG. 5. Specific heat of the superconducting transition of RbEu(Fe 1 -x Ni ) As x 4 4 for different doping levels. Alinear normal state background has been subtracted. The step height of the jump in C/T decreases with increasing Ni concentration from 0.23 J/molK 2 down to 0.12 J/molK . The observed 2 upturn of C/T at around 23K is due to magnetic fluctuations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-533fa6c5-a43e-4d05-92de-d279c6519c7a.jpg)

ing Fe by Ni induces essentially no change in the lattice specific heat and very little in the normal state electronic specific heat. The two clear features in C/T are (i) a step signaling the superconducting phase transition at T c and (ii) a kink at the magnetic transition temperature T m while the superconducting transition temperature is clearly suppressed-both in temperature and in step-height-upon increasing Ni doping, the shape and transition temperature T m of the magnetic transition change only weakly. Even for the largest Ni concentrations investigated here, the Eu moments order below the superconducting transition. Consistent with the resistivity data, there is no additional feature in the specific heat between 2 K and room temperature that would indicate a spin-density wave (on the Fe site) as is seen in most 122 parent compounds 3-9,22,25,39 . Figure 5 shows the superconducting contribution to the specific heat of the four samples obtained by subtracting a linear extrapolation of the normal state data above T c . The transition temperature decreases from T c = 36.9K for the undoped compound to T c = 26.7K for 4% Ni doping. While the step height is reduced from 0.23 J/mol K 2 to 0.12 J/mol K . 2 The evolution of the step size with Nidoping can be seen in inset of figure 4. In single-band weak-coupling BCS theory, the specific heat anomaly is given as ∆ C/T c = 1 43 . γ n , where γ n is the Sommerfeld coefficient of the electronic specific heat. Thus, the rapid suppression of ∆ C/T c with Ni-doping would indicate a strong reduction in the density of states, contrary to the transport data (see above). Alternatively, however, a strong dependence of ∆ C/T c on T c has been observed in a wide variety of iron-based superconductors, a phenomenon known as the BNC scaling 40 . It is believed that this behaviour arises in multi-band superconductors with strong inter-band pairing interaction 41 and thus might be

1.8

FIG. 6. C/T versus temperature for the three different Ni concentrations x = 0 02 (a), . x = 0 03 (b) and . x = 0 04 (c). . Fields of up to 9T were applied along the c axis suppressing the superconducting transition. The broadening of the transitions in field is considerably larger for higher Ni doping. The magnetic transition is transformed into a broad hump that moves to higher temperatures with fluctuations extending well above the superconducting transition. The insets (d) and (f) show a 2% and 4% doped RbEu(Fe 1 -x Ni ) As x 4 4 single crystal, respectively. Inset (e) shows C/T (normal state background subtracted) measured in 3 T for 0, 2% and 3% Ni doping. In the pristine sample the additional step below T c (marked by a dotted line) indicates the vortex lattice melting. This is not visible anymore for the higher doped samples.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-b8524a75-fd35-42a0-9852-103c12e1d532.jpg)

responsible for the strong reduction in ∆ C/T c .

Figure 6 shows C/T of the three Ni-doped samples in fields of up to 9T applied along the c -axis. In analogy to the magnetotransport data in Fig. 3, with increasing fields, the superconducting transition is suppressed and broadened. The broadening in field is considerably larger for the higher doped samples. At the same time, the magnetic transition transforms from a kink into a crossover that moves to higher temperatures in increasing field. The associated specific heat contribution due to strong magnetic fluctuations persists far above the superconducting transition. Previous specific heat measurements on pristine RbEuFe As revealed a step in 4 4 C/T slightly below the superconducting transition 38 that lines up with the steep drop in the resistance and is a signature of the vortex lattice melting transition. Inset (e) of Fig. 6 shows C/T (normal state background subtracted) measured in 3 T for 0, 2% and 3% Ni doping. In the pristine sample the additional step below T c (marked by the dotted line) indicates the vortex lattice melting. Upon Ni-doping, this feature is suppressed due to increased impurities and pinning sites in the sample which is consistent with the evolution of the magnetoresistance.

We have extracted the superconducting phase boundaries of the superconducting transition using an entropy conserving construction 42 and plotted it together with those obtained from resistive mid-points and magnetization, see Fig. 7. The agreement between resistivity and specific heat derived data is generally good. The difference in T c for the 3% sample can be attributed to measurements on crystals from different batches having slightly different doping levels. The upper critical field slopes dH c 2 /dT for both in-plane and out-of-plane fields are suppressed with increasing Ni doping. Taking the resistivity data, we observe that the c -axis upper critical field slope decreases from around 4 T/K for the undoped sample to 3.3 T/K for 4% doping, whereas the in-plane direction changes from 6.3 T/K ( x = 0 00) . to 5 T K ( / x = 0 04). . This results in a

H [T]

FIG. 7. Phase boundary of pure and Ni doped RbEuFe As 4 4 single crystals for in and out-of-plane fields. Data are obtained from resistivity (open symbols and dashed lines) and specific heat (closed symbols and solid lines). The superconducting transition obtained from susceptibility measurements is indicated as well (stars).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-2d032e10-9a3f-400a-9b29-11f0af672cd4.jpg)

9

FIG. 8. Specific heat near the magnetic transition for 0% (a), 2% (b), 3% (c), and 4% (d) Ni doping. A 9 T background curve has been subtracted. The transition temperature T m = 15 K, here the peak in C , appears unchanged upon Ni doping as does the shape of the transition. Small fields along the c axis shift the transition to lower temperatures while fields in-plane transform the transition into a cross-over that moves to higher temperatures. e) shows the phase diagram extracted from specific heat, magnetization and resistivity for the magnetic and superconducting transitions in zero field upon increasing Ni-doping.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-c17f377b-59ad-410c-996e-b9ccc1c06ce7.jpg)

small reduction of the superconducting anisotropy from Γ = ( dH ab c 2 /dT / dH ) ( ) c c 2 /dT ) = 1 6 to Γ = 1 5. . . Naively, one would expect that with increasing disorder, as is suggested by the strong increase of the residual resistivity, the upper critical fields H c 2 and dH c 2 /dT increase as the coherence length decreases. However, T c decreases also with Ni-doping. Since our data are taken close to T c , we present simple estimates using the Ginzburg-Landau formalism giving H c 2 ( T ) = Φ 0 / 2 πξ 2 ( T ). In the clean limit, ξ (0) &lt; l = v F τ with l , v F and τ the electron mean free path, the Fermi velocity and scattering time, respectively, the coherence length is given by the BCS coherence length ξ (0) = 0 74 . ξ 0 = 0 76 . glyph[planckover2pi1] v F / ∆(0) ≈ glyph[planckover2pi1] v F / 2 38 . k B T c . This results in H c 2 (0) ∝ T 2 c 2 /v 2 F and an upper critical field slope dH c 2 /dT c ∝ T /v c 2 F . Thus, assuming that changes in the electronic structure are not dominant (see above), the upper critical field and the slope of it indeed decrease with decreasing T c , that is, with increasing Nidoping. The same conclusion would also be reached in the dirty limit, when the coherence length is given as ξ (0) = 0 855( . ξ 0 l ) 1 / 2 resulting in H c 2 (0) ∝ T /v c 2 F τ . Since ρ ∝ 1 /τ is proportional to T c (with an offset) one again finds to leading order H c 2 (0) ∝ T 2 c 2 and a slope that decreases with increasing Ni-doping.

In order to gain insight into the nature of the magnetic transition, we performed detailed specific heat measurements in the close vicinity of the transition. In zero field, the transition is characterized by a kink in the specific heat for all doping levels, a non-singular behavior that was already observed for the undoped compound and identified as a BKT transition 28 . As noted above, the transition temperature of T m = 15 K is largely unaffected by the Ni doping. Applying small fields along the c axis shifts the magnetic transition to lower temperatures while applying the field in the plane replaces the transition with a crossover. The very similar fieldresponse of the specific heat of the doped samples (Fig. 8) and of the undoped samples (Ref. 28 ), leads us to conclude that Ni-doped RbEuFe As , just as the stoichio4 4 metric compound, belongs to the universality class of the anisotropic 2D XY model, with the Eu-moments having an easy-plane anisotropy, strong ferromagnetic intralayer coupling and very weak interlayer interactions leading to the BKT transition of the moments in the plane near T m = 15 K. This transition will be suppressed to lower temperatures for fields along the c axis, whereas it will be replaced with a broad crossover at higher temperatures for fields in the plane since any field in the plane destroys the rotational symmetry necessary for a BKT transition.

## D. Magnetization

We study the magnetic state of RbEu(Fe 1 -x Ni ) As x 4 4 by measuring the zero field-cooled (ZFC) and field cooled (FC) temperature dependence of the magnetization in different applied fields as well as the field dependence of the magnetization at different temperatures with the

FIG. 9. The ZFC and FC magnetic susceptibility of RbEu(Fe 1 -x Ni ) As x 4 4 single crystals with H = 1mT along the [001] direction. With increasing doping, the onset T c shifts lower while the transitions remain parallel, indicating homogenous doping. In this orientation, the diamagnetism masks most of the Eu magnetism.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-7f7bec58-6199-46d7-93e3-9a4dab339b3c.jpg)

field applied perpendicular or parallel to the c axis. Fig. 9 shows ZFC and FC measurements of magnetic susceptibility, χ = M/H , for x = 0 00, . 0 02, . 0 03, . and 0 04 . with H = 10G applied along [001]. Slight differences in demagnetization factors may cause variations in the absolute magnitude of susceptibility. As doping increases, the onset of diamagnetism-taken as T c -monotonically decreases from 36 5K to 27 5K, in agreement with the re-. . sistivity and the specific heat data. Due to the plate like sample geometries, magnetization data for H ‖ [001] contains large superconducting contributions, making the Eu magnetism difficult to discern in this geometry.

Applying magnetic fields along [110] decreases the contribution of diamagnetism, allowing for a more pro- nounced development of the Eu magnetism, as shown in Fig. 10. In the magnetization measurements, we find that T m , defined as the peak in the ZFC susceptibility below T c , is slightly suppressed from 15 3K to 15 0K by low Ni . . doping. This is in rough agreement with the specific heat data showing that the magnetic ordering temperature is essentially independent of Ni-doping.

In higher fields, the magnetic susceptibility for FC measurements at all doping levels is qualitatively very similar to the undoped compound 27 as shown in Fig. 11. In 0.1 T we observe a large anisotropy in the low temperature susceptibility with χ ab glyph[greatermuch] χ c under FC conditions for which vortex pinning effects are small. Thus, the strong easy-plane anisotropy of the Eu-moments observed in the undoped compound is preserved with Nidoping. At 1 T and above 50 K the data are well described by a Curie-Weiss law χ T ( ) = χ 0 + C/ T ( -Θ ), c yielding Θ c near 23 K and an effective moment µ eff = 2 837 . C 1 / 2 close to the expected Eu 2+ effective moment of µ eff = gµ B √ S S ( +1) = 7 94 . µ B / Eu (with g = 2 and S = 7/2) suggesting that these low levels of Ni doping do not substantially affect the Eu magnetic interactions as expected in numerical simulations 36 . Fit values are shown in table I. The higher value of µ eff found for x = 0 03 is likely due to surface EuFe . 2 As 2 phase on the sample, which would contribute twice as much Eu per volume as the 1144 phase would. The positive CurieWeiss temperature values suggests predominantly ferromagnetic interactions between the Eu-moments consistent with what has been seen in the undoped compound. Within a 2D Heisenberg model that includes in-plane versus out-of-plane spin anisotropy 28 , the anisotropy of the Curie-Weiss temperatures between in and out-of-plane fields is a measure of the anisotropy of the Eu spins. It is essentially unchanged upon Ni doping consistent with the data obtained at low temperatures.

[]

FIG. 10. The ZFC and FC magnetic susceptibility of RbEu(Fe 1 -x Ni ) As x 4 4 single crystals with H = 1 mT applied along the [110] direction. Solid circles mark ZFC, hollow triangles, FC. Hollow squares represent cooling in zero field, applying 10 mT, warming to 20 K, and then measuring FC. As doping increases, the peak of the magnetic transition is almost unchanged showing that magnetism is essentially unaffected by Ni doping.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-0f2737e5-a251-4867-a642-c4bfe7a14d17.jpg)

1/(

-

) []

FIG. 11. Temperature dependence of field-cooled magnetic susceptibility measured upon warming in RbEu(Fe 1 -x Ni ) As x 4 4 single crystals: (a) x = 0 02. . (b) x = 0 03. . (c) x = 0 04. . For all panels, colored symbols represent H ‖ [110], black symbols represent H ‖ [001]. For µ H 0 = 0 1T, the large in- to out-of-plane anisotropy is present at all doping levels. . The right-hand axis in each panel plots the inverse susceptibility measured for fields of 10 T; the low-temperature behavior is more isotropic, and the high-temperature ( &gt; 50K) data can be well described by Curie-Weiss behavior.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-19d9a6d2-effc-4f93-8199-757973c5ca29.jpg)

FIG. 12. Magnetization of the Eu-sublattice vs applied field H ‖ [110] (a-c) and along H ‖ [001] (d-f) for different doping levels and at various temperatures. The insets show the as-measured magnetization hysteresis loops. At high fields, the saturation magnetization approaches ≈ 300 emu/cm , approximately that of the pristine material. 3 In all measurements, the slope of the magnetization increases as T m is approached from above, and decreases for T &lt; T m as the magnetic lattice becomes stiffer.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-1f17b729-f0e7-4b19-ad1a-8d522398320b.jpg)

TABLE I. Properties of single crystals of RbEu(Fe 1 -x Ni ) As x 4 4 derived from x-ray diffraction, resistivity, specific heat and magnetization measurements. Θ xy C denotes the Curie temperature for in-plane fields and Θ z C the one for fields along the c axis.

|   x [ - ] |   c [A] ˚ |   T C [K] |   ρ @ T c [ µ Ω cm] |   ∆ C/T c [J/molK 2 ] |   µ eff [ µ B /Eu] |   Θ xy C [K] |   Θ z C [K] |   Θ xy C / Θ z C [ - ] |   M sat [ µ B /Eu] |
|-----------|-----------|-----------|---------------------|-----------------------|--------------------|--------------|-------------|------------------------|--------------------|
|      0    |   13.3047 |      36.9 |                  22 |                  0.23 |               7.75 |        23    |       21.4  |                  1.075 |                7.2 |
|      0.02 |   13.2911 |      32.5 |                  54 |                  0.19 |               7.87 |        23.18 |       21.94 |                  1.057 |                6.4 |
|      0.03 |   13.2825 |      30.7 |                  94 |                  0.14 |               8.89 |        22.53 |       20.87 |                  1.078 |                8.1 |
|      0.04 |   13.2733 |      26.9 |                 139 |                  0.12 |               7.97 |        22.84 |       21.04 |                  1.086 |                6.5 |

The insets of Figs. 12 (a)-(f) show magnetization hysteresis loops measured at 5 K for in and out-of-plane fields. The superposition of a hysteretic superconducting signal and a ferromagnetic-like signal is clearly seen, especially for H ‖ [001]. The large sample cross-section and high critical current density for H ‖ [001] makes this an expected behavior. In comparison to the x = 0 02 . and x = 0 04 samples, the . x = 0 03 sample displays a . lower superconducting hysteresis corresponding to weak pinning which suggests less disorder. Assuming that the superconducting hysteresis is symmetric about the equilibrium magnetization curve and that the magnetic hysteresis of the Eu spins is negligible as was found in the nonsuperconducting parent compound EuFe As 2 2 43 we determine the magnetization curve M(H) of the Eu spin sublattice as M = ( M + + M -) / 2 where M + ( M -) is the magnetization measured in increasing (decreasing) applied field. The results are shown for multiple temperatures in the main panels of Figs. 12 (a)-(f). Above T m , a Brillouin-like response is observed. The slopes of the magnetization curves grow as T m is approached, then become slightly lower as T decreases and the magnetic lattice stiffens. For H ‖ [110] and H ‖ [001], the satu-

M /

FIG. 13. Comparison of the symmetrized component, representing the ferromagnetic component, of the magnetic hysteresis of Ni-doped single crystals of RbEu(Fe 1 -x Ni ) As x 4 4 with H ‖ [110] (circles) and H ‖ [001] (triangles). Increasing the doping does not change the field necessary to collapse the Eu magnetic lattice.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00773v1-figure-41c3f994-49a0-4cb2-aad6-ff44012b86a8.jpg)

ration magnetization reaches ≈ 300 emu/cm 3 (compare table I) corresponding to 6.4 µ B /Eu, comparable to the expected full moment. Some excess moment was found for the x = 0 03 sample, that can again be attributed to . an impurity EuFe As 2 2 phase. Additionally, differences from the free moment value may arise from the approximation of the hysteresis model.

The 5 K Eu-sublattice magnetization curves for the three doping levels and a pristine sample in both field orientations are plotted together in Fig. 13. In order to minimize the effect of variations from varying degrees of hysteresis and from uncertainties in sample volume, we normalize each curve to its value at 6 T. There is no systematic change with doping of the saturation field in either orientation, with saturation fields of H ab sat ≈ 0 21T . and H c sat ≈ 0 8-1T for all doping levels, consistent with . the pronounced easy-plane magnetic anisotropy.

## IV. CONCLUSION

We have studied the effect of up to 4% Ni doping on the Fe-site in the magnetic superconductor RbEuFe As 4 4 in specific heat, resistivity and magnetization measurements. We observe a clear suppression of the superconducting transition temperature T c . Upon Ni-doping, the resistivity curve shifts up in a parallel fashion indicating a strong increase of the residual resistivity due to scattering by charged dopant atoms while the electronic structure appears largely unchanged. The observed step ∆ C/T c at the superconducting transition decreases strongly for increasing Ni doping in agreement with BNC scaling. The upper critical field slopes are reduced upon Ni doping for inas well as out-of-plane fields leading to a small reduction in the superconducting anisotropy. Specific heat measurements of the magnetic transition reveal the same BKT behavior close to the transition temperature T m for all doping levels with the transition temperature essentially unchanged and the in- to out-of-plane anisotropy of Eu-magnetism observed at small magnetic fields is unaffected. All of these observations indicate a decoupling of the Eu magnetism from superconductivity and essentially no influence of Ni doping on the Eu magnetism in this compound. These findings are in line with recent first principles calculations 36 indicating that the coupling between the Eu-moments is mediated by the RKKY interaction via the Fe-d z 2 orbitals which were found to be remarkably insensitive to Ni-doping.

- 1 Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. Am. Chem. Soc. 130 , 3296 (2008).
- 2 M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 101 (2008), 10.1103/PhysRevLett.101.107006.
- 3 J. W. Lynn and P. Dai, Physica C: Superconductivity 469 , 469 (2009).
- 4 D. C. Johnston, Advances in Physics 59 , 803 (2010).
- 5 G. R. Stewart, Rev. Mod. Phys. 83 , 1589 (2011).
- 6 K. Tanabe and H. Hosono, Japanese Journal of Applied Physics 51 , 010005 (2011).
- 7 H. Hosono and K. Kuroki, Physica C 514 , 399 (2015), superconducting Materials: Conventional, Unconventional and Undetermined.
- 8 Q. Si, R. Yu, and E. Abrahams, Japanese Journal of Applied Physics 1 , 13017 (2016).
- 9 Y. Bang and G. R. Stewart, Journal of Physics: Condensed Matter 29 , 123003 (2017).
- 10 N. Kurita, M. Kimata, K. Kodama, A. Harada, M. Tomita, H. S. Suzuki, T. Matsumoto, K. Murata, S. Uji, and T. Terashima, Phys. Rev. B 83 , 214513 (2011).
- 11 F. Hardy, N. J. Hillier, C. Meingast, D. Colson, Y. Li, N. Bariˇ si´, G. Yu, X. Zhao, M. Greven, c and J. S. Schilling, Phys. Rev. Lett. 105 , 167002 (2010).
- 12 C. Meingast, F. Hardy, R. Heid, P. Adelmann, A. B¨hmer, o P. Burger, D. Ernst, R. Fromknecht, P. Schweiss, and T. Wolf, Phys. Rev. Lett. 108 , 177004 (2012).
- 13 F. Hardy, A. E. B¨hmer, L. de' Medici, M. Capone, G. Gioo vannetti, R. Eder, L. Wang, M. He, T. Wolf, P. Schweiss, R. Heid, A. Herbig, P. Adelmann, R. A. Fisher, and C. Meingast, Phys. Rev. B 94 , 205113 (2016).
- 14 W. T. Jin, Y. Xiao, S. Nandi, S. Price, Y. Su, K. Schmalzl, W. Schmidt, T. Chatterji, A. Thamizhavel, and T. Br¨ckel, Phys. Rev. B u 100 , 014503 (2019).
- 15 A. Iyo, K. Kawashima, T. Kinjo, T. Nishio, S. Ishida, H. Fujihisa, Y. Gotoh, K. Kihou, H. Eisaki, and Y. Yoshida, Journal of the American Chemical Society 138 , 3410 (2016).
- 16 K. Kawashima, T. Kinjo, T. Nishio, S. Ishida, H. Fujihisa, Y. Gotoh, K. Kihou, H. Eisaki, Y. Yoshida, and A. Iyo, Journal of the Physical Society of Japan 85 , 064710 (2016).
- 17 Y. Liu, Y.-B. Liu, Q. Chen, Z.-T. Tang, W.-H. Jiao, Q. Tao, Z.-A. Xu, and G.-H. Cao, Science Bulletin 61 , 1213 (2016).
- 18 Y. Liu, Y.-B. Liu, Z.-T. Tang, H. Jiang, Z.-C. Wang, A. Ablimit, W.-H. Jiao, Q. Tao, C.-M. Feng, Z.-A. Xu, and G.-H. Cao, Phys. Rev. B 93 , 214503 (2016).
- 19 W. R. Meier, T. Kong, U. S. Kaluarachchi, V. Taufour, N. H. Jo, G. Drachuck, A. E. B¨hmer, S. M. Saunders, o A. Sapkota, A. Kreyssig, M. A. Tanatar, R. Prozorov, A. I. Goldman, F. F. Balakirev, A. Gurevich, S. L. Bud'ko, and P. C. Canfield, Phys. Rev. B 94 , 064501 (2016).

## V. ACKNOWLEDGMENT

This work was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division. K. W. acknowledges support from the Swiss National Science Foundation through the Early Postdoc.Mobility program.

- 20 W. Liu, L. Cao, S. Zhu, L. Kong, G. Wang, M. Papaj, P. Zhang, Y. Liu, H. Chen, G. Li, F. Yang, T. Kondo, S. Du, G. Cao, S. Shin, L. Fu, Z. Yin, H.-J. Gao, and H. Ding, arXiv:1907.00904 (2019).
- 21 H. S. Jeevan, D. Kasinathan, H. Rosner, and P. Gegenwart, Physical Review B 83 , 054511 (2011).
- 22 J. Koo, J. Park, S. Kook Cho, K. Duk Kim, S.-Y. Park, Y. Hee Jeong, Y. Jun Park, T. Yeong Koo, K.-P. Hong, C.-H. Lee, J.-Y. Kim, B.-K. Cho, K. Bong Lee, and H.-J. Kim, Journal of the Physical Society of Japan 79 , 114708 (2010).
- 23 S. Zapf and M. Dressel, Reports on Progress in Physics 80 , 016501 (2017).
- 24 S. Nandi, W. T. Jin, Y. Xiao, Y. Su, S. Price, D. K. Shukla, J. Strempfer, H. S. Jeevan, P. Gegenwart, and T. Br¨ckel, u Phys. Rev. B 89 , 014512 (2014).
- 25 Y. Xiao, Y. Su, M. Meven, R. Mittal, C. M. N. Kumar, T. Chatterji, S. Price, J. Persson, N. Kumar, S. K. Dhar, A. Thamizhavel, and T. Brueckel, Phys. Rev. B 80 , 174424 (2009).
- 26 J.-K. Bao, K. Willa, M. P. Smylie, H. Chen, U. Welp, D. Y. Chung, and M. G. Kanatzidis, Crystal Growth &amp; Design , Crystal Growth &amp; Design 18 , 3517 (2018).
- 27 M. P. Smylie, K. Willa, J.-K. Bao, K. Ryan, Z. Islam, H. Claus, Y. Simsek, Z. Diao, A. Rydh, A. E. Koshelev, W.-K. Kwok, D. Y. Chung, M. G. Kanatzidis, and U. Welp, Phys. Rev. B 98 , 104503 (2018).
- 28 K. Willa, R. Willa, J.-K. Bao, A. E. Koshelev, D. Y. Chung, M. G. Kanatzidis, W.-K. Kwok, and U. Welp, Phys. Rev. B 99 (2019), 10.1103/PhysRevB.99.180502.
- 29 M. P. Smylie, K. Willa, J.-K. Bao, K. Ryan, Z. Islam, H. Claus, Y. Simsek, Z. Diao, A. Rydh, A. E. Koshelev, W.-K. Kwok, D. Y. Chung, M. G. Kanatzidis, and U. Welp, Phys. Rev. B 98 , 104503 (2018).
- 30 K. Iida, Y. Nagai, S. Ishida, M. Ishikado, N. Murai, A. D. Christianson, H. Yoshida, Y. Inamura, H. Nakamura, A. Nakao, K. Munakata, D. Kagerbauer, M. Eisterer, K. Kawashima, Y. Yoshida, H. Eisaki, and A. Iyo, Phys. Rev. B 100 , 014506 (2019).
- 31 Y. Liu, Y.-B. Liu, Y.-L. Yu, Q. Tao, C.-M. Feng, and G.-H. Cao, Phys. Rev. B 96 , 224510 (2017).
- 32 K. Kawashima, S. Ishida, K. Oka, H. Kito, N. Takeshita, H. Fujihisa, Y. Gotoh, K. Kihou, H. Eisaki, Y. Yoshida, and A. Iyo, Journal of Physics: Conference Series 969 , 012027 (2018).
- 33 D. E. Jackson, D. VanGennep, W. Bi, D. Zhang, P. Materne, Y. Liu, G.-H. Cao, S. T. Weir, Y. K. Vohra, and J. J. Hamlin, Phys. Rev. B 98 , 014518 (2018).
- 34 S. Tagliati, V. M. Krasnov, and A. Rydh, Rev. Sci. Instrum. 83 , 055107 (2012).
- 35 K. Willa, Z. Diao, D. Campanini, U. Welp, R. Divan, M. Hudl, Z. Islam, W.-K. Kwok, and A. Rydh, Review

- of Scientific Instruments 88 , 125108 (2017).
- 36 C. Xu, Q. Chen, and C. Cao, Communications Physics 2 , 2399 (2019).
- 37 U. B. Paramanik, D. Das, R. Prasad, and Z. Hossain, Journal of Physics: Condensed Matter 25 , 265701 (2013).
- 38 A. E. Koshelev, K. Willa, R. Willa, M. Smylie, J.-K. Bao, D. Y. Chung, M. G. Kanatzidis, W.-K. Kwok, and U. Welp, (2019), arXiv:1906.10236.
- 39 P. J. Hirschfeld, M. M. Korshunov, and I. I. Mazin, Rep. Prog. Phys. 74 , 124508 (2011).
- 40 S. L. Bud'ko, N. Ni, and P. C. Canfield, Phys. Rev. B 79 , 220516 (2009).
- 41 Y. Bang and G. R. Stewart, New Journal of Physics 18 , 023017 (2016).
- 42 K. Willa, R. Willa, K. W. Song, G. D. Gu, J. A. Schneeloch, R. Zhong, A. E. Koshelev, W.-K. Kwok, and U. Welp, Phys. Rev. B 98 (2018), 10.1103/PhysRevB.98.184509.
- 43 S. Jiang, Y. Luo, Z. Ren, Z. Zhu, C. Wang, X. Xu, Q. Tao, G. Cao, and Z. Xu, New Journal of Physics 11 , 025007 (2009).
K. Willa M. P. Smylie Y. Simsek J. -K. Bao D. Y. Chung M. G. Kanatzidis W. -K. Kwok U. Welp
Magnetic and superconducting anisotropy in Ni-doped RbEuFe$_4$As$_4$ single crystals
We investigate the effect of Ni doping on the Fe-site in single crystals of
the magnetic superconductor RbEuFe$_4$As$_4$ for doping concentrations of up to
4%. A clear suppression in the superconducting transition temperature is
observed in specific heat, resistivity and magnetization measurements. Upon
Ni-doping, the resistivity curves shift up in a parallel fashion indicating a
strong increase of the residual resistivity due to scattering by charged dopand
atoms while the shape of the curve and thus the electronic structure appears
largely unchanged. The observed step $\Delta C/T_c$ at the superconducting
transition decreases strongly for increasing Ni doping in agreement with
expectations based on a model of multi-band superconductivity and strong
inter-band pairing. The upper critical field slopes are reduced upon Ni doping
for in- as well as out-of-plane fields leading to a small reduction in the
superconducting anisotropy. The specific heat measurements of the magnetic
transition reveal the same BKT behavior close to the transition temperature
$T_m$ for all doping levels. The transition temperature is essentially
unchanged upon doping. The in to out-of-plane anisotropy of Eu-magnetism
observed at small magnetic fields is unaltered as compared to the undoped
compound. All of these observations indicate a decoupling of the Eu magnetism
from superconductivity and essentially no influence of Ni doping on the Eu
magnetism in this compound.