## Abstract

We theoretically investigate the quantum-coherence properties of the cathodoluminescence (CL) emission produced by a temporally modulated electron beam. Specifically, we consider the quantum-optical correlations of CL produced by electrons that are previously shaped by a laser field. Our main prediction is the presence of phase correlations between the emitted CL field and the electron-modulating laser, even though the emission intensity and spectral profile are independent of the electron state. In addition, the coherence of the CL field extends to harmonics of the laser frequency. Since electron beams can be focused to below 1 Å, their ability to transfer optical coherence could enable the ultra-precise excitation, manipulation, and spectrally resolved probing of nanoscale quantum systems.

## INTRODUCTION

Inelastic electron scattering constitutes the basis of various powerful spectroscopy and spectrally selective imaging techniques (*1*). Cathodoluminescence (CL) and electron energy-loss spectroscopy (EELS) are related approaches for harnessing the spectral density of spontaneous interaction processes (*1*, *2*). In CL, the photon emission process involves no external field other than that provided by the electron, so there is no reference to compare the emission phase with. CL emission has spatial and temporal coherences (*3*–*9*) that represent phase correlations in space and time, respectively, with an absolute phase that may differ for each realization. In EELS, the incident and inelastically scattered electron states lose their coherence through the random phase associated with the resonant excitation. A stimulated counterpart for EELS interactions has been established in the form of photon-induced near-field electron microscopy (PINEM) (*10*–*12*) or electron energy-gain spectroscopy (EEGS) (*13*, *14*). In these techniques, the external excitation of a particular mode increases its interaction probability and thus selectively enhances the sensitivity for probing or imaging (*15*, *16*). The coherence of the electron-energy states is evident in a transverse (*17*–*19*) or longitudinal structuring of the electron beam as attosecond pulses (*20*–*23*). Recent work (*24*) has used semi-classical arguments to conclude a dependence of the excitation probability of two-level systems initially in the ground state by an individual electron on its wave function, although such dependence disappears in a full quantum treatment of the system (*25*). For CL, a quantum description of its properties should adhere to our current understanding, in which a point-particle description of the electron is sufficient. By addressing the quantum nature of the electron, one can ask how and to what extent would properties of CL, such as its intensity, coherence characteristics, and radiation pattern (*26*–*28*), be affected by the incident state of the electron.

Here, we directly address these questions in a rigorous theoretical framework, making predictions for the quantum state of radiation produced by phase- and density-modulated electron states. Establishing that an electron beam can coherently stimulate optical excitations, we introduce the notion of “electron-mediated coherence transfer.” Specifically, our results show that the temporal shaping of an electron beam has profound and measurable consequences for inelastic electron-light scattering. In the single-electron limit, we demonstrate that the coherence properties and the quantum-optical correlations strongly depend on the details of the electron state, while the CL intensity and spectral profile remain unaffected by the electronic wave function. In particular, we show that CL emission produced by PINEM-modulated electrons can exhibit mutual coherence with a replica of the PINEM-driving optical field or with its harmonics. We propose interferometric measurements for the extraction of phase information in CL by homodyne tomography of the radiation quantum state. The results are readily applicable not only to generated radiation such as CL but also to polarization created in nonradiative states such as deep-subwavelength excitation of higher-order multipoles (*29*). This concept defines a way of transferring optical polarization carried by electrons, which allow for sub-nanometer precision. Our study could also illuminate the role of single electrons in multi-electron bunches, which, according to several recent proposals (*24*, *30*–*34*), may be used for accessing and manipulating individual quantum systems.

## RESULTS AND DISCUSSION

Figure 1 represents a conceptual system to investigate coherent CL, which can be implemented within an electron microscope. Optical phase information from a PINEM-driving laser field is imprinted on and carried by the electron over a distance *z*, resulting in a coherent CL emission by interaction with an out-coupling sample system. Such electron transfer of optical coherence can be detected by an interferometric setting that targets either the linear or the nonlinear response of free electrons (Fig. 1, B and C, respectively).

We first consider comb-like electron-energy superposition states combined with the radiation vacuum ∣0〉*c _{j}* are the complex probability amplitudes for electron states with energy

*E*, and the index

_{j}*j*runs over electron-energy levels. The coefficients

*c*are normalized as ∑

_{j}*∣*

