Reciprocity between local moments and collective magnetic excitations in the phase diagram of BaFe2(As1−xPx)2

Unconventional superconductivity arises at the border between the strong coupling regime with local magnetic moments and the weak coupling regime with itinerant electrons, and stems from the physics of criticality that dissects the two. Unveiling the nature of the quasiparticles close to quantum criticality is fundamental to understand the phase diagram of quantum materials. Here, using resonant inelastic x-ray scattering (RIXS) and Fe−Kβ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rm{Fe}}-{{\rm{K}}}_{\beta }$$\end{document} emission spectroscopy (XES), we visualize the coexistence and evolution of local magnetic moments and collective spin excitations across the superconducting dome in isovalently-doped BaFe2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}$$\end{document}(As1−x\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{1-x}$$\end{document}Px\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{x}$$\end{document})2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}$$\end{document} (0.00 ≤\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \le $$\end{document} x ≤0.52\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\le 0.52$$\end{document}). Collective magnetic excitations resolved by RIXS are gradually hardened, whereas XES reveals a strong suppression of the local magnetic moment upon doping. This relationship is captured by an intermediate coupling theory, explicitly accounting for the partially localized and itinerant nature of the electrons in Fe pnictides. Finally, our work identifies a local-itinerant spin fluctuations channel through which the local moments transfer spin excitations to the particle-hole (paramagnons) continuum across the superconducting dome. Understanding the effect of the subtle interplay between strong and weak coupling regimes on the properties of quasiparticles is key to understanding unconventional superconductivity. Here, X-ray spectroscopies reveal the partially localized and itinerant magnetic character of quasiparticles in a doped iron pnictide material.

