Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Nano-photoluminescence of natural anyon molecules and topological quantum computation

## Abstract

The proposal of fault-tolerant quantum computations, which promise to dramatically improve the operation of quantum computers and to accelerate the development of the compact hardware for them, is based on topological quantum field theories, which rely on the existence in Nature of physical systems described by a Lagrangian containing a non-Abelian (NA) topological term. These are solid-state systems having two-dimensional electrons, which are coupled to magnetic-flux-quanta vortexes, forming complex particles, known as anyons. Topological quantum computing (TQC) operations thus represent a physical realization of the mathematical operations involving NA representations of a braid group Bn, generated by a set of n localized anyons, which can be braided and fused using a “tweezer” and controlled by a detector. For most of the potential TQC material systems known so far, which are 2D-electron–gas semiconductor structure at high magnetic field and a variety of hybrid superconductor/topological-material heterostructures, the realization of anyon localization versus tweezing and detecting meets serious obstacles, chief among which are the necessity of using current control, i.e., mobile particles, of the TQC operations and high density electron puddles (containing thousands of electrons) to generate a single vortex. Here we demonstrate a novel system, in which these obstacles can be overcome, and in which vortexes are generated by a single electron. This is a ~ 150 nm size many electron InP/GaInP2 self-organized quantum dot, in which molecules, consisting of a few localized anyons, are naturally formed and exist at zero external magnetic field. We used high-spatial-resolution scanning magneto-photoluminescence spectroscopy measurements of a set of the dots having five and six electrons, together with many-body quantum mechanical calculations to demonstrate spontaneous formation of the anyon magneto-electron particles (eν) having fractional charge ν = n/k, where n = 1–4 and k = 3–15 are the number of electrons and vortexes, respectively, arranged in molecular structures having a built-in (internal) magnetic field of 6–12 T. Using direct imaging of the molecular configurations we observed fusion and braiding of eν-anyons under photo-excitation and revealed the possibility of using charge sensing for their control. Our investigations show that InP/GaInP2 anyon-molecule QDs, which have intrinsic transformations of localized eν-anyons compatible with TQC operations and capable of being probed by charge sensing, are very promising for the realization of TQC.

## Introduction

Methods of topological quantum computation (TQC) are based on topological quantum field theories, which show that two-dimensional (2D) electron systems having magnetic-flux-vortexes can be modeled by a Hamiltonian whose eigenstates correspond to the states of quantum error correcting code1. This also means that they are described by an effective Lagrangian containing a non-Abelian (NA) topological term2. According to the theory, the quantum computation process in these systems will be protected from environmental distortion at the physical level and will be intrinsically fault-tolerant. Topological quantum-gate (TQG) processing represents a physical realization of mathematical operations involving NA representations of a braid group Bn generated by a set of n localized particles (anyons) having zero-energy excitations, which are Majorana zero modes (MZMs), each representing a qubit state. A MZM qubit involves magnetic-flux-vortexes having opposite directions, representing a particle-antiparticle pair. The operation of the gate includes moving particles around each other (braiding) and pairing (fusion) together with the generation of the particles and the control of the resulting state1,3. Up to now, all efforts to realize TQC have been focused on finding an appropriate system having anyons with MZMs4,5,6,7,8,9. Although it is not clear how the gate operation will be physically realized using the specific systems investigated4,10, it should generally involve some local potential perturbations for particle trapping and moving, i.e., “tweezing” 11,12,13. This tweezing procedure, changing the TQG qubit state, is the equivalent of the resonant electro-magnetic pulses that induce Rabi oscillations of a two-level qubit state in the conventional schemes of quantum computing (QC)14. However, while the MZM-qubit and tweezing are supposed to be topologically protected, two-level qubits and Rabi oscillations are not, and they require adding redundant qubits to permit error-correcting code processing. This processing should provide an extremely low error-probability (fidelity) threshold to make QC operations fault-tolerant15. Such a threshold has been demonstrated only for a single qubit in the major QC platforms developed so far, which are Josephson tunnel junctions16 and electron-spin qubits in the different solid state environment or ion traps configuration17,18,19,20,21, but for two or more qubits the fidelity is poorer22. Thus, the realization of TQC is of high demand.

TQC relies on the existence of anions, which are composite particles consisting of electrons and a few magnetic-flux-quanta vortexes23. The anyon has fractional charge and its wave-function can have an arbitrary phase after interchange. Multi-dimensional NA representations of Bn are formed by a set of n “coupled” anyons having a strongly degenerate ground state24. The wavefunction in this case is a vector depending on the position and the quantum number of each particle; particle exchange then gives a matrix, i.e. a NA transformation of this vector, which is topologically protected. Theoretically, NA anyons can be formed with half-vortexes in p-wave superconductors, which have MZM in the core 25,26,27.

The richest and most investigated anyon system is a 2D-electron semiconductor heterostructure in a perpendicular magnetic field in which anyon states are formed at fractional Landau level (LL) fillings and provide a dissipationless skipping-type edge current in conjunction with a transverse conductance plateau known as the fractional quantum Hall effect (FQHE)28. More than 50 FQHE anyon states were observed in the lowest Landau level (LLL)29 and about 14 states in the second Landau level (SLL)30. NA anyons have been suggested and probed for the state ν = 5/2 having half-filled p-wave state at SLL31,32,33. For the hypothetical SLL NA stateν = 12/5 the braids that yield a universal set of quantum gates were found3. FQHE anyons are probed by the edge currents involving hundreds of mobile particles, which make difficult direct TQC.

Another type of structures suggested for TQC are hybrid heterostructures, in which half-vortexes are expected to form using the proximity effect between a conventional s-wave superconductor and various topological structures/materials, such as semiconductor nanowires having strong a spin–orbit coupling5, chains of ferromagnetic atoms6, topological insulators7, anomalous Hall insulator–superconductor8, and an insulating 2D quantum magnet9. The key difference between these systems and the FQHE one is that a single vortex is formed by up to a few thousands of electrons forming normal core at the center of a superconducting current loop7, and because of this “massive” structure, tweezing seems to be hard to realize in these systems. The last two systems use edge currents and have the same difficulties as FQHE system.

