A brief review on muSR studies of unconventional Fe and Cr-based superconductors

Muon spin relaxation/rotation (muSR) is a vital technique for probing the superconducting gap structure, pairing symmetry and time reversal symmetry breaking, enabling an understanding of the mechanisms behind the unconventional superconductivity of cuprates and Fe-based high-temperature superconductors, which remain a puzzle. Very recently double layered Fe-based super-conductors having quasi-2D crystal structures and Cr-based superconductors with a quasi-1D structure have drawn considerable attention. Here we present a brief review of the characteristics of a few selected Fe- and Cr-based superconducting materials and highlight some of the major outstanding problems, with an emphasis on the superconducting pairing symmetries of these materials. We focus on muSR studies of the newly discovered superconductors ACa2Fe4As4F2(A = K, Rb, and Cs), ThFeAsN, and A2Cr3As3(A = K, Cs), which were used to determine the superconducting gap structures, the presence of spin fluctuations, and to search for time reversal symmetry breaking in the superconducting states. We also briefly discuss the results of muSR investigations of the superconductivity in hole and electron doped BaFe2As2.


I. INTRODUCTION
In 1911, Kamerlingh Onnes measured the electrical conductivity of numerous metals and discovered the abrupt disappearance of the resistance of a solid mercury wire in liquid helium 1 . Subsequently a wide range of superconductors have been found with increasing high transition temperatures (T c ), and it is hoped that a roomtemperature superconductor can ultimately be realized. Following the discovery of superconductivity it took a long time for a microscopic theoretical understanding of the phenomenon. Such a microscopic theory (BCS theory) was proposed in 1957 by John Bardeen, Leon Cooper, and John Robert Schrieffer 2,3 , where electrons with opposite momenta condense into Cooper pairs, which are bound via an attractive interaction mediated by lattice vibrations (phonons). Superconductors which are adequately characterized by BCS theory, where the Cooper pairs are bound by the electron-phonon interaction, are known as conventional superconductors, which make up the vast majority of known superconducting materials to date. In these systems, the Cooper pairs condense in a s-wave pairing state with zero orbital angular momentum, which leads to a gap in the single particle excitation spectrum across the whole Fermi surface 4 . For phonon mediated BCS superconductivity, the transition temperature is comparatively low, with an expected maximum of 30-40 K with strong coupling (at ambient pressure) 2,3,5 . A dramatic breakthrough in the field of superconductivity came from the discovery of a La-Ba-Cu-O system with a T c of around 30 K by Bednorz and Müller 6 , for which they were awarded the Nobel Prize in Physics 1987. The cuprates materials are quite distinct from previously known superconductors, consisting of doped CuO 2 layers. Soon after, cuprates with significantly higher values of T c were found, including the first material to become superconducting above the boiling point of liquid nitrogen, YBa 2 Cu 3 O 7−δ , with a T c of about 93 K 7 . The largest T c reported at ambient pressure is at around 135 K in Hg 2 Ba 2 Ca 2 Cu 3 O 10+δ 8 , which can be increased even further under high applied pressures 9 .
The iron pnictide superconductors (FeSCs) were the second family of high−T c superconductors to be discovered, by Kamihara et al. in 2006 13 , which was even more exciting given the long standing notion that ferromagnetic elements such as iron are detrimental to superconductivity. Prior to this finding, superconductivity had already been discovered with T c = 4 K in LaOFeP 14 , but the subsequent study of La(O 1−x F x )FeAs revealed a clearly enhanced value of T c = 26 K 13 . It was further found that upon switching La with other rare earth elements and with either fluorine doping or introducing oxygen deficiencies, T c could be further increased 15,16 . Much like the cuprates, it was quickly apparent that these classes of materials are also unconventional superconductors, with a pairing mechanism not mediated by the electron-phonon interaction of BCS theory [17][18][19][20][21][22] . Here the essential component for superconductivity appears to be the FeAs layers, as opposed to the Cu-O planes of the cuprates. One similarity to the cuprates is the tunabil-arXiv:1811.12677v1 [cond-mat.supr-con] 30 Nov 2018 ity of the properties upon carrier doping, starting from a non-superconducting parent compound. In the cuprates this parent compound corresponds to a Mott insulator, whereas for the Fe-based superconductors these are generally metallic or semimetallic with a spin-density wave ground state 19,23,24 . The crystal structures typically consist of arrangements of layers of FeAs and spacer layers, and a number of varieties are found. The structures of the 11, 111, 122,1111 and 32522-type materials are shown in Figs. 1(a-f) 19 , along with that of the recently discovered ACa 2 Fe 4 As 4 F 2 (A = K, Rb and Cs, 122442-family) superconductors. One of the striking features found in some iron based superconductors is the observation of electronic nematicity below a structural transition, [25][26][27] the origin of which is not yet resolved. In the nematic state the four fold rotational symmetry of the high temperature tetragonal structure is broken, and this appears to be electronically driven, corresponding to either spin or orbital nematic ordering [28][29][30][31][32][33][34][35][36][37] . Whether the nematic order in all the Fe-based systems has a common origin is another open question, particularly in the case of FeSe where there appears to be a lack of magnetic order in the nematic state, which has been accounted for by orbital ordering originating from the Aslamazov-Larkin vertex correction 38 . In general, nematic order appears to compete with superconductivity 39 , although when the s-wave and d-wave states are nearly degenerate, it has been suggested that the nematic order can lift the frustration and hence lead to an increase of T c 26,40,41 .
Another example of unconventional superconductivity is in the heavy fermion superconductors 42 , such as CeCu 2 Si 2 43 , CeCoIn 5 44 , CePt 3 Si 45 , UPt 3 46 and URu 2 Si 2 47 , which in spite of the much smaller values of T c , show a number of similarities to the high temperature superconductors. These similarities can be seen more clearly in Figs. 2(a)-(c), which show the schematic phase diagrams of three of the most investigated classes of unconventional superconductors, including (a), a cuprate superconductor (YBCO) from Leyraud et al. 10 , (b), a heavy fermion superconductor (CeRhIn 5 ) from Park et al. 11 , and (c) Fe-based (doped Ba-122) superconductors from Liu et al. 12 A common feature is that the application of a non-thermal tuning parameter such as doping or pressure suppresses a magnetically ordered phase, but leads to a superconducting dome, which may encompass a quantum critical point near to where T c is largest. Meanwhile, a plot showing the historical progression of realized values of T c is shown in Fig. 2(d).
Cooper pairs can be described by a wave function, consisting of spatial and spin parts, corresponding to particular orbital anglar momentum L and spin S. Spin-singlet wave functions have S = 0 and L = 0, 2, ..., whereas the spin-triplet wave functions have S = 1 and L = 1, 3, .... One salient feature of unconventional pairing states is that the Cooper pairs do not only form in the s-wave state of BCS theory, but can condense in higher angular momentum states, and the gap can be highly anisotropic, disappearing at points or lines on the Fermi surface. sults are most consistent with d-wave symmetry 48,49 , and the Cooper pairs are thought to mediated by antiferromagnetic spin fluctuations 50 . The gap corresponding to such a d-wave state is highly anisotropic, and it changes sign and has line nodes on the Fermi surface, as shown in Fig. 3(b). In fact a sign changing superconducting order parameter is generally required for the realization of spinfluctuation mediated superconductivity 50 , and therefore in analogy with the cuprates, one may also have anticipated d-wave symmetry in iron-arsenide superconductors. However, experimental probes on iron-based materials have in most cases revealed nodeless superconductivity, and this has often been accounted for by an s ± −state 51 similarly mediated by antiferromagnetic spin fluctuations, where there is a different sign of the order parameter between disconnected Fermi surface pockets, as shown in Fig. 3(f). Similar to the cuprates, Fe-based superconductors have also been much studied with inelastic neutron scattering (INS). These studies have revealed a spin resonance peak below T c in the 1111, 122, 111, and 11 families of superconductors, where the energy of the resonance peak has a linear relationship with T c 52 . The presence of such a spin resonance suggests a change of sign of the superconducting pairing function between regions of the Fermi surface separated by the ordering wave vector of the antiferromagnetic parent compound. Therefore in many cases, the observation of a spin resonance is in good agreement with that expected for the s ± state, since the AFM wave vector often connects the hole pocket at the zone center, and the electron pockets at the zone edge 52 .
There has been considerable debate over to what ex-tent the superconductivity of iron pnictides can be universally described by an s ± state. This is particularly the case for overdoped 122 FeAs-based systems such as KFe 2 As 2 which lack electron pockets at the zone edge 53 , and FeSe-based materials such as K x Fe 2−y Se 2 54 and (Li 1−x Fe x )OHFeSe 55,56 , where the hole pockets at the zone center are absent. For KFe 2 As 2 it was suggested that there is a crossover from a nodeless s ± state in optimally doped Ba 1−x K x Fe 2 As 2 to a nodal d-wave state in KFe 2 As 2 upon hole doping [57][58][59][60] . Before reaching the d-wave phase, the system would be expected to pass through an s + id state at intermediate doping, in which time reversal symmetry is broken 61 . On the other hand, other results suggest that even in the heavily hole doped region, the superconducting state remains s-wave 62,63 . ARPES measurements report evidence that KFe 2 As 2 has s-wave symmetry with a nodeless gap on the innermost hole Fermi surface, an octet line node structure on the middle surface, and a very small gap on the outer one 63 . For the FeSe based materials with no hole pockets at the zone center, there have been proposals for a different s ± state compatible with the electronic structure 61,64 . Another proposed alternative is an orbital selective s × τ 3 state, where the interband and intraband pairing functions have different d-wave symmetries, resulting in a sign change of the intraband pairing function between the electron pockets, but the superconducting gap remains fully open 65 .
Unconventional superconductors often exhibit the breaking of additional symmetries in the superconducting state, as well as the broken global gauge symmetry common to all superconductors. One possibility is that time reversal symmetry (TRS) may be broken, which leads to the spontaneous appearance of magnetic fields below T c . Two representative examples showing such behavior are the unconventional superconductors Sr 2 RuO 4 and UPt 3 , where evidence for TRS breaking was initially found from zero-field µSR 66,67 , which were later corroborated by optical Kerr rotation experiments 68,69 (although in Ref. 70 conflicting results were found for µSR measurements of UPt 3 ). The broken TRS is well accounted for by the leading theories of triplet superconductivity for these systems, where there are two-component order parameters corresponding to two-dimensional irreducible representations of the crystallographic point group 71,72 . However, in other cases such as LaNiC 2 and LaNiGa 2 , the broken TRS is ascribed to non-unitary triplet pairing [73][74][75] , where the Cooper pairs have an intrinsic spin polarization 76 . Meanwhile, it has also been proposed that loop Josephson-currents can account for the TRS breaking signals in some otherwise apparently conventional superconductors 77 . In the case of iron pnictides, µSR measurements were performed on Ba 1−x K x Fe 2 As 2 to look for the evidence of TRS breaking s+id state predicted to occur between the s ± and d-wave phases 61 . However, these yielded conflicting results, with one study finding no signatures of broken TRS 78 , while evidence was found from measurements of ion irradiated samples 79 . As such, performing µSR measurements to look for the presence of broken TRS is an important method for identifying the pairing state of new unconventional superconducting materials.
Recent discoveries have led to a particular interest in superconductivity of Cr-based materials. For example, CrAs, exhibits helimagnetic order at T N = 273 K at ambient pressure and becomes superconducting with T c = 1.5 K in an applied pressure of 1.09 GPa 80,81 . This spurred the synthesis of new Cr-based superconductors, which resulted in the discovery of superconductivity in A 2 Cr 3 As 3 (A = Na, K, Rb, and Cs), which have a quasi-1D crystal structure [82][83][84][85] . Subsequently, superconductivity was also observed in isostructural A 2 Mo 3 As 3 (A = K, Rb and Cs) 86,87 . As the elements Cr and Mo belong to the same group of transition metals with a similar electronic configuration, this poses the question of whether these two families of superconductors share some common origins for the superconductivity 86,87 . The A 2 Cr 3 As 3 superconductors have been intensively investigated both experimentally and theoretically 82-85,88-100 , as they are strong candidates for multiband triplet pairing. Proposals of either triplet p z or f -wave pairing states arise from theoretical calculations 92,93 , which may also account for the evidence of line nodes in the superconducting gap 90,96,100 . Various techniques can be used to investigate the gap structure of superconductors including measuring the temperature and field depen- dence of thermodynamic properties, angle-resolved measurements, and spectroscopic techniques. To completely understand the gap structure and pairing state, it is often necessary to consider a variety of measurements, and here we primarily discuss µSR investigations of unconventional superconductors. The temperature and field dependence of TF-µSR measurements yields information about the flux-line lattice and magnetic penetration depth, which can reveal the structure of the superconducting gap. Furthermore, TF-µSR measurements on single crystals can allow for the gap anisotropy to be directly probed 101,102 . µSR in zero-field can be used to determine whether time reversal symmetry is broken in the superconducting state, as well as the magnetic properties of the system. Detailed reviews about µSR measurements on superconductors and other condensed matter systems can also be found in Refs [103][104][105][106]108 . In this review, we principally concentrate on using µSR to probe the superconducting pairing states of newly discovered double layered FeSCs, as well as a family of Crbased one-dimensional materials. In the first part, a brief introduction to the muon-spin rotation/relaxation techniques is given, while the second part reviews the superfluid density of ACa 2 Fe 4 As 4 F 2 (A = K, Cs, and Rb), ThFeAsN, doped BaFe 2 As 2 and A 2 Cr 3 As 3 (A = K, and Cs). Finally, we close with concluding remarks where we discuss the results and suggest future directions. We note that this review only includes a small subsection of the extensive research on high-temperature superconductors, and further detailed reviews can be found in Refs. 19,21,52,108,[111][112][113] .

