Multi-particle interference in an electronic Mach-Zehnder interferometer

The recent development of dynamic single-electron sources makes it possible to observe and manipulate the quantum properties of individual charge carriers in mesoscopic circuits. Here, we investigate multi-particle effects in an electronic Mach-Zehnder interferometer driven by dynamic voltage pulses. To this end, we employ a Floquet scattering formalism to evaluate the interference current and the visibility in the outputs of the interferometer. An injected multi-particle state can be described by its first-order correlation function, which we decompose into a sum of elementary correlation functions that each represent a single particle. Each particle in the pulse contributes independently to the interference current, while the visibility (determined by the maximal interference current) exhibits a Fraunhofer-like diffraction pattern caused by the multi-particle interference between different particles in the pulse. For a sequence of multi-particle pulses, the visibility resembles the diffraction pattern from a grid, with the role of the grid and the spacing between the slits being played by the pulses and the time delay between them. Our findings may be observed in future experiments by injecting multi-particle pulses into an electronic Mach-Zehnder interferometer.


Introduction
Quantum-coherent circuits based on mesoscopic conductors [1] combined with dynamic single-electron emitters [2,3] have paved the way for experiments on high-frequency quantum transport [4][5][6][7][8][9] and are holding great promises for future quantum technologies. Advances in nanotechnology have made it possible to fabricate highly pure samples and cool them to sub-Kelvin temperatures, where the phase coherence of the charge carriers is preserved over large enough distances to exploit and control their quantum behavior. By utilizing the quantum Hall effect in a strong magnetic field, the electrons can be forced to move along chiral edge states, as if traveling on rail tracks [10][11][12][13]. In addition, several dynamic single-electron emitters have been developed. Single electrons can now be injected into a quantum-coherent circuit using driven mesoscopic capacitors [5,6,14], dynamic quantum dots [15][16][17] or by applying Lorentzian-shaped pulses to a contact [7,8] as first envisioned by Levitov and co-workers [18][19][20]. When operated in the giga-Hertz regime, these setups make it possible to investigate and observe the quantum behavior of a single or a few electrons in a circuit.
Interference of particles is one central prediction of quantum mechanics, which distinguishes it from classical physics. Interference patterns can develop because the states of quantum particles are described by wave functions. The quantum nature of photons was clarified in pioneering work such as the Hanbury Brown-Twiss [21,22] and Hong-Ou-Mandel [23] interference experiments and by the quantum mechanical theory of light by Glauber based on correlation functions [24][25][26]. This understanding led to the development of quantum optics, which uses beam splitters, mirrors, and other optical components to investigate the statistics and coherence of photons, and eventually paved the way for quantum technologies such as quantum communication with photons [27].
In this article, we investigate theoretically the interference of multi-particle excitations emitted into an electronic Mach-Zehnder interferometer building on earlier works on either static voltages [113], periodically modulated interferometers [114], or dynamic charge emitters [115]. We focus in particular on the injection of clean multi-particle pulses into the interferometer, and we show how the visibility measured in the outputs can be related to the excess correlation function of the incoming pulse, which can be further decomposed into elementary contributions from the individual charges making up the pulse. Based on this understanding, we can interpret an observed Fraunhofer-like diffraction pattern as arising due to the interference of the excess correlation functions of various elementary single-electron components of the multi-particle pulse. These findings may be observed in future experiments on electron quantum optics with dynamic charge emitters and electronic Mach-Zehnder interferometers.
The rest of the paper is organized as follows: In Sec. 2, we review the theoretical description of quantum transport in dynamically driven mesoscopic conductors based on the Floquet scattering formalism, and we introduce the notion of correlation functions in electron quantum optics. In Sec. 3, we describe the electronic Mach-Zehnder interferometer and show how the current in the outputs can be divided into a classical contribution and an interference term, which vanishes as the temperature is increased. We also show how the visibility of the interferometer can be expressed in terms of Figure 1. Illustration of a periodically-driven mesoscopic scatterer. Periodic voltage pulses, V(t), are applied to the left electrode, so that electrons are emitted towards the central scatterer, described by the scattering matrix S(E). There, they are either reflected back into the source electrode or transmitted to the right drain electrode, where the electric current is measured. The fermionic operatorsâ † 0 andâ 0 describe electrons in equilibrium at temperature T and chemical potential µ, before the application of the voltage pulses. The operatorsâ † andâ describe electrons after the application of the voltage pulses, while the operatorsb † andb describe electrons that have been transmitted through the scatterer. The different operators are related via the Floquet scattering matrix according to Eqs. (9) and (10). The electric current before and after the scatterer are denoted as I in (t) and I out (t), respectively. the excess correlation function of the injected charge pulses. Section 4 is devoted to our results and analysis of the visibility of the interferometer, which we interpret in terms of interferences between the different single-particle components of the multi-charge pulse, giving rise to a Fraunhofer-like diffraction pattern. Finally, in Sec. 5, we describe the main conclusions of our work.

