Arxiv:2009.11768V1 [Cond-Mat.Mes-Hall] 24 Sep 2020

Ground-state quantum geometry in superconductor–quantum dot chains

R. L. Klees,1 J. C. Cuevas,2 W. Belzig,1 and G. Rastelli1, 3, 4 1Fachbereich Physik, UniversitätKonstanz, D-78457 Konstanz, Germany 2Departamento de F´ısica Teórica de la Materia Condensada and Condensed Matter Physics Center (IFIMAC), Universidad Autónomade Madrid, E-28049 Madrid, Spain 3Zukunftskolleg, UniversitätKonstanz, D-78457 Konstanz, Germany 4INO-CNR BEC Center and Dipartimento di Fisica, Universitàdi Trento, I-38123 Povo, Italy (Dated: September 9, 2021) Multiterminal Josephson junctions constitute engineered topological systems in arbitrary synthetic dimensions defined by the superconducting phases. Microwave spectroscopy enables the measurement of the quantum geometric tensor, a fundamental quantity describing both the quantum geometry and the topology of the emergent Andreev bound states in a unified manner. In this work we propose an experimentally feasible and scalable multiterminal setup of N quantum dots connected to N + 1 superconducting leads which allows to deterministically study nontrivial topology in terms of the Chern number of the noninteracting ground state. An important result is that the nontrivial topology in a linear chain appears beyond a threshold value of the nonlocal proximity-induced pairing potential which represents the novel theoretical key ingredient of our proposal. Moreover, we generalize the microwave spectroscopy scheme to the multiband case and show that the elements of the quantum geometric tensor of the noninteracting ground state can be experimentally accessed from the measurable oscillator strengths at low temperature.

I. INTRODUCTION mented by simply increasing the number of superconducting leads and that building blocks can be conven- The concept of topology plays an important role in tional materials, although also topological superconduc- modern branches of theoretical physics by explaining the tors hosting Majorana zero modes have been studied in quantized nature of certain physical phenomena. The this context [27, 32–34]. Moreover, it has been recently pursue to find novel types of topologically nontrivial sys- suggested to use microwave spectroscopy to measure the tems providing quantized observables classified by inte- more fundamental quantum geometric tensor of ABS, gers has created a new research field called topological which provides both the information about the geome- band theory [1]. Since then, a huge variety of topological try of the state manifold and the topological information quantum matter has been predicted and discovered such contained in the Berry curvature [30]. While the Chern as topological insulators [2,3], topological semimetals [4], number follows directly from the local Berry curvature topological superconductors [5], as well as non-Hermitian by a simple integration, the local metric tensor carries (e.g., open/dissipative) topological systems [6]. Out of important information on, for instance, energy fluctu- these, the superconducting systems are also promising ations and the noise spectral function [36, 37], as well platforms for topologically protected quantum computa- as quantum phase transitions [38], wave packet dynam- tion that relies on non-Abelian exchange statistics, i.e., ics [39, 40], the superfluid weight [41–43], energy correc- braiding, of Majorana zero modes appearing at the Fermi tions of excitons [44], the orbital magnetic susceptibility energy [7–9] and their appearance was recently verified [45–47], and magnetic exchange constants [48]. Finally, experimentally [10]. Furthermore, qubits based on the it is also possible to strongly drive topologically trivial Andreev bound states (ABS) appearing inside a Joseph- Josephson systems to eventually generate nontrivial (Flo- son junction have also been proposed [11, 12]. These quet) topology [49, 50]. bound states can be experimentally accessed and coher- So far, many experiments to measure quantum geome- ently manipulated by microwave [13–16], tunneling [17], try (i.e., the Berry curvature and/or the quantum metric) and supercurrent spectroscopy [18] if the junction is em- and the related topology have already been successfully bedded in an rf superconducting quantum interference carried out in, e.g., ultracold fermionic atoms [51, 52], device (SQUID). superconducting qubits and qutrits [53, 54], supercon- arXiv:2009.11768v2 [cond-mat.mes-hall] 8 Sep 2021 More recently, multiterminal Josephson junctions ducting transmon qubits [55], superconducting quantum (MJJs) consisting of many superconducting terminals circuits [56], qubits in NV centers in diamond [57], high- have been theoretically investigated and shown to exhibit finesse planar microcavities [58], and coupled fiber loops topologically nontrivial physics [19–35]. In such multiter- [59], the last two representing examples of the emerging minal systems, topology emerges in the synthetic space field of topological photonics [60]. of superconducting phases and the integer-valued Chern Additional interesting directions are topological number can manifest itself in a quantized transconduc- fermion-parity [61] and Cooper-pair pumps [62–66], as tance between two terminals [22–28, 32, 34]. The ad- well as Josephson tunnel circuits [67–69], in which Weyl vantage of these systems is that, in principle, an ar- singularities in parameter space are responsible for a bitrary number of synthetic dimensions can be imple- quantized adiabatic transfer of single charges (charge 1e) 2

