# Non-perturbative series expansion of Green’s functions: The Anatomy of

Resonant Inelastic X-Ray Scattering in Doped Hubbard Model

###### Abstract

We present a non-perturbative, divergence-free series expansion of Green’s functions using effective operators. The method is especially suited for computing correlators of complex operators as a series of correlation functions of simpler forms. We apply the method to study low-energy excitations in resonant inelastic x-ray scattering (RIXS) in doped one- and two-dimensional single-band Hubbard models. The RIXS operator is expanded into polynomials of spin, density, and current operators weighted by fundamental x-ray spectral functions. These operators couple to different polarization channels resulting in simple selection rules. The incident photon energy dependent coefficients help to pinpoint main RIXS contributions from different degrees of freedom. We show in particular that, with parameters pertaining to cuprate superconductors, local spin excitation dominates the RIXS spectral weight over a wide doping range in the cross-polarization channel.

###### pacs:

Our understanding in many aspects of physics has benefited tremendously from perturbation theory. One of its triumphs lies in the practical calculation of various scattering cross sections and correlation functions in many branches of physics ranging from quantum electrodynamics to condensed matter. Although widely successful, such perturbative approaches require the “perturbing part” to be small to ensure convergence. This is not always a given, especially for systems that involve different degrees of freedom with comparable energy scales. The divergence of series expansions in quantum electrodynamics has been discussed early on Dyson (1952). Fuelled by recent advances in numerical diagrammatic methods aimed to tackle the many-body problem in regimes with competing interactions Toschi et al. (2007); Rubtsov et al. (2008); Kozik et al. (2010); Rohringer et al. (2013); Ayral and Parcollet (2015); Gull et al. (2011), it has been realized that divergent series need to be treated with care Gull et al. (2011); Lani et al. (2012); Schäfer et al. (2013); Janiš and Pokorný (2014); Kozik et al. (2015); Stan et al. (2015); Ribic et al. (2016); Rossi et al. (2016); Tarantino et al. (2017); Gunnarsson et al. (2017). We here present an alternative expansion scheme for the calculation of correlation functions that does not depend on small parameters and converges for any Hamiltonian. The difference between our projection method and standard perturbation theory is much like the difference between a Fourier series that can describe functions with poles and a Taylor series with finite range of convergence Sup ().

We apply our theory to describe the transition process of RIXS, which has emerged in recent years as a versatile tool for studying the energy-momentum structure of charge, spin, orbital, and lattice excitations in solids Ament et al. (2011). RIXS is a second-order process, described by a four-point two-particle Green’s function. In the first step, a core electron is promoted into an empty valence state by absorbing a high-energy photon. The created core hole interacts with its environment, and eventually decays into an excited state after certain lifetime by emitting a photon. Detailed information of various elementary excitations is then encoded in the change of energy, momentum, and polarization between the incident and scattered photons. The high photon flux provided by modern synchrotron x-ray sources and its energy-momentum window make RIXS an appealing complement to more established methods such as inelastic neutron scattering (INS). There are clear parallels between the two—both processes can be cast into the form of a momentum space correlation function

(1) |

While the results of INS can be directly communicated as the dynamic spin correlation function by identifying the operator with the spin operator , the complex structure of the RIXS operator originated from its non-trivial intermediate state prevents a direct interpretation of the measured spectra.

