# Topological supercurrents interaction and fluctuations in multiterminal Josephson effect

###### Abstract

We study Josephson effect in the multiterminal junction of topological superconductors. We use symmetry-constrained scattering matrix approach to derive band dispersions of emergent sub-gap Andreev bound states in a multidimensional parameter space of superconducting phase differences. We find distinct topologically protected band crossings that serve as monopoles of finite Berry curvature. Particularly, in a four-terminal junction the admixture of and periodic levels leads to an appearance of a finite energy Majorana-Weyl nodes. This topological regime in the junction can be characterized by a quantized nonlocal conductance that measures the Chern number of the corresponding bands. In addition, we calculate current-phase relations, variance, and cross-correlations of topological supercurrents in multiterminal contacts and discuss universality of these transport characteristics. At the technical level these results are obtained by integrating over the group of a circular ensemble that describes the scattering matrix of the junction. We briefly discuss our results in the context of observed fluctuations of the gate dependence of critical current in topological planar Josephson junctions and comment on the possibility of parity measurements from the switching current distributions in multiterminal Majorana junctions.

## I Introduction

The universality of conductance fluctuations (UCF) is the hallmark of mesoscopic physics BLA-JETPLett85 (); BLA-DEK-JETPLett86 (); Lee-Stone-PRL85 (); Lee-Stone-Fukuyama-PRB87 (). This phenomenon emerges from the quantum coherence of electron trajectories and is sensitive to changes in external magnetic field or gate voltage. At temperatures below the Thouless energy, , which is related to the inverse dwell time for an electron to diffuse across the sample , root-mean-square value of conductance fluctuations saturates to the universal value of order conductance quantum as long as characteristic sample size is smaller than dephasing length . Interaction effects in normal metals barely change the magnitude and universality of conductance fluctuations although they are crucially important in determining temperature dependence of dephasing effects and in particular Aleiner-Blanter-PRB02 (). Robustness of UCF can be rooted to random matrix theory description of Wigner-Dyson statistics of electron energy levels in disordered conductors Dyson (). Indeed, in the Landauer picture of transport across a mesoscopic sample, conductance is given by times the number of single particle levels within the energy strip of the width of Thouless energy. While the average number of such levels depends on the dimensionality, random matrix theory predicts that their fluctuation is universally of the order of one Altshuler-Shklovskii-JETP86 (); Mehta ().

When superconductivity is induced at the boundary of the mesoscopic sample via the proximity effect universality of fluctuations remains intact Exp-UCF-SN-1 (); Exp-UCF-SN-2 (). Indeed, the magnitude of sample-to-sample conductance fluctuations changes only by a numerical factor of the order of unity whose value depends on the underlying symmetry Brouwer-PhD (); Beenakker-RMP (). Interestingly, universality of fluctuations extends beyond conductance as it also manifests in the Josephson current of SNS bridge. Indeed, extending the original ideas of Altshuler and Spivak BLA-BZS-JETP87 (), who argued that random shifts of sub-gap energy levels with superconducting phase difference would alter the current, Beenakker showed Beenakker-PRL91 () that in short junctions, , where is the superconducting coherence length, root-mean-square value of critical current fluctuations saturates to a universal bound determined only by superconducting energy gap in the leads. Further a complete characterization of the supercurrent variance as a function of phase across the point contact Josephson junction was computed by Chalker and Macêdo Chalker (). In long junctions, , supercurrent fluctuations cease to be universal and scale with . However a remarkable property of these fluctuations is that there is a regime where the entire critical current through the junction can be determined by the mesoscopic contribution when the average current is suppressed.

