Importance of dxy\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{d}}_{{{xy}}}$$\end{document} orbital and electron correlation in iron-based superconductors revealed by phase diagram for 1111-system

The structural flexibility at three substitution sites in LaFeAsO enabled investigation of the relation between superconductivity and structural parameters over a wide range of crystal compositions. Substitutions of Nd for La, Sb or P for As, and F or H for O were performed. All these substitutions modify the local structural parameters, while the F/H-substitution also changes band filling. It was found that the superconducting transition temperature Tc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{\text{c}}$$\end{document} is strongly affected by the pnictogen height hPn\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h_{Pn}$$\end{document} from the Fe-plane that controls the electron correlation strength and the size of the dxy\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{xy}$$\end{document} hole Fermi surface (FS). With increasing hPn\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h_{Pn}$$\end{document}, weak coupling BCS superconductivity switches to the strong coupling non-BCS one where electron correlations and the dxy\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{xy}$$\end{document} hole FS may be important.

In this study, we extended our previous work on LaFe(As,P)(O,F/H) 29 by covering a wider range of structural parameters through three site substitutions (Nd for La, Sb or P for As, and F or H for O). All these substitutions modify the local structural parameters, while the F/H-substitution also changes band filling. (In other iron-based superconductors, AFe 2 As 2 systems (A = Ba, Sr, and Ca), the phase diagrams and T c have been investigated from the viewpoints of local structural parameters and carrier doping level 30 .) Precise measurements of structural parameters and resistivity for all samples revealed that the d xy hole FS is crucially important in a wide doping range for the high h Pn compounds. Although there are apparently three SC phases, two of them (SC1 and SC3) turn out to be of the same origin. It is likely that with increasing h Pn electron correlation becomes strong, which switches the FS nesting based weak coupling superconductivity (SC2) to other type of superconductivity such as the one derived from orbital-selective Mott systems (SC1/SC3). T-linear resistivity was commonly observed near the T c -maximum composition in each phase.