One recognizes the challenges in finding an effective theory for RIXS when realizing that the intermediate state in RIXS is the final state in x-ray absorption spectroscopy (XAS). For XAS, it is well known that core-valence interactions lead to excitons or resonances with asymmetric line-shapes, e.g. edge singularities Anderson (1967); Zaanen et al. (1985); Benjamin et al. (2013), combined with multiplets and charge-transfer shake-up excitations Haverkort et al. (2014); De Groot and Kotani (2008). No small parameters are present in the intermediate state, and previously proposed series expansions based on the assumption that the core-hole lifetime is small Luo et al. (1993); van den Brink and van Veenendaal (2006); Ament et al. (2007) turned out to be non-convergent Jia et al. (2016); Ament and van den Brink (2010). The lack of concrete understanding of RIXS cross section greatly limits its potential as a true alternative to INS. One acute example is the ongoing debate on the nature of the observed dispersing low-energy RIXS feature in overdoped cuprates Le Tacon et al. (2011, 2013); Dean et al. (2013); Ishii et al. (2014); Lee et al. (2014); Minola et al. (2015); Huang et al. (2016); Minola et al. (); Jia et al. (2014); Benjamin et al. (2014); Kanász-Nagy et al. (2016). On one hand, it was interpreted as damped collective magnetic excitations Le Tacon et al. (2011, 2013); Dean et al. (2013); Ishii et al. (2014); Lee et al. (2014); Minola et al. (2015, ); Jia et al. (2014); Huang et al. (2016), supported empirically by its similar momentum-space structure to that of dynamic spin structure factor Le Tacon et al. (2011); Jia et al. (2014). On the other hand, it was argued to be originated from incoherent particle-hole excitations Benjamin et al. (2014); Kanász-Nagy et al. (2016), as canonical Fermi-liquid behavior is expected at high doping levels.

In this paper, we present an alternative expansion method that does not depend on small parameters and is free from divergences. The RIXS operator is expanded into a series of effective operators composed of spin, density, and current operators in the proximity of the core-hole site. Good agreement with the exact solution can be reached within the first few orders, whereby the exact rate of convergence is model and material dependent. This approach provides an unbiased survey of all possible low-energy excitations that couple to the RIXS process with explicit considerations of light polarizations and incident photon energy. The result is an exact mapping of the scattering cross section to a set of intrinsic correlation functions. Such a procedure is generally applicable for a wide class of models.

Following Ref. Haverkort (2010), we start by rewriting Eq. (1), whose imaginary part gives the RIXS intensity as

(2) |

where denotes the initial state of the system with energy determined by the Hamiltonian . and are the momentum and energy transfer between the incoming () and outgoing () photons with polarizations and , respectively. The -dependent RIXS operator is obtained by a Fourier transform of the local operator

(3) |

which describes the RIXS process where a core hole is created at site by a dipole transition operator and subsequently annihilated locally by after some lifetime . A generally valid assumption is imposed here that the core hole is non-dispersive, as the overlap of core-electron wave-functions is usually negligible between different sites.

To understand what excitations can be probed by RIXS, it is favorable to express as a sum of effective operators that will consequently turn Eq. (2) into a series of simpler correlation functions as

(4) |

where the operator brings a state to . In cases where and are states in a complete orthonormal Hilbert space, the expansion is trivial and is readily given as Sup (). Expansion over a complete basis set is not practical due to the exponentially large size of the Hilbert space. We want to expand on operators that couple exponentially many states sharing certain property to exponentially many states having property in common. Two sufficient conditions for obtaining the expansion coefficients in Eq. (4) for such operators are:

(5) |

Using these relations we obtain for any given wave function by multiplying Eq. (4) on the right by and left by ,

(6) |

from which it follows that

(7) |

where we used the freedom to choose , and .

For RIXS (or any core-level spectroscopy in general), one should bear in mind its local nature defined by the immobile core hole and its short lifetime, which means that the expansion of will have a good convergence rate by expanding over operators that act in the proximity of the core-hole site . Let us define the Fock space of a subsystem , where denotes the core-hole site () or its up to -th-nearest neighbors (). is the single-site Fock space at site . Operators satisfying Eq. (5) can then be defined as , where and are orthonormal basis states of and is the identity operator acting on .

In the following we consider the single-band Hubbard model, although we note that generalization to multi-orbital case is straightforward. A multi-orbital expansion restricted to local spin operators can be found in Ref. Haverkort (2010). The single-site Fock space in this case is spanned by states with zero, single and double occupations {, , , }. In the zeroth-order () expansion, there are in total operators , among which only particle-number-conserving ones couple to the RIXS process, namely