II. INTRODUCTION TO THE µSR TECHNIQUE
Muon spin rotation and relaxation are very sensitive local probes which can be used to investigate magnetic and superconducting materials (for further details see Refs. [103][104][105][106]108 ). The muons utilized in µSR facilities are generated by colliding a high energy proton beam with a graphite target, which produces pions, and each positive (negative) pion decays into a positive (negative) muon and a (anti)neutrino. As a result of parity violation, the muon beam is spin-polarized, and these muons are implanted into the sample and come to rest via electrostatic interactions which preserve the muon spin direction. For µSR measurements, positive muons µ + are usually utilized, which decay with a half-life of 2.2 µs into a pair of neutrinos and a positron. Since the positrons are more likely to be emitted along the axis of the muon spin, the angular distribution of the emitted positrons reflects the distribution of the spins of the muons. In the presence of any local magnetic fields at the muon stopping site B loc , the µ + spins precess at the Larmor frequency ω µ = γ µ B loc , where γ µ /2π = 135.5 MHz T −1 is the muon gyromagnetic ratio. The time dependence of the polarization of the implanted muons is given by P µ (t) = G(t)P µ (0), where the function G(t) corresponds to the µ + spin autocorrelation function, which is determined by the internal magnetic field distribution 103 . A schematic illustration of a µSR experiment (adopted from Ref. 108 ) is shown in Fig. 4(a). A commonly measured quantity in µSR experiments is the time dependent asymmetry A(t) which is proportional to P µ (t) and is given by where N F (t) and N B (t) are the number of positrons detected in the forward and backward positions respectively, and α is a calibration constant. µSR measurements performed in zero external magnetic field are very sensitive to weak internal fields in a sample, down to the order of 0.1 G. In studies of superconductors this can be utilized to look for time reversal symmetry (TRS) breaking in the superconducting state. TRS breaking corresponds to the spontaneous appearance of magnetic fields below T c , which in turn leads to a more rapid muon spin depolarization. While a number of expressions can be utilized to analyze ZF-field µSR spectra, the data are often analyzed using 66,73 where λ is the Lorentzian relaxation rate, A 0 is the initial asymmetry of the relaxing component, and A bg is the background term. G KT (t) is the Kubo-Toyabe function given by where σ KT is the Gaussian relaxation rate. Here G KT characterizes the decay of the asymmetry due to a Gaussian distribution of magnetic fields which are static on the timescale of the muon lifetime (including those arising from nuclear moments), while the Lorentzian decay can correspond to either a Lorentzian distribution of static fields or rapidly fluctuating fields, and hence is often termed the electronic relaxation rate. The breaking of TRS in the superconducting state can be inferred from either an increase of λ 66,73 , or σ KT 74 below T c . The application of a longitudinal field (LF) parallel to the initial muon-spin serves to decouple the muon from the local magnetic fields in the sample. This decoupling occurs more rapidly for static fields than for dynamic fields, and therefore LF-µSR reveals information about both the magnitude and dynamics of the internal fields 104 . On the other hand, µSR measurements in applied transverse magnetics field (TF), perpendicular to the initial muon spin direction, are important for probing the flux line lattice 106 . In TF-µSR, the precession of the muon spins leads to oscillations in A(t), which decay due to the presence of a broadened field distribution P (B). The shape of the field distribution corresponding to the flux line lattice is shown in Fig. 4(b) 108 , and P (B) can be quantitatively utilized to extract the magnetic penetration depth λ and the coherence length ξ 106,116,117 . The decay of the oscillations in A(t) is often modelled using a Gaussian function and therefore A(t) can be analyzed using where ν s and ν bg are the frequencies of the muon precession from the sample and background, respectively with a phase offset φ. The Gaussian decay rate in the first term σ has two contributions below T c , where σ = (σ 2 sc + σ 2 nm ), and σ sc is the contribution from the vortex lattice and σ nm corresponds to the nuclear dipolar relaxation rate, which can be obtained from spectra measured in the normal state. σ sc can be directly related to the effective penetration depth λ ef f via σ sc /γ µ = 0.0609Φ 0 /λ 2 ef f , where Φ 0 is the magnetic flux quantum. This relation between σ sc and λ ef f is valid for 0.13/κ 2 (H/H c2 ) 1, for a Ginzburg-Landau parameter κ=λ/ξ 70 117 . Since λ(T ) is sensitive to low energy excitations, this in turn can provide information about the superconducting gap structure. For instance, if the superconducting gap is fully open across the whole Fermi surface, λ(T ) shows exponentially activated behaviour, but if there are line or point nodes in the gap, λ(T ) shows a linear or quadratic temperature dependence in the clean and dirty limits, respectively 118 . In particular, since σ sc is proportional to the superfluid density n s , it can be modelled using 96,119,120 , The temperature dependence of the superconducting gap function is approximated by the expression δ(T /T c ) = tanh{1.82[1.018(T c /T − 1)] 0.51 } while g(φ) describes the azimuthal angle dependence of the gap, given by 1 and |cos(2φ)| for an isotropic s−wave gap and a d−wave gap with line nodes, respectively 121,122 . A range of other potential gap structures have been used to analyze the results from iron-based superconductors, including g(φ) = 1 + αcos2φ for an extended s-wave state 123 , and ∆(T, φ)=[(∆ 1 (0)) 2 +(∆ 2 (0)sin(φ)) 2 ] 1/2 for an s×τ 3 state, where components corresponding to two different d-wave symmetries are summed in quadrature 65,199 . Therefore by analyzing σ sc with various gap models, the magnitude and structure of the superconducting gap can be probed, which in turn allows for the identification of the superconducting pairing symmetry.

