# Quantum tomography of electrical currents

R. Bisognin, A. Marguerite, B. Roussel, M. Kumar, C. Cabart, Camille Chapdelaine, Mohammad-Djafari Ali, Jean-Marc Berroir, Erwann Bocquillon, Bernard Plaçais, et al.

► To cite this version:

R. Bisognin, A. Marguerite, B. Roussel, M. Kumar, C. Cabart, et al. Quantum tomography of electrical currents. Nature Communications, 2019, 10 (1), pp.3379. 10.1038/s41467-019-11369-5. hal-02318674

HAL Id: hal-02318674

https://hal.science/hal-02318674v1

Submitted on 6 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

https://doi.org/10.1038/s41467-019-11369-5

OPEN

# Quantum tomography of electrical currents

R. Bisognin<fnref n="1" />,⁶, A. Marguerite ID <fnref n="1" />,⁶, B. Roussel²,³, M. Kumar<fnref n="1" />, C. Cabart², C. Chapdelaine⁴,

A. Mohammad-Djafari⁴, J.-M. Berroir<fnref n="1" />, E. Bocquillon<fnref n="1" />, B. Plaçais<fnref n="1" />, A. Cavanna ID ⁵, U. Gennser⁵, Y. Jin⁵,

P. Degiovanni ID ² & G. Fève<fnref n="1" />

In quantum nanoelectronics, time-dependent electrical currents are built from few elementary excitations emitted with well-defined wavefunctions. However, despite the realization of sources generating quantized numbers of excitations, and despite the development of the theoretical framework of time-dependent quantum electronics, extracting electron and hole wavefunctions from electrical currents has so far remained out of reach, both at the theoretical and experimental levels. In this work, we demonstrate a quantum tomography protocol which extracts the generated electron and hole wavefunctions and their emission probabilities from any electrical current. It combines two-particle interferometry with signal processing. Using our technique, we extract the wavefunctions generated by trains of Lorentzian pulses carrying one or two electrons. By demonstrating the synthesis and complete characterization of electronic wavefunctions in conductors, this work offers perspectives for quantum information processing with electrical currents and for investigating basic quantum physics in many-body systems.

<footnote n="1">¹Laboratoire de Physique de l'École normale supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris-Diderot, Sorbonne Paris Cité, Paris 75005, France. ²Univ Lyon, Ens de Lyon, Université Claude Bernard Lyon 1, CNRS, Laboratoire de Physique, F-69342 Lyon, France. ³European Space Agency —Advanced Concepts Team, ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands. ⁴Laboratoire des signaux et systèmes, CNRS, Centrale-Supélec — Université Paris-Saclay, Gif-sur-Yvette F-91190, France. ⁵Centre de Nanosciences et de Nanotechnologies (C2N), CNRS, Univ. Paris-Sud, Université Paris-Saclay, 91120 Palaiseau, France. ⁶These authors contributed equally: R. Bisognin and A. Marguerite. Correspondence and requests for materials should be addressed to G.F. (email: gwendal.feve@ens.fr)</footnote>
In the field of quantum technologies, controlling elementary excitations such as single photons¹, single atoms², or single ions³ is a resource for encoding quantum information as well as a way to develop our understanding of basic quantum physics in complex many-body problems. In quantum electronics, the availability of on-demand single-electron sources⁴⁻⁷ offers the possibility to generate time-dependent electrical currents carrying a controlled number of electron and hole excitations of a degenerate electronic fluid. At low temperatures, phase coherence is preserved, such that describing these excitations in terms of well-defined wavefunctions is meaningful. In addition, by implementing electron sources in ballistic low-dimensional conductors⁴,⁷,⁸, the elementary excitations can be guided along one-dimensional channels and used as flying qubits⁹,¹⁰ carrying information encoded in their quantum state. However, despite the development of a rich experimental toolbox to generate and propagate electronic states in a controlled way, very few tools are currently available to characterize these states. In particular, measuring electron or hole wavefunctions embedded within a quantum electrical current has, so far, been out of reach.

This absence of a universal tomography protocol in the fermionic case may seem peculiar, considering that such protocols are now commonly implemented to reconstruct the state of bosonic fields¹¹,¹². However, there are important differences between bosonic and fermionic fields. First, bosonic tomography protocols involve the use of a classical field¹², which has no counterpart for fermions. Second, the vacuum of a fermionic system being a Fermi sea, the electron and hole excitations are thus defined by the addition and removal of a particle. Quantum-state reconstruction in the fermionic case can be illustrated by the sketch of Fig. 1. In a one-dimensional conductor, a T-periodic source generates a time-dependent current consisting of periodic pulses labeled by the index $l \in \mathbb{Z}$. To define unambiguously the electron and hole excitations, we take the conductor at chemical potential $\mu = 0$ and temperature $T_{el} = 0$ K as a reference. The electronic excitations correspond to the filling of the states above the Fermi sea (energy $\hbar\omega \ge 0$) and the hole excitations to the emptying of the states below the Fermi sea ($\hbar\omega \le 0$). We introduce in Fig. 1 the emitted time-translated electron (e) and hole (h) wavefunctions $\varphi_{l,i}^{(\alpha)}(t) = \varphi_i^{(\alpha)}(t - lT)$, and their emission probabilities $p_i^{(\alpha)}$, where $\alpha = e$ or h labels the electron or hole states and $i$ runs from 1 to $N_\alpha$, the total number of electron ($N_e$) and hole ($N_h$) wavefunctions emitted per period. These emitted electron and hole wavefunctions form a set of mutually orthogonal states: $\langle \varphi_{l'}^{(\alpha)} | \varphi_{l,i}^{(\alpha)} \rangle = \delta_{i,i'} \delta_{\alpha,\alpha'} \delta_{l,l'}$.

Here, by combining two-particle interferometry¹³ with signal processing¹⁴, we demonstrate a quantum current analyzer, which

extracts the emitted wavefunctions $\varphi_{l,i}^{(\alpha)}(t)$ and their emission probabilities $p_i^{(\alpha)}$ from any periodic electrical current. For benchmarking, we first apply our analyzer on sinusoidal currents to validate our extraction method in the general case, when several excitations are emitted with non-unit probability. Sinusoidal drives are well suited to test the robustness of our procedure by comparing our results with parameter-free theoretical predictions. We then apply our technique to trains of Lorentzian pulses carrying an integer charge $q = -e$ and $q = -2e$ and extract their full content in terms of single-electron wavefunctions. At zero temperature, Lorentzian pulses of integer charge $-eN_e$ are predicted to generate an integer number $N_e$ of excitations exclusively above the Fermi sea⁶,¹⁵⁻¹⁷. By extracting all the emitted wavefunctions, we observe that thermal effects lead to the generation of a statistical mixture between the expected zero-temperature wavefunctions and additional undesired states. From the measurement of the emission probability of each generated wavefunction, we provide a quantitative analysis of the purity of the generated electronic states. By identifying specific single-electron and hole wavefunctions and determining their emission probabilities for various types of time-dependent currents, our work opens the way to a precise and systematic characterization of quantum information carried by electrical currents.

## Results

**Electronic coherence and Wigner distribution.** The main difficulty behind the extraction of the electron and hole wavefunctions from an electrical current lies in the explicit connection between the wavefunctions and a measurable physical quantity. So far, most of the characterizations of the excitations generated by electronic sources have been limited to the measurements of the average electrical current $I(t)^{4,18}$ and electronic distribution function $f(\omega)^{7,19}$. They provide information on the time and energy distributions, but cannot access the phase of electronic wavefunctions, which requires the use of interferometry techniques.

In analogy with optics, all interference effects are encoded in the first-order electronic coherence $G_{\rho,x}^{(e)}(t, t')^{20,21}$ defined as the time correlations of the fermion field $\hat{\Psi}(x, t)$, which annihilates an electron at position $x$ and time $t$ of the one-dimensional conductor: $G_{\rho,x}^{(e)}(t, t') = \langle \hat{\Psi}^{\dagger}(x, t') \hat{\Psi}(x, t) \rangle_{\rho}$. To simplify the notations in the rest of the paper, we suppress the superscript (e) and the dependence on the position $x$ and on the many-body density operator $\rho$ in the expression of the electronic coherence, which is written as $G(t, t')$. More generally, $G(t, t')$ contains all the information on the single-particle properties of the many-body electronic state. Electronic coherence being a priori a complex function, it is more convenient to use the electronic Wigner distribution$^{22,23}$ $W(t, \omega)$ obtained from $G(t + \tau/2, t - \tau/2)$ by Fourier transform along the time difference $\tau$. $W(t, \omega)$ is a real function of marginal distributions $I(t)$ and $f(\omega)$ obtained by respectively integrating $W(t, \omega)$ over energy $\omega$ and time $t$, thereby demonstrating that they only provide partial information.

Subtracting the reference contribution characterized by the zero-temperature Fermi distribution $\Theta(-\omega)$ (where $\Theta$ is the Heaviside function) defines $\Delta_0 W(t, \omega) = W(t, \omega) - \Theta(-\omega)$ (or equivalently its Fourier transform $\Delta_0 G(t, t')$). $\Delta_0 W(t, \omega)$ and $\Delta_0 G(t, t')$ are the key quantities that we explicitly connect to the wavefunctions $\varphi_{l,i}^{(\alpha)}$ and emission probabilities $p_i^{(\alpha)}$. This connection is trivial in the pure-state single-body case, that is, when a single excitation (either an electron or a hole) of wavefunction $\varphi$ is emitted with unit probability. In this simple limit, $\Delta_0 W(t, \omega) = W_{\varphi}(t, \omega)$, where $W_{\varphi}(t, \omega)$ is the Wigner representation$^{24}$ of the

<figure><img src="image_2.png" /><figcaption>**Fig. 1** Sketch of the elementary excitations generated by a periodic current $l$ ($t$). Electron and hole states are generated at each period above and below the Fermi sea. The time-translated wavefunctions $\varphi_{l,i}^{(\alpha)}(t) = \varphi_i^{(\alpha)}(t - lT)$, for $\alpha = e$ (electron) or h (hole) and $1 \le i \le N_\alpha$ are emitted at each period with probability $p_i^{(\alpha)}$</figcaption></figure>
wavefunction $\varphi$:

$$W_{\varphi}(t, \omega) = \int d\tau \varphi\left(t + \frac{\tau}{2}\right) \varphi^{*}\left(t - \frac{\tau}{2}\right) e^{i\omega\tau}. \quad (1)$$

In a recent experiment²⁵, Jullien et al. performed the first reconstruction of $W(t, \omega)$ in the case of a periodic train of single-electron Lorentzian pulses. Assuming that the single-body limit was valid, they extracted the electronic wavefunction $\varphi$ using $\Delta_0 W(t, \omega) = W_{\varphi}(t, \omega)$. However, the single-body limit can never be completely achieved, due to the presence of thermal excitations, to the periodic emission from the source or to deformations of the current pulse associated with imperfections of the voltage drive, or due to more fundamental effects, such as the Coulomb interaction. In addition, for multi-electron states, such as Lorentzian pulses carrying an integer number of excitations $N_e > 1^{16,26,27}$, going beyond the single-body limit to extract the electronic wavefunctions is absolutely required.

In the more complex case where one has to consider several wavefunctions generated with arbitrary probabilities, only specific sets of drives²⁸⁻³⁰ have been theoretically investigated. Furthermore, the connection between experimentally accessible quantities and the emitted wavefunctions was missing. Following the work of ref. ¹⁴, we explicitly connect the emitted electron $\varphi_{l,i}^{(e)}$ and hole $\varphi_{l,i}^{(h)}$ wavefunctions to the electronic coherence $\Delta_0 G$ by diagonalizing $\Delta_0 G$ in the subspace of electron and hole states (see "Methods"). As a result of the diagonalization procedure, $\Delta_0 G$ can be decomposed on the basis of electron and hole states $\varphi_{l,i}^{(\alpha)}$ by introducing the matrix elements

$$g_{ij}^{(\alpha\beta)}(l) = \langle \varphi_{l,i}^{(\alpha)} | \Delta_0 G | \varphi_{0,j}^{(\beta)} \rangle = \int dt \, dt' \varphi_{l,i}^{(\alpha)}(t)^* \Delta_0 G(t, t') \varphi_{0,j}^{(\beta)}(t').$$

$$\Delta_0 G(t, t') = \sum_{i=1}^{N_e} \sum_{(l, l') \in \mathbb{Z}^2} g_i^{(ee)}(l-l') \varphi_{l,i}^{(e)}(t) \varphi_{l',i}^{(e)}(t')^* \\ - \sum_{i=1}^{N_h} \sum_{(l, l') \in \mathbb{Z}^2} g_i^{(hh)}(l-l') \varphi_{l,i}^{(h)}(t) \varphi_{l',i}^{(h)}(t')^* \\ + \sum_{i,j} \sum_{(l, l') \in \mathbb{Z}^2} \left( g_{ij}^{(eh)}(l-l') \varphi_{l,i}^{(e)}(t) \varphi_{l',j}^{(h)}(t')^* \right. \\ \left. + g_{ji}^{(he)}(l-l') \varphi_{l,j}^{(h)}(t) \varphi_{l',i}^{(e)}(t')^* \right). \quad (2)$$