Results
First, we present structural data for several typical compositions. (Data for all remaining compositions are provided in the "Supplementary Information".) Figure 1a,b,f,g show the lattice constants a and c as functions of P/ Sb content (x) for RFeAs 1−x (P/Sb) x O 1−y (F/H) y with y = 0.1 and 0.3 (R = La and Nd). Part of the dataset is taken from our previous papers [27][28][29] . The P/Sb contents (x) were determined by energy-dispersive X-ray spectroscopy (EDX), while the F/H compositions (y) are nominal values. Here we see that the lattice constants vary linearly with x according to Vegard's law, demonstrating the successful substitution of P or Sb for As. The result indicates that P-substitution has the effect of lattice compression, while the Sb-substitution has the opposite effect. More detailed structural parameters, such as the pnictogen (Pn) height from the Fe-plane ( h Pn ), the Pn-Fe-Pn bond angle ( α ) and the Fe-Pn bond distance ( d Fe−Pn ) are presented in Fig. 1c-e,h-j. These parameters also vary linearly with x. With increasing Sb-content, thereby expanding the lattice, h Pn and d Fe−Pn increase, while α decreases. In Nd-systems, although the lattice constants are smaller, h Pn is larger and α is smaller than in the La-systems. As will be discussed later, among all these parameters, h Pn and/or α plays the most crucial role in determining the electronic state. From the present results, we find that the common primary effect of P/As(Sb), O/F(H) and La/ Nd substitutions is to increase h Pn .
To observe the change in carrier interaction, we measured the temperature (T) dependence of resistivity ( ρ ). In Fig. 2, the T-dependence of ρ − ρ 0 is plotted on a log-scale, where ρ 0 is residual resistivity. In all systems, with increasing P-content x, ρ − ρ 0 changed from non-Fermi liquid-like T-linear behavior to T 2 -behavior. This implies that As-rich compositions are in the strong interaction regime, while with increasing P-content, the interaction is weakened, and the system becomes a Fermi liquid. (It should be noted that, although the T 2 -behavior of resistivity could be induced by disorders in the strong coupling superconductors 31 , residual resistivity in the present compounds systematically decreases with x and reaches the lowest value at x = 1.0 where ρ ~ T 2 . This means that T 2 -behavior of resistivity is not due to disorder in the present case.) Accordingly, as seen in the insets, T c decreases with increasing P-content. www.nature.com/scientificreports/ In Fig. 3, we plot T c and the power n in ρ(T) = ρ 0 + AT n for all the P-doped samples in Fig. 2. Figure 3a,b show the results for y = 0.1 . In both La-and Nd-systems, at P100% composition ( x = 1.0 ), T c is very low (about 5 K) and n ∼ 2 . With decreasing x, T c gradually increases, while the power n monotonically decreases towards unity. This clear correlation between T c and n demonstrates that the key interaction for superconductivity is strengthened with decreasing x, as we have previously reported [27][28][29] . When x decreases below x = 0.4 , T c increases further in the Nd-system, while it is almost unchanged in the La-system. Our previous studies demonstrated    www.nature.com/scientificreports/ that there is a critical change in the electronic state at around x = 0.4 . This was evidenced not only by transport properties such as resistivity and Hall coefficient, but also by direct observation of the band structural change through the angle-resolved photoemission experiment 13 . We distinguish these two superconducting regions, naming the As-rich region "SC1" and the P-rich region "SC2". According to the theoretical calculations by Kuroki et al. 26 , when h Pn increases, the d xy hole FS around (π, π) becomes larger. Therefore, when h Pn decreases with x, as indicated in Fig. 1, the d xy hole FS is expected to shrink and eventually disappear as illustrated in Fig. 3e. We believe this Fermi surface change drives the phase change from SC2 to SC1. The latter can be defined as the phase in which the d xy hole FS contributes to superconductivity, while this is not expected in SC2. (Recently we observed d xy band in NdFeAs 1−x P x O 0.9 F 0.1 by angle-resolved photoemission spectroscopy. It was experimentally confirmed that the d xy band approaches the Fermi level at x < 0.2, forming an incipient band 32 .) For heavy electron-doping ( y = 0.3 ), the n-value is much larger than one at x = 1.0 , while it rapidly decreases when approaching x = 0 . (See Fig. 3c,d.) At x = 0 where the high T c of over 40 K is measured, n is close to 1, namely T-linear resistivity is observed. It is a common feature in the iron-based superconductors that while the P-rich compound is a weakly-interacting (Fermi-liquid like) system, the interaction becomes stronger with increasing As-content, which induces superconductivity and increases T c . A significant difference from the lightly-doped case ( y = 0.1 ) is that superconductivity disappears with P-substitution. As depicted in Fig. 3f, we expect that, with heavy electron doping, the electron FS becomes much larger than the hole FS at the zone center. This imbalance of hole-and electron-FS must prevent FS nesting-derived superconductivity in the P-rich compositions. Therefore, superconductivity in the P-rich compounds can be understood within the FS nesting scenario over a wide doping (y) range.
Since the d xy hole FS at the zone corners shrinks and eventually disappears with increasing P-doping level x 33 , we conclude that the d xy orbital is crucially important for superconductivity for heavily electron doping. In the Nd-system, T c remains high up to a larger x-value than in the La-system. This can be attributed to the larger d xy hole FS due to the higher h Pn in the Nd-system compared to the La-system.
We distinguish this heavily electron-doped region from the lightly doped region (SC1) by labelling it SC3, because the T c is strongly suppressed in the intermediate range of y in the La-system. However, in the Nd-system, there is no such T c suppression, and high T c is maintained over a wide y-range. Considering that the difference between the La-and Nd-systems is caused by the greater pnictogen height ( h Pn ), we enlarged h Pn by Sb substitution for As in the La-system. An interesting example at y = 0.14 is shown in Fig. 4. In the Nd-system (Fig. 4b), the x-dependence of T c at this doping level is qualitatively the same as that for y = 0.10 . However, in the La-system (Fig. 4a), T c is strongly depressed near x = 0 , while it recovers with increasing Sb-content, probably because of the enlargement of the d xy FS via the increase in h Pn . This suggests that SC1 and SC3 have the same origin. Figure 5a,b are contour plots illustrating how T c changes with x and y for the La-and Nd-systems, respectively. In these phase diagrams, we distinguish three superconducting phases (SC1, SC2 and SC3), particularly in the La-system. Compared to the change from SC1 to SC2 with the disappearance of the d xy FS, the change from SC1 to SC3 is not clearly defined. In the La-system, although these two phases are separated at x = 0 , when Sb is substituted, they merge with each other, as in the Nd-system, owing to the enlargement of the d xy FS by the increase in h Pn . We therefore think that the superconductivity mechanisms in phases SC1 and SC3 have the same origin. According to this scenario, the T c suppression near x = 0 for the intermediate doping range ( 0.14 < y < 0.25 ) is not caused by the bad FS nesting condition, but by the loss of d xy FS due to electron doping. The lost d xy FS and suppressed T c can be recovered by Sb substitution and/or electron doping. www.nature.com/scientificreports/ To visualize the relation between the local structural parameters and T c , Fig. 6 presents a contour plot of T c in the α-d Fe−Pn plane. With electron doping, the bond angle ( α ) decreases. Moreover, T c increases in the order of SC2, SC1 and SC3 as the structural parameters approach those for the regular tetrahedral structure.