_{j}*c*∣

_{j}^{2}= 1, and their phases vary with electron propagation in vacuum according to the free-particle dispersion. The quantum-optical and coherence properties of the CL from an electron state with a temporally modulated density depend on the interaction of the electron with an emitter. For simplicity, we consider a single nondegenerate optical band into which CL is emitted. General expressions for the quantum properties of the CL are derived in Materials and Methods, including the radiation continuum, as well as states with definite momentum, energy, and polarization. The photon frequency represents a good quantum number within the interaction bandwidth because it is a single-valued function of the longitudinal momentum. Transverse deflections can be neglected under the nonrecoil approximation (

*34*), which is valid for the examples that we discuss here. The quantum-optical properties of the CL emission can be described by the scattering operator,

*20*,

*36*). For every frequency, the scattering operator allows for annihilation and creation of a photon, marked by the operators

*g*

_{ω}accounts for the electron-photon coupling at the angular frequency ω, where the photon spectral density of the CL is given by ∣

*g*

_{ω}∣

^{2}. Similarly, ∣

*g*

_{ω}∣

^{2}would be the EELS spectral density in the absence of competing loss mechanisms (e.g., bulk plasmons and incoherent emissivity). The phase of

*g*

_{ω}is arbitrary, and

*g*

_{ω}can be chosen as a non-negative real-valued spectral function. However, in examples such as the radiation into normal modes of a fiber, it is convenient to impose a flat spectral phase on the photonic modes and place the spectral degree of freedom as a complex coupling function. A rigorous derivation of Eq. 2 and the conditions for which it applies are given in Materials and Methods (“Evolution operator” section).

The final quantum state *g*_{ω}∣^{2} ≪ 1, for which the final state is

The weak interaction has a small probability amplitude to create a photon with angular frequency ω, represented by state ∣1_{ω}〉, accompanied by a corresponding electron-energy loss. Notably, the CL intensity is unaffected by the specific electron superposition state*t* = 0 is

The proportionality makes it convenient to represent the field with the ladder operator, * _{f}*〉 in Eq. 3, we find

We assume that the PINEM-driven state is a comb separated by the photon energy, *37*) and write *n* is the energy exchange in terms of a harmonic of the fundamental ladder separation,

Since *c _{j}* are the amplitudes of the energy states

*E*, they are proportional to the Fourier coefficients of the electron wave function, that is,

_{j}*c*≡ 〈

_{j}*E*∣ψ〉 ∝ ∫ ψ(

_{j}*t*)

*e*

^{ijω0t}

*dt*, which can be used to simplify the CL field

Incidentally, the temporal electron-probability amplitude, ψ(*t*), can be represented spatially along the propagation axis *v*. Equation 8 is a central result of this paper, representing a general property of CL from a structured electron state (see detailed calculation in the “Mean electric field after interaction” section). We emphasize that only the expectation value of the field follows the electron density, whereas the mean photon number,

Obviously, for a coherent state of light (*38*), the degree of coherence is unity. Furthermore, the above expression remains unchanged for a continuous electron-energy spectrum as expressed in the “Degree of coherence” section.

Let us further explore the predictions of Eq. 10. First, although the DOC is defined as a local property of the CL emission, it depends purely on the properties of the luminescing electron, and not on the material or the geometry of the electron-light coupler. The local electromagnetic energy density in the denominator of Eq. 10 establishes the DOC as an invariant property of a given electron density ∣ψ(*t*)∣^{2}, regardless of the spatial profile of the CL emission. Since spatially varying electron profiles still require temporal density modulations to have a nonzero DOC, we refrain from discussing the complexities that spatial and spatio-temporal effects may add. Second, we have DOC(ω = 0) = 1, since the electron wave function is normalized. Third, the CL field is proportional to *39*)]. We note that the coherent CL field emitted from such a pulsed electron must be pulsed as well, with an equal duration; thus, the FWHM of the optical intensity is shorter by a factor of