(8) | ||||||

which are the linear combinations of local spin operators , , and density operators , , and defining the complete set of on-site excitations. Among the other 10 operators, 8 (2) change the particle number by 1 (2), which can combine with operators on neighboring sites to form number-conserving ones in higher-order expansions. An expansion including two (three) sites generates 256 (4096) operators, and 70 (924) of them are number conserving Sup (). Although this is a rapidly growing series, the expansion converges within the first few orders ^{1}^{1}1The Fermi velocity is m/s in a good metal, thus the distance an electron can travel in the intermediate state is Å, corresponding to only a few lattice spacings. See also Ref. Savary and Senthil ().
In addition, the number of operators can be further reduced by considering the local symmetry and light polarizations. For systems with local inversion symmetry, only has non-zero coefficient in the expansion for the cross-polarization channel ( and assuming ), or and for the parallel-polarization () one.

The Hamiltonian of the model we study reads

(9) |

which is defined on a one-dimensional (1D) chain or a two-dimensional (2D) square lattice. The first line of Eq. (Non-perturbative series expansion of Green’s functions: The Anatomy of Resonant Inelastic X-Ray Scattering in Doped Hubbard Model) describes the hoppings of the “”-electrons between nearest () and next-nearest () neighboring sites with local Coulomb interaction , and the second line accounts for the repulsion between the states and the core states “” with spin-orbit coupling . The numerical results will be evaluated with parameters given in the Supplemental Material Sup (). For both 1D and 2D, the calculation is performed on 12-site clusters with periodic boundary conditions using the many-body package Quanty Haverkort (2016); Lu et al. (2014).

We first show the results for the 1D case. The exact -integrated RIXS cross sections in the parallel-polarization channel for the hole-doped case with occupation are plotted in Fig. 1, in comparison to the approximate ones with effective operators expanded up to next-nearest neighbors (). The incident photon energy is tuned to the maximum of the XAS spectra (“resonance”) Sup () for Fig. 1(a) and also “detuned” at eV higher for Fig. 1(b). All spectra are broadened by a Lorentzian function with full width at half maximum of for plotting.

Within the expansion, regardless of doping and incident photon energy (see Fig. 1 and Supplemental Material Sup ()), the approximate cross sections only resemble the exact ones on the higher-energy end, while the low-energy spectral weight is suppressed. This highlights the importance of non-local charge fluctuations as pointed out by earlier studies Benjamin et al. (2014); Jia et al. (2016). The missing spectral weight can be largely restored by including the neighboring sites in the expansion, and the approximation reproduces almost fully the exact solution. The convergence rate can be quantified by comparing the spectral moments between the exact and approximate RIXS cross sections. A good convergence of the spectral weight and its center of mass is reached for with different values both before and after the resonance Sup ().

To address the relation between RIXS and the charge response function, we performed a partial expansion in the subspace using only the density operator . The resultant dynamic charge structure factor is shown in Fig. 1(a). The large discrepancy between the RIXS cross section and suggests that a direct association between the two should be discouraged. Their fundamental difference is not difficult to understand. In the presence of large on-site repulsion, the local charge fluctuation is associated with an energy scale of and thus strongly suppressed. The RIXS process, on the other hand, involves other low-energy excitations as shown by the additional effective operators such as ^{2}^{2}2The subscripts are the site indices, where 0 denotes the core-hole site, and the bonding combination of the neighborings sites with lattice spacings to the core-hole site., which corresponds to a two-spin excitation and contributes to the lower-energy RIXS spectral weight.

