Magnon interactions in a moderately correlated Mott insulator

Quantum fluctuations in low-dimensional systems and near quantum phase transitions have significant influences on material properties. Yet, it is difficult to experimentally gauge the strength and importance of quantum fluctuations. Here we provide a resonant inelastic x-ray scattering study of magnon excitations in Mott insulating cuprates. From the thin film of SrCuO2, single- and bi-magnon dispersions are derived. Using an effective Heisenberg Hamiltonian generated from the Hubbard model, we show that the single-magnon dispersion is only described satisfactorily when including significant quantum corrections stemming from magnon-magnon interactions. Comparative results on La2CuO4 indicate that quantum fluctuations are much stronger in SrCuO2 suggesting closer proximity to a magnetic quantum critical point. Monte Carlo calculations reveal that other magnetic orders may compete with the antiferromagnetic Néel order as the ground state. Our results indicate that SrCuO2—due to strong quantum fluctuations—is a unique starting point for the exploration of novel magnetic ground states.

One of the most studied electronic models is that of a two-dimensional square lattice [1].
The physics of these systems is well captured by the Hubbard model, from which an effective Heisenberg model can be generated to describe the AF ground state.In this strong-coupling limit, the spin exchange interactions are realized through virtual hopping processes.Upon down-tuning the interaction strength, the AF Mott state remains a theoretical ground state.However, in this limit, small perturbations (for example doping) can trigger new magnetic or metallic ground states [2,5].Pushing the Mott state into this soft regime is therefore of great interest.When potential and kinetic energy scales are comparable, quantum fluctuations enter the problem.Many experimental and theoretical studies have addressed this intermediate region of U/t.Upon doping, high-temperature superconductivity has been found [5] and spin-liquid states are predicted in certain models [6].
In spin-wave theory [7], quantum fluctuations describe the occupations of the bosonic modes as perturbations to the Néel state [8].The ground state can be viewed as a Néel state with finite boson density.The corresponding elementary excitation should therefore be considered as a magnon renormalized by its higher-order expansions, i.e., the magnon-magnon interactions.On an experimental level, it has however been difficult to measure or gauge the strength of quantum fluctuations.As mentioned above, magnon excitations [9] of the Néel state should be influenced by quantum fluctuations.The magnon dispersion is described by ℏω = Z c (k)ϵ k , where ϵ k is the "bare" magnon dispersion set by potential and kinetic energy scales, and Z c (k) = Z 0 c (1 + f k ) is the renormalization factor stemming from quantum fluctuations with f k being a momentum dependent function.In the strong coupling limit (U/t → ∞) quantum fluctuations are suppressed implying f k → 0 and Z c (k) ≈ 1.18 is essentially momentum independent [10,11].As quantum fluctuations grow stronger with gradually moderate values of U/t, the renormalization factor Z c (k) increases and ac-quires momentum dependence through a non-negligible f k .This limit governed by quantum fluctuations is interesting as it may provide physics beyond the antiferromagnetic Néel state.
Conceptually, this moderate U/t limit is complicated due to a multitude of comparable magnetic exchange interactions.Nearest and next-nearest neighbour exchange interactions are given by J 1 = 4t 2 /U and J 2 = 4t 4 /U 3 , in the projection of the Hubbard onto the Heisenberg model.In the moderate or weak interaction limit, higherorder exchange interaction terms are gaining prominence.The ring-exchange interaction J □ = 80t 4 /U 3 becomes a significant fraction of J 1 (J □ /J 1 = 20t 2 /U 2 ) and manifests by a magnon zone-boundary dispersion [12][13][14][15][16][17].In this limit, higher-order hopping integral t ′ , can introduce new magnetic interaction term J ′ = 4t ′2 /U that further adds to enhance the zone boundary (ZB) dispersion.Within the Hubbard-Heisenberg model, the zone boundary dispersion, quantum fluctuations, and Z c correlate in the U/t and t ′ /t parameter space.In fact, the renormalization factor Z c gains its momentum dependence from the higher-order exchange interactions.Enhanced quantum fluctuations may thus introduce new magnetic ground states and with that exotic magnonic quasiparticles [18,19].It is thus interesting to study materials with significant higher-order exchange couplings.In the cuprates, ACuO 2 with A = Sr, Ca has been studied with electron spectroscopy and resonant inelastic xray scattering (RIXS) [17,18,20] due to its large ring exchange interaction.Yet, no experiments have demonstrated the importance of magnon-magnon interactions and quantum fluctuations through direct measurements of a momentum dependent magnon renormalization factor.
Here we provide a RIXS study of SrCuO 2 (SCO) realized in thin-film format, which demonstrates a Mott-insulating nature by electron spectroscopy measurements [3,4].Analysis of the RIXS spectra led us to derive the single-and bi-magnon dispersions.Starting from an effective Heisenberg representation of the Hubbard model, we show that the observed single-magnon dispersion is inconsistent with a constant Z c for reasonable values of kinetic energy scales.We thus conclude that, in SrCuO 2 , quantum fluctuations are significantly influencing the magnon dispersion.Further, our analysis shows that the observed magnon dispersion is well described when introducing significant momentum dependence to Z c by evaluating the magnon-magnon interactions.
This finding is further supported by comparing cuprate compounds with different correlation strengths.Our results thus provide a gauge for quantum fluctuations which are getting increasingly important as U/t is reduced.Possible exotic magnetic ground states, emerging from quantum fluctuations are here explored by classical Monte Carlo calculations.