III. MULTIBAND SUPERCONDUCTIVITY IN ACa2Fe4As4F2
A. Crystal structure and superconductivity Recently, Wang et al. discovered superconductivity at 29-33 K in ACa 2 Fe 4 As 4 F 2 (A = K, Rb, Cs, 12442-family) 125,126 . Figures 1 (e)-(f) show the tetragonal crystal structure with space group I4/mmm (No. 139, Z = 2) in which the Fe 2 As 2 layers are sandwiched between A atoms on one side and Ca 2 F 2 on the other side, leading to two distinct As sites above and below the Fe plane. This crystal structure is similar to the CaAFe 4 As 4 -family (A = K, Rb, Cs, 1144 family) and since the superconducting properties are closely associated with the nature of the spacer layers between adjacent conducting layers 125,126 , it is of great interest to study the superconductivity of compounds with these double Fe 2 As 2 −layers. Furthermore, it was noted that for the 12442-family, T c has an inverse correlation with the lattice parameters a and c 125,126 . This is in contrast to the case of optimally hole-doped AFe 2 As 2 systems, where the maximum T c increases with the lattice parameters 125,126 . More specifically, for the 12242 materials T c decreases with increasing d intra (the spacing of Fe planes within the double Fe 2 As 2 −layers), yet T c increases almost linearly with the distance between the Fe planes separated by the Ca 2 F 2 layers, d inter . In FeSCs, T c has often been found to be related to both the As-Fe-As bond angle, and the As height from the Fe plane 127 , and for the 12442-family there are two sets of these parameters due to the asymmetric Fe 2 As 2 layers leading to inequivalent As sites 125,126 . In addition, unlike many other FeSCs, no structural phase transitions or SDW transitions have been reported at low temperatures 125,126 . These materials are intrinsically near optimal doping with hole conduction (0.25 holes/Fe 2+ ) 128 . Electronic structure calculations for KCa 2 Fe 4 As 4 F 2 reveal a multiband character and suggest that the system is close to a stripe antiferromagnetic ground state 129 . It is suggested that the magnetic ground state is suppressed by the self-doping, leading to the observed high temperature superconductivity.
B. Superconducting gap structure and time reversal symmetry