FIG. 1. Model of the superconductor-quantum dot chain. The chain consists of N quantum dots (each with onsite energy εj , j = 1,...,N) which are connected to N + 1 s-wave superconducting leads (each with a pairing phase ϕk, k = 0,...,N). The couplings between the quantum dots are wl,l+1 (l = 1,...,N − 1), while the couplings between the quantum dots and the superconducting leads are tj,j−1 and tjj . and Cooper pairs (charge 2e). on the topological phases of this system. In Sec.III, we Although recently a diffusive three-terminal MJJ in discuss the application of microwave spectroscopy in such a double-SQUID configuration has been both experi- systems. Importantly, we find that, although the quan- mentally and theoretically investigated on its topologi- tum geometry of the individual ABS is inaccessible at cal real-space properties [70] followed by a detailed theo- zero temperature, the topology and the quantum geom- retical study [71], it is still challenging to create MJJs etry of the noninteracting GS formed by many occupied [17, 72]. Furthermore, despite the fact that first ex- fermionic states can be accessed even in the multiband periments towards ballistic MJJs have been performed case. We discuss our findings and conclude our work in [73–75], the more general theoretical proposals of MJJs Sec.IV. Further details of derivations can be found in discussed above are expected to be quite challenging to the AppendixesA,B, andC. realize since the necessary space for experimental control gates decreases as the number of superconducting terminals increases. Such circular arrangements of the II. EFFECTIVE LOW-ENERGY MODEL AND terminals will eventually have an adverse effect on the TOPOLOGY scalability of these systems. In this work, we propose a linear design of an MJJ as an As described above, we start by considering a lin- alternative and scalable device to enable the determinis- ear chain of N + 1 BCS superconducting leads con- tic study quantum geometry and topology in topological nected to N quantum dots as shown in Fig.1. Each Josephson matter. This setup consists of a chain of quan- of the quantum dots is spin-degenerate with an onsite tum dots which are coupled to superconducting leads, as energy εj (j = 1,...,N) and each s-wave superconduc- presented in Fig.1. We think that this system is easier tor has a pairing potential ∆k with a pairing phase ϕk to realize and manipulate in experiments since there is (k = 0,...,N). The dot j is coupled to two neighbor- enough space from below for additional necessary control ing superconducting leads j − 1 and j via the hoppings gates and it could be implemented with already available tj,j−1 and tjj, respectively. To be more general, we in- ingredients such as carbon nanotubes or semiconduct- clude the effect of a direct tunneling between the dots ing nanowires contacted by superconductors [17, 76–80]. without passing through the superconductors, which is Moreover, we expect that it will be straightforward to described by a direct coupling wl,l+1 (l = 1,...,N − 1) extend the length of the chain and, hence, increase the between the dots l and l + 1. By integrating out the sub- number of dots and leads at will. We will show that such system of the superconducting leads and by considering a system allows for the emergence of topologically non- the limit of large superconducting gaps ∆k → ∞, we ob- trivial ABS and, therefore, represents an ideal platform tain an effective low-energy Hamiltonian describing the to experimentally study quantum geometry and topol- quantum dots proximity coupled to the superconductors. ogy in a deterministic way, in contrast to prior proposals The derivation of the effective Hamiltonian of our model based on statistical methods [22, 24, 26]. We will also is presented in AppendixA. show how to apply microwave spectroscopy in such multi- The resulting Hamiltonian has the form Heff = † ˆ (N) band systems to measure the quantum geometry of the d Heff d, where the 2N × 2N block-tridiagonal matrix noninteracting ground state (GS) as a whole, which is im- Hamiltonian is given by portant for (quantized) quantum transport phenomena  ˆ ˆ  due to the collection of occupied single-particle states. H1 V12 0 ··· 0 0 Vˆ Hˆ Vˆ ··· 0 0  The rest of the paper is organized as follows: In Sec.II,  12 2 23   0 Vˆ Hˆ ··· 0 0  we introduce and investigate the effective low-energy ˆ (N)  23 3  Heff =  ......  , (1) Hamiltonian of a linear chain of N dots coupled to N + 1  ......  superconductors. We also define the topological invariant    0 0 0 ··· Hˆ Vˆ  of the noninteracting GS, i.e., the Chern number, and an-  N−1 N−1,N  ˆ ˆ alyze the effects of the various parameters of this model 0 0 0 ··· VN−1,N HN 3 which is written in the basis defined by N = 3 dots in Fig. 2(a), which also shows that this the 2N-dimensional Nambu spinor d† = system allows for the appearance of Weyl singularities, † † † (†) i.e., points of degeneracy ϕ ∈ 3 at which two energy (d1↑, d1↓, d2↑, d2↓, . . . , dN↑, dN↓), with djσ being the W T annihilation (creation) operator of an electron with spin bands cross [Fig. 2(b)]. In this case, we have set the en- σ ∈ {↑, ↓} on the j-th dot. The 2 × 2 Nambu blocks are ergies of the quantum dots equal in magnitude with an given by alternating sign. However, Weyl singularities will also appear in more general situations. As we will see below, ˆ εj zj ˆ wl,l+1 zl,l+1 Weyl points for which εA,±1(ϕW) = 0 are responsible for Hj = ∗ , Vl,l+1 = ∗ , (2) zj −εj zl,l+1 −wl,l+1 topological phase transitions between regions of different values of the Chern number of the GS of the system. where the effective local and nonlocal proximity-induced The superconducting phases ϕ define N synthetic U(1) pairing potentials on the dots are defined as gauge fields for which we define gauge connection 1- (n) PN (n) (j−1) iϕj−1 (j) iϕj forms A = aα dϕ for each of the ABS |ψ i, zj = Γ e + Γ e , (3a) α=1 α n j j (n) where aα = ihψ |∂ ψ i is the (Abelian) Berry connec- (l) iϕl n α n zl,l+1 = Γ e , (3b) l,l+1 tion [83] and ∂α ≡ ∂/∂ϕα. The gauge-invariant cur- (n) (k) (l) vature 2-form F of the n-th single-particle state is 2 (n) (n) with Γj = πN0tjk and Γl,l+1 = πN0pltlltl+1,l, respec- defined as the exterior derivative of A , i.e., F = tively, where N0 is the normal density of states at the dA(n) = (1/2) PN f (n)dϕ ∧ dϕ , where f (n) = Fermi energy. In passing by, we modified the nonlocal α,β=1 αβ α β αβ ∂ a(n) − ∂ a(n) is the (Abelian) Berry curvature [84]. pairings by dimensionless factors pl ∈ [0, 1] motivated α β β α q Finally, the first Chern numbers of an ABS |ψ i are cal- by experiments in which typically Γ(l) < Γ(l)Γ(l) n l,l+1 l l+1 culated as [78, 79]. This modification allows us to treat the nonlocal couplings as independent parameters from the local Z 2π Z 2π ones. In general, the parameters p modifying the nonlo- (n) 1 (n) l c = f dϕα dϕβ , (4) (l) αβ αβ 2π 0 0 cal couplings Γl,l+1 depend on the geometrical details of the contacts as well as the coherence length of the Cooper pairs [81, 82]. which are defined by an integration in the (ϕα, ϕβ) sub- Finally, as in a system of N + 1 superconducting leads plane of the N-torus TN . Note that, due to PH sym- only N superconducting phases are independent, gauge- (−j) (j) (−j) metry, we have aα = −aα and, consequently, f = invariance allows us to set one of the superconducting αβ −f (j) and c(−j) = −c(j). We use the numerical algorithm phases to zero. We will consider ϕ0 = 0 in the fol- αβ αβ αβ lowing as a reference and the remaining N supercon- developed by Fukui, Hatsugai, and Suzuki in Ref. [85] for N ducting phases ϕ = (ϕ1, . . . , ϕN ) ∈ [0, 2π) define an a stable and gauge-independent calculation of the Chern N-dimensional periodic compact space TN (N-torus) in numbers. analogy to the quasimomenta in a periodic crystal with While the topological properties of an ABS |ψni in a first Brillouin zone. (n) terms of the Chern numbers cαβ is encoded in its local For N quantum dots coupled to N +1 superconducting (n) leads in the low-energy limit, there are 2N ABS. Defining Berry curvature fαβ , the full quantum geometry of each the sets of labels D± = {±1,..., ±N} for occupied (D−) of these ABS is described by the local quantum geometric and empty (D+) states at zero temperature, the 2N ABS tensor (QGT) [86] − + ˆ (N) |ψni (n ∈ D ∪ D ) satisfy Heff |ψni = εA,n|ψni, where (n) εA,n is the energy of the n-th ABS. In the following, we qαβ = h∂αψn|(1 − |ψnihψn|)|∂βψni . (5) will consider the ABS being ordered as εA,−N ≤ ... ≤ εA,−1 ≤ 0 ≤ εA,1 ≤ ... ≤ εA,N without loss of generality. The 2N ABS come in N particle-hole (PH) symmet- In particular, we obtain the (Fubini-Study) metric tensor ric pairs and fulfill |ψ i = P (N)|ψ i (j ∈ D+), where (n) (n) −j j gαβ = Re(qαβ ) from the real part of the QGT that pro- (N) (N) P = i(11 ⊗τ2)K is the operator of PH symmetry, K vides a measure of distance ds2 = PN g(n)dϕ dϕ (N) α,β=1 αβ α β is complex conjugation, 11 is the N × N identity ma- on the n-th ABS manifold of physical states, while the trix, and τ2 is the second Pauli matrix in Nambu space. Berry curvature f (n) = −2 Im(q(n)) is obtained from the Furthermore, the energies of a PH-symmetric pair are αβ αβ (N) imaginary part of the QGT that describes the geometric related by εA,−j = −εA,j since P satisfies the anti- phase of a physical quantum state. commutation relation {Hˆ (N),P (N)} = 0. eff The noninteracting GS of the system is defined by the Due to the large number of parameters, we choose (k) collection of all occupied single-particle states below the equal local couplings Γ = ΓL, equal nonlocal couplings j Fermi energy, i.e., the N ABS with energies εA,j ≤ 0 for (l) − Γl,l+1 = ΓNL, and equal interdot couplings wl,l+1 = w j ∈ D . As shown in AppendixB, the metric tensor for simplicity. As an example, we show the ABS for and the Berry curvature of the noninteracting GS can be 4

(a) (c) 1.0 0.9 0.8 1 0.7 0.66 0.6 0 0.5 0.33 0.4 -1 0.3 (b) 3 3 0

2 -3

-6

1 -9

-12 0 012 012012