Results
The crystal field environment around the copper site in SrCuO 2 is shown schematically in Fig. 1a.In contrast to, for example, La 2 CuO 4 (LCO), no apical oxygen is present in SrCuO 2 .Examples of Cu L-edge RIXS spectra, covering magnetic and dd excitations, are shown in Fig. 1b.The absence of the apical oxygen in SrCuO 2 pushes the d z 2 excitation well below the t 2g excitationsas previously established in CaCuO 2 [17,22].The two excitations at 1.43 and 1.73 eV are assigned to the d xy and degenerated d xz /d yz states, respectively.The origin of the additional peak at ∼2.06 eV, which has also been observed in CaCuO 2 [17,22], has been attributed to the incoherent component of the d xz /d yz orbital excitations due to the coupling to magnetic excitations [23].
Magnetic excitations have been recorded systematically along the (h, 0), (h, h), and zone boundary (azimuthal ϕ rotation with a constant in-plane momentum amplitude Q // ) directions with both linear vertical (σ) and linear horizontal (π) incident light polarisations.A single-magnon excitation manifests clearly in the π channel (see Fig. 1).When switching to σ polarisation, the single-magnon is suppressed as expected [24], and an excitation at higher energy appears.We interpret this as a magnetic continuum that in SrCuO 2 (and CaCuO 2 ) has a structure-sometimes referred to as a bi-magnon excitation [25][26][27].In what follows, we extract the singlemagnon and bi-magnon dispersions along the high symmetry directions.
We analyze the low-energy part of the RIXS spectra by considering four components that include elastic scattering (grey shaded area-enhanced in the σ channel [24]), single-and bi-magnon (orange shaded area) excitations, and a smoothly varying background (grey dashed line).Elastic scattering is mimicked by a Gaussian function centred at zero energy loss.The energy width is slightly larger than the instrumental resolution due to unresolved phonon modes [28,29].The single-and bi-magnon excitations are described respectively by a damped harmonic oscillator convoluted with the instrumental resolution and a Gaussian function.Background is modeled by a second-order polynomial.In grazing-exit geometry, RIXS cross section from the magnon (bi-magnon) is generally enhanced when using the π (σ) polarized incident lights-as shown in Fig. 1c,d and Supplementary Fig. 1.We fit globally across the two light polarisations to extract the two magnetic contributions.The resulting single-and bi-magnon dispersions are plotted on top of the inelastic RIXS spectral weight (π polarisation) in Fig. 2a-c.Consistent with previous reports on CaCuO 2 [17,18], a large zone boundary dispersion ) of the single magnon excitation is observed with an essentially non-dispersive section along the (h, h) direction.Away from the Brilliouin zone centre, the bi-magnon dispersion ℏω bm roughly mimics the single magnon dispersion ℏω sm .At the zone boundary position ( 14 , 1 4 ), ω bm /ω sm ≈ 2. This ratio however varies significantly along the high symmetry directions.
The fitting of the single-and bi-magnon excitations also provides information about spectral weight and quasiparticle lifetime.For most of the Brillouin zone, the energy width of the single-magnon is resolution-limited.However, around (0.5, 0) spectral weight suppression and shorter single-magnon lifetimes are observed consistently with what has previously been reported in La 2 CuO 4 and CaCuO 2 [14,17,18].

