Joule overheating poisons the fractional ac Josephson effect in topological Josephson junctions

Topological Josephson junctions designed on the surface of a 3D-topological insulator harbor Majorana bound states among a continuum of conventional Andreev bound states. The distinct feature of these Majorana bound states lies in the 4π-periodicity of their energy-phase relation that yields a fractional ac Josephson effect and a suppression of odd Shapiro steps under radio-frequency irradiation. Yet, recent experiments showed that a few, or only the first, odd Shapiro steps are missing, casting doubts on the interpretation. Here we show that Josephson junctions tailored on the large bandgap 3D-topological insulator Bi2Se3 exhibit a fractional ac Josephson effect acting on the first Shapiro step only. With a modified resistively shunted junction model, we demonstrate that the resilience of higher order odd Shapiro steps can be accounted for by thermal poisoning driven by Joule overheating. Furthermore, we uncover a residual supercurrent at the nodes between Shapiro lobes, which provides a direct and novel signature of the current carried by the Majorana bound states. Our findings showcase the crucial role of thermal effects in topological Josephson junctions and lend support to the Majorana origin of the partial suppression of odd Shapiro steps. The search for experimental evidence of Majorana modes is an area of intense research in condensed matter and quantum physics and uncovering clear evidence is complicated. The authors investigate the impact of Joule heating which can influence the analysis of experimental features related to Majorana bound states in topological Josephson junctions.

Another key approach to substantiate the very existence of MBS relies on the fractional AC Josephson effect [13][14][15] that develops in topological Josephson junctions 3,15 . Theory predicts that MBSs shall emerge in such junctions as a peculiar, spinless Andreev bound state (ABS). Contrary to the conventional ABSs whose energy level varies 2π-periodically with the phase difference ϕ between the junction electrodes, the MBS is 4π-periodic and crosses zero-energy for a phase π (see Fig. 1a), yielding a fractional AC Josephson effect at frequency f J /2 = eV/h (e is the electron charge, V the voltage drop across the junction and h the Planck constant), that is, half the Josephson frequency f J 3, 15 .
Yet, revealing such a 4π-periodic contribution has proven challenging in dc transport experiments due to the presence of often prevailing, conventional ABSs [16][17][18] . Moreover, poisoning processes-stochastic parity-changes of the quasiparticle occupation number-may obscure the MBS contribution by limiting its lifetime 19,20 . Measurement schemes probing at time scales shorter than this lifetime are thus essential. The Shapiro effect comes forth with the combined advantages of a radio-frequency (f rf ) excitation of the phase that can be faster than the poisoning dynamics 15,21,22 , and the ease of dc current-voltage (I-Vs) characteristics measurements.
The immediate consequence of the fractional AC Josephson effect is an unusual sequence of Shapiro voltage steps ΔV ¼ hf rf e in the I-V characteristics, twice of that of conventional Shapiro steps hf rf 2e 13,21-26 , providing direct evidence for the MBS 4π-periodicity. First experiments performed on InSb nanowires 27 , on strained HgTe 3D TI 28 , and on Bi 1−x Sb x 29 junctions however showed surprises in the sequences of Shapiro steps. In all these cases, only the n = ±1 steps were absent in a given range of radiofrequency (rf) power and frequency (n is the integer index of the Shapiro steps), an absence which was described as an (incomplete) signature of the fractional a.c. Josephson effect. More recent measurements on 2D HgTe quantum wells showed the absence of odd steps up to n = 9 30 and Josephson radiation at half the Josephson frequency 31 , though without the demonstration of time-reversal symmetry breaking that is required to induce MBSs in quantum spin-Hall edge channels 15,32 . While the latter observations advocate more clearly for the existence of 4π-periodic Andreev modes, the fact that the fractional AC Josephson effect acts only on some odd Shapiro steps depending on the system remains unclear. Whether it provides a signature of the Majorana mode is a central question for identifying topological superconductivity in a variety of implementations.
In this work we report on the observation and understanding of the partial fractional AC Josephson effect in Josephson junctions designed on exfoliated flakes of the 3D-topological insulator Bi 2 Se 3 . Our devices exhibit an anomalous sequence of Shapiro steps with (only) the first step absent at low rf power and frequency. The 4π-periodic contribution to the supercurrent is directly identified as a residual supercurrent at the first node of the critical current when the rf power is increased. To shed light on our findings, we develop a two-channel resistively shunted a Energy-phase spectrum of the Andreev bound states for a topological Josephson junction at the surface of a 3D-topological insulator. 4π-periodic, spinless, Majorana bound states coexist with conventional 2π-periodic Andreev bound states (ABS). The maximum energy E ABS is lower than the quasiparticle continuum at Δ in case of imperfect interface transparency or the presence of a magnetic barrier 3,15,43 . b Optical image of device. Scale bar is 9 μm. c Schematic of the Josephson junction geometry showing the superconducting electrodes (S) in orange that contact the top of the Bi 2 Se 3 flake (in blue). d Current-voltage characteristic of the junction indicated by the arrow in c. Measurements were carried out at 0.05 K junction (RSJ) model that includes the quasiparticle overheating induced by Joule effect [33][34][35] , and a thermally activated poisoning of the MBS. We show that Joule overheating suppresses the parity lifetime of the MBS and thus terminates the 4π-periodic contribution to any higher index Shapiro steps, accounting for the observed suppression of the first Shapiro step only.