I t is now well established that unconventional superconductivity (SC) originates from a unique and augmented manifestation of the electronic correlation properties that arises as the system is driven towards a quantum critical region via various tuning parameters, such as doping and pressure [1][2][3][4][5][6][7][8][9][10][11] . On one side of the superconducting dome, the correlation strength is strongly enhanced, leading to diverse quantum manybody effects, such as the Mott insulating state, non-Fermi-liquid behaviour and magnetic orders (see Fig. 1a, b) 3,4,[12][13][14] . On the other side of the dome, the correlation strength is often substantially suppressed, and the low-energy physics can be described by a more conventional Fermi-liquid theory 11 . The transition region between these two limits of correlation holds a quantum critical point (QCP) barely understood. Generally, SC is optimised in this intermediate region where the cooperation of a strongly enhanced non-Fermi-liquid behaviour, and the presence of a quantum criticality leads to unconventional and not understood physics 11 . This generic phase diagram suggests that intertwined electronic and magnetic instabilities, arising from the intermediate correlation strength close to this critical region, induce strong Cooper pairing 3,4,15 .
The physics that determines the evolution of the spin excitations is crucially dependent on the interaction strength of the electrons compared with their bandwidth. We explain this phenomenon schematically in Fig. 1c-g. The non-interacting density of states (DOS) of electrons ( Fig. 1c) is renormalised by the interaction in the two extreme limits as (i) the bands become sharper in the weak-coupling quasiparticle picture (Fig. 1d), and (ii) split into two Mott-like bands characterised by local moments in the strong-coupling limit (Fig. 1e). (i) In the weak-coupling limit, the low-energy spin excitations (χ ii , see Fig. 1d) are very fragile and mix with the particle-hole continuum, failing to form a localised moment. (ii) In the strong-coupling limit, the spin excitations across the Mott bands (χ ll , see Fig. 1e) feature a gapped behaviour on the order of the on-site energy U without the particle-hole continuum 13,14 . However, as the correlation strength is tuned to the intermediate coupling region (Fig. 1f), the correlated DOS co-hosts the quasiparticle DOS at low energy and Mott states at high energy. In this way, a new spin-excitation channel (χ il see Fig. 1g) appears, through which the local moment can now decay to the particle-hole (paramagnon) channels across the magnetic QCP. Fe pnictides (FePns) can be placed in this intermediate region where the interplay of local and itinerant electronic states leads to high-temperature SC 1,2,11,15 . An important open question in this context is how local magnetic moments and collective spin excitations are evolving across the superconducting dome in FePns.
FePns have a layered structure that reduces dimensionality, and the parent compounds exhibit a spin-density wave with collinear antiferromagnetic (AF) order (see Fig. 1b), which gives way to SC via doping as outlined in Fig. 1a [1][2][3][4][5][6][7][8][9][10][11] . The origin of the magnetism is poised between being itinerant, as in Cr, and localised, as in cuprates or heavy fermion materials 1,2,[5][6][7][8][9][10][11]16 producing an uncommon behaviour that has important consequences for the properties of FePns. Theoretical models proposed that the pairing interaction leading to the superconducting phase is provided by residual AF fluctuations persisting upon doping 3,4,15 . Yet, owing to the contribution of the AF fluctuations from local and itinerant states, a complex interplay between them arises, which is believed to play a crucial role in shaping the superconducting dome. Thus, the experimental study of magnetism across the phase diagram of high-temperature superconductors is of vital importance to provide a solid basis for testing these theories.
Inelastic neutron scattering (INS) is the traditional technique of choice to study magnetism, being able to detect magnetic fluctuations in the full Brillouin zone (BZ) 8,9 . Recently, resonant inelastic X-ray scattering (RIXS) has emerged as a complementary technique to INS by detecting spin excitations in FePns close to the Γ point as summarised in Refs. [17][18][19][20][21][22] . The detection of spin excitations is enabled in RIXS thanks to the spin-orbit coupling of the intermediate state mixing the quantum numbers L and S, thereby  1 Phase diagram, structure and excitations of BaFe 2 (As 1−x P x ) 2 . a Phase diagram of BaFe 2 (As 1−x P x ) 2 . The black squares represent the doping levels and temperature measured in our work. As blue line, we schematically depict the expected behaviour of the fluctuating moment. b Schematic real space magnetic structure of BaFe 2 As 2 . c-f Density of states (DOS) as a function of quasiparticle energy. Non-interacting c, weak-coupling random-phase approximation (RPA) d, strong coupling e and intermediate coupling momentum-resolved density-fluctation theory (MRDF) f. g Excitation spectrum of the spin excitations arising from weak, strong and intermediate coupling activating a channel for the detection of magnetic excitations 23,24 . An important consideration when comparing RIXS with INS is the portion of BZ probed by the two techniques, close to the Γ point in the case of RIXS and at the AF wave vector for the case of INS. Depending on the case, these two positions in momentum space can be equivalent or not. Moreover, it is hard to estimate and compare the absolute weight of magnetic excitations in these two regions of BZ, but it is generally accepted that the intensity at the AF wave vector is higher than close to the Γ point.
In this article, we use RIXS to systematically unveil the persistence and gradual hardening of the spin excitations in isovalently doped BaFe 2 (As 1−x P x ) 2 across the phase diagram (see Fig. 1a for a graphical description of the doping levels). Upon doping and without nominally injecting charge carriers, the spectral weight of the spin excitations increases. Our RIXS measurements are complemented with theory placing the Febased superconductors in the intermediate region of correlations and well describing the evolution of the spin excitations with the spin susceptibility (χ s ). In addition, an investigation of the local fluctuating magnetic moment unravels a decrease in the local moment (μ bare ) as a function of doping. This apparent dichotomy implies that there is a transfer of magnetic spectral weight from localised to itinerant as imposed by sum rule relations. We argue that the balance between localised and itinerant states is the key to achieve SC and plays an important role for the physics of criticality in BaFe 2 (As 1−x P x ) 2 . Fig. 2, we show a selection of raw RIXS spectra and fitting of the elastic line, fluorescence background and spin excitations for all the doping levels. The RIXS spectra display a low-energy mode ascribed to spin excitations dispersing as a function of in-plane momentum transfer (q == ) [17][18][19][20][21] . The bandwidth, better appreciated at high momentum transfer (top row in Fig. 2), gradually increases with doping along both (H, 0) and (H, H) as summarised for the highest q points in Fig. 3. The doping increases the width (damping) of the spin excitations, leading to broader dispersing modes at high doping levels. This is expected as the hardening of the spin excitations coupled to electron-hole pair excitations naturally leads to further broadening of the spin excitations 25 . Our main results are presented in Fig. 4a, where we report the doping dependence of the dispersion of the spin excitations extracted from RIXS experiments as dots with error bars, and overlaid on the calculated dynamical spin susceptibility (χ s ). All the doping levels investigated display dispersive spin excitations as shown in Fig. 2 and 4a and in Supplementary Fig. 6. Upon doping the spin excitations, observed in our RIXS experiments, harden in energy by about 45 meV between the parent (x = 0.0) compound and the most doped compound (x = 0.52) at (0.44, 0). The hardening for the same doping levels at (0.32, 0.32) is on the order of 65 meV as summarised in Fig. 3.

