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.

# Engineering topological states in atom-based semiconductor quantum dots

## Abstract

The realization of controllable fermionic quantum systems via quantum simulation is instrumental for exploring many of the most intriguing effects in condensed-matter physics1,2,3. Semiconductor quantum dots are particularly promising for quantum simulation as they can be engineered to achieve strong quantum correlations. However, although simulation of the Fermi–Hubbard model4 and Nagaoka ferromagnetism5 have been reported before, the simplest one-dimensional model of strongly correlated topological matter, the many-body Su–Schrieffer–Heeger (SSH) model6,7,8,9,10,11, has so far remained elusive—mostly owing to the challenge of precisely engineering long-range interactions between electrons to reproduce the chosen Hamiltonian. Here we show that for precision-placed atoms in silicon with strong Coulomb confinement, we can engineer a minimum of six all-epitaxial in-plane gates to tune the energy levels across a linear array of ten quantum dots to realize both the trivial and the topological phases of the many-body SSH model. The strong on-site energies (about 25 millielectronvolts) and the ability to engineer gates with subnanometre precision in a unique staggered design allow us to tune the ratio between intercell and intracell electron transport to observe clear signatures of a topological phase with two conductance peaks at quarter-filling, compared with the ten conductance peaks of the trivial phase. The demonstration of the SSH model in a fermionic system isomorphic to qubits showcases our highly controllable quantum system and its usefulness for future simulations of strongly interacting electrons.

## Main

Superconductivity, magnetism12, low-dimensional electron transport13, topological phases14 and other exotic phases of matter arise owing to the presence of strongly interacting particles within crystals15. However, the complexity of simulating such large quantum systems becomes intractable using classical computing methods16. A promising solution is to build a physical system at the same scale so that we can simulate these interacting fermionic systems17,18 directly, known as analogue quantum simulation19,20. The Su–Schrieffer–Heeger (SSH) model is the prototypical example of topological matter that describes a single electron hopping along a one-dimensional dimerized lattice with staggered tunnel couplings, v and w, as shown in Fig. 1a21. The SSH model has been experimentally simulated in physical systems of varying dimensions from Rydberg atoms (about 10 μm) to mechanical systems (about 10 mm) (Table 1). The coupling strengths of the various simulators lie in the nanoelectronvolt to microelectronvolt range, limiting their ability to reach the fully coherent regime. Importantly, these systems can be readily solved classically as they do not simulate many-body interactions. Only recently has the interacting many-body SSH model been observed using Rydberg atoms with an effective infinite on-site interaction (hardcore bosons)10. The ability to control the interaction strength, however, is critical for investigating fermionic systems22,23.

Semiconductor quantum dots are an emerging platform for the quantum simulation of strongly correlated electron systems4,5,20. Interacting electrons confined to quantum dots have been described by the Hubbard model16 involving Coulombic interactions that describe the energy required to add electrons to the same (on-site, U) or neighbouring (intersite, V) quantum dot24. Here intersite hopping is governed by the tunnel coupling, t, between quantum dots, and each dot can be tuned using electrostatic gates to raise or lower their energy levels, $${\epsilon }$$ (ref. 4). Phosphorus donors in silicon in particular have been proposed as promising candidates for simulators as they are nanoscale in size with very strong on-site energies (U ≈ 25 meV) and can be engineered to have strong intersite (V ≈ 5 meV) and hopping (t ≈ 5 meV) energies, while operating with a low thermal energy of kBT ≈ 0.02 meV (ref. 25), where kB is the Boltzmann constant and T is the temperature, reaching a range of U/t ≈ 1−100 while remaining in the low-temperature limit with t/kBT > 10 (ref. 23). The ability to reach the low-temperature, strongly interacting regime allows for a number of coveted quantum phases, such as superconductivity26 and antiferromagnetism27, to be simulated28. Despite the promise of semiconductor simulators, significant challenges have remained to simulate full quantum systems. These relate to the ability to precisely engineer, and tune, both the large on-site interaction energies and tunnel couplings to allow for the formation of a well defined coherent state across the system. In particular, for 10 quantum dots, we require precision control across 110 different experimental parameters related to U, V, t and $${\epsilon }$$.

In this paper, we utilize the atomic-precision placement accuracy of the scanning tunnelling microscope (STM) to engineer quantum dots with large on-site energies (U ≈ 25 meV) and uniform size to realize a homogeneous linear array for reliable simulation accuracy. If the quantum dots are too big, the capacitive coupling between individual dots becomes too large to independently control them. Conversely, if the they are too small, then a small change in the number of phosphorus donors within the quantum dot can substantially change the on-site energy, leading to randomness in the array. Importantly, our subnanometre-precision capability allows us to change the values of v and w with millielectronvolt resolution so that we can reliably enter both topologically trivial and topologically non-trivial regimes. Finally, a substantial challenge for gate-defined quantum dot architectures is that they require electrostatic gates to create the quantum dot potential and control the tunnel couplings requiring a minimum of about two gates per quantum dot29,30,31,32. With donor-based dots, we do not require these additional confinement gates and require only six electrostatic gates to control a ten-quantum dot array, thereby avoiding unnecessary cross-talk between gates. To ensure the creation of a well defined quantum state across the array, we designed an iterative maximum-current-alignment procedure to align the quantum dot energy levels within approximately 0.5 meV. The quantum state formed is then measured using bias spectroscopy via the planar source and drain leads. Having determined the necessary conditions to form the desired state, we simulate the one-dimensional topological phases associated with the interacting SSH model21.

The SSH model is one of the simplest known instances of topological quantum systems. The eigenenergies (Fig. 1b) of the SSH model give rise to two distinct phases dependent on the ratio of tunnel couplings with a topological phase transition at v = w. For v > w, there is a topologically trivial phase, where the lattice acts as a bulk insulator with the electron delocalized across the array and an energy gap between the upper and lower bulk states. For v < w, there exists a topologically non-trivial phase—a symmetry-protected topological phase—which gives rise to two zero-energy edge states where the electron is localized at the edge sites of the lattice33.

In previously measured instances of the SSH model, particle–hole symmetry was conserved owing to a lack of intersite electron–electron interactions21. However, quantum dots in semiconductors are affected by the intersite Coulomb interaction, Vi,j which is the change in energy of quantum dot j owing to the addition of an electron on quantum dot i. These long-range electron–electron interactions break the particle–hole symmetry, leading to non-degenerate electron and hole states4. As a consequence, it is important to control both the electron filling of the array and changes in the electrostatic environment to ensure correct simulation results. This requires us not only to have independent control of each quantum dot potential but also to alter the energy levels in unison—a tremendous technological challenge for such a small array with strong tunnel couplings and on-site energies (Supplementary Section II). As a result, we consider the full Hamiltonian of the extended (spinless) Hubbard model for a linear array of N quantum dots, given by

$${H}_{U}=\mathop{\sum }\limits_{i=1}^{N}{{\epsilon }}_{i}{n}_{i}+\mathop{\sum }\limits_{i}^{N}{U}_{i}{n}_{i}({n}_{i}-1)+\mathop{\sum }\limits_{i}^{N-1}{t}_{i,i+1}({c}_{i}^{\dagger }{c}_{i+1}+{\rm{h}}{\rm{.c}}.)+\mathop{\sum }\limits_{i,j}^{N}{V}_{i,j}{n}_{i}{n}_{j},$$
(1)

where $${{\epsilon }}_{i}$$ is the energy levels of the ith dot of the array, ni is the dot occupation operator, ti,i+1 is the tunnel coupling between nearest-neighbour ith and ith + 1 dots, Ui is the on-site Coulomb interaction term, Vi,j is the intersite Coulomb interaction terms between the ith and jth dots, $${c}_{i}^{\dagger }$$ (ci) are the creation (annihilation) operator of a particle on site i and h.c. indicates Hermitian conjugate.

In Fig. 1c, d, we show the 1,024 calculated multi-electron energy levels based on the two experimental arrays fabricated in the trivial and topological phases respectively, ordered by the number of electrons in the quantum dot array at quarter-filling. For ten quantum dots, in which each dot can host two electrons, quarter-filling corresponds to where there are five electrons shared across the dots. Here, if the array is in the trivial phase, the electrons will simply arrange themselves into each dimer such that they are spread evenly across the array. In contrast, if the array is in the topological configuration at quarter-filling, four electrons arrange themselves in the middle of the array (similar to the trivial phase). However, the fifth electron cannot occupy a dimer. We focus on quarter-filling of the array as the interacting topological states involve this fifth electron becoming localized to each end of the array. In the trivial phase (v > w), there is a singly degenerate ground state (box in Fig. 1c) delocalized across the entire array with a large energy gap of about 3.85 meV separating the five-electron state from the four- and six-electron states. The topological phase (v < w) exhibits a nearly four-fold degenerate (about 0.1 meV) ground state (box in Fig. 1d) involving four, five (two-fold degenerate) and six electrons. The ground state corresponds to four electrons on the inner eight quantum dots of the array with either zero, one or two additional electrons localized at the two edge quantum dots.

The experimental realization of two such arrays is shown in Fig. 1e, f, where the quantum dots are tunnel-coupled to source and drain leads to perform bias spectroscopy through the array. The regions outlined in orange and blue are the quantum dots, and the two regions at the start and end of the array outlined in white are the source and drain leads. The tunnel couplings ti,i+1 are engineered via the interdot separation, di,i+1 and follow an inverse exponential dependence, $${t}_{i,i+1}\propto \exp (-2{d}_{i,i+1}/3)$$ (ref. 34). By staggering the quantum dot array, we have ensured that non-nearest-neighbour tunnelling is exponentially suppressed with an estimated ti,i+2/ti,i+1 ≈ 0.01, ensuring that electron transport occurs in series through the array, while maximizing the differential lever arms to the dots (Supplementary Section I). The SSH model requires that the tunnel couplings have alternating strengths to observe the different topological phases while being simultaneously large enough to allow for a measurable transport current for bias spectroscopy. The quantum dot size is critical as the confinement potential experienced by the last outer electron, and hence the wavefunction overlap to the neighbouring quantum dots, depends on the number of donors comprising the quantum dot. We therefore fabricate quantum dots with an area of about 25 nm2 (about 25 phosphorus donors per site35). The nanoscale size of the quantum dots allows us to achieve large on-site energies of about 25 meV, but, importantly, where a small change in the size of the quantum dot does not significantly alter U, V, or t. The small separation of ≤10 nm also allows us to achieve large tunnel couplings t ≈ 1–10 meV for these quantum dot sizes35. Device I (Fig. 1e) is designed to be in the trivial phase with average staggered quantum dot distances, dv = 7.7 ± 0.1 nm and dw = 10.1 ± 0.2 nm, corresponding to $$\langle v/w\rangle =2.08$$. Device II (Fig. 1f) is designed to be in the topological phase where we now engineer an average staggered donor distance, dv = 9.6 ± 0.4 nm and dw = 7.8 ± 0.6 nm, corresponding to $$\langle v/w\rangle =0.265$$. These values highlight the subnanoscale accuracy that we can engineer devices with using STM lithography so that we can change $$\langle v/w\rangle$$ between 0.265 and 2.08 (ref. 36).

Figure 2a shows an STM image of the full device I. Here the outlined lighter regions show the lithographic hydrogen mask with six capacitively coupled control gates (G1 to G6), crucial to independently control the energy levels of the quantum dots (device II gate structure is nominally identical). Owing to the unique geometry of the device, the total lever arms of all gates linked together to each quantum dot are engineered to be consistent with a variation of less than 2.5%. This small variation means that we can also raise the global energy level of the whole quantum dot array for bias spectroscopy to measure the different phases of the SSH model. To align the energy levels of the quantum dots, we used a maximum-current-alignment scheme, in which the individual gates are tuned as outlined in Fig. 2b. This is achieved by initially setting the gate voltages at a conductance peak determined by sweeping G1, G2 and G3, against G4, G5 and G6, while measuring the current through the array. While positioned at this conductance peak, each gate is then individually swept around a set value, whereas all other gates were kept constant, as illustrated in Fig. 2c. After sweeping all six gates about their set voltage, the largest current peak is found and the corresponding gate is updated to the voltage at the centre of the current peak (G5 in the first iteration shown in Fig. 2d). All gates are then swept again, repeating this process, updating a single gate at a time as shown in Fig. 2d. Figure 2e shows the maximum current measured on each gate sweep per iteration for a constant source–drain bias, VSD. When the maximum current plateaus, VSD is reduced further and the entire process is repeated to increase the alignment accuracy.

Once the energy levels of the quantum dots are aligned, we perform a stability diagram measurement by shifting the energy levels of all quantum dots to investigate the electron occupation of the array. The stability diagram allows us to determine the electron occupation of the array as a function of the energy levels of the quantum dots. At zero source–drain bias (dashed white line in Fig. 2b), there is only enough energy to add a single electron to the array at a time. In this regime, we simulate the SSH model in the two devices. Figure 3 shows a comparison of the experimental and theoretical results of the two devices. Figure 3a shows the calculated normalized zero-bias conductance as a function of the ratio of tunnel couplings (v/w) of the array. The design of device I in the trivial phase is given by the dashed red line $$(\langle v/w\rangle =2.08)$$, whereas device II in the topological phase is given by the dashed blue line $$(\langle v/w\rangle =0.265)$$.

Figure 3b shows the zero-bias conductance as a function of the combined voltages on all the gates obtained from the trivial phase (device I) with the theoretical calculation shown in red. There are ten conductance peaks corresponding to a change in the total number of electrons on the array (Fig. 3c). We control the electron filling of the array by adjusting the gate voltages to tune the electron number from m to m + 10 (half-filling). At a quarter-filling (m + 5), there is a gap in the energy spectrum (7.9 meV, approximately twice the gap of 3.85 meV in Fig. 1c) corresponding to the single ground state of the SSH model for the trivial phase (Fig. 1c). From the estimate of Vi,j from electrostatic modelling, and fitting the magnitude of the tunnel coupling (Methods), we model the array to obtain the width of the different electron number regions, Sk (width of the m + k → m + k + 1 region). Figure 3d shows the width of the experimentally measured stability regions, obtained by determining the energy between adjacent conductance peaks, compared with the theoretical calculations based on electrostatic modelling with a tunnel coupling ratio $$\langle v/w\rangle =2.08$$. We find excellent agreement between the experimental and theoretical values. Small discrepancies ($$\lesssim$$1 meV) are most likely due to small misalignments of the quantum dot energy levels, which gives rise to on-site disorder, causing small shifts in the conductance peaks such that the peak structure is no longer symmetric around zero (Supplementary Section III).

We now look at the topological phase of the SSH model (device II) shown in Fig. 3e. The blue line in Fig. 3e represents the theoretical fit to the experiment, with a tunnel coupling ratio $$\langle v/w\rangle =0.265$$. We show a similar voltage range scan as for device I but here we observe only two sets of closely spaced peaks at zero gate voltage and at 85.5 mV corresponding to the average on-site energies across the array, $$\langle U\rangle$$ = 22.0 ± 3.2 meV. The conductance peaks from the states away from quarter-filling are not visible as they are now delocalized within the bulk of the array with a low probability of existing at the edge quantum dots. As a result, tunnelling between these bulk-like states and the source and drain leads is significantly suppressed. In the topological phase, the quarter-filling gap almost disappears completely with a sharp transition from the m + 4 to m + 6 states given by only two conductance peaks separated by about 0.2 meV, shown in Fig. 3f. These electron states correspond to where there are no electrons (m + 4 electrons) on the edge quantum dots to where both dots are occupied (m + 6 electrons). Importantly at exactly quarter-filling (m + 5 electrons), there is a non-zero probability that either four (P ≈ 0.05), five (P ≈ 0.30 or P ≈ 0.60 owing to the two-fold degeneracy) or six (P ≈ 0.05) electrons exist on the array at the same time owing to the nearly four-fold-degenerate ground state, as illustrated in Fig. 3g. This remarkable observation, that at zero gate bias there is a superposition of the number of electrons on the edge quantum dots, is a result of the near-zero energy of the topological states of the array and is a distinctive property of the many-body SSH model. As these topological states are localized at the edge quantum dots, the current flowing through the array corresponds to an electron moving from one side of the array to the other without occupying any of the inner quantum dots. This unique property is a direct consequence of the topology embedded within the SSH model as confirmed by the double conductance peak in Fig. 3f.

In conclusion, we have observed clear signatures of the topological states of the interacting SSH model in semiconductor quantum dots. To achieve this, we have precision engineered two devices with subnanometre resolution consisting of a linear array of ten donor-based quantum dots, with staggered nearest-neighbour tunnel couplings. The minimal gate design of our epitaxial devices allows for both individual and global alignment of all the energy levels of the quantum dots such that bias spectroscopy can probe the topological phases of the array. Importantly, we confirm the existence of the one-dimensional topological phase at quarter-filling $$(\langle v/w\rangle =0.265)$$ with the observation of two overlapping peaks in the zero-bias conductance corresponding to the near four-fold degeneracy of the many-body SSH model. For the trivial phase $$(\langle v/w\rangle =2.08)$$, we observe ten zero-bias conductance peaks corresponding to delocalized states across the entire array with an energy gap around quarter-filling owing to the large interdot coupling of about 3.5 meV. The low-gate-density design and low noise37 exhibited in these all-epitaxial devices offers a promising platform for simulating strongly interacting electron systems for quantum chemistry applications19. Future work will focus on extending the size of the quantum dot arrays, incorporating charge sensors and extending the simulations to engineered two-dimensional lattices38.

## Methods

### Device fabrication

The devices were fabricated on a 1–10 Ω cm p-type natural silicon substrate. The substrate was prepared via a series of high-temperature anneals, up to 1,100 °C, followed by a controlled cooling to 330 °C, resulting in a (2 × 1) reconstructed surface. The substrate was then terminated with atomic hydrogen, which can be selectively desorbed using the STM tip, leaving a hydrogen lithographic mask representing the device. After STM lithography, the substrate was exposed to a phosphine (PH3) precursor gas, in which phosphorus was absorbed and incorporated, at 350 °C, into the exposed areas. The device was then encapsulated with 40 nm of natural silicon using molecular beam epitaxy at a rate of about 0.125 nm min−1. A more detailed description of STM hydrogen lithography and the further device processing to electrically contact the device can be found in refs. 42,43,44.

### Experimental measurement setup

The electrical measurements were performed at millkelvin temperatures inside a dilution refrigerator. The two devices were measured in different fridges. Device I was measured in a fridge with a base temperature of about 10 mK, and device II was measured in a fridge with a base temperature of about 100 mK. Extended Data Fig. 1 shows a schematic of the electrical connections to the device. To perform the measurements, voltage sources were connected to the source and gates, to control their chemical potentials, and the transport current was measured through the drain. The drain current was amplified and converted to a voltage signal using a FEMTO variable-gain low-noise current amplifier DLPCA-200 with a low-pass 10-Hz filter. The filter signal was then digitized using a National Instrument data acquisition box (NIDAQ). The voltage sources used for the two devices were different. For device I, the source voltage was generated by the auxiliary of a Stanford Research Systems SR830 DSP lock-in amplifier and 1/50 room-temperature resistive-voltage divider, whereas the gate voltages were generated by two NIDAQs. Device II was measured on a different dilution refrigerator and had a different experimental setup. The source voltage was generated by a NIDAQ and a 1/50 room-temperature resistive-voltage divider. Gates 1, 2 and 3 voltages were also generated by the NIDAQ, whereas gates 4, 5 and 6 voltages were each generated by a Yokogawa 7651 programmable d.c. source. The different voltage sources used to measure the different devices should not affect the measurement results in a fundamental way.

### Calculation of the parameters of the extended Hubbard model

The quantum dot array can be described by the extended Hubbard model with the Hamiltonian given in equation (1). To theoretically solve the extended Hubbard model and calculate the parameters of the array, the data were fit using an open-source python package QmeQ45. The quantum dot Hamiltonian is described by the single-particle states $${\sum }_{i=1}^{N}{{\epsilon }}_{i}{n}_{i}+{\sum }_{i}^{N-1}{t}_{i,i+1}({c}_{i}^{\dagger }{c}_{i+1}+{\rm{h}}{\rm{.c}}.)$$ and the Coulomb matrix elements, $${\sum }_{i,j}^{N}{V}_{i,j}{n}_{i}{n}_{j}.$$ The Hamiltonian is constructed in the Fock basis45,46, for example, $$|0010101110\rangle$$, where 1 indicates an electron at that quantum dot and 0 indicates the quantum dot is unoccupied, and diagonalized exactly to obtain the exact many-body eigenstates $$|a\rangle$$, $$H=\sum {E}_{a}|a\rangle \langle a|$$. The transport current through the quantum dot array is then calculated numerically using the Pauli master equation. The array is assumed to be weakly coupled to the source and drain leads at a temperature, T, with a density of states that follows a Fermi distribution, f(E).

Analogous to the non-interacting energy-level diagram shown in Fig. 1b, we also calculated the energy-level diagram of the interacting system. As with Fig. 3a, we assume that the intersite Coulomb terms follow a d−1.5 dependence47, and that the energy levels of the quantum dots are tuned via $${{\epsilon }}_{i}=-{\sum }_{j}^{N}{V}_{i,j}$$. Extended Data Fig. 2 shows the energy spectrum for an array of ten sites for the m + 4, m + 5 and m + 6 electron states at quarter-filling, showing the more complex excited-state spectrum of the many-body states of the array.

For v/w > 1, at quarter-filling (m + 5 electrons), there is a large gap between the ground m + 5 and the lowest m + 4 and m + 6 electron states with a single ground state (labelled ‘bulk states’ analogous to Fig. 1b). By changing the tunnel couplings such that v/w < 1, the ground state becomes doubly degenerate with m + 5 electrons (labelled ‘edge states’ analogous to Fig. 1b). Here the quarter-filling gap is greatly reduced resulting in nearly degenerate states with differing electron numbers. The small energy gap at v/w < 1, is observed in the conductance trace in Fig. 3f and reflects the intersite Coulomb interactions present in the system.

In Fig. 3b and Fig. 3e, we show the theoretical calculation of the conductance through the quantum dot array fitted to the experimental data for the trivial and topological phases, respectively. In the theoretical calculation, we consider a spinless ten-dot array in the regime where there can be at most one electron on each of the quantum dots. Electron transport is restricted to sequential tunnelling through the array, as a result of the engineering of the device design, and only a single electron can tunnel through the quantum dots. No higher-order co-tunnelling events are allowed, as the source–drain bias is sufficiently small, making the choice of the Pauli master equation valid. All intersite Coulomb interaction terms, Vi,j, are included and accounted for by tuning the energy levels of the dots such that $${{\epsilon }}_{i}=-{\sum }_{j}^{N}{V}_{i,j}$$.

The tunnel coupling and intersite Coulomb interaction parameters of the extended Hubbard model are obtained by fitting to the measured conductance peaks in the trivial phase, and using the measured distances from the STM and the electrostatic modelling. The trivial phase was used for the fit, as in this phase all ten conductance peaks can be observed, whereas for the topological phase only two conductance peaks are observed. The conductance peaks in the trivial phase in Fig. 3b each correspond to transitions between different particle number ground states, which are separated in gate voltage by Sk (the voltage separation between the k and k + 1 particle number). These peak separations, Sk, are dominated by the tunnel coupling strengths, v and w, and the intersite Coulomb interaction strengths, Vi,j. Using the extracted values of Sk from the experimental results, we fit the overall magnitude of the tunnel coupling and intersite Coulomb interactions (while keeping the theoretically determined trends as a function of distance) to find the parameters shown in Supplementary Tables 4–6.

The tunnel coupling, t, is engineered via the interdot donor separation, d, and follows an exponential dependence, t = tMexp(−2d/3) eV, where the tunnel coupling magnitude tM = 0.1742 for a 1P–2P quantum dot system34. We anticipate a similar dependence of t as a function of interdot separation distance; however, tM is now used as a fitting parameter as it is known to depend on the crystallographic orientation of the quantum dots with v (about $$\langle 100\rangle$$) and w (about $$\langle 120\rangle$$) having different scaling factors owing to their different interdot axes. Supplementary Table 6 shows the distances and tunnel couplings for the two devices, with $$\langle v/w\rangle =2.08$$ (1.702, 2.460) for the trivial phase and $$\langle v/w\rangle =0.265$$ (0.146, 0.539) for the topological phase. The best fit to the zero-bias conductance peaks are achieved when the w tunnel couplings are approximately twice the v couplings for the same distance (tM,w ≈ 2tM,v). This angular dependence on the tunnel coupling between individual donors is a well known consequence of the silicon crystal lattice48,49,50; however, it has now been observed directly for quantum dots of larger size. This additional knowledge will assist in future experimental designs for finer engineering of the tunnel couplings.

The extracted tunnel-coupling parameters give excellent fits to the experimental data as seen in Fig. 3b, e, with Fig. 3d showing a comparison of the peak separations, Sk, for the experimental and theoretical data. Small variations between the experimental and theoretical data can be attributed to errors in the electrostatic modelling and small offsets in the alignments of the quantum dots.

The presence of intersite Coulomb interactions gives rise to an approximately 0.2 meV splitting observed in the topological states as seen in Fig. 3f. Extended Data Fig. 3 shows the calculated conductance for the topological phase for varying intersite Coulomb interaction strengths, Vi,j, with values given in Supplementary Tables 5, 6. With no intersite Coulomb interactions, a single conductance peak is observed (Extended Data Fig. 3a), whereas the inclusion of intersite Coulomb interactions gives rise to two conductance peaks, in which the separation between the peaks increases with increasing Vi,j. In addition to the presence of intersite Coulomb interactions, a small splitting (about 0.04 meV) in the topological states would arise owing to the non-zero v/w ratio and the finite length of the chain. Here the topological states, exponentially localized at opposite ends of the chain, have a finite overlap. This small splitting, as presented in Extended Data Fig. 3a, is, however, too small to be observed experimentally.

## Data availability

The data pertaining to this study are available from the corresponding author upon reasonable request.

## References

1. Bardeen, J., Cooper, L. N. & Schrieffer, J. R. Theory of superconductivity. Phys. Rev. 108, 1175–1204 (1957).

2. Osheroff, D. D., Richardson, R. C. & Lee, D. M. Evidence for a new phase of solid He3. Phys. Rev. Lett. 28, 885–888 (1972).

3. Klitzing, K. V., Dorda, G. & Pepper, M. New method for high-accuracy determination of the fine-structure constant based on quantized Hall resistance. Phys. Rev. Lett. 45, 494–497 (1980).

4. Hensgens, T. et al. Quantum simulation of a Fermi–Hubbard model using a semiconductor quantum dot array. Nature 548, 70–73 (2017).

5. Dehollain, J. P. et al. Nagaoka ferromagnetism observed in a quantum dot plaquette. Nature 579, 528–533 (2020).

6. Atala, M. et al. Direct measurement of the Zak phase in topological Bloch bands. Nat. Phys. 9, 795–800 (2013).

7. Meier, E. J., An, F. A. & Gadway, B. Observation of the topological soliton state in the Su–Schrieffer–Heeger model. Nat. Commun. 7, 13986 (2016).

8. St-Jean, P. et al. Lasing in topological edge states of a one-dimensional lattice. Nat. Photon. 11, 651–656 (2017).

9. Chaunsali, R., Kim, E., Thakkar, A., Kevrekidis, P. & Yang, J. Demonstrating an in situ topological band transition in cylindrical granular chains. Phys. Rev. Lett. 119, 024301 (2017).

10. de Léséleuc, S. et al. Observation of a symmetry-protected topological phase of interacting bosons with Rydberg atoms. Science 365, 775–780 (2019).

11. Cai, W. et al. Observation of topological magnon insulator states in a superconducting circuit. Phys. Rev. Lett. 123, 080501 (2019).

12. Nielsen, E. & Bhatt, R. N. Nanoscale ferromagnetism in nonmagnetic doped semiconductors. Phys. Rev. B 76, 161202 (2007).

13. Le, N. H., Fisher, A. J. & Ginossar, E. Extended Hubbard model for mesoscopic transport in donor arrays in silicon. Phys. Rev. B 96, 245406 (2017).

14. Le, N. H., Fisher, A. J., Curson, N. J. & Ginossar, E. Topological phases of a dimerized Fermi–Hubbard model for semiconductor nano-lattices. npj Quantum Inf. 6, 24 (2020).

15. Imada, M., Fujimori, A. & Tokura, Y. Metal–insulator transitions. Rev. Mod. Phys. 70, 1039–1263 (1998).

16. Qin, M., Schäfer, T., Andergassen, S., Corboz, P. & Gull, E. The Hubbard model: a computational perspective. Annu. Rev. Condens. Matter Phys. 13, 275–302 (2022).

17. Feynman, R. Simulating physics with computers. Int. J. Theor. Phys. 21, 46–488 (1982).

18. Cirac, J. I. & Zoller, P. Goals and opportunities in quantum simulation. Nat. Phys. 8, 264–266 (2012).

19. Argüello-Luengo, J., González-Tudela, A., Shi, T., Zoller, P. & Cirac, J. I. Analogue quantum chemistry simulation. Nature 574, 215–218 (2019).

20. Altman, E. et al. Quantum simulators: architectures and opportunities. PRX Quantum 2, 017003 (2021).

21. Su, W. P., Schrieffer, J. R. & Heeger, A. J. Solitons in polyacetylene. Phys. Rev. Lett. 42, 1698–1701 (1979).

22. Tasaki, H. The Hubbard model—an introduction and selected rigorous results. J. Phys Condens. Matter 10, 4353–4378 (1998).

23. Salfi, J. et al. Quantum simulation of the Hubbard model with dopant atoms in silicon. Nat. Commun. 7, 11342 (2016).

24. Wang, X., Yang, S. & Das Sarma, S. Quantum theory of the charge-stability diagram of semiconductor double-quantum-dot systems. Phys. Rev. B 84, 115301 (2011).

25. Dusko, A., Delgado, A., Saraiva, A. & Koiller, B. Adequacy of Si:P chains as Fermi–Hubbard simulators. npj Quantum Inf. 4, 1 (2018).

26. Maier, T., Jarrell, M., Pruschke, T. & Keller, J. d-wave superconductivity in the Hubbard model. Phys. Rev. Lett. 85, 1524–1527 (2000).

27. Hirsch, J. E. & Tang, S. Antiferromagnetism in the two-dimensional Hubbard model. Phys. Rev. Lett. 62, 591–594 (1989).

28. Arovas, D. P., Berg, E., Kivelson, S. & Raghu, S. The Hubbard model. Annu. Rev. Condens. Matter Phys. 13, 239–274 (2022).

29. Nichol, J. M. et al. High-fidelity entangling gate for double-quantum-dot spin qubits. npj Quantum Inf. 3, 3 (2017).

30. Mills, A. R. et al. Shuttling a single charge across a one-dimensional array of silicon quantum dots. Nat. Commun. 10, 1063 (2019).

31. Yang, C. H. et al. Operation of a silicon quantum processor unit cell above one kelvin. Nature 580, 350–354 (2020).

32. MacQuarrie, E. R. et al. Progress toward a capacitively mediated CNOT between two charge qubits in Si/SiGe. npj Quantum Inf. 6, 81 (2020).

33. Asbóth, J. K., Oroszlány, L. & Pályi, A. The Su–Schrieffer–Heeger (SSH) Model (Springer, 2016).

34. Wang, Y., Chen, C.-Y., Klimeck, G., Simmons, M. Y. & Rahman, R. All-electrical control of donor-bound electron spin qubits in silicon. Preprint at https://arxiv.org/abs/1703.05370 (2017).

35. Watson, T. F. et al. Transport in asymmetrically coupled donor-based silicon triple quantum dots. Nano Lett. 14, 1830–1835 (2014).

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

37. Kranz, L. et al. Exploiting a single-crystal environment to minimize the charge noise on qubits in silicon. Adv. Mater. 32, 2003361 (2020).

38. Wang, X. et al. Quantum simulation of an extended Fermi-Hubbard model using a 2D lattice of dopant-based quantum dots. Preprint at https://arxiv.org/abs/2110.08982 (2021).

39. Kim, E. et al. Quantum electrodynamics in a topological waveguide. Phys. Rev. 11, 011015 (2021).

40. Lin, S., Zhang, L., Tian, T., Duan, C.-K. & Du, J. Dynamic observation of topological soliton states in a programmable nanomechanical lattice. Nano Lett. 21, 1025–1031 (2021).

41. Kruk, S. et al. Edge states and topological phase transitions in chains of dielectric nanoparticles. Small 13, 1603190 (2017).

42. Rueß, F. J. et al. Realization of atomically controlled dopant devices in silicon. Small 3, 563–567 (2007).

43. Broome, M. A. et al. Two-electron spin correlations in precision placed donors in silicon. Nat. Commun. 9, 980 (2018).