Figure 2 shows the 1D exact and approximate RIXS cross sections in the cross-polarization channel. Similar to the parallel-polarization results, an overall improvement of the approximation can be achieved by extending the expansion range. Contrary to the former case where the inclusion of non-local excitation is essential, the majority of the spectral weight is already captured here at all dopings by a local expansion Sup (). This observation indicates that the RIXS spectrum consists mainly of dynamic spin structure factor not only in undoped magnetic insulators Haverkort (2010) but also up to a high doping range. This seemingly surprising result can be explained by examining the relevant non-local effective operators. The cross-polarization channel couples to those operators that break inversion symmetry such as a spin-pair exchange between the core-hole site and the -th-nearest-neighbor sites. Such processes are strongly prohibited compared to a local spin excitation when large local repulsion is present. Nonetheless, these non-local contributions can give rise to a nontrivial dependence of the RIXS spectra. While the energy of local spin excitations is independent, as shown by the constant values of spectra [or equivalently ] Sup (), the non-local spin-flip process will move to higher energy with increasing . The difference of dependence is crucial for pinpointing the nature of the experimentally observed low-energy feature of the cross-polarization RIXS in doped cuprates Benjamin et al. (2014); Minola et al. (2015); Huang et al. (2016); Minola et al. ().

Figure 3(a) shows the coefficient for and for . While the resonant energy of coincides with the XAS, that of is located at eV higher. The relative weight of non-local excitations thus increases by detuning the incident photon to higher energies. Therefore, while the dominant spectral weight in the cross-polarization channel still originates from the local spin excitations at moderate doping levels (Fig. 2), the relative increase of non-local contributions at high detuning energies may give rise to a fluorescence-like behavior Sup (). Figs. 3(b) and 3(c) show the expansion coefficients of operators and that appear in higher-order expansions. The near exponential convergence rate as a function of lattice spacing confirms the assumption that the most important excitations happen locally around the core-hole site.

In the last part we address the relation between the cross-polarization RIXS and in hole-doped 2D Hubbard model using the presented method. Figure 4 shows the -integrated and RIXS cross sections together with obtained by expansion. The higher order expansions are shown in the supplemental material Sup (). We find that the convergence rate for the 2D Hubbard model is slightly slower than for the 1D case. At half filling, as shown in Fig. 4(a) and (b), the cross-polarised RIXS probes nearly exact the local spin excitations at low energy, as evidenced by their almost identical spectral weight and line shape. Upon hole doping, while the discrepancy between the two increases due to growing non-local charge and spin fluctuations, the RIXS spectral weight originated from local excitations remains dominant. Even for highly overdoped case [Fig. 4(c) and (f)], which was deemed to be fully describable by Fermi liquid theory Benjamin et al. (2014), the local excitations still constitute about 60 percent of the total spectral weight, in line with recent experiments Minola et al. (). This finding supports earlier “paramagnon” interpretations Le Tacon et al. (2011, 2013); Minola et al. (2015) of the low-energy feature in cross-polarization RIXS on doped cuprates and suggests that considerations based purely on quasiparticles cannot account for the full spectral weight.

In summary, we presented a method to generate a divergence-free series expansion of Green’s functions using effective operators. We applied the method to expand the RIXS operator into a sum of spin, density, and current operators. The result is an exact mapping of the scattering cross section to a set of intrinsic correlation functions, independent of the model and parameters used. The coefficients of the effective operators encode the energy and polarization dependence of RIXS and help to identify the main excitations contributing to the RIXS spectral weight. A quantitative connection between the RIXS cross section and intrinsic correlation functions is provided. Using realistic models tailored to specific materials may help to resolve confusions in the understanding of current measurements and guide future experimental works.

We thank L. J. P. Ament, T. P. Devereaux and K. Wohlfeld for discussion.

## References