FIG. 2. (a) ABS energies εA,−3, . . . , εA,3 (from bottom to top) for N = 3 dots. Energy bands with the same color represent PH-symmetric ABS with energies εA,±j . Parameters: ε1 = −ε2 = ε3 = −ΓL/2, ϕ0 = 0, ϕ3 = π, w = 0.2 ΓL,ΓNL = 0.8 ΓL. 3 (b) Locations of the four Weyl points ϕW ∈ T for which εA,±1(ϕW) = 0 for the same parameters as in panel (a). The color indicates the topological charge of the Weyl points. The two blue Weyl points in the ϕ3 = π plane are visible as zero-energy (GS) band crossings in panel (a). (c) Topological phase diagrams of the ground state Chern number C12 for N = 3 dots for different values of the direct tunneling w between the dots (columns) as a function of the phase ϕ3. The upper row shows the phase diagram as a function of the nonlocal coupling ΓNL for fixed ε2 = ΓL/2, while the lower row shows it as a function of the central dot level ε2 for fixed ΓNL = 0.8 ΓL. All other parameters are the same as in panel (a). The dashed line with blue and red dots in the phase diagrams for w = 0.2 ΓL mark the locations of the planes with fixed ϕ3 in which the Weyl points of the same color appear in panel (b). constructed from the single-particle contributions as numbers from, e.g., 0 to ±1, the sum of them entering C(GS) in Eq. (7) will still be unchanged. (GS) X (j) X αβ G = g − h∂ ψ |ψ ihψ |∂ ψ i , (6a) αβ αβ α j k k β j In order to obtain topologically nontrivial states in − − j∈D j,k∈D (GS) j6=k terms of Cαβ , the minimal number of quantum dots (GS) X (j) in this system is N = 3 such that there are three tun- Fαβ = fαβ . (6b) j∈D− able phase differences, as has already been discussed in Refs. [21–24, 28, 30]. To analyze and discuss the stability (GS) of the topological phases for this case, we show topologi- While the result for the Berry curvature Fαβ is already well-known, the counter-intuitive result for the metric cal phase diagrams in Fig. 2(c) in terms of the GS Chern (GS) (GS) tensor G shows that it is not a simple sum over the number C12 . We see that there are large and stable αβ (GS) metric tensors of the individual single-particle states. regions for which nontrivial regions of C12 6= 0 appear. As a direct consequence of Eqs. (4) and (6b), the Chern This stability is due to the separation of Weyl points number of the noninteracting GS is simply given by the ϕW of different topological charge [Fig. 2(b)]. These sum of all individual Chern numbers of the occupied Weyl points move continuously in T3 while the param- single-particle states, i.e., eters of the system are changed, making these regions robust against small fluctuations. As long as two Weyl (GS) X (j) points of different topological charge do not meet and Cαβ = cαβ. (7) j∈D− annihilate, topological regions will be present. As shown in Fig. 2(c), the first important result is that (GS) Note that Cαβ only changes if there is a zero-energy the nontrivial topology in a linear chain appears only in crossing between the two PH-symmetric ABS |ψ−1i and the presence of a nonlocal pairing potential which rep- |ψ1i. Although there are of course Weyl points for which resents the novel theoretical key ingredient of our pro- other energy bands cross changing the individual Chern posal. A nonzero Chern number appears only beyond a 5 certain threshold value of the nonlocal pairing coupling ΓNL. Interestingly, we see that it is not necessary that the quantum dots are directly coupled (w = 0) for the Chern number to be nontrivial. For a small coupling w = 0.2 ΓL we still find topologically nontrivial regions. Remarkably, a larger interdot coupling of w = 2 ΓL can also enhance the size of topological regions in parameter space, although the two regions of finite Chern number are interchanged. However, if the interdot coupling becomes too strong (w ΓL, not shown), topological regions completely disappear since the effect of the superconductors will be negligible. In addition, we see that the nonlocal couplings ΓNL can be considerably smaller than the local couplings ΓL and that there is a lot of freedom in experimental tunability of, e.g., the central dot level ε2. In particular, the lower right plot in Fig. 2(c) shows that it is not required for the energy levels εj of the quantum dots to have the same magnitude with alternating sign in order to obtain nontrivial topological phases. For w = 2 ΓL and ΓNL = 0.8 ΓL, ε2 can be de- tuned even up to ε2 ≈ −13 ΓL with topological regions (1m) FIG. 3. Different contributions of the elements qαβ from still being present. (1) other ABS with energies εA,m (m 6= 0, 1) to the QGT qαβ of the ABS with energy εA,1 according to Eq. (8). III. QUANTUM GEOMETRY AND MICROWAVE SPECTROSCOPY ϕβ → ϕβ + 2A cos(ωt − γ)/~ω to another phase (α 6= β) In the following, we aim at applying the method of mi- with a fixed phase difference γ between the two mod- (nm) crowave spectroscopy described in Ref. [30] in the pres- ulations, to measure the off-diagonal elements qαβ of ence of more than one nondegenerate pair of ABS. The (γ) the QGT. The resulting transition rates are Rn→m,αβ = single-particle QGT in Eq. (5) of each ABS |ψ i can be n r(γ) δ(ε − ε − ω) with the oscillator strength decomposed into the sum n→m,αβ m n ~ [30, 87] N (n) X (nm) q = q (8) 2 αβ αβ (γ) 2πA n (nm) (nm) m=−N rn→m,αβ = qαα + qββ m6=0,n ~ iγ (nm) −iγ (nm)o (nm) + e qαβ + e qβα , (10) of nonadiabatic transition elements qαβ = h∂ ψ |ψ ihψ |∂ ψ i. α n m m β n which now depends on the relative phase difference γ be- In general, Eq. (8) implies a priori that we need to tween the two drives. Details of the derivation of Eqs. (9) measure all possible nonadiabatic transitions between a and (10) are presented in AppendixC. By comparing certain state |ψ i and all the other states |ψ i (n 6= m) n m Eqs. (9) and (10) with the definition of the QGT in by microwave spectroscopy to construct the QGT of this Eq. (8), we obtain the relations particular single-particle state |ψni, as sketched in Fig.3. We will first describe the general procedure and later N X discuss possible experimental difficulties. g(n) = ~ r , (11a) αα 2πA2 n→m,αα On the one hand, in order to measure the diago- m=−N (nm) m6=0,n nal elements qαα of the QGT, we need to apply a N small perturbation (A/ ω 1) to one phase ϕα → (n) X (0) (π) ~ g = ~ r − r , (11b) ϕα + 2A cos(ωt)/~ω, resulting in the transition (photon αβ 8πA2 n→m,αβ n→m,αβ m=−N absorption) with a rate Rn→m,αα = rn→m,αα δ(εm −εn − m6=0,n ~ω) with the oscillator strength [30, 87] N (n) X (π/2) (−π/2) f = ~ r − r . (11c) 2π αβ 4πA2 n→m,αβ n→m,αβ 2 (nm) m=−N rn→m,αα = A qαα (9) ~ m6=0,n according to Fermi’s golden rule. On the other hand, This allows us to measure the metric tensor by linear mi- we assume to simultaneously apply two small pertur- crowave spectroscopy, i.e., modulating one phase or two bations, ϕα → ϕα + 2A cos(ωt)/~ω to one phase and phases with γ = 0, π, while the Berry curvature follows 6 from circular microwave spectroscopy, i.e., modulating ≤ -5 0 ≥ 5 two phases with γ = ±π/2. 2 Although this seems to be straightforward to imple- ment, the experimental difficulty at very low temperature is to measure the transitions between two occupied or two empty states, respectively. To overcome this ob- stacle, one could either perform the experiment at suf- 1 ficiently high temperatures, for which previously filled (empty) states get partially depopulated (populated), or one could use pulses to empty/fill a certain state before measuring the transition. 0 Fortunately, for low-temperature quantum transport, 2 the experimentally relevant quantity is the QGT of the (GS) (GS) (GS) GS Qαβ = Gαβ − iFαβ /2 which, in the noninteracting case, can be constructed from the elements in Eq. (11) by means of Eq. (6). The resulting Berry curvature and 1 (GS) (GS) metric tensor of the GS, Fαβ and Gαβ , respectively, are still accessible via microwave spectroscopy even at zero temperature and the result simply reads (see Ap- pendixC for details) 0 0120 1 2

