Isotropic Pauli-Limited Superconductivity in the Infinite Layer Nickelate Nd 0.775 Sr 0.225 NiO 2
Bai Yang Wang, Danfeng Li, Berit H. Goodge, Kyuho Lee, Motoki Osada, Shannon P. Harvey, Lena F. Kourkoutis, Malcolm R. Beasley, Harold Y. Hwang
IIsotropic Pauli-Limited Superconductivity in the Infinite-Layer Nickelate Nd Sr NiO Bai Yang Wang,*
Danfeng Li,
Berit H. Goodge, Kyuho Lee,
Motoki Osada , Shannon P. Harvey , Lena F. Kourkoutis,
Malcolm R. Beasley , and Harold Y. Hwang* Department of Physics, Stanford University, Stanford, CA 94305, United States. Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory, Menlo Park, CA 94025, United States. Department of Applied Physics, Stanford University, Stanford, CA 94305, United States. School of Applied and Engineering Physics, Cornell University, Ithaca, NY 14853, United States. Department of Materials Science and Engineering, Stanford University, Stanford, CA 94305, United States. Kavli Institute at Cornell for Nanoscale Science, Cornell University, Ithaca, NY 14853, United States. *e-mail: [email protected]; [email protected].
Abstract
The recent observation of superconductivity in thin film infinite-layer nickelates offers fresh terrain for investigating superconductivity in layered oxides . A wide range of perspectives , emphasizing single- or multi-orbital electronic structure, Kondo or Hund’s coupling, and analogies to cuprates, have been theoretically proposed. Clearly further experimental characterization of the superconducting state is needed to develop a foundational understanding of the nickelates. Here we use magnetotransport measurements to probe the superconducting anisotropy in Nd Sr NiO . We find that the upper critical field is surprisingly isotropic at low temperatures, despite the layered crystal structure. We deduce that under magnetic field, superconductivity is strongly Pauli-limited, such that the paramagnetic effect dominates over orbital de-pairing. Underlying this isotropic response is a substantial anisotropy in the uperconducting coherence length, at least four times longer in-plane than out-of-plane. A prominent low-temperature upturn in the upper critical field indicates the likelihood of an unconventional ground state. Main
A key characteristic of layered superconductors is an upper critical field H c2 often exhibiting substantial anisotropy g º H c2,|| / H c2, ^ , where H c2,|| and H c2, ^ indicate magnetic fields in the a-b plane and along the crystal c -direction, respectively. In layered oxides such as cuprates , ruthenates , and cobaltates , the anisotropy reflects the two-dimensional (2D) nature of both the crystal and electronic structure. The case for the recently discovered infinite-layer nickelates is a priori interesting: although it shares the same crystal structure as infinite-layer cuprates, calculations indicate both a quasi-2D hole band and three-dimensional (3D) electron bands are present in the nickelates for a broad range of electron interactions . The relative importance of these bands is a subject of much interest with respect to the superconducting state. More generally, nickelates provide a new system within which to explore longstanding debates on the role of dimensionality for superconductivity. With these motivations, we have studied the superconducting anisotropy through resistive measurements of H c2 in thin-film Nd Sr NiO . This composition has the highest crystallinity in our recent investigations, as well as a superconducting transition temperature T c among the peak values observed across the superconducting dome . As previously reported, our optimal samples are capped with SrTiO , and are bounded in thickness to below ~10 nm; above this we observe mixed-phase formation or decomposition. A high-angle annular dark-field canning transmission electron microscopy (HAADF-STEM) image is shown in Fig. 1a. In good agreement with prior observations, c -axis oriented crystalline coherence is observed throughout the film, together with extended defects primarily based on vertical Ruddlesden-Popper (RP) type faults. Figure 1b, c shows the temperature-dependent resistivity r ( T ) in the low temperature region, with varying magnetic field applied along the c -axis ( H ^ ) and the a-b plane ( H || , transverse geometry). While minimal magnetoresistance is observed in the normal state, the superconducting transition is suppressed with increasing magnetic field, allowing us to extract the temperature dependence of H c2 . Figure 2 shows the main result of this study: the magnetic field-temperature phase diagram of H c2, ^ and H c2,|| , down to 0.5 K. Here T c and H c2 are identified as the resistive transition reaching 50% of the normal state, as discussed below. Overall, the most striking feature is that H c2 is largely isotropic at low temperatures (as emphasized in Fig. 2a inset), despite the layered crystal structure, the quasi-2D nature of the ! " $% hole band and the thin film geometry. Does this reflect a rather isotropic electronic structure, and superconducting coherence length, as unusually found in (K,Ba)FeAs (ref. )? Or does it have a different fundamental origin? Our aim is to address this central question. As we argue below, H c2 is dominated by the isotropic spin response, rather than orbital effects, over much of the phase diagram. Such strongly Pauli-limited superconductivity in nickelates appears unique among the layered oxides. We start by examining the magnetoresistive transition observed in Fig. 1b, c. The visible broadening can have contributions from static compositional inhomogeneity, superconducting fluctuations, and vortex motion in this type-II superconductor . In terms of inhomogeneity, the P-type faults likely give rise to a spread in T c . Regarding superconducting fluctuations, we fit the onset of superconductivity following the 2D Aslamazov-Larkin (AL) formulism , with representative fits shown in the insets of Fig. 1b, c. From this we obtain a corresponding T cAL , above which superconducting fluctuations significantly impact the resistive transition. We note that this regime also includes effects arising from static inhomogeneities. As for vortex motion, we investigate its contribution in scaled Arrhenius plots, as shown in Fig. 3a, b. Between the normal state and the measurement noise floor, exponentially decreasing resistivity is seen for both orientations across several orders of magnitude. This is well described by thermally-activated vortex motion : r ( T , H ) = r ( H ) exp( U / k B T ) where r is the resistivity prefactor, U is the activation energy and k B is Boltzmann’s constant. Figure 3c shows U extracted using this form. Similar to cuprates , U for H ^ and H || exhibits regions of power-law dependence on magnetic field. Note that the success of this fitting implies that vortex motion is relevant for H || , and the large difference in U values between H ^ and H || indicates substantial vortex-motion anisotropy. We define T cVM as the temperature above which resistivity deviates from thermally-activated behavior, and attribute the resistive transition below T cVM predominantly to vortex motion. Summarizing the onset of fluctuations and vortex motion, we construct the phase diagrams in Fig. 4a, b. In addition to H cVM ( T ) and H cAL ( T ), we plot representative candidate proxies for H c2 : H c90% ( T ), H c50% ( T ) and H c1% ( T ), corresponding to the resistivity reaching 90%, 50%, and 1% of the normal state. Between H cVM ( T ) and H cAL ( T ) (extrapolating to 27% and 80% of the resistive transition for H = 0 T) there is a regime where the magnetoresistance can serve as a reliable ndicator of H c2 , free from vortex creep or superconducting fluctuations which tend to under/over-estimate H c2 , as well as distort its functional form . Hence, we use H c50% to measure H c2 ; our conclusions are robust to variations within this regime aside from specific numerical values. We next focus on magnetotransport near T c , taking advantage of the simplification in the critical regime arising from diverging length scales. Here the central observation is a T -linear dependence of H c2, ^ and a T dependence of H c2,|| (Fig. 4c). We first consider purely orbital de-pairing (superconductivity suppressed by the energy cost of the field-induced supercurrent), with the characteristic T dependence of H c2,|| arising from 2D confinement. This is commonly observed for superconductors thinner than the Ginzburg-Landau (GL) coherence length , and is a reasonable possibility here, given the thickness limits for stabilizing Nd Sr NiO . In such geometries, H c2 is well described by the linearized GL form (generalized to an anisotropic, c -axis thickness-limited superconductor) as & '(,* (-) = 0 (0) 81 − -- ' ;, (1) & '(,∥ (-) = √120 (0)! 81 − -- ' ; >( , (2) where f is the flux quantum, x ab (0) is the in-plane zero-temperature coherence length, and d is the superconducting thickness. Using these forms, we extract x ab (0) = 42.6 ± Å and d = 225 ± Å . In previous studies on other material systems, d derived from such analysis was in good quantitative agreement with the physical thickness of the film , except for cases with dominant spin-orbit scattering , for hich the thickness is underestimated. Here we deduce a thickness that is 2.3 – 3 times larger than the observed thickness of the film (Table S1, Supplementary Information), which we take to be an unphysical discrepancy. Furthermore, this analysis presumes confinement preventing vortex entry for H || , inconsistent with evident vortex motion in this orientation (Fig. 3c). Finally, 2D confinement of an orbital-limited superconductor usually leads to a dramatic enhancement of H c2,|| over H c2, ^ that is inversely proportional to thickness, unlike the response observed here (insets to Fig. 2a, b). These observations strongly indicate a different origin for the behavior of H c2,|| . In fact, if the superconductor is Pauli-limited, the paramagnetic effect can also give a T dependence near T c (refs ), arising from the Zeeman energy gain of the normal state (~ µ H ) competing against the superconducting condensate energy (~ H /8 p ). In this case, fitting the T dependence allows us to extract the electronic effective magnetic moment near T c to be 2.37 ± . Additionally, as the purely orbital-limited H c2,|| necessarily exceeds the measured (Pauli-limited) H c2,|| , the observation of vortex motion below T cVM provides only an upper bound estimate of x c (0) £ Å , or 3.3 NiO layers (Supplementary Information), corresponding to x c (0)/ x ab (0) £ To quantitatively describe both orientations and lower temperatures beyond the critical regime, both orbital and paramagnetic effects clearly need to be considered, as captured in Werthamer–Helfand–Hohenberg (WHH) theory . Recognizing the unusual nature of the prominent upturn at lower temperatures, we focus first on T > 4 K and fit the data (independently for each orientation) o the single-band WWH form, constrained to a second-order superconducting transition for all temperatures. We use single-band WHH given that no signatures of multi-band superconductivity were observed near T c (ref. ). The three fit parameters allowed to vary (Table S1, Supplementary Information) are the Maki parameter a (characterizing the ratio between the orbital-limited and Pauli-limited upper critical fields), the spin-orbit scattering parameter l so , and the effective mass m eff . A number of consistency checks (Supplementary Information) confirm that: the spin-orbit scattering time is much longer than the elastic-transport scattering time (assumed in WHH); a reasonable electronic diffusion coefficient is deduced; x ab (0) is sensibly related to the clean-limit coherence length and the transport mean free path; and a is consistent with estimates from the critical regime near T c . Such overall consistency and the high quality of the fits allow us to draw some qualitative conclusions robust to the fitting details. The most prominent feature is the high values of a for both orientations, indicating the strong presence of the paramagnetic effect. In particular, the exceptionally large a || supports the purely Pauli-limited scenario for H c2,|| . Such an interpretation also explains the decreasing anisotropy of H c2 , as the apparently isotropic paramagnetic effect plays an increasingly important role approaching 4 K. This transition is visible through the substantial suppression of H c2, ^ below the purely orbital-limited response with decreasing temperature, as illustrated via the pair-breaking picture (Supplementary Information). The large a values also indicate that superconductivity here is in the dirty limit, in agreement with transport estimates . Finally, we note that the enhanced susceptibility independently deduced from Fig. 4c ( H c2,|| ) is consistent with the scale for the Pauli limit relevant for the WHH fit (with enhancements associated with l so ). Below 4 K, the anomalous upturns of H c2, ^ and H c2,|| require considerations beyond conventional WHH theory . Similar upturns, of much debated origin, have been observed for a number of layered superconductors . We are not in a position to suggest a most likely scenario, but discuss a few briefly here. The behavior is tantalizingly similar to expectations for the Fulde-Ferrell-Larkin-Ovchinnikov state ; however, the fact that we are in the dirty limit is inconsistent with anticipated requirements. The upturn could indicate the onset of two-band superconductivity , but the second-band contribution would also be expected to saturate to an isotropic, Pauli-limited value. Finally, we discuss the possible onset of a first-order superconducting transition . This tendency is expected for a strongly Pauli-limited superconductor at low temperatures, although with different curvature . If we remove the second-order restriction to our WHH fit, a comparable description of the data can be obtained in the high-temperature region, although with substantially smaller l so . Therefore, l so should be considered under-constrained by our present data, whereas the other conclusions we draw from WHH are robust to the possibility of a first-order transition, and quite reproducible (see Supplementary Information). In summary, our studies of superconductivity in Nd Sr NiO indicate that H c2 is largely isotropic at low temperatures because it is strongly Pauli-limited, although x c (0)/ x ab (0) £ H c2,|| extrapolates to the bare Pauli limit, uggesting that the anomalous low-temperature upturn and enhanced susceptibility may have a common origin. ethods Thin film synthesis.
The samples were grown on SrTiO substrates using pulsed-laser deposition under growth and reduction conditions previously reported . HAADF-STEM characterization.
Electron transparent STEM samples were prepared on a Thermo Fischer Scientific Helios G4 X focused ion beam (FIB) using the standard lift-out method, as was previously reported . Samples were thinned to <30 nm with 2 kV Ga ions, followed by a final polish at 1 kV to reduce effects of surface damage. HAADF-STEM characterization of the sample was then obtained on an aberration-correction FEI Titan Themis at an accelerating voltage of 300 kV with a convergence angle of 30 mrad and inner collection angle of 68 mrad.
Transport measurements.
The samples were contacted using wire-bonded aluminum wires in a 6-point Hall bar geometry. The 50% resistive criterion is applied to both temperature-sweep and field-sweep data to determine T c / H c2 . eferences
1. Li, D. et al.
Superconductivity in an infinite-layer nickelate.
Nature et al.
Superconducting dome in Nd – x Sr x NiO infinite layer films. Preprint at http://arxiv.org/abs/2003.08506 (2020). 3. Zeng, S. et al. Phase diagram and superconducting dome of infinite-layer Nd – x Sr x NiO thin films. Preprint at http://arxiv.org/abs/2004.11281 (2020). 4. Norman, M. R. Entering the nickel age of superconductivity. Physics
85 (2020). 5. Lee, K. W. & Pickett, W. E. Infinite-layer LaNiO : Ni is not Cu . Phys. Rev. B and CaCuO and implications for superconductivity. Phys. Rev. X et al.
Robust d " $% -wave superconductivity of infinite-layer nickelates. Phys. Rev. B . Phys. Rev. Lett. T c in superconducting nickelates. Phys. Rev. Res. Phys. Rev. B et al.
How to make superconducting-anisotropy least in high- T c cuprate superconductors. Phys. C RuO . Phys. C et al.
Thermodynamic and transport measurements of superconducting Na
CoO ·1.3H O single crystals prepared by electrochemical deintercalation.
Phys. Rev. Lett. et al.
Nearly isotropic superconductivity in (Ba,K)Fe As . Nature
Sov. Phys. Solid State
Z. Phys.
Phys. Rev. B et al.
Flux-creep dissipation in epitaxial YBa Cu O d film: Magnetic-field and electrical-current dependence. Phys. Rev. B
Phys. Rep.
Phys. Rev. B et al. Two-dimensional normal-state quantum oscillations in a superconducting heterostructure.
Nature
Phys. Rev. δ -doped SrTiO heterostructures. Phys. Rev. B
Phys. Rev. Lett. /LaAlO interface: a magnetotransport study. Phys. Rev. Lett.
Appl. Phys. Lett. Phys. Rev. Lett. Prog. Theor. Phys.
Phys. Rev. Lett. A . Phys. Rev. Lett. H c2 . III. Electron spin and spin-orbit effects. Phys. Rev.
Phys. Rev. B
J. Phys. Soc. Jpn. et al.
Calorimetric evidence for a Fulde-Ferrell-Larkin-Ovchinnikov superconducting state in the layered organic superconductor κ-(BEDT-TTF) Cu(NCS) . Phys. Rev. Lett.
Phys. Lett. et al.
Aspects of the synthesis of thin film superconducting infinite-layer nickelates.
APL Mater. et al. Doping evolution of the Mott-Hubbard landscape in infinite-layer nickelates. Preprint at http://arxiv.org/abs/2005.02847 (2020).
Acknowledgements
We thank R. L. Greene, A. Kapitulnik, S. A. Kivelson, P. B. Littlewood, B. Maiorov, G. A. Sawatzky, H. Takagi, R. Thomale, A. Viswanathan, and Y.-H. Zhang for discussions. The work at SLAC and Stanford was supported by the US Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering, under contract number DE-AC02-76SF00515, and the Gordon and Betty Moore Foundation’s Emergent Phenomena in Quantum Systems Initiative through grant number GBMF4415 (synthesis equipment). B.H.G. and L.F.K. acknowledge support by the Department of Defense Air Force Office of Scientific Research (No. FA 9550-16-1-0305). This work made use of the Cornell Center for Materials Research (CCMR) Shared Facilities, which are supported through the NSF MRSEC Program (No. DMR-1719875). he FEI Titan Themis 300 was acquired through No. NSF-MRI-1429155, with additional support from Cornell University, the Weill Institute, and the Kavli Institute at Cornell. The Thermo Fisher Helios G4 X FIB was acquired with support from the National Science Foundation Platform for Accelerated Realization, Analysis, and Discovery of Interface Materials (PARADIM) under Cooperative Agreement No. DMR-1539918.
Author Contributions
D.L., K.L., and M.O. grew the nickelate films, and conducted the reduction experiments and structural characterization. B.H.G. and L.F.K. conducted electron microscopy. B.Y.W. performed the transport measurements, and analysis with M.R.B. B.Y.W., S.P.H., and H.Y.H. wrote the manuscript with input from all authors.
Figure 1: Structure and magnetotransport properties of thin film Nd Sr NiO . a, Cross-sectional HAADF-STEM image of a Nd Sr NiO thin film. b, c, r ( T ) below 40 K at H = 0, 0.2, 0.4, 0.6, 0.8, 1, 2, 3, 4, 6, 8, 10, 12 T for magnetic field along the c -axis ( b ) and transverse in the a-b plane ( c ). The insets show representative resistive transitions with 2D Aslamazov-Larkin fits shown as dashed lines for field along the c -axis ( b ) and in the a-b plane ( c ). Figure 2: H c2 versus T phase diagram for Nd Sr NiO . a, b, H c2 – T phase diagrams for magnetic field along the c -axis ( a ) and in the a-b plane ( b ). The solid green lines are Werthamer–Helfand–Hohenberg (WHH) fits of the H c2 data above the low temperature upturn, constrained to 2 nd order superconducting transitions at all temperatures. The solid red line ( a ) is a 3D Ginzburg-Landau fit of the H c2, ^ data near T c . The black arrows indicate the free-electron Pauli-limited field H Pauli, µ =1 µ B = 1.86 ´ T c ( H = 0 T). The red arrows indicate the effective Pauli-limited field H Pauli, µ =2.4 µ B , corresponding to the enhanced magnetic susceptibility deduced near T c . The inset of ( a ) overlays the H c2 data and the corresponding WHH fits for both field orientations. The inset of ( b ) overlays the H c2,|| data of two samples with different thickness. The solid lines are guides to the eye. Figure 3: Transport signatures of vortex motion, and the magnetic field dependence of the activation energy. a, b, r ( T ) in scaled Arrhenius plots at field H = 0.2, 0.4, 0.6, 0.8, 1, 2, 3, 4, 6, 8, 10, 12 T for magnetic field along the c -axis ( a ) and in the a-b plane ( b ). The x -axes are scaled by the activation energies U obtained through fits of the resistive transition in the thermally activated dissipation region according to: r ( T , H ) = r ( H ) exp( U / k B T ). c, The extracted U as a function of magnetic field in double logarithmic scale for both field orientations. The solid lines are power law fits of U below 2 K for field along the c -axis and in the a-b plane, while the dashed lines are power law fits of U above 5 K. Figure 4: H c versus T phase diagram and temperature dependence of H c2 near T c . a, b, H c – T phase diagrams for magnetic field along the c -axis ( a ) and in the a-b plane ( b ). The regions separated by H c1% , H cVM , H cAL , and H c90% indicate superconductivity (SC), vortex motion, superconducting fluctuations, and the normal state. c, Temperature dependence of H c2 (for which we use H c50% ) and ( H c2 ) near T c for field along the c -axis and in the a-b plane, respectively. The dashed lines are linear fits of H c2 and ( H c2 ) above 85% of T / T c ( H = 0 T). These fits are also shown in ( a ) and ( b ). upplementary Information for “Isotropic Pauli-Limited Superconductivity in the Infinite Layer Nickelate Nd Sr NiO ” Bai Yang Wang,*
Danfeng Li,
Berit H. Goodge, Kyuho Lee,
Motoki Osada , Shannon P. Harvey , Lena F. Kourkoutis,
Malcolm R. Beasley , and Harold Y. Hwang* Department of Physics, Stanford University, Stanford, CA 94305, United States. Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory, Menlo Park, CA 94025, United States. Department of Applied Physics, Stanford University, Stanford, CA 94305, United States. School of Applied and Engineering Physics, Cornell University, Ithaca, NY 14853, United States. Department of Materials Science and Engineering, Stanford University, Stanford, CA 94305, United States. Kavli Institute at Cornell for Nanoscale Science, Cornell University, Ithaca, NY 14853, United States. *e-mail: [email protected]; [email protected]. . Reproducibility
The measurement results and conclusions drawn in the main text are reproducibly observed. We show in Fig. S1 the temperature dependent resistivity of a second sample, Sample H c2 versus T phase diagrams of Sample T behavior of H c2,|| near T c , a decreasing anisotropy of H c2 at lower temperatures (shown also in the inset of Fig. S1a), and anomalous upturns of H c2 approaching 0 K. A summary of sample and fitting parameters is given in Table S1. FIG. S1. r ( T ) for Sample H = 0, 0.5, 1, 2, 3, 4, 5, 6, 7, 8, 10, 12, 14 T for magnetic field along the c -axis ( a ) and the a-b plane ( b ). FIG. S2. a, b, H c2 – T phase diagrams of Sample c -axis ( a ) and the a-b plane ( b ). The solid green lines are WHH fits of the H c2 data above the low temperature upturn. The solid red line ( a ) is a 3D Ginzburg-Landau fit of the H c2, ^ data near T c . The black arrows correspond to the free-electron Pauli-limited field H Pauli, µ =1 µ B = 1.86 ´ T c ( H = 0 T). The red arrows correspond to the effective Pauli-limited field H Pauli, µ =2.2 µ B , corresponding to the enhanced magnetic susceptibility deduced near T c . The inset of ( a ) overlays the H c2 data and the corresponding WHH fits for both field orientations. I. Estimation of the Upper Bound of x c (0) Considering the 3D GL anisotropic formula : ! !",∥ = % &' (0)& ! (0) *1 − -- ! . an upper bound for x c (0) can be provided by the observation of vortex motion dissipation for H || . While we expect a relatively uniform superconducting order parameter for a thickness-limited ( d /1.8 < x c ) 2D superconductor with no vortex motion , we observe a substantial region of thermally activated dissipation below T cVM (Fig. S3, dashed orange line). Consequently, our observation requires x c to be in the 3D regime ( x c < d /1.8) below T cVM . This provides an upper bound of the zero-temperature c -axis coherence length: x c (0) £
10 Å (9 Å for Sample c -axis Ginzburg-Landau coherence length x c = x c (0) × (1 – T / T c ) –1/2 , for x c (0) ranging from 1 Å to 15 Å. The vertical dashed orange line indicates T cVM ( H = 0 T), below which thermally activated dissipation is observed for the H || orientation. The top horizontal dashed line indicates the 2D-3D thickness crossover for vortex motion, corresponding to d /1.8 where d is the sample thickness. The bottom horizontal dashed line indicates the c -lattice constant of the sample. II. Werthamer–Helfand–Hohenberg (WHH) Fitting
WHH fit parameter self-consistency checks:
Firstly, WHH theory relies on the spin-orbit scattering time being much longer than the elastic transport scattering time. To estimate the transport scattering rate, we performed separate Hall effect measurements of the same sample in H ^ orientation. For the undoped infinite layer nickelate, the broad consensus of electronic structure calculations is that both hole and electron bands are present in the system . As Sr doping is increasingly introduced, the electron pockets deplete. Although it is unclear when precisely the electron bands are fully depleted, we chose the highest superconducting doping level of 22.5%, where the low-temperature Hall coefficient is positive, to maximize the likelihood of a single hole-band scenario. Assuming simple Boltzmann transport and m eff , the obtained transport scattering time is t tr » ´ –16 s, which is indeed much shorter than the long spin-orbit scattering time t so » ´ –13 s derived from the fitted l so for the two field orientations, via l so = h /6 p k B T c t so , where h is Planck’s constant. Similarly, we can estimate the normal state in-plane electronic diffusion coefficient D || based on the transport scattering time and the Fermi velocity v f of the hole-band estimated from band calculations : 0.6 cm /s. This result agrees well with the fitted electronic diffusion coefficient 0.4 cm /s, derived from the fitted a ^ , according to a = 3 h /4 p m eff v f2 t tr . In addition, we note that the GL fit deduced in-plane coherence length x ab (0) ~ 42.6 ± 0.1 Å matches the geometric mean of the clean-limit coherence length x = 0.18 ´ hv f /2 p k B T c ~ 921 Å and the transport mean free path l mfp ~ 3.3 Å , as expected. Finally, the Maki parameter a can also be estimated from linear fits of H c2 near T c (ref. ): a || = 20.6 , a ^ = 1.1, in reasonable agreement with those obtained from the WHH fits (see Table S1). Based on the above estimates, it is worth noting that for our system, k f l mfp ~ 3.8 (where k f is the Fermi momentum ) and l mfp << x ab (0), reiterating that here superconductivity is in the dirty limit. (A similar k f l mfp ~ 4.2 results assuming a 2D free lectron gas.) This raises the possibility of localization effects on superconductivity , with regards to the anomalous upturn of H c2 . However, it is unclear how localization effects are manifested in a strongly Pauli-limited superconductor. nd order versus 1 st order WHH fits: Due to the unconventional nature of the anomalous upturns of H c2 at low temperatures , we limit our conventional WHH-theory-based fits to the high temperature regime. However, this leaves ambiguous, and experimentally under-constrained, the detailed WHH form at low temperatures. Therefore, we present two extremal cases here: 1) constraining the WHH fit to 2 nd order superconducting transitions at all temperatures (as in main text); and 2) allowing for 1 st order superconducting transitions to occur at low temperatures. The obtained WHH fits are shown in Fig. S4 as solid green and orange lines and the corresponding fit parameters are shown in Table S1. Although there is substantial variation in the fitted l so , which largely determines the order of the low-temperature phase transition, the magnitude of the obtained a and m eff are fairly independent of the fit choices, allowing robust qualitative conclusions to be drawn from these parameters. Note that the WHH form assumes a fixed g -factor of 2, relevant for comparison to the susceptibility enhancement deduced near T c . FIG. S4. a, b, H c2 – T phase diagrams of Sample c -axis ( a ) and the a-b plane ( b ). The solid green lines are WHH fits of the H c2 data above the low temperature upturn, constrained to 2 nd order superconducting transitions at all temperatures, and the orange lines are WHH fits allowing 1 st order transitions to occur at low temperatures. The dashed section of the orange line below the tricritical temperature T corresponds to the supercooling field where the superconducting transition becomes 1 st order . V. Paramagnetic versus Orbital Effects for H ^ Orientation The relative contributions from paramagnetic effects and orbital effects to the temperature dependence of the upper critical field can be identified using the pair-breaking form : ln( T / T c ) = y (1/2) – y (1/2 + b /2 p k B T c ), where y is the digamma function and b is the pair-breaking process energy. The b values corresponding to paramagnetic and orbital effects with perpendicular magnetic field are t so e hH /4 p m and D || eH / c , respectively. Here e is the electronic charge, c is the speed of light, and m is the free electron mass. The b values of the two effects can be added independently when considered together . As shown in Fig. S5, the pair-breaking formula considering only orbital effects (solid red line) fits H c2, ^ above 7 K satisfactorily. The extracted electronic diffusion coefficient D || » /s matches that obtained from both the WHH fit and the transport measurements discussed above. Below 7 K, however, H c2, ^ clearly deviates below the pair-breaking formula fit, suggesting an increasing role of the paramagnetic effect. Shown as a solid green line in Fig. S5, the pair-breaking formula fit considering both paramagnetic and orbital effects describes the temperature dependence of H c2, ^ accurately down to the lower temperature upturn. This supports the identification that H c2, ^ is orbital-limited near T c but crosses over to become Pauli-limited at lower temperatures. FIG. S5. Upper critical field H c2 versus T phase diagram of Sample H ^ orientation. Blue dots are measured H c2, ^ . Solid red and green lines are pair-breaking formula fits considering only orbital effects, and considering both orbital and paramagnetic effects, respectively. . TABLE S1. Parameters of the sample, superconducting phase, and WHH fits for Sample a as the superconductivity becomes strongly Pauli-limited ( a >> 1). We further note that m eff here is an effective parameter following the starting assumptions of WHH – i.e. of BCS superconductivity in a Fermi liquid with a spherical Fermi surface, thus not directly comparable to the band-structure or quasiparticle effective mass. T c (K) Sample thickness (Å) Extracted thickness d (Eq. 2) (Å) Extracted x ab (0) (Å) Effective moment ( µ B ) Maki parameter a (2 nd order/ 1 st order) Spin-orbit scattering parameter l so (2 nd order/ 1 st order) Effective mass m eff (2 nd order/ 1 st order) Sample H ^ H || H ^ H || eferences
1. Tinkham, M.
Introduction to superconductivity . (McGraw-Hill, Inc., 1996). 2. Maki, K. & Tsuneto, T. Pauli paramagnetism and superconducting state.
Prog. Theor. Phys.
Phys. Rev. B : Ni is not Cu . Phys. Rev. B and CaCuO and implications for superconductivity. Phys. Rev. X et al.
Robust d ( ! )* ! -wave superconductivity of infinite-layer nickelates. Phys. Rev. B
Phys. Rev. Res. et al. Model construction and a possibility of cuprate-like pairing in a new d . Preprint at https://arxiv.org/abs/1909.00060 (2019). 9. Nomura, Y. et al. Formation of a two-dimensional single-component correlated electron system and band engineering in the nickelate superconductor NdNiO . Phys. Rev. B . Phys. Rev. Lett. H c2 . III. Electron spin and spin-orbit effects. Phys. Rev.
Prog. Theor. Phys. Suppl.
Phys. Rev. B
Phys. Rep.
Phys. Rep.238,