Preprint
Article

Majorana Excitons in a Kitaev Chain of Semiconductor Quantum Dots in a Nanowire

Altmetrics

Downloads

82

Views

26

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

04 July 2023

Posted:

06 July 2023

You are already at the latest version

Alerts
Abstract
We present here a theory of Majorana excitons, photo-excited conduction electron-valence band hole pairs, interacting with Majorana Fermions in a Kitaev chain of semiconductor quantum dots embedded in a nanowire. Using analytical tools and exact diagonalisation methods we identify the presence of Majorana Zero Modes in the nanowire absorption spectra.
Keywords: 
Subject: Physical Sciences  -   Condensed Matter Physics

1. Introduction

There is currently interest in realizing synthetic topological quantum matter with topologically protected quasiparticles at its edges [1,2,3], with potential application in topological quantum computation [4,5,6,7,8,9]. Haldane fractional spin quasiparticles in a spin one chain and Majorana Fermions in topological superconductors are good examples [9,10,11]. To realize Majorana Fermions Kitaev proposed [11,12] a chain of quantum dots on a p-wave superconductor that carries such non-local zero energy Majorana Fermions localized on its two ends, the Majorana zero modes (MZMs). Since then there have been numerous proposals to realize the Kitaev chain [13,14,15,16,17,18,19]. In all cases, experimental confirmation of the presence of the MZMs has proved to be a non-trivial and challenging task [20,21,22,23,24,25,26,27,28].
Recent progress in semiconductor quantum dots in nanowires [29,30,31,32,33,34,35,36] opens the possibility of realizing Kitaev chains and optical detection of their Majorana zero modes. In this work, we consider such an array of InAsP quantum dots embedded in an InP nanowire as the material system [29,30,31,32,33,34,35,36] for realization of MZM, and study its signature in light-matter interaction. As the schematic in Figure 1(a) shows, we combine a semiconductor nanowire with p-wave superconductor [37,38,39,40,41,42]. The p-wave pairing in this system is introduced by proximity effect among electrons that are spin-polarized by an external magnetic field, making sure that Cooper pairs can only form between electrons in the conduction band (CB) of adjacent dots. We will show that one can tune the system parameters into topological regime, where two MZMs appear at the two ends of the chain. With semiconductor quantum dots, light can generate a hole in the valence band (VB) and an electron in the conduction band. The electron adds to an existing gas of Majorana Fermions while the hole then interacts with all the quasiparticles of the Kitaev chain, including MZMs, to form composite objects similar to excitons and trions in the Fermi Edge Singularity problem [43,44,45,46]. This leads to a structure in the absorption spectrum of the chain as a function of photon energy. Here we present a theory for the signatures of the MZMs in the optical spectra of the semiconductor nanowire.
After describing the model in Section 2, in Section 3 we introduce the exact diagonalization (ED) method, and introduce Majorana and bond Fermion representation of the Kitaev Hamiltonian. Next, in Section 4 we describe exciton-Majorana Fermion complexes and predict the absorption spectrum. We focus discussion on the optical signature of the MZM in the absorption spectrum. Finally, in Section 5 we conclude by summarising our results and discuss potential experiments detecting Majorana Fermions in a semiconductor Kitaev chain.

2. Kitaev chain in a semiconductor nanowire

Figure 1(a) shows a schematic representation of the Kitaev chain we are considering. It consists of a hexagonal InP nanowire with an array of embedded InAsP quantum dots in the proximity of a p-wave superconductor [37,38,39,40,41,42], in the presence of applied external magnetic field. Such arrays have been extensively investigated [29,30,31,32,33,34,35,36], including their excitonic complexes [29,36]. As Figure 1(b) shows, in our model we include the lowest conduction spin level of each dot and the highest spin valence band level, which are effectively both spin-polarized due to the external magnetic field. Consequently, in the presence of superconductivity, only the electrons from the adjacent conduction levels can pair up, as there is only one conduction level available in each dot. The Kitaev Hamiltonian H e in Equation (1a) describes the hopping and pairing of electrons in conduction band levels. The chemical potential is tuned to bring the chain to near half filling in the absence of superconductivity. Therefore, in the equilibrium the valence levels are full and the system is described by the Kitaev Hamiltonian. However, when a photon with energy close to the band gap of InAsP illuminates the dots, it generates a hole in VB and an electron in CB. The electron becomes one of the electrons in CB and decomposes into quasiparticles of the superconducting state. The hole then forms a bound state with the quasiparticles of the electronic system that is in a collective superconducting state. These possible bound states generate peaks in the absorption spectrum of the system, among which there is the signature of MZM, as we shall show below.
The hole is described by a simple tight binding Hamiltonian H h in Equation (1b). We also consider electron-hole interaction, H i n t in Equation (1c), which is strongest when both conduction electron and valence band hole are on the same quantum dot. Hence, we write the full Hamiltonian of the system as
H = H e + H h + H i n t ,
H e = t i = 1 N 1 c i + 1 c i + h . c . + Δ i = 1 N 1 c i + 1 c i + h . c . μ i = 1 N c i c i ,
H h = τ i = 1 N 1 h i + 1 h i + h . c . + η i = 1 N h i h i ,
H i n t = V i = 1 N n i e n i h ,
where c i ( h i ) is the normal Fermionic creation operator of an electron(hole) in dot i, t ( τ ) is hopping between adjacent conduction(valence) levels, Δ is pairing energy between adjacent conduction levels, μ is the chemical potential measured from the conduction energy level, and η is the CB to VB energy gap in each dot. In the interaction term Equation (1c), V is the Coulomb attraction energy between electrons and holes, where we also introduced n i e = c i c i ( n i h = h i h i ) , the electron(hole) number operator in dot i. Figure 1(b) schematically shows different terms of Equation (1) between two adjacent dots.
Next, before describing the absorption experiment, we start with a brief discussion of Kitaev Hamiltonian.

3. Majorana and bond Fermions in Kitaev Hamiltonian

The Kitaev Hamiltonian H e in Equation (1a), originally introduced in Ref. [11], supports two MZMs localized on the two ends of the chain, when the Hamiltonian is in topological regime. For a finite chain, the topological region is centered on parameters Δ = t and μ = 0 , which is our focus throughout this work. Here, after describing the exact diagonalization (ED) method for normal Fermions, following Kitaev [11], we show how using Majorana Fermions reveals the usefulness of a new set of Fermions we refer to as bond Fermions. Next, after matching energy spectra obtained by ED in both normal and bond Fermion bases, we shall use the bond Fermion basis for the rest of the paper.

3.1. Exact diagonalization in normal Fermion basis