Currently FQHE anyons are considered theoretically to be quasi-particle excitations of the host ground-state liquid 34 and the quasi-particle picture is the basis for topological quantum field theories describing specific NA states and their use for building qubit and gate operations3,8. We have shown recently, however, that the FQHE anyon can exist as a single localized particles, which do not involve many-body interactions for their formation nor require quasi-particle concept for their description. Such fractionally charge particles were proposed to explain magneto-photoluminescence (magneto-PL) measurements of quasi-2D InP/GaInP2 single-electron islands (quantum dots) having Wigner–Seitz radius rs ~ 4 35. In this particle, which we called a magneto-electron (eν), a corresponding number of magnetic-quantum-flux vortexes (k = 1/ν) are self-generated.

Here we report the observation of molecular structures of eνs using magneto-PL measurements of QDs having about six-electrons and rs ~ 2 and show that these eν-anyon molecules (eν-AM) represent a novel system, which can be used for the realization of TQGs. The measurements include the imaging of the emission area of individual PL lines, together with quantum mechanical calculations and analysis of their electronic structure in a magnetic field. Using these we demonstrate the self-formation of eν-AMs and report observation of the eν-AMs having ν ~ 3/5–1/4, corresponding to a built-in magnetic field 6–12 T. We also observe a transformation of the eν arrangement in a 3 e2/7-AM under photo-excitation demonstrating fusion and braiding of the anyons. The observed transformation reveals a significant redistribution of the fractional charge within the dot, which suggests the use of single electron transistor charge sensing to control the TQC operations.

## Experimental details and data

### InP/GaInP2 quantum dots

We have studied nine single self-assembled InP/GaInP2 QDs36 having a lateral size of D ~ 120 nm and five-to-six electrons using near-field scanning optical microscope (NSOM) at temperature 10 K, external magnetic field range Be = 0–10 T and spatial resolution up to 25 nm (see “Methods”).

The dots were selected using NSOM measurements of a set of about 200 dots in variety of InP/GaInP2 QDs structures in previous studies36,37,38. Specifically, nine dots investigated were obtained in the sample X3141 (see Ref.36) and all of them have “anomalous” PL spectral features revealed in Refs.36,37,39 in comparison with the “normal” dots typically having more electrons and smaller size (see below). Two of such dots, also discussed here, were measured in the same sample.

All dots (except one normal) are denoted by Dkx, where the letter D indicates the initial identification of the shell filling (see the dot D5 in Ref.36, which is labeled D4e dot here), k enumerates the dots and x = a, b, c … indicates the NSOM run. One normal dot is denoted F1e. The dots are located within an area ~ 2 × 2 µm2 for the same run/scan (see below) and within an area ~ 2 × 2 mm2 for different runs.

For anomalous dots the parameters measured include the AM type, the main peak energy E0, the s-p splitting ΔEsp, the size D, the Wigner–Seitz radius rs, the built-in magnetic field Bbi, ν and the eν configuration. The AM type, E0 and ΔEsp, were measured from intensity distribution, position and energy splitting of the PL spectral components. The D values for most of the dots were measured/estimated directly from scanning experiments and used for calculation of rs = DN-0.5/(2aB*), where aB* ~ 8 nm35. Bbi and ν was estimated from a complex analysis of the whole data set including the dependence of the PL spectra on Be, the NSOM maps and theoretical calculations/analysis (see below). The eν configuration was suggested from analysis of the whole data set. These parameters apply to the photo-excited state (PS), which for all dots (except D1e) have N* = N + 1 = 6 electrons, where N is number of electrons in initial (IS) state. For D1e N* = 7.

Table 1 summarizes the dot parameters measured. The parameters, important for further discussion, are rs, which changes from 2.2 to 2.6 making up 20% variations and Bbi, which changes from 6 to 12 T making up two times (200%) variations.

### Theoretical description

The analysis of the experimental data was done using a phenomenological description within a framework of a general theory based on exact quantum mechanical calculations involving Fock–Darwin (FD), Hartree–Fock (HF) and configuration interaction (CI) approaches, developed for a few 2D electrons confined in circular potential in a perpendicular magnetic field about 20 years ago 40,41,42.

The many-body HF and CI methods were used to calculate dependence of the energy structure together with the shell energy and electron distribution on magnetic field for specific QD studied taking into account their shape, and, thus, a circular symmetry breaking effects (see “Methods”). These were used for analysis and comparison with corresponding experimental data.

Single-particle FD spectrum for circular symmetric dots43 was used in the form Ek,l(B) = Ek,l(Be + Bbi), where Ek,l(B) = Ek,lFD(Ω) = ħΩ(k + 1) − 1/2ħωcl, and Ω2 = ω02 + ωc2/4, where ω0 and ωc = e*B/m* are the quantum confinement and cyclotron frequencies, respectively, and k and l are the principal (k = 0, 1, ..) and azimuthal (l = -k, -k + 2, …, k-2) quantum numbers. We used an approximate expression for E0,0(B0) $$\approx$$ N*Vc derived in Ref.38 for N = 1, where νc=0.07ωc(m*Ω)−0.5 is a matrix element of the Coulomb interaction.

The FD spectrum was used to fit the experimental shifts of PL lines versus external magnetic field Be to estimate the electron charge e* and Bbi.

### Photo-luminescence spectra

#### Anomalous spectral features and intensity distributions