Figure 2 presents a sinusoidal modulation of the electron through PINEM and the resulting properties of the DOC for the emitted radiation. In particular, Fig. 2A illustrates a PINEM-driving laser field that imprints an oscillatory phase (orange) on the electron, without an immediate change in the electron density (gray). The electron-phase curvature represents a varying velocity, Δ*v*, which acts as temporal lensing that modulates the electron probability density after propagation. CL (e.g., into an aligned waveguide) should be locked to the phase of the modulated electron and, therefore, to the phase of the driving field. The electron arrival time does not affect the locking to the phase of the PINEM-driving laser, as illustrated by the few electron replicas sketched in Fig. 2A as ellipsoids. Figure 2B shows the DOC for CL from a PINEM-modulated electron, as a function of the propagation distance, *z*, to the CL emitter. For this numerical simulation, we used a standard electron beam acceleration voltage of 200 keV, a wavelength of 800 nm *z* = 0 is given by *e*^{−2∣β∣i cos ω0t}. In particular, Fig. 2B presents the sum ∑* _{j}*(

*c*(

_{j}*z*))*

*c*

_{j+n}(

*z*), which is nonzero only for harmonics of the modulating laser, as predicted by Eq. 7. Here, the initial amplitudes of the electron state are the Bessel functions of the first kind (

*20*,

*36*),

*c*(

_{j}*z*= 0) =

*J*(2∣β∣), corresponding to the Fourier transform of the oscillatory electron phase. The coefficients

_{j}*c*evolve in vacuum as

_{j}*c*(

_{j}*z*) =

*c*(0)

_{j}*e*

^{ikzz}, where

*E*and ℏ is the reduced Planck constant (see the “Degree of coherence for a PINEM-modulated electron” section for an explicit representation). We note that, in some of the literature, β (occasionally marked

_{j}*g*) is referred to as “the coupling coefficient” in PINEM experiments, as it couples adjacent electron-energy levels by stimulated emission and absorption of laser photons. The coupling in this work,

*g*

_{ω}, is the coefficient for exchanging electron-energy loss with spontaneous photon emission, assuming negligible losses. These two quantities are linked, as described previously (

*36*,

*37*): if one populates our conceptual CL emitter with a coherent state having a mean of 〈

*n*〉 photons at frequency ω and use it to drive electrons, the effective PINEM parameter would be

*z*= 6.43 mm (marked with a white dashed line in Fig. 2B), where the CL spectrum is widest. Aside from the

*z*dependence, the peak value of the DOC is fixed and located at exact integer multiples of ω

_{0}. The distribution in the vicinity of the harmonics is determined by the duration of the pre-structured electron, as plotted in Fig. 2C for Gaussian electron pulses. Taking Fig. 2C to two of its extremes, if the electron pulse duration is infinitely long, the DOC is nonzero only at an infinitesimal band for each harmonic. In contrast, if the electron is infinitely compressed to a point-like particle, the DOC of the various harmonics merge, and the CL is fully coherent. For a PINEM-driven modulation, the CL field (proportional to

To offer a quantitative example, we consider CL into a parallel dielectric waveguide (*36*, *40*) produced by a 200-fs-long Gaussian electron pulse that is modulated by PINEM. Our interest in this particular system is explained at a later stage. The waveguide parameters are chosen such that the phase velocity of the optical mode at the frequency of the PINEM-driving laser equals the electron group velocity. For electrons accelerated to 200 keV(corresponding to 69% of the speed of light in vacuum), we find that for electrons passing near the surface of a cylindrical silicon-nitride waveguide, the coupling is most efficient for a diameter of 345.8 nm. Figure 3 describes such an experiment. The CL field, *E _{CL}*, can be collected in a waveguide and mixed with a replica of the driving laser (red),

*E*, using a 50/50 beam splitter. Figure 3B shows the coupling ∣

_{R}*g*

_{ω}∣ alongside the emitted field

*I*

_{1}and

*I*

_{2}(see the “On the calculation of quantum averages of time-dependent field operators” section), is proportional to the CL field

*41*). In addition, since the noise of the reference pulse is suppressed, there is no strict requirement that the reference pulse matches the electron pulse duration. The reference simply needs to overlap with the CL pulse. The next nonzero term of the noise (see the “Intensity and noise in a balanced detection experiment for finite electron and laser pulses” section) can be smaller by orders of magnitude. Consequently, the limiting factor in an experiment would be the buildup of noise originating from the imperfect or imbalanced quantum efficiency of the two detectors and from intensity fluctuations of the reference laser field (

*42*).