We start off by introducing the exact diagonalization method (ED) for finding the energy spectrum of the Kitaev Hamiltonian. In ED we span the Hilbert space of the system by configuration basis [47]. For our electronic system being made of N spinless orbitals, there are N 0 + N 1 + . . . + N N 1 + N N = 2 N possible configurations, which we construct as
α 1 α N = i = 1 N ( c i ) α i 0 ,
where 0 is the vacuum of electrons, α i = 1 or 0 , corresponds to having(1) or not having(0) electron in orbital i.
For a given number of electrons M we generate electron configurations p M . But as Kitaev Hamiltonian, being a Hamiltonian for a superconductor, does not conserve particle number, its eigenstates are coherent linear combinations of electronic configurations with different electron numbers as
ψ ν = M , p M C M , p M ν M , p M ,
where we are populating N sites with M = 0 , 1 , , N electrons. To solve for coefficients C M , p M ν , we apply the Hamiltonian on this state, and by using the orthogonality of the configurations we obtain the eigenvalue equation
p M , M q M , M H e p M , M C M , p M ν = E ν C M , q M ν .
However, since the Kitaev Hamiltonian H e in Equation (1a) only changes particle number in pairs, the matrix element q M , M H e p M , M is non-zero only if M and M have the same parity, i.e. if they are both even or odd. This parity symmetry allows us to break the Hilbert space into two decoupled subspaces of even and odd configurations. In Appendix A we explicitly show the configurations and the Hamiltonian matrix q M , M H e p M , M in each of these subspaces, for the case of N = 3 .

3.2. Bond Fermions

We now express the Kitaev Hamiltonian in Equation (1a) in terms of Majorana and bond Fermions. First, as schematically shown in Figure 2, we write each electron operator , c and c + in terms of two Majorana Fermion operators γ 1 and γ 2 as
c j = 1 2 ( γ j , 1 + i γ j , 2 ) , c j = 1 2 ( γ j , 1 i γ j , 2 ) ,
where the γ ’s are Majorana Fermion operators. Majorana Fermions satisfy a slightly different anti-commutation relation than the ordinary Fermions, { γ i , α , γ j , β } = 2 δ i j δ α β .
Using Equation (5) and Majorana anti-commutation relations, the Hamiltonian H e can be written in terms of Majorana Fermions as
H e = i 2 ( t + Δ ) j = 1 N 1 γ j , 1 γ j + 1 , 2 + ( t Δ ) j = 1 N 1 γ j + 1 , 1 γ j , 2 μ j = 1 N ( γ j , 1 γ j , 2 i ) .
The form in Equation (6) shows the pairing between Majoranas of different types in adjacent sites. But it also shows that the Hamiltonian is not diagonal in Majorana Fermions, they are not quasiparticles of the Kitaev Hamiltonian. Following Kitaev [11], as shown in Figure 2, we define a new set of Fermionic operators, bond Fermions, which are made of two Majoranas of different types from adjacent sites as
a j = 1 2 ( γ j , 1 + i γ j + 1 , 2 ) = 1 2 ( c j + c j + c j + 1 c j + 1 ) ,
a N = 1 2 ( γ N , 1 + i γ 1 , 2 ) = 1 2 ( c N + c N + c 1 c 1 ) ,
where we also defined a N , to which we refer as the zero mode, out of the two unpaired Majoranas at the two ends of the chain, as shown in Figure 2. Then, the Hamiltonian in terms of bond Fermion operators is
H e = 1 2 ( t + Δ ) j = 1 N 1 ( 2 a j a j 1 ) + ( t Δ ) j = 1 N 1 ( a j + 1 a j 1 + a j + 1 a j 1 + h . c . ) μ j = 1 N ( 1 + ( a j a j 1 + a j a j 1 + h . c . ) ) ,
where in the second and the third sum one should identify a 0 a N . Note that in topological regime, when t = Δ and μ = 0 , the bond Fermions diagonalize the Hamiltonian in Equation (8) and reduce it to
H e = t j = 1 N 1 ( 2 a j a j 1 ) ,
which implies a set of N 1 quasiparticles with energy 2 t , and one non-local quasiparticle a N with zero energy, hence the name zero mode. In this case, since bond Fermions are the quasiparticles of Kitaev Hamiltonian, their configurations are the eigenstates of the system.
In this spirit, we also use bond Fermion configurations for exact diagonalization of Kitaev Hamiltonian. In the same fashion as in Equation (2) we define bond Fermion configurations as
α 1 α N ¯ = i = 1 N ( a i ) α i 0 a ,
where 0 a is the vacuum of bond Fermions, and we used the overline to distinguish these configurations from the normal Fermion configurations. Next, an equation similar to the equation in Equation (3) can be written for the eigenstates of the Hamiltonian in terms of bond Fermion configurations, where now M , p M would represent the p M configuration of having M bond Fermions. And similar to the case of normal Fermions, since H e in Equation (8) conserves the parity of bond Fermion numbers too, we can split the Hilbert space into even and odd subspaces. In Appendix A we explicitly show the bond Fermion configurations and the Hamiltonian matrix of Equation (8) in each of these subspaces, for the case of N = 3 .

3.3. Energy spectrum

To demonstrate the usefulness of bond Fermion basis, we now describe the energy spectrum of a chain of N = 3 quantum dots, obtained both in the normal and bond Fermion basis.
Figure 3 shows the energy spectrum for the case of Δ = t < 0 . Throughout the work, we consider t < 0 , as it is the case for conduction bands hopping integrals. As we mentioned above, in this case, the configurations of bond Fermions are also the eigenstates of the system. Being in topological regime, with these parameters the system has a doubly degenerate ground state, one in the odd subspace GS = 111 ¯ with all bond Fermions, and the other in the even subspace GS ¯ = 110 ¯ , which is missing the zero energy bond Fermion a N a 3 . Next, we have the singly excited states, missing one non-zero bond Fermion, with excitation energy 2 | t | , from which we have two in each subspace, a 1 = 011 ¯ and a 2 = 101 ¯ in the even subspace, and a 1 ¯ = 010 ¯ and a 2 ¯ = 100 ¯ in the odd subspace. Finally, in each subspace, there is one doubly excited state, missing two non-zero bond Fermion with excitation energy 4 | t | , a 1 a 2 ¯ = 000 in the even subspace, and a 1 a 2 = 001 in the odd subspace. Table 1 summarizes the description of the spectrum in terms of bond Fermions.

4. Kitaev chain and a light induced valence hole