Discussion
This figure indicates that the local crystal structure strongly affects the T c value in Fe-based superconductors. In SC3 phase where h Pn is close to 1.38 Å, T c is very high and insensitive to the band filling. Note that band filling should be another important parameter in weak coupling superconductivity where bosonic fluctuation via FS nesting acts as a pairing glue. In fact, in Fig. 6, the T c suppression due to electron doping is visualized by the white bands extending from upper left to lower right. Therefore, in SC2 phase (P-rich compositions), the FS-nesting condition is crucially important, which supports a weak coupling regime [34][35][36] . However, in the large h Pn region, there is no white area, namely, the FS-nesting condition does not affect T c at all.
One of the theoretical scenarios is derived from strong electron correlations near the d 5 Mott insulator regime [16][17][18] . In this scenario, the total electron correlation is strong enough to induce antiferromagnetic order in LaFeAsO. Then, superconductivity is realized when the carriers are doped in this strongly-correlated parent compound. The electron correlation is enhanced with increasing h Pn . Because h Pn can be controlled by As(Sb)/P, and Nd/La substitutions, the T c enhancement with these substitutions can be roughly interpreted as a result of increase in electron correlation. O/F(H) substitution has two effects, namely, carrier doping and increase in h Pn . The former weakens the electron correlation, while the latter strengthens it. The apparent separation of SC1 and SC3 in the La-system (Fig. 5a) is a result of the competition of these two effects. In Fig. 6, we can see that the electron correlation is enhanced towards the optimal point. T-linear resistivity observed in SC3 (Fig. 3) is a common property of strongly-correlated superconductors. Another strength of this model is that the electron  www.nature.com/scientificreports/ correlation depends on the orbital, and the d xy orbital plays a crucial role. This is consistent with the present results for SC1 and SC3.
There is another theoretical model for iron-based superconductors, shape resonance model [19][20][21][22][23] , where T c shows a clear enhancement around Lifshitz transition. In NdFeAs 1−x P x (O,F), we observed the change of FS topology at Brillouin zone corner around x = 0-0.2, and that at zone center around x = 0.2-0.4 by angle-resolved photoemission spectroscopy measurement 13 , which may correspond to Lifshitz transition. However, the T c value changes gradually with x but does not show any abrupt change at a particular composition. So far, it is unclear whether this model can explain the observed complicated phase diagram in Fig. 5. Further experimental and theoretical studies are required to clarify that this model can be applied to the superconductivity mechanism of iron-based superconductors.
The evidence for the deviation from a weak coupling BCS regime can be seen in the so-called Uemura plot 37 . This famous plot demonstrates that many exotic superconductors, including cuprates, show a linear correlation between T c and muon relaxation rate σ that is proportional to the superconducting carrier density. The important fact here is that these superconductors are characterized by low carrier density and high T c , which suggests that the coupling interaction is extremely strong in these compounds, even beyond that of strong coupling superconductivity described by Eliashberg's theory. Namely, they are beyond BCS superconductors, and may be close to BCS-BEC crossover. Shortly after the discovery of Fe-based superconductors, σ was determined for several RFeAs(O,F) (R = La, Nd and Ce) compounds 38 . All the data points for the measured compounds lie on the Uemura line, indicating that these compounds are categorized as strong coupling exotic superconductors beyond BCS theory.
On the other hand, σ for the low T c ( ∼ 5 K) compound LaFePO was found to be comparable to that for other high T c ( ∼ 30 K) compounds 38 . This means that the LaFePO data point in the T c -σ plot is far from the Uemura line. Recently, we performed systematic measurements of σ-values for various compounds of LaFeAs 1−x P x (O,F), and found that with increasing As-content the data point approaches the Uemura-line 39 . This result implies that the coupling strength perhaps originating from electron correlation increases with increasing As-content, as we expect from this study, and as Misawa et al. predicted in their calculations 16,18 . Therefore, in the SC2 phase, the system changes from a weak coupling BCS superconductor to a strong coupling non-BCS one with increasing As-content, while SC1 and SC3 are in the latter regime. Note that the BCS-BEC crossover in multiband superconductors must exhibit properties different from the single band case, as has been theoretically studied for the application to iron-based superconductors 40,41 .