Figure 3 (C to E) depicts the amplitudes of the coupling and the field for various propagation lengths of the electron in the near field of the waveguide, alongside the emitted CL field in Fig. 3 (F to H). We note that for a long propagation path along the fiber, the degree of coherence (blue) can be wider than the CL coupling *g*_{ω} (red), leading to a lower signal-to-noise ratio (SNR) and distortions in the temporal shape of the CL field, as shown in Fig. 3H. In such a case, the CL-field emission is optimal for longer electron pulses that match better to the bandwidth of the coupling.

Generally, all of our results are applicable to scenarios with multiple modes, including free-space radiation from surfaces and nano-structures, which are typically characterized by CL. However, choosing CL into a parallel waveguide as the main example in this work has a few instructional, conceptual, and practical merits: (i) An optical fiber can carry a single transverse mode that allows for a straightforward and mode-matched interferometric detection of the field. (ii) The coupling has a simple form, similar to a sinc function. (iii) The coupling can be made strong by macroscopic phase matching between the free electron and the radiation, which could enable the detection of CL from even a single electron. A similar logic was applied to radiation-driven electrons by elongated phase matching near dielectric surfaces (*43*, *44*). (iv) The interaction length provides for a continuous tuning knob to control the coupling bandwidth. (v) Dielectric waveguides can be decoupled from the radiation they host and maintain photon entanglement over many kilometers (*45*). Thus, the CL photons maintain quantum correlations with the electron (*36*), which may be lost if the CL is mediated by lossy polaritonic modes such as plasmons (*6*). (vi) The low dimensionality of a parallel waveguide allowed us to focus this work on temporally modulated electron wave functions. More generally, the consequences of imprinting temporal phases on electrons are extendable to scenarios that spatio-temporally manipulate the electron beam (*17*–*19*). The CL field emitted at optical frequencies from a static transverse structuring of the electron wave functions [e.g., bipolar or vortex beams (*46*, *47*)] has spatial correlations, but the expectation value of its complex amplitude vanishes [see Eq. 8]. (vii) The emission into a dielectric waveguide is instantaneous and allows us to discuss the basic coherence-transfer mechanism. Additional complexity is added if CL is emitted into a cavity with a comparable or longer lifetime than the electron pulse duration as, for instance, a dielectric microresonator. The collected coherent CL from such objects could carry their temporal and phase signatures in addition to that of the electron.

### Conclusions

In conclusion, we have shown that even though the emission rate of CL is independent of the temporal structure of the electron, information embedded on a modulated electron wave function can be retrieved by homodyne mixing with a reference field. Thus, the characterization of the CL emission can reach the shot-noise limit. For a particular emitter, a deviation of the degree of coherence from the calculation shown here can quantify the portion of the incoherent CL emission experimentally. Conceptually, the proposed experimental system is a Mach-Zehnder–like interferometer, mixing the reference light from one arm with the CL-mediated optical coherence carried by the electron in the other arm. We exemplify luminescence into a parallel dielectric waveguide, where CL can be sensitive to even a single electron, and additionally exhibit a narrow linewidth. However, equivalent electron-mediated Mach-Zehnder interferences could be observed using collection systems for CL mounted on scanning and transmission ultrafast electron microscopes, with additional fiber-based or free-space optics. The CL reveals the extreme nonlinear nature of coherence transfer by free electrons, which is capable of emitting a broadband spectrum with intricate correlations, including a degree of coherence for harmonics of the driving field. This concept holds potential for sensitive spectroscopy, as well as for the coherent manipulation of heterostructures and individual quantum systems with femtosecond optical resolution at the atomic scale.

## MATERIALS AND METHODS

### Evolution operator

We consider an electron beam that is well described by a wave packet localized in momentum space around a wave vector **k**_{0} and that the photon energies involved are small compared with the electron relativistic energy *37*)