In recent years the interest in Josephson physics has shifted towards junctions whose elements either include topological materials Sacepe (); Veldhorst (); Bestwick (); Mason (); Finck-PRX (); Kurter-PRB (); Kurter-NC (); Sochnikov (); Stehno () or where topological properties are enabled by a specific design of the hybrid-junction with otherwise conventional materials Frolov (); Shtrikman (); Marcus (); Moler (); Kouwenhoven (). These possibilities and advances motivate our work to investigate how universal mesoscopic effects manifest in topological Josephson junctions that in particular host Majorana states (see review Beenakker-RMP-MF () and references therein). We carry out this analysis in the context of multiterminal devices that were brought into the spotlight of recent theoretical attention with the observation that they can emulate topological matter Riwar (); Eriksson (); Houzet (); HYX1 (); HYX2 (); Qi (); Nazarov (); Deb (); Meyer (), which triggered experimental efforts in realizing these systems in various proximitized circuits Giazotto (); Heiblum (); Finkelstein (); Manucharyan (); Pribiag ().

The rest of the paper is organized as follows. In Sec. II we briefly review symmetry-constrained scattering matrix transport formalism in application to Josephson effect in multiterminal circuits. In Sec. III we apply these methods to two-terminal junction as a benchmark, and then extend our analysis to three- and four-terminal devices where we compute emergent band structure of sub-gap states, investigate their topology, and derive transport characteristics such as transconductance and supercurrent. In Sec. IV we focus our attention on the statistical properties of topological supercurrents and obtain analytical results for variance that takes a universal form and also inherits periodicity of Majorana Josephson effect.

## Ii Scattering matrix formalism

Consider a Josephson junction (JJ) where superconducting (S) terminals are connected through the common normal (N) region, thus forming a multiterminal SNS contact. To keep the presentation simple, we assume that each superconducting lead is coupled by only a single conducting channel in the normal region and both time-reversal and chiral symmetries are broken (unconventional classes D and C Beenakker-RMP-MF ()). Formation of the sub-gap bound states in the JJs is the result of coherent Andreev reflections that describe electron-to-hole conversion at the superconductor-normal interface. In -terminal junctions an elastic scattering event at energy is characterized by a scattering matrix where “” denotes the particle-hole degrees of freedom. In what follows we assume that all leads have the same superconducting gap and normalize all energies in units of . The particle-hole (PH) symmetry is represented by

(1) |

where the antiunitary PH transform falls into two categories . For example, for -wave paring, () in spin-nondegenerate case, and () in spin-degenerate case, where are the Pauli matrices acting in particle-hole space and denotes the complex conjugation. The Andreev bound state energies are determined by the determinant equation Beenakker-PRL91 ()

(2) |

Here is the scattering matrix of the normal region, and is the scattering matrix describing Andreev reflections, where is the diagonal matrix of superconducting phases. We set owing to global gauge invariance. Due to the PH symmetry Eq. (1) these scattering matrices take the block-diagonal forms

(3) |

where . The determinant in Eq. (2) simplifies further to a degree- characteristic polynomial of ,

(4) |

which is (anti)palindromic . Importantly, from Eq. (4) we observe that, for a fixed normal-region scattering matrix , the Andreev bands of symmetry classes are dual via the relation

(5) |

Previously, we have extensively discussed the scenario in application to three- and four-terminal junctions HYX1 (); HYX2 (). In this work we primarily focus on the Andreev spectrum of junctions that can support zero-energy Majorana modes. In what follows, we also assume energy-independent scattering matrices that corresponds to, for example, weak links where length of the junction is small compared to the superconducting coherence length, , so that retardation effects of traveling quasiparticles can be neglected. We note that the existence of Majorana zero modes does not depend on this assumption.

To study energy spectra of emergent states in junctions with terminals, we introduce the scattering matrix at ,

(6) |

that belongs to the circular real ensemble since . Via Eq. (2) the zero-energy Majorana modes are determined by the determinant equation of an antisymmetric matrix ,

(7) |

From here we draw important properties. (i) For , Eq. (7) is generally satisfied for any scattering matrices and phases . This implies that Andreev-Majorana zero modes are present at any phases and robust to elastic scattering and superconducting order parameter nonuniformity. These nondispersive flat bands do not contribute to Josephson currents. (ii) For , the Andreev bands cross at zero energy at phases determined by the Pfaffian equation

(8) |

Based on our study on two- and four-terminal junctions, we conjecture that there always exist a pair of Majorana zero-modes modes on an -dimensional hypersurface in the space described by Eq. (8). Next we reveal the energy spectrum of the junction for several concrete forms of the scattering matrix.