- Dyson (1952) F. J. Dyson, Phys. Rev. 85, 631 (1952).
- Toschi et al. (2007) A. Toschi, A. A. Katanin, and K. Held, Phys. Rev. B 75, 045118 (2007).
- Rubtsov et al. (2008) A. N. Rubtsov, M. I. Katsnelson, and A. I. Lichtenstein, Phys. Rev. B 77, 033101 (2008).
- Kozik et al. (2010) E. Kozik, K. V. Houcke, E. Gull, L. Pollet, N. Prokof’ev, B. Svistunov, and M. Troyer, EPL 90, 10004 (2010).
- Rohringer et al. (2013) G. Rohringer, A. Toschi, H. Hafermann, K. Held, V. I. Anisimov, and A. A. Katanin, Phys. Rev. B 88, 115112 (2013).
- Ayral and Parcollet (2015) T. Ayral and O. Parcollet, Phys. Rev. B 92, 115109 (2015).
- Gull et al. (2011) E. Gull, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, M. Troyer, and P. Werner, Rev. Mod. Phys. 83, 349 (2011).
- Lani et al. (2012) G. Lani, P. Romaniello, and L. Reining, New Journal of Physics 14, 013056 (2012).
- Schäfer et al. (2013) T. Schäfer, G. Rohringer, O. Gunnarsson, S. Ciuchi, G. Sangiovanni, and A. Toschi, Phys. Rev. Lett. 110, 246405 (2013).
- Janiš and Pokorný (2014) V. Janiš and V. Pokorný, Phys. Rev. B 90, 045143 (2014).
- Kozik et al. (2015) E. Kozik, M. Ferrero, and A. Georges, Phys. Rev. Lett. 114, 156402 (2015).
- Stan et al. (2015) A. Stan, P. Romaniello, S. Rigamonti, L. Reining, and J. A. Berger, New J. Phys. 17, 093045 (2015).
- Ribic et al. (2016) T. Ribic, G. Rohringer, and K. Held, Phys. Rev. B 93, 195105 (2016).
- Rossi et al. (2016) R. Rossi, F. Werner, N. Prokof’ev, and B. Svistunov, Phys. Rev. B 93, 161102 (2016).
- Tarantino et al. (2017) W. Tarantino, P. Romaniello, J. A. Berger, and L. Reining, Phys. Rev. B 96, 045124 (2017).
- Gunnarsson et al. (2017) O. Gunnarsson, G. Rohringer, T. Schäfer, G. Sangiovanni, and A. Toschi, Phys. Rev. Lett. 119, 056402 (2017).
- (17) See appended Supplemental Material for the details of parameter choice, coefficients of expansions, and numerical results for other doping levels.
- Ament et al. (2011) L. J. P. Ament, M. van Veenendaal, T. P. Devereaux, J. P. Hill, and J. van den Brink, Rev. Mod. Phys. 83, 705 (2011).
- Anderson (1967) P. W. Anderson, Phys. Rev. Lett. 18, 1049 (1967).
- Zaanen et al. (1985) J. Zaanen, G. A. Sawatzky, J. Fink, W. Speier, and J. C. Fuggle, Phys. Rev. B 32, 4905 (1985).
- Benjamin et al. (2013) D. Benjamin, D. Abanin, P. Abbamonte, and E. Demler, Phys. Rev. Lett. 110, 137002 (2013).
- Haverkort et al. (2014) M. W. Haverkort, G. Sangiovanni, P. Hansmann, A. Toschi, Y. Lu, and S. Macke, EPL 108, 57004 (2014).
- De Groot and Kotani (2008) F. De Groot and A. Kotani, Core Level Spectroscopy of Solids (CRC Press, Boca Raton, FL, 2008).
- Luo et al. (1993) J. Luo, G. T. Trammell, and J. P. Hannon, Phys. Rev. Lett. 71, 287 (1993).
- van den Brink and van Veenendaal (2006) J. van den Brink and M. van Veenendaal, EPL 73, 121 (2006).
- Ament et al. (2007) L. J. P. Ament, F. Forte, and J. van den Brink, Phys. Rev. B 75, 115118 (2007).
- Jia et al. (2016) C. Jia, K. Wohlfeld, Y. Wang, B. Moritz, and T. P. Devereaux, Phys. Rev. X 6, 021020 (2016).
- Ament and van den Brink (2010) L. J. P. Ament and J. van den Brink, (2010), arXiv:1002.3773 .
- Le Tacon et al. (2011) M. Le Tacon, G. Ghiringhelli, J. Chaloupka, M. M. Sala, V. Hinkov, M. W. Haverkort, M. Minola, M. Bakr, K. J. Zhou, S. Blanco-Canosa, C. Monney, Y. T. Song, G. L. Sun, C. T. Lin, G. M. De Luca, M. Salluzzo, G. Khaliullin, T. Schmitt, L. Braicovich, and B. Keimer, Nat. Phys. 7, 725 (2011).
- Le Tacon et al. (2013) M. Le Tacon, M. Minola, D. C. Peets, M. Moretti Sala, S. Blanco-Canosa, V. Hinkov, R. Liang, D. A. Bonn, W. N. Hardy, C. T. Lin, T. Schmitt, L. Braicovich, G. Ghiringhelli, and B. Keimer, Phys. Rev. B 88, 020501 (2013).
- Dean et al. (2013) M. P. M. Dean, G. Dellea, R. S. Springell, F. Yakhou-Harris, K. Kummer, N. B. Brookes, X. Liu, Y.-J. Sun, J. Strle, T. Schmitt, et al., Nat. Mater. 12, 1019 (2013).
- Ishii et al. (2014) K. Ishii, M. Fujita, T. Sasaki, M. Minola, G. Dellea, C. Mazzoli, K. Kummer, G. Ghiringhelli, L. Braicovich, T. Tohyama, K. Tsutsumi, K. Sato, R. Kajimoto, K. Ikeuchi, K. Yamada, M. Yoshida, M. Kurooka, and J. Mizuki, Nat. Commun. 5 (2014).
- Lee et al. (2014) W. S. Lee, J. J. Lee, E. A. Nowadnick, S. Gerber, W. Tabis, S. W. Huang, V. N. Strocov, E. M. Motoyama, G. Yu, B. Moritz, H. Y. Huang, R. P. Wang, Y. B. Huang, W. B. Wu, C. T. Chen, D. J. Huang, M. Greven, T. Schmitt, Z. X. Shen, and T. P. Devereaux, Nat. Phys. 10, 883 (2014).
- Minola et al. (2015) M. Minola, G. Dellea, H. Gretarsson, Y. Y. Peng, Y. Lu, J. Porras, T. Loew, F. Yakhou, N. B. Brookes, Y. B. Huang, J. Pelliciari, T. Schmitt, G. Ghiringhelli, B. Keimer, L. Braicovich, and M. Le Tacon, Phys. Rev. Lett. 114, 217003 (2015).
- Huang et al. (2016) H. Y. Huang, C. J. Jia, Z. Y. Chen, K. Wohlfeld, B. Moritz, T. P. Devereaux, W. B. Wu, J. Okamoto, W. S. Lee, M. Hashimoto, Y. He, Z. X. Shen, Y. Yoshida, H. Eisaki, C. Y. Mou, C. T. Chen, and D. J. Huang, Sci. Rep. 6, 19657 (2016).
- (36) M. Minola, Y. Lu, Y. Y. Peng, G. Dellea, H. Gretarsson, M. W. Haverkort, Y. Ding, X. Sun, X. J. Zhou, D. C. Peets, L. Chauviere, P. Dosanjh, D. A. Bonn, R. Liang, A. Damascelli, N. B. Brookes, F. Yakhou, J. Pelliciari, M. Dantz, X. Lu, T. Schmitt, L. Braicovich, G. Ghiringhelli, B. Keimer, and M. L. Tacon, Phys. Rev. Lett. (accepted) arXiv:1708.00663 .
- Jia et al. (2014) C. J. Jia, E. A. Nowadnick, K. Wohlfeld, Y. F. Kung, C. C. Chen, S. Johnston, T. Tohyama, B. Moritz, and T. P. Devereaux, Nat. Commun. 5, 3314 EP (2014).
- Benjamin et al. (2014) D. Benjamin, I. Klich, and E. Demler, Phys. Rev. Lett. 112, 247002 (2014).
- Kanász-Nagy et al. (2016) M. Kanász-Nagy, Y. Shi, I. Klich, and E. A. Demler, Phys. Rev. B 94, 165127 (2016).
- Haverkort (2010) M. W. Haverkort, Phys. Rev. Lett. 105, 167404 (2010).
- (41) The Fermi velocity is m/s in a good metal, thus the distance an electron can travel in the intermediate state is Å, corresponding to only a few lattice spacings. See also Ref. Savary and Senthil ().
- Haverkort (2016) M. W. Haverkort, JPCS 712, 012001 (2016).
- Lu et al. (2014) Y. Lu, M. Höppner, O. Gunnarsson, and M. W. Haverkort, Phys. Rev. B 90, 085102 (2014).
- (44) The subscripts are the site indices, where 0 denotes the core-hole site, and the bonding combination of the neighborings sites with lattice spacings to the core-hole site.
- (45) L. Savary and T. Senthil, arXiv:1506.04752 .