Conclusion
The electronic phase diagram was extensively investigated for RFeAs 1−x (P/Sb) x O 1−y (F/H) y (R = La and Nd) over a wide composition range. It was observed that the pnictogen height ( h Pn ) increases with increasing As/Sbcomposition and electron doping, and that h Pn in the Nd-system is larger than that in the La-system. The changes in the superconducting transition temperature ( T c ) with composition can be understood in terms of the contribution of the d xy hole Fermi surface (FS) and electron correlation strength, both of which are controlled by h Pn . The apparently-separated two superconducting phases (SC1 and SC3) in the La-system merged with Sb-substitution, suggesting that the origins of these two phases are the same. While FS nesting controlled by band filling plays an important role in the weak coupling BCS superconductivity for the P-rich compositions (SC2 phase), as h Pn increases with increasing As/Sb-content, superconductivity regime switches to strong coupling exotic one with d xy orbital playing an important role (SC1/SC3), which may support the orbital-selective Mott scenario.

Methods
Polycrystalline RFeAs 1−x (P/Sb) x O 1−y F y (R = La and Nd) were synthesized using solid state reaction methods. For R = La, a stoichiometric-ratio mixture of LaAs, LaSb (or LaP), As (or P), Fe 2 O 3 , Fe and LaF 3 powder was pressed into a pellet in a pure Ar-filled glove box and heated at 1100-1250 • C for 40-60 h in an evacuated silica tube. For R = Nd, a mixture of NdAs, NdP, Fe 2 O 3 , Fe and FeF 2 powder was used for synthesis. As the solubility limit of fluorine is low ( y < 0.15-0.20), hydrogen was used to synthesize heavily electron-doped samples ( y = 0.20 -0.40). H-substituted samples were synthesized under high pressure. Among the above powders, LaF 3 and FeF 2 were replaced by LaH 2 and NdH 2 for the La-and Nd-systems, respectively. A mixture of all the chemicals with the appropriate stoichiometric ratio was pressed into a pellet and heated at 1100 • C for 2 h under 4 GPa.
All the samples were characterized using high-resolution X-ray diffraction with beam energy of 11.5 keV and 15 keV at room temperature in BL-8A/8B of Photon Factory in KEK, Japan. The lattice constants, pnictogen (Pn) height from the Fe-plane ( h Pn ), Pn-Fe-Pn bond angle ( α ), and Fe-Pn bond distance were calculated from the experimental data by Rietveld analysis 42 .
As, P, and Sb concentrations were determined by energy-dispersive X-ray (EDX) spectroscopy. In all samples, the estimated As-, P-, and Sb-concentrations were almost the same as the nominal values. We could not determine the F and H concentrations by EDX, because F and H are light elements.
Magnetic susceptibility measurements were performed using a Magnetic Property Measurement System (MPMS), with an applied field of 10 Oe. Electrical resistivity was measured using a standard four-probe method. Most of the T c -values presented here were those determined from zero resistivity.