## Iii Multiterminal Josephson effect

### iii.1 Two-terminal junctions

We first study two-terminal junctions as a benchmark. We parametrize the unitary matrix by four independent parameters,

(9) |

where representing the normal-region transmission and scattering phases . The sub-gap spectrum of excitations is determined by the characteristic polynomial (4) via equation , where the -function take the form with and . The two branches of dispersive solutions are given by

(10) |

where for comparison we remind the results for the conventional junctions. For finite transmission , the Pfaffian equation (8) reduces to so that a Majorana crossing occurs at with . The zero-temperature Josephson current takes the form

(11) |

Typical energy dispersion and supercurrent-phase relation are shown in Fig. 1. We note that for double-degenerate Majorana states emerge at and the energy and supercurrent exhibit periodicity in . In addition, the bound states are detached from the continuum with a minimal gap at and . Equations (10) and (11) are consistent with the prior results, e.g. Ref. Kwon ().

### iii.2 Three-terminal junctions

For the spectrum of localized states is determined by the palindromic polynomial and composed of three bands,

(12) |

Adopting the same parametrization of the scattering matrix as in Ref. HYX1 () the -function can be found in the closed analytical form

(13) |

Consequently, there are only six independent parameters of the scattering matrix that enter the spectrum of Andreev bound states (ABS). Furthermore, scattering phases only shift the phases of the leads .

Depending on the choice of scattering matrix parameters we find rich behavior of the energy bands. For a special case and the spectrum exhibits nontrivial topology, as shown in Fig. 2. Zero-energy Weyl points appear at for [Fig. 2(b)-(d)]. As shown in Fig. 2(a), the Chern number of corresponding band structure exhibits a sign jump for . We also note that the other topological phase transitions for and are related to the gap closing/reopening at the Andreev band edge . Figure 2(e) displays Josephson currents in two terminals when the system is tuned to the nodal gapless states. In Fig. 2(f) the series of one-dimensional cuts in either phase or show how Josephson currents change as one tunes to the vicinity of nodal points. We observe that moving across the node currents exhibit discontinuous jumps. Note that the Majorana flat band does not contribute to the Josephson current.

### iii.3 Four-terminal junctions

The energy spectrum of four-terminal junctions can host Majorana zero modes and Weyl nodes simultaneously. The four Andreev bands determined by the palindromic equation are given explicitly by the following expressions

(14) |

where the - and -functions are defined by

(15) |

Here we have used short-hand notations for phases , , permutations , and denotes the real part of a complex number. Additionally parameters and are functions of the scattering matrix elements . Specifically

(16) |

and

(17) |

with . The scattering matrix is parametrized by sixteen real parameters as in Ref. Dita (), where and . An inspection of these expressions reveals that despite the fact that we need ten independent phases to parametrize the scattering matrix only six effective angles , , , , , affect the Andreev spectrum in Eq. (14). The zero-energy states are determined by the Pfaffian equation (8),

(18) |

Via the unitary condition of , Eq. (18) implies that where and are real functions of . Most importantly, this determines a Majorana-crossing surface in space given by .

As a practical example, we study the energy bands of this model for the choice of incommensurate parameters: , , , , , , , , , and . The energy spectrum, corresponding Chern number, and Josephson currents are shown in Fig. 3. We observe that the lower bands exhibit -periodicity due to the Majorana crossings described by Eq. (18). Moreover, at [Fig. 3(b)] and [Fig. 3(c)] finite-energy Weyl nodes form at and respectively between one of the higher and lower bands. The appearance of these nodal points is signaled by a change of the Chern number [Fig. 3(a)]. At this point we comment that it has been recently shown that such Majorana-Weyl crossings occur in a different model of a four-terminal junction formed between the end-states of one-dimensional topological superconductors (TS) of class D Meyer (). It has been pointed out that a finite Chern number in this regime is associated with a quantized transconductance . We confirm this result in our scattering matrix model and remark that the extra phase transitions in Fig. 3(a) are related to gap closing/reopening at the band edge that may not be stable since the higher bands can strongly hybridize with the continuum .