Supplemental Material for

“Non-perturbative series expansion of Green’s functions: The Anatomy of

Resonant Inelastic X-Ray Scattering in Doped Hubbard Model”

## Appendix A Convergence of series expansions

Perturbative solutions of the eigenvalues, eigenvectors, or Green’s functions for a given Hamiltonian written as the sum of a bare Hamiltonian and perturbing Hamiltonian can be written as a power series in the perturbing Hamiltonian with expansion coefficients given as derivatives of with respect to . This is a feature perturbative methods have in common with a Taylor series expansion of a function. The correspondence between a Taylor series and more general perturbative expansions implies that similar convergence criteria apply. For a realistic interacting Hamiltonian, perturbative solutions such as diagrammatic series expansions (e.g. Dyson series) have a tendency to diverge, even with moderate interaction strength. In this letter we provide a non-divergent alternative. Instead of creating a power series in the perturbing interaction, we use a projection method to expand the Green’s functions. This is analogous to approximating a function with a Fourier series instead of a Taylor expansion. The convergence criteria of the former are in general more easily satisfied. In particular, poles, discontinuities or divergences in the complex plane of the original function do not pose problems for the Fourier expansion.

In the following two subsections we discuss the convergence of our projection method for operators. In the first subsection we illustrate the impact of a divergent series by comparing the Taylor and Fourier expansions of the function . In the second subsection we prove that the projection method we propose in the main text leads to a convergent series.