Floquet scattering theory
Mesoscopic conductors driven by a periodic voltage can be described using Floquet scattering theory [112]. We thus consider systems governed by a time-periodic Hamiltonian H(x, t) = H(x, t + T ) with period T and frequency Ω = 2π/T . According to the Floquet theorem, a complete set of solutions to the time-dependent Schrödinger equation can be written as [112] We can expand the periodic function φ (x, t) in a Fourier series, such that the wave functions become The wave functions are invariant with respect to the shift E → E + mhΩ, with m ∈ Z. Therefore, the energy E is a quasi-energy (Floquet energy), which is defined only up to integer multiples ofhΩ. In the following, we consider a mesoscopic scatterer that is connected by ballistic leads to two electronic reservoirs, and a periodic voltage V(t) = V(t + T ) is applied to one of them, see Fig. 1. Electrons are then emitted from the driven reservoir towards the scatterer from where they are either reflected back into the source electrode, or they are transmitted through the scatterer and eventually reach the drain electrode. An applied voltage can be treated as a spatially uniform potential [112]. The wave function of electrons leaving the reservoir can then be derived from the Schrödinger equation, where ψ 0 (x) is the solution without the voltage pulses, and J(t) is a time-dependent phase factor, Owing to the periodicity of the time-dependent voltage, we can expand J(t) in a Fourier series as with the Fourier coefficients given as As a result, the wave function of electrons leaving the driven contact becomes which is obviously a Floquet wave function as in Eq. (2). Now, based on Eq. (7), the possible energies of an electron after being excited by the periodic voltage are E n = E + nhΩ, where n is a positive or negative integer and E is the energy of the electron before interacting with the voltage pulses.
In what follows, we consider quasi-one-dimensional transport as in the chiral edge states of a quantum Hall system. We assume that the propagating electrons have energies close to the Fermi energy µ, and approximate their dispersion relation as linear, where k F is the Fermi wavevector. Thus, the electrons are assumed to propagate at the Fermi velocity v F . Using second quantization, the operatorsâ † 0,j (E) andâ 0,j (E) describe electrons with energy E in the reservoir labeled j connected to the incoming lead, before being excited by the voltage pulses. These electrons are in thermal equilibrium, allowing us to write their quantum-statistical average as where |0 is the state of the undisturbed Fermi sea, and f j (E) = 1/(e β(E−µ) + 1) is the Fermi-Dirac distribution with β = 1/(k B T) being the inverse electronic temperature. We also define the creation and annihilation operatorsâ † j andâ j for electrons incident on the scatterer from lead j, after being excited by the voltage. They are related to the equilibrium ones aŝ Next, we describe the transmission and reflection of charge pulses on the central scatterer. The operators b † i (E) and b i (E) for scattered electrons, with i labeling the reservoirs, are related to the incident ones asb where S(E) is the scattering matrix. Then, by combining Eqs. (9) and (10), we can express the operators for the scattered electrons aŝ whereby we can identify the product S(E n )J n as the Floquet scattering matrix, for a static conductor driven by periodic voltage pulses. The Floquet scattering matrix contains the probability amplitudes for an electron with energy E to change its energy to E n = E + nhΩ by exchanging n modulation quanta of sizehΩ with the voltage and to be transmitted through the scatterer. For our purposes, it is useful that the effects of the voltage drive and the scatterer can be factorized as in Eq. (12), which would not be the case if the scatterer itself was modulated in time.
Finally, the electric current generated by a time-dependent voltage reads [104,112,116] Using Eqs. (8) and (11), the time-dependent current can be recast as where the Fourier components of the voltage pulses are defined in Eq. (6).