KCa2Fe4As4F2
KCa 2 Fe 4 As 4 F 2 with T c = 33.36 K exhibits multigap superconductivity with line nodes as reported by Smidman et al. 130 . TF-µSR data reveal that the temperature dependence of the in-plane penetration depth λ −2 ab does not have the form of a s-wave superconductor, since it does not become flat and saturate at low temperatures as shown in Fig. 5(b). Furthermore, λ −2 ab cannot be described by either single band models or a multigap s-wave expression, but are well described by two-gap models with line nodes. Here good fits are obtained for both a s + d model where one gap is fully open and the other is nodal, as well as a d + d scenario where both gaps have line nodes. The respective fitted values of the larger of the two gaps are 3.52 k B T c and 5.08 k B T c , which are significantly higher than the weakly-coupled gap values (1.76 k B T c and 2.14 k B T c for s-and d-wave, respectively), suggesting strongly coupled superconductivity. An interesting comparison can be made with the ARPES studies of two other nodal iron-pnictide superconductors BaFe 2 (As 0.7 P 0.3 ) 2 and KFe 2 As 2 63,131 . For BaFe 2 (As 0.7 P 0.3 ) 2 it was found that the superconducting gap nodes only exist on a small region of the Fermi surface 131 , and it was suggested that this relatively limited nodal area can account for the relatively high T c of 30 K. On the other hand, the ARPES study of KFe 2 As 2 revealed a much more extended nodal area of the Fermi surface with an octet line node structure, and in turn this compound has a much lower T c of around 3 K 63 . In KCa 2 Fe 4 As 4 F 2 the superfluid density does not drop as dramatically with increasing temperature as would for instance be expected for a weakly coupled d-wave superconductor, which also suggests a relatively small nodal region of the Fermi surface as in BaFe 2 (As 0.7 P 0.3 ) 2 , which is consistent with the similarly large T c .
show the results from ZF-µSR, which indicate that TRS is preserved in the superconducting state. In addition, a temperature dependence of the relaxation rate well above T c suggests the presence of spin fluctuations in this system. It is of interest to understand the origin of nodal superconductivity in this system and in particular how changing the spacer layers in these asymmetric FeSC tunes the properties of the pairing state. In comparison, the 1144 superconductor CaKFe 4 As 4 exhibits two-gap nodeless superconductivity [134][135][136] , and therefore the exchange of a Ca layer for Ca 2 F 2 can tune the pairing state from nodeless to nodal in these materials 132 . Since it is difficult to account for this change along similar lines to other nodal FeSCs, this suggests that this material represent a different means of realizing nodal iron-based superconductors. Furthermore, from the TF−µSR study superconducting parameters of an in-plane penetration depth λ ab (0) = 229.5 nm, superconducting carrier density n s = 1.39 × 10 27 m −3 , and carrier's effective-mass m * = 2.59 m e are estimated using Θ D = 366 K, and λ e−ph = 1.59.

RbCa2Fe4As4F2
Another member of the 12442 family RbCa 2 Fe 4 As 4 F 2 with T c = 29.19 K was investigated by Adroja et al. 109 . For this compound, the magnetic susceptibility and specific heat reveal the presence of bulk superconductivity with a T c of around 29 K. Figure 6(a) displays the TF−µSR results which show that λ −2 ab cannot be accounted for by either nodal or nodeless single-gap functions. The λ −2 ab data are well described by both an s + s model with two nodeless gaps, as well as an s + d model where one of the gaps has line nodes. From the TF−µSR study, superconducting parameters of an in-plane penetration depth λ ab (0) = 231.5 nm, superconducting carrier density n s = 1.29 × 10 27 m −3 , and carrier's effectivemass m * = 2.44 m e are estimated using Θ D = 352 K, and λ e−ph = 1.45. Additionally, the ZF−µSR study does not reveal a clear sign of TRS breaking below T c , but the temperature dependence of the relaxation rate of the

CsCa2Fe4As4F2
Kirschner et al. 133 reported a TF-and ZF-µSR study on CsCa 2 Fe 4 As 4 F 2 with T c = 28.3 K to understand its superconducting and magnetic properties. They find that λ ab (T ) −2 of CsCa 2 Fe 4 As 4 F 2 does not plateau at low temperature, indicating nodal superconductivity as shown in Fig. 6(b). At intermediate temperatures, there is an inflection point in the temperature dependence of the superfluid density from which it is concluded that CsCa 2 Fe 4 As 4 F 2 is a nodal multigap superconductor, which can be fitted with a s + d-wave model. ZF−µSR measurements again show two components in the relaxation of the asymmetry, and both the Gaussian and a Lorentzian relaxation rates show a weak temperature dependence between 2 K and 50 K, without any clear sign of TRS breaking at T c , which is very similar to the results for KCa 2 Fe 4 As 4 F 2 and RbCa 2 Fe 4 As 4 F 2 109,130 . This compound contains the largest alkali atom in this family of superconductors, and λ ab (0) = 244 nm, n s = 1.16 × 10 27 m −3 , and m * = 2.44m e are estimated using Θ D = 344 K, and λ e−ph = 1.44. It is interesting to note that the multigap behavior observed in ACa 2 Fe 4 As 4 F 2 (A = K, Rb, and Cs) is very common in other FeSC materials (see Table-I). For example, recent µSR measurements on CaKFe 4 As 4 revealed nodeless two-gap behavior, with gap sizes of 2.5 and 8.6 meV 134 , which is consistent with the results from tunnel-diode resonator and scanning tunneling microscopy measurements 135 , as well as a high-resolution ARPES study 136 . Meanwhile NMR results reveal spin singlet superconductivity from the decrease of the Knight shift below T c , as well as nearly temperature independent behavior in the normal state 132 . Polarized Raman spectroscopy reveals no evidence for Pomeranchuk-like electronic nematic fluctuations 137 , suggesting that the occurence of electronic nematicity is not a requirement for realizing high−T c superconductivity in FeSCs. On the other hand, the analysis of the superfluid density of KCa 2 Fe 4 As 4 F 2 and CsCa 2 Fe 4 As 4 F 2 indicate that at least one gap has line nodes, which is very different to the fully gapped behavior of CaKFe 4 As 4 , while the analysis for RbCa 2 Fe 4 As 4 F 2 cannot distinguish between nodeless and nodal two-gap models. Why swapping the Ca layer in CaKFe 4 As 4 for Ca 2 F 2 induces nodal behavior needs to be determined by future studies since it does not correspond to a pressure effect, reduced height of the As atoms, or hole doping 130 .