(GS) ~ X X Gαα = rj→k,αα, (12a) 2πA2 (GS) j∈D− k∈D+ FIG. 4. Ground state metric tensor G12 and Berry curva- (GS) (GS) X X (0) (π) ture F12 for N = 3 dots at ϕ3 = 0.4 π in the topological ~ (GS) Gαβ = (rj→k,αβ − rj→k,αβ), (12b) 8πA2 phase with C12 = −1 [cf. Fig. 2(c)]. The elements of the j∈D− k∈D+ QGT are strongly peaked for values of ϕ close to a Weyl point. (GS) ~ X X (π/2) (−π/2) The coordinates of the Weyl point indicated in this plot are F = r − r , (12c) ϕ = (ϕ , ϕ , ϕ ) ≈ (1.603, 0.781, 0.384)π [cf. lower αβ 4πA2 j→k,αβ j→k,αβ W W,1 W,2 W,3 j∈D− k∈D+ red dot in Fig. 2(b)]. The visible range of values of the elements of the QGT has been limited from −5 to +5 to make the overall structure for small values visible. In fact, the val- which shows that, indeed, only transitions from occu- ues of the elements of the QGT at the point (ϕW,1, ϕW,2) pied (D−) to empty states (D+) are needed. Both the are exceeding ±800. Parameters: ΓNL = 0.8 ΓL, w = 0.2 ΓL, GS Berry curvature and the GS metric tensor can there- ε1 = −ε2 = ε3 = −ΓL/2, ϕ0 = 0. fore be simply obtained from all the experimentally accessible transitions by using microwave spectroscopy. In Fig.4, we show, as an example of Eq. (12), the metric IV. DISCUSSION AND CONCLUSION tensor and the Berry curvature of the GS for N = 3 dots in the topological region for ϕ3 = 0.4 π in which (GS) In this work we have proposed an experimentally fea- we find a Chern number C12 = −1 [cf. Fig. 2(c)]. A local measurement of the elements of the QGT by mi- sible system consisting of a linear chain of N quantum crowave spectroscopy will not only reveal the topological dots connected to N +1 s-wave superconducting leads, as phase in terms of the Chern number by integrating the shown in Fig.1, to study quantum geometry in topologi- Berry curvature, it is also possible to locate Weyl points cal Josephson matter. This linear setup does not require in phase space TN for which zero-energy states are ex- a direct connection between the different superconduct- pected. In fact, Fig.4 shows that there are regions where ing terminals, as opposed to the system considered in the elements of the QGT are significantly different from Ref. [30], which might take a lot of engineering effort to zero, signaling the location of a Weyl point located at realize in practice. Furthermore, such quantum-dot mod- els allow for the deterministic study of the local quantum ϕW/π ≈ (1.603, 0.781, 0.384), with the ϕ1 and ϕ2 coordinates indicated in the plots. A hint to another Weyl geometry in topological Josephson matter in contrast to point is visible as a blurred point at (ϕ1, ϕ2) ≈ (0.9, 1.5)π prior proposals based on statistical methods [22, 24, 26]. which is actually in a different ϕ3 = const. plane [cf. We have derived an effective low-energy Hamiltonian Fig. 2(b)] and all of the Weyl points can be located by of the chain, presented in Eq. (1), which reveals large many measurements for different fixed ϕ3. While possi- topological regions in terms of the integer-valued Chern ble zero-energy states might be useful on its own, let us number and, therefore, leaves large experimental freedom recall here that the Chern number will be responsible for for the concrete values of the various system parameters a quantized transconductance across two terminals [22]. to study quantum geometry and topology in the syn- 7 thetic space of N independent superconducting phases was supported by the Deutsche Forschungsgemeinschaft [cf. Fig.2]. In particular, as discussed in Fig. 2(c), it through SFB 767 and Grant No. RA 2810/1. J.C.C. is also not necessary for the dots to be directly coupled acknowledges the support via the Mercator Program of (w = 0) and it turns out that a larger value of the di- the Deutsche Forschungsgemeinschaft in the frame of the rect coupling can even lead to an enhancement of the SFB 767. robustness of the topological phase (w = 2 ΓL). We furthermore believe that this or similar engineered systems also provide useful and easily scalable condensed-matter Appendix A: Effective Hamiltonian for N dots realizations to study higher-dimensional topological invariants such as, e.g., higher-dimensional Chern numbers Here, we derive the effective Hamiltonian of N dots [88–99] or the so-called Dixmier-Douady invariants [100– connected to N+1 superconducting terminals as sketched 104], since the synthetic dimension N of the space of in Fig.1. Due to the inclusion of superconductivity, we N superconducting phases T scales with the number of start by defining a set of Pauli matrices (τ1, τ2, τ3) in superconducting leads which can be, in principle, arbi- Nambu space. The starting Hamiltonian reads trary. Moreover, we have worked out the insightful connec- H = HD + HS + HDS. (A1) tion between the QGT of single-particle ABS and the The N dots and their couplings are described by HD = † QGT of the corresponding noninteracting GS. At this d HˆDd , where the 2N ×2N matrix Hamiltonian is given point we want to stress that this connection is relevant by to any physical system hosting multiple energy bands.   As only the geometrical property of the GS is of particu- ε1 w12 0 ··· 0 0 lar importance for low-temperature experiments, we have w12 ε2 w23 ··· 0 0   0 w ε ··· 0 0  shown that, although the QGT of an individual ABS can- ˆ  23 3  HD =  ......  ⊗ τ3 not be accessed without further preparation, the QGT of  ......  the GS as a whole as presented in Eq. (12) is still ac-    0 0 0 ··· εN−1 wN−1,N  cessible by microwave spectroscopy from the oscillator 0 0 0 ··· wN−1,N εN strengths of multiple absorption lines. In this regard, we (A2) have generalized the previously discussed application of microwave spectroscopy to measure the QGT of topolog- with εj being the onsite energy of dot j, wj,j+1 being ical Josephson matter in Ref. [30] to the case of the pres- the hopping energy between dots j and j + 1, and d† = ence of multiple ABS. We emphasize that this method † † † (d1↑, d1↓, d2↑, d2↓, . . . , dN↑, dN↓) is the Nambu spinor con- provides an alternative way to measure the GS topology sisting of electronic annihilation (creation) operators d(†) beyond the already mentioned transconductance mea- jσ of spin σ =↑, ↓ on dot j. surements originally proposed in Ref. [22]. Furthermore, The Hamiltonian describing the N + 1 BCS super- microwave spectroscopy provides a tool to gain access to P † ˆ † the local QGT and makes it possible to localize the po- conducting leads reads HS = k ckHS,kck, with ck = † † † sition of zero-energy states through the measurement of (c0k↑, c0(−k)↓, c1k↑, c1(−k)↓, . . . , cNk↑, cN(−k)↓) containing Weyl points in phase space [cf. Fig.4]. (†) electronic annihilation (creation) operators cjkσ of quasi- Since it is based on quantum dots, our system provides momentum k and spin σ in lead j, and the 2(N + 1) × an ideal starting point to study the effect of Coulomb 2(N + 1) block-diagonal matrix Hamiltonian interactions on the geometry and topology of the GS. Futhermore, it may be relevant for concrete experiments HˆS,k = diag(HˆS,0k, HˆS,1k,..., HˆS,Nk). (A3) to study the effect of a finite superconducting gap, which ˆ iϕj τ3 can be of the order of the tunneling [78, 79], and the influ- Here, HS,jk = ξjkτ3 + ∆je τ1, where ∆j and ϕj are ence of the continuum [28, 105] on the quantum geometry the superconducting pairing potential and phase, respec- of the GS, which is not captured in our low-energy theory tively, and ξjk is the normal state dispersion in each lead. at infinite gap. Finally, the general role of a finite tem- Finally, the coupling between the dots and neighboring perature and dissipation on both quantum geometry and superconducting leads is given by topology remains an open question, both theoretically X H = (d†Vˆ c + h.c.), (A4) and experimentally [106–109], and it will be interesting DS DS k to apply these concepts to the proposed system in the k future. with the 2N × 2(N + 1) matrix coupling

t10 t11 0 0 ··· 0 0  0 t21 t22 0 ··· 0 0 ACKNOWLEDGMENTS   ˆ  0 0 t32 t33 ··· 0 0  VDS =   ⊗ τ3 (A5)  ......  We thank Hannes Weisbrich, Hugues Pothier, and  ......  Cristian Urbina for useful discussions. This work 0 0 0 0 ··· tN,N−1 tNN 8 and the hopping energy tjk between dot j and lead k. Appendix B: Quantum geometry of general We use the Green’s function (GF) technique to derive fermionic n-particle states an effective Hamiltonian of the subsystem of the dots. ˆ −1 Fromg ˆS,jk = (ε − HS,jk) , we find the 2 × 2 GF of the For n ∈ N, let |ψji, j ∈ {1, . . . , n}, be single-particle j-th lead as states with their individual single-particle quantum geometric tensors (QGTs) ε + ξ τ + ∆ eiϕj τ3 τ gˆ = jk 3 j 1 . (A6) S,jk ε2 − ξ2 − ∆2 (j) jk j qαβ = h∂αψj|(1 − |ψjihψj|)|∂βψji, (B1) Since the dots and hoppings do not depend on quasi- their Berry curvatures f (j) = −2 Im(q(j)) and their momentum k, we sum over all quasimomenta and define αβ αβ P (j) (j) gˆS,j = gˆS,jk. Turning the summation into an energy (Fubini-Study) metric tensors gαβ = Re(qαβ ). The dif- k P R integration in the wide-band limit via k → N0 dξjk, ferent single-particle states shall be orthonormalized, i.e., with N0 being the normal state density of states at the hψj|ψki = δjk. Furthermore, let Fermi energy, we obtain 1 X h i |Ψ i = √ (−1)p(P) P |ψ i ⊗ · · · ⊗ |ψ i (B2) iϕj τ3 n 1 n ε + ∆ e τ ∆j →∞ j 1 iϕj τ3 n! gˆS,j = −πN0 q −→ −πN0e τ1, (A7) P ∆2 − ε2 j be the corresponding general antisymmetric noninteracting n-particle fermionic state, in which we sum over all where we assume that the pairing potentials are larger n! different permutations with the permutation operator than all relevant energy scales. Finally, the bare 2(N + P and the number of transpositions p(P) of a particular 1) × 2(N + 1) matrix GF of the subsystem of the super- permutation. How are the single-particle QGTs related conducting leads is given by the block-diagonal matrix to the corresponding noninteracting n-particle QGT