Voltage pulses
We now specify the types of voltage pulses that we will consider. A current that carries the charge of one electron per period can be created by applying a periodic voltage V(t), which satisfies the condition e 2 h t+T t V(t )dt = e, where T is the period of the drive [20]. For almost any voltage drive, such a current is accompanied by neutral electron-hole pairs, and the resulting charge pulses consist of more than just a single electron. However, by applying Lorentzian-shaped pulses, clean single-electron excitations can be created [18][19][20]. Experimentally, such single-electron pulses were first realized by Dubois and coworkers, who named them levitons [7,117]. A Lorentzian voltage pulse has the form where Γ is the half-width of the pulse, andn controls the average charge per voltage pulse, which can be either integer or non-integer. Technically, we treat single voltage pulses, or a train of a finite number of pulses, within the Floquet scattering formalism by considering a period that is much longer than all other time scales in the problem. We will also consider single Gaussian voltage pulses of the form as well as sinusoidal voltage pulses reading Finally, a sequence of m pulses is given as where V 1 is the voltage for a single pulse, and w is the separation between the pulses.

Excess correlation function
Electron quantum optics experiments aim to realize ideas from quantum optics, however, with electrons in solid-state architectures instead of photons in vacuum or in optical fibres. Therefore, it is relevant to adapt the theoretical toolbox from quantum optics to condensed matter systems. Among the useful tools in quantum optics are the correlation functions introduced by Glauber [24][25][26]. Correlation functions, also known as Green's functions, involve quantum statistical averages of products of field operators. To facilitate the analysis of electron quantum optics experiments, the ideas of Glauber have recently been extended to electronic systems and mesoscopic conductors [63,118,119]. An important difference in the electronic case compared to quantum optics is that, unlike photons which often propagate in vacuum, electrons in mesoscopic conductors propagate on top of the Fermi sea. For this reason, we need to use excess correlation functions, which subtract the effects of the underlying Fermi sea, leaving only correlations between the propagating electrons close to the Fermi level.
The first-order excess correlation functionG (1) contains all information about the single-particle states of the emitted electrons [119]. It also provides some information about multi-particle states, although their full characterization requires higher-order correlation functions [120]. For electrons injected by a voltage drive, the first-order excess correlation function is defined as [63] G (1) where the term with the subscript "on" contains field operators for the case with a time-dependent voltage V(t), and the term with the subscript "off" corresponds to the case where there is no voltage drive. Due to the linear dispersion relation, we can evaluate the correlation function anywhere along a lead, and we generally leave out the spatial dependence in these definitions as well as in wave functions and other related quantities. If the injected electrons are transmitted through a scatterer, we must separately define the correlation functions before and after the scatterer. However, we will focus on the correlation function before the scatterer, which can be related to measurable quantities such as the interference current and the visibility in a Mach-Zehnder interferometer as we will see. To evaluate the correlation function, we first define the electronic field operators that create an electron at position x in lead j aŝ where ψ (in/out) E,j (x, t)/ √ hv F are the incident and scattered wave functions for electrons in lead j with energy E. The prefactor 1/ √ hv F ensures that the wave functions carry a unit flux of particles. Without a voltage applied to the reservoir, the creation and annihilation operators in the lead are simply the same as the reservoir operators,â † 0,j ,â 0,j . However, when a time-dependent voltage is applied to the reservoir, the operators in the lead are related to the reservoir operators according to Eq. (9). To keep the notation simple, we only consider the correlation function for operators in the same lead and henceforth omit the indices of the leads. The correlation function in Eq. (19) then reads Since the wave functions in the lead are plane waves, ψ E (x, t) = e ikx e −iEt/h , we moreover find Using Eq. (9), we can insert equilibrium operators in this expression, which allows us to evaluate the quantum-statistical averages based on Eq. (8). We can perform one of the integrations, which yields This expression holds for any temperature and driving frequency. However, we still need to evaluate the last integral. To this end, we employ an approximation which is valid at low temperatures and driving frequencies compared to the Fermi energy, which is reasonable for mesoscopic conductors. Generally, we need to evaluate energy integrals of the form Changing the integration variable to ε = β(E + nhΩ − µ) allows us to rewrite the integral as where we have assumed that both the temperature and the driving frequency are small compared to the Fermi energy, namely, µβ 1 and µ nhΩ, where n is among the largest numbers of excitation quanta that an electron can emit or absorb because of the time-dependent voltage. As a result, we can extend the lower limit of the integral to minus infinity and evaluate the resulting integral using contour integration in the complex plane, leading to the simple expression Finally, by combining Eq. (23) and Eq. (26), we arrive at the useful approximatioñ To proceed, we note that the sums in Eq. (27) exactly correspond to the Fourier expansions of J * (t 1 ) and J(t 2 ) in Eq. (5). Moreover, by using Eq. (4), we find such that the correlation function can be written as where we have introduced the envelope correlation function [121] The form of the correlation function in Eq. (29) is particularly useful, as it separates the effects of the Fermi energy, the temperature, and the voltage [47,61]. The effects of the Fermi energy enter as the phase factor e iµ(t 1 −t 2 )/h , while the temperature enters via the function χ, which decreases from χ = 1 at zero temperature to χ 0 at large temperatures, where quantum coherence is lost. Finally, the effects of the time-dependent voltage are contained in G (1) through the phase factor Φ(t 1 , t 2 ). It is easy to add the other factors, if needed, and we therefore focus on the envelope correlation function.