IV. MULTIBAND SUPERCONDUCTIVITY IN ThFeAsN
The study of stoichiometric iron pnictide superconductors has been of particular interest, due to the advantage of having reduced disorder, which can be useful for the development of realistic theoretical models for understanding the unconventional superconductivity in the Febased materials. Recently, Wang et al. 138 140,141 , and as a result this material presents the scenario of a stoichiometric material near optimal T c , which has an electronic structure much like most other FeSCs.

A. Crystal structure and Tc
Wang et al. 138 reported that ThFeAsN crystallizes in the ZrCuSiAs-type tetragonal crystal structure (space-group P 4/nmm, No. 129, Z= 2) with lattice parameters a = 4.0367 and c=8.5262Å. In this structure there are separate layers of Th and N perpendicular to c-axis. The As and Fe layers are halfway along the c-axis, and the Fe and As ions form tetrahedrons with two As-Fe-As bond angles α ∼ 107.0 • and β ∼ 114.5 • at 300 K. The layered structure of ThFeAsN is very similar to others in the 1111 family of iron pnictide superconductors, such as LaFeAsO 19 . The ρ(T ) data shown in Fig. 7(a) shows a superconducting transition onsetting below T c = 30 K. As shown in Figs. 7 (b) and (e), bulk superconductivity below 30 K is revealed by measurements of the magnetic susceptibility and specific heat.  Fig. 7(g) is fitted better with models with two gaps, than either the single-gap isotropic s-wave or d-wave models 146 . It can be seen that similarly good fits are obtained with either an isotropic s + s wave model, or an s + d wave model which has one nodal gap as as well as a nodeless one. However as shown in Fig. 7(f), a nonlinear field-dependence of the specific heat can be observed in the low temperature limit. This is evidence for a nodal superconducting gap, suggesting that the (s + d)−wave model best explains the gap structure of ThFeAsN 146 . Furthermore, the estimated superconducting parameters were λ L (0) = 375 nm, n s = 4.97×10 27 m −3 , and m * = 2.48 m e 146 . Shiroka et al. 147 , also investigated ThFeAsN using TF-and ZF-µSR, as well as NMR measurements, which detect the presence of magnetic fluctuations, which weaken upon entering the superconducting state. This was taken as evidence for competition between superconductivity and magnetism, which is in contrast to the often observed scenario of microscopic coexistence 52,148 . They also report that the superconductivity is best accounted for by pairing states without gap nodes, with both a two-gap s-wave and anisotropic s-wave models accounting for the data. Moreover, their ZF-µSR study suggests that time reversal symmetry is preserved in the superconducting state, since there is only a very weak change of the relaxation rate with decreasing temperature.

V. MULTIBAND SUPERCONDUCTIVITY IN DOPED AFe2As2
The observation of high-temperature superconductivity in doped AFe 2 As 2 (A = Ba, Sr, Ca) with T c as high as 38 K in Ba 0.6 K 0.4 Fe 2 As 2 attracted significant interests in these materials 127,[157][158][159] . In particular, the research activities were sparked due to their simple ThCr 2 Si 2 -type crystal structure, the easy growth of large single crystals, and the similarity of their generic phase diagram with that of the high-T c cuprates. With a semimetallic ground state the parent AFe 2 As 2 compounds exhibit an itinerant antiferromagnetic spin density wave transition accompanied by a tetragonal to orthorhombic structural distortion and become superconducting upon suppression of the SDW transition by partial chemical substitutions at either the A, Fe or As sites or by application of external pressure. The emergence of superconductivity upon the suppression of long-range antiferromagnetic ordering is qualitatively similar to the observation of super-conductivity in high-T c cuprates where the parent compounds are antiferromagnetic insulators 127,[157][158][159][160][161] . Understanding the superconducting gap structure and the microscopic pairing mechanism is a key problem for FeSCs, and in general for the understanding of high-T c superconductivity. While two-gap superconductivity (with one small gap and one large gap) is generally believed to occur for the doped-122 superconductors, different superconducting gap structures have been proposed. For optimally doped Ba 1−x K x Fe 2 As 2 , Ba 1−x Rb x Fe 2 As 2 , Ba(Fe 1−x Co x ) 2 As 2 and Ba(Fe 1−x Ni x ) 2 As 2 a nodeless isotropic gap structure seems more appropriate [162][163][164][165][166][167][168][169] .
On the other hand signatures of nodal superconducting gap structures have been reported for BaFe 2 (As 1−x P x ) 2 and Ba(Fe 1−x Ru x ) 2 As 2 170-173 , as well as in overdoped Ba 1−x K x Fe 2 As 2 and Ba 1−x Rb x Fe 2 As 2 62,174 . The Fermi surface structure plays an important role in determining the superconducting gap structure. In optimally doped Ba 1−x K x Fe 2 As 2 , the Fermi surface consists of both electron and hole pockets, whereas the Fermi surface of the end point compound KFe 2 As 2 consists of only hole pockets 62,175,176 .

A. Crystal structure and superconductivity
The parent AFe 2 As 2 compounds form in the familiar ThCr 2 Si 2 -type body-centered tetragonal (bct) structure (space group I4/mmm, No. 139, Z = 2) which consists of alternating layers of A +2 and [Fe 2 As 2 ] −2 slabs stacked along the c-axis. In this structure each of the A, Fe and As atoms occupy only one crystallographic site. The A atoms form a bct sublattice and the Fe atoms form square-planar layers. Fe and As form FeAs 4 tetrahedra which are edge shared to form the three-dimensional networks of [Fe 2 As 2 ] −2 ; four of the six edges of FeAs 4 are shared. In doped-AFe 2 As 2 the T c is suggested to be related to both the As-Fe-As bond angle, and the distance (height) between the Fe layer and the adjacent pnictogen layer 127 . For the K-doped BaFe 2 As 2 , the highest T c appears for the As-Fe-As bond angle of ≈ 109.47 • within an undistorted FeAs 4 tetrahedron 127 .