Discussion
The single-magnon dispersion of SrCuO 2 features two peculiar characteristics.A steep zone boundary dispersion is followed by a non-dispersive section along the (h, h) direction.Magnon excitations of layered copperoxides have been discussed via a Heisenberg Hamiltonian derived from the Hubbard model [2].In the simplest form, the nearest-neighbour exchange interaction J 1 = 4t 2 /U is described through the Coulomb interaction U and nearest-neighbour hopping integral t.The magnon dispersion is, in this limit, isotropic-given by ℏω  [12], however, revealed a zone boundary dispersion indicating the importance of higher-order exchange interaction terms.To account for this zone boundary term, a ring exchange interaction J □ ∼ t 4 /U 3 was included to satisfactorily describe the observed magnon dispersion [12,14].Later more detailed studies [21,30] included higher-order hop-ping terms, i.e., next, and next-next nearest-neighbour hopping intergrals t ′ and t ′′ .This extended model, yields a magnon dispersion ℏω = Z c (U, t, t ′ , t ′′ )ϵ k (U, t, t ′ , t ′′ ), where Z c is a momentum dependent quantum renormalization factor and ϵ k (U, t, t ′ , t ′′ ) = A 2 k − B 2 k is the bare magnon dispersion with A k and B k determined by U , t, t ′ , and t ′′ , as described in refs.21 and 30.In the U/t → ∞ limit, the magnon-magnon renormalization factor Z c is momentum independent, and ϵ k (U, t, t ′ , t ′′ ) has an analytic expression [16,31].This expression has been used to fit the magnon dispersion of La 2 CuO 4 , with realistic values of U , t, t ′ , and t ′′ .In particular, ratios of t ′ /t ∼ −0.4 and t ′′ /t ∼ 0.2 are found consistent with density functional theory calculations [32] and angle-resolved photoemission spectroscopy experiments [16,33].
For SrCuO 2 , however, the constant Z c solution does not provide a satisfactory description of the observed single-magnon dispersion (see blue dashed lines in Fig. 3g-l).As La 2 CuO 4 and SrCuO 2 share similar square lattice structures, similar values of t ′ /t are expected.However, an unbiased fit yields unphysical values for the hopping parameters and physical sensible values provide poor fits.We are thus led to reject the initial ansatz in a numerical self-consistent fashion.For La 2 CuO 4 , this methodology confirms that Z c is roughly constant (see Fig. 3a-c and Table I) with marginal changes to U, t, t ′ , t ′′ (compared to the constant Z c model [16]).However, for SrCuO 2 , an entirely new solution emerges.Values of t ′ /t and t ′′ /t comparable to La 2 CuO 4 and a smaller U/t now describes the magnon dispersion-see solid lines in Fig. 3 and Table I.We stress that this new solution describes the observed dispersion using fewer fitting parameters as Z c is now given by U, t, t ′ , t ′′ .The moderate value of U/t implies a magnon-magnon renormalization factor that is strongly momentum dependent (see Fig. 3g-i Our results thus indicate that quantum fluctuations have a significant impact on the magnon dispersion in SrCuO 2 . In Fig. 4, we plot the constant χ 2 (goodness-of-fit) contour lines that encircle solutions that are within 10% of the minimum of χ 2 .With the currently available data, a rather broad set of parameters describe the magnon dispersion of La 2 CuO 4 .For SrCuO 2 , however, χ 2 has a unique well-defined minimum confined to a narrow region of the parameter space.In Table I, the fitting values to describe the magnon dispersions for La 2 CuO 4 and SrCuO 2 are listed.For SrCuO 2 , these values represent a minimum in the χ 2 function.For La 2 CuO 4 , we further constrain the solutions by fixing t ′ /t = −0.4.In Fig. 4, the modeled zone boundary dispersion is plotted as a function of J □ /J 1 and J ′ /J 1 .Within the same parameter space, the Brillouin zone average Zc is shown.Generally, the model displays a correlation between magnon ZB dispersion and Zc .The Heisenberg-Hamiltonian projection from Hubbard model is breaking down in the limit where Zc ≫ Z 0 c -that is when (J □ + J ′ )/J 1 is large.Compared to La 2 CuO 4 , SrCuO 2 displays a larger ring exchange coupling J □ .In fact, the fitting parameters obtained for SrCuO 2 are close to the limit where the Heisenberg representation of the Hubbard model breaks down.This limit is characterized by a complete suppression of the staggered magnetization and the imaginary solution of the magnon dispersion (Fig. 4).It has been reported that the absence of apical oxygen in cuprates leads to a decrease in the electronic correlation strength [34,35], which agrees with the observations here.To corroborate our findings, we analysed the magnon dispersion of CaCuO 2 (CCO), which is another infinitelayer cuprate compound with a large ring-exchange coupling [17].As observed in ref. [17], the one-band Hubbard model with underestimated quantum renormalization fails to describe the magnon dispersion and yields an unrealistically small U (U/t = 4.9).As shown in Fig. 4 and Supplementary Figs. 2 and 3, the magnon dispersions in both compounds are only well described when in- and SCO/GSO (blue).The solid lines are corresponding fits using a Hubbard model including higher-order terms (see text).The parameters extracted from the fits are listed in Table I.Blue dashed lines mark the magnon dispersion obtained assuming a constant Zc = 1.219 [16,21], with U = 2.15 eV, U/t = 6.25, t ′ /t = −0.4,and t ′′ /t ′ = −0.5.Error bars indicate one standard deviation.The in-plane momentum amplitude Q // takes 0.461, 0.444, and 0.463 for LCO/STO, LCO/LSAO, and SCO/GSO, respectively.Data on LCO/STO and LCO/LSAO are taken from ref. [16].Top panels (a-c, g-i) display the momentum dependence of the quantum fluctuation factor Zc obtained from fitting with the Hubbard model.
cluding substantial quantum corrections generated from magnon-magnon interactions.We thus demonstrate that SrCuO 2 -a moderately correlated Mott insulator-hosts strong quantum fluctuations that can potentially stabilize ground states beyond the AF ordered Néel state.Enhancing further J □ /J 1 or J ′ /J 1 would be of great interest to explore new quantum matter ground states.
To gain insight into the possible magnetic ground states when the Néel order breaks down by the increase of t and t ′ , we perform Monte Carlo calculations using a Heisenberg Hamiltonian including the first-, second-, and third-nearest-neighbour, as well as a four-spin ring exchange coupling (see Methods).To compare with the experimental results on SrCuO 2 , we fixed U = 2.66 eV and t ′′ /t ′ = −0.5.Incommensurate magnetic orders characterized by a quartet of magnetic Bragg peaks around (0.5, 0.5), i.e., with wave vectors Q M = (0.5 ± δ, 0.5) and (0.5, 0.5 ± δ), or Q M = (0.5 ± δ/ √ 2, 0.5 ± δ/ √ 2) are found in the parameter space between the antiferromagnetic Néel and columnar orders (see Supplementary Fig. 4 for examples of the calculated spin structure factor).We plot in Fig. 4d the distance δ between Q M and the Néel wave vector (0.5, 0.5) as a function of t and t ′ /t.While SrCuO 2 is in the Néel AF ordered state, it is located not far from incommensurate magnetic ordered phases which can be reached by increasing t.Further increase of t ′ /t enhances the next-nearest neighbour coupling and stabilizes the columnar antiferromagnetic order.This tuning can be potentially realized by strain application with different substrates [16,36].Explorations on possible strain or pressure induced quantum critical behaviour would provide more insights into the nature of the magnetic ground states.Note that recent theoretical and numerical works on the two-dimensional Hubbard model have shown that under small doping, the Néel order becomes unstable and replaced by other magnetic orders [37][38][39][40][41].A spin-charge stripe state with an incommensurate ordering wave vector has also been experimentally established in underdoped cuprates [42][43][44].We point out that the classical Monte Carlo simulations do not capture the quantum nature of the problem, but they show the parameter space, where the Néel state is expected to break down.When antiferromagnetic order is suppressed by the strong quantum fluctuationsstemming from magnon-magnon interactions-near the phase boundary, superconductivity could be potentially enhanced.It would therefore be of great interest to  12 , 0), renormalized staggered magnetization as a function of the exchange interactions J □ /J1 and J ′ /J1.(c) Average renormalization factor Zc across the Brillouin zone for the same parameter space as in (a,b).Green, purple, yellow, and blue dashed circles indicate constant χ 2 contour lines with solutions that are within 10% of the minimum of χ 2 for LCO/STO, LCO/LSAO, CCO, and SCO, respectively.Grey, green, and purple filled symbols denote parameters determined from the Hubbard model with constant Zc on respectively bulk LCO [12], LCO/STO, and LCO/LSAO [16].Magnon dispersion of CCO is adapted from ref. [17].Details of the analysis are described in the Supplementary Information.Empty areas in (a,c) and (b) indicate where the magnon dispersion becomes imaginary and the magnetization negative, respectively.(d) Magnetic ground state structure as a function of t and t ′ /t obtained from Monte Carlo calculations with U = 2.66 eV and t ′′ /t ′ = −0.5 fixed.White pentagram marks the position of SCO.
study how the different magnetic ground states influence superconductivity upon doping.
Future theoretical studies including extended dynamical mean-field theory calculations on the spin excitations [45][46][47][48][49]-beyond the scope of the current work-could offer more insights into the relationship between the quantum fluctuations and the degenerate ground states near the critical region.On the experimental front, such quantum effects could also be addressed by comparative RIXS measurements at a higher temperature (T ∼ 0.1J), which call for future exploration.