The current operator *c*_{k}, and the relativistic electron group velocity **v** = ℏ*c*^{2}**k**_{0}/*E*_{0}. *m* and *e* are the electron’s mass and electric charge, and *c* is the speed of light in vacuum. *E*_{0} is the electron’s initial energy and *k*_{0} = ∣**k**_{0}∣. Here, α = *x*, *y*, *z* and *v*_{α} are the projections on axis α of the vector-potential operator at **r** and the electron velocity, respectively. *V* is the volume of the quantization box of the electron. **r**. Therefore, the scattering operator governing the evolution of the system from time *t* = −∞ to *t* = ∞ is given by Di Giulio *et al.* (*34*)

*G*_{j,i}(**r**, **r**′, ω) is the electromagnetic Green tensor at frequency ω, with *i*, *j* = *x*, *y*, *z* satisfying the equation**r** − **r**′) is Dirac Delta distribution function. For electrons traveling parallel to the *z* axis, the Green tensor is rewritten as *G*_{z, i}(**R**_{0}, *z*, **r**′, ω) to split the coordinate **r** into its longitudinal component *z* and the transverse radius-vector **R**_{0}. Im{ϵ(**r**′, ω)} is the imaginary part of the local dielectric constant for an optical frequency ω. *k _{z}* is the electron wave-vector projection on the propagation axis,

*(*

_{i}**R**)∣

^{2}≈ δ(

**R**−

**R**

_{0}), and

^{EELS}(ω) being the real-valued electron energy-loss probability (

*35*). The operator

*48*) and it is not of interest for this study. In addition, the previous operators satisfy the following commutation relations

*g*

_{ω}. In the Results and Discussion section, there is a preferred selection of the optical-mode phase, and hence,

*g*

_{ω}can be complex. In the example of waveguide modes, as in Fig. 3H, the side lobes are a direct result of a sign flip of

*g*

_{ω}within the bandwidth of the coherent emission.

### Initial joint electron-sample state

If we consider the incoming electron being in a superposition of energy states and the sample starting from the ground state (i.e., zero photons), the electron-sample initial state in the interaction picture can be written as* _{kz}* is the coefficient for the momentum component

*k*

_{z}in the electron wave function.

### Mean number of excitations after interaction

By evolving the initial state to the final state

### Mean electric field after interaction

The quantum average (or the expectation value) of the positive-energy electric field operator, defined as (*49*)

We have introduced the incoming electron wave function in the interaction picture *L* denoting the quantization length. Interestingly, this expression can be recast by using the relation

This result is equivalent to Eq. 8, and thus

### Degree of coherence

To quantify the coherence transmitted from the modulated electron to the sample excitations, we define the frequency dependent degree of coherence

This quantity can be calculated from Eq. 16 together with the expression

### Degree of coherence for a PINEM-modulated electron

It is interesting to study the degree of coherence for an electron that is modulated by free propagation after PINEM interaction with a laser field. Here, we assume a PINEM laser frequency ω_{0} and a PINEM coupling parameter β (*1*). In particular, beyond the PINEM interaction region, the electron wave function in the Schrödinger picture reduces to*L* is the quantization length (*50*). Therefore, the Fourier components of the wave function are given by

We have expanded the wave vector up to second order in 𝓁 and we have used the fact that the spatial region of interaction is smaller than the Talbot distance *z _{T}*, namely,

*dz*e

^{i(k−k′)z}=

*L*δ

_{k, k′}. We can now move to the interaction picture by multiplying

^{i[ε0 + v(k−k0)]t}. Additionally, by going back to real space, we obtain the electron wave function in the interaction picture

*vt*by the propagation distance from the PINEM interaction region

*d*. We are now ready to calculate the Fourier transform of the electron density, which reads

From Eq. 17, we see that the integral vanishes unless ω = *n*ω_{0}, where *n* is an integer. With this assumption, we obtain

All the quantities proportional to the Fourier transform in Eq. 18 are calculated for a given harmonic order *n* = ω/ω_{0}.

### On the calculation of quantum averages of time-dependent field operators

According to Glauber’s prescription (*38*), to calculate averages of time-dependent operators to then compute measurable quantities (e.g., light intensities and correlation functions), the operators have to be understood in the Heisenberg picture. For instance, the time-varying light intensity at a given point in space **r** for polarization α has to be calculated as*H* stands for Heisenberg picture. Now, by using an adiabatic switching of the interaction, which provides the connection between the interaction and Heisenberg pictures

From the previous expression, it is clear that the scattering matrix *t* = ∞, which is not the quantity in which one is usually interested. However, a simplification that can be used in this study consists of considering the time-dependent field observables at large times, thus neglecting the few-femtosecond transient period in which the electron is still interacting with the sample and producing nonzero quantum averages of the electric field and light intensity. By doing so, we can extend the final time in the scattering operator to infinity, which leads to

We lastly remark that all the time-dependent quantities in this work imply this assumption.

### Intensity and noise in a balanced detection experiment for finite electron and laser pulses

In this section, we compute the signal and noise in a balanced detection experiment. We assume that a replica of the PINEM-driving laser field is the reference field and that the weak CL emanating from the interaction of the PINEM-modulated electrons with a sample placed downstream is the signal. These two fields are mixed in a symmetric beam splitter, as depicted in Fig. 1. The beam splitter is characterized by reflection and transmission coefficients *R* and *T*, which we assume to be independent of the optical frequency or polarization. Subsequently, the total signal is collected by two ideal detectors, labeled D1 and D2, respectively. We write the reference electric field operator as

Intuitively, the operator *T* and *R* coefficients. The fluence (optical energy per unit area) impinging on each detector is

The difference between the two detectors, which is our signal, is given by

We redefined here the average 〈·〉 symbol such that it includes a continuous mode coherent state ∣{α^{R}}〉, where ^{R}(ω) is the frequency profile of the reference light pulse. One can choose, for example, Fresnel coefficients such as

The noise of such a measurement can be calculated similarly. The square of the noise is defined by the variance as

Since we consider a strong reference laser field, ^{R} dominates Eq. 21. For the choice ∣*R*∣ = ∣*T*∣, the terms for the noise variance scaling as ∣α^{R}∣^{4} and ∣α^{R}∣^{3} vanish. Thus, the next leading order is substantially smaller and comprises combinations in which *t* = 0 is implicitly understood in the CL operators. The contribution of this leading term to the noise is smaller by the ratio of the reference and CL fields, which could be many orders of magnitude smaller than the shot noise of the reference on each detector. Thus, this result implies a theoretical limit to the noise floor, independent of the intrinsic noise of the detectors. We conclude that the corresponding ideal SNR in this system is SNR = S/N.

### Explicit derivation of Eq. 8

One can separate the initial electron-photon state in Eq. 1 as the product of the electron part ∣ψ* _{e}*〉 and the vacuum of the radiation, ∣ψ

*〉 = ∣ψ*

_{in}*〉 ⊗ ∣0〉, where ∣ψ*

_{e}*〉 = ∑*

_{e}*∣*

_{j}c_{j}*E*〉. The ladder coefficients can be derived through simple algebra

_{j}We can now write the sum in Eq. 7 as

Here, we assume that the energy states are identical, aside from their energy difference, and that ∣ψ* _{e}*〉 is fully spanned by the discrete and complete set of ∣

*E*〉 states, that is

_{j}In the sample region, we can approximate the electron dispersion as linear, and thus, the wave function can be written as a function of time, ψ(*t*), relying on the relation *z* − *vt* = *constant*. Using this temporal electron wave function ψ(*t*), one finds

### Generalization for strong electron-photon coupling

Starting from the scattering operator in Eq. 2, we are interested in the expectation value for a specific frequency Ω, that is,

Since

The expectation value for

The commutation relation above allows us to write

### Higher-order correlations

The moment of order *N* for the quantum correlations of the CL is

The equation marked with (*) uses the commutation in Eq. 23. We can now substitute

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**

**Funding:**This work was funded in part by the Deutsche Forschungsgemeinschaft (DFG) in the Collaborative Research Center “Atomic Scale Control of Energy Conversion” (DFG-SFB 1073, project A05) and in the Priority Program “Quantum Dynamics in Tailored Intense Fields” (DFG-SPP 1840), ERC (Advanced Grant 789104-eNANO), the Spanish MINECO (MAT2017-88492-R and SEV2015-0522), the European Commission (project no. 101017720-EBEAM), and the Humboldt Foundation. O.K. acknowledges the Max Planck Society for funding from the Manfred Eigen Fellowship for postdoctoral fellows from abroad. V.D.G. acknowledges support from the EU (Marie Skłdowska-Curie Grant 713729).

**Author contributions:**O.K., V.D.G., F.J.G.d.A., and C.R. conceived the concept. O.K. constructed the paper and prepared the figures. O.K. and V.D.G. wrote the paper with contributions from F.J.G.d.A. and C.R.

**Competing interests:**The authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

- Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).