B. Superconductivity in hole doped Ba1−xKxFe2As2 and Ba1−xRbxFe2As2
The hole doping induced by partial substitution of Ba by K in BaFe 2 As 2 , which exhibits structural and SDW transitions near 140 K 177,178 , suppresses the SDW transition leading to the development of superconductivity in Ba 1−x K x Fe 2 As 2 , with a maximum T c ≈ 38 K for x = 0.4, i.e. in Ba 0.6 K 0.4 Fe 2 As 2 179 . The µSR investigations on Ba 1−x K x Fe 2 As 2 have suggested the coexistence of magnetism and superconductivity as well as phase separation of superconducting (paramagnetic) and antiferromagnetic regions 162,164,[180][181][182] . The µSR study on Ba 0.6 K 0.4 Fe 2 As 2 by Hiraishi et al. 164 shows evidence for s-wave superconductivity where the σ sc (T ) data in the superconducting state could be fitted either by a singlegap structure with 2∆(0)/k B T c = 5.09(4) or by a two-gap structure with 2∆ 1 (0)/k B T c = 7.3 and 2∆ 2 (0)/k B T c = 4.1 (2). The large value of 2∆(0)/k B T c compared to the weak coupling BCS value of 3.53 clearly demonstrates strong coupling superconductivity in Ba 0.6 K 0.4 Fe 2 As 2 . A µSR study on Ba 0.5 K 0.5 Fe 2 As 2 by Goko et al. 162 found a nearly linear temperature dependence of σ in the superconducting state [see Fig. 8(a), figure from Ref. 162 ] which could not be explained by a single isotropic energy gap. Instead the result was suggested to represent either an anisotropic superconducting gap with line nodes or a multi-gap structure. The ARPES study on Ba 0.6 K 0.4 Fe 2 As 2 favours two band superconductivity with strong coupling 163,165,166 . Based on ARPES data, the presence of two isotropic and nodeless superconducting gaps and an inter-band scattering scenario for the pairing mechanism have been suggested for superconductivity in Ba 0.6 K 0.4 Fe 2 As 2 .
Attempts were also made to probe the time reversal symmetry breaking in Ba 1−x K x Fe 2 As 2 , however the µSR investigations on polycrystalline Ba 1−x K x Fe 2 As 2 (0.5 ≤ x ≤ 0.9) could not detect any evidence for spontaneous magnetic fields associated with TRS breaking 78 . A recent µSR study on high quality ion-irradiated Ba 0.27 K 0.73 Fe 2 As 2 single crystals reported evidence for TRS breaking 79 . It has been suggested that the broken TRS state exists only for a very narrow region of x from 0.7 to 0.8 in the phase diagram of Ba 1−x K x Fe 2 As 2 79 . The end point compound KFe 2 As 2 was proposed to exhibit a nodal d-wave pairing 59 , and a crossover from a nodeless gap structure to a nodal state upon heavy hole doping is expected to occur near x ≈ 0.8 62 . The change from nodeless to nodal gap structure is believed to be related to a Lifshitz transition i.e., the disappearance of the electron pockets in the Fermi surface, since while the Fermi surface of optimally doped Ba 0.6 K 0.4 Fe 2 As 2 possesses both electron and hole pockets, the Fermi surface of KFe 2 As 2 has only hole pockets 62,175,176 . The partial substitution of Ba by Rb in BaFe 2 As 2 also leads to superconductivity in Ba 1−x Rb x Fe 2 As 2 . It is observed that with Rb in place of K, the Fe-sublattice AFM can be completely suppressed, consequently, there is no phase separation of superconducting and antiferromagnetic regions in optimally doped Ba 1−x Rb x Fe 2 As 2 . The µSR study on Ba 1−x Rb x Fe 2 As 2 also finds evidence for two-gap superconductivity 167 . The temperature dependence of the penetration depth λ obtained from the µSR measurements on Ba 0.65 Rb 0.35 Fe 2 As 2 is shown in Fig. 8 The partial substitution of Fe by Co in BaFe 2 As 2 also suppresses the SDW transition and leads to the emergence of superconductivity in Ba(Fe 1−x Co x ) 2 As 2 as a result of effective electron doping, with a maximum T c ≈ 23 K for x ≈ 0.07 [184][185][186] . The µSR investigations on Ba(Fe 1−x Co x ) 2 As 2 report the coexistence of superconductivity and magnetism/SDW fluctuations even in the optimally doped case 168,169,[187][188][189][190][191] . A µSR study on Ba(Fe 0.926 Co 0.074 ) 2 As 2 showed evidence for two-gap superconductivity with 2∆ 1 (0)/k B T c = 3.77 and 2∆ 2 (0)/k B T c = 1.57 168,169 . The temperature dependence of the penetration depth λ is shown in Fig. 8(d) (figure from Ref. 168 ) as 1/λ 2 vs. T along with the fit with two-gap model. The parameter 2∆(0)/k B T c suggests that the superconductivity in Ba(Fe 0.926 Co 0.074 ) 2 As 2 is not in a strong coupling limit as was the case for Ba 0.6 K 0.4 Fe 2 As 2 , though they both show a two-gap superconductivity. On the other hand µSR investigations on SrFe 1.75 Co 0.25 As 2 reveal the two gaps to be 2∆ 1 (0)/k B T c = 7. Following the discovery of superconductivity under pressure in CrAs 80,81 , a new family of isostructural superconductors A 2 Cr 3 As 3 (A = Na, K, Rb, and Cs) were discovered [82][83][84][85] . These were found to exhibit superconductivity at ambient pressure, with higher values of T c of 8.6, 6.1, 4.8 and 2.2 K for A= Na, K, Rb, and Cs respectively. K 2 Cr 3 As 3 crystallizes in a hexagonal crystal structure [space group P6m2 (No. 187)] with lattice parameters a = 9.9832Å, c = 4.2304Å. 82 . Here the atoms do not form a planar hexagonal arrangement, but instead lie within two planes offset along the c-axis by half a unit cell, as illustrated in Figs. 9(a). The Cr and As atoms form quasi-one dimensional double-walled chains of [(Cr 3 As 3 ) 2− ] ∞ , consisting of face sharing octahedra with the outer wall being formed from As and with Cr atoms on the inside, as shown in Figs. 9(b)-(c). Within the Cr octahedra, the bond lengths lie between 2.61 and 2.69Å, suggesting that metallic bonding between Cr atoms is expected. The double-walled chains run along the crystallographic c-axis, and are separated by the K + cations, and hence giving rise to a quasi-onedimensional crystal structure.