The spectra of anomalous dots presented in Fig. 1a have a set of sharp peaks, which are a main zero-energy e0-peak and about five e1- e5 peaks, related to anti-Stokes components (ASCs) having a splitting of ~ 0.5 meV. This is significantly different from the spectra of normal dots presented in the insert (see figure caption for their parameters), which reveal two or three ASC peaks, a few times larger splitting (3–5 meV), and order-of-magnitude larger peak broadening. While, in the normal dots, the zero-energy and ASCs peaks are related to the occupied degenerate s- p-, d-, … electron (e) shells (including a photo-excited e), in the anomalous ones they are related to single spin-polarized es and manifest Bbi, as we will show below. The energy splitting between the e0 and the e3 peaks ΔEsp, shown by a horizontal arrow, is related to the s-p-shell splitting, which, according to HF calculations, determines a quantum confinement energy ħω0* = 0.7ΔEsp. Besides the ASCs, Stokes components (SCs) denoted as nw0 (n = 1–3, ħω0 $$\approx$$ ΔEsp) are observed for some dots (D1d, D2b, D4e and D1e) and they are related to a center-of-mass vibration.

While the anomalous spectral features are nearly the same for all dots, the intensity distribution of ASCs shows significant variations. These variations reveal three types denoted by AM5,1, AMm and AM6,0. For the AM5,1 type (see two lowest spectra in Fig. 1a) the ASCs are an order of magnitude weaker than e0. For the AMm type (four middle spectra) the intensity of the ASCs increases a few times reaching a value of up to half of e0. For the AM6,0 type (see three upper spectra) the intensity of the ASC peaks increases further and becomes nearly the same as the intensity of the peak e0.

### Magnetic field dependence

In the magnetic fields Be = 0–10 T (see Fig. 1b) the AM5,1-type D1d dot reveals a very weak diamagnetic shifts of PL spectral lines (about 0.5 meV for 10 T) and a strong change in relative intensity versus Be resulting in the appearance of the type AMm at Be = 4–5 T and type AM6,0 at 8–10 T. The emergence of these two types is accompanied by the appearance of an additional peak e6 between peaks e1 and e2. For the AM6,0 type a SC at 2ω0 appears, similar to the dot D1e (see Fig. 1a).

For normal dots the shifts (not shown here) are an order of magnitude stronger (up to 1 meV/T) having anti-crossings and paramagnetic regions at fields < 2 T , while the intensity distributions do not show significant changes, as will be discussed elsewhere.

### Spatially-resolve spectra and imaging

The dots reveal anomalous size of the emission area (EA) of the single-e PL peaks, which is seen from spatially resolved data for dots D1e and D4e (see Fig. 1c,d and Fig. 2a–e). A variation of the intensity of e0-e5 peaks down to ~ 0.4% per nm is seen in spectra in Fig. 1c,d. This corresponds to an EA size dEA down to ~ 25 nm as it is seen in the PL intensity maps in Fig. 2a,b, giving dEA = 40 ± 20 nm.

The maps also show the location of the individual es, and the combined maps presented in Fig. 2c,d display molecular structures. The molecular structure is absent in the normal D0e and F1e dots, the EA of which corresponds to the dot size of ~ 90 and 110 nm, respectively. This is seen in the combined map of the area containing these four dots and presented in Fig. 2e. In the map the images of IS are presented for D1e and D4e dots. The map shows that these are neighboring dots arranged in a liner chain with separation ~ 200 nm.

The molecular structure of the QD D4e in PS in Fig. 2c is a compressed hexagon having a size ~ 120 × 70 nm and a bond length ~ 40 nm. The IS structure (see Fig. 2e) is a pentagon, in which e0 and e1 electrons are shifted by 10–20 nm along the (x,y) diagonal, as follows from ω0 EAs shape.

A complicated quasi-1D molecular structure oriented along the (x,y)-diagonal and having a size ~ 150 × 70 nm2 is observed for D1e (see Fig. 2d). In this structure e1, e2, e3 and e6 EAs overlap and e0,1 and e5,6 EAs have two locations, revealing degenerate 2e states, which gives a 7e PS of this dot. The overlapped peaks are located in the upper right corner of the map; the rest three peaks are located along the (x,y)-diagonal below at the distance ~ 70 nm and separated from each other by 40 nm.

The SCs maps correspond to a ~ 60 nm vertical up shift of e4 (ω0 map), a ~ 60 nm down shift e3 (2ω0 map) and a ~ 40 nm down shift of e5 along the (x,-y)-diagonal (2ω0 map). This results in e0-e4 and e3-e5 pairing and in the IS the molecular structure has nearly equilateral triangle arrangement of paired es at the vertexes, having sides ~ 80 nm and 90 nm and a bond length ~ 60 nm.

In the combined maps in Fig. 2c,d a dark–light color code of EAs corresponds to large-small separation of es from the photo-excited hole, which allows determine its location as shown in the maps. Thus, specific spectral shape observed is related to a specific distribution of such separations (see below).

### Analysis of the data

#### Classification of energies and states in magnetic field

A general theory of circular dots having N electrons in a magnetic field distinguishes four contributions to the total energy Etot, first discussed for 2e in Ref.44, which are the kinetic energies of localization induced by the confinement potential Econf(N) and by the magnetic field, i.e., cyclotron motion, Ecycl(N,B), respectively, and the Coulomb energies of electrons center-of-mass ECoul,c.m(N) and relative ECoul,rel(N,B) motion, respectively. The B-independent parts are Econf(N) = K(N)ħω0* and ECoul,c.m(N) rsEconf(N). The B-dependent cyclotron term Ecycl(N,B) = ħωcN*/2 is dominant for large fields when ωc > ω0. The B-dependent term ECoul,rel(N,B) is non-trivial and, in spite of its extremely small value ~ 0.01ECoul,c.m (see Fig. 3d below), provides ground-state transitions having discrete total angular momentum values Lz(B) = $${\sum }_{i}^{N}{l}_{zi}$$, where lzi is the angular momentum of single-es, known as magic numbers (LzMN). These appear because of the necessity of matching the spatial structure of the e wave-functions and the spatial symmetry of the e arrangement45. The e state of the specific LzMN is related to a filling factor ν = L0/LzMN, where L0 = N(N − 1)/2, linking it to FQHE46. The set of values B(LzMN) = Bν represent a “spectrum”, which for specific N value, is scaled as ~ 1/rs.

