Abstract
We present a fully quantum model for the excitation of a bound electron based on the “free-electron bound-electron resonant interaction” (FEBERI) scheme. The bound electron is modeled as a quantum two-level system (TLS) at any initial quantum (qubit) state, and the free electron is presented as a pre-shaped quantum electron wavepacket (QEW). In the case that the QEW is short or modulated at optical frequency, the TLS quantum state may be coherently controlled with multiple modulation-correlated QEWs. For this case, we derive the transition probability of the TLS due to interaction with a multi-particle beam based on an analytical approximate solution of the Schrodinger equation that amounts to using Born’s probabilistic interpretation of the quantum electron wavefunction. We verify the credibility of the analytical model at its validity ranges using a fully quantum density matrix computation procedure. It is shown that the transition probability can grow quadratically with the number of correlated QEWs and exhibit Rabi oscillation. The study indicates a possibility of engineering the quantum state of a TLS by utilizing a beam of shaped QEWs.
Introduction
Recent technological advances enable the shaping of free electron wavefunctions in the transverse [1–3] and longitudinal [4–6] dimensions in an electron microscope setup. In the longitudinal propagation direction of the electron (or in time dimension), the energy and density expectation values of the electron wavefunction can be modulated by a coherent laser beam at optical frequencies [7–9], utilizing the scheme of photo-induced near-field electron microscopy (PINEM) [10]. The modulation is made possible by a multiphoton emission/absorption process in the near field of a nanostructure [7], a foil [11], or a laser-beat (pondermotive potential) [12, 13]. After the PINEM interaction, the energy spectrum of a single quantum electron wavepacket (QEW) exhibits discrete energy sidebands spaced apart by the laser photon energy . It was also shown that due to the nonlinear energy dispersion of electrons in free space drift, the discrete energy modulation of the QEW is converted into tight periodic density modulation (bunching) at attosecond short levels [14, 15]. Furthermore, it was shown that single QEWs can be shaped and compressed to femtosecond pulse duration by means of a chirping and streaking process with coherent THz or Infra-Red (IR) beams [16, 17]. The reality of the finite-sized shape and the density modulation of the QEW in interaction with light was asserted and analyzed in the context of stimulated Smith-Purcell radiation [18, 19]. It has been confirmed experimentally by the acceleration/deceleration interaction of the pre-shaped/pre-modulated QEW with another coherent laser beam synchronous with the modulating laser [20, 21]. The reality of the QEW sculpting and modulation features in stimulated radiative interaction and superradiance has been asserted also in the case of multiple modulation-correlated electron wavepackets [22, 23], as an extension of the classical case of a pre-bunched particle beam [18].
Shaping and transfer of coherence and quantum properties from light to free electrons wavefunction have received recently significant research attention within an emerging new research field of “quantum electron optics”. It has been shown that the coherence and incoherence features of light, and even the quantum statistical states of light, can be transferred to the free-electron quantum electron wavefunction by means of the PINEM process [24, 25]. This paves the way for new applications of these emerging technologies for coherent control of quantum systems of light and matter using pre-shaped and optically modulated QEWs.
A simple example of coherent control of quantum states by electrons is the effect of free-electron bound-electron resonant interaction (FEBERI) proposed in Ref. [26]. In this process, a pre-shaped or pre-modulated beam of QEWs interacts with a bound electron modeled by a quantum two-level system (TLS). Such a TLS model is valid because, in general, in the linear response regime, the transition amplitude responds linearly to all possible transitions of quantum levels in matter with frequency-dependent amplitude. Therefore, targeting a single two-level transition does not reduce the generality of the model, even if the targets have multiple levels. The QEWs induce excitations of the TLS when passing in the vicinity of the TLS target. It has been suggested that a beam of probability-density modulation-correlated QEWs can interact resonantly with the TLS when the harmonic of its optical frequency modulation matches the TLS quantum energy level transitions , where is the periodic temporal modulation frequency of the QEW density distribution and is the quantum energy gap of the TLS. This assertion has raised a debate [27–29], but also a stream of numerous recently published papers relating to different aspects of this effect and its potential applications in electron microscopy, diagnostics, and coherent control of qubits [30–33].
Using a fully quantum-mechanical analysis (both free and bound electrons quantized) of the FEBERI interaction with a single arbitrarily shaped QEW [34, 35], we showed that the FEBERI effect can be applied for coherent control and interrogation of the qubit state of a target TLS. However, because the FEBERI effect is practically very weak for a single QEW and single TLS, it is necessary to consider the interaction of the TLS with multiple modulation-correlated QEWs. Such a beam of density modulation-correlated QEWs can be produced if all electrons are energy modulated in the PINEM process by the same coherent laser beam and then develop the same phase-correlated density bunching after drifting to the FEBERI interaction point. In this case, as proposed in the early paper of Gover and Yariv [32], based on a semi-classical model, a coherent build-up of the quantum state of the TLS may be possible, including full Rabi oscillation between the TLS quantum levels. Such a Rabi oscillation process that initially scales quadraticlly with the number of interacting electrons N, would be the analog of spontaneous superradiance (in the sense of Dicke [36]) by a density modulated electron beam, which is characterized by an emission rate proportional to where is the number of interacting electrons. This process is well established for bunched classical particle electron beams [18, 37], but was also shown to take place when the electrons are density-modulated on the level of the quantum wavefunction (in expectation value) and modulation-correlated [22, 38] Likewise, the semiclassical analysis of the FEBERI effect results in quadratic growth of the transition from the ground state to the upper state of the TLS at resonance condition [39].
In this paper, we study the dynamics of a TLS excited by multiple finite size or modulation-correlated QEWs under the assumption that the interaction period is smaller than the relaxation time and the decoherence time of the TLS. The analytical model is based on a semiclassical approximation of the Schrodinger equation for the entangled free and bound electrons. We present an approximation in which the free-electron quantum recoil is neglected. This amounts to using a probabilistic model for the free electron location, relying on Born’s interpretation of the quantum electron wavefunction. In parallel, we present for the first time a fully quantum model for single and multiple electrons FEBERI interactions, and use it to compute numerically the transition probability dynamics of the TLS. This is employed to evaluate the applicability range of the approximate analytical model, and specifically, to rigorously affirm the N2 scaling of the FEBERI effect predicted in the semiclassical theory [26] which has been under debate [27–29].
The paper is organized as follows: In the first section we present the model and the theoretical framework of the FEBERI setup. In the subsequent two sections we apply the probabilistic model to interpret the excitation of a TLS with a single near point-particle (short size) QEW and with a beam of such QEWs. In the next two sections we present the excitation of a TLS with a single density modulated QEW and then, with a beam of such modulation correlated QEWs. All approximate analytical expressions are compared to the results of the quantum density matrix computations. Conclusions are presented in the final section.
System Model and Theoretical Framework
The setup of our system model is shown in Figure 1, in which a thin free-electron QEW propagates in proximity to a Hydrogen-like atom that is modeled as a TLS. The interaction of the free electron and bound electron is considered to be the Coulomb interaction. We denote the joint wavefunction of the free and bound electrons by . Then the dynamics of the considered system are governed by the Schrödinger equation,where is the unperturbed Hamiltonian of the free electron and bound electron. represents the interaction Hamiltonian. In order to apply the analysis also to relativistic electrons, we use a “relativistic” Hamiltonian for the free electron of energy which we have derived by second-order iterative approximation of the Klein-Gordon equation, neglecting the spin effect [41,48]:
FIGURE 1
This corresponds to second-order expansion of the relativistic energy dispersion of a free electron / , where and . This “relativistic Hamiltonian” has been derived recently also directly from the Dirac equation [40] without the quadratic term.
For the simplified model where the spin is neglected, we assume that the free and bound electrons do not overlap spatially. Therefore, there are no exchange energy or spin-orbit interaction effects, and we can avoid the intricate second quantization of the many-body interaction theory [41]. For the Coulomb interaction in the near field and neglecting retardation [42], the interaction Hamiltonian is
Here we used Feynman’s expression for the Coulomb potential [43] in order to keep the analysis valid in the relativistic regime. A more accurate form would be to use the Darwin potential for relativistic Coulomb interaction between moving charged particles [44, 45]. However, for the parameters of the cases delineated here, the corrections due to this model are negligible.
Having the bound electron modeled as a TLS, the eigenfunction solutions of the Hamiltonian satisfywhere . Then the general wavefunction of the bound electron can be represented as , where . The wavefunction solution of the free electron in zero-order is taken to be a general wavepacket
The shape of the wavepacket is determined by the coefficients in momentum dimension - . Here we consider two cases. The first case is a finite size QEW of Gaussian shape, arriving at the FEBERI interaction point after a free drift length from its waist position. In 1-D momentum space it is presented as:
The second case is a PINEM phase-modulated Gaussian wavepacket arriving at the FEBERI interaction point after a free drift length from the PINEM modulation point:
Such a wavepacket represents in spatiotemporal space a density modulated finite-size QEW [14]. Here is the PINEM coupling strength parameter, and are the frequency and reference phase of the laser-beam-induced field at the PINEM interaction point, and is the QEW momentum spread. For simplicity, we assume that the QEW reaches the interaction point at a time at its longitudinal waist so that the axial coordinate spread of the QEW is (the expansion of the QEW during the short interaction time is negligible [19]).
During the interaction process, the expansion coefficients of the QEW get entangled with the coefficients of the bound electron . Then the generally combined wavefunction of the free electron and bound electron during the interaction can be represented in terms of the eigenfunctions:
After substituting this expression into Eq. () and canceling out the no-interaction terms, we are left with
By multiplying by and integrating over space, we reach an integrodifferential equation that needs to be solved as a function of timewhere , and we have used the ortho-normality relation and defined the self-interaction term . If , the integration over in can be carried out independently of , and for the interaction Hamiltonian () we havewhere is the dipole transition matrix element.
Probabilistic Model for the Excitation of TLS With Single Free Electron
In order to describe the TLS dynamics analytically, we present in this section an iterative approach to solving the source equation (). Substitutingon the RHS of () allows calculation of the development in time of the TLS, neglecting the recoil dynamics of the QEW. Multiplying Eq. () by the complex conjugate of the free electron wavefunction Eq. () and integrating over space, one obtainsWith , we have
This presentation is reminiscent of interaction with an unperturbed point-particle that arrives at time at the interaction point with Born’s quantum wavefunction probability . It should be stressed that is not well determined for a single electron. We assume that it is possible to solve Eq. () with substitution of its expectation value - , and the solution will then represent the result of interaction with an ensemble of identical QEWs.
The probability distribution of a single electron QEW of narrow width is:where is normalized over time. Then Eq. () can be simplified towhere is the weighed interaction strength, which for a Gaussian can be calculated from the wavepacket (). By neglecting the dynamics of the QEW around the interaction time also on the LHS, we can turn Eq. () into coupled differential equations for the TLS:and after integration
For a single Gaussian wavepacket at its longitudinal waist , we have
Normalizing time to the transit time parameter , where , and defining , the weighed interaction strength can be recast into a convolution relationwith the parameter being the ratio of the wavepacket duration and the transit time . Figure 2 shows the weighed interaction strength (20) between the TLS and free-electron for the longitudinal and vertical dipole moment orientations with different electron wavepacket sizes quantified by parameter . The maximum interaction strength decreases with the increase of the wavepacket size. The longitudinal component of the electric field induced by the traversing electron at the position of the dipole reverses the sign. This is the reason for the antisymmetric shape of the interaction strength as a function of time when the dipole moment is oriented longitudinally (Figure 2A), and its symmetry when it is oriented vertically (Figure 2B).
FIGURE 2
Defining Fourier transform , the TLS transition amplitude during the interaction (17) turns out to be proportional to the Fourier transform of the weighted function (19) at the transition resonance frequency:where
For a Gaussian QEW with , . Then the incremental transition probability amplitude in Eq. () reads
Therefore, the transition probability of TLS after the interaction iswhere , , and . For a finite-size QEW, the transition probability to the upper-level quantum state for excitation of the TLS from its ground state ( and ) iswhere , and the first-order term vanishes. In the case of excitation of the TLS from a superposition state, the second-order transition term is the same as Eq. (), while the first-order transition term in Eq () is:where with being the phase of the quantum state of the TLS.
The probabilistic model approximation is presumed to apply in the very short QEW regime , corresponding to the limit of short interaction time . Eqs. (, ) manifest through the parameter the wavepacket size dependence of the transition probabilities for excitation of the TLS from ground state or superposition state in the near-point-particle parameters regime. Note that for excitation from a superposition state, the incremental transition probability also depends on the relative phase of the quantum state of the TLS. The phase-match timing dependence on the dipole oscillation phase of the short interaction impulse at the QEW arrival time to the interaction point, is predictive of a coherent interaction enhancement by multiple electrons with correlated arrival timing, as discussed in the following section.
In order to check the validity of the analytical approximations, we have developed numerical computation codes for solving the FEBERI problem of interaction between a single finite-size QEW and a TLS at any initial state, starting from the Schrödinger equation (). The computation examples of the FEBERI effect were performed for a model of a Gaussian QEW and were studied as a function of its size , in order to examine the claimed dependence of the interaction on the wavepacket shape. In all the current examples the dipole polarization was taken to be transverse. The parameters used in the examples are typical of electron microscope PINEM-kind experiments [14], shown in Table 1. Note that the Hamiltonian of the combined system is time-independent so that the solution of the Schrödinger equation () can be represented by the density matrixwhere the evolution operator is and the initial state of the combined system is . Because we are interested now in the excitation of the TLS, the state of the TLS can be obtained by tracing out the free electron states
TABLE 1
| Physical Parameters | Typical Values |
|---|---|
| Beam Energy | |
| Free electron impact parameter | |
| TLS energy gap (transition frequency) | |
| Dipole moment |
The typical values of the simulation parameters.
The transition probability to the TLS upper quantum state is defined as . This procedure makes it possible to generalize the simulation of the FEBERI to multiple QEWs, in which case a train of QEWs sequentially interacts with the bound electron. After each interaction of an electron with the TLS, we trace out the free electron states to find the state of TLS Eq. (), which is used then as the initial state for the interaction with the next electron. We then repeat the calculation of Eq. (). The computation algorithm details are given in the appendix of Ref. [39]. This algorithm provides a fully quantum multiparticle interaction computation formulation, valid for rigorous evaluation of the TLS final state, and we employ it in the coming sections to examine the range of validity of the analytical expressions and the semiclassical FEBERI theory.
We first apply the quantum numerical computation to check the validity limits of the analytical expressions for the TLS transition probability for a single QEW (Eq. 24–Eq. 26). Figures 3A,B display the upper-level probability after interacting with a single QEW as a function of wave-packet size for the initial TLS in the ground state and superposition state, respectively. Numerical simulation results show that the transition probability is independent of the wave-packet size if the TLS starts from the ground state (blue curve), which contradicts the analytical result (Eq. ) of exponential decay (red curve). This discrepancy is expectable, because our probabilistic model approximation, neglecting the free-electron recoil, is presumed to apply only in the very short-size QEW regime. Indeed, the quantum numerical computation result of finite wavepacket size-independent transition probability from the ground state, seems to be more agreeable than the analytical result even from the philosophical point of view of Born’s probability interpretation of the electron wavefunction: when is large, the probability of the point-particle arrival to the TLS location is spread over a longer time, but it always happens at some time during the passage of the QEW, and must exhibit the same inelastic scattering, but the phase of the TLS dipole moment oscillation, undefined initially, is random after an interaction. On the other hand, when the initial state of the coherently pre-excited TLS is a quantum superposition state of a well-defined dipole oscillation phase, the quantum-mechanical numerical simulation result of the post-interaction probability is consistent with the analytical approximation expression of near point particle QEW with well-defined phase relative to the TLS dipole oscillation phase. In this case, the transition probability strongly depends on the wavepacket size, decaying with the increase of the wave-packet size in either model.
FIGURE 3
Excitation of TLS With a Bunched Electron Beam
The excitation of a TLS with multiple QEWs is theoretically an intricate multi-particle quantum interaction problem that involves the entanglement of the free electron wavefunction with the states of TLS. Here, we resort again to the simple analytic approximate probabilistic model, in which we extend Eq. () to multiple particleswhere . We then solve for the cumulative incremental transition probability for the case of periodically injected near-point-particle QEWs. Under the assumption that the relaxation time of the TLS is much longer than the duration of the QEWs pulse, we substitute the particles probability function in Eq. (), and changing order of integration in and results inWith change of variable ,
The incremental probability amplitude in this equation averages to zero for random , except when , where , namely, when the QEWs arrive at the interaction point at a rate that is a sub-harmonic of the transition frequency . Then, with the approximation of a small change in the amplitude , the amplitude of the upper level is
For a Gaussian QEW (), the transition probability to upper level for the QEWs case is
As we learned from the previous chapter, this approximate result may not be rigorous in the initial stage of the multiple electrons transition buildup from the ground state, when the phase of the dipole moment oscillation is not well-defined. We conjecture that when is large enough, the phase of the TLS gets established by the first near-point-particle QEWs of the train, and the subsequent QEWs then continue to build up the transitions in-phase.
This case of a periodically spaced train of near-point-particle QEWs may be realistic for low (microwave or THz) frequency TLS transitions, where classical Klystron-kind electron current modulation is available. It has thus been termed a “Quantum Klystron” in [30]. It can be comprehended as the quadratic approximation of the scaling of a Rabi oscillation process with Rabi frequency , and it is the analog of the classical bunched-particles beam superradiance effect [18] Note that in the classical point particle limit and low (microwave) frequencies [30] high current density of the electron beam is allowed (with the limitations of beam quality and space charge effect) and there may be then multiple electrons per period. We also point out that the case of multiple periodic trains of QEWs, is closely related to the earlier studied effect of “pulsed beam scattering” [46, 47].
It is instructive to compare the quadratic dependence of () on the number of QEWs - to the same dependence in the case of superradiance [18, 23]. In this comparison, the exponential decay factor that originates from the finite size of the Gaussian QEW (), is the quantum limit of the “bunching coefficient” in a bunched point-particle beam superradiance [18]. Note that when the QEWs arrival times t0k are random, then the second term in Eq. 31 averages to zero, but its averaged absolute value squares is proportional to N. This case is analogous to “shot-noise” spontaneous radiation emission by point-particle beam radiators.
Figure 4 shows a simulation of the buildup of the TLS upper-level probability with particles arriving at times , which is a random number. For the electrons arriving to the interaction point in-phase with the modulating laser modulo the bunching period , the probability growth is evidently quadratic, as claimed. For comparison, we show in the figure by the magenta curve the case that is taken to be entirely random. The growth rate is linear and the upper-level probability of arriving at the same value requires particles [29].
FIGURE 4
Excitation of a TLS With a Modulated Single Free Quantum Electron Wavepacket
Here we extend our Born’s probability interpretation of the electron wavefunction to model the case of a density modulated QEW. In Eq. () we model the density expectation value of a single electron wavepacket aswhere is the envelope function of the drifting QEW and the modulation function is periodic, and therefore can be expressed in terms of a Fourier series . The coefficients were derived in Ref. [22] for the case of the wavefunction of a modulated Gaussian QEW, which is the modulation phase, is determined by the modulating laser beam.
We next derive the incremental excitation probabilities in Eq. (). For a modulated QEW, the weighed probability function becomes
We substitute this probability distribution of a modulated QEW in Eq. () and change the order of integrationsWith a change of variables ,Then the first-order incremental transition probability is.
If the envelope Gaussian distribution is a wide function -, then the spectral function is a narrow function around a harmonic that is resonant with the transition frequency . In such a case, only one harmonic - can excite resonantly the transition:under the condition that the QEW modulation phase is matched to the TLS initial dipole moment oscillation phase .
Likewise, the second-order incremental transition probability can be calculated.
It can be seen that both first order and second-order expressions of the incremental probabilities are dependent on the QEW shape and modulation features and display resonant excitation characteristics around the condition , which would manifest the QEW modulation characteristics in a properly set experiment. Note that in a modulated QEW the QEW envelope is necessarily longer than the modulation period: . Comparing the first and second-order transition probability expressions of the modulated QEW (, ) to the corresponding terms of the unmodulated finite-size QEW (, ), we find out that at this limit the latter decay to zero, but the former (modulated QEW expressions) do not decay, as long as a harmonic of the modulation frequency is synchronous with the transition frequency –. This indicates a possibility for measuring the modulation features of the QEW. However, quite remarkably, for single modulated QEWs, there is no enhancement of the transition probability even at resonance. At resonance, Eqs. (, ) reduce to Eqs. (, ), except for a Fourier series component coefficient.
We check the results of the probabilistic model approximation by comparing it to the results of the quantum density matrix model numerical computation that is based on the solution of Schrodinger equation (), (described in the appendix of Ref. [39]). While the analytical expressions provide only post-interaction transition probabilities, the numerical computation, interestingly enough, lets us follow the dynamics of the TLS quantum transition probability during the interaction time (the QEW transit time at proximity to the TLS) and the final transition probability after passage. Figures 5B,C display computation results of the dynamics of quantum transition to the upper state of the TLS for two cases of pre-shaped QEW distributions (see Figure 5A): an unmodulated (broad) QEW () (red curve) and a modulated QEW with the same size envelope ( ) (blue curve). The density modulation was pre-evaluated assuming a PINEM process [20] with the beam parameters of Table 1 and .
FIGURE 5
Figures 5B,C present the dynamic buildup of the transition probability starting from the ground state and from a coherent superposition (qubit) state correspondingly. The bunching frequency was set to synchronize with the transition frequency at the fundamental harmonic . The dynamics of transitions with the modulated QEW suggest gradual “quantum jumps” any time a sub-bunch of the QEW arrives at the FEBERI interaction point. This seems to be consistent with the Born picture of the probability of point particle arrival. However, these are only probabilities of interaction events, the measurable post-interaction transition probability of the TLS from the ground state (Figure 5B) generated by a single passing-by resonantly modulated QEW is the same as with a single unmodulated near-point-particle QEW. This is in partial agreement with Eqs. (, ) except for the exponential decay factors in these two expressions that should be set equal to 1 because the analytical model, neglecting quantum recoil, is not rigorous for a long QEW (see also blue curve in Figure 3A). On the other hand, when starting from a coherent superposition state of the TLS (Figure 5C) the computation confirms the analytical expressions () (for modulated QEW at resonance) and () (for unmodulated long QEW): high transition probability in the first case and diminished transition probability in the latter (note that in this case, the simulation result reflects practically only the contribution of the dominant first-order incremental probability term because the second-order contribution is minute for a long unmodulated QEW).
Excitation of a TLS With a Beam of Modulation-Correlated Quantum Electron Wavepackets
In this section, we extend our Born’s probability interpretation analytical model to the case of multiple modulation-correlated QEWs. Consider the case of multiple long-size QEWs, all phase- (energy-) modulated at a PINEM interaction point at the level of their quantum wavefunctions [14] by the same coherent laser beam of frequency and phase . Passing all the same drift lengths to the FEBERI interaction point, the expectation values of their density are modulation phase correlated. We extend the probability distribution expression of a single electron QEW Eq. () to the modulation-correlated multiple particles by substitution in Eq. ()where are the centroid arrival times of the envelopes of the modulated QEWs, and the modulation function, common to all QEWs is periodic in time . For modulation-correlated QEWs, the weighed probability distribution function isThen, substitution in Eq. (), changing the integration order of and , results in.
Again, for broad Gaussian distribution , the spectral function is a narrow function around a harmonic that is resonant at the transition frequency . Take (upper and lower levels), then with the approximation of a small change in the amplitude , we haveThis averages to zero for random arrival times of the wavepacket centroids, except at the resonance case , where.And independently of the arrival times , the transition probability to the upper level is
This expression explicitly manifests the scaling buildup of the upper quantum level probability in the case of multiple modulation-correlated QEWs, similar to the case of periodically modulated point particles and in analogy to the superradiance of bunched particles [18].
To check the result of the probabilistic model approximation, we compare it to the results of the quantum model (Schrodinger equation) numerical computations (see details in appendix of Ref. [39]). For simulating the dynamics of multiple QEWs in interaction with the TLS, we consider a train of incoming QEWs, sequentially interacting with the TLS, and repeat the calculation for each electron passing by. Equation 7 is used to represent the initial quantum electron wavefunction of each electron upon arrival at the FEBERI interaction point. In Eq. 7 we use a random value (replacing ) corresponding to the random arrival of the QEW envelope centroids to the interaction point. The initial state of the TLS in each interaction is obtained by tracing out the free electron states in the post-interaction density matrix of the prior interaction. Possible overlap of QEWs is neglected, assuming a sparse electron beam. In Figure 6, we show the results of the simulation for the parameters given in Table 1. The modulation density of each QEW in the spatiotemporal domain is shown in Figure 5A. The wavepacket envelope size is . The simulation results, depicted in Figure 6 display the quadratic buildup (blue curve) of the upper quantum level probability in the case of multiple modulation-correlated QEWs at the resonant fundamental harmonic FEBERI interaction condition . Even though the QEW centroids arrive at random times , the phase (in Eq. 7) of all QEWs is the same, determined by the phase of the coherent laser beam that pre-modulates them all by a PINEM process. The noisy linear buildup curve (red) corresponds to the case of PINEM interaction with an incoherent light, where the phase (replacing in Eq. 7) is random, as well as the envelope centroid arrival times .
FIGURE 6
Conclusion
We studied the excitation of a bound electron in an arbitrary initial quantum (qubit) state by a free electron beam, where the bound electron is modeled as a quantum two-level system (TLS) and the free electrons are represented as quantum electron wavepackets (QEWs). We developed an analytical model for the FEBERI interaction based on Born’s probabilistic interpretation of the quantum electron wavefunction and derived expressions for the transition probability of the TLS subject to interaction with a shaped QEW (finite size gaussian), optically modulated QEWs, and a beam of such QEWs. We tested the analytical results against the simulation results of an accurate quantum model of the FEBERI effect, based on the density matrix solution of the Schrodinger equation for the entangled free-bound electron wavefunctions.
The accurate quantum model simulations show quadratic growth of the TLS quantum transition probability as a function of the number of electrons in a train of modulation-correlated QEWs, when a harmonic of the QEW modulation frequency is resonant with the TLS quantum transition frequency, and the modulation initial phase of all electrons is the same (e.g., they are modulated by the same coherent laser beam in a PINEM process). This result confirms the prediction of the analytical approximation model based on Born’s probabilistic interpretation of the electron wavefunction and the earlier semiclassical analysis of the FEBERI effect [26].
The quadratic scaling of the TLS excitation build-up in the FEBERI process can be explained as the second-order expansion of the sinus-squared scaling of a Rabi oscillation process [26], where the optical frequency near-field of the beam of modulation-correlated QEWs play the same role as a laser beam in a conventional Rabi oscillation experiment. Another instructive observation that came out of the numerical simulations is that in the case of modulation-correlated QEWs beam, the temporal intervals between the probability-density micro bunches of the QEWs (Figure 5) are spaced apart at resonance an integral number of times the TLS transition period ( ) even though the centroids of the QEWs arrive at random. The coherent quadratic buildup can be then viewed in Born’s probabilistic interpretation as a result of quasi-periodic arrival times of a train of quantum- probability-determined “point-particles” in-phase with the TLS dipole moment oscillation at the quantum transition frequency. This is a bridge to the semiclassical case of FEBERI interaction with a point-particle density-modulation beam (quantum klystron [30]) and an analog of the radiative process of bunched electron beam superradiance [18].
To test the surprising result of coherent transition probability buildup with a modulation-correlated electron beam independently of the random arrival times of the electrons to the FEBERI point, we show in Figure 6 (red curve) a simulation of interaction with a beam of modulated QEWs that are not modulation phase-correlated (random ). In this case, that corresponds to electron beam PINEM modulation with an incoherent laser beam, the FEBERI scaling is a noisy linear curve, indicating uncorrelated excitation of the TLS. This observation is in line with recent observations that PINEM modulation carries the coherence properties and even the quantum statistical fingerprint of the modulating light beam and transfers it to the cathodoluminescence spectrum [25] and to the EELS spectrum [24].
We point out that coherent control and demonstration of Rabi oscillation of a single TLS with the proposed FEBERI effect is attractive in comparison to such operations with a laser beam, because of the atomic-scale resolution of an electron beam. However, the experimental realization of this scheme is challenging due to the very small value of the interaction coupling factor in typical TLS targets (). Enhancement of the interaction may be possible in consideration of high dipole moment targets, such as multiple TLS in a Dicke (“super-qubit”) state [35].
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
AG, AY, JS, and DR conceived the concept. BZ and DR performed the theoretical derivations and prepared the figures. BZ, DR, and AG wrote the paper with contributions from RI and AF All authors reviewed and discussed the manuscript and made significant contributions to it.
Funding
This work was supported by the Israel Science Foundation (ISF) under Grant No. 00010001000, the National Natural Science Foundation of China under Grant No. 12104068, and the Natural Science Foundation of Chongqing under Grant No. cstc2021jcyj-msxmX0684.
Acknowledgments
We acknowledge helpful communications with O. Kfir, Y. M. Pan, and F. J. Garcia de Abajo. DR acknowledges support by the PBC program of the Israel Council of Higher Education.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
VerbeeckJTianHSchattschneiderP. Production and Application of Electron Vortex Beams. Nature (2010) 467(7313):301–4. 10.1038/nature09366
2.
Voloch-BlochNLereahYLilachYGoverAArieA. Generation of Electron Airy Beams. Nature (2013) 494(7437):331–5. 10.1038/nature11840
3.
ShilohRLereahYLilachYArieA. Sculpturing the Electron Wave Function Using Nanoscale Phase Masks. Ultramicroscopy (2014) 144:26–31. 10.1016/j.ultramic.2014.04.007
4.
BaumP. Quantum Dynamics of Attosecond Electron Pulse Compression. J Appl Phys (2017) 122(22):223105. 10.1063/1.5006864
5.
MorimotoYBaumP. Single-cycle Optical Control of Beam Electrons. Phys Rev Lett (2020) 125(19):193202. 10.1103/physrevlett.125.193202
6.
VanacoreGMMadanICarboneF. Spatio-temporal Shaping of a Free-Electron Wave Function via Coherent Light-Electron Interaction. Riv Nuovo Cim (2020) 43(11):567–97. 10.1007/s40766-020-00012-5
7.
PiazzaLLummenTTQuiñonezEMurookaYReedBWBarwickBet alSimultaneous Observation of the Quantization and the Interference Pattern of a Plasmonic Near-Field. Nat Commun (2015) 6(1):6407–7. 10.1038/ncomms7407
8.
ReinhardtOKaminerI. Theory of Shaping Electron Wavepackets with Light. ACS Photon (2020) 7(10):2859–70. 10.1021/acsphotonics.0c01133
9.
YaluninSVFeistARopersC. Tailored High-Contrast Attosecond Electron Pulses for Coherent Excitation and Scattering. Phys Rev Res (2021) 3(3):L032036. 10.1103/physrevresearch.3.l032036
10.
BarwickBFlanniganDJZewailAH. Photon-induced Near-Field Electron Microscopy. Nature (2009) 462(7275):902–6. 10.1038/nature08662
11.
VanacoreGMMadanIBerrutoGWangKPomaricoELambRJMcGroutherDKaminerIBarwickBGarcía de AbajoFJCarboneF. Attosecond Coherent Control of Free-Electron Wave Functions Using Semi-infinite Light fields. Nat Commun (2018) 9(1):2694–11. 10.1038/s41467-018-05021-x
12.
KozákMSchönenbergerNHommelhoffP. Ponderomotive Generation and Detection of Attosecond Free-Electron Pulse Trains. Phys Rev Lett (2018) 120(10):103203. 10.1103/physrevlett.120.103203
13.
KozákMEcksteinTSchönenbergerNHommelhoffP. Inelastic Ponderomotive Scattering of Electrons at a High-Intensity Optical Travelling Wave in Vacuum. Nat Phys (2018) 14(2):121–5. 10.1038/nphys4282
14.
FeistAEchternkampKESchaussJYaluninSVSchäferSRopersC. Quantum Coherent Optical Phase Modulation in an Ultrafast Transmission Electron Microscope. Nature (2015) 521(7551):200–3. 10.1038/nature14463
15.
PriebeKERathjeCYaluninSVHohageTFeistASchäferSet alAttosecond Electron Pulse Trains and Quantum State Reconstruction in Ultrafast Transmission Electron Microscopy. Nat Photon (2017) 11(12):793–7. 10.1038/s41566-017-0045-8
16.
BaumP. Quantum Dynamics of Attosecond Electron Pulse Compression. J Appl Phys (2017) 122(22):223105. 10.1063/1.5006864
17.
MorimotoYBaumP. Single-cycle Optical Control of Beam Electrons. Phys Rev Lett (2020) 125(19):193202. 10.1103/physrevlett.125.193202
18.
GoverAIanconescuRFriedmanAEmmaCSudarNMusumeciPet alSuperradiant and Stimulated-Superradiant Emission of Bunched Electron Beams. Rev Mod Phys (2019) 91(3):035003. 10.1103/revmodphys.91.035003
19.
GoverAPanY. Dimension-dependent Stimulated Radiative Interaction of a Single Electron Quantum Wavepacket. Phys Lett A (2018) 382(23):1550–5. 10.1016/j.physleta.2018.03.049
20.
BlackDSNiedermayerUMiaoYZhaoZSolgaardOByerRLet alNet Acceleration and Direct Measurement of Attosecond Electron Pulses in a Silicon Dielectric Laser Accelerator. Phys Rev Lett (2019) 123(26):264802. 10.1103/physrevlett.123.264802
21.
SchönenbergerNMittelbachAYousefiPMcNeurJNiedermayerUHommelhoffP. Generation and Characterization of Attosecond Microbunched Electron Pulse Trains via Dielectric Laser Acceleration. Phys Rev Lett (2019) 123(26):264803. 10.1103/physrevlett.123.264803
22.
PanYGoverA. Spontaneous and Stimulated Radiative Emission of Modulated Free-Electron Quantum Wavepackets-Semiclassical Analysis. J Phys Commun (2018) 2(11):115026. 10.1088/2399-6528/aae2ec
23.
PanYGoverA. Spontaneous and Stimulated Emissions of a Preformed Quantum Free-Electron Wave Function. Phys Rev A (2019) 99(5):052107. 10.1103/physreva.99.052107
24.
DahanRGorlachAHaeuslerUKarnieliAEyalOYousefiPSegevMArieAEisensteinGHommelhoffPKaminerI. Imprinting the Quantum Statistics of Photons on Free Electrons. Science (2021) 373(6561):eabj7128. 10.1126/science.abj7128
25.
KfirODi GiulioVde AbajoFJGRopersC. Optical Coherence Transfer Mediated by Free Electrons. Sci Adv (2021) 7(18):eabf6380. 10.1126/sciadv.abf6380
26.
GoverAYarivA. Free-Electron-Bound-Electron Resonant Interaction. Phys Rev Lett (2020) 124(6):064801. 10.1103/PhysRevLett.124.064801
27.
de AbajoFJG. Comment on “Free-Electron–Bound-Electron Resonant Interaction”. Phys Rev Lett (2021) 126(1):019501. 10.1103/physrevlett.126.069902
28.
de AbajoFJG. Erratum: Comment on “Free-Electron–Bound-Electron Resonant Interaction”[Phys. Rev. Lett. 126. Phys Rev Lett (20212021) 126(6):069902:019501. 10.1103/physrevlett.126.069902
29.
GoverAYarivA. Gover and Yariv Reply. Phys Rev Lett (2021) 126(1):019502. 10.1103/physrevlett.126.019502
30.
RätzelDHartleyDSchwartzOHaslingerP. Controlling Quantum Systems with Modulated Electron Beams. Phys Rev Res (2021) 3(2):023247. 10.1103/physrevresearch.3.023247
31.
ZhaoZSunX-QFanS. Quantum Entanglement and Modulation Enhancement of Free-Electron-Bound-Electron Interaction. Phys Rev Lett (2021) 126(23):233402. 10.1103/physrevlett.126.233402
32.
García de AbajoFJDi GiulioV. Optical Excitations with Electron Beams: Challenges and Opportunities. ACS Photon (2021) 8(4):945–74. 10.1021/acsphotonics.0c01950
33.
RuimyRGorlachAMechelCRiveraNKaminerI. Toward Atomic-Resolution Quantum Measurements with Coherently Shaped Free Electrons. Phys Rev Lett (2021) 126(23):233403. 10.1103/physrevlett.126.233403
34.
ZhangBRanDIanconescuRFriedmanAScheuerJYarivAet alQuantum Wave-Particle Duality in Free-Electron-Bound-Electron Interaction. Phys Rev Lett (2021) 126(24):244801. 10.1103/physrevlett.126.244801
35.
ZhangBRanDIanconescuRFriedmanAScheuerJYarivAet alQuantum States Interrogation Using a Pre-shaped Free Electron Wavefunction. Phys. Rev. Research (Forthcoming 2022).
36.
DickeRH. Coherence in Spontaneous Radiation Processes. Phys Rev (1954) 93(1):99–110. 10.1103/physrev.93.99
37.
GoverA. Superradiant and Stimulated-Superradiant Emission in Prebunched Electron-Beam Radiators. I. Formulation. Phys Rev Spec Topics-Accelerators Beams (2005) 8(3):030701. 10.1103/physrevstab.8.030701
38.
García de AbajoFJDi GiulioV. Optical Excitations with Electron Beams: Challenges and Opportunities. ACS Photon (2021) 8(4):945–74. 10.1021/acsphotonics.0c01950
39.
GoverAZhangBRanDIanconescuRFriedmanAScheuerJet alResonant Interaction of Modulation-Correlated Quantum Electron Wavepackets with Bound Electron States (2020). arXiv preprint arXiv:2010.15756.
40.
Di GiulioVKociakMde AbajoFJG. Probing Quantum Optical Excitations with Fast Electrons. Optica (2019) 6(12):1524–34. 10.1364/optica.6.001524
41.
BartlettRJ. Many-body Perturbation Theory and Coupled Cluster Theory for Electron Correlation in Molecules. Annu Rev Phys Chem (1981) 32(1):359–401. 10.1146/annurev.pc.32.100181.002043
42.
García de AbajoFJ. Optical Excitations in Electron Microscopy. Rev Mod Phys (2010) 82(1):209–75. 10.1103/revmodphys.82.209
43.
Feynman RichardPMainly Electromagnetism and Matter, Feynman Lectures on Physics, Vol. 2. 33. Available at: https://www.feynmanlectures.caltech.edu/II_26.html.
44.
BreitG. Dirac's Equation and the Spin-Spin Interactions of Two Electrons. Phys Rev (1932) 39(4):616–24. 10.1103/physrev.39.616
45.
LhuillierCFarouxJP. Hamiltonian of a many-electron Atom in an External Magnetic Field and Classical Electrodynamics. J Phys France (1977) 38(7):747–55. 10.1051/jphys:01977003807074700
46.
RobicheauxFNoordamLD. Coherent Scattering with Pulsed Matter Beams. Phys Rev Lett (2000) 84(17):3735–9. 10.1103/physrevlett.84.3735
47.
PindzolaMSWitthoeftMRobicheauxF. Coherent Electron Beam Excitation of a Model Hydrogen Atom. J Phys B: Mol Opt Phys (2000) 33(21):L839–L846. 10.1088/0953-4075/33/21/107
Summary
Keywords
quantum electron wavepackets, superradiance, coherent control, wavepackets size, electron-matter interaction
Citation
Ran D, Zhang B, Ianconescu R, Friedman A, Scheuer J, Yariv A and Gover A (2022) Coherent Excitation of Bound Electron Quantum State With Quantum Electron Wavepackets. Front. Phys. 10:920701. doi: 10.3389/fphy.2022.920701
Received
14 April 2022
Accepted
23 May 2022
Published
26 July 2022
Volume
10 - 2022
Edited by
Giovanni Maria Vanacore, University of Milano-Bicocca, Italy
Reviewed by
Dawei Wang, Zhejiang University, China
Jiteng Sheng, East China Normal University, China
Updates
Copyright
© 2022 Ran, Zhang, Ianconescu, Friedman, Scheuer, Yariv and Gover.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Avraham Gover, gover@eng.tau.ac.il
†These authors have contributed equally to this work and share first authorship
This article was submitted to Optics and Photonics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.