B. Physical properties: K2Cr3As3 and Cs2Cr3As3
Specific heat measurements of K 2 Cr 3 As 3 reveal the presence of bulk superconductivity, where the large jump in the specific heat at the transition of 2.2 γT c is considerably greater than that expected from weak coupling theory, indicating that the superconductivity is strongly coupled 82 . From probing the normal state, strong electronic correlations were revealed by a large electronic specific heat coefficient and non-Fermi liquid transport behavior 82,88 , which are a common feature of unconventional superconductors. A lack of BCS pairing was also suggested from the very large and anisotropic upper critical fields (H c2 ), where H c2 perpendicular to the Cr−chains is significantly larger than that corresponding to the parallel direction. Moreover, paramagnetic pair breaking appears to be absent for the perpendicular direction, but parallel to the Cr−chain H c2 appears still to have paramagnetic limiting, but with an enhanced Pauli field 88,89 .
Measurements of the magnetic susceptibility show that superconductivity occurs below 2.1 K and the isothermal M (H) curve suggests type-II superconductivity in Cs 2 Cr 3 As 3 84,100 . The reported value of the upper critical field H c2 = 6.45 T is higher than the Pauli limit, H P = 1.84T c = 3.86 T, much like the other isostructural materials, which again indicates that the superconductivity is unconventional. In zero field, ρ(T ) shows metallic behavior with a linear temperature dependence from 50 K down to just above T c , indicating non-Fermiliquid behavior, suggesting the importance of spin fluctuations 84,100 .
C. Superconducting gap structure and time reversal symmetry

K2Cr3As3
Adroja et al. 96 reported a TF-and ZF-µSR study on K 2 Cr 3 As 3 . The temperature dependence of σ sc obtained from the TF-µSR measurements could be fitted to both an isotropic s-wave model, as well as a d-wave model with line nodes 96 . However, the goodness-of-fit values indicate that the data fit better to the d-wave model than the swave model [see Fig. 9(e)]. Therefore the µSR results are more consistent with having line nodes than a fully open superconducting gap, which is consistent with the findings of penetration depth measurements using a tunnel diode oscillator technique, where a linear temperature dependence is revealed at low temperatures in the highest quality samples 90 . The ZF-µSR study reveals some evidence for the appearance of weak internal fields below T c , suggesting that the superconducting state may be unconventional, but these results cannot ascertain whether there is spontaneous broken time reversal symmetry. Moreover, from the TF-µSR study the estimated superconducting parameters are λ L (0) = 432(4) nm, n s = 2.7 × 10 27 m −3 , and m * = 1.75 m e , respectively. The NMR study by Zhi et al. 95 on K 2 Cr 3 As 3 also supported the conclusion of unconventional non s-wave superconductivity from the absence of a Hebel-Slichter coherence peak in the spin-lattice relaxation rate 1/T 1 just below T c . The measurements of 1/T 1 also revealed enhanced Cr spin fluctuations in the normal state above T c , where the power law dependence of 1/T T 1 ∼ T −0.25 was suggested to be compatible with a Tomonaga-Luttinger liquid. Signatures of Tomonoga-Luttinger liquid physics were also found from ARPES measurements by Watson et al.,97 , and the relationship between the low dimensional spin fluctuations and the unconventional superconductivity requires further exploration. This results indicates that λ shows an increase below T c , suggesting that the superconductivity of Cs 2 Cr 3 As 3 is associated with the spin fluctuation mechanism, and is hence unconventional.

Cs2Cr3As3
Zhi et al. 98 studied the NMR and 75 As nuclear quadrupole resonance to probe the normal state of Cs 2 Cr 3 As 3 . Their Knight shift measurements corresponding to the Cs1 site show an increase of the uniform spin susceptibility upon decreasing the temperature from 295 K to around 60 K. Upon further lowering the temperature there is a slight decrease before flattening out below about 10 K. On the other hand, the Knight shift for the Cs2 is very small, indicating little contribution to the exchange interactions. From a comparison of the low temperature spin-lattice relaxation rates with K 2 Cr 3 As 3 and Rb 2 Cr 3 As 3 , it is suggested that substituting increasingly large alkali atoms both suppresses the Cr spin fluctuations and the superconducting T c 95,98,99 . This suggests a relationship between spin fluctuations and the superconductivity in ACr 3 As 3 , again pointing towards an unconventional pairing mechanism.
VII. UEMURA CLASSIFICATION SCHEME: In this section we focus upon the Uemura classification scheme 193,194 which is based on the correlation between the superconducting T c and the effective Fermi temperature, T F , determined from µSR measurements of the superconducting penetration depth. Within this scheme strongly correlated "exotic" superconductors, i.e. high T c cuprates, heavy fermions and organic materials, all lie within a particular part of the diagram, which is indicative of a degree of universal scaling 194 of T c with T F such that 1/10 ≥ (T c /T F ) ≥ 1/100. For conventional BCS superconductors 1/1000 ≥ T c /T F . Considering that the value of T c /T F = 30/4969.4 = 0.006 for ThFeAsN [see Fig.11], this material can be classified as lying close to the proposed limit 193,194 . Taking into account the value of T c /T F = 33.36/1285 = 0.026 for KCa 2 Fe 4 As 4 F 2 , T c /T F = 29.19/1260 = 0.023 for RbCa 2 Fe 4 As 4 F 2 , T c /T F = 28.31/1130 = 0.025 for CsCa 2 Fe 4 As 4 F 2 (see Fig. 11), these materials can be designated as exotic superconductors. For K 2 Cr 3 As 3 , T c /T F = 6.0/4689 = 0.0013 and Cs 2 Cr 3 As 3 = T c /T F = 2.1/1651 = 0.0012. Moreover, we have also included the data of ACa 2 Fe 4 As 4 F 2 (A = K, Rb and Cs) and A 2 Cr 3 As 3 (A = K and Cs) to Uemura plot as shown in Fig. 11, which shows that the former compounds belong to the exotic class, while the latter are away from this region.