Absorption of a photon injects an electron-hole pair into the system. Therefore, the relevant optically excited states live in the subspace of all configurations with one hole. Here, after studying the energy spectrum of the full Hamiltonian in Equation (1) with one hole in the configuration space of bond Fermions, we discuss the absorption spectrum of the chain and the optical signature of the MZM.

4.1. Exact diagonalization of electron-hole system

Having demonstrated the benefit of bond Fermion basis, we now study the Hamiltonian with one hole in the configuration basis of bond Fermions and one hole. Using M , p M ; m to refer to M bond Fermions being in their p M configuration, and the hole being at site m, we can find the spectrum by solving an equation similar to Equation (4), but considering the full Hamiltonian H in Equation (1) rather than H e , and in the configuration basis of bond Fermions and one hole.
For instance, for N = 3 dots, following the convention we introduced in Section 3.3 and Table 1, we can list these configurations as
In this subspace, the hole Hamiltonian H h in Equation (1b) amounts to a constant η and mixes states with the same electronic configurations and different locations of the hole by hopping matrix element τ . Therefore, with the ordering in Table 2, the full Hamiltonian with one hole for the example of N = 3 dots has a structure like
H = H 1 τ 0 τ H 2 τ 0 τ H 3 + η ,
where each block is a 4 × 4 matrix, τ is the identity matrix times τ , and the diagonal blocks are given by the matrix elements of H e + H i n t in Equation (1) over the configurations in Table 2. The interaction term H i n t in Equation (1c) for each of the diagonal blocks H j is V n j e , and it mixes up different bond Fermion configurations as we have
n j e = 1 2 + 1 2 ( a j 1 a j + a j 1 a j + h . c . ) , 1 < j N ,
n 1 e = 1 2 + 1 2 ( a N a 1 + a N a 1 + h . c . ) , j = 1 ,
which implies that when the hole is not at the two ends of the chain then the interaction mixes up two non-zero bond Fermions, and when it is at one of the two ends, the interaction mixes the zero mode with one of the non-zero ones. For instance, for the operator n 1 e and n 2 e in the even configuration basis in Table 2 we have
n 1 e = 1 2 1 1 0 0 1 1 0 0 0 0 1 1 0 0 1 1 , n 2 e = 1 2 1 0 0 1 0 1 1 0 0 1 1 0 1 0 0 1 ,
while n 3 e , in a similar fashion to n 1 e , mixes GS ¯ with a 2 , and a 1 with a 1 a 2 ¯ .

4.2. Energy spectrum of the electron-hole system

Figure 4 shows the energy spectrum of a chain of length N = 3 dots in the even subspace and for Δ = t , and μ = 0 , as the electron-hole interaction V increases; for a localized hole ( τ = 0 ) on the left panel, and for a mobile hole with τ = 0 . 3 | t | on the right panel. Both cases show branching into two groups, pertaining to bonding and antibonding pairs of states, mixed by the interaction V.
The case of localized hole allows us to understand the spectrum better. There are four states associated with each dot, and since the two end dots are geometrically the same, the spectrum always shows four pairs of doubly degenerate states. As we show in the Appendix B, and it can be seen from Equation (13), for the two end dots two of these four states are mixtures of GS ¯ ; 1 ( 3 ) and a 1 ( 2 ) ; 1 ( 3 ) that give us visible peaks at E ± described in Section 4.3.1, and also indicated on the plot. The two other pairs of degenerate levels are mixtures of a 2 ( 1 ) ; 1 ( 3 ) and a 1 a 2 ¯ ; 1 ( 3 ) , which do not get excited by absorbing a photon. For the middle dot, as can be seen from Equation (13), one pair of states are mixture of a 1 ; 2 and a 2 ; 2 , where only the bonded state gets excited by absorbing a photon (see Appendix B), resulting in the peak E 0 ; also described in Section 4.3.1 and shown on the plot. Finally, the last pair of states, which also do not get excited by absorbing a photon, are mixtures of GS ¯ ; 2 and a 1 a 2 ¯ ; 1 ( 3 ) .
As can be seen in the right panel of Figure 4, for a mobile hole, when τ 0 , we still have four pairs of doubly degenerate states as a result of the chain’s spatial symmetry. At V = 0 for the case of localized hole there is an extra triple degeneracy because of the non-dispersive nature of the localized hole band. But for a mobile hole, it can be seen on the right panel of Figure 4, that at V = 0 the degenerate levels split into sets of triples, corresponding to the three propagating modes of the hole band. More importantly, in this case, since the three dot subspaces are connected by hole hopping (see Equation (11)), the above described pairs of states mix up by τ , and the ones that are closer in energy mix more. As a result of this mixture, more peaks arise in the absorption spectrum, as we discuss in the next section.

4.3. Absorption spectrum

As the schematic in Figure 1(a) shows, in absorption experiment a photon probes the chain along the nanowire. InAsP dots having significantly smaller bandgap than InP bulk of the nanowire [48], guaranties that the photon can only be absorbed by the dots. For calculating the absorption spectrum of the chain, we assume that the photon creates an electron-hole pair with uniform probability along the nanowire, and so define the polarization operator as
P = 1 N i = 1 N c i h i = 1 N i = 1 N P i ,
where we also introduced the local electron-hole pair operator P i = c i h i .
We are assuming that one can also setup the system to create the electron-hole pair on a chosen specific dot i [35,49,50,51], i.e. acting with the operator P i on the chain, rather than P. As we discuss, having access to such spatially resolved spectrum is important in detecting the optical signature of the MZM.
The polarization operator P ( i ) – we use this notation to simultaneously refer to P and P i – takes the ground state of the system to an excited state with one hole and one electron in CB. Since the ground state can be degenerate, as it is when Δ = t and μ = 0 , the absorption spectrum has an even and an odd part pertaining to each ground state
A ( i ) ( E ) = | β even | 2 ϕ odd | ϕ odd | P ( i ) GS even | 2 δ ( E E ϕ odd + E GS ) + | β odd | 2 ϕ even | ϕ even | P ( i ) GS odd | 2 δ ( E E ϕ even + E GS ) = | β even | 2 A ( i ) even ( E ) + | β odd | 2 A ( i ) odd ( E ) ,
where ϕ even ( odd ) are the eigenstates of the one hole subspace and the corresponding electron parity, and we used the notation A ( i ) to simultaneously refer to the regular absorption spectrum A, and A i the spatially resolved absorption spectrum coming from dot i.

4.3.1. Analytic result for localized hole