The one-dimensional cut of the spectrum in Fig. 3(b) along at is shown in Fig. 3(d). The Josephson currents as functions of corresponding to the spectrum in Fig. 3(b) are shown in Fig. 3(e), where the hedgehog-like singularities are present at the Weyl nodes and . We note that the other two nodal points at and do not induce current singularities since the higher- and lower-band contributions cancel each other. This is can be observed in Fig. 3(f) which displays along the cut at .

## Iv Fluctuations in topological junctions

In the previous section we studied Josephson current for the given realization of the scattering matrix. As alluded in the introduction, this current is expected to display reproducible sample-to-sample fluctuations and it is thus of interest to study its statistical properties. We primarily focus on its variance and also on the cross-correlation function that can be experimentally accessed in the multiterminal devices. As is known from quantum transport theoretical approaches, statistical transport properties of phase-coherent mesoscopic systems can be conveniently computed by means of averaging over a random-matrix that describes the system. In open systems, the averaging is done over the scattering matrix and one typically considers two models of junctions: chaotic cavities or disordered contacts. The former case is more suitable for the model considered in this work. We thus follow classical works by Baranger and Mello Baranger (), and Jalabert, Pichard, and Beenakker Jalabert () who studied conduction through a chaotic cavity on the assumption that the scattering matrix is uniformly distributed in the unitary group, restricted only by symmetry. This is the circular ensemble, introduced by Dyson, and shown to apply to a chaotic cavity by Blumel and Smilansky Smilansky (). In other words we consider SNS junction where normal region is a chaotic quantum dot Whisler ().

The probability density , an invariant Haar measure, of the -matrix parameters is given by

(19) |

The distribution function of an observable , defined as , is in practice calculated by the characteristic function

(20) |

where denoting the CUE ensemble average.

For benchmark we first study the statistics of the Josephson current in the two-terminal junctions and take as the units of in the following discussion. From Eqs. (9) and (19) we obtain a constant invariant measure . We remind that this simplicity is specific to unitary case, for instance in orthogonal symmetry probability density is not flat in even for a single channel limit. All the moments as well as the distribution function of the Josephson current can be obtained analytically. The -moment is given by the expression

(21a) | ||||

(21b) |

where and is the hypergeometric function. Therefore, the variance for reads

(22a) | ||||

whereas for the non-topological case is has a different look | ||||

(22b) |

as depicted in Fig. 4(a). In the topological regime variance inherits periodicity and has remarkably simple form. In the non-topological regime, our result is similar to that of Chalker-Macêdo Chalker () albeit the different numerical coefficients as they considered multi-mode disordered junction model where averaging is done over the Dorokhov distribution of transmission eigenvalues. Finally the Josephson-current distribution function takes the form for

(23a) | |||

and for | |||

(23b) |

where and being the critical currents, is the Heaviside step function, and the function . As shown in Fig. 4(b), in both classes the relation is satisfied. (i) For , is a linear function of for which and the slope is defined by the phase . In particular for with . (ii) For , is smaller for larger current amplitude and for with .

We proceed to study the Josephson-current statistics of the three-terminal junctions. From Eq. (19) we obtain the probability density of the effective parameters ,

(24) |

where is the normalization constant. The numerical results of the expectation value, variance, and covariance of are shown in Fig. 5, where the covariance is defined as . The general relations and are satisfied. We observe that the variances and covariances distinguish the junctions.

For three-terminal junctions, by integrating and terminal we can construct an effective two-terminal S-TS junction Ioselevich () which supports a single channel on one lead (topological), with phase , and two channels on the other (conventional), with phase . Defining and in Eqs. (12) and (13), we obtain the Josephson current through the two leads

(25) |

that depends on five independent parameters . We present the numerical results of the statistical properties of for such a configuration in Fig. 6. The expectation and variance of as functions of are shown in Fig. 6(a). The characteristic function for various is shown in Fig. 6(b). We observe that, for , with the exponent being -independent. We calculate the -averaged characteristic function and fix the exponent . In particular, this analysis enables us to extract the asymptotic behavior of the full distribution function which is found to exhibit a universal power-law scaling in the limit .