For a circular symmetry analog of the dot D1d, for which Etot(0) = 90 meV, Ecycl(B) = 4.2/B (meV/T) and ΔLz = 5 or 6 46, respectively, the molecular states are spin-polarized, ν < 1-states and are developed at B > 4 T, after maximum density droplet formation 47. In the range B = 5–14 T, embracing ν = 1/2–3/10, the dominant configuration is (5,1)46.

### Molecular configurations and PL spectra intensity distributions

In the QD D1d the mixing of the configurations takes place, due to non-circular shape, as is seen in electron density distributions, calculated using CI, in Fig. 3a. Nine density distributions shown for ν from 1 to 3/14 occupying B range 2.7–16 T reveal the onset of molecular structure formation near ν = 5/7 (B ~ 5 T) similar to the circular symmetry. A “pure” (6,0) and (5,1) configurations appear for ν = 5/7 and ¼ and for ν = 2/7 and 3/14, respectively. For other values of ν two mixed configurations are appeared. In the one (see ν = 1/2) a central e maximum is shifted ~ 10–20 nm down from the center, and in the other (see ν = 1/3) a pentagon with adjacent single-e maxima at the left is formed. Note, that at ν = 1/4 the calculated configuration is a compressed hexagon, exactly the same as observed experimentally in D1e, which implies Bbi ~ 12 T.

Molecular configurations in Fig. 3a can be assigned to specific AM types spectra observed in Fig. 1a,b. They, together with the maps in Fig. 2c,d, show that the appearance of different AM types spectra reflects a competition between (6,0) and (5,1) isomers for a non-circular dot shape and different Bbi or B. In these the former has nearly equal AM6,0 type ASC intensity, due to the nearly equal separation of electrons from the hole, while the latter has one dominant peak AM5,1 type owing to its location closer to the hole and the AMm-type has intermediate, mixed configurations and intensity distributions.

### Magneto-electrons and spontaneous anyon molecule formation

The analysis of the calculated B-dispersion of HF and FD shell energies, Lz, Etot and their comparison with experimental data of the dot D1d shown in Fig. 3b together with the FD fit to these data in Fig. 3c and ECoul,rel and Etot, shown in Figs. 3d,e, respectively, reveal the formation of the fractionally charged eνs.

The Etot curve in Fig. 3b shows a nearly dispersionless Econf(N) + ECoul,c.m(N) contribution for B < 3 T, a linear increase from Ecycl(N,B) and B > 4 T, and weak oscillations from ECoul,rel(B) at Bν (see also plot of ECoul,rel(B) in Fig. 3d) over the entire B-range. The LzMN states of the dot D1d are clearly visible on Lz(B) curve as plateaus and weak kinks at Lz = 5, 9, 15, 21, 25 30, 35 and 45. The corresponding Bν-spectrum has a ν values set from 3 to 1/3 and is very close to that of circular dot46. The plateaus are also seen in the region ν > 1corresponding to the integer (ν = 3) and SLL fractional (ν = 5/3) QHE states, demonstrating its deep connection to the localized es states in QDs.

This connection is also revealed for FD states, which appears as a matching of the crossings between the (1,0) and (0,l) FD levels to LzMN transitions at ν = 3 to 5/3 and to 1, corresponding to a total spin transitions from Sz = 0 to 1 to 2. The HF energies, accounting for spin and px-py circular symmetry distortion splittings, are shown in the figure in corresponding regions and approximately match the FD energies, neglecting the splitting. Both HF and FD energies at large fields come out to the LLL line ħωc/2, with a nearly order-of-magnitude reduction of inter level splitting, which evolves to zero in the limit B $$\to \infty$$.

From Fig. 3b we can see, that the critical discrepancy between the calculated FD/HF energies and the measured ASC shifts is that the latter have negligible B-dispersion and an order of magnitude larger level splitting over the entire 0–10 T magnetic field range. This implies a reduction of ωc, as can be revealed from the FD fit, shown in Fig. 3c. The FD fit gives a general matching with experimental shifts and splittings of e0-e5 peaks for Bbi ~ 6 T and a three-fold reduction of ωc. The total internal field B = Bbi + Be thus has the range B = 6–16 T, which includes ν = ½ and ¼.

The matching, however, does not account for the appearance of the additional peak e6 and the ~ 0.2 meV variations of the experimental B-dispersion, taking place near B ~ 7, 10 and 14 T. These features correspond to the peculiarities of ECoul,rel(B) in Fig. 3d at ν = ½ , 1/3 , 2/7 and ¼ and, thus, can be interpreted as the experimental signatures of these states. The appearance of the e6 peak could indicate crossing of single-e lz levels near steps in LzMN.

Reduction of ωc can be interpreted as fractional charge and manifests a self-generation of the magnetic-quanta-flux vortexes by single es forming eνs 35. The vortex self-generation arises because of the reduction of Etot parts ECoul,c.m + Ecycl caused by the decrease in charge and the dissipationless (superconducting) motion of es occupying quantum confined states for rs > 2. Such a reduction of Etot for eνs compared to es is shown in Fig. 3e. It reveals a ~ 10 meV minimum at ν = 5/3 and a gradual decrease from 100 meV for ν~ 5/7 to 50 meV for ν ~ 1/4, which corresponds to increase of the energy drop from 10 to 80 meV.

### Size dependence of molecular structure