### a.1 Taylor v.s. Fourier expansion of functions

To illustrate the difference between a divergent and a convergent series expansion, we compare the Taylor and Fourier expansions of the function

(S1) |

The function is shown in the left panel of Fig. S1. It is continues and infinitely differentiable on the real axis. Being an even function, only even-order terms arise in its expansion. The Taylor and Fourier series are given as:

(S2) | |||||

with the expansion coefficients

(S3) |

and

(S4) |

As has poles in the complex plane at , the Taylor series only converge for . This can be seen in the middle panel of Fig. S1, where the series is truncated at and . The expansion does not approximate well to the original function with increasing outside of the convergence range . A Fourier series, on the other hand, nicely converges to the original function, as shown in the right panel of Fig. S1.

### a.2 Convergence of the effective operator series

In Eq. 4 of the main article (), we define an exact mapping of the RIXS operator to a set of effective operators. The operators are defined on a complete set of orthonormal states with and as . The expansion coefficients of operator expanded on this complete set are given as

(S5) |

which, as stated in the main text can be derived by multiplying Eq. 4 from the left and right with and .

Absolute point-wise convergence of this series expansion can be proven by realizing that the RIXS intensity is bounded. In a RIXS process, the RIXS operator promotes a given initial state into as

(S6) |

As the RIXS intensity is bounded by some constant (), it follows that , which is a sufficient condition for the series to converge absolutely.

