Dynamical scattering of single-electron wave packets

Characterizing distinct electron wave packets is a basic task for solid-state electron quantum optics with applications in quantum metrology and sensing. A important circuit element for this task is a non-stationary potential barrier than enables backscattering of chiral particles depending on their energy and time of arrival. Here we solve the quantum mechanical problem of single-particle scattering by a ballistic constriction in an fully depleted quantum Hall system under spatially uniform but time-dependent electrostatic potential modulation. The result describes electrons distributed in time-energy space according to a modified Wigner quasiprobability distribution and scattered with an energy-dependent transmission probability that characterizes constriction in the absence of modulation. Modification of the incoming Wigner distribution due to external time-dependent potential simplifies in case of linear time-dependence and admits semiclassical interpretation. Our results support a recently proposed and implemented method for measuring time and energy distribution of solitary electrons as a quantum tomography technique, and offer new paths for experimental exploration of on-demand sources of coherent electrons.


I. INTRODUCTION
Electron quantum optics is a relatively new field, aiming to reproduce quantum optics experiments with electron wave packets instead of photons [1][2][3][4]. It offers the prospects of probing the interaction between just a few electrons, as well as studying phenomena on the scale of electron coherence time. The field has potentially promising applications in signal processing [5] and quantum sensing [6].
One of the crucial ingredients that has made the investigation of single electron excitations in ballistic waveguides possible is the advent of devices that emit ordered streams of electrons with sufficient separation between individual particles [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]. Techniques to characterise quantum properties of electrical current have been adapted from photon quantum optics. Statistical properties of the source and the exchange statistics of the particles can be probed using intensity interferometers such as the Hanbury-Brown and Twiss interferometer [22][23][24][25], while coherence, entanglement, and the wave-like nature of particles can be probed using amplitude interferometers such as the Mach-Zehnder interferometer [26][27][28].
However, most of the existing electron quantum optics experiments have focused on single-electron sources that emit electrons close to the Fermi sea [3]. These excitations can be accessed by perturbing the Fermi sea through the application of periodic gate potentials allowing for quantum state reconstruction through correlation measurements between the unknown and the reference signal [29][30][31]. Such an approach is not possible in the case of "high-energy" electrons [19,[32][33][34][35][36] emitted far above (tens of meV) the Fermi sea, because they do not overlap with perturbations around the Fermi energy. These excitations provide a new kind of quasi-particle whose coherence properties are largely unexplored.
A recent work by Fletcher et al. [37] addresses the characterization of these "high-energy" electrons by using an en- * Corresponding author: slava@latnet.lv FIG. 1. A sketch of the setup: two counterpropagating chiral edge states (orange) in a quantum Hall system pass through a constriction where backscattering between the edges is possible. Individual electronic excitations are launched from the source (left) and get either transmitted or reflected. The scattering region is electrostatically gated and the energy-dependent transmission probability T (E) can be modulated by the voltage on the gate (shaded shape). It is assumed that the potential V (t) created by the gate voltage modulation is spatially uniform in the scattering region where backscattering between the edge channels takes place. Tuning the time-and energydependence of the transmission function T (E − V (t)) by appropriately designed time-dependence of the gate potential V (t) enables tomography -the measurement of joint energy-time content -of the wavepacket. ergy barrier that is tuned to match the energy of the quasiparticles. The authors propose and implement an electron tomography protocol that reconstructs the joint energy-time quasiprobability distribution p(E, t) of the incoming electrons by measuring the transmitted charge through an energy barrier whose height is varied linearly in time. The experimental setup can be described by a model shown schematically in Fig.  1, where the time-dependent energy barrier is modelled with a static constriction subject to a (locally) uniform gate voltage V (t) = V 0 + αt. The two parameters of the linear modulation (the offset in energy V 0 and the slope α in the time-energy plane) provide a two-dimensional map of the charge Q(V 0 , α), that is used to infer p(E, t) using inversion techniques from tomographic image processing [37]. The key relation connecting the properties of the incoming wave-packets with the measured signal can be written as where T (E) is the energy-dependent transmission probability of the constriction in the absence of the gate voltage modulation. * Equation (1) has been previously derived classically [36] by considering scattering of a statistical ensemble of electrons with simultaneously well-defined energy and time of arrival, impinging on a scattering barrier with a transmission probability T (E, t). Due to gauge invariance under uniform modulation, shift in voltage is equivalent to shift in energy, giving Eq. (1). In this classical picture p(E, t) is simply the joint probability density characterizing the electrons emitted by the source.
In this work, we derive a general expression for the transmitted charge for an arbitrary time-dependence of V (t) and show that the distribution p(E, t) in a general case needs to be replaced by a suitably modified Wigner function. We also show that if the gate voltage V (t) has a linear dependence on time, no modification of the Wigner function is necessary and Eq. (1) is valid if we take p(E, t) to be the Wigner function of the incoming electron ρ in (E, t). This confirms that the protocol implemented in [37] can be used not only for classical but also for quantum tomography, i.e., quantum state reconstruction. We argue that if there are deviations from a linear time dependence of the gate potential, the modified Wigner function instead of the actual incoming Wigner function will be observed. We note that the Wigner distribution function has previously been found to be a useful concept in electron quantum optics [30,38,39]. Here it appears naturally as a quantum counterpart of the classical probability density.
Our manuscript is organized as follows: In Sec. II we give a precise definition of our system and of the approximations involved, as well as introduce the Wigner distribution function. Our main result, a quantum version of Eq. (1), is presented and derived in Sec. III. In that Section we also discuss special cases, such as a gate voltage with a linear time dependence (for which Eq. (1) is exact), a gate voltage with sharp edges, and the limit of a slowly varying gate voltage. Section IV illustrates our results, by comparing the "classical" and "quantum" expressions for a few characteristic examples. In Section V we discuss the implication of our results for the analysis and improvement of electron tomography experiments. We conclude in Sec. VI.

II. MODEL
We consider a constriction with two counterpropagating quantum Hall edge channels. The coordinate x measures the * Strictly speaking, in Eq. (1) V (t) should be replaced by the integral of the electric field for electrons traveling along the edge, which need not be equal to V (t) for a time-dependent gate potential, see Sec. III.
where v is the electron velocity,σ 3 is the Pauli matrix, and the 2 × 2 matrix structure refers to the two edges. The second termĤ bs describes backscattering in the constriction, and is a gate potential. The backscattering Hamiltonian H bs is characterized by its transmission probability in the absence of the gate voltage V (x, t), where τ (E) is the transmission amplitude. Since backscattering takes place locally at the contact, the backscattering term has finite support which we define to be |x| < x b . We furthermore assume that This assumption is justified if backscattering exclusively takes place near the center of the constriction, which is where the two counterpropagating edge states are closest, and where the externally controlled part of the potential is spatially uniform. Although, strictly speaking, the condition (5) is stronger than the assumption that a (time-independent) shift of the gate voltage V is equivalent to a shift in energy, the relation T (E, V ) = T (E − V ) (which is central to the tomography experiment [37]) does not hold for a generic scatterer if the condition (5) is not satisfied.
Finally, we set the gate potential to zero sufficiently far away from the constriction, We consider a wavepacket incident on the constriction, such that it reaches the center of the constriction for t ≈ 0. For times t → −∞ this initial state is described by the wavefunction where the function ϕ in (t) is peaked near t = 0. The initialstate wavefunction is normalized such that Similarly, for times t → ∞ the wavepacket is described by the wavefunction where ϕ out (t − x/v) represents the transmitted part of the wave packet, and ϕ out (t + x/v) represents the reflected part. Instead of the function ϕ in (t) it is convenient to use the Wigner transform ρ in (E, t), which is defined as The classical limit of the Wigner distribution function ρ(E, t) is the joint probability density p(E, t) of energy E and time t and corresponds to a Dirac delta distribution peaked around the classical trajectory E(t). Correspondingly, the spread of the Wigner distribution around the classical trajectory is a manifestation of the quantumness of the wave packet [39]. The interpretation of ρ(E, t) as a joint probability distribution is not rigorous, though, since the Wigner distribution can also take on negative values. Integrating ρ(E, t) with respect to one of its arguments, however, one obtains a positive function, which is the probability density (marginal distribution) of the other argument, i.e.,

A. Derivation of the main result
We now proceed with the calculation of the charge Q transmitted through the constriction, For this calculation it is sufficient to consider the right-moving edge only. We drop the spinor notation of the previous Section, use the scalar wavefunction ψ(x, t) to refer to the wavefunction component at the right-moving edge, and write In the absence of the gate voltage, V (x, t) = 0, the result can be easily expressed in terms of the energy-dependent transmission probability T (E), which is a special case of the classical equation (1), with p(E, t) replaced by ρ in (E, t).
We now consider the general case of time-and energydependent scattering, specializing to the geometry described in Sec. II, for which the time dependence comes from the gate potential V (x, t).
Since the potential V (x, t) = 0 for x < −x g , the initialstate solution (7) is valid for all x < −x g . Similarly, the expression (9) for the transmitted wavepacket is valid for all x > x g . In the first step of our calculation we solve the timedependent Schrödinger equation to calculate the wavefunction at all positions In a second step we perform the gauge transformation where V (t) is the value of the gate voltage at the center of the constriction, see Eq. (5). For −x b < x < x b the transformed wavefunctionψ satisfies a Schrödinger equation without gate potential, so that where is the Fourier transform of the constriction's transmission amplitude τ (E). The choice of the offsets in the time arguments in Eq. (16) ensures that, for a perfectly transmitting constric- -all the information about coherent dynamics inside the modulated scattering region is encoded in the matrix element τ (t) of the time evolution operator which is time-translation invariant in the gauge expressed by Eq. (15). This is valid as long as the modulation is spatially flat as expressed by condition (5). Matching Eqs. (14) and (16) Transmitted charge (12) can be written in the form that closely resembles the classical result (1) and the limiting case (13), but with the functioñ Here is the integral of the electric field along the electron's trajectory. The expressions (20)- (22) are the key results of this article. The modified Wigner functionρ(E, t) can be interpreted as the Wigner representation (10) of a modified incoming asymptotic state,φ in (t) = ϕ in (t) e −i t V(t )dt / . The effective gate potential V(t) represents the relevant modification of the applied gate potential V (t) due to finite spatial extent of the gate. Althoughφ in (t), and henceρ(E, t), does not in general represent the actual quantum state of the electron at any time instant,ρ(E, t) can still be measured as the outcome of a tomographic reconstruction, as we argue in Sec. V.

B. Special cases
Below we consider a number of special cases of V (x, t) and discuss the corresponding forms of V(t) andρ(E, t).
a. Time-independent potential. One verifies that for a time-independent potential ∂V ( where V is the value of the spatially uniform potential in the scattering region. We see that the time-independent V just changes the energy reference level in Eq. (13), as expected from gauge invariance.
Typically, u(x) is expected to be a smooth sigmoid function, in which case the convolution (25) limits the sharpness of time-dependent features in V(t) compared to V (t). c. Gate potential with sharp edges. In the limit of a sharp edge, is applicable, and the correction of V (t) becomes just the shift of the argument, We see that the modification of the wave-packet happens at the edge of the gate-affected region. d. Linear time dependence. The experiment of Ref. 37 features a potential with a linear time dependence, which yields a linear V(t), In this case we find that the functionρ(E, t) can be expressed directly in terms of the Wigner distribution ρ in (E, t) of the incoming wavepacket, Upon comparing Eqs. (20) and (29) with Eq. (1) we conclude that for a potential with linear time dependence the quantummechanical theory reproduces the classical approximation (1). Note that V(t) is not in general equal to the value of the potential V (x = 0, t) = V (t) at the centre of the gate at the same time instant t, but rather at an earlier time corresponding to the effective position of the gate edge. We can define the latter explicitly by rewriting Eq. (28) as e. Slowly varying potential. If the rate of change of the gate voltage is slow compared to the velocity of the wave packet, it is reasonable to expand Truncating after the quadratic term, this leads to the formal expressioñ where Equation (31) contains the result (29) for a gate potential with a linear time dependence and is a good starting point for an expansion in a small second derivative V 2 (x) of the timedependent gate potential. f. Linear time dependence superimposed on an arbitrary gate potential V (x, t). If on top of a linear time dependence the potential V (x, t) also contains arbitrary additional terms V (x, t), the modified Wigner functionρ(E, t) can be expressed as where V(t) is given by Eq. (28) andρ (E + V(t), t) is calculated with respect to V (x, t) only. The relevance of this result will be discussed in Sec. V. (Note that Eq. (34) simplifies to Eq. (29) for the special case V (x, t) = 0.)

IV. EXAMPLES
In this Section we show explicit results for three examples: A gate voltage with linear time dependence, a gate voltage with an abrupt step-like time dependence, and a gate voltage with parabolic time dependence. We take the spatial profile of the gate voltage V (x, t) to be a spatially uniform potential with sharp edges at x = ±x g . We also shift the time origin by −x b /v, so that the integral V(t) of the electric field can be replaced directly by the potential V (t) at the center of the constriction, see Eq. (26). For the incoming wavepacket we take the uncorrelated Gaussian form, The scale σ t sets the characteristic width in the time domain. The energy of the wavepacket is centered around E 0 , with fluctuations of order /σ t . The Wigner distribution ρ in (E, t) corresponding to the wavefunction (35) is For the examples in the this and the following sections, we model the energy-dependence of the transmission probability through a static constriction as where δ gives the width of the energy window in which the transmission changes from 0 to 1 and the zero of energy is chosen to coincide with half-transmission point of the constriction. For the numerical calculations we will take the idealized limit δ → 0, corresponding to a point contact that perfectly selects the sign of the electron's energy.

A. Linear ramp
We first consider the case of a gate voltage with a linear time dependence, Gate voltages with different offsets V 0 are related by a delay time t d = −V 0 /α. For this case the modification of  the Wigner distribution can be obtained from Eq. (29), which gives:ρ The modified Wigner distribution for different values of the ramp rate α and the offset V 0 is shown in Figure 3, along with the non-transformed Wigner distribution. In comparison to the original Wigner distribution ρ in , the modified distributioñ ρ has a shifted center energy, determined by the offset V 0 , and it is stretched along a straight line with slope −α. Figure 4 shows the transmitted charge as a density plot in the case of a linear ramp as a function of the parameters α and V 0 . This transformation gives the basis for interpreting p(E, t) in Eq.
(1) as the ρ in (E, t), as discussed in the discussed in Sec. III. B.
Step-like ramp As mentioned in the introduction, performing the detection measurement with an energy-independent transmission probability that abruptly switches on or off allows one to measure the electron probability distribution as a function of time. Such an instantaneous switching on or off of the transmission function of the constriction requires an increase of the gate voltage V (t) by an amount much larger than the energy uncertainty /σ t of the incoming wavepacket. Richer information about time and energy distribution can be obtained if the gate voltage jumps by a finite amount at time t = t 0 , For a Gaussian incoming wavepacket, the corresponding modified Wigner distribution can be calculated from Eqs. (21) and (26), where we abbreviated A comparison between the initial-state Wigner function ρ in (E, t), the modified Wigner functionρ(E, t) for two different parameter choices, and the semiclassical expectation † † We call ρ in (E + V (t), t) "semiclassical" because it can be interpreted as a quasi-probability density for electrons that have a well-defined energy E + V (t) at every time instant t. Fig. 5. Because of the abrupt time dependence at t = t 0 , the semiclassical expectation is a poor approximation, as can be seen in the figure. In particular, the true modified distribution function is a continuous function of E and t, whereas the semiclassical expectation has a discontinuous jump by V f − V i at t = t 0 . The two lower panels of the figure show that in contrast to a linear ramp, a step-like ramp transforms the Wigner distribution such that it assumes negative values. In Fig. 6 we show the transmitted charge at a fixed step height V f − V i = 3 /σ t as a function of the center potential result (20), and the top right panel the semiclassical approximation (1). Their difference is shown in the bottom panel. As expected, if |V | is sufficiently large, the transmitted charge is independent of t 0 and approaches 0 or 1 (corresponding to maximum possible transmittance). For intermediate values of V , there is a transition from 0 to 1 transmitted charge as t 0 goes from large negative to large positive values. Although the exact and the semiclassical results both reproduce the correct limits for large |V | and |t 0 |, the behavior for smallV and t 0 shows qualitative differences (such as shape of the median Q = e/2 curve), as well as quantitative differences (more than 0.1e, see bottom panel).
We can check explicitly that a large-amplitude sudden step in V (t) will sample the time distribution. For this we need to prove that for −V i , V f /σ t , δ, finite E 0 , andρ(E, t) given by Eq. (41), the transmitted charge Of the two terms in Eq. (41), the contribution of the second one to the time integral in Eq. (43) vanishes the limit of large V f − V i due to fast oscillations, while the first term contributes only at t < t 0 , when integrand is non-zero in the vicinity of E ≈ −V i . Taking the limit V i → −∞, this gives where z = E √ 2σ t / + i(t 0 − t)/( √ 2σ t ) and w(z) is Faddeeva function w(z) = (i/π) e −ξ 2 /(z − ξ) dξ (Im z > 0).  (20), the center column shows the semiclassical approximation (29), and right column shows the slow-potential approximation, obtained by truncating the exponential in Eq. (31) at first order in W .
The first term in curly brackets in Eq. (44) gives the expected result. Straightforward integration confirms that the second term, proportional to w(z), contributes zero.

C. Parabolic ramp
The third example is the case of a parabolic ramp, which may serve as an approximation of a periodic (sinusoidal) gate potential near the maximum or minimum of the potential. We consider a gate potential of the form The modified Wigner function is where Ai(z) is the Airy function of the first kind. Figure 7 contains a comparison between the exact result for the modified Wigner function, the semiclassical approximation (29) and a slow-potential approximation obtained by FIG. 8. Transmitted charge Q for a Gaussian wavepacket incident on a constriction subject to a gate potential with the parabolic time dependence (45) with W = 0.2 /σ 3 t (left) and W = 2 /σ 3 t (right). The top row is obtained from the full quantum-mechanical expression Eq. (20), the center row uses the slow-potential approximation (31) with the exponential truncated after first order in the second derivative of the potential. The bottom row shows the difference between the slow-potential approximation and the full quantummechanical result.
truncating the exponential in Eq. (31) at first order in the second derivative of the potential, for different values of W and for V 0 = E 0 + /σ t , t 0 = 0. The figure confirms that both the semiclassical and the slow-potential approximation are good approximations for the modified Wigner functionρ(E, t) if |W | /σ 3 t , whereby the slow-potential approximation also faithfully reproduces some of the fringes at larger values of W , which are absent from the semiclassical approximation.
In Fig. 8 we show the transmitted charge for two different values of the the curvature parameter W , with separate panels for the exact result, the semiclassical approximation, and the slow-potential approximation. Their difference is shown in the bottom panels. The transmitted charge goes from zero for values of V 0 − E 0 far below the line V 0 − E 0 = 1 2 W t 2 0 to 1 for values of V 0 − E 0 far above this line. The differences between the exact and the slow-potential solution grows with W , remaining below 1% for W = 0.2 /σ 3 t and going up to 10% for W = 2 σ 2 t .

V. CONNECTION TO EXPERIMENT
In this section we discuss the relevance of our results to the recently proposed and implemented tomography protocol for solitary electrons [37].
A. Scattering by a linearly modulated barrier as a Radon transform As mentioned in the introduction, the tomography experiment provides a map of the charge Q(V 0 , α) transmitted through barrier that is controled by a gate voltage V (t) = αt + V 0 with a linear time dependence as a function of the energy offset V 0 and the energy-time slope α = α 0 tan θ, where α 0 is a scale factor adapted to the characteristic time and energy scales of the experiment. The quantity represents the Radon transform [40] of the incoming Wigner distribution ρ in (E, t) convolved with the energy derivative of the transmission coefficient dT /dE. In Eq. (47) we have used the fact that, according to Eq. (29), for a linear ramp the modified Wigner functionρ(E, t) can be expressed in terms of the incoming Wigner function ρ in (E, t) simply by adding V (t) to the energy. The tomography protocol of Ref. [37] consists of measuring (47) for a sufficiently wide range of V 0 and θ to enable numerical computation of the the inverse Radon transform using the standard filtered back-projection algorithm [40]. A sufficiently sharp T (E), such that the derivative dT /dE can be approximated as a delta function, ensures that the inversion accurately represents the unknown ρ in (E, t).
To illustrate the Radon transform, in Fig. 9 we have shown the quantity ∂Q(V 0 , θ)/∂V 0 calculated for the incoming Wigner distribution ρ in (E, t) of a Gaussian wave packet, exactly like the one depicted in the top left panel of Fig. 3. The Radon transform is the line integral of ρ in (E, t) along the straight line E = V 0 + α 0 t tan θ. The non-zero values in Fig. 9 are concentrated around θ = 0. A correlation between energy and time [39] that has a slope β in the time-energy plane would introduce a shift by arctan(β/α 0 ) along the θ axis [37]. Notice that the width of the non-zero regions along the V 0 axis at θ = 0 corresponds to the quantum-limited energy width /(2σ t ) of the Gaussian distribution, as δ → 0 in this example. with ∂T /∂E in (47) can be interpreted as the Wigner representation of an effective density matrix where ϕ E (t) = ϕ in (t)e −iEt/ is the incoming state ϕ in (t) shifted in energy by E, and Π(t) = e −iEt/ dT (E). In the limit of δ → 0, dT /dE becomes a delta-function, Π(t − t ) becomes 1, and ρ eff (E, t) in Eq. (48) reduces to the pure-state Wigner function (10). For finite δ though, the characteristic temporal width /δ of the Fourier transform Π(t − t ) in the time-domain representation [Eq. (50)] sets the upper limit on the coherence time of the wave-packet that can be resolved by the tomographic reconstruction. For δ > 0, electron partitioning at the barrier introduces shot noise that reduces the purity of the reconstructed effective mixed stateρ. We illustrate this by computing the convolution of the uncorrelated Gaussian wavepacket (35) and a Gaussian dT /dE from Eq. (37), Quantum purity, defined as γ = tr ρ 2 , is straightforward to compute in the Wigner representation, γ = h ρ eff (E, t) dE dt [37]. Using (51) we obtain the effective purity of tomographic reconstruction by a finite-width barrier for a quantum-limited Gaussian wave-packet, We see that γ in the case of an uncorrelated Gaussian wavepacket is equal to the ratio between the ideal Heisenberg uncertainty product ( /2) and the product of time and energy widths in the effective distribution (51) which is broadened by a finite resolution δ of the energy detector. Note that in general, the uncertainty product alone is insufficient to distinguish the incoherent broadening (e.g. δ-dependent terms in Eq. (51)] from coherent broadening (e.g. the energy-time correlation created in the source [37,39]). Hence a full tomography technique is essential for evaluating the quality a source with respect to the quantum limit of localizing a single particle in energy and time.
The dependence of δ on the parameters of the constriction depends on the details of the scattering interactionĤ bs and goes beyong the scope of this work. Nevertheless, we can use the simple sketch in Fig. 2 for an order of magnitude estimate. For a single-channel model with the backscattering amplitude distributed spatially over a length x ∈ [−x b . . . + x b ], the temporal width of τ (t) (and hence Π(t)) is limited by 2 x b /v. This gives an estimate for the achievable energy resolution of the tomography technique within a specific model of the scattering barrier depicted schematically in Fig. 2. This bound can also be interpreted as x b ∆p where ∆p = δ/v is the Fourier (diffraction) limit on the momentum resolution achievable over the characteristic length x b of the backscattering region for a particle propagating with velocity v.
One can expect from Eq. (53) that the effective length of the spatially uniform part of the scattering region is an important parameter for designing a barrier suitable for high-resolution tomography. On the other hand, as we have seen in Sec. III B, the length and the shape of the gate edge region from −x g to −x b is not important as long as the linear-in-time modulation condition (27) is fulfilled and the velocity dispersion can be ignored.

C. Measurement of the modified Wigner function
As mentioned in Sec. III, the modified Wigner functioñ ρ(E, t) is not a physical observable and does not represent the Wigner function of the wave packet at any stage. However, it can be experimentally obtained as the outcome of tomographic reconstruction. To see this, we consider a gate voltage V (t) that is the sum of a linear part V 0 + α t, as in the experiment of Ref. [37], and an additional perturbation V (t). This scenario is discussed theoretically at the end of Sec. III B, see Eqs. (33) and (34). Analogously to Eq. (47), the modified Wigner function corresponding to the perturbation V (t) only, can be expressed as where we have assumed ideal energy resolution, inserted V(t) = V 0 + α 0 t tan θ, and absorbed the correction to V (t) due to a possibly non-sharp edge of the gate into V 0 . This result means that, if we perform the inverse Radon transform of ∂Q(V 0 , θ)/∂V 0 , the outcome would be the modified Wigner distribution calculated for the non-linear part V (t) of the potential V (t) only. For example, the functionρ(E, t) of Eq. (41) can be measured using the linear tomographic reconstruction protocol of Ref. [37] while applying a gate voltage that is the sum of Eq. (38) and Eq. (40), i.e., By introducing such a sharp voltage "kick" on top of a linear time dependence of the gate potential, it would be possible to gather direct evidence of coherence of the electron source, since negative values in the measurement ofρ (E, t) arise from interference between phase-coherent parts of the incoming wave packet. This is relevant for the investigation of the properties of single-electron sources that emit solitary electrons high above the Fermi energy, where coherence is key to potential interferometric applications [3].
Our results also provide a way to describe a different kind of problem: sampling the potential V (t) with single electrons whose incoming distribution is well known [41]. In this case, the deviations in the measurement of ρ in (E, t) from the actual incoming distribution may allow one to measure the waveform V (t) beyond the classical resolution limit. In this context, Eq. (25) can elucidate the influence of the potential shape on the bandwidth of such a "quantum oscilloscopic" measurement.

VI. CONCLUSIONS
Inspired by recent experiments [36,37], we have constructed a fully quantum-mechanical description of a dynamical scattering problem of electron wave packets in a onedimensional chiral channel passing through a constriction subject to a time-dependent gate potential V (t). We have shown that the expression for the transmitted charge in this system is analogous to the corresponding classical expression, with a modified Wigner distribution functionρ being the quantum analog of the classical probability distribution. In particular, if the gate voltage time dependence is linear then the modified Wigner distribution is obtained from the Wigner function of the incoming electron by a time-dependent shift of the energy E. In this case the full quantum-mechanical theory agrees with the classical-limit estimate (1) justifying the quantum tomography proposal put forward in Ref. 37.
A finite energy resolution δ of the scattering barrier adds shot noise to the tomographically reconstructed Wigner distribution, limiting the maximal coherence time that can be probed by /δ. Small δ requires large spatial extent x b of the uniformly modulated part of the constriction, δ ≥ v/x b . We have shown that neither the length of the wave-packet nor the shape of the accelerating edge of the gate potential perturbs the linear tomography protocol as long as the modulation of the scattering region remains spatially uniform.
Our fully quantum-mechanical theory provides the tools to analyze experiments with a general time dependence of the gate voltage, relevant for wavepackets that are stretched in the time domain. The results presented here can be used to describe existing and propose new quantum measurement protocols. They can also serve to analyze a reverse problem: determining the gate potential V (t) when the incoming electron distribution is known. We have shown that additional modification of the Wigner function by non-linear gate modulation can in principle be revealed by quantum tomography within the same device. This opens a possibility to probe coherence of single-electron sources with explicit signatures of quantum interference (e.g., negativity of the modified Wigner function) which does not require multiple-path layouts.