Since the vortexes have fixed radius ~ aB* 35, the resulting charge, i.e., ν, Bbi, and AM configuration, are determined mostly by rs. The size of the eν area dν increases by a fraction ~ 0.1 per vortex, as can be found from simple geometrical drawing (see the corresponding cartoons for ν = ½, 1/3 and ¼ eν areas in Fig. 3e). The dν value of ~ 40 nm comes from the calculations for ν = 1/4 is in agreement with the measured dEA in the dot D4e. Thus, this dot can be considered as 6e1/4 and 5e1/4 AM in PS and IS, respectively. For the dot D1e a pairing of EAs in IS indicates 3e2/7 AM, which corresponds to Bbi ~ 11 T. For Bbi ~ 6 T of the dot D1d the PS state is 2e3/5 AM.

In Fig. 3f we present expected schematic arrangements of eνs in these three AMs obtained by placing vortexes onto the maxima of their e distribution. The e3/5 and e2/7 are realizations of the complex eνs, in which ν = n/p and n > 1, consisting of n-es fused by a “coupling” vortex (see below). The figure accounts for different e distribution size, i.e., rs, and demonstrates matching of ν and rs. This justifies our analysis, which, being extended to the rest of the dots, results in the values ν and Bbi and the eν-AM configurations in the PS given in Table 1. The analysis shows factor-of-two increase of Bbi for a 20% increase of rs in agreement with the experiment (see Table 1).

### Magneto-electrons and topological quantum computations

While the AM configurations observed and analyzed in 6e states show symmetric, homogenous molecular structures compatible with distorted five- and sixfold circular symmetry and with the calculated electron density distributions, for 7e state they are not. The observed anomalous 1D composite configuration for the 7e state in the dot D1e (see Fig. 2d) can be assigned to a decomposition of a symmetrical 3e2/7 AM of a 6e state after an extra e is added. In this case, the additional e2/7 needed for symmetrical arrangement, cannot be generated, since it requires pair of es. Thus, we can suppose that the additional e creates e1/3, which results to a transition from 21 to 24 vortexes state. The corresponding configuration could then consists of a 3e1/3 AM and a single e4/15. The corresponding arrangements of the eνs are presented in Figs. 4a,b. In the figures the coupling vortexes are shown by dashed circle.

Comparison of the arrangements and particle displacements in the scheme of Fig. 4a,b reveals transformations and interchange of eνs, which correspond to elementary topological quantum computing operations (TQCOs). These involve first, two unfusions of (e0e4) and (e3e5) e2/7 anyon pairs resulting in four e1/3 single anyons; second, two braids of (e0,e4) and (e1,e3) pairs; and, third, two fusions, which are a fusion of a e1/3 (e3,e6) pair forming a e2/7 anyon and a fusion of a e2/7 (e1e2) and a (e3e6) pair forming e4/15 anyon.

The above configuration transformation involves spatially separated eν-anyons localized on the specific levels of QD confinement potential and a full description of such a transformation should specify a localization energy, i.e., include the energy coordinate. Along these lines we add an energy coordinate to x–y plots and present the energy/space diagrams in Fig. 4c,d overlaid on the confinement potential. The diagrams outline a few meV difference in the confinement energy of coupled eνs.

They show that the coupling vortex should be presented by a pair of half-flux vortexes synchronously generated in gapped eν-states. These vortexes should have twice size (not shown in the figure) and thus embrace one of the neighboring vortexes. While such a representation of complex eνs seems unusual in the framework of conventional theories of FQHE states, assuming degenerate states and zeros of many-electron wave-functions for the vortex description34, it is supported by our data and by the experimental observations of 2/5 anyons by Aharonov-Bhom interferometry in quantum Hall bars48, giving independent evidence of their existence.

TQCOs resulting from the PL process in the dot D1e, presented in the Fig. 4a,b, are described by a B7 braid group diagram, i.e., time-position world lines, shown in Fig. 4e. The braid of e4e0 and e4e0 e1/3-anyons should result in a $$\Pi$$/3 phase change 23, which according to the diagram, is encoded in the unfusion of the corresponding e2/7 anyon and accompanied by a local charge redistribution. This allows one to use local charge sensing to control TQCOs, which can be provided with high sensitivity by a single-e transistor (SET)49. Since control of phase is a key step in TQCOs and a technology and methodology of SET fabrication and measurements are well developed50, using SETs opens up a new possibility for physical realization of TQGs.

The design of the TQG based on InP/GaInP2 eν-AMs will include few SETs and nano-electrode gates having sizes ~ 50 nm adjacent to individual QDs formed by semiconductor manufacturing technology using e-beam lithography 51. The arrays of QDs (see Fig. 2e)] can be used to realize circuits. Since we demonstrated that eν-AM structure is sensitive to small variation of rs, using nano-gates changing local potentials near QD gives possibility of the fine tuning of the initial state.

In Fig. 4a,b we draw three SETs and a gate near D1e QD, which can be used to realize a prototype of the device for electrical testing of TQCOs shown in Fig. 4e. In the device RF SETs acting simultaneously and utilizing a single RF line by means of carrier multiplexing 52 will be used to detect "charge triangulation". Each of the devices is biased at the slope of SET transfer characteristic and employs charge cancellation technique53 to minimize effects of direct capacitive coupling to the pulsing gate. In this way signals obtained from SETs (Usi) under variation of rs or injection/removing of the electrons by the gate voltage (UG), i.e., Usi(UG)-functions, will be used for designing of topological quantum computing processing.

Finally, we should point out, that the QDs considered do not directly involve NA anyons and MZMs, as suggested in initial TQC proposals. However, our magneto-PL measurements and the preliminary analysis of N ~ 8 InP/GaInP2 dots having rs ~ 1.5 (similar to D04 in the insert in Fig. 1a), which will be published elsewhere, reveal Bbi ~ 2 T close to ν ~ 5/2. This can indicate e1/4 anyon supporting MZMs, similar to that discussed for corresponding SLL FQHE state. The possibility of forming of the corresponding state in an appropriate QD naturally follows from eν-Etot of in Fig. 3e for ν> 1. Moreover, for such dot types we also observed a quasi 1D shape of EA36, similar to that in Fig. 2b, which can be a signature of TQCOs. We already demonstrated Coulomb blockade control of N in these dots in the range 8–18 having ΔUG ~ 0.3 meV per electron and the possible formation of a spin-polarized state for N > 15 54.