RIXS. In
The hardening of the spin excitations has been observed by INS for BaFe 2 (As 0:7 P 0:3 ) 2 26 . At the zone boundary, INS observed a The persistence of spin excitations along the superconducting dome was previously observed in both electron-and hole-doped FePns but with different effects on the energy of these modes. In hole-doped Ba 1Àx K x Fe 2 As 2 , the spin excitations soften upon doping due to the increase in the electronic correlations as demonstrated by theoretical calculations and accurate measurements of the Sommerfeld constant [29][30][31] . In electron-doped FePns, a different behaviour has been detected with the spin excitations being unaffected by doping in their bandwidth but with a decrease in spectral weight 29,32 . In our present work, we uncover that the energy of the spin excitations increases upon isovalent doping and does not directly correlate with the critical temperature (T c ). This indicates that T c is likely connected to other microscopic details aside the effective exchange constant. Thus, the case of isovalent BaFe 2 (As 1−x P x ) 2 reveals an unprecedented behaviour, in stark contrast to hole-and electron-doped BaFe 2 As 2 , adding new features to the large diversity in the doping effects and magnetism in FePns.
The persistence of spin excitations in unconventional superconductors outside the antiferromagnetic parent compound is also a hallmark of the cuprates [33][34][35][36][37] , where paramagnons have been observed across all the phase diagrams. In cuprates, the energy evolution of the spin excitations upon doping is different than the case of FePns: in hole-doped cuprates, there is a more or less constancy of the bandwidth along the antinodal direction, whereas electron-doped cuprates display hardening of the magnetic excitations [35][36][37] .
Momentum-resolved density fluctuation theory. One approach able to describe the regime of intermediate coupling is momentum-resolved density fluctuation (MRDF) theory. MRDF dissects the electronic spectrum into strongly renormalised, incoherent 'local' states, and low-energy, itinerant Bloch states as displayed in Fig. 1f. The back reaction of the quantum fluctuations to the electronic states leads to momentum, energy and orbital-dependent renormalisation as well as lifetime broadening. Within the quantum field theory, these effects are captured by the real and imaginary part of the self-energy correction 38 . MRDF theory self-consistently computes the dynamical correlation functions and the corresponding self-energy in the full momentum and energy space 15,38 . This theory captures the momentum and energy-dependent evolution of the localised states as well as the dispersive quasiparticle states, resulting in an anisotropic pairing symmetry, compatible with what is observed experimentally 3,4 .
To describe the spin excitations probed by RIXS, the dynamical susceptibility (χ) has been calculated and decomposed into the spin (χ s ) and charge (χ c ) channels with both the random-phase approximation (RPA) and MRDF method. The intensity of χ s dominates over χ c in agreement with a previous work 20 . Density functional theory (DFT)-RPA calculations fail to distribute the spectral weight between the local and itinerant states and overestimate the peak energy of the spin excitations compared with our experimental data (see Supplementary Figs. 10 and 11 and Supplementary Notes 1 and 2). The reason for this is that DFT-RPA underestimates the electronic correlation strength placing the spin excitations at too high energy. Moreover, these RPA calculations cannot account for the hardening of the dispersion of spin excitations observed upon doping. In Fig. 1g, we show an idealised scheme of the excitation pattern achieved in the intermediate coupling region comprising local-local (χ ll ), local-itinerant (χ il ) and itinerant-itinerant excitations (χ ii ). The low-energy χ ii excitations have an energy in the order of the spin gap (%1-10 meV) and are not observable by RIXS due to current limitations of the energy resolution but are observed in INS 8,9 . The high-energy χ ll appears at energies of 900-1500 meV and are not detectable by RIXS due to the intense fluorescence background. The χ il excitations are the (para)-magnon excitations and can be qualitatively compared with our RIXS measurements as a function of doping.
In Fig. 4b, c, we depict the doping dependence of χ s directly extracted from the MRDF calculations as well as the RIXS spectra at (0.44, 0). To better visualise the renormalisation in energy of the magnetic excitations, we also take the difference between the selected doping levels and the parent compound and present the results in Fig. 4d for both theory and experiments. The agreement between theory and experiment is remarkable individuating the MRDF, and more in general the intermediate coupling approach as appropriate to describe the magnetism and electronic structure of FePns. In previous works, DFT-RPA has been employed to successfully describe the spin excitations in overdoped cuprates [39][40][41] . The failure of DFT-RPA to account for the spin excitations in BaFe 2 (As 1−x P x ) 2 implies that the FePns are not weakly correlated systems and cannot be compared with overdoped cuprates with reduced electronic correlations but rather need to be placed in the family of multi-orbital correlated systems similarly to heavy fermion materials 16,39 .
X-ray emission spectroscopy. To complement our measurements of the spin excitations and assess the local magnetism of The error bars are defined by the uncertainty of locating the zero-energy position in the raw RIXS spectra. We took conservatively 30 meV (one pixel of our detector), which is much larger than the error resulting from the fitting analysis BaFe 2 (As 1−x P x ) 2 , we employed X-ray emission spectroscopy (XES)-a classical technique that has been established as a sensitive probe of the local magnetic moment (μ bare ) 19,[42][43][44][45][46][47][48] . XES is sensitive to the local fluctuating magnetic moment and does not require a net magnetisation or ordering, such as X-ray magnetic circular dichroism (XMCD) or neutron diffraction, but detects directly the paramagnetic moment [44][45][46] . In this technique, a photon (hν ¼ 7140 eV) excites an Fe-1s core electron into the continuum, creating a core hole that is filled by a Fe-3p electron with the consequent emission of a photon (hν ¼ 7040-7065 eV). The final Fe-3p 5 state has a wave function partly overlapped with the Fe-3d orbitals, and is, thereby, affected by the spin of the valence band through the exchange interaction 19,[42][43][44][45][47][48][49] . The femtosecond timescale of this technique allows for the measurement of fast fluctuations of μ bare 5,42,43,49 , and by means of normalisation and calibration [carried out with FeCrAs (0 μ B ) and BaFe 2 As 2 (set arbitrarily to one to define a relative scale)] μ bare can be determined 42,43,49 .
In Fig. 5a, b, we show XES spectra for BaFe 2 As 2 , and BaFe 2 (As 1−x P x ) 2 (x = 0.52), and FeCrAs and the respective difference spectra. Clearly, a gradual decrease in μ bare is inferred from the difference spectra depicted in Fig. 5c. The values of μ bare , presented in Fig. 5d, e are continuously reduced, (e.g., μ bare ¼ 1.0 ± 0:1 for x ¼ 0:00, μ bare ¼ 0:6 ± 0:1 for x ¼ 0:28 and μ bare ¼ 0:4 ± 0:1 for x ¼ 0:52), but not fully quenched by doping, despite the complete disappearance of the ordered magnetic moment observed by neutron scattering 50,51 . This evidence is remarkable in light of the constancy of the Fe oxidation state observed in XAS (see Supplementary Fig. 2) with isovalent doping driving the antiferromagnetically long-range ordered BaFe 2 As 2 into a paramagnetic phase.
We compare the doping dependence of the local moment with the strength of the Fermi surface (FS) nesting at the antiferromagnetic wavevector. In Fig. 5d, the theoretical data are presented as Δχ(x) = χ(x) -χ(0.6), where χ(x) is the computed static susceptibility at the antiferromagnetic wavevector at doping x [we normalise Δχ (x = 0) = 1]. Note that this comparison is only qualitative, and a self-consistent estimation of the static magnetic moment is computationally expensive. Similar doping dependence of the magnetic instability arising from the FS nesting and the observed local moment indicates an electronic mechanism of the magnetic ground state in this system.