RIXS experiments
Cu L 3 -edge RIXS experiments were carried out at the ADRESS beamline [53,54], of the Swiss Light Source (SLS) synchrotron at the Paul Scherrer Institut.All data were collected at base temperature (∼20 K) under ultrahigh vacuum (UHV) conditions, 10 −9 mbar or better.RIXS spectra were acquired in grazing-exit geometry with both linear horizontal (π) and linear vertical (σ) incident light polarisations with a fixed scattering angle 2θ = 130 • .The two-dimensional nature of the system ensures that the out-of-plane dependence of the magnon dispersion is negligible, as confirmed in a recent RIXS study on CaCuO 2 [18].The energy resolution, estimated by the full-width-at-half-maximum of the elastic scattering from an amorphous carbon sample, is 118.5 meV at Cu L 3 edge (∼931.5 eV).Momentum transfer is expressed in reciprocal lattice units (r.l.u.) using a pseudo-tetragonal unit cell with a = b = 3.97 Å and c = 3.4 Å for SrCuO 2 , a = b = 3.8 Å and c = 13.2Å for La 2 CuO 4 .RIXS intensities are normalised to the weight of dd excitations [55].
The staggered magnetization is calculated by , where N is the number of spin on the square lattice.Our code is available upon request.

Monte Carlo simulations.
Monte Carlo simulations were carried out using the Heisenberg Hamiltonian projected from the Hubbard model, taking into consideration the leading contribu- where For each set of input values, (U, t, t ′ , t ′′ ), the simulation ran for 10 7 Monte Carlo steps with random starting configurations.The Monte Carlo calculation code is available upon request.

