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

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_{int}$ 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

$$\begin{array}{cc}\hfill & H=H_e+H_h+H_{int},\hfill \end{array}$$

$$\begin{array}{cc}\hfill & H_e=t\sum _{i=1}^{N-1}{c}_{i+1}^{\mathsf{†}}{c}_i+h.c.+\mathrm{\Delta }\sum _{i=1}^{N-1}{c}_{i+1}^{\mathsf{†}}{c}_i^{\mathsf{†}}+h.c.-\mu \sum _{i=1}^N{c}_i^{\mathsf{†}}{c}_i,\hfill \end{array}$$

(1a)

$$\begin{array}{cc}\hfill & H_h=-\tau \sum _{i=1}^{N-1}{h}_{i+1}^{\mathsf{†}}{h}_i+h.c.+\eta \sum _{i=1}^N{h}_i^{\mathsf{†}}{h}_i,\hfill \end{array}$$

(1b)

$$\begin{array}{cc}\hfill & H_{int}=-V\sum _{i=1}^N{n}_i^e{n}_i^h,\hfill \end{array}$$

(1c)

where $c_i^{\mathsf{†}}\left(h_i^{\mathsf{†}}\right)$ is the normal Fermionic creation operator of an electron(hole) in dot i, $t\left(\tau \right)$ is hopping between adjacent conduction(valence) levels, $\mathrm{\Delta }$ is pairing energy between adjacent conduction levels, $\mu$ is the chemical potential measured from the conduction energy level, and $\eta$ 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^{\mathsf{†}}{c}_i(n_i^h=h_i^{\mathsf{†}}{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 $\mathrm{\Delta }=t$ and $\mu =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 $\left(\genfrac{}{}{0pt}{}{N}{0}\right)+\left(\genfrac{}{}{0pt}{}{N}{1}\right)+...+\left(\genfrac{}{}{0pt}{}{N}{N-1}\right)+\left(\genfrac{}{}{0pt}{}{N}{N}\right)=2^N$ possible configurations, which we construct as

$$\left|{\alpha }_1\cdots {\alpha }_N\right.〉=\prod _{i=1}^N{\left(c_i^{\mathsf{†}}\right)}^{{\alpha }_i}\left|0\right.〉,$$

(2)

where $\left|0\right.〉$ is the vacuum of electrons, ${\alpha }_i=1\mathrm{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

$$\left|{\psi }^{\nu }\right.〉=\sum _{M,p_M}{C}_{M,p_M}^{\nu }\left|M,p_M\right.〉,$$

(3)

where we are populating N sites with $M=0,1,\cdots ,N$ electrons. To solve for coefficients $C_{M,p_M}^{\nu }$, we apply the Hamiltonian on this state, and by using the orthogonality of the configurations we obtain the eigenvalue equation

$$\sum _{p_M,M}\left.〈q_{M^{\prime }},M^{\prime }\right|H_e\left|p_M,M\right.〉C_{M,p_M}^{\nu }=E^{\nu }{C}_{M^{\prime },q_{M^{\prime }}}^{\nu }.$$

(4)

However, since the Kitaev Hamiltonian $H_e$ in Equation (1a) only changes particle number in pairs, the matrix element $\left.〈q_{M^{\prime }},M^{\prime }\right|H_e\left|p_M,M\right.〉$ is non-zero only if M and $M^{\prime }$ 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 $\left.〈q_{M^{\prime }},M^{\prime }\right|H_e\left|p_M,M\right.〉$ 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 ${\gamma }_1$ and ${\gamma }_2$ as

$$\begin{array}{ccc}\hfill & c_j=\frac{1}{2}({\gamma }_{j,1}+i{\gamma }_{j,2}),\hfill & \hfill c_j^{\mathsf{†}}=\frac{1}{2}({\gamma }_{j,1}-i{\gamma }_{j,2}),\end{array}$$

(5)

where the $\gamma$’s are Majorana Fermion operators. Majorana Fermions satisfy a slightly different anti-commutation relation than the ordinary Fermions, $\{{\gamma }_{i,\alpha },{\gamma }_{j,\beta }\}=2{\delta }_{ij}{\delta }_{\alpha \beta }$.

Using Equation (5) and Majorana anti-commutation relations, the Hamiltonian $H_e$ can be written in terms of Majorana Fermions as

$$H_e=\frac{i}{2}\left((t+\mathrm{\Delta })\sum _{j=1}^{N-1}{\gamma }_{j,1}{\gamma }_{j+1,2}+(t-\mathrm{\Delta })\sum _{j=1}^{N-1}{\gamma }_{j+1,1}{\gamma }_{j,2}-\mu \sum _{j=1}^N({\gamma }_{j,1}{\gamma }_{j,2}-i)\right).$$

(6)

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

$$\begin{array}{cc}\hfill & a_j=\frac{1}{2}({\gamma }_{j,1}+i{\gamma }_{j+1,2})=\frac{1}{2}(c_j^{\mathsf{†}}+c_j+c_{j+1}-c_{j+1}^{\mathsf{†}}),\hfill \end{array}$$

(7a)

$$\begin{array}{cc}\hfill & a_N=\frac{1}{2}({\gamma }_{N,1}+i{\gamma }_{1,2})=\frac{1}{2}(c_N^{\mathsf{†}}+c_N+c_1-c_1^{\mathsf{†}}),\hfill \end{array}$$

(7b)

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

$$\begin{array}{cc}\hfill H_e=\frac{1}{2}& \left((t+\mathrm{\Delta })\sum _{j=1}^{N-1}(2a_j^{\mathsf{†}}{a}_j-1)+(t-\mathrm{\Delta })\sum _{j=1}^{N-1}(a_{j+1}^{\mathsf{†}}{a}_{j-1}+a_{j+1}{a}_{j-1}+h.c.)\right.\hfill \\ \hfill & \left.-\mu \sum _{j=1}^N(1+(a_j^{\mathsf{†}}{a}_{j-1}+a_j{a}_{j-1}+h.c.))\right),\hfill \end{array}$$

(8)

where in the second and the third sum one should identify $a_0\equiv a_N$. Note that in topological regime, when $t=\mathrm{\Delta }$ and $\mu =0$, the bond Fermions diagonalize the Hamiltonian in Equation (8) and reduce it to

$$H_e=t\sum _{j=1}^{N-1}(2a_j^{\mathsf{†}}{a}_j-1),$$

(9)

which implies a set of $N-1$ quasiparticles with energy $2t$, 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

$$\left|\overline{{\alpha }_1\cdots {\alpha }_N}\right.〉=\prod _{i=1}^N{\left(a_i^{\mathsf{†}}\right)}^{{\alpha }_i}\left|0_a\right.〉,$$

(10)

where $\left|0_a\right.〉$ 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 $\left|M,p_M\right.〉$ 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 $\mathrm{\Delta }=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 $\left|\mathrm{GS}\right.〉=\left|\overline{111}\right.〉$ with all bond Fermions, and the other in the even subspace $\left|\overline{\mathrm{GS}}\right.〉=\left|\overline{110}\right.〉$, which is missing the zero energy bond Fermion $a_N\equiv a_3$. Next, we have the singly excited states, missing one non-zero bond Fermion, with excitation energy $2\left|t\right|$, from which we have two in each subspace, $\left|a_1\right.〉=\left|\overline{011}\right.〉$ and $\left|a_2\right.〉=\left|\overline{101}\right.〉$ in the even subspace, and $\left|\overline{a_1}\right.〉=\left|\overline{010}\right.〉$ and $\left|\overline{a_2}\right.〉=\left|\overline{100}\right.〉$ in the odd subspace. Finally, in each subspace, there is one doubly excited state, missing two non-zero bond Fermion with excitation energy $4\left|t\right|$, $\left|\overline{a_1{a}_2}\right.〉=\left|000\right.〉$ in the even subspace, and $\left|a_1{a}_2\right.〉=\left|001\right.〉$ 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 $\left|M,p_M;m\right.〉$ 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 $\eta$ and mixes states with the same electronic configurations and different locations of the hole by hopping matrix element $\tau$. 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=\left[\begin{array}{ccc}{H}_1& -\tau & 0\\ -\tau & H_2& -\tau \\ 0& -\tau & H_3\end{array}\right]+\eta ,$$

(11)

where each block is a $4\times 4$ matrix, $\tau$ is the identity matrix times $\tau$, and the diagonal blocks are given by the matrix elements of $H_e+H_{int}$ in Equation (1) over the configurations in Table 2. The interaction term $H_{int}$ 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

$$\begin{array}{ccc}\hfill n_j^e& =\frac{1}{2}+\frac{1}{2}(a_{j-1}^{\mathsf{†}}{a}_j+a_{j-1}^{\mathsf{†}}{a}_j^{\mathsf{†}}+h.c.),\hfill & \hfill 1<j\le N,\end{array}$$

(12a)

$$\begin{array}{ccc}\hfill n_1^e& =\frac{1}{2}+\frac{1}{2}(a_N^{\mathsf{†}}{a}_1+a_N^{\mathsf{†}}{a}_1^{\mathsf{†}}+h.c.),\hfill & \hfill j=1,\end{array}$$

(12b)

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

$$\begin{array}{ccc}\hfill & n_1^e=\frac{1}{2}\left[\begin{array}{cccc}1& -1& 0& 0\\ -1& 1& 0& 0\\ 0& 0& 1& -1\\ 0& 0& -1& 1\end{array}\right],\hfill & \hfill n_2^e=\frac{1}{2}\left[\begin{array}{cccc}1& 0& 0& 1\\ 0& 1& -1& 0\\ 0& -1& 1& 0\\ 1& 0& 0& 1\end{array}\right],\end{array}$$

(13)

while $n_3^e$, in a similar fashion to $n_1^e$, mixes $\left|\overline{\mathrm{GS}}\right.〉$ with $\left|a_2\right.〉$, and $\left|a_1\right.〉$ with $\left|\overline{a_1{a}_2}\right.〉$.

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 $\mathrm{\Delta }=t$, and $\mu =0$, as the electron-hole interaction V increases; for a localized hole ($\tau =0$) on the left panel, and for a mobile hole with $\tau =0.3\left|t\right|$ 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 $\left|\overline{\mathrm{GS}};1\left(3\right)\right.〉$ and $\left|a_{1\left(2\right)};1\left(3\right)\right.〉$ that give us visible peaks at $E_{\pm }$ described in Section 4.3.1, and also indicated on the plot. The two other pairs of degenerate levels are mixtures of $\left|a_{2\left(1\right)};1\left(3\right)\right.〉$ and $\left|\overline{a_1{a}_2};1\left(3\right)\right.〉$, 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 $\left|a_1;2\right.〉$ and $\left|a_2;2\right.〉$, 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 $\left|\overline{\mathrm{GS}};2\right.〉$ and $\left|\overline{a_1{a}_2};1\left(3\right)\right.〉$.

As can be seen in the right panel of Figure 4, for a mobile hole, when $\tau \ne 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 $\tau$, 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=\frac{1}{\sqrt{N}}\sum _{i=1}^N{c}_i^{\mathsf{†}}{h}_i^{\mathsf{†}}=\frac{1}{\sqrt{N}}\sum _{i=1}^N{P}_i,$$

(14)

where we also introduced the local electron-hole pair operator $P_i=c_i^{\mathsf{†}}{h}_i^{\mathsf{†}}$.

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_{\left(i\right)}$ – 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 $\mathrm{\Delta }=t$ and $\mu =0$, the absorption spectrum has an even and an odd part pertaining to each ground state

$$\begin{array}{cc}\hfill A_{\left(i\right)}\left(E\right)& =|{\beta }_{\mathrm{even}}{|}^2\sum _{{\varphi }_{\mathrm{odd}}}|\langle {\varphi }_{\mathrm{odd}}{|P_{\left(i\right)}\left|{\mathrm{GS}}_{\mathrm{even}}\right.〉|}^2\delta (E-E_{{\varphi }_{\mathrm{odd}}}+E_{\mathrm{GS}})\hfill \\ \hfill & +|{\beta }_{\mathrm{odd}}{|}^2\sum _{{\varphi }_{\mathrm{even}}}|\langle {\varphi }_{\mathrm{even}}{|P_{\left(i\right)}\left|{\mathrm{GS}}_{\mathrm{odd}}\right.〉|}^2\delta (E-E_{{\varphi }_{\mathrm{even}}}+E_{\mathrm{GS}})\hfill \\ \hfill & =|{\beta }_{\mathrm{even}}{|}^2{A}_{\left(i\right)}^{\mathrm{even}}\left(E\right)+{|{\beta }_{\mathrm{odd}}|}^2{A}_{\left(i\right)}^{\mathrm{odd}}\left(E\right),\hfill \end{array}$$

(15)

where $\left|{\varphi }_{\mathrm{even}\left(\mathrm{odd}\right)}\right.〉$ are the eigenstates of the one hole subspace and the corresponding electron parity, and we used the notation $A_{\left(i\right)}$ 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 $\tau =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\left(E\right)=\frac{1}{N}\sum _{i=1}^N{A}_i\left(E\right).$$

(16)

At the heart of topological regime when $\mathrm{\Delta }=t$ and $\mu =0$, as depicted graphically in Figure 5, and expressed in Equation (A6) , an electron created at site i by $c_i^{\mathsf{†}}$, decomposes into a superposition of creation and annihilation operators of the two bond Fermions on its two sides, $a_i^{(\mathsf{†})}$ and $a_{i-1}^{(\mathsf{†})}$. 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\left(E\right)=\frac{1}{2}\left\{\begin{array}{cc}\delta (E-E_0)\hfill & 1<i<N\hfill \\ A_{-}\delta (E-E_{-})+A_{+}\delta (E-E_{+})\hfill & i=1,N\hfill \end{array}\right.,$$

(17)

where

$$\begin{array}{cc}\hfill & E_0=\eta +2\left|t\right|-V,\hfill \end{array}$$

(18a)

$$\begin{array}{cc}\hfill & E_{\pm }=\eta +\left|t\right|-\frac{V}{2}\pm \sqrt{t^2+{\left(\frac{V}{2}\right)}^2},\hfill \end{array}$$

(18b)

$$\begin{array}{cc}\hfill & A_{\pm }=\frac{1}{2}\left(1\mp \frac{V}{\sqrt{4t^2+V^2}}\right),\hfill \end{array}$$

(18c)

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_{\pm }$ are present at the two ends of the chain. As we show in the Appendix B, the two peaks $E_{\pm }$ 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\gg \left|t\right|$. 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 $\tau =0.1\left|t\right|$ (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_{\pm }$. Moreover, note how in the full spectrum $\overline{A}$ (top left) there is only one visible peak at high energy near $E_{+}$, and how the same peak is large in ${\overline{A}}_1$ (top middle) and faint in ${\overline{A}}_2$ (top right), pertaining to the localized nature of the MZM that $E_{+}$ carries. In plotting Figure 6 we used ${\overline{A}}_{\left(i\right)}=\frac{1}{2}(A_{\left(i\right)}^{\mathrm{even}}+A_{\left(i\right)}^{\mathrm{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\left|t\right|$ while changing $\tau$. When $\tau \to 0$ we approach the idealized case of localized hole, where the subspaces of having the hole on each dot are decoupled. Growing $\tau$ 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 $\tau =0.3\left|t\right|$ 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 (${\overline{A}}_2$) until around $\tau =0.1\left|t\right|$. In contrast, we can observe in ${\overline{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 $\tau$’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.