## Conclusion

We have introduced a novel system, which can be used for the realization of TQG. These are eν-anyon molecules (eν-AM) naturally formed in self-organized many electron InP/GaInP2 QD structures. We demonstrated their promising TQC prospects using high-spatial-resolution magneto-PL measurements of five-to-seven electron QDs having rs ~ 2. Our data, together with the quantum mechanical calculations of their electronic structure in a magnetic field and analysis demonstrate a self-formation of eν AMs, having ν ~ 3/5–1/4, corresponding to a built-in magnetic field of 6–12 T. Our measurements of PL spectra and the mapping of the intensity of individual PL lines gave a direct imaging of molecular configurations and allow the observation of a transformation of a eν arrangement in a e2/7-AM under photo-excitation involving fusion and braiding of the anyons. The observed transformation reveals a significant redistribution of the fractional charge within the dot, which suggests the use of single electron transistor charge sensing to control TQC operations. Our investigations show that InP/GaInP2 AM QDs having intrinsic anyon localization at zero external magnetic field combined with charge sensing control of anyon’s states open up novel directions for the realization of TQC.

## Methods

### Magneto-PL measurements and data processing

Spatially-resolved magneto-PL spectra were measured using NSOM operating at 10 K and magnetic fields of up to 10 T and using optical fiber probes having an aperture size of 50–300 nm in a collection-illumination mode. The spectra were excited by the 514.5 nm Ar-laser line and measured using a CCD (multi-channel) detector together with a 280 mm focal length monochromator. The excitation power measured before the fiber coupler was ~ 5 µW, which provided a power density of ~ 0.5 W/cm2. The spectral resolution of the system is 0.2–0.4 meV.

For the Lorenz contour deconvolution of PL spectra, we used a multi-peak fitting procedure from Origin 8.0 graphic software.

The spatially-resolved PL intensity at the selected wavelengths (image) was generated using the spectra taken in a square grid having a mesh of 50 nm. We plotted the experimental data using a contour plot option of Origin 8.0 and the division of the intensity data into 20 levels. The size of the emission area was estimated as a size of a NSOM image at the emission intensity level of 0.9, i.e. two upper contour plot levels. NSOM image of ASC are related to PS Stokes components images correspond to a difference of electron positions in the IS and PS36.

### Configuration interaction and Hartree–Fock many-body calculations

The states of $$N$$ confined electrons (N ~ 6 and rs ~ 2.5) in a strong perpendicular magnetic field $${B}_{z}=B$$ are modeled theoretically by configuration interaction (CI), using a method described in more detail in Ref. 36. We assume a quasi-2D effective-mass approximation with Hamiltonian (in atomic units, $$4\pi {\epsilon }_{0}={m}_{e}=\hslash =\left|e\right|=1$$)

$$H=\sum_{i=1}^{N}\left\{\frac{1}{2{m}^{*}}{\left[{\mathbf{p}}_{i}-e\mathbf{A}\left({\mathbf{r}}_{i}\right)\right]}^{2}+{V}_{\mathrm{ext}}\left({\mathbf{r}}_{i}\right)+{g}^{*}{{\mu }_{B}B}_{z}{S}_{zi}\right\}+\sum_{i>j}^{N}\frac{{e}^{2}}{\kappa \left|{\mathbf{r}}_{i}-{\mathbf{r}}_{j}\right|} ,$$

where the magnetic vector potential is given by $$\mathbf{A}\left(x,y\right)=\left({B}_{z}/2\right)\left(-y,x,0\right)$$, and the effective mass $${m}^{*}=0.077$$ and dielectric constant $$\kappa =12.61$$ are taken to be the values for bulk InP. The effective 2D confining potential $${V}_{\mathrm{ext}}\left(\mathbf{r}\right)$$ in the $$x$$-$$y$$ plane is assumed to be harmonic near the center, with frequency parameter $$\hslash {\omega }_{0}$$, and to develop a “hard wall” near the physical boundary of the dot on the substrate. The confining potential also has a slight angular deformation, indicated by the potential contours in Fig. 3a, which was chosen to be typical of the dots synthesized experimentally. The spin Zeeman term $${g}^{*}{{\mu }_{B}B}_{z}{S}_{zi}$$ is usually not included explicitly. The first step in the CI calculation is to calculate a single-particle basis set $$|i\rangle$$ using spin-polarized Hartree–Fock (HF),

$$\left[\frac{1}{2{m}^{*}}{\left({\varvec{p}}-e\mathbf{A}\right)}^{2}+{V}_{\mathrm{ext}}+{V}_{\mathrm{HF}}\right]\left|i\right.\rangle ={\epsilon }_{i}\left|i\right.\rangle .$$

The lowest $$N$$ HF orbitals $$|a\rangle$$ (where $$a=1, \dots , N$$) are occupied, with one electron per orbital assuming complete spin polarization. The occupied orbitals contribute to the self-consistent HF potential $${V}_{\mathrm{HF}}$$, which is defined as

$$\left.\langle j\right|{V}_{\mathrm{HF}}\left|i\right.\rangle =\frac{{e}^{2}}{\kappa }\sum_{a}^{\mathrm{occ}}\left(\langle ja\left|{r}_{12}^{-1}\right|ia\rangle -\langle ja\left|{r}_{12}^{-1}\right|ai\rangle \right) ,$$