Mach-Zehnder interferometer
An electronic Mach-Zehnder interferometer can be used to investigate interference effects of electrons, similarly to an optical Mach-Zehnder interferometer for photons. The setup can be realized using chiral edge states in a Corbino disk geometry [35]. The interferometer is made up of two quantum point contacts that act as beam splitters with energy independent reflection and transmission amplitudes that we denote by r 1 , r 2 and t 1 , t 2 . There are two paths for electrons to travel from the first beam splitter to the next one, as illustrated in Fig. 2, and we denote the length of the upper and lower paths as L u and L d , respectively. In addition, a magnetic flux threads the area enclosed by the paths, causing electrons to acquire the additional phases φ u and φ d depending on the path. In the following, we analyze the current injected in the upper input and the current that reaches the upper output.

Injected current
We start by evaluating the time-dependent current injected in the input. In this case, the scattering matrix in Eq. (14) equals unity, and the integral over energy becomes nhΩ. We then find From Eq.(5), we recognize that this expression can be rewritten as Moreover, using that we arrive at the simple expression for the injected current showing that it is given simply by the conductance quantum, G 0 = e 2 /h, times the the applied voltage.

Output current
Next, we turn to the current in the upper output of the interferometer. The scattering amplitude is given by [35,122] where τ u,d = L u,d /v F are the times it takes for electrons to travel along the upper or lower arm. The two contributions above describe the propagation of electrons in each of the two arms. We now see that the product of scattering matrices can be written as a sum of "classical" and "interference" contributions, with and (38) where ∆τ = τ u − τ d is the time difference between traveling along each arm, and φ = φ u − φ d is the phase difference controlled by the magnetic flux. We have also used the unitarity of the scattering matrices of the two quantum point contacts to deduce that r j t * j = −r * j t j , which leads to r 1 where T j = |t j | 2 and R j = 1 − T j , j = 1, 2. Based on Eq. (14), we see that the current in the output can also be divided into a classical part and an interference contribution. Specifically, we find I out (t) = I cl (t) + I int (t), (39) with the two contributions given by the corresponding terms of the transmission in Eq. (36), For the classical part I cl , there is no energy dependence in the scattering amplitudes, so I cl can be derived in essentially the same way as the injected current in Eq. (34) with the only differences being the prefactors R 1 R 2 and T 1 T 2 and the exponential factors containing τ u and τ d . Hence, we find which is simply a sum of shifted and scaled injected currents, and this result can be understood without resorting to quantum mechanics. By contrast, quantum effects are important for the interference current. For ∆τ = 0, the scattering amplitudes are energy independent, and we find the interference current just as we found the classical one. On the other hand, for ∆τ = 0, the scattering amplitudes are energy dependent, and we find the interference curret using Eq. (26). The interference current then becomes with the envelope correlation function G (1) defined in Eq. (30), and the function χ given in Eq. (26).

Transferred charge
In addition to the output current, we can consider the transferred charge given by the integral which is taken over a period, which for a finite number of pulses is assumed to be much longer than any other time scale in the problem. The transferred charge can again be separated into classical and interference contributions as where Q cl and Q int are given by Eq. (43) by inserting the corresponding currents. From Eq. (41), we see that the classical contribution reads which is simply the charge of the injected current multiplied by R 1 R 2 + T 1 T 2 . The charge transferred due to the interference current follows from Eq. (42) and can be expressed as (46) in terms of the excess correlation function, where we have assumed that ∆τ = 0.