If τ = 0 and the hole is localized the full Hamiltonian becomes block diagonal (see Equation (11)), i.e. the subspaces of having the hole in each of the dots decouple. Consequently, we have
A ( E ) = 1 N i = 1 N A i ( E ) .
At the heart of topological regime when Δ = t and μ = 0 , as depicted graphically in Figure 5, and expressed in Equation (A6) , an electron created at site i by c i , decomposes into a superposition of creation and annihilation operators of the two bond Fermions on its two sides, a i ( ) and a i 1 ( ) . If the electron is created at one of the two ends, one of the bond Fermions is the zero mode a N . On the other hand, when the hole is at site i, the interaction V n i e n i h mixes the two bond Fermions, as shown in Figure 5, and expressed in Equations (A8) and (A12). As we show in the Appendix B, combining these two mechanisms, one can find an analytic expression for the spatially resolved absorption spectrum for a chain of arbitrary length N if the hole is created on site i as
A i ( E ) = 1 2 δ ( E E 0 ) 1 < i < N A δ ( E E ) + A + δ ( E E + ) i = 1 , N ,
where
E 0 = η + 2 | t | V ,
E ± = η + | t | V 2 ± t 2 + V 2 2 ,
A ± = 1 2 1 V 4 t 2 + V 2 ,
and then the full absorption spectrum is given by the simple sum in Equation (16).
The bottom row of Figure 6 shows the results in Equations (17) and (18) for the case of N = 3 . The peak E 0 is only present in the middle, while the peaks E ± are present at the two ends of the chain. As we show in the Appendix B, the two peaks E ± have a mixture of zero mode in them, while E 0 is purely made of non-zero bond Fermions. At V = 0 , E is purely made of zero mode while E + is purely made of non-zero bond Fermions. As we increase V, E + acquires more zero mode contribution while E mixes more with a non-zero bond Fermions. At the same time, by increasing V, the peak at E + diminishes, as can be seen from Equation (18c). If not too weak, E + peak is a better resolved optical signature for the MZM than E , as it is separated from the rest of the spectrum by V, and we expect to have V | t | . This presents an advantage over scanning tunnelling microscopy approach for detecting MZM [21]. Moreover, if one can perform spatially resolved absorption spectroscopy on the chain, the presence of the zero mode can be determined by the presence of a visible peak at high energy near E + when probing the end dots, and its absence when probing other dots.

4.3.2. Absorption for mobile hole

When the hole is mobile, there are N itinerant hole states with different energies. Therefore, one would expect N different transitions to each electronic state. More importantly, as can be seen in Equation (11), hopping hole mixes up different subspaces of having the hole in different dots. As a result more transitions become allowed leading to the emergence of more peaks in the absorption spectrum.
In Figure 6 we compare the absorption spectrum of a mobile hole with τ = 0 . 1 | t | (top row), and the analytic result of Equation (17) for localized hole (bottom row), for the case of N = 3 . It is evident how more peaks are visible for the case of mobile hole, while the major peaks are still close to the location of E 0 and E ± . Moreover, note how in the full spectrum A ¯ (top left) there is only one visible peak at high energy near E + , and how the same peak is large in A ¯ 1 (top middle) and faint in A ¯ 2 (top right), pertaining to the localized nature of the MZM that E + carries. In plotting Figure 6 we used A ¯ ( i ) = 1 2 ( A ( i ) even + A ( i ) odd ) , as for a mobile hole the even and odd parts of the absorption spectrum are not the same. But since there is no preference between the two ground states, one would expect to observe an average of the two.
The same logic is valid for a chain of any length, as the analytic result in Equation (17) is for general N. Figure 7 shows the absorption spectrum of a chain of length N = 9 . Here, we set V = 10 | t | while changing τ . When τ 0 we approach the idealized case of localized hole, where the subspaces of having the hole on each dot are decoupled. Growing τ mixes up the modes of different dots. Consequently, the zero mode starts leaking out of the two ends of the chain. On the first panel of Figure 7 we can see that at high energy there is still only one visible peak near E + , until about τ = 0 . 3 | t | a faint peak appears to the right of it. This makes E + a very robust signature for a relatively large range of hole hopping. Having access to spatially resolved spectrum we can further confirm that the peak is indeed coming from the two ends. It can be seen from the third panel of Figure 7 that the there is no visible high energy peak on the site next to the end dot ( A ¯ 2 ) until around τ = 0 . 1 | t | . In contrast, we can observe in A ¯ 2 that E , which also contains a large share of zero mode and is a stronger peak, starts leaking out of the end dot very quickly for small τ ’s.

5. Conclusion

We present here a theory of Majorana excitons, photo-excited conduction electron-valence band hole pairs, interacting with Majorana Fermions in a Kitaev chain of semiconductor quantum dots embedded in a nanowire. Using exact diagonalization techniques and Majorana and Bond Fermions we compute the energy spectra of the system. We confirm the existence of nonlocal bond Fermion, a superposition of Majorana Fermions at the two ends of the chain, with zero energy. We introduce a valence band hole and describe its interaction with Majorana fermions. We predict interband absorption spectra and discuss the signature of Majorana Zero Modes in the absorption spectra. We demonstrate how spatially resolved absorption spectrum can be used to confirm the localized character of the MZMs.
We hope this preliminary work motivates future theoretical and experimental work on hybrid nanowire semiconductor quantum dots /superconductor systems for the demonstration of Majorana Fermions.

Author Contributions

Conceptualization, P.H. and M.K.; methodology, P.H.; software, M.M., D.M., and H.A.; validation, M.M., D.M., H.A., and D.G.; formal analysis, H.A., M.M., and D.G.; investigation, P.H., M.M., and H.A.; resources, P.H., D.M.; data curation, H.A. and M.M.; writing—original draft preparation, H.A. and M.M; writing—review and editing, P.H., H.A., M.M., and D.M.; visualization, M.M. and H.A.; supervision, P.H.; project administration, P.H.; funding acquisition, P.H. and M.K.. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Quantum Sensors and Applied Quantum Computing Challenge Programs at the National Research Council of Canada, by NSERC Discovery Grant No. RGPIN- 2019-05714, and University of Ottawa Research Chair in Quantum Theory of Quantum Materials, Nanostructures, and Devices.

Data Availability Statement

The codes used to produce data presented in this study are openly available in "Kitaev Exciton" at reference [52].

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Exact diagonalization for chain of length three