where $$|i\rangle$$ and $$|j\rangle$$ are general states (occupied or unoccupied). The eigenvalues $${\epsilon }_{a}$$ of the occupied orbitals, at least for zero to moderate $${B}_{z}$$, may be expected to give approximations to the removal energies of electrons from the system, according to Koopman’s theorem from atomic and molecular physics54. In the calculated electron density distributions the dot size (Dc) was estimated as the size of the area containing 96% of electron density.

## References

1. 1.

Kitaev, A. Y. Fault-tolerant quantum computation by anyons. Ann. Phys. 303, 2–30 (2003).

2. 2.

Freedman, M. H. P/NP, and the quantum field computer. PNAS 95, 98–101 (1998).

3. 3.

Bonesteel, N. E., Hormozi, L., Zikos, G. & Simon, S. H. Braid topologies for quantum computation. Phys. Rev. Lett. 95, 140503–140504 (2005).

4. 4.

Sarma, S. D., Freedman, M. & Nayak, C. Topologically protected qubits from a possible non-abelian fractional quantum Hall state. Phys. Rev. Lett. 94, 166802–166804 (2005).

5. 5.

Mourik, V. et al. Signatures of Majorana fermions in hybrid superconductor-semiconductor nanowire devices. Science 336, 1003–1007 (2012).

6. 6.

Nadj-Perge, S. et al. Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor. Science 346, 602–6007 (2012).

7. 7.

Xu, J.-P. et al. Experimental detection of a Majorana mode in the core of a magnetic vortex inside. Phys. Rev. Lett. 114, 017001–017005 (2015).

8. 8.

He, Q. L. et al. Chiral Majorana fermion modes in a quantum anomalous Hall insulator–superconductor structure. Science 357, 294–299 (2017).

9. 9.

Kasahara, Y. et al. Majorana quantization and half-integer thermal quantum Hall effect in a Kitaev spin liquid. Nature 559, 227–231 (2008).

10. 10.

Sarma, S. D., Freedman, M. & Nayak, C. Majorana zero modes and topological quantum computation. npj Quantum Inf. 1, 15001–15013 (2015).

11. 11.

Wan, X., Hu, Z.-X., Rezayi, E. H. & Yang, K. Fractional quantum Hall effect at ν ¼ 5=2: Ground states, non-Abelian quasiholes, and edge modes in a microscopic model. Phys. Rev. B 77, 165316–165415 (2008).

12. 12.

Prodan, E. & Haldane, F. D. M. Mapping the braiding properties of the Moore-Read state. Phys. Rev. B 80, 115121–115214 (2009).

13. 13.

Macaluso, E., Compari, T., Mazza, L. & Carusotto, I. Fusion channels of non-abelian anyons from angular-momentum and density-profile measurements. Phys. Rev. Lett. 123, 266801–266807 (2019).

14. 14.

Ladd, T. D. et al. Quantum computers. Nature 464, 45–53 (2010).

15. 15.

Shor, P. W. Scheme for reducing decoherence in quantum computing memory. Phys. Rev. A 52, R2493–R2496 (1995).

16. 16.

Barends R.et al,. Superconducting quantum circuits at the surface code threshold for fault tolerance. Nature 508, 500–503 (2014).

17. 17.

Harty, T. P. et al. High-fidelity preparation, gates, memory, and readout of a trapped-ion quantum bit. Phys. Rev. Lett. 113, 220501–220505 (2014).

18. 18.

Veldhorst, M. et al. An addressable quantum dot qubit with fault-tolerant control-fidelity. Nat. Nanotechnol. 9, 981–985 (2014).

19. 19.

Watson, F. et al. A programmable two-qubit quantum processor in silicon. Nature 555, 633–637 (2015).

20. 20.

He, Y. et al. A two-qubit gate between phosphorus donor electrons in silicon. Nature 571, 371–375 (2019).

21. 21.

Rong, X. et al. Experimental fault-tolerant universal quantum gates with solid-state spins under ambient conditions. Nat. Commun. 6, 8748–8757 (2015).

22. 22.

Wright, K. et al. Benchmarking an 11-qubit quantum computer. Nat. Commun. 10, 5464–5466 (2019).

23. 23.

Wilczek, F. Magnetic flux, angular momentum, and statistics. Phys. Rev. Lett. 48, 1144–1146 (1982).

24. 24.

Moore, G. & Read, N. Nonabelions in the fractional quantum Hall effect. Nucl. Phys. B 360, 362–396 (1991).

25. 25.

Volovik, G. E. Fermion zero modes on vortices in chiral superconductors. JETP Lett. 70, 609–614 (1999).

26. 26.

Read, N. & Green, D. Paired states of fermions in two dimensions with breaking of parity and time-reversal symmetries and the fractional quantum Hall effect. Phys. Rev. B 61, 10267–10297 (2000).

27. 27.

Ivanov, D. A. Non-abelian statistics of half-quantum vortices in p-wave superconductors. Phys. Rev. Lett. 86, 268–271 (2001).

28. 28.

Tsui, D. C., Stormer, H. L. & Gossard, A. C. Two-dimensional magneto-transport in the extreme quantum limit. Phys. Rev. Lett. 48, 1559–1562 (1982).

29. 29.

Xia, J. S. et al. Electron correlation in the second Landau level: A competition between many nearly degenerate quantum phases. Phys. Rev. Lett. 93, 176809–176814 (2004).

30. 30.

Choi, H. C., Kang, W., DasSarma, S., Pfeiffer, L. N. & West, K. W. Activation gaps of fractional quantum Hall effect in the second Landau level. Phys. Rev. B 77, 081301–4(R) (2008).

31. 31.

Dolev, M., Heiblum, M., Umansky, V., Stern, A. & Mahalu, D. Observation of a quarter of an electron charge at the v = 5/2 quantum Hall state. Nature 452, 829–834 (2008).

32. 32.

Willett, R. L., Nayak, C., Shtenge, K., Pfeiffer, L. N. & West, K. W. Magnetic-field-tuned Aharonov-Bohm oscillations and evidence for non-abelian anyons at v =5/2. Phys. Rev. Lett. 111, 186401–186405 (2013).