Visibility
A common way to characterize interference is to measure the visibility of a physical observable. This technique can be applied to various situations involving wave interference, from classical electromagnetic waves to quantum mechanical particles. Typically, there is a physical observable which oscillates as a function of a control parameter. The usual definition of visibility in this context is the ratio of the amplitude of the oscillations over the average value of the oscillations. The charge visibility has so far been the main method for characterizing interference effects in Mach-Zehnder interferometers [35].
The charge transferred through a Mach-Zehnder interferometer oscillates as a function of the magnetic flux according to Eqs. (43,45,46). The standard definition of the charge visibility is then where we have to maximize and minimize the transferred charge over the phase difference φ. To evaluate the visibility from Eq. (47), we again decompose the charge as Q = Q cl + Q int and consider the two terms seperately. For ∆τ = 0, the interference term can be expressed as where we have introduced the phase α 1 (φ) = k F ∆L − φ, defined the function and then written G(τ u , τ d ) in polar form with the complex phase α 2 . As a result, the integral reduces to From Eq. (48), we see that Q int only depends on the phase φ through the cosine term, which takes values between -1 and 1. Thus, for the maximum value of Q int , we get while for the minimum, we have min(Q int ) = −max(Q int ).
For ∆τ = 0, we have and therefore together with Eq. (52). Moreover, since Q cl is independent of φ, we immediately get max(Q cl ) = min(Q cl ) = Q cl , and combining these results, we obtain the expression where we have defined the prefactor which is independent of the voltage pulses.

Results & analysis
We are now ready to analyze the charge transferred through the Mach-Zehnder interferometer and the corresponding interference effects. We have already seen how the charge visibility is captured by the absolute value of the time-averaged excess correlation function, cf. Eqs. (49) and (55). We now go on to show how the properties of the correlation function, in particular the possibility to decompose it into elementary components, can explain the interference patterns in the visibility.

Fraunhofer-like diffraction pattern
We start by considering a single voltage pulse, and in Fig. 3 we show the visibility given by Eq. (55) as a function of the difference of travel times ∆τ over the pulse witdh Γ. We consider Lorentzian pulses, Eq. (15), Gaussian pulses, Eq. (16), and sinusoidal pulses, Eq. (17). In each case, we see that the visibility decreases as a function of the difference of travel times. If the difference of travel times is large compared to the width, the parts of the wave function going through different arms do not interfere when they reach the second quantum point contact, leading to a suppression of the visibility. By contrast, for short differences, of travel times compared to the pulse width, a pulse can interfere with itself at the second quantum point contact.
The visibility for Gaussian and sinusoidal pulses are nearly identical and qualitatively similar to the Lorentzian pulses. The difference between the Lorentzian and the other pulses is mainly due to the long tails of the Lorentzian ones, which extend beyond |t| ≥ Γ. The similarity between the Gaussian and the sinusoidal pulses suggests that the exact shape of the pulses only has a small influence on the visibility. Owing to the similar behavior of the different pulses, we focus now on the Lorentzian ones, which can be interpreted in terms of elementary excitations that leave the Fermi sea unaltered.
Continuing with a single pulse, we now allow it to carry more than a single charge. The visibility is quite featureless forn = 1 according to Fig. 3. By contrast, more interesting structures appear as n is increased. Forn > 1, oscillations appear in the visibility, forming a Fraunhofer-like interference pattern (a peak with superimposed interference oscillations in the tails) as seen in Fig. 3, with the number of oscillations being proportional ton. To understand this Fraunhofer-like diffraction pattern, we analyze the correlation function G (1) in more detail.
All voltages in Fig. 3 are proportional ton, and we can write the phase in Eq. (28) as Φ(t 1 , t 2 ) = −nϕ(t 1 , t 2 ) with ϕ being the phase for a single pulse withn = 1. Hence, for integer values ofn, we can decompose the numerator of the correlation function as When combined with Eq. (30), we can then write the correlation function as with We note that g(t 1 , t 2 ) is the correlation function of a single pulse withn = 1. Consequently, Eq. (58) implies that the correlation function G (1) for a pulse with integern > 1 is a sum ofn elementary correlation functions, which differ from each other only by the relative phases kϕ with k = 0, . . . ,n − 1.
The visibility is proportional to the time-integrated correlation function G(τ u , τ d ) according to Eq. (49). We thus denote the correlation function G(τ u , τ d ) for a single pulse with chargen as Gn(τ u , τ d ). Using Eq. (58) and omitting the dependence on τ u and τ d to keep the notation simple, we find where we have defined the real numbers G k and θ k such that G k e iθ k = dtge ikϕ . In analogy with Eq. (58), we can now interpret this expression as a sum of elementary contributions, where each term of the sum has a different phase θ k − θ 0 and amplitude G k /G 0 . Hence, the interference effects can be understood as interference between the correlation functions of the elementary excitations within the pulse. We refer to this as interference of the second kind in contrast to conventional interference, where it is the wave function that interferes. We stress that while the wave function of a single particle can interfere with itself only, the correlation functions of various single particles can interfere with each other. Specifically, the Fraunhofer-like diffraction pattern in Fig. 3 can be captured by the expression  which is valid for ∆τ = 0 andn being an integer. In addition to the interference oscillations, we see that the overall decay of the visibility in Fig. 3 forn > 1 is similar to the casen = 1. Generally, the visibility decreases as a function of the difference of travel times, but the decay is faster for larger values ofn. Thus, the interference caused by the relative phases tend to reduce the visibility. From Fig. 3 we also conclude that to observe interference, the difference of travel times, ∆τ, and the width, 2Γ, should be of the same order. For ∆τ Γ or ∆τ Γ, there are no oscillations in the visibility.