Results
Partial even-odd effect in Bi 2 Se 3 Josephson junctions. Our samples are based on flakes of Bi 2 Se 3 crystals exfoliated with the scotch tape technique. Figure 1b shows a 30 nm thick flake of Bi 2 Se 3 contacted with multiple electrodes of vanadium enabling both magneto-transport and Josephson junction measurements on the same Bi 2 Se 3 crystal (see Appendix for fabrication and measurement details). Analysis of Shubnikov-de-Haas oscillations and Hall effect enable identification of three electronic populations contributing to the sample conductance: Bulk states with a charge carrier density of 4.5 × 10 19 cm −3 , and the top and bottom surface states with densities of 1 × 10 12 and 4 × 10 12 cm −2 respectively (see Supplementary Note 1). All three channels may thus carry supercurrent by proximity effect 36 .
We focus here our discussion on the Josephson junction of length L = 125 nm and width W = 2.25 μm (see geometry in Fig. 1c) indicated by the white arrow in Fig. 1c. Below the superconducting transition temperature of the electrodes (T c = 5 K), the proximity effect develops in the TI, leading to a dissipationless supercurrent in the I-Vs as shown in Fig. 1d (see Supplementary Note 2 for I-Vs at higher current bias showing an excess current due to Andreev reflections). The transition to the resistive state of the junction (R = 7.5 Ω) is hysteretic at 0.05 K with switching and retrapping currents of I sw = 7.3 μA and I r = 5.0 μA, respectively. Such a hysteresis is a common feature of most Josephson junctions made with metallic weak links and results from a quasiparticle overheating in the normal section of the junction 34 . As we will show below, the ensuing quasiparticle overheating is key for understanding the suppression of the n = ±1 Shapiro steps only.
The dc response of the Josephson junction to an rf irradiation is shown in the Shapiro maps of Fig. 2a, b that display the colorcoded differential resistance dV/dI versus rf power P rf and dc current I. For an rf frequency of 3.5 GHz, Fig. 2a, well-defined Shapiro steps develop in the I − V curves, two of which are shown in Fig. 2c, with voltage steps that match the standard value V n ¼ ± n hf rf 2e ¼ n 7:2 μV expected for a 2π-periodic current-phase relation. In the Shapiro map, the black areas indicate dV/dI = 0 and hence the position and amplitude of the Shapiro steps in the P rf − I plane. Two features standard for a usual Shapiro map are visible. First, on increasing P rf , the critical current continuously decreases till nearly full suppression at P rf = 9 dBm and then oscillates at higher P rf . Second, the sequence of appearance of the Shapiro steps with P rf is sorted by the Shapiro step index n, and, importantly, starts with the step n = 1.
The central experimental result of this work is displayed in Fig. 2b, where we show the rf response of the same junction at a lower frequency of 1 GHz. This Shapiro map exhibits distinct features that markedly differentiate it from the higher frequency map. On increasing P rf , the first Shapiro step n = 1 sets in at high P rf after the steps of higher index for both switching and retrapping currents. This anomaly, sometimes termed even-odd effect 21,22 , results in the conspicuous absence of the first Shapiro step in I-Vs picked out at low P rf , while steps of higher indexes already appears, see Fig. 2d (similar results obtained on another sample are shown in Supplementary Note 3). Our findings match those recently obtained on InSb nanowires 27 , strained HgTe 3D TI 28 and Bi 1−x Sb x alloy 29 , which were interpreted as a signature of a 4π-periodic MBS contribution. A second and new feature emerges at the first minimum of the critical current when the rf power is increased, i.e., at P rf ≃ 10.7 dBm. Contrary to a conventional Shapiro map where a complete supercurrent suppression is expected at what can be called a resistive node, a small supercurrent I c~1 90 nA remains. This appears clearly in the individual I-Vs of Fig. 2d, see for instance the blue curve there. Figure 2e displays a similar Shapiro map obtained at a slightly different f rf in the low-frequency regime, but zoomed on the resistive node where the critical current is expected to vanish but does not. We shall see in the following that this residual supercurrent provides a direct signature of the presence of a 4π-periodic mode in the junction.
Determination of the coherent transport regimes. Capturing the ABS spectrum of a Josephson junction on 3D TIs remains difficult as several conduction channels, including bulk and surfaces, intervene. Should all channels carry supercurrent, the nature of charge carriers in them may lead to virtually different regimes of coherent transport that we assess in the following. For bulk carriers, a rough estimate of the mean-free path (see Supplementary Note 1) gives a Thouless energy E th = ħD/L 2 ≃ 417 μeV, smaller than the superconducting gap of the vanadium electrodes Δ = 800 μeV. This channel thus belong to the class of long diffusive Josephson junctions 37 . In contrast, for the topological surface states, the spin texture of the Dirac electrons stemming from the spin-momentum locking leads to a very strong scattering anisotropy which promotes forward scattering. As a result, the transport time τ tr is expected to be significantly enhanced with respect to the elastic scattering time τ e , with ratio τ tr /τ e up to 60 depending on the disorder source 38 . Recent experiments on Bi 2 Se 3 flakes combining field effect mobility and quantum oscillations assessed a ratio τ tr /τ e ≳ 8 39 . Taking the latter value as a conservative estimate and the surface state elastic mean-free path l e ≃ 28 nm of our sample (see Supplementary Note 1) leads to a transport length l tr ≳ 225 nm. These considerations suggest that surface transport is ballistic with, importantly, a non-zero probability for straight electronic trajectories impinging both electrodes. This is also consistent with signatures of ballistic transport over 300 nm evidenced in Bi 2 Se 3 nanowires 40,41 .
Consequently, we consider the topological surface state channel as ballistic. As such, the relevant energy scale for the ABSs is ħv F /L = 2.8 meV with v F = 5.4 × 10 5 m/s the Fermi velocity 42 . It is greater than Δ, which shall lead to ABSs in the short ballistic limit. Theory then predicts that a 4π-periodic MBS exists, even in the case where the Fermi level is far from the Dirac point of the surface states, which is here the experimentally relevant regime 43 . This MBS corresponds to ballistic trajectories impinging perpendicularly the superconducting electrodes, all other incidence angles yielding conventional 2π-periodic ABSs 3,43 .
Eventually, observability of the 4π-periodicity theoretically implies a strong constraint on the Andreev spectrum: The MBSs must be decoupled from the quasiparticle continuum at ϕ = 0 and 2π to avoid direct transfer of quasiparticles into or from the continuum. Such a transition would indeed occur from the excited to the ground state every 2π, restoring an effective 2πperiodicity for the MBS. This detrimental effect can be remedied by adding a magnetic layer or magnetic field that break timereversal symmetry, and thus open a gap between the MBS and the quasiparticle continuum 15,21,23,43,44 , as sketched in Fig. 1a. In our samples, the vanadium that we use as superconducting electrodes is known to form magnetic dopants in Bi 2 Se 3 and eventually a ferromagnetic phase at large concentration 45 . Given that a smooth ion milling of the Bi 2 Se 3 surface is processed before vanadium deposition, favoring vanadium diffusion into the Bi 2 Se 3 crystal, there is presumably a magnetic layer or local magnetic moments at the superconducting interface as well as magnetic moments on the oxidized vanadium side surfaces of the electrodes. This singular configuration is likely to break timereversal symmetry on the scale of the junction, thus leading to the decoupling of the MBS from the continuum and to the ensuing observability of 4π-periodicity in our Shapiro maps.
Two-channel RSJ model. To understand our experimental findings, we consider a RSJ model comprising a pure Josephson junction in parallel with a shunt resistor R. In the usual scheme of a single Josephson channel with a critical current I c , the key parameter for the phase dynamics is the phase relaxation time which sets the typical time scale for the phase to adapt to a drive current change 33 . With a rf drive, the RSJ model thus acts on the phase as a low pass filter of cutoff frequency 1/τ J (see Supplementary Note 4). The regime of visibility of Shapiro steps is thus defined by f rf τ J < 1.
To phenomenologically capture the complex dynamics of a topological Josephson junction where a MBS lies within (a majority of) ABSs, we include two different Josephson junctions J1 and J2 in the RSJ model (see Fig. 3a) 24,26 . The first junction J1 stands for the conventional ABSs with a 2π-periodic current-phase relation, and the second one J2 represents the 4π-periodic MBSs.
Having both 2π and 4π-periodic contributions in the total supercurrent I s ðϕÞ ¼ I 2π c sinðϕÞ þ I 4π c sinðϕ=2Þ drastically changes the Shapiro steps sequence. The dynamics is now ruled by two different phase relaxation times τ 2π J ¼ h and I 4π c . Consequently, the 4π-periodic contribution will impact the junction dynamics only for drive frequencies f rf < 1=τ 4π J . This can be straightforwardly seen in the Shapiro maps that we obtained by numerically solving the RSJ equation together with the Josephson . In the low-frequency limit f rf τ 2π J < f rf τ 4π J < 1, all even Shapiro steps develop at rf power P rf lower than their neighboring odd steps, leading to the following appearance sequence |n| = {2, 1, 4, 3, 6, 5…} on increasing P rf (see Fig. 3b). Lowering f rf would enhance this even-odd effect and result ultimately in a quasi-suppression of odd Shapiro steps. Conversely, when the drive frequency is faster than 1=τ 4π J but still lower than 1=τ 2π J , that is, f rf τ 2π J < 1 < f rf τ 4π J , the 4π-periodic contribution is suppressed, restoring the regular sequence of Shapiro steps appearance |n| = {1, 2, 3, 4…} on increasing P rf , as shown in Fig. 3c. In a 2D space (f rf , P rf ), the ranges of existence of the even and odd Shapiro steps are sketched in Fig. 3d. Importantly, the even-odd effect is robust even if I 4π c sounds negligible compared to I 2π c , since the even-odd effect will always be present at low enough frequency as soon as f rf < 1=τ 4π J . This low-frequency observability of the even-odd effect furthermore excludes an explanation for the existence 4π-periodic contribution based on Landau-Zener transitions at a soft gap, which should be otherwise enhanced at high frequencies.
Let us now consider in the same computed maps the P rfdependence of the Shapiro steps amplitude. In contrast with the standard oscillatory behavior with a complete suppression at the resistive nodes, the even steps, including the supercurrent branch (n = 0), exhibit a non-vanishing amplitude at every two resistive node on increasing P rf , see Fig. 3b, c. This unusual feature has been predicted in a recent paper by Domínguez et al. 26 . It is a direct consequence of the presence of the 4π-periodic channel. It also conspicuously matches the residual supercurrent at the first node of the supercurrent branch in the experimental data of Fig. 2b, e.
We demonstrate below that the amplitude of this residual supercurrent can be quantitatively related to the 4π-periodic critical current I 4π c . Figure 4a displays the computed switching current I 0 of the Shapiro step n = 0, extracted from Fig. 3b as a function of P rf . This plot both highlights the oscillatory behavior of I 0 with P rf and enables us to identify the residual supercurrent of the first node that we note I k¼1 0 , with k the node index. To demonstrate the correlation between I k¼1 0 and I 4π c , we performed a numerical study of the dependence of I k¼1 0 on the relevant parameters I 4π c and f rf by systematically computing Shapiro maps for different sets of parameters. Figure 4b displays the calculated I k¼1 0 =I 4π c versus f rf τ 4π J for different I 4π c =I c ratios (we define here the total critical current I c ¼ I 2π c þ I 4π c ). All I k¼1 0 =I 4π c values collapse on a single curve which tends to saturate to~0.7 in the limit f rf τ 4π J ( 1. Conversely, I k¼1 0 vanishes when f rf τ 4π J $ 16, indicating that this residual supercurrent is visible at higher frequencies than the even-odd effect on the step appearance order. The remarkable consequence of the collapse is that I 4π c is uniquely defined for a fixed set of parameters I k¼1 0 , R and f rf . Using a polynomial fit of the data points in Fig. 4b ( Joule-induced poisoning of the MBS. The two-channel RSJ model, however, does not account for the experimental absence of the first Shapiro step only. Inclusion of a junction capacitance geometric (or instrinsic 47 ) through a two-channel RSCJ model mitigates the suppression of the n ≥ 3 odd steps 48 , but is however not relevant for our strongly overdamped Bi 2 Se 3 -based Josephson junctions (With a geometrical capacitance of the order of 10 aF estimated with a planar capacitor approximation for the superconducting electrodes, we obtain a damping parameter σ ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi h=2eI c R 2 C p $ 260. We also estimate an intrinsic capacitance C *~0 .1 pF predicted in ref. 47 which leads to σ = 2.5. This overdamped regime is in contradiction with the hysteresis seen in the I-Vs, consequently bearing out the electron overheating scenario.). Furthermore, in the case of strongly underdamped, hysteresis should also be sizable in the Shapiro steps 49 , which is not observed in our experiment.
We propose instead that the absence of the first step can be explained by including thermal effects resulting from Joule heating, and their impact on the 4π-periodic mode. Some of us recently showed that electron overheating must be taken into account with the RSJ model to capture Shapiro maps in conventional metallic Josephson junctions 35 . Assuming that quasiparticles in the junction form a thermal distribution with an effective temperature T qp different from the phonon bath temperature T ph , we included the temperature dependence of the critical current, I c (T qp ), and solved the RSJ equation selfconsistently together with the heat balance equation P ¼ IðtÞVðtÞ h i¼ ΣΩ T 5 qp À T 5 ph (Σ is the electron-phonon coupling constant and Ω the volume of the normal part, see Supplementary Note 6 for values) to extract T qp for each dc current. Solving such a thermal RSJ model (tRSJ), we obtain a significant raise of T qp , when a dc voltage drop sets in on the firstly developed Shapiro step 35 .
We then conjecture that the ensuing excess of non-equilibrium quasiparticles that we express in terms of an effective quasiparticle temperature poisons the 4π-periodic mode. Note that any non-thermal distribution would have the same consequences. Acting on a single mode, poisoning causes a stochastic paritychange 15,21,22,31,44 that switches in time the quasiparticle occupation from the excited to the ground state, as illustrated in Fig. 1a by the black dotted arrow. Within the two-channel RSJ model, we model this parity-change of the 4π-periodic contribution to the supercurrent as I 4π s ðtÞ ¼ ðÀ1Þ n sw ðtÞ I 4π c sin where n sw (t) is a random occupation number of characteristic time scale determined by a switching parity lifetime τ sw (see Supplementary Note 7). The exact microscopic processes that yield such a diabatic event can involve several transitions, including pair breaking, quasiparticle recombination with various rates, and possible coupling to the bulk states. Phenomenologically, but without loosing the generality of the foregoing, we follow the approach of Fu and Kane 15 and consider a thermally activated switching parity lifetime τ sw for the 4π-periodic mode: where Δ − E MBS is the minimal energy gap separating the 4πperiodic mode to the continuum (see Fig. 1b), τ 0 the characteristic time scale for the activation process, and k B the Boltzmann constant. When the quasiparticle effective temperature T qp is high, the switching time is small τ sw ≲ τ J and poisoning suppresses the 4π-periodic Josephson effect, and hence the even-odd effect, by switching frequently the current-phase relation (see Supplementary Note 7). For the opposite limit τ sw ) τ J , poisoning is irrelevant. We thus claim that understanding the partial even-odd effect in topological Josephson junctions relies on the interplay between activated poisoning and thermal effects in the two-channel tRSJ model. We show in Fig. 5a, b two Shapiro maps computed for the same rf frequency and junction parameters as those in Fig. 3b. Figure 5a is the result of the tRSJ model that includes a Tdependence of I 2π c ðTÞ given by the measured I sw (T) 35 and realistic parameters for the heat balance equation (see Supplementary Note 6 for an estimate of ΣΩ in our junctions). We assume here that most of the T-dependence of I sw (T) comes from the 2πperiodic modes. Compared to Fig. 3b, the electron overheating leads to a broadening of the resistive transitions at high P rf between the periodic oscillations of the Shapiro steps. Nevertheless, the full even-odd effect acting on all odd steps remains as in the absence of thermal effect, see Fig. 3b. Figure 5b is the main result of our theoretical analysis. It displays a Shapiro map computed with the same tRSJ than Fig. 5a, but including thermally activated poisoning of the 4π-periodic channel defined by Eqs. (1) and (2) which ensues from Joule overheating. The 4π-periodic channel now impacts only the steps n = ±2 by enhancing their amplitude, therefore inverting the appearance order on increasing P rf between step n = 1 and 2. The sequence of appearance of all higher order steps turns out to be regularized due to the suppression of the 4π-periodic contribution by poisoning. This finding, that is, the even-odd effect limited to the first Shapiro step only, is in full agreement with our experiment shown in Fig. 2b and with works on other systems [27][28][29] .
The thermally activated poisoning can be captured by inspecting the computed T qp and τ sw for two different P rf . Figure 6a displays the I-V curves corresponding to the black arrows in Fig. 5b. In Fig. 6b, we show the corresponding T qp versus I, which raises linearly once a dissipative voltage sets in. Accordingly, τ sw is exponentially suppressed, and becomes inferior to τ J on the Shapiro steps n ≥ 2 (Fig. 6c). This explains why the 4π-periodic component acts only on the appearance order of the n = 1 and 2 steps, leaving all other steps unaffected.
Furthermore, inspecting Fig. 5b we see that the residual supercurrent of only the first resistive node I k¼1 0 remains in presence of poisoning, as observed in the data of Fig. 2b. This confirms that the residual supercurrent I k¼1 0 is a robust feature that provides a new indicator for the presence of 4π-periodic modes, and enables a direct and quantitative determination of the corresponding 4π-periodic critical current as discussed above.
Within this approach of thermal poisoning by Joule overheating, it is interesting to compare the power dissipated between the available experiments on different TI systems. For instance, estimates of the dissipated power on the n = 2 Shapiro step give two orders of magnitude difference between different TI materials: P ' 2 hf rf 2e I $ 25 pW for our data and 17 pW for strained HgTe 28 . For InAs nanowires 27 and Bi 1−x Sb x 29 , P~2 pW. For HgTe quantum wells 30 , the dissipated power is the smallest P~0.2 pW. Comparing the strained HgTe to the HgTe quantum wells that share the same electron-phonon coupling constant, we estimate a power per unit of volume of~80 fW μm −3 and~2 fW μm −3 , respectively (~300 fW μm −3 in our Bi 2 Se 3 samples). Such a small dissipated power density in the experiment on HgTe quantum wells should result in a minimized amount of non-equilibrium quasiparticles and limited poisoning, therefore explaining the observed full suppression of not only the first but of several odd Shapiro steps in this system.

Discussion
Our theoretical approach, combining two Josephson channels in parallel, electronic overheating and quasiparticle poisoning, can be extended to more elaborate situations. For instance, recent theory works predict for the case of the quantum spin Hall regime an 8π-periodicicity due to either interactions 50 or quantum magnetic impurities 51,52 . Although this has not been reported so far in experiments, it would be interesting to study how multiple periodicities mix in the Shapiro response. Given the understanding of the two-channel RSJ model, we expect the 8π-periodicity to enhance every steps of index ±4n and significantly modify the beating pattern at high P rf . Other effects such as the voltage dependence of the phase relaxation time 21 could be included in our model and should enhance the effect of thermal poisoning on the partial suppression of the even-odd effect.
To conclude, our work elucidates the origin of the puzzling suppression of only the first Shapiro step in topological Josephson junctions. In our Bi 2 Se 3 Josephson junctions, this suppression is accompanied by a residual supercurrent that provides a new indicator of the 4π-periodic contribution to the supercurrent. Together, these observations can be captured by a two-channel thermal RSJ model in which Joule overheating activates poisoning of the 4π-periodic mode. The even-odd effect restricted to the first Shapiro step and the residual supercurrent do provide a clear signature of a 4π-periodic mode in the Andreev spectrum, conspicuously pointing to MBSs. Our phenomenological model illustrates a direct consequence of thermal poisoning on MBSs, signaling that dissipation must be scrutinized with attention in dc-biased measurements. Addressing a microscopic description of the enhanced poisoning in such non-equilibrium measurement schemes is a challenging task for theory that should lead to significant progress towards new devices for Majorana physics and possible MBS qubits.

Methods
Bi 2 Se 3 crystals were synthesized by melting growth method with high purity (5N) Bi and Se in an evacuated quartz tube. Crystals were analyzed by X-rays diffraction, and angle-resolved photoemission spectroscopy. Flakes of Bi 2 Se 3 were exfoliated from the bulk crystal on silicon wafer and systematically inspected by atomic force microscopy to ensure crystal quality. V/Au superconducting electrodes were patterned by e-beam lithography and deposited by e-gun evaporation after a soft ion beam etching.
Measurements were performed in a dilution refrigerator equipped with highly filtered dc lines that comprise room temperature feed-through Pi-filters, lossy custom-made coaxial cables and capacitors to ground on the sample holder. Radio-  Fig. 5 Poisoning of the even-odd effect. a Shapiro map computed with the two-channel thermal resistively shunted junction (RSJ) model displaying the differential resistance dV/dI versus dc current I normalized to the critical current I c and radio-frequency power P rf . For the sake of clarity, we took the same parameters as in Fig. 3b, that is, f rf τ 4π J ¼ 0:4 (see Supplementary Note 8), with a phonon bath temperature T ph = 0.1 K. The model includes the experimental T-dependence for I 2π c ðTÞ and an estimate of the electron-phonon coupling in Bi 2 Se 3 (see Supplementary Note 6). On each point a recursive algorithm solves successively the RSJ equation, the ensuing time dependent voltage and dissipated power, the raise of the quasiparticle temperature T qp , and then re-solves the RSJ equation with the new I 2π c T qp till convergence. b Shapiro map computed with the two-channel thermal RSJ model as in a but including the thermal poisoning. The implementation of a thermally activated poisoning suppresses the even-odd effect for the Shapiro step of indexes n ≥ 3 frequency are fed through a dedicated coaxial cable ending as an antenna that were adjusted in the vicinity of the devices. Shapiro map measurements were performed with standard lock-in amplifier technique.

Data availability
The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.  Fig. 5b at two different radio-frequency powers P rf indicated by the black arrows in Fig. 5b, showing the voltage V normalized to hf rf /2e (a), Quasiparticle temperature T qp (b), and the switching time τ sw (c) versus current I normalized to the critical current I c . At P rf = −112 dBm, the first Shapiro step is absent in the I − V. The increase of T qp when a voltage sets in on Shapiro steps leads to an exponential suppression of τ sw . When τ sw < τ 4π J (the phase relaxation time τ 4π J is indicated by the black dotted line in c), poisoning suppresses the contribution of the 4π-periodic channel to the Shapiro steps n > 2, and all odd Shapiro steps are present in the I − V, see blue curve in (a)