Knowing that these projection methods lead to a convergent series, Eq. 5 to Eq. 8 focus on the question on how to optimize the rate of convergence of this series.

## Appendix B Effective operators for expansion

Equation 8 shows the 6 effective operators used for the expansion of the RIXS operator locally on the core-hole site (). We here provide the 70 operators that arise in the expansion.

There are several ways to create the 70 operators needed for the expansion. One should note that these 70 operators define a (non-unique) basis, and any linear combination of these operators can be used as a basis set to expand the RIXS operator on. This freedom of choice of set of orthonormal operators to project to is much like the difference between a Fourier series and Chebyshev polynomial expansion for functions. In our calculation, we created the non-local operators as products of the complete sets of local operators acting on the core-hole site and its neighboring sites.

As a didactic introduction, we derive here the 70 possible operators following the same arguments as done in the main text for the case. The single-site Fock space containing all states with locally 0, 1 or 2 electrons is spanned by the 4 states {, , , }. The basis set for two sites are given as their outer product, which contains 16 states {, , , , , , , , , , , , }. There are thus linear independent operators possible within this Fock space. We label the 16 basis states from to . Considering the conditions Eq. 5, we define , and thus , , , , and . As the RIXS process is particle number conserving, the coefficients of operators like are zero in the expansion, leaving only 70 operators relevant for the expansion. The above procedure defines a valid and straightforward method to obtain all operators for expansion of any order. However, it does not always give a set of operators which provide the most physical insight into the RIXS process. From here on one can however linear combine them to the more familiar operators, such as the identity, spin, number, and current operators.

Table B shows one possible basis set of operators for the expansion. Subscripts 0 and 1 denote the core-hole site and the bonding linear combination of the nearest-neighbor sites, respectively. Here we define the number conserving local operators of each site as the identity , number , double occupation , and spin operators and . This is equivalent to Eq. 8, but gives a simpler form of the multi-site operators. Creation and annihilation operators with superscript such as is a shorthand for . Note that operators of the form and can be linearly combined as , with the () combination only coupling to the parallel (cross) polarization channel.

1 | 2 | 3 | 4 | 5 | |||||

6 | 7 | 8 | 9 | 10 | |||||

11 | 12 | 13 | 14 | 15 | |||||

16 | 17 | 18 | 19 | 20 | |||||

21 | 22 | 23 | 24 | 25 | |||||

26 | 27 | 28 | 29 | 30 | |||||

31 | 32 | 33 | 34 | 35 | |||||

36 | 37 | 38 | 39 | 40 | |||||

41 | 42 | 43 | 44 | 45 | |||||

46 | 47 | 48 | 49 | 50 | |||||

51 | 52 | 53 | 54 | 55 | |||||

56 | 57 | 58 | 59 | 60 | |||||

61 | 62 | 63 | 64 | 65 | |||||

66 | 67 | 68 | 69 | 70 |

## Appendix C XAS of 1D Hubbard model

Figure S2 shows the calculated XAS spectra for three different doping levels with occupation = 0.83 (hole-doped), 1.00 (half-filled), and 1.17 (electron-doped). The parameters are chosen most relevant to the Cu edge in cuprates, with eV, , and eV. The Coulomb repulsion is defined as , with the coordination number for 1D and for 2D. This on one hand defines the same magnetic energy scale for both cases with approximately , and on the other hand gives a commonly used value for 2D. The inverse core-hole life time is chosen as for the intermediate states.

## Appendix D Overview of exact and approximate RIXS cross sections for 1D Hubbard model

## Appendix E Zeroth and first moment of the exact and approximate RIXS cross sections for 1D Hubbard model

## Appendix F Range dependence of expansion coefficients and exponential rate of convergence

In Fig. 3 of the main text we show the expansion coefficient of and . Here we show the full spectra for . Figs. 3(b) and 3(c) plot the corresponding values of at .