gˆS = diag(ˆgS,0,..., gˆS,N ), (A8) (n) Qαβ = h∂αΨn|(1 − |ΨnihΨn|)|∂βΨni (B3) written in the basis c† = (c†,..., c† ) = P c† . 0 N k k and, in particular, what are the corresponding relations Formally, the bare 2N ×2N matrix GF of the dots with (n) (n) −1 to the Berry curvature F = −2 Im(Q ) and the met- their couplings is given byg ˆD = (ε−HˆD) . The explicit αβ αβ (n) (n) form ofg ˆD is not needed for the following derivation. ric tensor Gαβ = Re(Qαβ ) of the n-particle state |Ψni? ˆ (N) After a lengthy but straightforward calculation, the We define an effective matrix Hamiltonian Heff from the dressed GF of the 2N × 2N dot subsystem, i.e., final result reads ˆ ˆ (N) −1 GD = (ε − Heff ) , which is the solution of the Dyson n n ˆ ˆ ˆ (n) X (j) X equation GD =g ˆD +g ˆDΣDGD, with the self-energy Gαβ = gαβ − h∂αψj|ψkihψk|∂βψji , (B4a) ˆ ˆ ˆ † j=1 j,k=1 ΣD = VDS gˆS VDS. Using the formal definition ofg ˆD, we j6=k (N) obtain the solution Hˆ = HˆD +Σˆ D which describes the n eff (n) X (j) dots proximity coupled to the superconductors. Finally, Fαβ = fαβ , (B4b) the effective low-energy Hamiltonian is j=1

† (N) which is also presented in Eq. (6) in Sec.II of the main Heff = d Hˆ d = HD + HLP + HNLP, (A9) eff text for the GS of the system that contains N occupied with the local and nonlocal pairings (LP and NLP) orig- single-particle ABS. While the result for the metric ten- inating from the self-energy Σˆ D given by sor reveals a nontrivial and less intuitive relation, the result for the Berry curvature can be understood more N intuitively and matches the expected result. On an intu- X † (j−1) iϕj−1τ3 (j) iϕj τ3 HLP = dj Γj e + Γj e τ1dj, itive level, suppose that the single-particle states acquire iφj j=1 the Berry phases φj according to |ψji 7→ e |ψji after (A10a) a cyclic evolution in parameter space [83]. Then, the N−1 evolution of the corresponding n-particle state results in X (j) † iϕj τ3 † iϕj τ3 HNLP = Γ (d e τ1d + d e τ1d ), j,j+1 j j+1 j+1 j iΦn j=1 |Ψni 7→ e |Ψni (B5) (A10b) according to the definition in Eq. (B2), where the Berry Pn (k) 2 (k) phase Φn = j=1 φj is the sum of the individual Berry with Γ = πN0t ,Γ = πN0tjktj+1,k and the local j jk j,j+1 phases φj. This directly translates to the Berry curvature † † Nambu spinors dj = (dj↑, dj↓). Eq. (A9), together with by means of Stokes’ integral theorem. Eq. (A10), represents the effective Hamiltonian presented As a direct consequence of Eq. (B4b), the Chern num- in Eq. (1) in Sec.II of the main text. ber of the GS is simply given by the sum of all individual 9

(γ) Chern numbers of the N occupied single-particle states according to Fermi’s golden rule. Note that rn→m,αβ = as expected, i.e., (−γ) rm→n,αβ by deﬁnition. Using the identities −1 (GS) X (j) C = c , (B6) αβ αβ h∂αψn| ψmi = −hψn| ∂αψmi, (C3a) j=−N ˆ (N) hψn| (∂αHeﬀ )|ψmi = (εn − εm)h∂αψn| ψmi, (C3b) which is also the result presented in Eq. (7) in Sec.II of the main text. for n 6= m, we obtain

2 Appendix C: Microwave spectroscopy of the 2πA (nm) rn→m,αα = qαα , (C4a) quantum geometry of the ground state ~ 2 (γ) 2πA n (nm) (nm) (n) (n) rn→m,αβ = qαα + qββ For the diagonal elements, qαα = gαα , we need ~ o to apply a perturbation to one phase ϕα → ϕα + iγ (nm) −iγ (nm) + e qαβ + e qβα . (C4b) 2A cos(ωt)/~ω, resulting the transition rate Rn→m,αα = rn→m,αα δ(εm − εn − ~ω) with the oscillator strength [30, 87] which are exactly Eqs. (9) and (10) as presented in Sec.III of the main text. Finally, for α 6= β, this yields 2 ˆ (N) 2 2πA hψm|(∂αHeff )|ψni the useful identities rn→m,αα = 2 (C1) ~(εm − εn) Re(q(nm)) = ~ (r(0) − r(π) ), (C5a) according to Fermi’s golden rule. In order to measure αβ 8πA2 n→m,αβ n→m,αβ (n) the off-diagonal elements of the QGT qαβ for α 6= β, we (nm) ~ (π/2) (−π/2) Im(qαβ ) = − rn→m,αβ − rn→m,αβ . (C5b) need to simultaneously apply the perturbations ϕα → 8πA2 ϕα + 2A cos(ωt)/~ω and ϕβ → ϕβ + 2A cos(ωt − γ)/~ω At zero temperature, the GS is given by an (antisymmet- with a phase difference γ, resulting in the transition rates ric) combination of all N occupied single-particle states (γ) (γ) Rn→m,αβ = rn→m,αβ δ(εm − εn − ~ω) with the oscillator up to the Fermi energy, as generally defined in Eq. (B2). strength [30, 87] We now use the general result in Eq. (B4) to relate the measured oscillator strengths from microwave spec- 2 ˆ (N) iγ ˆ (N) 2 (GS) (γ) 2πA hψm| ∂αHeff + e ∂βHeff |ψni troscopy to the Berry curvature Fαβ and the metric rn→m,αβ = 2 (GS) ~(εm − εn) tensor Gαβ of the ground state. (C2) For the Berry curvature, we obtain

−1 −1 N −1 N (GS) (B4b) X (j) (8) X X (jk) (C5b) X X (π/2) (−π/2) F = f = −2 Im(q ) = ~ r − r αβ αβ αβ 4πA2 j→k,αβ j→k,αβ j=−N j=−N k=−N j=−N k=−N k6=0,j k6=0,j " −1 −1 N # −1 N X (π/2) (−π/2) X X (π/2) (−π/2) X X (π/2) (−π/2) = ~ r − r + r − r = ~ r − r , (C6) 4πA2 j→k,αβ j→k,αβ j→k,αβ j→k,αβ 4πA2 j→k,αβ j→k,αβ j,k=−N j=−N k=1 j=−N k=1 k6=j | {z } = 0, Eq. (C2) which is the one presented in Eq. (12c) in Sec.III of the main text. For the metric tensor, we obtain

−1 −1 −1 N −1 (GS) (B4a) X (j) X (8) X X (jk) X (jk) Gαβ = gαβ − h∂αψj|ψki hψk|∂βψji = Re(qαβ ) − qαβ j=−N j,k=−N j=−N k=−N j,k=−N j6=k k6=0,j k6=j −1 N −1 −1 N X X (jk) 1 X (jk) (kj) X X (jk) = Re(q ) − (q + q ) = Re(q ) . (C7) αβ 2 αβ αβ αβ j=−N k=−N j,k=−N | {z } j=−N k=1 k6=0,j k6=j (jk) 2 Re(qαβ )

In both cases we see that we only need to measure transitions between occupied and empty states. Finally, we get for 10 the diagonal and oﬀ-diagonal elements