As the electron and hole wavefunctions $\varphi_i^{(\alpha)}$ are extracted from the diagonalization of $\Delta_0 G$, it naturally implies that there are no quantum coherences in Eq. (2) between states $\varphi_{l,i}^{(\alpha)}$ and $\varphi_{l',i}^{(\alpha)}$ whenever $i \neq i'$: $g_{i \neq i'}^{(\alpha\alpha)} = 0$.

Each term of Eq. (2) can be separately interpreted. The first (second) term represents the contribution of electron (hole) wavepackets to the first-order coherence. For $l=l'$, the real numbers $0 \le g_i^{(ee)}(0) \le 1$ and $0 \le g_i^{(hh)}(0) \le 1$ represent the probability for emitting the electron ($\varphi_i^{(e)}$) and hole ($\varphi_i^{(h)}$) wavefunctions. Following the notation introduced at the beginning of the paper, we thus have $p_i^{(\alpha)} = g_i^{(\alpha\alpha)}(0)$. Compared with the simple picture sketched in Fig. 1, the T-periodicity of the source requires to consider also the complex numbers $g_i^{(ee)}(l-l')$ (resp. $g_i^{(hh)}(l-l')$) for $l \neq l'$ representing coherences between electronic (resp. hole) wavepackets emitted at different periods. The last two terms of Eq. (2) then represent the coherence between the electron and hole states $\varphi_{l,i}^{(e)}$ and $\varphi_{l',j}^{(h)}$ encoded in $g_{ij}^{(eh)}(l-l')$. It can only be nonzero when the electron and hole emission probabilities $p_i^{(e)}$ and $p_j^{(h)}$ are different from 0 or 1. It

then expresses the existence of a quantum superposition between the unperturbed ground state and the creation of the electron/hole pair built from the single-particle states $\varphi_i^{(e)}$ and $\varphi_j^{(h)}$. The coefficients $g_{ij}^{(eh)}(l-l')$ encode the modulus and phase of such a quantum superposition. In this description, the ideal emission of a quantized number of $N_e$ electrons and $N_h$ holes is characterized by $g_i^{(ee)}(l-l') = \delta_{l,l'}$ and $g_i^{(hh)}(l-l') = \delta_{l,l'}$ implying that $g_{ij}^{(eh)}(l-l') = 0$.

This formalism serves as the theoretical background for the extraction of the electron and hole wavefunctions from experimental measurements. Using two-particle interferences, we proceed to the measurement of the electronic coherence $\Delta_0 W$ and $\Delta_0 G$ for arbitrary electrical currents. We then implement an algorithm (see "Methods"), which identifies the emitted wavefunctions $\varphi_i^{(e)}$ and $\varphi_j^{(h)}$ from the diagonalization of $\Delta_0 G$ in the subspace of electron and hole states and recasts it in the form given by Eq. (2). This set of data describes completely the single-particle content of the electronic current and quantifies how far it deviates from the ideal emission regime.

**Experimental setup and protocol.** The experiment is performed in a high-mobility GaAs/AlGaAs two-dimensional electron gas placed in a strong perpendicular magnetic field so as to reach the quantum Hall regime at filling factors $v=2$ or $v=3$, where charge propagates along one-dimensional chiral-edge channels. We focus on the propagation on the outer-edge channel, which realizes a ballistic spin-polarized one-dimensional conductor. The electronic source is a metallic gate capacitively coupled to the edge channels, allowing us to shape any charge distribution³¹ by applying the proper time-dependent voltage to the gate. The resulting Wigner distribution $W_S(t, \omega)$ can be reconstructed³² by measuring two-electron interferences³³,³⁴, using an electronic Hong-Ou-Mandel³⁵ interferometer³⁶⁻³⁸. As shown in Fig. 2, the interferometer consists of a quantum point contact used as an electronic beam splitter partitioning the excitations propagating from inputs 1 and 2 with transmission probability $T$. Input 1 is connected to the source, whereas input 2 is connected to a voltage-driven ohmic contact that will generate a set of known reference states, called probe states, of Wigner distribution $W_{P_n}$ for $n \in \mathbb{N}$. For each probe state, we measure the excess noise $\Delta S_n$ at output 3 between the source being switched on and off¹³:

$$\Delta S_n = 2e^2 T(1-T) \int \frac{d\omega}{2\pi} \left[ \frac{\overline{\Delta W_S(t, \omega)^t}}{2} (1 - 2f_{\text{eq}}(\omega)) - \frac{\overline{2\Delta W_S(t, \omega)\Delta W_{P_n}(t, \omega)^t}}{\overline{\Delta W_{P_n}(t, \omega)^t}} \right] \quad (3)$$

where $\overline{\cdot \cdot \cdot t}$ denotes the average over time $t$, and $\Delta W_{S/P_n}$ are, respectively, the source and probe excess Wigner distributions with respect to the Fermi-Dirac distribution $f_{\text{eq}}(\omega)$ at temperature $T_{\text{el}} \neq 0$: $W_{S/P_n}(t, \omega) = f_{\text{eq}}(\omega) + \Delta W_{S/P_n}(t, \omega)$. The first term in Eq. (3) represents the classical random partition noise of the source. It is reduced by the second term in Eq. (3), which represents the antibunching between indistinguishable source and probe excitations colliding on the splitter. Their degree of indistinguishability is given by the overlap between $\Delta W_S$ and $\Delta W_{P_n}$. By properly choosing the set of probe states, Eq. (3) allows for the reconstruction of any unknown Wigner distribution²²,³².

A convenient set of probe states can be used to reconstruct each harmonic of the Fourier expansion of the excess source Wigner distribution

$$\Delta W_S(t, \omega) = \sum_{n \in \mathbb{Z}} \Delta W_{S,n}(\omega) e^{2\pi i n \int t}, \quad (4)$$
where $f = 1/T$ denotes the driving frequency. The $n=0$ harmonic represents the source excess electronic distribution function $\Delta f(\omega)$. All the time dependence of $\Delta W_S(t, \omega)$ is encoded in the $n \neq 0$ harmonics. To select the contribution from the $n$th harmonic in Eq. (3), we apply on the probe input a small ac signal at frequency $nf$ on top of a dc bias$^{32}$: $V_{P_n}(t) = V_{dc} + V_{P_n} \cos(2\pi n f t + \varphi)$. The resulting Wigner distribution $W_{P_n}$ (plotted in Fig. 2) evolves periodically in time at frequency $nf$. By measuring the output noise $\Delta W_{S,n}$ as a function of $\varphi$ and $V_{dc}$ (see “Methods”), the real and imaginary parts of $\Delta W_{S,n}(\omega)$ can be extracted.

**Electronic Wigner distribution of sinusoidal drives.** We first apply our quantum current analyzer to sinusoidal drives, $V_S(t) = V_S \cos(2\pi n f t)$ at various frequencies $f$. Figure 3a presents the measurements of the $n=0, 1, 2, 3$ harmonics of $\Re(\Delta W_{S,n})$ ($\Im(\Delta W_{S,n})=0$) of three sinusoidal drives of similar amplitudes ($V_S \approx 32$ µV). We first focus on the effect of frequency by comparing $\Delta W_{S,n}$ for $f=10$ MHz and $f=9$ GHz at $T_{el}=100$ mK. The $n=0, 2$, and $3$ harmonics are lower for $f=9$ GHz compared with $f=10$ MHz ($\Delta W_{S,n}=3$ even falls below our experimental resolution for $f=9$ GHz).

Indeed, in a photo-assisted description of electronic transport$^{39,40}$, the $n \neq 1$ harmonics are related to multiphoton absorption/emission processes, whose strength increases with the ratio $eV_S/hf$, which equals 800 for $f=10$ MHz compared with 0.8 for $f=9$ GHz. We then turn to the effect of temperature by comparing $\Delta W_{S,n}$ for the two drives at $f=9$ GHz, but at different temperatures. Decreasing the temperature from $T_{el}=100$ mK to $T_{el}=60$ mK leads to a narrowing of all the harmonics and to an increase of their amplitude. For the three drives, the agreement between the data and theoretical predictions (dashed lines) is excellent, showing the robustness of our reconstruction procedure. After measuring all relevant $\Delta W_{S,n}$ we can combine them in Eq. (4) to reconstruct $W_S(t, \omega)$.

The Wigner distributions are represented in Fig. 3b. Within experimental accuracy, the $f=10$-MHz case follows an equilibrium distribution function, $W_S(t, \omega) = f_{eq,\mu}(t)(\omega)$, with a time-varying

chemical potential following the ac drive: $\mu(t) = -eV_S \cos(2\pi n f t)$. This is expected as the $f=10$ MHz case corresponds to a quasi-classical current ($hf \ll k_B T_{el}$) characterized by bounded values of the Wigner distribution, $0 \le W(t, \omega) \le 1$, such that $W(t, \omega)$ can be interpreted as a time-dependent electronic distribution function and viewed as an adiabatic evolution of the stationary (dc) case$^{22}$. In contrast, $hf \ge k_B T_{el}$ corresponds to the quantum case, where the Wigner distribution can take negative or above one values. This is what we observe for the $f=9$-GHz drives, with a strong emphasis of these quantum features at the lowest temperature $T_{el}=60$ mK. Consequently, in the quantum regime, single-particle properties are no longer described in terms of a time-varying electronic distribution function. This is the case where $W(t, \omega)$ can be used to extract electron and hole wavefunctions.

**Electron/hole wavefunctions generated by sinusoidal drives.**

The second step of our analyzer extracts individual electronic wavepackets from the reconstructed Wigner distribution by implementing an algorithm (see "Methods"), which recasts our measurements in the form of Eq. (2). Figure 4 presents the result of this analysis on the experimental data obtained for the $f=9$-GHz sinusoidal drives. As the probability to emit more than one electron/hole is very small, the analysis can be limited to one electron $\varphi_1^{(e)}$ and one hole $\varphi_1^{(h)}$ wavefunction ($p_{i \to 1}^{(a)} \approx 10^{-3} \ll 1$). They are plotted in the Wigner representation in the case $T_{el}=60$ mK in Fig. 4a. The hole is shifted by half a period with respect to the electron and its energy distribution $|\varphi_1^{(h)}(\omega)|^2$ mirrors that of the electron's at positive energy. As a figure of merit of the procedure, we evaluate the state fidelity defined as the overlap between electron and hole wavefunctions extracted from the experimental data and the electron and hole wavefunctions extracted from numerical computations of the Wigner distribution using Floquet scattering theory (see Supplementary Note 1). The results are in excellent agreement with a fidelity >0.99 for all the extracted wavefunctions, demonstrating the accuracy of the state reconstruction.

<figure><img src="image_3.png" /><figcaption>**Fig. 2** Sketch of the experimental setup. The two-dimensional electron gas is represented in blue color and the edge channels as blue lines. A quantum point contact (red color) is used as an electronic beam splitter. A time-dependent voltage $V_S(t)$ is applied to a mesoscopic capacitor (gate in gold color capacitively coupled to the edge channel) placed at input 1 of the beam splitter and generates the unknown Wigner distribution $\Delta W_S(t, \omega)$. The probe signal $V_{P_n}(t)$, a low-amplitude sinusoidal drive at frequency $nf$ is generated at input 2. The corresponding probe Wigner distributions $\Delta W_{P_n}(t, \omega)$ are plotted for $n=0$ to $n=2$ (the frequency is $f=5$ GHz and the temperature $T_{el}=80$ mK). The current noise at the output of the splitter is converted to a voltage noise on the quantized resistance $R_\nu = h/(\nu e^2)$. $R_\nu$ is connected to an LC tank circuit used to shift the measurement frequency at the resonance $f_0 = 1.45$ MHz</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 3** Wigner distribution of sinusoidal drives. **a** Measured $\Delta W_{S,n}(\omega)$ for $n=0$ to $n=3$ for sinusoidal drives at frequency $f=10$ MHz and $f=9$ GHz. The observed parity: $\Delta W_{S,n}(-\omega) = (-1)^{n+1}\Delta W_{S,n}(\omega)$, directly stems from the electron/hole symmetry of the sinusoidal drive. Error bars are defined as standard error of the mean. **b** Time-energy representation of the Wigner distribution $W_S(t, \omega)$</figcaption></figure>

Figure 4b depicts the moduli of the inter-period coherences between emitted wavepackets for different temperatures (bars represent numerical simulations, circles represent data). When the temperature increases, the occupation probabilities which, in the present case, are very close due to electron/hole symmetry ($p_1^{(e)} \approx p_1^{(h)}$) increase from 0.17 (numerical calculation at $T_{el}=0$ K) to 0.25 (60 mK) and 0.27 (100 mK). Coherences between different periods ($g_1^{(ee)}(l \neq 0)$ and $g_1^{(hh)}(l \neq 0)$) also appear and extend over the thermal coherence time ($h/k_B T_{el} \approx 0.5$ ns at 100 mK). This reflects that at finite temperature, the electron and hole states $\phi_1^{(e/h)}$ have a finite probability to be occupied by thermal excitations. As the probability to emit the electron and hole differ from 1, we also observe nonzero electron/hole coherence: $g_{11}^{(eh)}(l-l') \neq 0$. Interestingly, these terms are suppressed by thermal fluctuations, reflecting the transition from a pure quantum state at $T_{el}=0$ K to a statistical mixture at higher temperature. At $T_{el}=0$ K, a single process occurs: the generation of the quantum superposition between the unperturbed ground state (with probability $1-p_1^{(e)}$) and the creation of the electron/hole pair (with probability $p_1^{(e)}$). Thermal fluctuations allow two additional processes: only the electron state, or only the hole state, can be generated. The resulting state at finite temperature is a

statistical mixture between these three processes. Our algorithm enables a quantitative description by computing a purity indicator, $P$, from the extracted inter-period coherences (see “Methods”). It quantifies the weight of coherent electron/hole processes with respect to all emitted excitations. By construction, $P=1$ at zero temperature, and from our experimental data, decreases to 0.71 at $T_{el}=60$ mK and to 0.58 at $T_{el}=100$ mK. Numerical evaluation of the same quantity calculated using Floquet scattering theory (see Supplementary Note 1) give 0.999 at zero temperature, 0.725 at 60 mK and 0.588 at 100 mK in very good agreement with the experimental data.

**Single-electron Lorentzian pulse ($q = -e$).** We now turn to the analysis of a current generated by periodic Lorentzian voltage pulses $V_S(t) = \sum_l -\frac{V_0}{1+(t-IT)^2/\tau^2}$ with $\tau = 42$ ps, $f=4$ GHz, and $V_0$ chosen such that each pulse carries exactly a single-electron charge: $\frac{e^2}{\hbar} \int_0^T V_S(u) du = -e$. The Lorentzian pulses are generated by calibrating the amplitude and phase of each harmonic of the current at the location of the beam splitter (see “Methods”). This ensures that phase and amplitude shifts caused by Coulomb interaction effects$^{26,41}$ during the propagation along the edge channels can be absorbed in the calibration process. This procedure allows us to demonstrate the proof of principle of our
<figure><img src="image_5.png" /><figcaption>**Fig. 4** Electron and hole wavefunctions generated by low-amplitude sinusoidal drives. **a** Wigner distribution functions $W_{f^{(c/h)}}(t, \omega) = \int d\tau \varphi_1^{(c/h)}(t + \tau/2) \varphi_1^{(c/h)*}(t - \tau/2) e^{i\omega\tau}$ for the dominant electronic $\varphi_1^{(e)}$ and hole $\varphi_1^{(h)}$ wavefunctions for $f = 9$ GHz and $T_{el} = 60$ mK ($\varphi_1^{(e/h)}$ obtained at $T_{el} = 100$ mK are almost identical). The panels in the margins of the color plots represent the time $|\varphi_1^{(c/h)}(t)|^2/f$ and energy $f|\varphi_1^{(c/h)}(\omega)|^2$ distributions obtained by integrating $W_{f^{(c/h)}}(t, \omega)$ over $\omega$ and t. **b** Moduli of the interperiod coherence $|g^{(ee)}(l)|$, $|g^{(hh)}(l)|$ and $|g^{(eh)}(l)|$ (colored bars correspond to numerical calculations and colored dots to the experimental data)</figcaption></figure>

quantum current tomography without having to consider Coulomb interaction effects.

The measured $n=0$ to $n=4$ harmonics of $\Delta W_S(t, \omega)$ are plotted in Fig. 5a. They are located on the positive energy side with maxima shifted by $nhf/2e$ for increasing $n$ and, indeed, take very small value for $|\omega| \le n\pi f$ showing that almost no hole excitation is emitted. Their energy width is imposed by the temporal width of the pulse $\tau$ and their amplitude reproduces the decrease of the harmonics of $I(t)$ (obtained by integration of $\Delta W_S$, $n(\omega)$). The overall agreement with theoretical predictions (dashed lines) is good (no fitting parameter). The resulting Wigner distribution $W_S(t, \omega)$ is plotted in Fig. 5b. Strong nonclassical features ($W_S(t, \omega) \approx 1.2$) are observed at the location of the electron excitation in the $(t, \omega)$ plane.

The wavefunctions extracted from our analysis are plotted in Fig. 5c. First, as expected for a single-electron Lorentzian pulse, this analysis quantitatively confirms that almost no hole excitation is emitted: $p^{(h)} = 0.03 \pm 0.01$. Second, contrary to the previous case (low-amplitude sine drive), two electronic wavefunctions $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$ contribute with emission probabilities $p_1^{(e)} = 0.83 \pm 0.01$ and $p_2^{(e)} = 0.18 \pm 0.01$. This means that finite temperature (depending on the ratio $k_B T_{el}/hf$, see Supplementary Note 2) leads to the generation of a statistical mixture$^{42}$ (of purity $P = 0.75$) between $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$ with weights $p_1^{(e)}$ and $p_2^{(e)}$ ($p_1^{(e)} + p_2^{(e)} = 1.01 \pm 0.01$). Note that such an emission of a statistical mixture between different single-electron wavefunctions could not be captured within the single-body limit considered in the previous analysis of Jullien et al.$^{25}$, which

underlines the need for developing the general approach we demonstrate here.

The wavefunctions $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$ can be compared with the expected ones generated by Lorentzian voltage pulses at zero temperature. A single Lorentzian pulse of temporal width $\tau$ and carrying an integer number $N_e$ of electrons is the Slater determinant built over the $1 \le n \le N_e$ electronic wavefunctions$^{16,26}$ $\varphi_n^{(\text{single})}(\omega) = \frac{1}{\sqrt{N}} \Theta(\omega) \exp(-\omega\tau) \mathcal{L}_{n-1}(2\omega\tau)$, where $N$ is a normalization constant, $\mathcal{L}_n$ is the $n$th Laguerre polynomial. The two first ones $\varphi_{n-1}^{(\text{single})}(\omega)$ and $\varphi_{n-2}^{(\text{single})}(\omega)$, plotted in red dashed lines in Fig. 5c, are very similar to $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$. However, they do not reproduce the discretization of the energy distribution in units of $hf$. This discretization is related to the periodicity of the single-electron emission, which is not captured by the expression of $\varphi_n^{(\text{single})}$. The $n=1$ wavefunction of the periodic train of Lorentzian pulses has been shown$^{14,43}$ to be given by $\varphi_{L,n=1}(\omega) = \frac{1}{\sqrt{N}} \Theta(\omega) \exp(-\omega_1 \tau)$ where $\omega_1 = 2\pi f \left[ \frac{\omega}{2\pi f} \right]$ has quantized steps $2\pi f$ related to the pulse periodicity. Here, we generalize this expression to the $n=2$ and $n=3$ wavefunctions by using the following ansatz: $\varphi_{L,n}(\omega) = \frac{1}{\sqrt{N}} \Theta(\omega) \exp(-\omega_n \tau) \mathcal{L}_{n-1}(2\omega_n \tau)$ where $\omega_n = 2\pi f \left( x_n + \left[ \frac{\omega}{2\pi f} \right] \right)$. We take $x_1 = 0$ following refs.$^{14,43}$ and then numerically deduce $x_2 = 0.33$ and $x_3 = 0.24$ from the constraint that the wavefunctions should be orthogonal: $\langle \varphi_{L,n} | \varphi_{L,n'} \rangle = \delta_{n,n'}$. Comparing the wavefunctions $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$ extracted from our experimental data to these theoretical predictions, we observe that $\varphi_1^{(e)}$ is very close
<figure><img src="image_6.png" /><figcaption>**Fig. 5** Electron wavefunction generated by a Lorentzian pulse, $q = -e$. **a** Measured $\Delta W_{S,n}(\omega)$ ($n=0$ to $n=4$) for the single-electron Lorentzian pulse at $f=4$ GHz, $\tau=42$ ps, and $T_{el}=50$ mK. Error bars are defined as standard error of the mean. **b** $W_S(t, \omega)$, the dashed line represents the voltage pulse $V(t)=h/e^2/(t)$, where $l(t)$ is obtained by integrating $W_S(t, \omega)$ on energy $\omega$. **c** Wigner representation of $\varphi_1^{(e)}(t)$ (left) and $\varphi_2^{(e)}(t)$ (right). The panels in the margins represent the time $|\varphi_1^{(e)}(t)|^2$ and energy $|\varphi_2^{(e)}(\omega)|^2$ distributions obtained by integrating $W_{\varphi^{(e)}}(t, \omega)$ over $\omega$ and $t$. $|\varphi_1^{(e)}(\omega)|^2$ is represented in log scale for better comparison with theoretical predictions with single shot $\varphi_n^{(\text{single})}$ (red dashed line) and periodic $\varphi_{L,n}$ (blue dashed line) Lorentzian wavepackets (with $\tau=42$ ps). As it can be seen from the very good agreement with the blue dashed line, the steplike behavior of $|\varphi_1^{(e)}(\omega)|^2$ comes from the periodicity of the drive</figcaption></figure>

to the expected wavefunction $\varphi_{L,n=1}$ (blue dashed line in Fig. 5c) with an overlap of 0.98. Interestingly, $\varphi_2^{(e)}$ which is emitted at higher energy strongly resembles $\varphi_{L,n=2}$ with an overlap of 0.93 leading to this simple interpretation of temperature effects: finite temperature leads to the emission of a statistical mixture between the expected $n=1$ Lorentzian wavefunction and the wavefunctions corresponding to higher excitations numbers $n>1$. Importantly, our observations are not related to imperfections of the emission drive or to errors of our extraction method, but only to thermal effects. This probabilistic description of the electron state stems from the nonzero entropy of the finite-temperature ground state, which reveals the statistical (non-quantum) fluctuations of the ground state. This interpretation is confirmed by numerical calculations. Applying our method on perfect periodic Lorentzian pulses calculated using Floquet scattering theory (see Supplementary Note 2), we recover that for $T_{el}=0$ K, $p_1^{(e)}=1$ and $p_2^{(e)}=0$ (pure state) but for $T_{el}=50$ mK, $p_1^{(e)}=0.84$ and $p_2^{(e)}=0.18$, in very good agreement with our experimental results.

**Two-electron Lorentzian pulse** ($q = -2e$). Finally, we analyze periodic trains of Lorentzian pulses carrying the charge of two electronic excitations: $V_S(t) = \sum_l -\frac{2V_0}{1+(t-IT)^2/\tau^2}$. Our

experimental results are plotted in Fig. 6. As in the single-electron case, thermal effects lead to the generation of a statistical mixture of different wavefunctions, one more than the number of emitted charges. The first wavefunction $\varphi_1^{(e)}$ is emitted with unit probability $p_1^{(e)}=1$, but $\varphi_2^{(e)}$ and $\varphi_3^{(e)}$ are emitted with probabilities smaller than one, $p_2^{(e)} = 0.69 \pm 0.02$ and $p_3^{(e)} = 0.24 \pm 0.02$, reflecting that the emitted state is a statistical mixture of purity $P = 0.68$.

Interestingly, and contrary to the $q = -e$, case, $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$ do not correspond to the expected Lorentzian wavefunctions $\varphi_{L,n=1}$ and $\varphi_{L,n=2}$ plotted in red dashed line in Fig. 6c. This can be understood by discussing first the zero-temperature case. At $T_{el} = 0$ K, the generated state is predicted to be described by the Slater determinant formed from the wavefunctions $\varphi_{L,n=1}$ and $\varphi_{L,n=2}$. However, any choice of basis obtained by a linear combination of $\varphi_{L,n=1}$ and $\varphi_{L,n=2}$ is equally valid to describe this Slater determinant (see Supplementary Note 3). At finite temperature, $T_{el} \neq 0$ K, this ambiguity in the choice of the two wavefunctions describing the electronic state is lifted as $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$ are no longer generated with the same probability: $p_1^{(e)} \neq p_2^{(e)}(\neq 1)$.

Searching for the basis of states which maximizes the overlap with $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$, we observe that finite temperature favors the
<figure><img src="image_7.png" /><figcaption>**Fig. 6** Electron wavefunction generated by a Lorentzian pulse, $q = -2e$. **a** Measured $\Delta W_{S,n}(\omega)$ ($n=0$ to $n=4$) for the two-electron Lorentzian pulse at $f = 4$ GHz, $\tau = 42$ ps, and $T_{el} = 50$ mK. Error bars are defined as standard error of the mean. **b** $W_S(t, \omega)$, the dashed line represents the voltage pulse $V(t) = h/e^2/(t)$, where $I(t)$ is obtained by integrating $W_S(t, \omega)$ on energy $\omega$. **c** Wigner representation of $\varphi_1^{(e)}(t)$ (left) and $\varphi_2^{(e)}(t)$ (middle) and $\varphi_3^{(e)}(t)$ (right). The panels in the margins represent the time $|\varphi_i^{(e)}(t)|^2$ and energy $|\varphi_i^{(e)}(\omega)|^2$ distributions obtained by integration of $W_{\varphi_i^{(e)}}(t, \omega)$ over $\omega$ and $t$. $|\varphi_1^{(e)}(\omega)|^2$ is represented in log scale. The red dashed line represent the theoretical predictions for the $n=1$ to $n=3$ wavefunctions of periodic trains of Lorentzian pulses $\varphi_{L,n}$. The blue dashed lines represent the theoretical predictions for $\varphi'_{L,n}$ obtained from linear combinations of the $\varphi_{L,n}$.</figcaption></figure>

emergence of two specific states obtained by the following linear combination of $\varphi_{L,n=1}$ and $\varphi_{L,n=2}$:

$$ \varphi'_{L,n=1} = \cos\left(\frac{\theta}{2}\right) \varphi_{L,n=1} + \sin\left(\frac{\theta}{2}\right) \varphi_{L,n=2} \quad (5) $$

$$ \varphi'_{L,n=2} = \sin\left(\frac{\theta}{2}\right) \varphi_{L,n=1} - \cos\left(\frac{\theta}{2}\right) \varphi_{L,n=2} \quad (6) $$

with $\theta \sim \frac{\pi}{2} \times 0.37$. The single-electron wavefunctions $\varphi'_{L,n=1}$ and $\varphi'_{L,n=2}$, plotted in blue dashed lines in Fig. 6c, have a very strong overlap with $\varphi_1^{(e)}$ and $\varphi_2^{(e)}$ (0.99 and 0.96). As $p_1^{(e)} = 1$ and $p_2^{(e)} = 0.69$, it means that with probability 0.69 the two electron state described by the Slater determinant formed from $\varphi'_{L,1}$ and $\varphi'_{L,2}$ is generated. This state is equivalent to the expected Slater determinant formed from $\varphi_{L,1}$ and $\varphi_{L,2}$. However, with probability $p_2^{(e)} = 0.24$ a different two electron state is generated corresponding to the Slater determinant formed from $\varphi'_{L,1}$ and $\varphi_3^{(e)}$.

The connection between $\varphi_3^{(e)}$ and a theoretically predicted wavefunction is less clear. $|\varphi_3^{(e)}(\omega)|^2$ resembles $|\varphi_{L,n=3}(\omega)|^2$, the energy distribution of the $n=3$ Lorentzian wavefunction (blue dashed lines in Fig. 6c). However the time distributions are

different, which is reflected by the relatively small overlap of 0.86 between the two states.

As for the $q = -e$ case, we can check using numerical computations that the generation of a statistical mixture between two different Slater determinants is caused by the finite temperature. The simulations presented in Supplementary Note 3 confirm that the probability to generate the Slater determinant formed from $\varphi_{L,1}$ and $\varphi_{L,2}$ decreases from 1 at zero temperature to 0.79 at $T_{el} = 50$ mK, while the probability to generate the Slater determinant formed from $\varphi'_{L,1}$ and $\varphi_3^{(e)}$ increases from 0 at zero temperature to 0.18 at $T_{el} = 50$ mK in reasonable agreement with our observations.

## Discussion

We have demonstrated a quantum tomography protocol for arbitrary electrical currents. Without any a priori knowledge on the electronic state, this protocol extracts all the electron and hole wavefunctions and their emission probabilities.

Our protocol is the tool of choice for characterizing single- to few-electron sources by extracting all the $N_e$ single-electron wavefunctions generated at each period. We have analyzed in this work the $N_e = 1$ and $N_e = 2$ Lorentzian voltage pulses and have extracted the wavefunctions of each generated electronic
excitation. The generation of Lorentzian electronic wavepacket in a ballistic one-dimensional channel, which we demonstrate here, is an important milestone for the development of time-resolved quantum electronics. Numerous recent proposals suggest the use of time-resolved charge or energy currents⁴⁴ carried by integer charge Lorentzian voltage pulses to probe the timescales of quantum coherent conductors⁴⁵ or to dynamically control their interference pattern⁴⁶.

The reconstruction of the quantum state of single electronic excitations is a currently active research field as illustrated by the very recent achievement of the state tomography of high energy electrons propagating along quantum Hall edge channels⁴⁷,⁴⁸. These recent works highlight the importance of characterizing the purity of the emitted states. Importantly, our protocol fully captures the differences between pure states and statistical mixtures and provides a quantitative evaluation of the purity. This ability to quantify the purity of quantum states generated by electronic sources is crucial for future applications of quantum electronics. More specifically, we find that for single-charge Lorentzian pulses, which are predicted to generate a pure single-electron wavefunction at zero temperature, poisoning by thermal excitations results in the emission of a mixture (with purity $P = 0.75$) of two different states, which correspond to the $n = 1$ and $n = 2$ Lorentzian wavepackets $\phi_{L,n}$. For two-electron Lorentzian pulses, we show that thermal effects lead to the generation of a mixture between the zero temperature Slater determinant formed from $\phi_{L,1}$ and $\phi_{L,2}$ and an undesired two electron state which we fully characterize.

The generation and characterization of multi-electron states in quantum conductors also opens the way to the study of correlations and interactions between a controlled number of excitations emitted on demand in the circuit, with applications to the controlled generation of entangled electron or electron/hole⁴⁹ pairs. In this context, this protocol can also be applied to identify single-particle wavefunctions generated in interacting conductors⁴¹ and supplemented by other measurements⁵⁰, to quantify the importance of interaction-induced quantum correlations.

Finally, it can establish a bridge between electron and microwave quantum optics⁵¹,⁵² by probing the electronic content of microwave photons injected from a transmission line into a quantum conductor.

## Methods

**Sample and noise measurements.** The sample is a GaAs/AlGaAs two-dimensional electron gas of charge density $n_s = 1.9 \times 10^{15} \text{ m}^{-2}$ and mobility $\mu = 2.4 \times 10^6 \text{ cm}^{-2} \text{ V}^{-1} \text{ s}^{-1}$. It is placed in a high magnetic field to reach the quantum Hall regime at filling factor $v = 2$ ($B = 3.7$ T) and $v = 3$ ($B = 2.6$ T). The measurements of the Wigner distribution of sinusoidal drives have been performed at $v = 2$. The measurements of the Wigner distributions of Lorentzian pulses have been performed at $v = 3$. The current noise at the output of the quantum point contact is converted to a voltage noise via the quantum Hall edge channel resistance $R_v = h/ve^2$ between output 3 and the ohmic ground (see Fig. 2). In order to move the noise measurement frequency in the MHz range, $R_v$ is connected to an LC tank circuit of resonance frequency $f_0 = 1.45$ MHz. The tank circuit is followed by a pair of homemade cryogenic amplifiers followed by room-temperature amplifiers. A vector signal analyzer measures the correlations between the voltages at the output of the two amplification chains in a 78-kHz bandwidth centered on $f_0$. The noise measurements are calibrated by measuring both the thermal noise of the output resistance $R_v$, and the partition noise of a d.c. bias applied on input 2 of the electronic beam splitter.

**Generation of Lorentzian current pulses.** The single-electron and two-electron periodic trains of Lorentzian current pulses are generated by applying the ac part of the signal $V_{ac}(t)$ on the mesoscopic capacitor placed at input 1 of the beam splitter, see Fig. 1, and the dc part of the signal $V_{dc}$ on the ohmic contact connected to input 1 of the beam splitter such that $V_S(t) = V_{ac}(t) + V_{dc} = \sum_{l} \frac{V_0}{1+(t-IT)^2/\tau^2}$. More precisely, $V_{ac}(t)$ is generated harmonic by harmonic (from $n = 1$ to $n = 5$):
$$V_{ac}(t) = \sum_{n=1}^{n=5} V_{ac,n} \cos(n\Omega t + \phi_n).$$
The careful calibration of the amplitude $V_{ac,n}$

and phase $\phi_n$ of each harmonic is performed at the location of the beam-splitter using low-frequency noise measurements. The calibration of each amplitude is performed by sending a single harmonic $V_{ac,n} \cos(2\pi n t + \phi_n)$ toward the splitter and measuring the low-frequency noise as a function of $V_{ac,n}$. The calibration of the relative phases $\phi_n$ is more difficult and involves two-particle interferences between two harmonics at two different frequencies. As an example, to calibrate the relative phase $\phi_n$ between the first and $n$th harmonic of the signal, we generate the voltage $V_S(t) = V_{ac,1} \cos(2\pi n t) + V_{ac,n} \cos(2\pi n t + \phi_n)$ at input 1 of the splitter. From the two-particle interference effect, the noise at the splitter output depends on the relative phase $\phi_n$ between the two harmonics. It is minimal (or maximal depending on the harmonics considered) when the two harmonics are in phase, allowing for an accurate calibration of the relative phase between the different harmonics. Note that for even harmonics, the two-particle interferences between the two harmonics vanish at zero bias voltage such that a small bias voltage of a few tens of microvolts needs to be added for their phase calibration. The fact that amplitude and phase of each harmonic are calibrated at the splitter location, is very important regarding the effect of Coulomb interaction during the propagation of single-electron excitations. Previous works have emphasized⁵³ the importance of these effects in quantum Hall edge channels and their impact on the relaxation and decoherence⁵⁴ of electronic excitations or on the fractionalization⁵⁵ of these excitations. In our experimental setup, single or two-electron Lorentzian wavepackets will eventually fractionalize²⁶ due to Coulomb interaction effect. However, this fractionalization will occur after the beam splitter, where we performed the tomography experiment. Indeed, these single-electron excitations result from the generation of voltage pulses propagating along the edge channels. For such types of electronics sources, which are qualitatively different from quantum dot emitters where electron emission involves the tunneling through a transmission barrier, Coulomb interaction can be taken into account by a renormalization²⁶ of the amplitude and phase of each harmonic of the voltage pulse. By calibrating these amplitudes and phases at the splitter location, it means that we absorb the effect of Coulomb interaction by accommodating the amplitudes and phases of the signal to reconstruct at the splitter a Lorentzian pulse that is only limited by our calibration accuracy. This allows us to neglect Coulomb interaction effects in this experiment contrary to previous experiments³⁸,⁵⁴ performed with ac-driven quantum dots.

**Reconstruction of $\Delta W_{S,n}(\omega)$ from noise measurements.** In this paper, we implement a reconstruction of the source Wigner distribution $W_S(t, \omega)$ from the measurement of the current noise $\Delta S_n$ that does not rely on any assumption on the electronic state generated by the source.

First, the excess electronic distribution function $\Delta W_{S,n=0}(\omega)$ can be obtained via the derivative of the noise $\Delta S_{n=0}$ with respect to the d.c. bias $\omega_{dc} = -eV_{dc}/\hbar$ applied on the probe port 2.

$$\Delta S_{n=0} = 2e^2 T (1-T) \int \frac{d\omega}{2\pi} \left[ \Delta \overline{W}_S^{-t} \left( 1 - 2f_{eq}(\omega - \omega_{DC}) \right) \right] \quad (7)$$

$$\Delta \overline{W}_{S,0} = - \frac{\pi}{2e^2 T} \frac{\partial \Delta S_{n=0}}{\partial \omega_{DC}} = \int d\omega \Delta W_{S,0}(\omega) \left( \frac{-\partial f_{eq}}{\partial \omega} \right) (\omega - \omega_{DC}) \quad (8)$$

As shown by Eq. (8), the experimental signal $\Delta \overline{W}_{S,0}$ does not provide directly $\Delta W_{S,0}$ but its convolution with the thermally broadened function $\left(\frac{-\partial f_{eq}}{\partial \omega}\right)$. Knowing the electronic temperature, one can reconstruct $\Delta W_{S,0}(\omega)$ from the measurement of $\Delta W_{S,0}$ using Bayesian deconvolution techniques presented in the next section.

We then turn to the higher-order terms, $\Delta W_{S,n=0}$, which encode all the time dependence. They are reconstructed from the measurement of the output noise $\Delta S_{n,\phi}$ as a function of the d.c. voltage $V_{DC}$ and the phase $\phi$ of the a.c. voltage $V_{P_n}(t)$ applied on the probe port. For a sinusoidal drive at frequency $nf$ of small amplitude $V_{P_n}$ on top of a d.c. drive $V_{DC}$, the probe excess Wigner distribution is given by²²
$$\Delta W_{P_n}(t, \omega) = - \frac{eV_{P_n}}{2\pi} \cos(2\pi n t f + \phi) h_n(\omega - \omega_{DC}), \quad \text{with } h_n(\omega) = (f_{eq}(\omega - n\pi f) - f_{eq}(\omega + n\pi f))/(2\pi n f) \text{ and } \omega_{DC} = -eV_{DC}/\hbar.$$ Inserting this expression for the probe Wigner distribution in Eq. (3) we can reconstruct the real and imaginary parts of $\Delta W_{S,n}$:

$$\Re(\Delta \overline{W}_{S,n}) = \frac{h}{8e^3 V_{P_n} T (1-T)} \left( \Delta S_{n,\phi=0} - \Delta S_{n,\phi=\pi} \right) \\ = \int d\omega \Im(\Delta W_{S,n}(\omega)) h_n(\omega - \omega_{DC}) \quad (9)$$

$$\Im(\Delta \overline{W}_{S,n}) = \frac{h}{8e^3 V_{P_n} T (1-T)} \left( \Delta S_{n,\phi=\frac{\pi}{2}} - \Delta S_{n,\phi=\frac{3\pi}{2}} \right) \\ = \int d\omega \Im(\Delta W_{S,n}(\omega)) h_n(\omega - \omega_{DC}) \quad (10)$$

As in the $n = 0$ case, the experimental signal is the convolution between $\Delta W_{S,n}$ and $h_n$. The real and imaginary parts of $\Delta W_{S,n}$ are thus reconstructed using deconvolution techniques (see next section). A specific difficulty arises for the $n \neq 0$ terms, as their reconstruction process requires the accurate knowledge of amplitude and phase of the probe signals for various values of $n$. The amplitude and phase calibration of all the probe signals $V_{P_n}(t)$ is performed similarly to the calibration
of the amplitude and phase of the harmonics of the Lorentzian voltage pulses. As a result of the phase calibration, we find that, as theoretically expected for the sine and Lorentzian drives, $\Im(\Delta\tilde{W}_{S,n}) = 0$ for all $n$ and all $\omega$. Note that, as mentioned above, Coulomb interaction only modifies the phase and amplitude of sinusoidal drives. As these are calibrated at the splitter location, it means that we can simply ignore Coulomb interaction effects on the probe signals. Finally, when measuring the $n \neq 0$ harmonics, we also systematically checked the linear dependence of the output noise with the probe amplitude in order to check the validity of the linear approximation relating $\Delta W_{P,n}(t, \omega)$ to $V_{P,n}(t)$.

**Bayesian deconvolution method.** The relation between $\Delta W_{S,n}(\omega)$ and $\Delta\tilde{W}_{S,n}(\omega)$ is given by the convolution product, see Eqs. (8)–(10):

$$ \Delta\tilde{W}_{S,n}(\omega) = (h_n * \Delta W_{S,n})(\omega) \quad (11) $$

$$ = \int d\omega' h_n(\omega' - \omega) \Delta W_{S,n}(\omega') \quad (12) $$

In order to estimate $\Delta W_{S,n}(\omega)$ based on $\Delta\tilde{W}_{S,n}(\omega)$, we need to implement a deconvolution algorithm. Since deconvolution is an ill-posed problem, simply performing a division in Fourier space leads to an estimation, which is not robust to measurement errors. The sensibility to errors, due to lost information, corresponds to zero or close to zero values of the Fourier transform of $h_n$. In order to find a more robust estimation with correct physical properties, we propose to add appropriate prior information on $\Delta W_{S,n}(\omega)$ thanks to a Bayesian framework$^{56-58}$.

The discretized forward model for convolution (12) can be expressed as

$$ \Delta\tilde{W}_{S,n} = H_n \cdot \Delta W_{S,n} + N_n, \quad (13) $$

where bold characters stand for vectors and matrices resulting from discretization. $\Delta\tilde{W}_{S,n}$ is the vector of data points, $H_n$ is the convolution matrix, and $\Delta W_{S,n}$ is the unknown quantity we are looking for. The term $N_n$ is added to take account for all the errors (measurement and discretization). It is modeled as Gaussian random vector, with a known covariance matrix $V_e$ with diagonal elements $V_{e,d}$ estimated, thanks to repeated experiments. This gives the expression of the probability distribution of $\Delta\tilde{W}_{S,n}$ knowing $\Delta W_{S,n}$ and $V_e$, which is called the likelihood

$$ p(\Delta\tilde{W}_{S,n} | \Delta W_{S,n}, V_e) \propto \exp\left(-\frac{1}{2} \|\Delta\tilde{W}_{S,n} - H_n \cdot \Delta W_{S,n}\|_{V_e}^2\right) \quad (14) $$

where $\|\mathbf{x}\|_{V_e}^2 = \sum_i \frac{v_e^2}{v_i^2}$. Finding the argument that maximizes the likelihood, is equivalent to performing a division in Fourier space, since the convolution matrix $H_n$ is diagonal in the Fourier basis. This argument is dominated by $H_n^{-1}N_n$ terms. In the Bayesian framework, by adding a prior information, we want to enforce physical properties such that $\Delta W_{S,n}(\omega)$ tends to zero when $|\omega|$ increases. For this purpose, we assign a Gaussian prior distribution on $\Delta W_{S,n}$:

$$ p(\Delta W_{S,n} | V_f) \propto \exp\left(-\frac{1}{2} \|\Delta W_{S,n}\|_{V_f}^2\right). \quad (15) $$

For variances $V_f$, we use the expression

$$ V_f(\omega) = v_f \exp\left(-\frac{\omega^2}{\omega^2}\right), \quad (16) $$

where $v_f$ and $\omega$ are parameters tuned to enforce limit condition when $|\omega|$ increases (the influence of the parameters $w$ and $v_f$ on the deconvoluted signal are presented in the Supplementary Note 4).

Applying Bayes' rule, the posterior probability distribution of $\Delta W_{S,n}$ combines likelihood (14) and prior distribution (15):

$$ p(\Delta W_{S,n} | \Delta\tilde{W}_{S,n}, V_e, V_f) = \frac{p(\Delta\tilde{W}_{S,n} | \Delta W_{S,n}, V_e) p(\Delta W_{S,n} | V_f)}{p(\Delta\tilde{W}_{S,n} | V_e, V_f)} \quad (17) $$

The argument that maximizes this posterior distribution (17), is the most likely estimate of $\Delta W_{S,n}$ knowing both the measurement results $\Delta\tilde{W}_{S,n}$, $V_e$ and prior information encoded in $V_f$. Indeed, this maximum a posteriori (MAP), which also in this case is the posterior mean, is robust to errors $N_n$. Because of model evidence, the term in the denominator of (17), $p(\Delta\tilde{W}_{S,n} | V_e, V_f)$ does not depend on $\Delta W_{S,n}$, MAP estimate becomes equivalent to the minimization of the criterion:

$$ J(\Delta W_{S,n}) = \frac{1}{2} \| \Delta\tilde{W}_{S,n} - H_n \cdot \Delta W_{S,n} \|_{V_e}^2 + \frac{1}{2} \| \Delta W_{S,n} \|_{V_f}^2 \quad (18) $$

The estimated $\Delta W_{S,n}$ has to comply with a box constraint given by Pauli exclusion principle and Cauchy–Schwartz inequality$^{22}$. Consequently, the implemented algorithm in Box 1 looks for a minimum of criterion (18) inside the box constraint, thanks to a Projected Gradient Descent method$^{59}$. In this algorithm, $M_n(\omega)$ denotes Cauchy–Schwartz inequality bounds $\forall n, \omega$ and is explicitly given in ref. $^{22}$.

## Box 1 | Detailed deconvolution algorithm

Compute Cauchy-Schwartz bounds $M_n(\omega)$

Choose amplitude $v_f$ and width of prior (16) for $V_f(\omega)$

Compute the minimum of criterion (18)

$$ \Delta W_{S,n} = (H^T V_e^{-1} H + V_f^{-1})^{-1} H^T V_e^{-1} \Delta\tilde{W}_{S,n} $$

Project the solution inside the box given by Cauchy-Schwartz bounds:

$$ \Delta W_{S,n}(\omega) = \min(\Delta W_{S,n}(\omega), M_n(\omega)) $$

$$ \text{and } \Delta W_{S,n}(\omega) = \max(\Delta W_{S,n}(\omega), -M_n(\omega)) $$

**repeat**

Compute the gradient of criterion (18)

$$ \nabla(\Delta W_{S,n}) = -H^T V_e^{-1} (\Delta\tilde{W}_{S,n} - H\Delta W_{S,n}) + V_f^{-1} \Delta W_{S,n} $$

Project the gradient $\nabla(\Delta W_{S,n})$ to stay in the box constraint if $|\Delta W_{S,n}(\omega)| \ge M_n(\omega)$ and $\Delta W_{S,n}(\omega) * \nabla(\Delta W_{S,n})(\omega) \le 0$ then $\nabla(\Delta W_{S,n})(\omega) = 0$

**else**

$$ \nabla(\Delta W_{S,n})(\omega) = \nabla(\Delta W_{S,n})(\omega) $$

**end if**

Compute the furthest displacement $d_\infty$ in the box along $\nabla(\Delta W_{S,n})$ direction

**For all** $\nabla(\Delta W_{S,n})(\omega) \neq 0$ **do**

$$ d_\infty := \min \left( \frac{M_n(\omega) - \Delta W_{S,n}(\omega)}{\nabla(\Delta W_{S,n})(\omega)}, d_\infty \right) $$

**end for**

Compute the optimum displacement $d_0$ along $\nabla(\Delta W_{S,n})$ direction

$$ d_0 = \| \nabla(\Delta W_{S,n}) \|^{-2} \left( \| H \nabla(\Delta W_{S,n}) \|_{V_e}^2 + \| \nabla(\Delta W_{S,n}) \|_{V_f}^2 \right) $$

Compute one projected descent gradient step

$$ \Delta W_{S,n} = \Delta W_{S,n} - \min(d_0, d_\infty) \nabla(\Delta W_{S,n}) $$

until $-\ln(p(\Delta W_{S,n} | \Delta\tilde{W}_{S,n}, V_e, V_f))$ is minimized

**Electron and hole wavefunction extraction.** The extraction of electron and hole wavefunctions from the experimental data for $\Delta W_S(t, \omega)$ relies on an algorithm that recasts any excess T-periodic single-electron coherence under the form given by Eq. (2) of the article. The algorithm is a generalization of the Kahrunen–Loève analysis$^{60}$ to electron quantum optics. It is based on an exact diagonalization of the projections of the single-electron coherence (represented by $\Delta_0 W_S(t, \omega)$) onto the electronic and hole quadrants with respect to the reference chemical potential (here $\mu = 0$). As explained in ref. $^{22}$, these projections are defined by decomposing the space of a single-particle state into positive (for electrons) and negative energy (for holes) states.

In practice, they are obtained through the following procedure: after deconvolution, the experimental data come as a finite set of real values $\Delta W_{S,n}(\omega_t)$ where $\omega_k$ are the discretized values of $\omega$ and $n = 0, \pm 1, \pm 2, \dots, \pm 5$. First, we add the thermal excess $f_{eq}(\omega) - \Theta(-\omega)$ of the equilibrium Fermi–Dirac distribution at temperature $T_a$ to $\Delta W_{S,n}=0(\omega)$ to obtain the experimental data set for $\Delta_0 W_S(\omega)$. Then, in order to extract a square matrix for the exact diagonalization, the next step is to interpolate the data on a grid well suited to the electronic and hole quadrants. These two quadrants are defined in the frequency domain as corresponding to the sectors, where purely electronic (resp. purely hole) excitations contribute to $\Delta W_S$, $n(\omega)$. For a periodically driven source, they correspond to $\omega \ge |n|\pi f$ for the electron quadrant and to $\omega \le -|n|\pi f$ for the hole quadrant$^{22}$. For each $n$, the data set $\Delta_0 W_S$, $n(\omega_k)$ is first interpolated using cubic splines to infer a new data set on a grid adapted to the electronic and hole quadrants (i.e., such that this grid intersects the boundaries $\omega \pm n\pi f = 0$ of the electronic and hole quadrants). This new data grid has a discretization step $\delta\omega$ such that $h\delta\omega \simeq 0.19 \mu eV$. This data set is then used to build the matrices corresponding to the projections on the electron and hole quadrants of this interpolated data for $\Delta_0 W_S(\omega)$.

Due to time periodicity of the single-electron Wigner function $\Delta_0 W_S(t, \omega) = \Delta_0 W_S(t + T, \omega)$, diagonalizing these two projections onto the electron and hole quadrants leads to electronic ($\alpha = e$) and hole ($\alpha = h$) probability spectral bands $g_i^{(\alpha\alpha)}(\nu)$ ($i$ being a band index) depending on quasi-pulsation interval $0 \le \nu \le 2\pi f$ associated with the time period $T$. The corresponding eigenvectors are the Floquet version of Bloch waves of solid-state physics.

Then, the electron and hole wavefunctions $\phi_{i,l}^{(\alpha)}$ generated at each time period$^{14,61}$ are the analogous of the Wannier functions$^{62}$. They consist of normalized single-particle wavepackets such that

$$ \phi_{i,l}^{(\alpha)}(t) = \phi_i^{(\alpha)}(t - IT) \quad (19) $$

$$ \langle \phi_{i,l}^{(\alpha)} | \phi_{i,l}^{(\alpha)} \rangle = \delta_{i,l} \delta_{\alpha,\alpha'} \delta_{l,l'} \quad (20) $$
where $\alpha$ denotes the electron or hole label and $i$ the band index. These electron and hole wavefunctions are therefore very well suited to describe the excitations generated by time-periodic electron beams. In solid-state physics$^{63}$, the Wannier wavefunctions are not uniquely defined since one can impose an arbitrary quasi-momentum dependent phase in front of each Bloch wave. Here, the same problem is present, and exactly as in solid-state physics, this ambiguity is lifted by minimizing their time spreading. This provides us electronic and hole atoms of signals that are maximally localized in the time domain. Finally, the electron coherence $g_i^{(ee)}(l)$ between electronic atoms of signals translated by $l \in \mathbb{Z}^*$ time periods as well as the hole coherence $g_{i,l}^{(hh)}(l)$ can be obtained from the probability spectra through Fourier transform. For $l \in \mathbb{Z}$

$$g_i^{(aa)}(l) = \int_0^{2\pi f} e^{i\nu l/f} g_i^{(\alpha)}(\nu) \frac{d\nu}{2\pi f} \qquad (21)$$

Note that there is no electronic coherence between atoms of a signal of the same type of excitation, but with a different band index. The electron/hole coherences $g_{i,l}^{(eh)}(l)$ are defined as the single-electron coherence between the electronic atom of signal $\varphi_{i,l}^{(e)}$ and the hole atom of signal $\varphi_{i,r}^{(h)}$:

$$g_{i,l}^{(eh)}(l) = \mathrm{Tr} \left( \psi |\varphi_{i,l}^{(e)}\rangle \psi^\dagger |\varphi_{0,r}^{(h)}\rangle \right) \qquad (22)$$

$$= \int dt d\tau \varphi_{i,l}^{(e)}(t) \ast \Delta_0 G(t, t') \varphi_{0,l}^{(h)}(t'). \qquad (23)$$

where $l \in \mathbb{Z}$, $i$ and $l'$ are possibly different. It is obtained from the electronic Wigner function using the explicit numerical data for the electronic and hole wavefunctions.

**Purity indicator.** The general expression of the state purity is beyond the scope of this paper. We focus here on the two limiting cases considered in this paper: the sinusoidal drives where one electron and one hole wavefunction are generated and the periodic train of single-electron Lorentzian pulses, where two or three electronic wavefunctions need to be considered (the probability for hole emission can be neglected).

Let us consider first the sinusoidal drive case with only one electron and one hole branch so that the branch index $i$ can be dropped out. Time periodicity implies that only coherences between Floquet-Bloch eigenvectors with the same quasi-pulsation do not vanish. As only one electron and one hole wavepackets are emitted, we can therefore consider, at each given quasi-pulsation $0 \le v < 2\pi f$, the reduced density matrix $\rho_{eh}$ in the occupation number basis of the electron and hole states: $|n_e n_h\rangle$, $n_e = 0$ or 1 and $n_h = 0$ or 1 are the occupation numbers of the corresponding single-particle states. As a result of the superposition rule that forbids quantum superposition between states with fermion numbers of different parity$^{64,65}$, $\rho_{eh}$ corresponds to a pure state in three situations: either the electronic and hole levels are both filled (state $|11\rangle$), or both empty (state $|00\rangle$), or populated in a coherent way (state $|u|01\rangle + |v|10\rangle$ with $|u|^2 + |v|^2 = 1$). Any deviation from purity thus reflects incoherent electron/hole processes. The purity indicator, which is defined as $\mathrm{Tr}(\rho_{eh}^2)$, is a good quantity for measuring the weight of coherent processes:

$$\mathrm{Tr}(\rho_{eh}^2) = 1 - 2A(\nu)(1 - A(\nu)) - 2B(\nu)(1 - B(\nu)) \qquad (24)$$

where

$$A(\nu) = g^{(ee)}(\nu)(1 - g^{(hh)}(\nu)) - |g^{(eh)}(\nu)|^2 \qquad (25)$$

$$B(\nu) = g^{(hh)}(\nu)(1 - g^{(ee)}(\nu)) - |g^{(eh)}(\nu)|^2 \qquad (26)$$

are computed in terms of the eigenvalues $g^{(ee)}(\nu)$ and $g^{(hh)}(\nu)$ obtained from our diagonalization algorithm and of the corresponding electron/hole coherences $g^{(eh)}(\nu)$.

Let us now discuss all the reduced density matrices $\rho_{eh}$ for all $0 \le v < 2\pi f$. When Wick's theorem is valid, there are no correlations between different quasi-pulsations $v_1 \ne v_2$. Consequently, we can take the infinite-dimensional product over all the Floquet-Bloch pairs of electron and hole modes for $0 \le v < 2\pi f$ and take the trace of its square, which is a formal infinite product over all $0 \le v < 2\pi f$. This quantity has to be regularized in the infrared by discretizing the quasi energies $v_n = 2\pi f n/N$ for $n = 0...N-1$ and taking the $1/N$th power of the result. This procedure leads to the resulting quantity

$$\mathbb{P} = \exp \left[ \int_{0}^{2\pi f} \ln(1 - 2A(\nu)(1 - A(\nu)) - 2B(\nu)(1 - B(\nu))) \frac{d\nu}{2\pi f} \right] \qquad (27)$$

which is equal to unity if and only if the many-body state is pure and obtained from a Fermi sea vacuum by adding on top of it coherent superposition of electron/hole pairs. When there are incoherent processes such as in the case of nonzero temperature, $\mathbb{P} < 1$. In the present situation, the condition $A(\nu) = B(\nu) = 0$, which ensures unit purity, corresponds to $|g^{(eh)}(\nu)|^2 = g^{(ee)}(\nu)(1 - g^{(hh)}(\nu)) = g^{(hh)}(\nu)(1 - g^{(ee)}(\nu))$, which implies that $g^{(ee)}(\nu) = g^{(hh)}(\nu)$. Then, the condition on $g^{(eh)}(\nu)$ ensures that, for each quasi-pulsation $v$, we have acted on the Fermi sea $|F\rangle$ through the coherent sum of the identity operator, and of the elementary electron/hole pair creation operator that is the product of a creation operator for the electron single-

particle state and of a destruction operator for the hole single-particle state. This is equivalent to putting each quasi-particle which, in $|F\rangle$, is in the hole state $\varphi_\nu^{(h)}$, into the linear combination $u(\nu)\varphi_\nu^{(h)} + v(\nu)\varphi_\nu^{(e)}$. The resulting many-body state is then pure and of the form

$$|\Psi\rangle = \prod_{0 \le v < 2\pi f} (u(\nu) + v(\nu)\psi^\dagger[\varphi_\nu^{(e)}]\psi[\varphi_\nu^{(h)}]) |F\rangle. \qquad (28)$$

This specific form was also obtained in ref. $^{30}$, which considered a conductor at zero temperature described by a single-particle time-dependent scattering matrix. It reduces to the form given by Vanevic et al.$^{29}$ in the case where the Floquet-Bloch spectrum is flat as a function of $v$ (in which case there are no interperiod electronic and hole coherences).

Equation (27) for the purity can be adapted to the Lorentzian case we study in the paper. Three electronic bands need to be considered at most for the Lorentzian pulse carrying two electrons. These bands are not coupled in our specific experimental situation (the term coupling the two bands are the electron-hole coherences $g_{ij}^{(eh)}$, which we measure to be negligible). In this case, the extension of Eq. (27) is straightforward:

$$\mathbb{P} = \mathbb{P}_1 \times \mathbb{P}_2 \times \mathbb{P}_3 \qquad (29)$$

$$\mathbb{P}_i = \exp \left[ \int_0^{2\pi f} \ln \left( 1 - 2g_i^{(ee)}(\nu)(1 - g_i^{(ee)}(\nu)) \right) \frac{d\nu}{2\pi f} \right] \qquad (30)$$

where we have used the simplified expressions of $A_i(\nu) = g_i^{(ee)}(\nu)$ and $B_i(\nu) = 0$ in the case where $g_i^{(hh)}(\nu) \approx g_{ij}^{(eh)}(\nu) \approx 0$. The pure state $\mathbb{P} = \mathbb{P}_i = 1$ is only

recovered for $g_i^{(ee)}(\nu) = 0$ or 1. In our experimental situation for the single-electron Lorentzian pulse, two wavefunctions are emitted with probabilities smaller than 1 and we find $\mathbb{P} = 0.75$, showing that the generated state is a mixture of two single-electron wavefunctions. For the two-electron Lorentzian pulse, we find a slightly smaller purity $\mathbb{P} = 0.68$.

## Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Received: 23 November 2018 Accepted: 4 July 2019

Published online: 29 July 2019

## References

1. Eisaman, M. D., Fan, J., Migdall, A. & Polyakov, S. V. Single-photon sources and detectors. Rev. Sci. Instrum. 82, 071101 (2011).

2. Barredo, D., de Léséleuc, S., Lienhard, V., Lahaye, T. & Browaeys, A. An atom-by-atom assembler of defect-free arbitrary two-dimensional atomic arrays. Science 354, 1021-1023 (2016).

3. Chou, C.-W. et al. Preparation and coherent manipulation of pure quantum states of a single molecular ion. Nature 545, 203-207 (2017).

4. Fève, G. et al. An on-demand coherent single electron source. Science 316, 1169-1172 (2007).

5. Leicht, C. et al. Generation of energy selective excitations in quantum Hall edge states. Semicond. Sci. Technol. 26, 055010 (2011).

6. Dubois, J. et al. Minimal-excitation states for electron quantum optics using levitons. Nature 502, 659-663 (2013).

7. Fletcher, J. D. et al. Clock-controlled emission of single-electron wave packets in a solid-state circuit. Phys. Rev. Lett. 111, 216807 (2013).

8. Roussely, G. et al. Unveiling the bosonic nature of an ultrashort few-electron pulse. Nat. Commun. 9, 2811 (2018).

9. Bertoni, A., Bordone, P., Brunetti, R., Jacoboni, C. & Reggiani, S. Quantum logic gates based on coherent electron transport in quantum wires. Phys. Rev. Lett. 84, 5912-5915 (2000).

10. Yamamoto, M. et al. Electrical control of a solid-state flying qubit. Nat. Nanotechnol. 7, 247-251 (2012).

11. Smithey, D. T., Beck, M., Raymer, M. G. & Faridani, A. Measurement of the Wigner distribution and the density matrix of a light mode using optical homodyne tomography: application to squeezed states and the vacuum. Phys. Rev. Lett. 70, 1244-1247 (1993).

12. Lvovsky, A. I. & Raymer, M. G. Continuous-variable optical quantum-state tomography. Rev. Mod. Phys. 81, 299-332 (2009).

13. Marguerite, A. Two-particle interferometry in quantum Hall edge channels. Phys. Status Solidi B 254, 1600618 (2017).

14. Roussel, B. Autopsy of a Quantum Electrical Current. PhD thesis, Université de Lyon (2017).

15. Levitov, L. S., Lee, H. & Lesovik, G. Electron counting statistics and coherent states of electric current. J. Math. Phys. 37, 4845-4866 (1996).
16. Keeling, J., Klich, I. & Levitov, L. S. Minimal excitation states of electrons in one dimensional wires. *Phys. Rev. Lett.* **97**, 116403 (2006).

17. Gabelli, J. & Reulet, B. Shaping a time-dependent excitation to minimize the shot noise in a tunnel junction. *Phys. Rev. B* **87**, 075403 (2013).

18. Kataoka, M. et al. Time-of-flight measurements of single-electron wave packets in quantum hall edge states. *Phys. Rev. Lett.* **116**, 126803 (2016).

19. Altimiras, C. et al. Non-equilibrium edge-channel spectroscopy in the integer quantum Hall regime. *Nat. Phys.* **6**, 34–39 (2010).

20. Bocquillon, E. et al. Electron quantum optics in ballistic chiral conductors. *Ann. Phys. (Berl.)* **526**, 1–30 (2014).

21. Haack, G., Moskalets, M. & Büttiker, M. Glauber coherence of single-electron sources. *Phys. Rev. B* **87**, 201302 (2013).

22. Ferraro, D. et al. Wigner function approach to single electron coherence in quantum Hall edge channels. *Phys. Rev. B* **88**, 205303 (2013).

23. Kashcheyevs, V. & Samuelsson, P. Classical-to-quantum crossover in electron on-demand emission. *Phys. Rev. B* **95**, 245424 (2017).

24. Wigner, E. On the quantum correction for thermodynamic. *Equilib., Phys. Rev.* **40**, 749–759 (1932).

25. Jullien, T. et al. Quantum tomography of an electron. *Nature* **514**, 603–607 (2014).

26. Grenier, C. et al. Fractionalization of minimal excitations in integer quantum Hall edge channels. *Phys. Rev. B* **88**, 085302 (2013).

27. Ferraro, D. et al. Hong-Ou-Mandel characterization of multiply charged Levitons. *Eur. Phys. J. Spec. Top.* **227**, 1345–1359 (2018).

28. Vanević, M., Nazarov, Y. V. & Belzig, W. Elementary events of electron transfer in a voltage-driven quantum point contact. *Phys. Rev. Lett.* **99**, 076601 (2007).

29. Vanević, M., Gabelli, J., Belzig, W. & Reulet, B. Electron and electron-hole quasiparticle states in a driven quantum contact. *Phys. Rev. B* **93**, 041416 (2016).

30. Yin, Y. Quasiparticle states of on-demand coherent electron sources, Quasiparticle states of on-demand coherent electron sources. *J. Phys. Condens. Matter* **31**, 245301 (2019).

31. Misiorny, M., Feve, G. & Splettstoesser, J. Shaping charge excitations in chiral edge states with a time-dependent gate voltage. *Phys. Rev. B* **97**, 075426 (2018).

32. Ch. Grenier, R. et al. Single-electron quantum tomography in quantum Hall edge channels. *New J. Phys.* **13**, 093007 (2011).

33. Liu, R. C., Odom, B., Yamamoto, Y. & Tarucha, S. Quantum interference in electron collision. *Nature* **391**, 263–265 (1998).

34. Neder, I. et al. Interference between two indistinguishable electrons from independent sources. *Nature* **448**, 333–337 (2007).

35. Hong, C. K., Ou, Z. Y. & Mandel, L. Measurement of subpicosecond time intervals between two photons by interference. *Phys. Rev. Lett.* **59**, 2044–2046 (1987).

36. Ol'khovskaya, S., Splettstoesser, J., Moskalets, M. & Büttiker, M. Shot noise of a mesoscopic two-particle collider. *Phys. Rev. Lett.* **101**, 166802 (2008).

37. Jonckheere, T., Rech, J., Wahl, C. & Martin, T. Electron and hole Hong-Ou-Mandel interferometry. *Phys. Rev. B* **86**, 125425 (2012).

38. Bocquillon, E. et al. Coherence and indistinguishability of single electron wavepackets emitted by independent sources. *Science* **339**, 1054–1057 (2013).

39. Reydellet, L.-H., Roche, P., Glattli, D. C., Etienne, B. & Jin, Y. Quantum partition noise of photon-created electron-hole pairs. *Phys. Rev. Lett.* **90**, 176803 (2003).

40. Rychkov, V. S., Polianski, M. L. & Büttiker, M. Photon-assisted electron-hole shot noise in multiterminal conductors. *Phys. Rev. B* **72**, 155326 (2005).

41. Cabart, C., Roussel, B., Fève, G. & Degiovanni, P. Taming electronic decoherence in one-dimensional chiral ballistic quantum conductors. *Phys. Rev. B* **98**, 155302 (2018).

42. Moskalets, M. Single-particle emission at finite temperatures. *Low. Temp. Phys.* **43**, 865–876 (2017).

43. Moskalets, M. First-order correlation function of a stream of single-electron wavepackets. *Phys. Rev. B* **91**, 195431 (2015).

44. Dashti, N., Misiorny, M., Kheadoud, S., Samuelsson, P. & Splettstoesser, J. Minimal-excitation single-particle emitters: A comparison of charge and energy transport properties. *Phys. Rev. B* **100**, 035405 (2019).

45. Burset, P., Kotilahti, J., Moskalets, M. & Flindt, C. Time-domain spectroscopy of mesoscopic conductors using voltage pulses. *Adv. Quantum Technol.* **2**, 1900014 (2019).

46. Gaury, B. & Waintal, X. Time-Domain Spectroscopy of Mesoscopic Conductors Using Voltage Pulses. *Nat. Commun.* **5**, 3844 (2014).

47. Fletcher, J. D. et al. Quantum tomography of solitary electrons. Preprint at https://arxiv.org/abs/1901.10985 (2019).

48. Locane, E., Brouwer, P. W., and Kashcheyevs, V., Dynamical scattering of single-electron wave packets. Preprint at https://arxiv.org/abs/1901.08940 (2019).

49. Hofer, P. P., Dasenbrook, D. & Flindt, C. On-demand entanglement generation using dynamic single electron sources. *Phys. Status Solidi B* **254**, 1600582 (2017).

50. Thibierge, É. et al. Two-electron coherence and its measurement in electron quantum optics. *Phys. Rev. B* **93**, 081302 (2016).

51. Grimsmo, A. L., Qassemi, F., Reulet, B. & Blais, A. Quantum optics theory of electronic noise in quantum conductors. *Phys. Rev. Lett.* **116**, 043602 (2016).

52. Virally, S., Simoneau, J. O., Lupien, C. & Reulet, B. Discrete photon statistics from continuous microwave measurements. *Phys. Rev. A* **93**, 043813 (2016).

53. Wahl, C., Rech, J., Jonckheere, T. & Martin, T. Interactions and charge fractionalization in an electronic Hong-Ou-Mandel interferometer. *Phys. Rev. Lett.* **112**, 046802 (2014).

54. Marguerite, A. et al. Decoherence and relaxation of a single electron in a one-dimensional conductor. *Phys. Rev. B* **94**, 115311 (2016).

55. Freulon, V. et al. Hong-Ou-Mandel experiment for temporal investigation of single-electron fractionalization. *Nat. Commun.* **6**, 6854 (2015).

56. Ayasso, H. & Mohammad-Djafari, A. Joint NDT image restoration and segmentation using Gauss–Markov–Potts prior models and variational Bayesian computation. *IEEE Trans. Image Process.* **19**, 2265–2277 (2010).

57. Mohammad-Djafari, A. & Dumitru, M. Bayesian sparse solutions to linear inverse problems with non-stationary noise with Student-t priors. *Digit. Signal Process.* **47**, 128–156 (2015).

58. Zhao, N., Basarab, A., Kouamé, D. & Tournetet, J. Y. Joint segmentation and deconvolution of ultrasound images using a hierarchical Bayesian model based on generalized Gaussian priors. *IEEE Trans. Image Process.* **25**, 3736–3750 (2016).

59. Bertsekas, D. P. *Nonlinear Programming* 3rd edn (Athena Scientific, Boston, MA, 2016).

60. Papoulis, A. & Pillai, S. U., Probability, random variables and stochastic processes 4th edn, (McGraw-Hill, New York 2002).

61. Roussel, B., Cabart, C., Fève, G., Thibierge, E. & Degiovanni, P. Electron quantum optics as quantum signal processing. *Phys. Status Solidi B* **254**, 16000621 (2017).

62. Wannier, G. H. The structure of electronic excitation levels in insulating crystal. *Phys. Rev.* **52**, 191–197 (1937).

63. Marzari, N. et al. Maximally localized Wannier functions: theory and applications. *Rev. Mod. Phys.* **84**, 1419–1475 (2012).

64. Wick, G. C., Wightman, A. S. & Wigner, E. P. The intrinsic parity of elementary particles. *Phys. Rev. B* **88**, 101–105 (1952).

65. Johansson, M. Comment on 'Reasonable fermionic quantum information theories require relativity. Preprint at https://arxiv.org/abs/1610.00539 (2016).

## Acknowledgements

This work has been supported by ANR grants "Ishot reloaded" (ANR-14-CE32-0017), the project EMPIR 17FUN04 SEQUOIA, and the ERC consolidator grant "EQuO" (No. 648236).

## Author contributions

Y.J. fabricated the sample on GaAs/AlGaAs heterostructures grown by A.C. and U.G. R.B., A.M. and M.K. conducted the measurements. Deconvolution algorithms were developed by CCh, R.B. and A.M.D. The algorithm for the extraction of electron and hole wavefunctions was developed by B.R., C.Ca, and P.D. R.B., A.M., M.K., J.M.B., E.B., B.P. and G.F. participated in the data analysis. All authors participated in the writing of the paper. G.F. supervised the project.

## Additional information

Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467-019-11369-5.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Peer review information:** Nature Communications thanks Géraldine Haack and other anonymous reviewer(s) for their contribution to the peer review of this work.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution

4.0 International License, which permits use, sharing, adaptation,

distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.