VIII. CONCLUDING REMARKS
Despite extensive research activities worldwide following the discovery of FeSCs, the pairing states and mechanisms of superconductivity in these materials are still not well understood. Nevertheless, µSR is a powerful technique for elucidating these mysteries, in particular, due to the ability to characterize the superconducting gap structure, which is vital for revealing the underlying pairing interactions. Although conventional phonon-mediated superconductors exhibit an isotropic s-wave gap with a 2∆/k B T c ≈ 3.53 for weak-coupling, so far no consensus has been reached on the gap character of FeSCs. In this review, we discuss TF and ZF-µSR measurements in the normal and superconducting states of ACa 2 Fe 4 As 4 F 2 (A = K, Rb and Cs), ThFeAsN, A 2 Cr 3 As 3 (A = K and Cs) and in doped AFe 2 As 2 . The superfluid densities of ACa 2 Fe 4 As 4 F 2 (A = K, Rb and Cs) can be modelled using two-gap models, which better accounts for the data than single gap s-wave or dwave models 109,130,133 . Here an s + d model can describe the data for all three members of this family of materials, while A = K can also be accounted for by a d + d model and A = Rb can also be fitted with a s + s model. A larger value of 2∆/k B T c (∼ 6.5) than the 3.53 expected from BCS theory have also been obtained from the µSR investigations of all these compounds, indicating strongly coupled superconductivity. The observation of two gaps is very similar to what is often found in FeSCs (see Table-

I).
The nodal superconductivity observed in KCa 2 Fe 4 As 4 F 2 and CsCa 2 Fe 4 As 4 F 2 130,133 is different to the case of the CaKFe 4 As 4 , where in the latter distinct signatures are found for a nodeless s± pairing state [134][135][136] , despite both systems being near optimal hole doping 136 . Moreover the electronic structure calculations by Wang et al. 129 are in line with many other Fe-based materials with hole pockets at the zone center and electron pockets at the edge, and therefore a nodeless s± state might be expected. Consequently the change in superconducting gap structure between CaKFe 4 As 4 and KCa 2 Fe 4 As 4 F 2 is not similar to the case of (Ba 1−x K x )Fe 2 As 2 , where the change from nodeless to nodal superconductivity is brought about by hole doping 63,163,196,197 . Guguchia et al. 183 proposed that pressure leads to a change from nodeless to nodal superconductivity in Ba 0.65 Rb 0.35 Fe 2 As 2 , while it was also suggested that nodal superconductivity occurs in Fe based materials when the pnictogen is less than 1.33Å above the Fe-layer 195 . Nevertheless, in KCa 2 Fe 4 As 4 F 2 the As atoms are at heights of 1.40 and 1.44Å 130 , well above this upper limit and indeed above the values found for CaKFe 4 As 4 (1.40 and 1.35 A) 198 . This then suggests that this represents a different scenario for the occurence of nodal superconductivity in FeSCs, and whether the asymmetry above and below the Fe-layers which is present in both CaAFe 4 As 4 and ACa 2 Fe 4 As 4 F 2 plays a role in this, is yet to be determined. Meanwhile for RbCa 2 Fe 4 As 4 F 2 , Adroja et al. 109 find that the nodeless s + s-and nodal s + d-wave mod-els can both fit the superfluid density well 109 . Since this study could not distinguish between the nodeless and nodal models, additonal measurements are important to determine whether the gap structure is similar to the other materials in this family. An outstanding puzzle in ThFeAsN is the origin of the superconductivity in a stoichiometric compound in the apparent absence of a nearby ordered phase 147 . Interestingly, band structure calculations predict a stripe antiferromagnetic ground state for ThFeAsN which is not observed 140,141 , but instead ZF−µSR reveals the presence of magnetic fluctuations 146,147 . Adroja et al. 146 determine that for ThFeAsN, σ sc (T ) can be modelled with a two-gap model with either an isotropic s+s or a s+d wave gap structure, instead of a single-gap isotropic s-wave, anisotropic s-wave, or d-wave models, which is in accordance with the common trend of multiband superconductivity in FeSCs. However, it was suggested that the field-dependent heat capacity can discriminate between the two-gap nodal and nodeless scenarios, and from this it was concluded that the (s + d)-wave model is a better description of the gap function of ThFeAsN 146 . Additionally, the gap value of 2∆/k B T c = 4.29 obtained from the superfluid density analysis suggests strongly coupled superconductivity in ThFeAsN 146 . On the other hand, Shiroka et al. 147 concluded that the TF-µSR results are best described by nodeless two-gap superconductivity. These results may be important for the proposal of theoretical models to explain the origin of unconventional superconductivity in ThFeAsN. It would be of particular interest to examine the nature of the superconductivity in the second superconducting dome in ThFeAsN 1−x O x (x ∼ 0.3) using µSR and neutron scattering measurements.
The µSR investigations of both hole-doped and electrondoped BaFe 2 As 2 have been very useful for understanding the nature of the superconducting order parameter and two-gap superconductivity is well supported in these doped 122 systems. µSR studies on hole doped Ba 1−x K x Fe 2 As 2 have shown two isotropic and nodeless superconducting gaps in optimally doped Ba 0.6 K 0.4 Fe 2 As 2 with a strong coupling 164 . Evidence for TRS breaking is also found in ion-irradiated Ba 0.27 K 0.73 Fe 2 As 2 79 . Furthermore, a crossover from a nodeless to nodal gap is proposed to be induced by heavy hole doping near x ≈ 0.8 62 , associated with the disappearance of the electron pockets in the Fermi surface at a Lifshitz transition in heavily hole doped materials. This kind of crossover was also seen in Ba 1−x Rb x Fe 2 As 2 , which is induced by both heavy hole doping as well as by hydrostatic pressure 167,174 . The µSR investigations on electron-doped Ba(Fe 1−x Co x ) 2 As 2 have also shown evidence for two-gap superconductivity, however, the superconductivity does not seem to be strongly coupled 168,169 .
An important puzzle for the A 2 Cr 3 As 3 (A = K, Rb, and Cs) family of superconductors is the nature of the superconducting pairing symmetry, especially given the clear evidence for line nodes in the superconducting gap from µSR and tunnel diode oscillator based measurements 90,96,100 . One of the leading pairing symmetries predicted from theoretical proposals for a broad range of parameters is a triplet p z -wave pairing state arising due to ferromagnetic fluctuations 92,93 . This state has a nodal plane at k z = 0, and is also able to account for the disappearance of paramagnetic limiting in the upper critical field perpendicular to the chain direction 88,89 . On the other hand, the anisotropy of the paramagnetic limiting fields in A 2 Cr 3 As 3 may also be a consequence of the strong antisymmetric spin-orbit coupling arising due to the noncentrosymmetric crystal structure 199 . Another possibility is an f -wave pairing state which has three line nodes within the ab-plane 92 , which is compatible with the observed in-plane six fold modulation of the upper critical field 200 . Meanwhile the presence of spin fluctuations apparently modified below T c is revealed from ZF-µSR, pointing towards a link between these fluctuations and the superconductivity 96,100 , in line with the above theoretical scenarios. With the findings of low dimensional fluctuations compatible with a Tomonaga-Luttinger liquid 95,97 , it is of particular importance to understand any relation between these and the unconventional superconductivity. Despite extensive efforts, there are still many unresolved questions related to the gap structure and pairing states of FeSCs. A thorough understanding is all the more important given the potential for novel applications, particularly those related to quantum computation 52 . For instance, it is important to determine the origin of nodal superconductivity in some systems, and whether these are manifestations of the same multigap sign-changing swave state, or if multiple pairing symmetries are realized. In addition to the various gap structure probes, inelastic neutron scattering investigations have been very important for elucidating the pairing states of FeSCs 52 , and therefore INS studies on all the Fe-based and Cr-based materials presented in this review are highly desirable. More generally, comprehensive studies using a broad variety of complementary methods are very important for a clear elucidation of the gap structures and pairing symmetries of FeSCs and related materials, as well as for the study of proximate magnetic and nematic phases and their relationship to superconductivity. To this end µSR is an important tool for such studies, in particular being a sensitive local probe which can be used to determine the superconducting gap structure. Together with other methods such as thermodynamic probes, ARPES, inelastic neutron scattering, resonant inelastic X-ray scattering, and nuclear magnetic resonance, this can provide vital information for revealing the microscopic origin of high-T c superconductivity.