−1 N (C4a) X X G(GS) = ~ r , (C8a) αα 2πA2 j→k,αα j=−N k=1 −1 N (GS) (C5a) X X (0) (π) G = ~ (r − r ) , (C8b) αβ 8πA2 j→k,αβ j→k,αβ j=−N k=1 which are the relations presented in Eqs. (12a) and (12b), respectively, in Sec. III of the main text.

[1] A. Bansil, H. Lin, and T. Das, Colloquium: Topological [16] M. Hays, G. de Lange, K. Serniak, D. J. van Woerkom, band theory, Rev. Mod. Phys. 88, 021004 (2016). D. Bouman, P. Krogstrup, J. Nyg˚ard, A. Geresdi, [2] C. L. Kane and E. J. Mele, Z2 Topological Order and the and M. H. Devoret, Direct Microwave Measurement Quantum Spin Hall Effect, Phys. Rev. Lett. 95, 146802 of Andreev-Bound-State Dynamics in a Semiconductor- (2005). Nanowire Josephson Junction, Phys. Rev. Lett. 121, [3] M. Z. Hasan and C. L. Kane, Colloquium: Topological 047001 (2018). insulators, Rev. Mod. Phys. 82, 3045 (2010). [17] J.-D. Pillet, C. H. L. Quay, P. Morfin, C. Bena, [4] N. P. Armitage, E. J. Mele, and A. Vishwanath, Weyl A. Levy Yeyati, and P. Joyez, Andreev bound states in and Dirac semimetals in three-dimensional solids, Rev. supercurrent-carrying carbon nanotubes revealed, Nat. Mod. Phys. 90, 015001 (2018). Phys. 6, 965 (2010). [5] M. Sato and Y. Ando, Topological superconductors: A [18] L. Bretheau, Ç. O.¨ Girit, C. Urbina, D. Esteve, review, Rep. Prog. Phys. 80, 076501 (2017). and H. Pothier, Supercurrent Spectroscopy of Andreev [6] K. Kawabata, K. Shiozaki, M. Ueda, and M. Sato, Sym- States, Phys. Rev. X 3, 041034 (2013). metry and Topology in Non-Hermitian Physics, Phys. [19] B. van Heck, S. Mi, and A. R. Akhmerov, Single fermion Rev. X 9, 041015 (2019). manipulation via superconducting phase differences in [7] A. Y. Kitaev, Unpaired Majorana fermions in quantum multiterminal Josephson junctions, Phys. Rev. B 90, wires, Phys. Usp. 44, 131 (2001). 155450 (2014). [8] A. Y. Kitaev, Fault-tolerant quantum computation by [20] C. Padurariu, T. Jonckheere, J. Rech, R. Mélin,D. Fein- anyons, Annals of Physics 303, 2–30 (2003). berg, T. Martin, and Y. V. Nazarov, Closing the prox- [9] L. Fu and C. L. Kane, Superconducting Proximity Effect imity gap in a metallic Josephson junction between and Majorana Fermions at the Surface of a Topological three superconductors, Phys. Rev. B 92, 205409 (2015). Insulator, Phys. Rev. Lett. 100, 096407 (2008). [21] T. Yokoyama and Y. V. Nazarov, Singularities in the [10] G. Yang, Z. Lyu, J. Wang, J. Ying, X. Zhang, J. Shen, Andreev spectrum of a multiterminal Josephson junc- G. Liu, J. Fan, Z. Ji, X. Jing, F. Qu, and L. Lu, Pro- tion, Phys. Rev. B 92, 155437 (2015). tected gap closing in Josephson trijunctions constructed [22] R.-P. Riwar, M. Houzet, J. S. Meyer, and Y. V. Nazarov, on Bi2Te3, Phys. Rev. B 100, 180501(R) (2019). Multi-terminal Josephson junctions as topological mat- [11] N. M. Chtchelkatchev and Y. V. Nazarov, Andreev ter, Nat. Commun. 7, 11167 (2016). Quantum Dots for Spin Manipulation, Phys. Rev. Lett. [23] H.-Y. Xie, M. G. Vavilov, and A. Levchenko, Weyl 90, 226806 (2003). nodes in Andreev spectra of multiterminal Josephson [12] A. Zazunov, V. S. Shumeiko, E. N. Bratus, J. Lantz, junctions: Chern numbers, conductances, and super- and G. Wendin, Andreev Level Qubit, Phys. Rev. Lett. currents, Phys. Rev. B 97, 035443 (2018). 90, 087003 (2003). [24] E. Eriksson, R.-P. Riwar, M. Houzet, J. S. Meyer, and [13] L. Bretheau, Ç. O.¨ Girit, H. Pothier, D. Esteve, and Y. V. Nazarov, Topological transconductance quantiza- C. Urbina, Exciting Andreev pairs in a superconducting tion in a four-terminal Josephson junction, Phys. Rev. atomic contact, Nature (London) 499, 312 (2013). B 95, 075417 (2017). [14] C. Janvier, L. Tosi, L. Bretheau, Ç. O.¨ Girit, M. Stern, [25] H.-Y. Xie, M. G. Vavilov, and A. Levchenko, Topolog- P. Bertet, P. Joyez, D. Vion, D. Esteve, M. F. Goff- ical Andreev bands in three-terminal Josephson junc- man, H. Pothier, and C. Urbina, Coherent manipulation tions, Phys. Rev. B 96, 161406 (2017). of Andreev states in superconducting atomic contacts, [26] J. S. Meyer and M. Houzet, Nontrivial Chern Num- Science 349, 1199 (2015). bers in Three-Terminal Josephson Junctions, Phys. Rev. [15] D. J. van Woerkom, A. Proutski, B. van Heck, Lett. 119, 136807 (2017). D. Bouman, J. I. Väyrynen, L. I. Glazman, [27] L. P. Gavensky, G. Usaj, and C. A. Balseiro, Topological P. Krogstrup, J. Nyg˚ard,L. P. Kouwenhoven, and phase diagram of a three-terminal Josephson junction: A. Geresdi, Microwave spectroscopy of spinful Andreev From the conventional to the Majorana regime, Phys. bound states in ballistic semiconductor Josephson junc- Rev. B 100, 014514 (2019). tions, Nat. Phys. 13, 876 (2017). [28] E. V. Repin, Y. Chen, and Y. V. Nazarov, Topological properties of multiterminal superconducting nanostruc- 11