Following Equation (2), for a chain of length N = 3 , the configurations of normal Fermions in the even and odd subspaces are
Even : { 000 , 110 , 101 , 011 } ,
Odd : { 100 , 010 , 001 , 111 } .
Computing the matrix elements of all terms in the Kitaev Hamiltonian H e in Equation (1a) between every pairs of configurations in Equation (A1a) and Equation (A1b) we can explicitly write the matrix q M , M H e p M , M in Equation (4), in each of the subspaces as
H c even = 0 Δ 0 Δ Δ 2 μ t 0 0 t 2 μ t Δ 0 t 2 μ ,
H c odd = μ t 0 Δ t μ t 0 0 t μ Δ Δ 0 Δ 3 μ ,
for which we used the same ordering as in Equation (A1).
Then, for ED in the bond Fermion basis, following Equation (10), the configurations are
Even : { 000 ¯ , 110 ¯ , 101 ¯ , 011 ¯ } ,
Odd : { 100 ¯ , 010 ¯ , 001 ¯ , 111 ¯ } ,
and one can use Equation (8) to find the corresponding Kitaev Hamiltonian matrices in each of the subspaces as
H a even = 1 2 2 ( t + Δ ) μ t Δ + μ μ t + Δ μ 2 ( t + Δ ) t Δ μ μ t + Δ t Δ + μ t Δ μ 0 μ μ t + Δ μ t + Δ μ 0 3 2 μ ,
H a odd = 1 2 0 μ t Δ μ Δ t μ μ 0 t Δ μ Δ t μ t Δ μ t Δ μ 2 ( t + Δ ) μ Δ t μ Δ t μ μ 2 ( t + Δ ) 3 2 μ ,
for which we used the same ordering as in Equation (A3). Notice how the two Hamiltonian matrices in Equation (A4) become diagonal when Δ = t and μ = 0 .

Appendix B. Analytic calculation of absorption spectrum for localized hole

Following the notation in Table 1, for a chain of arbitrary length N, and for Δ = t < 0 and μ = 0 , we express the two degenerate ground states of the system as
GS = j = 1 N a j | 0 a ,
GS ¯ = j = 1 N 1 a j | 0 a .
Here, we derive the result in Equations (17) and (18) using GS in Equation (A5a), and the procedure is the same for GS ¯ .
To start, first note that from Equation (7) we have
c j = 1 2 ( a j + a j + a j 1 a j 1 ) , 1 < j N
c 1 = 1 2 ( a 1 + a 1 + a N a N ) , j = 1 .
which means for 1 < j < N we have
P j GS = c j h j GS = 1 2 ( a j a j 1 ) GS ; j = 1 2 ( a j ; j a j 1 ; j ) ,
where we used the same notation as in Table 2 for the excited states. On the other hand, from Equation (12a) we have
n j e a j 1 a j = 1 2 1 1 1 1 a j 1 a j .
Recalling that both a j and a j 1 are eigenstates of H e with excitation energy 2 | t | , then in the basis { a j ; j , a j 1 ; j } , the full Hamiltonian H = H e + η V n j e is
H = 2 | t | + η V 2 1 1 1 1 ,
with the following two eigenstates
a j ; j = 1 2 a j ; j a j 1 ; j , E 0 = η + 2 | t | V
a j + ; j = 1 2 a j ; j + a j 1 ; j , E 1 = η + 2 | t | .
Now using Equation (A7) we have a j ; j | P j GS 2 = 1 2 , and a j + ; j | P j GS 2 = 0 , which give us the result in Equation (17) for the case of 1 < j < N .
Next, we show the second case for j = 1 , and the procedure is the same for j = N . In this case, using Equation (A6b) we have
P 1 GS = c 1 h 1 GS = 1 2 ( a 1 a N ) GS ; 1 = 1 2 ( a 1 ; 1 a N ; 1 ) ,
where a N is identical to the other ground state GS ¯ up to a global phase, hence its excitation energy is zero. Then from Equation (12b) we have
n 1 e a N a 1 = 1 2 1 1 1 1 a N a 1 .
Therefore, considering that a 1 and a N are eigenstates of H e with excitation energy 2 | t | and zero, respectively, the full Hamiltonian H = H e + η V n 1 e in the basis { a N ; 1 , a 1 ; 1 } is
H = η V 2 + 0 V 2 V 2 2 | t | ,
and its two eigenstates are given by
a 1 ; 1 = cos ( θ ) a 1 ; 1 sin ( θ ) a N ; 1 , E = η + | t | V 2 t 2 + V 2 2 ,
a 1 + ; 1 = sin ( θ ) a 1 ; 1 + cos ( θ ) a N ; 1 , E + = η + | t | V 2 + t 2 + V 2 2 ,
with
cos ( θ ) = 1 2 + t 2 4 t 2 + V 2 , sin ( θ ) = 1 2 t 2 4 t 2 + V 2 .
Now using Equation (A11) we have a 1 ; 1 P 1 GS 2 = 1 2 ( cos ( θ ) + sin ( θ ) ) 2 = 1 2 A , and a 1 + ; 1 P 1 GS 2 = 1 2 ( cos ( θ ) sin ( θ ) ) 2 = 1 2 A + , where we used Equation (A15) to express these matrix elements in terms of A ± in Equation (18c). From this follows the result in Equation (17) for the case of j = 1 , N .