FIG. 1 .
FIG. 1. RIXS on SrCuO2 with different incident light polarisations.(a) Crystal structure of SrCuO2 and schematic illustration of the RIXS scattering geometry.Incident light, either linear horizontal (π) or vertical (σ), is directed to the film with variable angle θ.(b) RIXS spectra at Q = (0.34, 0) measured with π (blue line) and σ (purple line) incident x-rays.Spectra are vertically shifted for clarity.Dotted lines mark the peak positions of the dd excitations, determined from fitting with Gaussian components denoted by the shaded areas.(c,d) Zooms of the low-energy part (within the grey dashed boxes) of the spectra in (b).Solid lines are the sum of a four-component fit.Each component is indicated by dashed lines and shaded areas-see text and Supplementary Information for more details.Vertical dashed lines indicate the zero energy loss.

hFIG. 2 .
FIG. 2.Magnon and bi-magnon dispersions and spectral weights observed in SrCuO2.(a-c) RIXS intensity as a function of momentum and energy loss measured with π polarized incident light, along three different directions in the reciprocal space as indicated by the solid color lines in the inset.Elastic and background scattering has been subtracted.Open circles and diamond points indicate respectively the magnon and bi-magnon pole position (see text for detailed description of the analysis).In (b) the in-plane momentum amplitude is Q // = 0.463.(d,e) Single magnon inverse lifetime Γ ∼ ℏ/τ along the (h, 0) and (h, h) directions.Horizontal dashed line indicates the applied energy resolution.(f-i) Normalized single-magnon (f,g) and bi-magnon (h,i) spectral weight along the (h, 0) and (h, h) high symmetry directions.Error bars are determined from the fitting uncertainty.