33. 33.

Banerjee, M. et al. Observation of half-integer thermal Hall conductance. Nature 559, 205–210 (2018).

34. 34.

Laughlin, R. B. Anomalous quantum Hall effect: An incompressible quantum fluid with fractionally charged excitations. Phys. Rev. Lett. 50, 1395–1398 (1983).

35. 35.

Mintairov, A. M. et al. Fractional charge states in the magneto-photoluminescence spectra of a single-electron quantum dot. Nanomaterials 11, 493–514 (2021).

36. 36.

Mintairov, A. M. et al. Control of Wigner localization and electron cavity effects in near-field emission spectra of In(Ga)P/GaInP quantum-dot structures. Phys. Rev. B 97, 195443–195449 (2018).

37. 37.

Mintairov, A. M., Merz, J. L. & Blundell, S. Molecular states of electrons: emission of single molecules in self-organized InP/GaInP quantum dots. in Fingerprints in the Optical and Transport Properties of Quantum Dots, InTech ISBN 978-953-51-0648. 126–152 (2012).

38. 38.

Mintairov, A. M., Kapaldo, J., Merz, J. L., Vlasov, A. S. & Blundell, S. A. Wigner molecules and charged excitons in near-field magnetophotoluminescence spectra of self-organized InP/GaInP2 quantum dots. Phys. Rev. B 95, 11544–11610 (2017).

39. 39.

Mintairov, A. M. et al. Nanoindentation and near-field spectroscopy of single semiconductor quantum dots. Phys. Rev. B 69, 155306–155312 (2004).

40. 40.

Chacraborty, T. Physics of artificial atoms: Quantum dots in a magnetic field. Comments Condens. Matter Phys. 16, 35–68 (1992).

41. 41.

Maksym, P. A., Imamura, H., Mallon, G. P. & Aoki, H. Molecular aspects of electron correlation in quantum dots. J. Phys. Condens. Matter 12, R299–R334 (2000).

42. 42.

Reimann, S. M. & Manninen, M. Electronic structure of quantum dots. Rev. Mod. Phys. 74, 1287–1343 (2002).

43. 43.

Jacak, L., Hawrylak, P. & Wojs, A. Quantum Dots 176 (Springer, 1998).

44. 44.

Phannkuche, D., Gernhardts, R. R., Maksym, P. A. & Gudmundson, V. Theory of quantum dot helium. Physica B 189, 6–15 (1993).

45. 45.

Ruan, W. Y., Liu, Y. Y., Bao, C. G. & Zhang, Z. Q. Origin of magic angular momenta in few-electron quantum dots. Phys. Rev. B 51, 7942–7945 (1995).

46. 46.

Yannouleas, C. & Landman, U. Unified microscopic approach to the interplay of pinned-Wigner-solid and liquid behavior of the lowest Landau-level states in the neighborhood of ν=1/3. Phys. Rev. B 84, 165327–165417 (2011).

47. 47.

Reimann, S. M., Koskinen, M., Manninen, M. & Mottelson, B. R. Quantum dots in magnetic fields: Phase diagram and broken symmetry at the maximum-density-droplet edge. Phys. Rev. Lett. 83, 3270–3273 (1999).

48. 48.

Venkatachalam, V., Yacoby, A., Pfeiffer, L. & West, K. Local charge of the ν= 5/2 fractional quantum Hall state. Nature 469, 185–188 (2011).

49. 49.

Aassime, A., Johansson, G., Wendin, G., Schoelkopf, R. J. & Delsing, P. Radio-frequency single-electron transistor as readout device for qubits: Charge sensitivity and backaction. Phys. Rev. Lett. 86, 3376–3379 (2001).

50. 50.

Crippa, A. et al. Electrical spin driving by g-matrix modulation in spin-orbit qubits. Phys Rev. Lett. 120, 137702–137705 (2018).

51. 51.

Stevenson, T. R., Pellerano, F. A., Stahle, C. M., Aidala, K. & Schoelkopf, R. J. Multiplexing of radio-frequency single-electron transistors. Appl. Phys. Lett. 80, 3012–3014 (2002).

52. 52.

Yong, T., Orlov, A. O., Snider, G. L. & Patrick, F. J. Radio frequency operation of clocked quantum-dot cellular automata latch. Appl. Phys. Lett. 95, 193109–193113 (2009).

53. 53.

Mintairov, A. M. et al. Piezo-electric fields and state-filling photo-luminescence in natural InP/GaInP2 Wigner molecule structures. Appl. Phys. Lett. 118, 121101 (2021).

54. 54.

Lindgren, I. & Morrison, J. Atomic Many-Body Theory 2nd edn. (Springer, 1986).

## Acknowledgements

A.M.M., D.V.L. and A.S.V acknowledge the financial support of the Russian Science Foundation (project #19-19-00246) of magneto-PL NSOM experiments. G.S., A.O. and A.M.M. acknowledges the financial support of the National Science Foundation Grant (#DMR-1904610) of development TQG proposal.

## Author information

Authors

### Contributions

A.M.M., D.V.L and A.S.V. conceived the magneto-PL NSOM experiment. G.L.S. and A.O.O. suggest and describe the proposal for TQG using SET. The CI and HF calculations were carried out by S.A.B., using codes developed by him. The manuscript was written by A.M.M. with the input from all authors and with the critical contributions from S.A.B., A.O.O., and G.L.S.

### Corresponding author

Correspondence to Alexander M. Mintairov.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

### Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Mintairov, A.M., Lebedev, D.V., Vlasov, A.S. et al. Nano-photoluminescence of natural anyon molecules and topological quantum computation. Sci Rep 11, 21440 (2021). https://doi.org/10.1038/s41598-021-00859-6

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s41598-021-00859-6