References

  1. Gyongyosi, L.; Imre, S. A survey on quantum computing technology. Computer Science Review 2019, 31, 51–71. doi:10.1016/j.cosrev.2018.11.002. [CrossRef]
  2. Field, B.; Simula, T. Introduction to topological quantum computation with non-Abelian anyons. Quantum Science and Technology 2018, 3, 045004. doi:10.1088/2058-9565/aacad2. [CrossRef]
  3. Campbell, E.T.; Terhal, B.M.; Vuillot, C. Roads towards fault-tolerant universal quantum computation. Nature 2017, 549, 172–179. doi:10.1038/nature23460. [CrossRef]
  4. Stern, A.; Lindner, N.H. Topological quantum computation—from basic concepts to first experiments. Science 2013, 339, 1179–1184. doi:10.1126/science.1231473. [CrossRef]
  5. Nayak, C.; Simon, S.H.; Stern, A.; Freedman, M.; Sarma, S.D. Non-Abelian anyons and topological quantum computation. Reviews of Modern Physics 2008, 80, 1083. doi:10.1103/revmodphys.80.1083. [CrossRef]
  6. Sarma, S.D.; Freedman, M.; Nayak, C. Majorana zero modes and topological quantum computation. npj Quantum Information 2015, 1, 1–13. doi:10.1038/npjqi.2015.1. [CrossRef]
  7. Das Sarma, S.; Freedman, M.; Nayak, C. Topological quantum computation. Physics Today 2006, 59, 32–38. doi:10.1063/1.2337825. [CrossRef]
  8. Freedman, M.; Kitaev, A.; Larsen, M.; Wang, Z. Topological quantum computation. Bulletin of the American Mathematical Society, 40, 31–38. doi:10.48550/arXiv.quant-ph/0101025. [CrossRef]
  9. Jaworowski, B.; Hawrylak, P. Quantum bits with macroscopic topologically protected states in semiconductor devices. Applied Sciences 2019, 9, 474. doi:10.3390/app9030474. [CrossRef]
  10. Haldane, F.D.M. Nonlinear field theory of large-spin Heisenberg antiferromagnets: semiclassically quantized solitons of the one-dimensional easy-axis Néel state. Physical Review Letters 1983, 50, 1153. doi:10.1103/PhysRevLett.50.1153. [CrossRef]
  11. Kitaev, A.Y. Unpaired Majorana fermions in quantum wires. Physics-Uspekhi 2001, 44, 131. doi:10.1070/1063-7869/44/10s/s29. [CrossRef]
  12. Kitaev, A.Y. Fault-tolerant quantum computation by anyons. Annals of Physics 2003, 303, 2–30. doi:10.1016/s0003-4916(02)00018-0. [CrossRef]
  13. Lutchyn, R.M.; Sau, J.D.; Sarma, S.D. Majorana fermions and a topological phase transition in semiconductor-superconductor heterostructures. Physical Review Letters 2010, 105, 077001. doi:10.1103/PhysRevLett.105.077001. [CrossRef]
  14. Mourik, V.; Zuo, K.; Frolov, S.M.; Plissard, S.; Bakkers, E.P.; Kouwenhoven, L.P. Signatures of Majorana fermions in hybrid superconductor-semiconductor nanowire devices. Science 2012, 336, 1003–1007. doi:10.1126/science.1222360. [CrossRef]
  15. Sau, J.D.; Sarma, S.D. Realizing a robust practical Majorana chain in a quantum-dot-superconductor linear array. Nature Communications 2012, 3, 964. doi:10.1038/ncomms1966. [CrossRef]
  16. Leijnse, M.; Flensberg, K. Parity qubits and poor man’s Majorana bound states in double quantum dots. Physical Review B 2012, 86, 134528. doi:10.1103/PhysRevB.86.134528. [CrossRef]
  17. Dvir, T.; Wang, G.; van Loo, N.; Liu, C.X.; Mazur, G.P.; Bordin, A.; Ten Haaf, S.L.; Wang, J.Y.; van Driel, D.; Zatelli, F.; others. Realization of a minimal Kitaev chain in coupled quantum dots. Nature 2023, 614, 445–450. doi:10.1038/s41586-022-05585-1. [CrossRef]
  18. Nadj-Perge, S.; Drozdov, I.K.; Li, J.; Chen, H.; Jeon, S.; Seo, J.; MacDonald, A.H.; Bernevig, B.A.; Yazdani, A. Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor. Science 2014, 346, 602–607. doi:10.1126/science.1259327. [CrossRef]
  19. Sun, H.H.; Jia, J.F. Detection of Majorana zero mode in the vortex. npj Quantum Materials 2017, 2, 34. doi:10.1038/s41535-017-0037-4. [CrossRef]
  20. Liu, D.E.; Baranger, H.U. Detecting a Majorana-fermion zero mode using a quantum dot. Physical Review B 2011, 84, 201308. doi:10.1103/physrevb.84.201308. [CrossRef]
  21. Jäck, B.; Xie, Y.; Yazdani, A. Detecting and distinguishing Majorana zero modes with the scanning tunnelling microscope. Nature Reviews Physics 2021, 3, 541–554. doi:10.1038/s42254-021-00328-z. [CrossRef]
  22. Pientka, F.; Romito, A.; Duckheim, M.; Oreg, Y.; von Oppen, F. Signatures of topological phase transitions in mesoscopic superconducting rings. New Journal of Physics 2013, 15, 025001. doi:10.1088/1367-2630/15/2/025001. [CrossRef]
  23. Pikulin, D.I.; van Heck, B.; Karzig, T.; Martinez, E.A.; Nijholt, B.; Laeven, T.; Winkler, G.W.; Watson, J.D.; Heedt, S.; Temurhan, M.; others. Protocol to identify a topological superconducting phase in a three-terminal device. arXiv preprint arXiv:2103.12217 2021. doi:10.48550/arXiv.2103.12217. [CrossRef]
  24. Liu, J.; Potter, A.C.; Law, K.T.; Lee, P.A. Zero-bias peaks in the tunneling conductance of spin-orbit-coupled superconducting wires with and without Majorana end-states. Physical Review Letters 2012, 109, 267002. doi:10.1103/PhysRevLett.109.267002. [CrossRef]
  25. Sarma, S.D.; Pan, H. Disorder-induced zero-bias peaks in Majorana nanowires. Physical Review B 2021, 103, 195158. doi:10.1103/PhysRevB.103.195158. [CrossRef]
  26. Rubbert, S.; Akhmerov, A. Detecting Majorana nonlocality using strongly coupled Majorana bound states. Physical Review B 2016, 94, 115430. doi:10.1103/physrevb.94.115430. [CrossRef]
  27. Aghaee, M.; Akkala, A.; Alam, Z.; Ali, R.; Ramirez, A.A.; Andrzejczuk, M.; Antipov, A.E.; Astafev, M.; Bauer, B.; Becker, J.; others. InAs-Al hybrid devices passing the topological gap protocol. arXiv preprint arXiv:2207.02472. doi:10.48550/arXiv.2207.02472. [CrossRef]
  28. Baldelli, N.; Bhattacharya, U.; González-Cuadra, D.; Lewenstein, M.; Graß, T. Detecting Majorana zero modes via strong field dynamics. ACS Omega 2022. doi:10.1021/acsomega.2c07169. [CrossRef]
  29. Cygorek, M.; Korkusinski, M.; Hawrylak, P. Atomistic theory of electronic and optical properties of InAsP/InP nanowire quantum dots. Physical Review B 2020, 101, 075307. doi:10.1103/physrevb.101.075307. [CrossRef]
  30. Manalo, J.; Cygorek, M.; Altintas, A.; Hawrylak, P. Electronic and magnetic properties of many-electron complexes in charged InAsxP1-x quantum dots in InP nanowires. Physical Review B 2021, 104. doi:10.1103/PhysRevB.104.125402. [CrossRef]
  31. Koong, Z.X.; Ballesteros-Garcia, G.; Proux, R.; Dalacu, D.; Poole, P.J.; Gerardot, B.D. Multiplexed single photons from deterministically positioned nanowire quantum dots. Physical Review Applied 2020, 14, 034011. doi:10.1103/PhysRevApplied.14.034011. [CrossRef]
  32. Dalacu, D.; Mnaymneh, K.; Lapointe, J.; Wu, X.; Poole, P.J.; Bulgarini, G.; Zwiller, V.; Reimer, M.E. Ultraclean emission from InAsP quantum dots in defect-free wurtzite InP nanowires. Nano Letters 2012, 12, 5919–5923. doi:10.1021/nl303327h. [CrossRef]
  33. Jaworowski, B.; Rogers, N.; Grabowski, M.; Hawrylak, P. Macroscopic singlet-triplet qubit in synthetic spin-one chain in semiconductor nanowires. Scientific Reports 2017, 7, 5529. doi:10.1038/s41598-017-05655-9. [CrossRef]
  34. Phoenix, J.; Korkusinski, M.; Dalacu, D.; Poole, P.J.; Zawadzki, P.; Studenikin, S.; Williams, R.L.; Sachrajda, A.S.; Gaudreau, L. Magnetic tuning of tunnel coupling between InAsP double quantum dots in InP nanowires. Scientific Reports 2022, 12, 5100. doi:10.1038/s41598-022-08548-8. [CrossRef]
  35. Northeast, D.B.; Weber, J.F.; Dalacu, D.; Phoenix, J.; Poole, P.J.; Aers, G.; Lapointe, J.; Williams, R.L. Optical fibre-based (plug-and-play) single photon source using InAsP quantum dot nanowires and gradient-index lens collection. arXiv preprint arXiv:2104.11197. doi:10.1038/s41598-021-02287-y. [CrossRef]
  36. Laferrière, P.; Yeung, E.; Korkusinski, M.; Poole, P.J.; Williams, R.L.; Dalacu, D.; Manalo, J.; Cygorek, M.; Altintas, A.; Hawrylak, P. Systematic study of the emission spectra of nanowire quantum dots. Applied Physics Letters 2021, 118, 161107. doi:10.1063/5.0045880. [CrossRef]
  37. Talantsev, E.; Iida, K.; Ohmura, T.; Matsumoto, T.; Crump, W.; Strickland, N.; Wimbush, S.; Ikuta, H. p-wave superconductivity in iron-based superconductors. Scientific Reports 2019, 9, 14245. doi:10.1038/s41598-019-50687-y. [CrossRef]
  38. Wang, W.S.; Zhang, C.C.; Zhang, F.C.; Wang, Q.H. Theory of chiral p-wave superconductivity with near nodes for Sr2RuO4. Physical Review Letters 2019, 122, 027002. doi:10.1103/PhysRevLett.122.027002. [CrossRef]
  39. Yuan, N.F.; Mak, K.F.; Law, K. Possible topological superconducting phases of MoS2. Physical Review Letters 2014, 113, 097001. doi:10.1103/PhysRevLett.113.097001. [CrossRef]
  40. Frigeri, P.; Agterberg, D.; Koga, A.; Sigrist, M. Superconductivity without Inversion Symmetry: MnSi versus CePt3Si. Physical Review Letters 2004, 92, 097001. doi:10.1103/PhysRevLett.92.097001. [CrossRef]
  41. Hardy, F.; Huxley, A. p-wave superconductivity in the ferromagnetic superconductor URhGe. Physical Review Letters 2005, 94, 247006. doi:10.1103/PhysRevLett.94.247006. [CrossRef]
  42. Ishida, K.; Mukuda, H.; Kitaoka, Y.; Asayama, K.; Mao, Z.; Mori, Y.; Maeno, Y. Spin-triplet superconductivity in Sr2RuO4 identified by 17O Knight shift. Nature 1998, 396, 658–660. doi:10.1038/25315. [CrossRef]
  43. Mahan, G. Many-Particle Physics; Physics of Solids and Liquids, 2012. doi:10.1007/978-1-4757-5714-9. [CrossRef]
  44. Wojs, A.; Hawrylak, P. Negatively charged magnetoexcitons in quantum dots. Physical Review B 1995, 51, 10880–10885. doi:10.1103/PhysRevB.51.10880. [CrossRef]
  45. Hawrylak, P. Excitonic effects in optical spectra of a quasi-one-dimensional electron gas. Solid State Communications 1992, 81, 525–527. doi:https://doi.org/10.1016/0038-1098(92)90605-9. [CrossRef]
  46. Hawrylak, P. Optical properties of a two-dimensional electron gas: Evolution of spectra from excitons to Fermi-edge singularities. Physical Review B 1991, 44, 3821–3828. doi:10.1103/PhysRevB.44.3821. [CrossRef]
  47. Weiße, A.; Fehske, H., Exact Diagonalization Techniques. In Computational Many-Particle Physics; Fehske, H.; Schneider, R.; Weiße, A., Eds.; Springer Berlin Heidelberg: Berlin, Heidelberg, 2008; pp. 529–544. doi:10.1007/978-3-540-74686-7_18. [CrossRef]
  48. Adachi, S., Energy-Band Structure: Energy-Band Gaps. In Properties of Group-IV, III-V and II-VI Semiconductors; John Wiley & Sons, Ltd, 2005; chapter 6, p. 116. doi:10.1002/0470090340.ch6. [CrossRef]
  49. Dalacu, D.; Poole, P.J.; Williams, R.L. Tailoring the geometry of bottom-up nanowires: application to high efficiency single photon sources. Nanomaterials 2021, 11. doi:10.3390/nano11051201. [CrossRef]
  50. Laferriere, P.; Yeung, E.; Giner, L.; Haffouz, S.; Lapointe, J.; Aers, G.C.; Poole, P.J.; Williams, R.L.; Dalacu, D. Multiplexed single-photon source based on multiple quantum dots embedded within a single nanowire. Nano Letters 2020, 20, 3688–3693. doi:10.1021/acs.nanolett.0c00607. [CrossRef]
  51. Laferriére, P.; Haffouz, S.; Northeast, D.B.; Poole, P.J.; Williams, R.L.; Dalacu, D. Position-controlled telecom single photon emitters operating at elevated temperatures. Nano Letters 2023, 23, 962–968. doi:10.1021/acs.nanolett.2c04375. [CrossRef]
  52. Allami, H. Kitaev Exciton. Available at https://github.com/hassan-allami/KitaevExciton.git, 2023.