hhFIG. 3 .
FIG.3.Magnon dispersion and momentum dependence of quantum renormalization factor Zc. Bottom panels (d-f, j-l) display the magnon dispersion along the indicated momentum trajectories for LCO/STO (green), LCO/LSAO (pink), and SCO/GSO (blue).The solid lines are corresponding fits using a Hubbard model including higher-order terms (see text).The parameters extracted from the fits are listed in TableI.Blue dashed lines mark the magnon dispersion obtained assuming a constant Zc = 1.219[16,21], with U = 2.15 eV, U/t = 6.25, t ′ /t = −0.4,and t ′′ /t ′ = −0.5.Error bars indicate one standard deviation.The in-plane momentum amplitude Q // takes 0.461, 0.444, and 0.463 for LCO/STO, LCO/LSAO, and SCO/GSO, respectively.Data on LCO/STO and LCO/LSAO are taken from ref.[16].Top panels (a-c, g-i) display the momentum dependence of the quantum fluctuation factor Zc obtained from fitting with the Hubbard model.

2 FIG. 4 .
FIG. 4.Quantum renormalization effect within the Hubbard model.(a,b) Zone boundary dispersion ratio EZB/E( 1 2 , 0), renormalized staggered magnetization as a function of the exchange interactions J □ /J1 and J ′ /J1.(c) Average renormalization factor Zc across the Brillouin zone for the same parameter space as in (a,b).Green, purple, yellow, and blue dashed circles indicate constant χ 2 contour lines with solutions that are within 10% of the minimum of χ 2 for LCO/STO, LCO/LSAO, CCO, and SCO, respectively.Grey, green, and purple filled symbols denote parameters determined from the Hubbard model with constant Zc on respectively bulk LCO[12], LCO/STO, and LCO/LSAO[16].Magnon dispersion of CCO is adapted from ref.[17].Details of the analysis are described in the Supplementary Information.Empty areas in (a,c) and (b) indicate where the magnon dispersion becomes imaginary and the magnetization negative, respectively.(d) Magnetic ground state structure as a function of t and t ′ /t obtained from Monte Carlo calculations with U = 2.66 eV and t ′′ /t ′ = −0.5 fixed.White pentagram marks the position of SCO.

2 .
The first three sums count respectively nearest, next-nearest, and next-next-nearest neighbour using indices ⟨i, j⟩, ⟨i, i ′ ⟩ and ⟨i, i ′′ ⟩.The last sum counts around the squares following the clockwise direction.Simulations were performed for a temperature T = 0.1 K on a sheet of 50 × 50 unit cells in the (a, b)-plane, i.e., 2500 magnetic sites with classical spin S = 1 2 .
. Early neutron scattering experiments on La 2 CuO 4