## V Discussion and outlook

In this work we applied methods of scattering matrix theory to study transport properties of multiterminal Josephson junctions of topological superconductors. We have examined spectrum of sub-gap states in two-, three-, and four-terminal configurations and determined that texture of resulting Andreev bands in the multidimendisonal parameter space of superconducting phases can produce nonvanishing fluxes of Berry curvature. These properties translate into the quantized nonlocal conductances of these devices. We have also studied current-phase relationships and interaction of supercurrents, as well as their mesoscopic statistical properties. In particular, we discussed universal regime of current fluctuations and computed supercurrent variance as well as current cross-correlation function in the topological regime. We close this work with few comments in relation to existing and future possible experiments where fundamental physics of multiterminal Josephsonic devices could be further explored.

Recently Josephson supercurrent and conductance were measured as a function of geometry, temperature, and gate voltage in proximitized planar junction devices comprised of superconductors and surface states of topological insulator (S-TI-S junctions) in order to determine the nature of the electronic transport in these systems. The supercurrent was found to exhibit a sharp drop as a function of gate voltage, see Figs. (2) and (6) of Ref. Kurter-PRB (), superimposed with reproducible noise whose magnitude was a fraction of the critical current. The systematic trend in the critical current dependence was explained by a mechanism related to the relocation of the topological surface state with respect to trivial conducting two-dimensional states formed by band-banding near the surface. In real space, a negative gating potential pushes the trivial state below the topological surface states, exposing the topological state to the disordered surface of the TI. As a result, the magnitude of the supercurrent changes sharply. The noise was attributed to the percolation effects as near the voltage threshold it is likely that local charge fluctuations cause the path of the supercurrent to be highly meandering. We wish to point out that there is possibly for an alternative picture as this noise could be of mesoscopic origin. This evidence is further supported by observed similar reproducible noise features in Fraunhofer magneto-oscillations of the critical current. While our model is not directly applicable to S-TI-S junctions we draw an observation the magnitude of current fluctuations is consistent with the expectations that disorder scattering causes observed mesoscopic effects.

In addition, we wish to comment that related statistical properties of supercurrents can be also studied by measuring switching current distributions. In particular, for topological Josephson devices, the critical current measurements can potentially enable determining the parity state of a Majorana fermion (pair) in a junction since supercurrent acquires an anomalous fractional component due to Majorana modes, , where the sign encodes the parity. The typical switching measurement is performed by ramping the bias current through the junction to detect the current value at which the junction jumps to the finite voltage state. By repeating this protocol many times and accumulating statistics of random supercurrent switching events (as previously successfully implemented in various mesoscopic proximity circuits, e.g. nanowires and graphene layers Coskun (); Murphy ()), one expects to reveal a bimodal distribution indicating the two parity states. If the separation of the two peaks in this distribution is wide enough, one can detect (with some fidelity) the parity state. Mutiterminal devices considered in this work can provide an actual hardware platform to conduct such experiments and our transport theory will be useful in modeling future measurements. In particular, knowledge of current-phase relationship is needed for determining the energy barrier of a phase slip that triggers the switching. Furthermore, these developments are also inspired by the potential application of multiterminal devices in design of protected superconducting qubits.

## Acknowledgment

We would like to thank Dale Van Harlingen for the discussions, especially on the topic of the gate dependence of the critical current fluctuations in topological Josephson junctions Kurter-PRB () and feasibility of switching current experiments in tri-junction configuration. The work of H.-Y. X was supported by National Science Foundation Grant DMR-1653661. The work of A.L was supported in part by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, under Award DESC0017888, and the Army Research Office, Laboratories for Physical Sciences Grant W911NF-18-1-0115. H.-Y. X. acknowledges the hospitality of the Kavli Institute of Theoretical Sciences (KITS) where this work was partially done.

## References