Discussion
Our XES results are puzzling, and assuming spectral weight sum rules, there has to be spectral weight transfer from localised moments into spin excitations as probed by RIXS. In Fig. 5e, we report the integrated spectral weight of the spin excitations determined from our RIXS spectra. Within error bars, the spin excitation intensities detected in the portion of BZ pertinent to RIXS possibly show gradual transfer of spectral weight upon doping. However, additional enhancement of spectral weight in  This spectral redistribution signals that SC is optimised by a balance of these types of interactions. In summary, we unveiled the persistence and hardening of magnetic excitations along the superconducting dome of BaFe 2 (As 1−x P x ) 2 with a concomitant decrease in local magnetic moments. The spectral weight of the spin excitations observed In order to visualise the reciprocal redistribution of local magnetic moments and spin excitations, we represent the superconducting dome as a purple shaded area in the background. The error bars on the intensity of the spin excitations have been quantified in ± 0.1 with respect to the value of the parent compound. The fitting error is actually smaller, but the sum procedure at different q points produces a propagation of the errors, so we set a slightly higher upper boundary with RIXS slightly increases, corroborating that SC emerges from a balance of localised and itinerant electronic interactions. The experimental results agree well with intermediate coupling calculations placing the FePns in the family of multi-orbital correlated systems and uncovering clearly the mixed electronic and magnetic interactions present with optimal SC.