Diffraction grid
Having established that the interference of the elementary components in a single pulse is revealed in the visibility as a Fraunhofer-like diffraction pattern, we now consider two separated pulses with arbitrary charge. In the case with only a single pulse, the interference pattern was fully determined by the charge of the pulse and the difference of the arm lengths. Now, the time delay w between the two pulses, see Eq. (18), offers another way to control the interference pattern. Figure 4 shows the effect of changing the delay w between two pulses with chargen = 2 (n = 1) in the left (right) panel. For w = 0, the pulses combine into a single pulse with the charge given by the total charge of the pulses. As the time delay in increased, the pulse splits into the two separate pulses. As we further decompose the multi-charged pulse into its elementary parts withn = 1, the number of oscillations decreases and the visibility is enhanced as can be seen by comparing the left and right panels. We also note that if we take w Γ, we completely split the pulses and they do not interfere. In Fig. 4, the oscillations are concentrated near w = ∆τ = 2Γ, and we should generally choose w close to those values to observe the clearest interference effects.
Under these conditions, we can further enhance the interference effects by increasing both the number of pulses m and their chargen. This leads to an increased number of oscillations as seen in Fig. 5, where we consider the splitting of a single pulse withn = 25 into five pulses withn = 5. To maximize the Fraunhofer-like pattern, the pulses are separated by the same time delay w. The magnitude of the visibility and its oscillations are reduced by the increased interference. Again, the visibility of a single pulse is given by its width Γ and its total chargen. Figure 5 shows how the visibility starts to rapidly oscillate with the time delay w as we decompose the large pulse into smaller ones. This behavior resembles the diffraction pattern generated by a regular spatial grid, where the smaller pulses here are acting like the slits of the grid, which are spaced by the delay time w. Thus, by tuning the separation between the pulses, one can directly control the resulting diffraction pattern.

Conclusions
We have theoretically investigated multi-particle interference in an electronic Mach-Zehnder interferometer driven by dynamic voltage pulses. To this end, we have described a Floquet scattering theory, which allows us to calculate the time-dependent current and the excess correlation function in a mesoscopic conductor driven by a time-dependent voltage. The current can be divided into a classical part and an interference term, with the classical contribution being independent of the temperature, while the interference term is gradually washed out by an increasing temperature. Moreover, the interference term determines the visibility that we have investigated for different pulse types.
For a single Lorentzian-shaped voltage pulse containing several charges, we have shown how oscillations in the visibility can be related to interferences between the individual charges making up the pulse. These interferences give rise to a Fraunhofer-like diffraction pattern consisting of interference oscillations superimposed on a central peak structure. Additional features appear as several pulses are injected into the interferometer, making it possible for different pulses to arrive simultaneously at the second quantum point contact of the setup. In that case, we observe an interference pattern that resembles what one would get with a regular spatial grid, however, with the spacing between the slits of the grid replaced by the delay time between the pulses. These predictions may be observed in future experiments by injecting multi-particle pulses into an electronic Mach-Zehnder inteferometer.