- (1) B. L. Al’tshuler, Fluctuations in the extrinsic conductivity of disordered conductors, JETP Lett. 41, 648 (1985).
- (2) B. L. Al’tshuler and D. E. Khmel’nitskii, Fluctuation properties of small conductors, JETP Lett. 42, 359 (1986).
- (3) P. A. Lee and A. D. Stone, Universal Conductance Fluctuations in Metals, Phys. Rev. Lett. 55, 1622 (1985).
- (4) P. A. Lee, A. D. Stone, and H. Fukuyama, Universal conductance fluctuations in metals: Effects of finite temperature, interactions, and magnetic field, Phys. Rev. B 35, 1039 (1987).
- (5) I. L. Aleiner and Ya. M. Blanter, Inelastic scattering time for conductance fluctuations, Phys. Rev. B 65, 115317 (2002).
- (6) F. J. Dyson, Statistical Theory of the Energy Levels of Complex Systems, J. Math. Phys. 3, 140 (1962).
- (7) B. L. Altshuler and B. I. Shklovskii, Repulsion of energy levels and conductivity of small metal samples, Sov. Phys. JETP 64, 127 (1986).
- (8) M. L. Mehta, Random Matrices, (Elsevier, 2004).
- (9) S. G. den Hartog, C. M. A. Kapteyn, B. J. van Wees, T. M. Klapwijk, W. van der Graaf, and G. Borghs, Sample-Specific Conductance Fluctuations Modulated by the Superconducting Phase, Phys. Rev. Lett. 76, 4592 (1996).
- (10) K. Hecker, H. Hegger, A. Altland, and K. Fiegle, Conductance Fluctuations in Mesoscopic Normal-Metal/Superconductor Samples, Phys. Rev. Lett. 79, 1547 (1997).
- (11) P. W. Brouwer, On the Random-Matrix Theory of Quantum Transport, PhD Thesis (Leiden, 1997).
- (12) C. W. J. Beenakker, Random-matrix theory of quantum transport, Rev. Mod. Phys. 69, 731 (1997).
- (13) B. L. Al’tshuler and B. Z. Spivak, Mesoscopic fluctuations in a superconductor-normal metal-superconductor junction, Sov. Phys. JETP 65, 343 (1987).
- (14) C. W. J. Beenakker, Universal Limit of Critical-Current Fluctuations in Mesoscopic Josephson Junctions, Phys. Rev. Lett. 67, 3836 (1991).
- (15) J. T. Chalker and A. M. S. Macêdo, Complete Characterization of Universal Fluctuations in Quasi-One-Dimensional Mesoscopic Conductors, Phys. Rev. Lett. 71, 3693 (1993).
- (16) B. Sacepe, J. B. Oostinga, J. Li, A. Ubaldini, N. J. G. Couto, E. Giannini, and A. F. Morpurgo, Gate-tuned normal and superconducting transport at the surface of a topological insulator, Nat. Commun. 2, 575 (2011).
- (17) M. Veldhorst, M. Snelder, M. Hoek, T. Gang, V. K. Guduru, X. L. Wang, U. Zeitler, W. G. Van der wiel, A. A. Golubov, H. Hilgenkamp, and A. Bririkman, Josephson supercurrent through a topological insulator surface state, Nat. Mater. 11, 417 (2012).
- (18) J. R. Williams, A. J. Bestwick, P. Gallagher, Seung Sae Hong, Y. Cui, Andrew S. Bleich, J. G. Analytis, I. R. Fisher, and D. Goldhaber-Gordon, Unconventional Josephson Effect in Hybrid Superconductor-Topological Insulator Devices, Phys. Rev. Lett. 109, 056803 (2012).
- (19) S. Cho, B. Dellabetta, A. Yang, J. Schneeloch, Z. Xu, T. Valla, G. Gu, M. J. Gilbert, and N. Mason, Symmetry protected Josephson supercurrents in three-dimensional topological insulators, Nat. Commun. 4, 1689 (2013).
- (20) A. D. K. Finck, C. Kurter, Y. S. Hor, D. J. Van Harlingen, Phase Coherence and Andreev Reflection in Topological Insulator Devices, Phys. Rev. X 4, 041022 (2014).
- (21) C. Kurter, A. D. K. Finck, P. Ghaemi, Y. S. Hor, and D. J. Van Harlingen, Dynamical gate-tunable supercurrents in topological Josephson junctions, Phys. Rev. B 90, 014501 (2014).
- (22) C. Kurter, A. D. K. Finck, Y. S. Hor, D. J. Van Harlingen, Evidence for an anomalous current-phase relation in topological insulator Josephson junctions, Nat. Commun. 6, 7130 (2015).
- (23) I. Sochnikov, L. Maier, C. A. Watson, J. R. Kirtley, C. Gould, G. Tkachov, E. M. Hankiewicz, C. Brüne, H. Buhmann, L. W. Molenkamp, and K. A. Moler, Nonsinusoidal Current-Phase Relationship in Josephson Junctions from the 3D Topological Insulator HgTe, Phys. Rev. Lett. 114, 066801 (2015).
- (24) M. P. Stehno, V. Orlyanchik, C. D. Nugroho, P. Ghaemi, M. Brahlek, N. Koirala, S. Oh, and D. J. Van Harlingen, Signature of a topological phase transition in the Josephson supercurrent through a topological insulator, Phys. Rev. B 93, 035307 (2016).
- (25) V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, L. P. Kouwenhoven, Signatures of Majorana fermions in hybrid superconductor-semiconductor nanowire devices, Science 336, 1003 (2012).
- (26) Anindya Das, Yuval Ronen, Yonatan Most, Yuval Oreg, Moty Heiblum, Hadas Shtrikman, Evidence of Majorana fermions in an Al - InAs nanowire topological superconductor, Nature Physics 8, 887 (2012).
- (27) D. Sherman, J. S. Yodh, S. M. Albrecht, J. Nygard, P. Krogstrup, C. M. Marcus, Hybrid Double Quantum Dots: Normal, Superconducting, and Topological Regimes, Nature Nanotechnology 12, 212 (2017).
- (28) Eric M. Spanton, Mingtang Deng, Saulius Vaitiekenas, Peter Krogstrup, Jesper Nygard, Charles M. Marcus, Kathryn A. Moler, Current-phase relations of few-mode InAs nanowire Josephson junctions, Nature Physics 13, 1177 (2017).
- (29) R. M. Lutchyn, E. P. A. M. Bakkers, L. P. Kouwenhoven, P. Krogstrup, C. M. Marcus, Y. Oreg, Realizing Majorana zero modes in superconductor-semiconductor heterostructures, Nat. Rev. Mater. 3, 52 (2018).
- (30) C. W. J. Beenakker, Random-matrix theory of Majorana fermions and topological superconductors, Rev. Mod. Phys. 87, 1037 (2015).
- (31) Roman-Pascal Riwar, Manuel Houzet, Julia S. Meyer, and Yuli V. Nazarov, Multi-terminal Josephson junctions as topological matter, Nature Communications 7, 11167 (2016).
- (32) Erik Eriksson, Roman-Pascal Riwar, Manuel Houzet, Julia S. Meyer, and Yuli V. Nazarov, Topological transconductance quantization in a four-terminal Josephson junction, Phys. Rev. B 95, 075417 (2017).
- (33) Julia S. Meyer and Manuel Houzet, Non-trivial Chern numbers in three-terminal Josephson junctions, Phys. Rev. Lett. 119, 136807 (2017).
- (34) Hong-Yi Xie, Maxim G. Vavilov, Alex Levchenko, Topological Andreev bands in three-terminal Josephson junctions, Phys. Rev. B 96, 161406(R) (2017).
- (35) Hong-Yi Xie, Maxim G. Vavilov, Alex Levchenko, Weyl nodes in Andreev spectra of multiterminal Josephson junctions: Chern numbers, conductances, and supercurrents, Phys. Rev. B 97, 035443 (2018).
- (36) Zhenyi Qi, Hong-Yi Xie, Javad Shabani, Vladimir E. Manucharyan, Alex Levchenko, and Maxim G. Vavilov, Controlled-Z gate for transmon qubits coupled by semiconductor junctions, Phys. Rev. B 97, 134518 (2018).
- (37) Janis Erdmanis, Arpad Lukacs, and Yuli V. Nazarov, Weyl disks: Theoretical prediction, Phys. Rev. B 98, 241105(R) (2018).
- (38) Oindrila Deb, K. Sengupta, and Diptiman Sen, Josephson junctions of multiple superconducting wires, Phys. Rev. B 97, 174518 (2018)
- (39) Manuel Houzet and Julia S. Meyer, Majorana-Weyl crossings in topological multi-terminal junctions, preprint arXiv:1810:09962.
- (40) E. Strambini, S. D’Ambrosio, F. Vischi, F. S. Bergeret, Yu. V. Nazarov, and F. Giazotto, The -SQUIPT as a tool to phase-engineer Josephson topological materials, Nat. Nanotechnol. 11, 1055 (2016).
- (41) Yonatan Cohen, Yuval Ronen, Jung-Hyun Kang, Moty Heiblum, Denis Feinberg, Regis Melin, Hadas Shtrikman, Non-local Supercurrent of Quartets in a Three-Terminal Josephson Junction, PNAS 115 (27), 6991 (2018).
- (42) Anne W. Draelos, Ming-Tso Wei, Andrew Seredinski, Hengming Li, Yash Mehta, Kenji Watanabe, Takashi Taniguchi, Ivan V. Borzenets, Francois Amet, Gleb Finkelstein, Supercurrent Flow in Multi-Terminal Graphene Josephson Junctions, preprint arXiv:1810.11632.
- (43) Natalia Pankratova, Hanho Lee, Roman Kuzmin, Maxim Vavilov, Kaushini Wickramasinghe, William Mayer, Joseph Yuan, Javad Shabani, Vladimir E. Manucharyan, The multi-terminal Josephson effect, preprint arXiv:1812.06017.
- (44) Gino Graziano, Mihir Pendharkar, Joon Sue Lee, Chris Palmstrøm, and Vlad S. Pribiag, Coherent transport in multi-terminal Josephson junctions, under preparation.
- (45) H.-J. Kwon, K. Sengupta, and V. M. Yakovenko, Fractional ac Josephson effects in p- and d-wave superconductors, Eur. Phys. J. B 37, 349 (2004).
- (46) P. Dita, Parametrisation of unitary matrices, J. Phys. A: Math. Gen. 15, 3465 (1982).
- (47) Harold U. Baranger and Pier A. Mello, Mesoscopic Transport through Chaotic Cavities: A Random S-Matrix Theory Approach, Phys. Rev. Lett. 73, 142 (1994).
- (48) R. A. Jalabert, J.-L. Pichard, and C. W. J. Beenakker, Universal Quantum Signatures of Chaos in Ballistic Transport, Europhys. Lett. 27, 255 (1994).
- (49) R. Blumel and U. Smilansky, Classical irregular scattering and its quantum-mechanical implications, Phys. Rev. Lett. 60, 477 (1988).
- (50) Colin M. Whisler, Maxim G. Vavilov, Alex Levchenko, Josephson currents in chaotic quantum dots, Phys. Rev. B 97, 224515 (2018).
- (51) P. A. Ioselevich, P. M. Ostrovsky, and M. V. Feigel’man, Josephson current between topological and conventional superconductors, Phys. Rev. B 93, 125435 (2016).
- (52) U. C. Coskun, M. Brenner, T. Hymel, V. Vakaryuk, A. Levchenko, and A. Bezryadin, Distribution of Supercurrent Switching in Graphene under the Proximity Effect, Phys. Rev. Lett. 108, 097003 (2012).
- (53) A. Murphy, P. Weinberg, T. Aref, U. C. Coskun, V. Vakaryuk, A. Levchenko, and A. Bezryadin, Universal Features of Counting Statistics of Thermal and Quantum Phase Slips in Nanosize Superconducting Circuits, Phys. Rev. Lett. 110, 247001 (2013).