44. Fricke, L. et al. Coherent control of a donor-molecule electron spin qubit in silicon. Nat. Commun. 12, 3323 (2021).

45. Kiršanskas, G., Pedersen, J. N., Karlström, O., Leijnse, M. & Wacker, A. QmeQ 1.0: an open-source python package for calculations of transport through quantum dot devices. Comput. Phys. Commun. 221, 317–342 (2017).

46. Lin, H. Q., Gubernatis, J. E., Gould, H. & Tobochnik, J. Exact diagonalization methods for quantum systems. Comput. Phys. 7, 400 (1993).

47. Hogg, M. Scalable Readout and Control of Phosphorus Donor Spin Qubits in Silicon. PhD thesis, UNSW Sydney Australia (2021).

48. Koiller, B., Hu, X. & Das Sarma, S. Exchange in silicon-based quantum computer architecture. Phys. Rev. Lett. 88, 027903 (2002).

49. Wellard, C. J. et al. Electron exchange coupling for single-donor solid-state spin qubits. Phys. Rev. B 68, 195209 (2003).

50. Joecker, B. et al. Full configuration interaction simulations of exchange-coupled donors in silicon using multi-valley effective mass theory. New J. Phys. 23, 073007 (2021).

## Acknowledgements

The research was supported by Silicon Quantum Computing Pty Ltd and the Australian Research Council Centre of Excellence for Quantum Computation and Communication Technology (project number CE170100012). M.Y.S. acknowledges an Australian Research Council Laureate Fellowship.

## Author information

Authors

### Contributions

S.K.G., Y.H. and M.Y.S. conceived the project. M.K., H.G., M.B.D., J.G.K. and Y.C. fabricated the devices. M.K. and S.K.G. performed the measurements, analysed the data and performed the theoretical calculations. The manuscript was written by M.K., S.K.G., J.G.K. and M.Y.S. with input from all authors. M.Y.S supervised the project.

### Corresponding author

Correspondence to M. Y. Simmons.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Peer review

### Peer review information

Nature thanks Ingmar Swart and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

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

## Extended data figures and tables

### Extended Data Fig. 1 Experimental measurement setup.

a, Schematic of the experimental set up for Device I, showing all electrical connections to the device. A STM image of the device I is shown. The gate voltages are controlled by a NIDAQ and three Yokogawas, and the source voltage by the NIDAQ. The drain current (green line) is amplified via a FEMTO low noise amplifier and acquired by the NIDAQ. b, Schematic of the experimental setup for Device II, showing all electrical connections to the device. A STM image of the 10 dot array is shown, taken before patterning of the control gates. The gate voltages are controlled by two NIDAQs, and the source voltage controlled by a SR830 DSP lock-in amplifier. The drain current (green line) is amplified via a FEMTO low noise amplifier and acquired by the NIDAQ.

### Extended Data Fig. 2 Energies of the many-body states in the interacting SSH model around quarter-filling.

The energies of m + 4, m + 5, and m + 6 electron states for an array of 10 sites as a function of the tunnel coupling ratio, v/w. The energies are calculated at quarter-filling (that is, the ground state is always m + 5 electrons). The ground state evolves from a singly degenerate state for v/w > 1 (bulk-like states) to a two-fold degenerate m + 5 electron state for v/w < 1 (edge-like states). The quarter-filling energy gap between the m + 4 and m + 6 electron ground states reduces as v/w < 1 resulting in nearly degenerate states, reflecting the almost zero-energy cost associated with adding an electron to the topological phase of the array, with a small energy gap due to the inter-site Coulomb interactions.

### Extended Data Fig. 3 Splitting of the topological phase due to the presence of inter-site Coulomb interactions.

A comparison of the theoretically calculated conductance traces (blue lines) with the experimental data observed (circles) as the inter-site Coulomb interaction strength is varied from, a, no Coulomb interactions, to b, 0.5 × Vi,j, c, Vi,j, and d, 2 × Vi,j. In the case of no Coulomb interactions only a single peak is observed in the theoretical conductance trace, which does not match the roughly 0.2 meV splitting observed in the experimental data. As the Coulomb interaction strength is increased two peaks evolve with the splitting between the peaks increasing, with the roughly 0.2 meV splitting matching the experimental data in c.

## Supplementary information

### Supplementary Information

This Supplementary Information file includes Supplementary Tables 1–6, Figs. 1–5 and References.

## Rights and permissions

Reprints and Permissions

Kiczynski, M., Gorman, S.K., Geng, H. et al. Engineering topological states in atom-based semiconductor quantum dots. Nature 606, 694–699 (2022). https://doi.org/10.1038/s41586-022-04706-0

• Accepted:

• Published:

• Issue Date:

• DOI: https://doi.org/10.1038/s41586-022-04706-0