Methods
Sample preparations. Single crystals of BaFe 2 (As 1−x P x ) 2 were grown either by stoichiometric melt 52 or Ba 2 As 3 /Ba 2 P 3 self-flux method 53 . The samples were characterised either with resistivity or magnetisation measurements. For x ¼ 0:15, 0.22, 0.27 and 0.48, crystals grown by the stoichiometric melt technique were used for the RIXS experiments, while for x ¼ 0:38 and 0.52 we use the crystals grown by the self-flux technique. Supplementary Fig. 1 shows the temperature dependence of the resistivity on crystals selected from the same batches. For x > 0:3 the spindensity wave (SDW) long-range order is completely suppressed 52,53 .
X-ray absorption spectroscopy and RIXS. X-ray absorption spectra (XAS) and RIXS experiments were performed at the ADRESS beamline of the Swiss Light Source, Paul Scherrer Institut, Villigen PSI, Switzerland 54,55 . In agreement with previous works [17][18][19][20]56 , the samples were mounted with the ab plane perpendicular to the scattering plane and the c axis lying on it (sketch in Supplementary Fig. 2) and post cleaved in situ at a pressure better than 2.0 × 10 −10 mbar. The reciprocal space directions studied are (0, 0)!(1, 0) and (0, 0)!(1, 1) according to the orthorhombic unfolded crystallographic notation. The values of in-plane momentum transferred are expressed as relative lattice units (RLU) (q == ¼ 2π/a). We use the convention of 1 Fe per unit cell. All the measurements were carried out at 13 K by cooling the manipulator with liquid helium. XAS spectra were measured in both total fluorescence yield (TFY) and total electron yield (TEY) modes. No difference was observed between TEY and TFY, which indicates the sample integrity across the thickness probed in our experiment. We measured Fe L 2;3 XAS spectra for all samples at 15°of incidence angle relative to the sample surface. All the XAS spectra are reported in Supplementary Fig. 3 and display the constancy of the iron oxidation state. There are small spectral differences at around 710 eV that are possibly due to a different covalency between the FeAs and FeP. The RIXS spectrometer was set to a scattering angle of 130°, and the incidence angle on the samples surface was varied to change the in-plane momentumtransferred (q == ) from (0, 0) to (0.44, 0) and from (0, 0) to (0.31, 0.31). All RIXS measurements in this paper were recorded in grazing-incidence configuration as depicted in Supplementary Fig. 2. The zero-energy loss of our RIXS spectra was determined by measuring spectra in σ polarisation. The total energy resolution was measured by employing the elastic scattering of carbon-filled acrylic tape and is around 110 meV.
RIXS spectra were normalised to unity, and the main emission line was fitted according to Refs. [17][18][19][20]56 employing the following formulas: with and In the first formula, the first part is a second-order polynomial function describing the emission line at low-energy loss, the second part is an exponential decay describing the emission line at high-energy loss. The two behaviours are swapped into each other by the g γ term. The third term is a Gaussian curve (G) observed at around -4.2-eV energy loss. An exemplary fitting of the full RIXS spectra is shown in Supplementary Fig. 4 for BaFe 2 As 2 at (0.44, 0.0). In Supplementary Fig. 5, we show RIXS spectra over an extended range of energy loss. Overall, the three spectra resemble a metallic system with the lineshape not being affected by doping. Due to changes in the band structure upon doping, small energy shifts and changes in the RIXS intensity are detected.
K β emission spectroscopy. We performed the Fe-K β X-ray emission (XES) experiments at BL11XU of SPring-8, Japan. The incoming beam was monochromatised by a Si (111) double crystal and a Si (400) secondary channel-cut crystal. The energy was calibrated by measuring X-ray absorption of a Fe foil and set to 7.140 keV with π polarisation. We employed three spherical diced Ge (620) analysers and a detector in Rowland geometry at 2-m distance. The total combined resolution was about 400 meV estimated from full width half maximum (FWHM) of the elastic line.
We scanned the absolute emission energy between 7.02 keV and 7.08 keV and normalised the intensity by the incident flux monitored by an ionisation chamber. We carried out the measurements at 15 K. The experimental geometry of the experiment is depicted in Supplementary Fig. 7. The determination of μ bare was performed by employing the integrated absolute difference (IAD) method described by Vanko et al. 45 . The areas of the XES spectra were normalised to unity as shown in Supplementary Fig. 8a-f. To overcome problems of alignment of the absolute energy between different samples, we estimated the centre of mass of the spectra and aligned the energy to have the same centre of mass as described by Glatzel et al. 44 . Then the difference with the spectrum of the FeCrAs reference was taken as displayed in Supplementary Fig. 8a-f and summarised in Supplementary Fig. 8g. The integration of this difference gives the IAD, which is directly proportional to μ bare 45 .
Doping-dependent band structure. Two extreme doping regimes of the BaFe 2 (As 1−x P x ) 2 samples, namely BaFe 2 As 2 (x ¼ 0:00) and BaFe 2 P 2 (x ¼ 1:00), have been studied within DFT 57 . Also accurate tight-binding models are available to reproduce the DFT band structure, consistent with angle-resolved photoemission spectroscopy (ARPES) data (after including renormalisation effects). However, P-doping effects on the As site, which do not change the carrier concentration, are neither trivial to accurately calculate within the DFT framework, nor within a simple rigid band shift technique that works reasonably well for electron-and holedoping cases 58 . Since P atom is smaller in size compared with As atoms, P doping is expected to decrease the lattice volume. Indeed, both X-ray 52 11 . Among all the FS pockets, the outermost hole pocket shows most dominant and drastic changes with doping: it is very close to cylindrical as a function of k z at x ¼ 0:00, but gradually becomes corrugated with increasing doping. The other pockets acquire comparatively less k z dispersion [59][60][61] . This could be quantified by calculating the nesting strength via the integration of the static susceptibility at the nesting wave vector. The nesting strength has been shown to decrease upon doping 11,58,60 . Projected orbital symmetries on each FS reveal that the hole pockets are dominated by the d xz /d yz orbitals, while the electron pockets are made of hybridisation between the d xz /d yz orbitals with the d xy orbitals. The d z 2 orbital has little contribution at x ¼ 0:00, while it grows with doping. d z 2 basically contributes to the strongly k z dispersive parts of the outermost hole pocket 60 . Based on the aforementioned experimental observations on the band structure evolution with doping, we construct an effective five-orbital tight-binding (TB) model to reproduce the low-energy dispersion and FS topology across the entire phase diagram. The 2 Fe unit cell incorporates two Fe sublattices producing a tenband model. We start with a five-orbital TB model as derived in Ref. 57 for BaFe 2 As 2 . The intra-orbital dispersions are defined by ξ i k ¼ ϵ i k þ Δ i À ϵ F , where k is the crystal momentum and i ¼ 1 À 5 is the orbital index. ϵ i k is the momentumdependent part of the dispersion that arises from the nearest, next and higher neighbour hoppings between the same orbitals, Δ i is the corresponding on-site potentials and ϵ F is the chemical potential. As the unit cell volume decreases monotonically with P doping, it is expected that the electron hopping amplitude increases monotonically with doping. We model this effect by a simple renormalisation factor λ as ϵ i k ! λϵ i k , where λ increases with doping. Similarly, due to the monotonic increase in the pnictogen coordinate and the height (z Pn , h Pn ), the on-site potential also changes. Interestingly, we find that an orbital-dependent modification of the on-site potential is required to properly reproduce the experimental behaviour of the k z dispersion and FS changes. We set Δ i ! Δ i þ δ for d xz /d yz orbitals, Δ i ! Δ i À δ for d x 2 Ày 2 and d 2 z orbitals and Δ i ! Δ i À 2δ for d xy orbital, where δðxÞ ¼ 0:08x, with x being the doping value. The chemical potential is adjusted at each doping to keep the same number of electrons. Such a model was also invoked earlier to model the FS topological transition in KFe 2 Se 2 under pressure 62 . The corresponding FS topologies at three representative dopings are shown in Supplementary Fig. 9.

Data availability
Data that support the findings of this study are available upon reasonable request from the corresponding authors.