tures: Effect of a continuous spectrum, Phys. Rev. B ric without Berry curvature, Phys. Rev. B 94, 134423 99, 165414 (2019). (2016). [29] T. Yokoyama, J. Reutlinger, W. Belzig, and [48] F. Freimuth, S. Blügel,and Y. Mokrousov, Geometrical Y. V. Nazarov, Order, disorder, and tunable gaps in contributions to the exchange constants: Free electrons the spectrum of Andreev bound states in a multitermi- with spin-orbit interaction, Phys. Rev. B 95, 184428 nal superconducting device, Phys. Rev. B 95, 045411 (2017). (2017). [49] L. P. Gavensky, G. Usaj, D. Feinberg, and C. A. Bal- [30] R. L. Klees, G. Rastelli, J. C. Cuevas, and W. Belzig, seiro, Berry curvature tomography and realization of Microwave Spectroscopy Reveals the Quantum Geomet- topological Haldane model in driven three-terminal ric Tensor of Topological Josephson Matter, Phys. Rev. Josephson junctions, Phys. Rev. B 97, 220505 (2018). Lett. 124, 197002 (2020). [50] B. Venitucci, D. Feinberg, R. Mélin,and B. Dou¸cot, [31] Y. Chen and Y. V. Nazarov, Spin Weyl quantum unit: A Nonadiabatic Josephson current pumping by chiral mi- theoretical proposal, Phys. Rev. B 103, 045410 (2021). crowave irradiation, Phys. Rev. B 97, 195423 (2018). [32] M. Houzet and J. S. Meyer, Majorana-Weyl crossings in [51] N. Fläschner, B. S. Rem, M. Tarnowski, D. Vogel, D.- topological multiterminal junctions, Phys. Rev. B 100, S. Lühmann,K. Sengstock, and C. Weitenberg, Exper- 014521 (2019). imental reconstruction of the Berry curvature in a Flo- [33] M. Trif and P. Simon, Dynamics of a Majorana Trijunc- quet Bloch band, Science 352, 1091 (2016). tion in a Microwave Cavity, Adv. Quantum Technol. 2, [52] L. Asteria, D. T. Tran, T. Ozawa, M. Tarnowski, 1900091 (2019). B. S. Rem, N. Fläschner, K. Sengstock, N. Gold- [34] K. Sakurai, M. T. Mercaldo, S. Kobayashi, A. Yamak- man, and C. Weitenberg, Measuring quantized circular age, S. Ikegaya, T. Habe, P. Kotetes, M. Cuoco, and dichroism in ultracold topological matter, Nat. Phys. Y. Asano, Nodal Andreev spectra in multi-Majorana 15, 449 (2019). three-terminal Josephson junctions, Phys. Rev. B 101, [53] X. Tan, D.-W. Zhang, Q. Liu, G. Xue, H.-F. Yu, Y.- 174506 (2020). Q. Zhu, H. Yan, S.-L. Zhu, and Y. Yu, Topological [35] E. V. Repin and Y. V. Nazarov, Weyl points in the Maxwell Metal Bands in a Superconducting Qutrit, multi-terminal Hybrid Superconductor-Semiconductor Phys. Rev. Lett. 120, 130503 (2018). Nanowire devices, arXiv:2010.11494. [54] M. D. Schroer, M. H. Kolodrubetz, W. F. Kindel, [36] M. Kolodrubetz, V. Gritsev, and A. Polkovnikov, Clas- M. Sandberg, J. Gao, M. R. Vissers, D. P. Pappas, sifying and measuring geometry of a quantum ground A. Polkovnikov, and K. W. Lehnert, Measuring a Topo- state manifold, Phys. Rev. B 88, 064304 (2013). logical Transition in an Artificial Spin-1/2 System, [37] T. Neupert, C. Chamon, and C. Mudry, Measuring the Phys. Rev. Lett. 113, 050402 (2014). quantum geometry of Bloch bands with current noise, [55] X. Tan, D.-W. Zhang, Z. Yang, J. Chu, Y.-Q. Zhu, Phys. Rev. B 87, 245103 (2013). D. Li, X. Yang, S. Song, Z. Han, Z. Li, Y. Dong, H.- [38] P. Zanardi, P. Giorda, and M. Cozzini, Information- F. Yu, H. Yan, S.-L. Zhu, and Y. Yu, Experimental Mea- Theoretic Differential Geometry of Quantum Phase surement of the Quantum Metric Tensor and Related Transitions, Phys. Rev. Lett. 99, 100603 (2007). Topological Phase Transition with a Superconducting [39] O. Bleu, G. Malpuech, Y. Gao, and D. D. Solnyshkov, Qubit, Phys. Rev. Lett. 122, 210401 (2019). Effective Theory of Nonadiabatic Quantum Evolution [56] P. Roushan, C. Neill, Y. Chen, M. Kolodrubetz, Based on the Quantum Geometric Tensor, Phys. Rev. C. Quintana, N. Leung, M. Fang, R. Barends, B. Camp- Lett. 121, 020401 (2018). bell, Z. Chen, B. Chiaro, A. Dunsworth, E. Jeffrey, [40] D. D. Solnyshkov, C. Leblanc, L. Bessonart, A. Nalitov, J. Kelly, A. Megrant, J. Mutus, P. J. J. O’Malley, J. Ren, Q. Liao, F. Li, and G. Malpuech, Quantum D. Sank, A. Vainsencher, J. Wenner, T. White, metric and wave packets at exceptional points in non- A. Polkovnikov, A. N. Cleland, and J. M. Martinis, Ob- Hermitian systems, Phys. Rev. B 103, 125302 (2021). servation of topological transitions in interacting quan- [41] S. Peotta and P. Törmä,Superfluidity in topologically tum circuits, Nature 515, 241 (2014). nontrivial flat bands, Nat. Commun. 6, 8944 (2015). [57] M. Yu, P. Yang, M. Gong, Q. Cao, Q. Lu, H. Liu, [42] A. Julku, S. Peotta, T. I. Vanhala, D.-H. Kim, and S. Zhang, M. B. Plenio, F. Jelezko, T. Ozawa, N. Gold- P. Törmä,Geometric Origin of Superfluidity in the Lieb- man, S.-L. Zhang, and J.-M. Cai, Experimental mea- Lattice Flat Band, Phys. Rev. Lett. 117, 045303 (2016). surement of the quantum geometric tensor using cou- [43] L. Liang, T. I. Vanhala, S. Peotta, T. Siro, A. Harju, pled qubits in diamond, Nat. Sci. Rev. 7, 254 (2019). and P. Törmä,Band geometry, Berry curvature, and [58] A. Gianfrate, O. Bleu, L. Dominici, V. Ardizzone, superfluid weight, Phys. Rev. B 95, 024515 (2017). M. De Giorgi, D. Ballarini, G. Lerario, K. W. West, [44] A. Srivastava and A. Imamoˇglu,Signatures of Bloch- L. N. Pfeiffer, D. D. Solnyshkov, D. Sanvitto, and Band Geometry on Excitons: Nonhydrogenic Spectra G. Malpuech, Measurement of the quantum geometric in Transition-Metal Dichalcogenides, Phys. Rev. Lett. tensor and of the anomalous Hall drift, Nature 578, 381 115, 166802 (2015). (2020). [45] Y. Gao, S. A. Yang, and Q. Niu, Field Induced Posi- [59] M. Wimmer, H. M. Price, I. Carusotto, and U. Peschel, tional Shift of Bloch Electrons and Its Dynamical Im- Experimental measurement of the Berry curvature from plications, Phys. Rev. Lett. 112, 166601 (2014). anomalous transport, Nat. Phys. 13, 545 (2017). [46] Y. Gao, S. A. Yang, and Q. Niu, Geometrical effects in [60] T. Ozawa, H. M. Price, A. Amo, N. Goldman, orbital magnetic susceptibility, Phys. Rev. B 91, 214405 M. Hafezi, L. Lu, M. C. Rechtsman, D. Schuster, J. Si- (2015). mon, O. Zilberberg, and I. Carusotto, Topological pho- [47] F. Piéchon, A. Raoux, J.-N. Fuchs, and G. Montam- tonics, Rev. Mod. Phys. 91, 015006 (2019). baux, Geometric orbital susceptibility: Quantum met- [61] P. Kotetes, M. T. Mercaldo, and M. Cuoco, Synthetic 12