Figure 1. (a) Schematic of the system and the light absorption experiment. (b) Schematic of the Hamiltonian terms between two adjacent dots according to Equation (1), where conduction(valence) levels are labeled by c i ( h i ) operators. The conduction level is the reference of energy, hence the downward arrow indicates negative μ .
Figure 1. (a) Schematic of the system and the light absorption experiment. (b) Schematic of the Hamiltonian terms between two adjacent dots according to Equation (1), where conduction(valence) levels are labeled by c i ( h i ) operators. The conduction level is the reference of energy, hence the downward arrow indicates negative μ .
Preprints 78616 g001
Figure 2. Schematic of Kitaev chain in the Majorana and bond representation, with non-zero bond Fermions in purple, and the nonlocal zero mode a N living on the two ends of the chain.
Figure 2. Schematic of Kitaev chain in the Majorana and bond representation, with non-zero bond Fermions in purple, and the nonlocal zero mode a N living on the two ends of the chain.
Preprints 78616 g002
Figure 3. Energy spectra of Kitaev chain in normal(left) and bond(right) basis, where Δ = t < 0 and μ = 0 . Energy is normalized to | t | .
Figure 3. Energy spectra of Kitaev chain in normal(left) and bond(right) basis, where Δ = t < 0 and μ = 0 . Energy is normalized to | t | .
Preprints 78616 g003
Figure 4. Energy spectra of the full Hamiltonian Equation (1) with one hole in the even subspace, as a function of electron-hole interaction V, for N = 3 dots, Δ = t and μ = 0 . (left) for the case of localized hole, τ = 0 , (right) for a mobile hole with τ = 0 . 3 | t | . The overlap of transparent markers makes the degenerate levels look darker. The peak energies E 0 and E ± discussed in Section 4.3.1 are also shown according to Equations (18a) and (18b).
Figure 4. Energy spectra of the full Hamiltonian Equation (1) with one hole in the even subspace, as a function of electron-hole interaction V, for N = 3 dots, Δ = t and μ = 0 . (left) for the case of localized hole, τ = 0 , (right) for a mobile hole with τ = 0 . 3 | t | . The overlap of transparent markers makes the degenerate levels look darker. The peak energies E 0 and E ± discussed in Section 4.3.1 are also shown according to Equations (18a) and (18b).
Preprints 78616 g004
Figure 5. An electron created by c i is a superposition of creation and annihilation operators of two bond Fermions a i ( ) and a i 1 ( ) , according to Equation (A6). The interaction V n i e n i h mixes up the two bond Fermions according to Equation (A8). Note that when i is one of the two ends, then one of the bond Fermions is the zero mode a N (see the Appendix B for more details).
Figure 5. An electron created by c i is a superposition of creation and annihilation operators of two bond Fermions a i ( ) and a i 1 ( ) , according to Equation (A6). The interaction V n i e n i h mixes up the two bond Fermions according to Equation (A8). Note that when i is one of the two ends, then one of the bond Fermions is the zero mode a N (see the Appendix B for more details).
Preprints 78616 g005
Figure 6. (left) The averaged absorption spectrum A ¯ ( E ) , and (middle and right) spatially resolved absorption A ¯ i ( E ) , for Δ = t , μ = 0 , and for N = 3 dots. (top) for a mobile hole with τ = 0 . 1 | t | , (bottom) for a localized hole, τ = 0 , according to the analytic results in Equations (17) and (18). The spectra are plotted against ( E η ) / | t | while changing V / | t | on the y-axis. The bright curves show the location of the peaks as V changes, and the color scale shows their heights. Gaussian profile was used for the peaks with the width σ = 0 . 025 | t | . The maximum value of each peak shows the magnitude of the corresponding matrix element.
Figure 6. (left) The averaged absorption spectrum A ¯ ( E ) , and (middle and right) spatially resolved absorption A ¯ i ( E ) , for Δ = t , μ = 0 , and for N = 3 dots. (top) for a mobile hole with τ = 0 . 1 | t | , (bottom) for a localized hole, τ = 0 , according to the analytic results in Equations (17) and (18). The spectra are plotted against ( E η ) / | t | while changing V / | t | on the y-axis. The bright curves show the location of the peaks as V changes, and the color scale shows their heights. Gaussian profile was used for the peaks with the width σ = 0 . 025 | t | . The maximum value of each peak shows the magnitude of the corresponding matrix element.
Preprints 78616 g006
Figure 7. Absorption spectrum for a chain of length N = 9 , and for Δ = t , μ = 0 , V = 10 | t | , and changing τ . (left) The full averaged spectrum A ¯ , (middle) the spatially resolved spectrum for the first dot A ¯ 1 , and (right) the spatially resolved spectrum for the second dot A ¯ 2 . The bright curves show the location of the peaks as τ changes, and the colorscale shows their heights. Gaussian profile was used for the peaks with the width σ = 0 . 025 | t | . The maximum value of each peak shows the magnitude of the corresponding matrix element.
Figure 7. Absorption spectrum for a chain of length N = 9 , and for Δ = t , μ = 0 , V = 10 | t | , and changing τ . (left) The full averaged spectrum A ¯ , (middle) the spatially resolved spectrum for the first dot A ¯ 1 , and (right) the spatially resolved spectrum for the second dot A ¯ 2 . The bright curves show the location of the peaks as τ changes, and the colorscale shows their heights. Gaussian profile was used for the peaks with the width σ = 0 . 025 | t | . The maximum value of each peak shows the magnitude of the corresponding matrix element.
Preprints 78616 g007
Table 1. Describing the spectra plotted in Figure 3. Configurations of bond Fermions are the eigenstates of Kitaev Hamiltonian when Δ = t and μ = 0 .
Table 1. Describing the spectra plotted in Figure 3. Configurations of bond Fermions are the eigenstates of Kitaev Hamiltonian when Δ = t and μ = 0 .
index 1 2 3 4 5 6 7 8
configuration 111 ¯ 110 ¯ 010 ¯ 011 ¯ 100 ¯ 101 ¯ 001 ¯ 000 ¯
label GS GS ¯ a ¯ 1 a 1 a ¯ 2 a 2 a 1 a 2 a 1 a 2 ¯
parity odd even odd even odd even odd even
excitation
energy
0 0 2 | t | 2 | t | 2 | t | 2 | t | 4 | t | 4 | t |
Table 2. Configurations of bond Fermions with one hole for N = 3 dots.
Table 2. Configurations of bond Fermions with one hole for N = 3 dots.
Even Odd
GS ¯ ; 1 a 1 ; 1 a 2 ; 1 a 1 a 2 ¯ ; 1
GS ¯ ; 2 a 1 ; 2 a 2 ; 2 a 1 a 2 ¯ ; 2
GS ¯ ; 3 a 1 ; 3 a 2 ; 3 a 1 a 2 ¯ ; 3
GS ; 1 a ¯ 1 ; 1 a ¯ 2 ; 1 a 1 a 2 ; 1
GS ; 2 a ¯ 1 ; 2 a ¯ 2 ; 2 a 1 a 2 ; 2
GS ; 3 a ¯ 1 ; 3 a ¯ 2 ; 3 a 1 a 2 ; 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated