Thilo Hahn

Diana Vaclavkova

Miroslav Bartos

Karol Nogajewski

Marek Potemski

Kenji Watanabe

Takashi Taniguchi

Pawel Machnikowski

Tilmann Kuhn

Jacek Kasprzak

Paweł Machnikowski

Daniel Wigger
email: daniel.wigger@pwr.edu.pl

Destructive photon echo formation in six-wave mixing signals of a MoSe 2 monolayer

Introduction

The spin echo [1] is an essential effect in nuclear magnetic resonance spectroscopy and the basis for all sorts of complex radio pulse sequences [2] that are routinely applied in medicine [START_REF] Hendee | Physics and applications of medical imaging[END_REF], chemistry [START_REF] Hore | Nuclear magnetic resonance[END_REF], or physics [START_REF] Žutić | Spintronics: Fundamentals and applications[END_REF][START_REF] Weichselbaumer | Echo trains in pulsed electron spin resonance of a strongly coupled spin ensemble[END_REF][START_REF] Debnath | Selfstimulated pulse echo trains from inhomogeneously broadened spin ensembles[END_REF]. While spin resonances are driven by radio frequencies, optical frequencies are required to resonantly excite typical charge transitions; in this regime the analogous phenomenon is called photon echo [START_REF] Kurnit | Observation of a photon echo[END_REF]. First demonstrations were performed on ruby crystals [START_REF] Kurnit | Observation of a photon echo[END_REF][START_REF] Abella | Photon echoes[END_REF], but the photon echo spectroscopy also has a long standing history in semiconductor optics. It has been used to study different types of exciton dynamics, ranging from exciton-exciton scattering [START_REF] Schultheis | Ultrafast phase relaxation of excitons via exciton-exciton and exciton-electron collisions[END_REF][START_REF] Wegener | Line shape of time-resolved four-wave mixing[END_REF][START_REF] Leo | Effects of coherent polarization interactions on time-resolved degenerate four-wave mixing[END_REF][START_REF] Lindberg | Theory of the semiconductor photon echo[END_REF] to exciton-phonon coupling [START_REF] Carter | Echo peak-shift spectroscopy of non-markovian exciton dynamics in quantum wells[END_REF][START_REF] Jakubczyk | Impact of phonons on dephasing of individual excitons in deterministic quantum dot microlenses[END_REF], and it has been applied to three dimensional bulk [START_REF] Schultheis | Ultrafast phase relaxation of excitons via exciton-exciton and exciton-electron collisions[END_REF][START_REF] Becker | Femtosecond photon echoes from bandto-band transitions in GaAs[END_REF][START_REF] Lohner | Coherent optical polarization of bulk GaAs studied by femtosecond photon-echo spectroscopy[END_REF], two dimensional quantum well [START_REF] Webb | Observation of time-resolved picosecond stimulated photon echoes and free polarization decay in GaAs/AlGaAs multiple quantum wells[END_REF][START_REF] Langer | Access to long-term optical memories using photon echoes retrieved from semiconductor spins[END_REF], one dimensional nanowire [START_REF] Wagner | Population dynamics and dephasing of excitons and electron-hole pairs in polytype wurtzite/zincblende inp nanowires[END_REF], and zero dimensional quantum dot structures [START_REF] Langbein | Microscopic measurement of photon echo formation in groups of individual excitonic transitions[END_REF][START_REF] Moody | Fifth-order nonlinear optical response of excitonic states in an InAs quantum dot ensemble measured with two-dimensional spectroscopy[END_REF][START_REF] Delmonte | Coherent coupling of individual quantum dots measured with phase-referenced twodimensional spectroscopy: Photon echo versus double quantum coherence[END_REF][START_REF] Kosarev | Accurate photon echo timing by optical freezing of exciton dephasing and rephasing in quantum dots[END_REF]. However, the technique is not restricted to solid state samples; it has also been applied to liquids [START_REF] Vöhringer | Time-gated photon echo spectroscopy in liquids[END_REF][START_REF] De Boeij | Ultrafast solvation dynamics explored by femtosecond photon echo spectroscopies[END_REF][START_REF] Jordanides | Solvation dynamics in protein environments studied by photon echo spectroscopy[END_REF]. In its classical form the photon echo is based on a four-wave mixing (FWM) process and therefore it constitutes a nonlinear process of third order (χ (3) ) in the low excitation limit [START_REF] Cho | Photon echoes and related four-wave-mixing spectroscopies using phase-locked pulses[END_REF]. Especially in single low-dimensional systems like quantum dots it requires large effort to detect such nonlinear optical signals [START_REF] Cundiff | Coherent spectroscopy of semiconductors[END_REF][START_REF] Langbein | Coherent optical spectroscopy of semiconductor nanostructures[END_REF]. Interestingly, due to their strong excitonic optical response, monolayers of transition metal dichalcogenides (TMDCs) show a remarkable signal strength in FWM spectroscopy [START_REF] Moody | Intrinsic homogeneous linewidth and broadening mechanisms of excitons in monolayer transition metal dichalcogenides[END_REF][START_REF] Jakubczyk | Radiatively limited dephasing and exciton dynamics in MoSe 2 monolayers revealed with fourwave mixing microscopy[END_REF][START_REF] Hao | Neutral and charged inter-valley biexcitons in monolayer MoSe 2[END_REF][START_REF] Jakubczyk | Coherence and density dynamics of excitons in a single-layer MoS 2 reaching the homogeneous limit[END_REF]. So far the direct correspondence between inhomogeneous spectral broadening and photon echo duration has been used to map the inhomogeneity of TMDC monolayers [START_REF] Moody | Intrinsic homogeneous linewidth and broadening mechanisms of excitons in monolayer transition metal dichalcogenides[END_REF][START_REF] Jakubczyk | Radiatively limited dephasing and exciton dynamics in MoSe 2 monolayers revealed with fourwave mixing microscopy[END_REF][START_REF] Jakubczyk | Coherence and density dynamics of excitons in a single-layer MoS 2 reaching the homogeneous limit[END_REF][START_REF] Boule | Coherent dynamics and mapping of excitons in single-layer MoSe 2 and WSe 2 at the homogeneous limit[END_REF].

We here exploit the strong excitonic optical response of TMDC monolayers, and explore six-wave mixing (SWM) signals from a MoSe 2 monolayer in the low excitation limit, which here represents the χ (5) -regime [START_REF] Tominaga | Fifth order optical response of liquid CS 2 observed by ultrafast nonresonant six-wave mixing[END_REF][START_REF] Bolton | Demonstration of sixth-order coulomb correlations in a semiconductor single quantum well[END_REF][START_REF] Voss | Biexcitonic effects in the coherent control of the excitonic polarization detected in six-wave-mixing signals[END_REF][START_REF] Zhang | Opening four-wave mixing and six-wave mixing channels via dual electromagnetically induced transparency windows[END_REF]. We find that the signal dynamics exhibit a characteristic temporary suppression depending on the delay between the pulses. Supported by a theoretical model based on the local field effect describing exciton-exciton interaction, we explain that this suppression can be understood as the formation of a destructive photon echo.

Results

Sample and four-wave mixing dynamics We perform multiwave mixing experiments on an hBN/MoSe 2 /hBN heterostructure as schematically depicted in Fig. 1. By radio frequency modulating the incoming beams, the different pulses are labeled by phases φ j [START_REF] Langbein | Heterodyne spectral interferometry for multidimensional nonlinear spectroscopy of individual quantum systems[END_REF]. After heterodyning the emitted light with a specific N -wave mixing (NWM) phase, that is a particular phase combination of the form φ NWM = j a j φ j (a j ∈ Z, j |a j | = N -1, j a j = 1), we retrieve the corresponding nonlinear signal from the investigated monolayer [START_REF] Mukamel | Principles of nonlinear optical spectroscopy[END_REF]. Details are given in the experimental section. We use the same sample as investigated in Ref. [START_REF] Boule | Coherent dynamics and mapping of excitons in single-layer MoSe 2 and WSe 2 at the homogeneous limit[END_REF] where the echo formation in two-pulse FWM signals with φ FWM = 2φ 2φ 1 was used to study the inhomogeneity of the structure. The well known photon echo appears due to the dephasing impact of the structural inhomogeneity on the exciton's coherence as schematically depicted via Bloch vectors in Fig. 2 [START_REF] Shah | Ultrafast spectroscopy of semiconductors and semiconductor nanostructures[END_REF]. For the sake of simplicity in the illustration we show a combination of a π/2 and a π pulse, which results in a photon echo in the full coherence. However, when considering the FWM coherence characterized by the phase φ FWM the photon echo appears for any combination of pulse areas. The first laser pulse in Fig. 2(a) having a pulse area of θ 1 = π/2 and arriving at the time t = -τ generates an exciton coherence. Because of the presence of different transition energies, originating from strain and dielectric variations, the coherences generated by the first laser pulse oscillate with different frequencies resulting in a dephasing of the total coherence. As depicted in Fig. 2(b), after a delay τ , i.e., at the time t = 0, a second laser pulse with the pulse area θ 2 = π inverts all Bloch vectors. This is followed by a rephasing of the different coherences. The rephasing takes the same time as the dephasing and therefore the FWM signal is significantly enhanced due to constructive interferences at the time precisely given by the delay time t = τ between the two pulses.

In Fig. 3(a) we show the measured FWM signal from our sample as a function of the real time t after the second pulse and the delay τ . The pulse alignment is such that the pulse with phase φ 1 arrives at t = -τ and the pulse with φ 2 at t = 0. We find that the amplitude is not concentrated along the diagonal, which would represent a photon echo. Instead, starting from its maximum around (τ, t) = (0, 0) it basically decays to positive delays and times. This demonstrates that the investigated sample position is virtually homogenous because the photon echo is not present.

dp dt = i(1 -2n)[Ω(t) + V p] -(β + iω 0 )p , (1a)
dn dt = 2Im[Ω * (t)p] -Γn -λ(n -n ) , (1b)
dn dt = -Γn -λ(n -n) . (1c)
Here, n and n are the occupations in the optically coupled and uncoupled valleys, respectively, Ω(t) is the time dependent Rabi frequency describing the optical excitation by cocirclarly polarized pulses, β and Γ are the dephasing and the decay rate, respectively, ω 0 is the exciton energy in the absence of the local field coupling, and V is the strength of the local field coupling. Compared to Ref. [START_REF] Rodek | Local field effects in ultrafast lightmatter interaction measured by pump-probe spectroscopy of monolayer MoSe 2[END_REF] we do not find a significant impact of excitation induced dephasing as discussed in more detail in the supporting information (SI). As recently studied in Ref. [START_REF] Rodek | Local field effects in ultrafast lightmatter interaction measured by pump-probe spectroscopy of monolayer MoSe 2[END_REF] we additionally take an intervalley scattering contribution with the rate λ into account. In the special case λ = 0, the system reduces to a two-level system. We have recently investigated the local field model in the context of nonlinear optical signals focussing on FWM [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF] and pumpprobe spectroscopy [START_REF] Rodek | Local field effects in ultrafast lightmatter interaction measured by pump-probe spectroscopy of monolayer MoSe 2[END_REF], showing that the handy description produces resonant optical spectroscopy signals that are consistent with experiment. In a nutshell, the local field effect leads to energy shifts of the exciton depending on the its occupation.

In the limit of ultrafast optical pulses Eqs. ( 1) can be solved analytically relating the coherence p + and the occupation n + after the pulse to the respective values p -and n -immediately before the pulse

p + = p -cos 2 θ 2 + i 2 sin(θ)(1 -2n -)e iφ + sin 2 θ 2 p - * e i2φ ≈ p -+ i θ 2 (1 -2n -)e iφ , (2a)
n + = n -+ sin 2 θ 2 (1 -2n -) + sin(θ)Im p -e -iφ ≈ n -+ θ Im p -e -iφ , ( 2b
)
where θ is the pulse area and φ its phase. The approximations in Eqs. (2) describe the light-field induced contributions of a single pulse in first order of the pulse area, which is sufficient for reproducing the contributions of the order V 2 to the SWM signal, discussed in the main text of this work (other contributions, involving terms up to the second order in the pulse areas, are discussed in the SI.

Starting from the excitonic ground state characterized by n = n - 1 = 0, p = p - 1 = 0 the first pulse in the linear order creates

p + 1 ≈ i θ 1 2 e iφ1 , (3a)
n + 1 ≈ 0 , (3b)
which shows that relevant exciton occupations will only be created by second laser pulse from p + 1 . Our goal is to derive the nonlinear optical response in the lowest order of the pulse area θ because the experiments are performed with low pulse powers. In the case of the SWM signal this is the fifth order O(θ 5 ) (χ (5) -regime).

Once optical pulses have generated coherences p 0 and occupations n 0 the corresponding free propagation [Ω = 0, in Eqs. (1)] is governed by pure dephasing of p, exciton decay and intervalley scattering of n ( ) , and local field coupling between p and n. The corresponding dynamics can also be calculated analytically. Focusing on the optically addressed occupation first, its time-dependence reads

n(t) = n 0 + n 0 2 e -Γt + n 0 -n 0 2 e -(Γ+2λ)t , (4a)
which leads to a balanced occupation between the two valleys n and n on the timescale 1/(2λ) and a decay of both occupations with the rate Γ. As known from literature [START_REF] Rodek | Local field effects in ultrafast lightmatter interaction measured by pump-probe spectroscopy of monolayer MoSe 2[END_REF][START_REF] Schmidt | Ultrafast coulomb-induced intervalley coupling in atomically thin WS 2[END_REF][START_REF] Wang | Colloquium: Excitons in atomically thin transition metal dichalcogenides[END_REF] and as considered in this work the intervalley scattering is typically significantly faster than the decay, i.e., Γ 2λ. As a simplifying approximation for the sake of interpreting the results we can therefore assume that the occupations are balanced rapidly after an optical pulse, i.e.,

n(t) ≈ n 0 + n 0 2 e -Γt . (4b)
This step might lead to slight deviations for short delays in the range τ ≈ 1/(2λ) which are however hardly visible for the parameters chosen here as shown in the SI. Further we can only optically manipulate the occupation n, while n remains unchanged by the applied laser pulses. According to Eq. (3a) these pulses moreover add phase labels φ. As explained at the beginning of this section and as practically applied below, we are only interested in specific phase combinations that describe the considered wave-mixing signal. Consequently, any change of the occupation n is labeled by phase factors which do not apply to the other valley n . Therefore, the latter is irrelevant for the final optical signal and we can consider n 0 = 0 in Eq. (4b) resulting in

n(t) ≈ n 0 2 e -Γt . (4c)
We want to remark that in the opposite limit of a very slow scattering rate 2λ Γ the occupation dynamics would directly be given by n(t) = n 0 e -Γt . In this case all later discussions would work in the same way. By replacing V → 2V in all following derivations one can even directly retrieve the corresponding equations.

Based on this approximation for the occupation dynam-ics we can calculate the coherence dynamics

p(t) = p 0 exp   -i 2V Γ t 0 dt n(t )   e (iV -β)t ≈ p 0 exp -i 2V Γ n 0 2 1 -e -Γt e (iV -β)t ≈ p 0 e iV n0t e (iV -β)t ≈ p 0 1 -iV n 0 t - 1 2 (V n 0 t) 2 e (iV -β)t . (5)
In the first approximation step we have used the approximated occuption from Eq. (4c). As the exciton decay is much slower than the dephasing Γ β, in the second step we take Γ → 0. Finally, we keep the local field-induced contributions up to the second order in the exciton occupation, i.e., considering

V n 0 t
1. Below we will see that these are the terms which contribute to the SWM signal in the χ (5) -regime. In the SI we show that only terms up to O(V 2 ) contribute to the χ (5) -regime of the SWM signal, while contributions with higher powers in V appear in higher orders of the optical field. The last equation tells us that the local field induced mixing of the occupation n 0 with the coherence p 0 into the coherence p(t) in the lowest order grows linearly in time with a rate given by n 0 and the local field strength V . We will use this argument later on.

We simulate the FWM signal with the phase combination 2φ 2φ 1 by calculating the coherence dynamics p(t) following a two-pulse sequence and filter this quantity with respect to the required phase combination as described in Refs. [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF][START_REF] Rodek | Local field effects in ultrafast lightmatter interaction measured by pump-probe spectroscopy of monolayer MoSe 2[END_REF]. With this we find the signal dynamics |S FWM | ∼ 2φ2-φ1 p 2 (t, τ ) depicted in Fig. 3(b). To achieve the excellent agreement with the experiment in Fig. 3(a) we used a local field strength of V = 100 ps -1 and pulse areas of θ 1 = θ 2 /2 = θ = 0.02π, a Gaussian pulse duration of ∆t = 70 fs (standard deviation), a dephasing rate of β = 3 ps -1 , an intervalley scattering rate λ = 4 ps -1 , and a decay rate of Γ = 0.6 ps -1 . Note, that the depicted signal was calculated numerically because we considered a non-vanishing pulse duration and we did not employ the approximations mentioned in Eqs. ( 2) - [START_REF] Žutić | Spintronics: Fundamentals and applications[END_REF]. To set the strength of the local field coupling into context, typical values for GaAs quantum wells were reported in the range of a few meV [START_REF] Kim | Unusually slow temporal evolution of femtosecond four-wave-mixing signals in intrinsic GaAs quantum wells: Direct evidence for the dominance of interaction effects[END_REF][START_REF] Mayer | Evidence of biexcitonic contributions to four-wave mixing in GaAs quantum wells[END_REF], which is at least one order of magnitude smaller than considered here.

Six-wave mixing dynamics

In the present study we go one step further in multi-wave mixing and consider one of the possible SWM signals generated by three laser pulses, namely φ SWM = 2φ 3 -2φ 2 + φ 1 . In principle there are two delays in this pulse sequence that could be varied but, as schematically shown in Fig. 1 (right), we set the delay between the first two pulses to τ 12 = 0 and only vary the second one τ 23 = τ . The impact of a non-vanishing τ 12 is discussed in the SI. In the pure two-level system this configuration probes the polarization, i.e., it contains the same information as the FWM signal discussed before as shown in the SI. The measured SWM dynamics as a function of the real time t after the third pulse and the delay τ are shown in Fig. 4(a). Here, the two pulses with φ 1 and φ 2 arrive at t = -τ and the pulse with φ 3 at t = 0.

The signal consists of a strong maximum at small t ≈ 0.5 ps and τ ≈ 0. Moving to negative delays τ < 0 (pulse 3 is arriving before 1 and 2), the signal is strongly damped. For positive delays τ > 0 it decays much slower on the same timescale as the FWM signal in Fig. 3. We find a remarkable depression of the signal that stretches along a curved diagonal given by Eq. (9) (dashed black line), as will be derived below. In correspondence with the previously described constructive signal enhancement in the photon echo, we call this pronounced signal reduction a destructive photon echo. Later, we will discuss criteria justifying the labeling of this feature as an echo effect.

To identify the origin of this peculiar dynamical feature we model the SWM signal within the local field model described above. We consider the same system parameters as for the FWM signal but choose equal pulse area for all three pulses

θ 1 = θ 2 = θ 3 = θ = 0.02π, in agreement with the experiment.
In the SI it is shown that the exact choice of θ in the low excitation regime does not change the SWM signal dynamics. The SWM signal is extracted via the phase combination 2φ 3 -2φ 2 + φ 1 and the signal dynamics |S SWM | ∼ 2φ3-2φ2+φ1 p 3 (t, τ ) are depicted in Fig. 4(b). The retrieved signal agrees very well with the respective experiment in Fig. 4(a) showing the same characteristic suppression of the signal, i.e., the destructive photon echo. We give an overview regarding the impact of the different parameters on the SWM signal in the SI. Most importantly, it is shown in the SI, that the SWM dynamics do not exhibit the suppression for small local field strengths V demonstrating that this feature is a result of exciton-exciton interaction. Other specific features in the signal's dynamics ) and occupations (n, red). Green waved arrows show pulse induced, violet dotted ones local field induced, black ones free dynamics without, and blue ones with dephasing. The phase differences are defined as ∆ nm = φ nφ m . The flow chart holds for any order of the optical field. However, in the linear response regime considered here it is

0 n + 1 = 0 n + 2 = 0 n - 3 = 0.
like the asymmetric decay between positive and negative delays or between τ and t were studied in detail in Ref. [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF] and behave similarly in SWM. The destructive photon echo effect should also be present in an inhomogeneously broadened system, where also a traditional constructive photon echo develops. As the constructive echo selects only a specific time interval of the emitted signal around t = τ , it leads to a suppression of the entire SWM signal for all other times. Consequently, the SWM signal and therefore the destructive photon echo, which bends away from the diagonal (discussed later), is only visible in the vicinity of t = τ . This aspect is discussed in more detail in the SI.

Note, that the simulation shown in Fig. 4(b) takes the non-vanishing pulse duration into account and is therefore performed numerically. In the limit of ultrafast laser pulses we can find analytical expressions for the SWM signal. Given that the experiment is carried out with weak pulse powers, we restrict the following studies on the lowest order in the optical field which is the χ (5) -regime. In this order we have already eight different contributions as derived in the SI. From these we will focus on the ones with the strongest local field contribution which is O(V 2 ), i.e., we omit all terms O(V 1 ) and O(V 0 ). The reason for this is the absence of the destructive echo for small V . In Ref. [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF] we have derived a flow chart representation for the construction of nonlinear wave mixing signals. In Fig. 5 we employ this procedure to disentangle the origin of the different contributions to the SWM signal. The flow chart only shows coherences p (blue) and occupations n (red) with phase combinations (green, given as left indices) relevant for the final signal; corresponding flow charts for the contributions with O(V 1 ) and O(V 0 ) are provided in the SI. Conveniently, the O(V 2 )-contributions can be derived with the appoximations given in Eq. ( 2). Note, that we introduce phase differences ∆ nm = φ nφ m here. The right lower index refers to the pulse number, while the upper -(+) indicates times immediately before (after) that pulse. Starting from the excitonic ground state with n = 0, p = 0 the first pulse generates the occupation 0 n + 1 and the coherence φ1 p + 1 . Note that in the scheme we do not restrict ourselves to the lowest order contributions in the light field and thus include also the occupation which is of second order in the pulse amplitude. The second pulse arrives at the same time (τ 12 = 0) and creates two occupations 0 n + 2 and |∆21| n + 2 and the coherence φ2 p + 2 . During the following propagation for the time of the delay τ two relevant things happen: On the one hand all coherences experience dephasing (blue arrows). On the other hand φ2 p + 2 is mixed with |∆21| n + 2 via the local field coupling resulting in the coherence 2φ2-φ1 p - 3 before the third pulse. Note, that this contribution carries the FWM phase φ FWM = 2φ 2φ 1 which we will come back to below. After the final third pulse we have five relevant terms 1 : The coherence φ3 p + 3 and the three occupations |∆32| n + 3 , |∆32-∆21| n + 3 , and |∆21| n + 3 . The polarization φ1 p + 1 is not affected by the second and third pulse and just evolves into φ1 p + 3 . During the remaining propagation step in real time t the three relevant SWM contributions are generated by local field mixing processes according to Eq. ( 5).

Considering the contribution on the right first, we have to mix φ1 p + 3 with |∆32| n + 3 twice. This results in the phase combination φ 1 + 2(φ 3φ 2 ) which is the SWM phase combination. According to Eq. ( 5) each of these local field mixing processes contributes with an amplitude of V t resulting in the amplitude (V t) 2 . In addition the amplitude is damped due to the dephasing happening during the delay. Following the two paths in the diagram back to this propagation step we find that φ1 p + 2 → φ1 p - 3 and φ2 p + 2 → φ2 p - 3 contribute with a dephas-1 Note, that the contribution to 2φ2-φ1 p - 3 generated directly by the first and second pulse is of lower order in V and therefore neglected here but considered in the SI. However, it is fully included in the numerical solution.

ing term ∼ e -βτ . The latter one is used twice due to the double local field mixing resulting in the total damping rate of e -3βτ .

Moving on to the left contribution in Fig. 5 we find that in the last propagation the coherence φ3 p + 3 is local-field mixed with |∆32| n + 3 and |∆21| n + 3 once, resulting in the SWM phase combination φ 3 + (φ 3φ 2 ) -(φ 2φ 1 ). Each of these processes contributes an amplitude of V t again resulting in (V t) 2 . The difference to the first term is that here only |∆32| n + 3 stems from a coherence, namely φ2 p + 2 which experiences a dephasing during the delay. Consequently, the entire contribution is damped by e -βτ . This already shows that this contribution is more important for the final SWM signal than the right term discussed first.

The final main contribution is the middle one in Fig. 5. Here, in the last real time propagation the local field mixing happens between φ3 p + 3 and |∆32-∆21| n + 3 which contributes with an amplitude of V t. Following the path back, we find that the occupation is created from the polarization 2φ2-φ1 p - 3 which itself is produced by a local field mixing step between φ2 p + 2 and |∆21| n + 2 . The mixing process lasts for the delay time and therefore contributes with a factor V τ to the final SWM amplitude. This is also the step where the only dephasing happens resulting in a final amplitude of V 2 τ te -βτ . Directly comparing the three contributions in Fig. 5 we find that the right one is smaller than the other two due to the stronger dephasing during the delay. We will therefore disregard this term from now on. The other two terms are of the order V 2 and exhibit the same dephasing with e -βτ . From the full derivation given in the SI we find that these two terms in Fig. 5 carry opposite signs. We will explain these signs later when introducing an effective Bloch vector picture to illustrate the destructive echo formation. Finally, adding the two terms we get

SWM p 3 (t, τ ) ≈ V 2 (τ t -t 2 )e -β(τ +t) , (6)
where we also added the dephasing rate for the propagation in t after the last pulse. We find that the two terms exactly compensate each other for t = τ . This renders our first step towards the understanding of the destructive echo. We identified two paths (signal contributions) with the same damping but with different magnitudes depending on the delay τ that act destructively in the total SWM signal. To illustrate the interplay between these two terms in Fig. 6 we schematically plot the applied laser pulses and the absolute value of the final SWM signal retrieved from the previous derivations. We also include the intermediate growth of the FWM coherence lasting for the delay τ . The full SWM signal in the χ (5) -regime is derived in the SI and reads

SWM p 3 (t, τ ) = θ 2 5 i + V (τ -4t) - i 2 (V t) 2 e -2βτ + iV 2 t(τ -t) e -β(t+τ ) . (7)
Interestingly, we also find a suppression of the signal when only considering the terms O(V 1 ), namely for t = τ /4. We obviously do not find such a feature in our measurement, which shows that the linear oder in V does not have a significant contribution. One advantage of the spectral interferometry in the applied approach is the possibility to detect -besides the amplitude -also the phase of the SWM signal [START_REF] Kasprzak | Coherent coupling between distant excitons revealed by twodimensional nonlinear hyperspectral imaging[END_REF]. We see that for t < τ the positive contribution ∼ τ t dominates while for t > τ the negative one ∼ -t 2 is larger. Therefore we expect a phase jump when crossing the destructive echo in time t < τ → t > τ . To confirm this in Fig. 7(a) we plot the measured (solid blue) and calculated (dashed red) SWM signal amplitude as a function of time t for τ = 0.35 ps. We slightly adjusted the pulse duration to ∆t = 80 fs to achieve the good agreement with this experiment. The finding that the simulated signal does not drop to zero shows the impact of the noncompensating contributions ∼ V 2 , like the ordinary SWM signal in Eq. ( 7), and the influence of the considered non-vanishing pulse duration. In Fig. 7(b) we show the corresponding SWM phase from S SWM = |S SWM | e iϕ . We indeed find a jump of approximately π at the destructive echo as would be expected for the two dominant terms. However, all other contributions, which are naturally present in the numerical simulation, lead to a reduction of the phase jump and a further distortion of the destructive echo dynamics. We again find that the full depression of the signal happens at t ≈ 0.4 ps which is slightly later than t = τ = 0.35 ps. This slight delay of the effect mainly stems from the non-vanishing decay rate of the exciton Γ = 0.

To include the exciton decay we have to go back to the full dynamics in Eq. ( 5) and directly consider V n 0 1e -Γt /Γ 1 without taking the limit Γ → 0 leading to

n(t) = n 0 e -Γt , (8a)
p(t) ≈ p 0 1 -i V Γ n 0 1 -e -Γt - V 2 2Γ 2 n 2 0 1 -e -Γt 2 . ( 8b
)
Now we consider the different paths in Fig. 5 to determine the decay's impact on the different SWM contributions. The term ∼ t 2 stems from the delay step |∆21| n + 2 → |∆21| n - 3 , contributing a factor e -Γτ , and the final double wave mixing step gives 1e -Γt 2 . The ∼ τ t contribution has a local-field mixing step during the delay leading to 1e -Γτ and only a single local field mixing after the last pulse, i.e., 1e -Γt . Finally, the two contributions compensate each other for

1 -e -Γt 2 e -Γτ = 1 -e -Γt 1 -e -Γτ , ⇒ τ = 1 Γ ln 2 -e -Γt . (9)
We find that the decay additionally acts on the local field mixing processes because it dynamically reduces the occupation during the mixing. This affects the t 2 -contribution in a different way than the part ∼ τ t and the destructive echo gets delayed with respect to t = τ . Equation ( 9) does not depend on V as confirmed in the SI. Consequently, it also holds in both limiting cases of small and large intervalley scattering as discussed below Eq. (4c) and can be expected to be correct also for intermediate λ-values. Note, that the curve describing the dynamics of the SWM signal depression in Fig. 4 follows Eq. ( 9) and is additionally shifted by the pulse width ∆t to shorter times to compensate the non-vanishing pulse duration. The curve almost perfectly follows the distribution of the destructive echo.

Quasi-Bloch vector and phase inhomogeneity As discussed in the previous section, the FWM-polarization constitutes the basis for the SWM signal. Therefore, it is instructive to introduce the concept of the quasi-Bloch vector to explain the wavemixing origin for the FWM signal in the next section. Building on this picture we will be able to explain the destructive photon echo formation in the SWM signal in the following section.

Four-wave mixing quasi-Bloch vector

The ordinary echo in FWM is most instructively visualized within the Bloch vector picture in Fig. 2 where the inhomogeneity leads to different rotation speeds of the various Bloch vector realizations. In that case each exciton energy is represented by a single Bloch vector. Here, we are dealing with a slightly different situation. The wave mixing experiment is repeated numerous times with successively different phase combinations (φ 1 , φ 2 , φ 3 ) for each repetition of the measurement. Therefore we have to represent each run of the experiment, i.e., each phase combination, by a single Bloch vector. To find the realized Bloch vectors we have a look at the coherence and occupation immediately after the second pulse [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF] p

+ 2 = i θ 2 e iφ1 + e iφ2 -θ 2 e i(2φ2-φ1) + O(θ 5 ) , (10a)
n + 2 = θ 2 2 [1 + cos(φ 2 -φ 1 )] + O(θ 4 ) , (10b)
where we consider all terms up to the third order in the optical field (χ (3) -regime) because our measurements are performed with low excitation powers. Following the flow chart in Fig. 5 we find that the two terms φ2 p + 2 and |∆21| n + 2 perform local field induced mixing during the following propagation before the third pulse (delay step). Exactly this term was identified as a contribution of the FWM signal in Ref. [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF]. Therefore, we will study this process in more detail. Our goal is to isolate a set of Bloch vectors that allows us to directly extract the final wave mixing signal by integrating over the entire set. To achieve this we have to already filter the coherence with the respective wave mixing phase factor. In this case we are interested in the FWM contribution, which means that we have to consider the phase factor e iφFWM = e i(2φ2-φ1) and consequently the filtered polarization φ2 p + 2 e -iφFWM . This brings us to the two relevant terms

p = φ2 p + 2 e -iφFWM = i θ 2 e -i(φ2-φ1) , (11a)
ñ = |∆21| n + 2 = θ 2 2 cos(φ 2 -φ 1 ) . (11b)
We use these expressions to define a quasi-Bloch vector for this FWM contribution via

ṽFWM (t = 0) =    Re(p) Im(p) ñ    = θ 2    sin(∆ 21 ) sin(∆ 21 + π/2) θ sin(∆ 21 + π/2)    , (12)
with the previously introduced phase difference ∆ 21 = φ 2 -φ 1 . We call this quantity quasi-Bloch vector because it does not represent the entire density matrix, as the full Bloch vector does. When varying ∆ 21 as it is done in experiment and simulation, we already see that this set of quasi-Bloch vectors follows a three dimensional Lissajous curve with a frequency relationship of 1:1:1 as depicted in Fig. 8(a). It appears as a circle tilted around the Re(p)-axis and its projections on the different planes of the coordinate system form Lissajous curves with 1:1 frequency relations. This results in a distribution of occupations ñ, whose spread, represented by the tilt angle, is given by the considered pulse areas.

R e ( p)

Im ( p) The final integrated FWM coherence FWM p is marked as blue cross and its propagation is given by V t.

Starting from this distribution of quasi-Bloch vectors directly after the second laser pulse their dynamics is governed by the local field mixing between p and ñ described by [see Eq. ( 5)]

dṽ dt = ṽ ×    0 0 V ñ   . (13)
These dynamics are a rotation of the coherence p depending on the occupation ñ. Particularly important is that the rotation flips its direction for opposite signs of ñ. At the same time ñ remains unaffected. Finally, from the distribution of quasi-Bloch vectors at a given time ṽ(t) we can directly retrieve the final wave mixing contribution by integrating over all possible phase combinations, i.e., over ∆ 21 . We directly see that the

For the initial distribution of quasi-Bloch vectors in Eq. ( 12) we find that the final FWM coherence vanishes because all vectors are equally distributed on the tilted circle in Fig. 8(a). However, the following dynamics lead to a nonvanishing FWM contribution. The rotation of p with frequencies proportional to ñ has two crucial consequences: (i) The presence of a distribution of ñ results in different frequencies. As already discussed in Ref. [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF] this means that each run of the experiment results in a slightly different emission energy of the exciton and consequently leads to a broadened FWM spectrum. Because the distribution of ñ stems from the variation of the applied laser pulse phase combinations we call this spectral broadening local field induced phase inhomogeneity. (ii) According to Eq. ( 13) the sign of ñ determines the rotation direction of the quasi-Bloch vectors. The initial distribution of quasi-Bloch vectors in Fig. 8(a) is a circle that is tilted around the Re(p)-axis and all points on this axis are not affected by the local field induced rotation (ñ = 0). This axis remains a symmetry axis also for the dynamics of the quasi-Bloch vectors: one half of the circle rotates clockwise, the other counterclockwise in exactly the same way. Consequently, when integrating over the relative phase ∆ 21 the opposing Im(p) contributions compensate each other and the final FWM coherence is real, i.e., Re(p) = 0. Because only Re(p) will later contribute to the FWM signal, the crucial quantity determining the propagation direction of Re(p) is ñ Im(p) as can be seen in Eq. ( 13). In Fig. 8(b) we show the quasi-Bloch vector dynamics in the Re(p),Im(p)-plane after the second pulse. The bright open circles represent the initial homogeneous distribution and the filled dark circles the distribution for a time t > 0. Looking at the initial distribution in Fig. 8(a) we find that the semicircle with Im(p) < 0 also has ñ < 0 and the other semicircle with Im(p) > 0 has ñ > 0. This means that the former rotate clockwise, while the latter rotate counter-clockwise as indicated by the curved black arrows in Fig. 8(b). Consequently, the weight of all quasi-Bloch vectors, which is depicted by the blue cross, moves into the positive Re(p)-direction. The speed of this movement is given by V and the tilt angle of the initial quasi-Bloch vector distribution. This illustrates that the FWM contribution grows in time and that the local field induced phase inhomogeneity governs this process.

Six-wave mixing quasi-Bloch vector

Equipped with the quasi-Bloch vector picture we can retrieve the two most important SWM signal contributions by carrying out an analogue discussion. According to Fig. 5, from all possible contributions after the third laser pulse, we only need the coherence φ3 p + 3 and the occupations |∆32| n + 3 , |∆32-∆21| n + 3 , and |∆21| n + 3 . Omitting all other terms we get

p → i θ 2 e iφ3 , ( 15a
)
n → θ 2 2 cos(∆ 32 ) + θ 2 2 cos(∆ 21 ) - θ 4 8 V τ sin(∆ 32 -∆ 21 ) . ( 15b
)
To generate the SWM contribution ∼ τ t we need to local-field mix the coherence filtered with respect to the SWM phase combination 2φ 3 -2φ 2 + φ 1 with the last term of the occupation in Eq. (15b), i.e.,

φ3 p + 3 e -iφSWM = i θ 2 e -i(∆32-∆21) , ( 16a
) |∆32-∆21| n + 3 = - θ 4 8 V τ sin(∆ 32 -∆ 21 ) . (16b)
From this we can directly read the respective initial quasi-Bloch vector

ṽSWM,τt (t = 0) = θ 2    sin(α) sin(α + π/2) -1 4 θ 3 V τ sin(α)    , (17)
with α = ∆ 32 -∆ 21 . We see that this is the same Lissajous curve structure as in the previously discussed FWM case. The only differences are that the phase variation is now given by ∆ 32 -∆ 21 , which still uniformly covers all angles, a tilt around the Im(p)-axis, and a tilt angle such that the occupation distribution ñ has an amplitude of -θ 4 V τ /8. The subsequent free dynamics are again governed by Eq. ( 13).

Before discussing the shape and dynamics of this quasi-Bloch vector set we consider the other relevant SWM contribution ∼ t 2 . Therefore, we have a look at

φ3 p + 3 e -iφSWM = i θ 2 e -i(∆32-∆21) , ( 18a
) |∆32| n + 3 + |∆21| n + 3 = θ 2 cos ∆ 32 -∆ 21 2 × cos ∆ 32 + ∆ 21 2 . ( 18b
)
We are finally only interested in the coherence SWM p 3 describing the SWM signal. So, at this point we can already R e ( p)

Im ( p) ṽSWM,t 2 ñ ṽSWM,τt
d(∆ 32 + ∆ 21 )θ 2 cos ∆ 32 -∆ 21 2 cos ∆ 32 + ∆ 21 2 = 4θ 2 cos ∆ 32 -∆ 21 2 . ( 19
)
This leaves us with the quasi-Bloch vector for this SWM contribution reading

ṽSWM,t 2 (t = 0) = θ 2    sin(α) sin(α + π/2) 8θ sin (α/2 + π/2)    , (20)
with α = ∆ 32 -∆ 21 , which is obviously a three dimensional Lissajous curve with the frequency ratio 2:2:1.

The two initial distributions of quasi-Bloch vectors given in Eqs. ( 17) and ( 20) are depicted in Fig. 9 in blue and violet, respectively. We vividly see the Lissajous curves leading to the distributions in ñ-direction. While the three dimensional perspective of the curves is appealing, we will consider the projections to the three different planes of the coordinate system when discussing the dynamics of the quasi-Bloch vectors, which are depicted in pale colors. These projections are shown in Fig. 10, where different times are sorted in rows (increasing from top to bottom) and different perspectives in columns.

First considering the ∼ τ t contribution (blue) we find essentially the same situation as for the intermediate FWM contribution discussed in Fig. 8. The only important difference is that for the SWM contribution the circle is tilted around the Im(p)-axis. Following the same arguments as before the quasi-Bloch vector realizations move into the positive Im(p)direction. This movement is shown in the second column from the right in Fig. 10. The initial homogeneous distribution is depicted by bright empty circles and for the considered times (written on the left) as filled dark circles. The movement is marked by the curved black arrow and the final SWM contribution after integration over all quasi-Bloch vector realizations as a dark blue cross. The speed of this movement is again constant and given by V τ , which is proportional to the spread of ñ and therefore given by the tilt angle of the circle.

For the ∼ t 2 contribution (violet) the situation is more involved because the Lissajous curve is less trivial. In this situation the crossing point of the two loops lies at ṽSWM,t 2 = (0, Im(p), 0). Because of the vanishing occupation this point will not move in the following propagation and the Im(p)-axis is again the symmetry axis of the dynamics. Consequently, the opposing Re(p) values compensate each other, such that the final SWM contribution is again purely imaginary. Therefore, we have to consider the movement into the direction governed by the product ñ Re(p). When looking at the left column in Fig. 10 the distribution reaches values with ñ Re(p) > 0 (large dots) and ñ Re(p) < 0 (small dots). Following Eq. ( 13) the large dots move into the negative Im(p)-direction while the small ones move into the positive Im(p)-direction as indicated by the black horizontal arrows in the Im(p),ñ-plane at t < τ . In the complex p-plane this leads to two counter movements marked as arrows in the right column. The small dots move up (dashed curved arrows) and the large dots move down (solid curved arrows). Initially both of the these movements start with the same velocity because large and small dots have the same absolute values |ñ Re(p)|, i.e., the Lissajous curve is symmetric. Therefore, the integrated SWM signal starts with a vanishing velocity for small t ≈ 0. However, during the following dynamics the quasi-Bloch vectors with positive Im(p)-velocity (small dots) rotate collectively in such a way that this part of the curve obtains smaller occupations. This can be seen by focussing on the bottom-left and upper-right corners of the panels in the left column. Accordingly, those with negative Im(p)velocity (large dots) obtain larger |ñ| as can be seen in the second column [Im(p),ñ-plane] marked by the horizontal black arrows. In the left Lissajous curve this means that the large dots move towards the top left and bottom right corner of the plot, while the small dots move inwards. Consequently, the negative Im(p)-velocity component gains speed because |ñ Re(p) > 0| (large dots) grows, while the positive component slows down as |ñ Re(p) < 0| (small dots) shrinks. In the complex p-plane (right column) this means that the upwards movement of the quasi-Bloch vectors represented by small dots (dashed arrows) slows down for increasing times t while the downwards movement of the large dots (solid arrows) speeds up. In summary, for the integrated SWM coherence, marked by the violet cross, we find an acceleration towards negative Im(p)-values given by

-2 -1 0 1 2 t = 0 Re(p), ñ Im(p), ñ Re(p), Im(p) ṽSWM,τt ṽSWM,t 2 -1 0 1 SWM p(τ t) SWM p(t 2 ) -2 -1 0 1 2 t < τ t < τ -1 0 1 t < τ 10× ×10 -2 -1 0 1 2 t = τ t = τ -1 0 1 t = τ 10× -2 -1 0 1 2 -1 0 1 t > τ -1 0 1 t > τ -1 0 1 -1 0 1 t > τ -1 0 1 ñ/(2θ 2 ) 2 Im(p)/θ ñ/(2θ 2 ) 2 Im(p)/θ ñ/(2θ 2 ) 2 Im(p)/θ ×10 ñ/( 2θ
-V 2 t 2 .
Exactly at t = τ in the third row the two integrated SWM contributions marked in the third and forth column have the same distance from the center as indicated by the dashed horizontal lines which show that they compensate each other. At this time the destructive echo appears due to the destructive

Conclusion

To conclude we again compare the traditional photon echo in FWM with the newly discovered destructive photon echo in SWM. As summarized in Tab. 1 we define three criteria that characterize the echo formation. Firstly, the timing of the photon echo formation is basically determined by the chosen delay τ between the two laser pulses. Slight deviations from the exact t = τ timing are used to learn about internal dynamics of the studied system [START_REF] Poltavtsev | Photon echo transients from an inhomogeneous ensemble of semiconductor quantum dots[END_REF]. The same holds for the destructive echo in SWM. As discussed in detail the two main contributions compensate each other exactly at t = τ . Nonetheless, in the detected and simulated full signal dynamics in Fig. 4 the depression happens at t > τ . The reason is that exciton decay dynamically changes the slope of the destructive echo. Secondly, the photon echo appears because all considered Bloch vectors form a constructive interference. To explain this effect only the conventional FWM signal has to be considered, which would appear as a single path in a flow chart like the one in Fig. 5 (see also Ref. [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF]). In the case of the destructive echo we have shown that we need two paths which interfere destructively to explain the novel echo effect. Finally, the fundamental source of the photon echo is any sort of inhomogeneity, which might appear in space by a locally varying exciton energy [START_REF] Koch | Subpicosecond photon-echo spectroscopy on GaAs/AlAs short-period superlattices[END_REF], or in time via external noise that acts on the optical transition energy [START_REF] Langbein | Coherent optical spectroscopy of semiconductor nanostructures[END_REF]. In the case of the destructive echo we have shown that the local field together with a variation of the applied laser pulse phases results in a spectral broadening. As discussed in Ref. [START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF] this local field induced phase inhomogeneity can already be detected in FWM spectra. This summary shows that the criteria defining the traditional photon echo can also be applied to the destructive echo effect. Therefore, we conclude that we really found a new photon echo effect.

In addition we have developed a quasi-Bloch vector picture to illustrate the generation of the destructive echo inspired by the instructive Bloch vector image of the photon echo. We found that the SWM-relevant quasi-Bloch vectors are distributed along Lissajous curves which adds surprising aesthetics to the involved Bloch vector dynamics in SWM.

Following this proof-of-principle demonstration of a new photon echo effect the natural next tasks will be to further explore the destructive photon echo's application possibilities. One obvious step will be to measure photon echos at various spots of the sample and thereby get a measure for the spatial distribution of the local field coupling strength. We know that different pulse sequences in FWM probe different observables [START_REF] Mermillod | Dynamics of excitons in individual InAs quantum dots revealed in four-wave mixing spectroscopy[END_REF]. Therefore we hope that novel pulse sequences in higher wave-mixing spectroscopy will be designed that should allow to access parameters governing the nonlinear lightinduced dynamics in semiconductors, like the strength of interaction mechanisms, that are not easily accessible with traditional techniques. These developments are particularly attractive for materials with a large nonlinearity promoting strong wave mixing responses like TMDCs.

Experimental Section

In the experiment, we employ a mode-locked laser (Tsunami Femto provided by Spectra Physics) generating femtoscond pulses at a repetition rate of 76 MHz. The wavelength is tuned to 752 nm, which corresponds to the A exciton transition in the investigated hBN/MoSe 2 /hBN heterostructure at T = 5 K. The sample is kept in an optical He-flow cryostat and microspectroscopy is performed using an Olympus microscope objective (LCPLN50XIR, numerical aperture of 0.65), which is installed on an XYZ piezo stage from Physik Instrumente. To perform the multi-pulse, heterodyne experiments, the initial pulse train from the laser source is split into three and each of these beams is focussed into a separate acousto-optic modulator (AOM). The AOMs are driven at distinct radio-frequencies of Ω 1 = 79, Ω 2 = 80.77, and Ω 3 = 81 MHz, such that the deflected beams acquire corresponding frequency upshifts. Next, the time delays between the first two (τ 12 ) and the last two beams (τ 23 ) are introduced by a pair of mechanical delay stages. In addition, a grating-based pulse shaper is employed to correct the temporal chirp of the pulses, when passing through optical elements of the setup, especially AOMs and the microscope objective. The such prepared pulse sequence is then recombined into the same spatial mode and focussed at the heterostructure sample reaching a diffraction limited spot diameter of around 0.8 µm. The nonlinear optical signal from the sample is retrieved in the back-reflectance geometry. To isolate the desired wave-mixing response another AOM is used, operating at the heterodyne frequency generated by a home-made analogue radio-frequency mixing electronics, assembled from individual components provided by Mini Circuits. For example, the studied SWM signal is detected at (Ω SWM = 2Ω 3 -2Ω 2 + Ω 1 )=81.54 MHz. After being deflected from the AOM the signal is frequency downshifted by Ω SWM and the unique wave-mixing component under consideration carries the original frequency from the laser source. The signal can now interfere with a reference beam from the same laser source, which propagates in the vicinity of the driving beams. At the same time all other optical response components present in the signal reflected from the sample are still modulated in the MHz range and thus average out completely during a single acquisition time of a few ms. Importantly, the mixing AOM operates in a Bragg configuration. Hence, simultaneously the reference beam gets deflected onto the signal, receiving a frequency upshift by Ω SWM . This again allows for interference with the signal, which was not deflected by the AOM (operating at Ω SWM ). In total we generate two beams in which the considered wave-mixing signal is interfering with the reference pulse, but with the opposite phase, as required by energy conservation. The spectrally-resolved interference is obtained with an imaging spectrometer (Princeton Instruments with a focal length of 750 mm) and detected with an CCD camera (PIXIS from Princeton Instruments, with an eXcelon coating). A background free, shot-noise limited detection is achieved by a balanced detection and multi-acquisition. Namely, we exploit the Bragg configuration and subtract the two mixed beams impinging at different positions on the CCD. In addition, the phase of the mixing AOM is cycled between 0 and π to overcome any classical noise from the CCD camera. The wave-mixing signal amplitude and phase are obtained from the interferogram by performing spectral interferometry.

Again, we have omitted all terms which will not contribute to the SWM signal in the χ (5) -regime. Note, that due to intervalley scattering, the occupation n - 3 is split between the valleys as discussed in the main text. The third pulse then transforms polarization and occupation into

p + 3 = i θ 2 e iφ1 e -βτ + e iφ2 e -βτ + e iφ3 -i θ 3 4 e i(φ3-φ2+φ1) - θ 3 8 e i(2φ2-φ1) (i -V τ ) e -βτ -i θ 3 8 e i(2φ3-φ1) + e i(2φ3-φ2) e -βτ + θ 5 32 e i(2φ3-2φ2+φ1) (i + V τ ) e -βτ , (S5a)
n + 3 = θ 2 2 cos(φ 2 -φ 1 ) + cos(φ 3 -φ 1 )e -βτ + cos(φ 3 -φ 2 )e -βτ - θ 4 8 cos(φ 3 -2φ 2 + φ 1 )e -βτ - θ 4 8 V τ sin(φ 3 -2φ 2 + φ 1 )e -βτ . (S5b)
Note, that we have again omitted all terms that are irrelevant for the final SWM signal, especially constant terms for the occupation that do not carry any phase information. Next, we consider the time-evolution after the third pulse, where we can already sort for the SWM phase 2φ 3 -2φ 2 +φ 1 . In Fig. 5 in the main text we schematically depict the origin of the contributions with two local-field mixing processes for O(V 2 ). We disentangle the different contributions with O(V 1 ) and the SWM polarization from the pure TLS without local-field mixing with O(V 0 ) by the help of the flow charts in Fig. S1.

Starting with the polarization in (S5a) that already carries the SWM phase, their signal contribution reads

SWM p/θ 5 = i 32 e -βτ + 1 32 V τ e -βτ . (S6a)
The first part represents the SWM signal from the pure TLS as it is independent of V . Its origin is depicted in Fig. S1 on the left side (i) and we directly see that the coherence is detected when scanning the delay between the second and third pulse as mentioned in the main text. The second term (ii) is proportional to V τ , therefore it stems from the local-field mixing between φ2 p + 2 and |∆21| n + 2 during the delay propagation, depicted is the most right path in Fig. S1.

From the FWM-polarizations created by the third pulse, we get the two SWM contributions (iii) and (iv) by a single local field mixing process in the real time propagation. According to Eq. ( 5) we have to calculate p(t) = -iV tp 0 n 0 for single and From the analytical solution, we find the destructive echo only in the highest possible order, i.e., V 2 . To understand, how the local field strength influences the destructive echo, in Fig. S3 we show the SWM dynamics for different values of V increasing from (a) to (d). While for V = 2 ps -1 (a) and V = 5 ps -1 (b) no destructive echo is visible, for V = 10 ps -1 (c), the depression is already visible. For comparison Fig. S3(d) shows the case with V = 100 ps -1 from the main text. This again confirms, that the V 2 -order has to dominate the signal to observe the destructive echo. For V = 10 ps -1 in Fig. S3(c) we see that the destructive echo follows the dynamics of the minimum approximated by Eq. ( 9) shown as dashed curved line. This confirms our finding that this behavior is not affected by the strength of the local field coupling V . In addition to the local field coupling which shifts the transition energy depending on the exciton occupation, we can also take excitation induced dephasing (EID) into account [1], where the dephasing grows with the occupation β total = β + W (n + n ). For an EID on the same order of magnitude as the local field the destructive echo effect is strongly suppressed as illustrated in Fig. S4( In Fig. S5 we plot the same SWM amplitude dynamics as in Fig. 7(a) in the main text but in the simulation consider a EID of W = 10 ps -1 , which is much smaller than the applied local field coupling of V = 100 ps -1 . We see that the destructive photon echo minimum is significantly less pronounced when taking the EID into account. Therefore, to maintain the good agreement between theory and experiment, we have excluded EID from the discussion. We can conclude that the EID should not have a significant impact in our system.

p(t) = -(V t) 2 p 0 n 2 0 /

S3 Impact of excitation induced dephasing

S4 Impact of intervalley scattering

In the context of the approximation in Eq. (4b) we estimated that the intervalley scattering should only have an impact for small delays. To confirm this we calculate the SWM amplitude dynamics for λ = 0 as depicted in Fig. S6 In the main text we have considered a vanishing delay between the two pulses labeled with φ 1 and φ 2 . As a reference the corresponding SWM dynamics are shown in Fig. S7(a) as functions of τ 23 . To identify the impact of a non-vanishing delay between these two pulses, in Fig. S7(b) we plot the SWM dynamics for τ 12 = 0.35 ps. As depicted in the inset, for this delay the two pulses are almost entirely separated. The corresponding SWM amplitude dynamics are identical to the ones with τ 12 = 0 in Fig. S7(a) but have an overall smaller amplitude. This demonstrates that the destructive echo is not caused by pulse overlap effects. The reason for the decreased amplitude is, that the coherence φ1 p + 1 after the first arriving pulse declines due to the considered dephasing. Consequently, the following processes, i.e., optical excitations and local-field mixing, are the same but start from a smaller coherence. Therefore, the final SWM signal is weaker.

Figure 1 :

 1 Figure 1: Sample structure consisting of a multilayer hBN, monolayer MoSe 2 , multilayer hBN stack. (Left) Four-wave mixing (FWM) generated by two laser pulses with tunable delay τ . (Right) Six-wave mixing (SWM) with three excitation pulses having a tunable delay τ between pulses 2 and 3. The FWM and SWM dynamics are measured in real time t.

Figure 2 :

 2 Figure 2: Schematic Bloch vector image of the photon echo process. (a) 1st laser pulse excitation with following dissipation. (b) 2nd excitation inverting all coherences and following rephasing. Laser pulse rotations are depicted in green, initial and final Bloch vectors in grey and blue, respectively, and the inhomogeneity induced coherence dynamics in red.

Figure 3 :

 3 Figure 3: Four-wave mixing dynamics as a function of real time t and delay τ . (a) Experiment and (b) simulation.

Figure 3 (

 3 b) shows the corresponding simulated FWM amplitude |S FWM | within a few-level system where we take the exciton transfer into the optically uncoupled valley as well as the exciton-exciton interaction in terms of the local field model into account. It describes the optically induced dynamics of the exciton's coherence p and its occupations n and n via[START_REF] Wegener | Line shape of time-resolved four-wave mixing[END_REF][START_REF] Hahn | Influence of local fields on the dynamics of four-wave mixing signals from 2D semiconductor systems[END_REF][START_REF] Rodek | Local field effects in ultrafast lightmatter interaction measured by pump-probe spectroscopy of monolayer MoSe 2[END_REF]

Figure 4 :

 4 Figure 4: Six-wave mixing dynamics as a function of real time t and delay τ . (a) Experiment, (b) simulation, with the curved dashed line depicting Eq. (9).

Figure 5 :

 5 Figure5: Flow chart for the three main contributions to the SWM signal listing intermediate phase-filtered coherences (p, blue) and occupations (n, red). Green waved arrows show pulse induced, violet dotted ones local field induced, black ones free dynamics without, and blue ones with dephasing. The phase differences are defined as ∆ nm = φ nφ m . The flow chart holds for any order of the optical field. However, in the linear response regime considered here it is 0 n + 1 = 0 n + 2 = 0 n - 3 = 0.

Figure 6 :

 6 Figure 6: Schematic picture of the pulse sequence and signal development from the two main contributions.

Figure 7 :

 7 Figure 7: Amplitude and phase dynamics of the destructive echo. (a) SWM amplitude dynamics for a delay of τ = 0.35 ps. Experiment as solid blue and simulation as dashed red line. (b) Respective dynamics of the phase of the SWM signal.

Figure 8 :

 8 Figure 8: Quasi-Bloch vector of the FWM contribution. (a) Lissajous curve of the initial distribution of quasi-Bloch vectors. (b) Dynamics of the quasi-Bloch vectors. Initial dis-tribution as bright open circles and for a time t > 0 after the second pulse as filled dark circles. The final integrated FWM coherence FWM p is marked as blue cross and its propagation is given by V t.

integral over ñ

  always vanishes and we are left with the final FWM coherence FWM p(t) = 2π 0 d∆ 21 p(t; ∆ 21 ) .

Figure 9 :

 9 Figure9: Quasi-Bloch vector distributions after the laser pulse. The ∼ t 2 contribution (violet) and ∼ τ t (blue) form Lissajous curves with frequency ratios 2:2:1 and 1:1:1, respectively.

2 )Figure 10 :

 210 Figure10: Projections of the quasi-Bloch vector representation of the destructive echo formation. Columns show projections of the quasi-Bloch vector defined in Eq. (14) as listed on the top, the rows show different real times as given on the left. The contribution ∼ tτ is shown as blue, ∼ t 2 as violet dots. The size of the violet dots represents the product ñ Re(p), positive values are larger and negative ones smaller. The situation for t = 0 is shows as bright circles for all other times. The black arrows indicate the movement of the quasi-Bloch vectors, solid (dashed) for positive (negative) ñ Re(p). The final SWM coherences are marked as crosses in the two right columns.

2 Figure S2 :- 1 Figure S3 :

 2S21S3 Figure S2: SWM dynamics for different pulse areas θ being halved (a) or doubled (b) compared to Fig. 4(b) in the main text.

  5 0 0.5 1 1.5 2 -0.5 0 0.5 1 1.5 2 delay τ (ps) time t (ps) a W = 50 ps -1 time t (ps) b W = 0 ps -1

Figure S4 :

 S4 Figure S4: Impact of a non-vanishing excitation induced dephasing with W = 50 ps -1 (a) and (b) is the example from Fig. 4(b) in the main text with W = 0.

  a) for W = 50 ps -1 . As reference in Fig. S4(b) we show the case without EID (W = 0) from the main text [Fig. 4(b)]. Obviously, the EID prohibits the development of a destructive photon echo.

Figure S5 :

 S5 Figure S5: Impact of a non-vanishing excitation induced dephasing with W 10 ps -1 on the SWM amplitude dynamics for τ = 0.35 ps. Experiment as solid blue line and simulation as dashed red line. The measurement is the same as in Fig. 7(a).

Figure S6 :

 S6 Figure S6: Impact of a vanishing intervalley scattering λ = 0 on the SWM amplitude dynamics. (a) SWM amplitude dynamics. (b) Difference of the SWM amplitude with respect to the case with λ = 4 ps -1 from Fig. 4(b) in the main text.

Figure S7 :

 S7 Figure S7: Impact of a non-vanishing pulse delay τ 12 on the SWM amplitude dynamics. (a) Normalized SWM amplitude dynamics for τ 12 = 0 from Fig. 4(b) in the main text. (b) SWM amplitude dynamics for τ 12 = 0.35 ps normalized to the maximum value of (a). The insets schematically show the timing of the pulses with φ 1 and φ 2 .

Figure S8 :

 S8 Figure S8: Impact of a non-vanishing inhomogeneous broadening on the SWM amplitude dynamics. (a) SWM amplitude dynamics with destructive echo from Fig. 4(b) in the main text. (b) SWM amplitude dynamics for V = 0 but with an inhomogeneous broadening of σ = 4 ps -1 resulting in a constructive photon echo marked by the dotted black lines. (c) Combination of constructive and destructive photon echo. The signal depression is only visible in the vicinity of the constructive echo, while the entire signal is suppressed for all other times.

Table 1 :

 1 Comparison of criteria defining the photon echo in FWM and the destructive photon echo in SWM.

		echo	destructive echo
		(FWM)	(SWM)
	timing	delay + internal dynamics	delay + internal dynamics
	interference	constructive (1 path)	destructive (2 paths)
	source	inhomogeneous broadening	local field induced phase inhomogeneity
	quantum interference of all depicted quasi-Bloch vectors. For
	even larger times t < τ (bottom row) the accelerated ∼ t 2 contribution is stronger than the one with constant velocity ∼ τ t and the SWM is non-vanishing again.

Acknowledgements

T.H. thanks the German Academic Exchange Service (DAAD) for financial support (No. 57504619). D.W. thanks the Polish National Agency for Academic Exchange (NAWA) for financial support within the ULAM program (No. PPN/ULM/2019/1/00064). T.H., D.W., T.K., and P.M. acknowledge support from NAWA under an APM grant. M.P, D.V., and M.B thank the EU Graphene Flagship project. M.P, D.V, K.N, and J.K acknowledge support from CNRS via IRP "2D materials". M.P. and K.N acknowledge support from TEAM programme of the Foundation for Polish Science, co-financed by the EU within the ERDFund. M.P. is supported by the ESF under the project CZ.02.2.69/0.0/0.0/20_079/0017436. K.W. and T.T. acknowledge support from the Elemental Strategy Initiative conducted by the MEXT, Japan (Grant Number JPMXP0112101001) and JSPS KAKENHI (Grant Numbers 19H05790 and JP20H00354).

Destructive photon echo formation in six-wave mixing signals of a MoSe 2 monolayer (Supplementary Information)

Thilo Hahn, 1,2 Diana Vaclavkova, 3 Miroslav Bartos, 3,4 Karol Nogajewski, 5 Marek Potemski, 3,5 Kenji Watanabe, 6 Takashi Taniguchi, 7 Paweł Machnikowski, 2 Tilmann Kuhn, 1 Jacek Kasprzak, 8 Daniel Wigger mixing the three-pulse FWM term ∆32+φ1 p + 3 = φ3-φ2+φ1 p + 3 and

The second one (iv) is created by the two-pulse FWM term 2φ3-φ2 p + 3 mixed with |∆21| n + 3 leading to

Finally, contribution (v) is created by local field mixing

For completeness we also give the three SWM contributions with O(V 2 ) here. As illustrated in Fig. 5 in the main text they are retrieved by double local field mixing between φ1 p + 3

between φ3 p + 3 , ∆32 n + 3 , and

where the factor 2 stems from the sum of the two occupations, and finally the single local field mixing between φ3 p + 3 and ∆32-∆21 n

Collecting all the expressions finally yields

for the SWM signal in the χ (5) -regime.

From Eq. ( 7)/(S8) we see that in the lowest order of the light-matter coupling, the signal amplitude scales with the fifth power of the pulse area. Therefore, operating in the low excitation regime the SWM dynamics do not depend on the choice of θ. To confirm this in Fig. S2 we plot the same SWM dynamics as in Fig. 4(b) in the main text but with halved (a) and doubled (b) pulse area. There are no obvious differences between the two simulations. This confirms that the chosen pulse area