Weyl Points and Chiral Anomaly in Majorana Devices Rev. Lett. 91, 057005 (2003). with Nonstandard Andreev-Bound-State Spectra, Phys. [78] L. Hofstetter, S. Csonka, J. Nyg˚ard, and Rev. Lett. 123, 126802 (2019). C. Schönenberger, Cooper pair splitter realized in [62] L. J. Geerligs, S. M. Verbrugh, P. Hadley, J. E. Mooij, a two-quantum-dot Y-junction, Nature 461, 960 H. Pothier, P. Lafarge, C. Urbina, D. Estève, and (2009). M. H. Devoret, Single Cooper pair pump, Z. Phys. B: [79] L. G. Herrmann, F. Portier, P. Roche, A. Levy Yey- Condens. Matter 85, 349 (1991). ati, T. Kontos, and C. Strunk, Carbon Nanotubes [63] J. P. Pekola, J. J. Toppari, M. Aunola, as Cooper-Pair Beam Splitters, Phys. Rev. Lett. 104, M. T. Savolainen, and D. V. Averin, Adiabatic 026801 (2010). transport of Cooper pairs in arrays of Josephson [80] G. Fülöp,F. Dom´ınguez,S. d’Hollosy, A. Baumgartner, junctions, Phys. Rev. B 60, R9931(R) (1999). P. Makk, M. H. Madsen, V. A. Guzenko, J. Nyg˚ard, [64] M. Aunola and J. J. Toppari, Connecting Berry’s phase C. Schönenberger, A. Levy Yeyati, and S. Csonka, and the pumped charge in a Cooper pair pump, Phys. Magnetic Field Tuning and Quantum Interference in Rev. B 68, 020502(R) (2003). a Cooper Pair Splitter, Phys. Rev. Lett. 115, 227003 [65] R. Leone and L. Lévy, Topological quantization by con- (2015). trolled paths: Application to Cooper pairs pumps, Phys. [81] P. Recher, E. V. Sukhorukov, and D. Loss, Andreev Rev. B 77, 064524 (2008). tunneling, Coulomb blockade, and resonant transport [66] P. A. Erdman, F. Taddei, J. T. Peltonen, R. Fazio, and of non-local spin-entangled electrons, Phys. Rev. B 63, J. P. Pekola, Fast and accurate Cooper pair pump, Phys. 165314 (2001). Rev. B 100, 235428 (2019). [82] D. Feinberg, Andreev scattering and cotunneling be- [67] V. Fatemi, A. R. Akhmerov, and L. Bretheau, Weyl tween two superconductor-normal metal interfaces: the Josephson circuits, Phys. Rev. Research 3, 013288 dirty limit, Eur. Phys. J. B 36, 419 (2003). (2021). [83] M. V. Berry, Quantal phase factors accompanying adi- [68] L. Peyruchat, J. Griesmar, J.-D. Pillet, and Ç. O.¨ Girit, abatic changes, Proc. R. Soc. A 392, 45 (1984). Transconductance quantization in a topological Joseph- [84] M. Nakahara, Geometry, Topology, and Physics (Insti- son tunnel junction circuit, Phys. Rev. Research 3, tute of Physics Publishing, Bristol, 2003). 013289 (2021). [85] T. Fukui, Y. Hatsugai, and H. Suzuki, Chern Numbers [69] T. Herrig and R.-P. Riwar, A “minimal” topological in Discretized Brillouin Zone: Efficient Method of Com- quantum circuit, arXiv:2012.10655. puting (Spin) Hall Conductances, J. Phys. Soc. Jpn. 74, [70] E. Strambini, S. D’Ambrosio, F. Vischi, F. S. Bergeret, 1674 (2005). Y. V. Nazarov, and F. Giazotto, The ω-SQUIPT as a [86] M. Kolodrubetz, D. Sels, P. Mehta, and A. Polkovnikov, tool to phase-engineer Josephson topological materials, Geometry and non-adiabatic response in quantum and Nat. Nanotechnol. 11, 1055 (2016). classical systems, Phys. Rep. 697, 1 (2017). [71] M. Amundsen, J. A. Ouassou, and J. Linder, Ana- [87] T. Ozawa and N. Goldman, Extracting the quantum lytically determined topological phase diagram of the metric tensor through periodic driving, Phys. Rev. B proximity-induced gap in diffusive n-terminal Joseph- 97, 201117 (2018). son junctions, Sci. Rep. 7, 40578 (2017). [88] S.-C. Zhang and J. Hu, A Four-Dimensional General- [72] S. R. Plissard, I. van Weperen, D. Car, M. A. Verheijen, ization of the Quantum Hall Effect, Science 294, 823 G. W. G. Immink, J. Kammhuber, L. J. Cornelissen, (2001). D. B. Szombati, A. Geresdi, S. M. Frolov, L. P. Kouwen- [89] Y. E. Kraus, Z. Ringel, and O. Zilberberg, Four- hoven, and E. P. A. M. Bakkers, Formation and elec- Dimensional Quantum Hall Effect in a Two- tronic properties of InSb nanocrosses, Nat. Nanotech- Dimensional Quasicrystal, Phys. Rev. Lett. 111, nol. 8, 859 (2013). 226401 (2013). [73] A. W. Draelos, M.-T. Wei, A. Seredinski, H. Li, [90] H. M. Price, O. Zilberberg, T. Ozawa, I. Carusotto, Y. Mehta, K. Watanabe, T. Taniguchi, I. V. Borzenets, and N. Goldman, Four-Dimensional Quantum Hall Ef- F. Amet, and G. Finkelstein, Supercurrent flow in multi- fect with Ultracold Atoms, Phys. Rev. Lett. 115, 195303 terminal graphene Josephson junctions, Nano Lett. 19, (2015). 1039 (2019). [91] T. Ozawa, H. M. Price, N. Goldman, O. Zilberberg, [74] N. Pankratova, H. Lee, R. Kuzmin, K. Wickramasinghe, and I. Carusotto, Synthetic dimensions in integrated W. Mayer, J. Yuan, M. G. Vavilov, J. Shabani, and photonics: From optical isolation to four-dimensional V. E. Manucharyan, Multiterminal Josephson Effect, quantum Hall physics, Phys. Rev. A 93, 043827 (2016). Phys. Rev. X 10, 031051 (2020). [92] C. Chan and X.-J. Liu, Non-Abelian Majorana Modes [75] E. G. Arnault, T. Larson, A. Seredinski, L. Zhao, H. Li, Protected by an Emergent Second Chern Number, Phys. K. Watanabe, T. Taniguchi, I. V. Borzenets, F. Amet, Rev. Lett. 118, 207002 (2017). and G. Finkelstein, The Multi-terminal Inverse AC [93] S. Sugawa, F. Salces-Carcoba, A. R. Perry, Y. Yue, and Josephson Effect, arXiv:2012.15253. I. B. Spielman, Second Chern number of a quantum- [76] A. Levy Yeyati, J. C. Cuevas, A. López-Dávalos, and simulated non-Abelian Yang monopole, Science 360, A. Mart´ın-Rodero, Resonant tunneling through a small 1429 (2018). quantum dot coupled to superconducting leads, Phys. [94] M. Lohse, C. Schweizer, H. M. Price, O. Zilberberg, and Rev. B 55, R6137(R) (1997). I. Bloch, Exploring 4D quantum Hall physics with a 2D [77] M. R. Buitelaar, W. Belzig, T. Nussbaumer, B. Babić, topological charge pump, Nature 553, 55 (2018). C. Bruder, and C. Schönenberger, Multiple Andreev Re- [95] O. Zilberberg, S. Huang, J. Guglielmon, M. Wang, flections in a Carbon Nanotube Quantum Dot, Phys. K. P. Chen, Y. E. Kraus, and M. C. Rechtsman, Pho- tonic topological boundary pumping as a probe of 4D 13

quantum Hall physics, Nature 553, 59 (2018). Monopoles with a Superconducting Qudit, Phys. Rev. [96] I. Petrides, H. M. Price, and O. Zilberberg, Six- Lett. 126, 017702 (2021). dimensional quantum Hall effect and three-dimensional [104] M. Chen, C. Li, G. Palumbo, Y.-Q. Zhu, N. Gold- topological pumps, Phys. Rev. B 98, 125431 (2018). man, and P. Cappellaro, Experimental characterization [97] C. H. Lee, Y. Wang, Y. Chen, and X. Zhang, Electro- of the 4D tensor monopole and topological nodal rings, magnetic response of quantum Hall systems in dimen- arXiv:2008.00596. sions five and six and beyond, Phys. Rev. B 98, 094434 [105] M. Zgirski, L. Bretheau, Q. Le Masne, H. Pothier, (2018). D. Esteve, and C. Urbina, Evidence for Long-Lived [98] L. Lu, H. Gao, and Z. Wang, Topological one-way fiber Quasiparticles Trapped in Superconducting Point Con- of second Chern number, Nat. Commun. 9, 5384 (2018). tacts, Phys. Rev. Lett. 106, 257003 (2011). [99] H. Weisbrich, R. L. Klees, G. Rastelli, and W. Belzig, [106] J. Dittmann and G. Rudolph, On a connection govern- Second Chern Number and Non-Abelian Berry Phase in ing parallel transport along 2 × 2 density matrices,J. Topological Superconducting Systems, PRX Quantum Geom. Phys. 10, 93 (1992). 2, 010310 (2021) [107] O. Viyuela, A. Rivas, and M. A. Martin-Delgado, [100] M. K. Murray, Bundle Gerbes, J. London Math. Soc. Uhlmann Phase as a Topological Measure for One- 54, 405 (1996). Dimensional Fermion Systems, Phys. Rev. Lett. 112, [101] G. Palumbo and N. Goldman, Revealing Tensor 130401 (2014). Monopoles through Quantum-Metric Measurements, [108] O. Viyuela, A. Rivas, S. Gasparinetti, A. Wallraff, Phys. Rev. Lett. 121, 170401 (2018). S. Filipp, and M. A. Martin-Delgado, Observation [102] G. Palumbo and N. Goldman, Tensor Berry connec- of topological Uhlmann phases with superconducting tions and their topological invariants, Phys. Rev. B 99, qubits, npj Quantum Inf. 4, 10 (2018). 045154 (2019). [109] C. Gneiting, A. Koottandavida, A. V. Rozhkov, and [103] X. Tan, D.-W. Zhang, W. Zheng, X. Yang, S. Song, F. Nori, Unraveling the topology of dissipative quantum Z. Han, Y. Dong, Z. Wang, D. Lan, H. Yan, S.- systems, arXiv:2007.05960. L. Zhu, and Y. Yu, Experimental Observation of Tensor