Quantum transport and localization in 1d and 2d tight-binding lattices

Particle transport and localization phenomena in condensed-matter systems can be modeled using a tight-binding lattice Hamiltonian. The ideal experimental emulation of such a model utilizes simultaneous, high-fidelity control and readout of each lattice site in a highly coherent quantum system. Here, we experimentally study quantum transport in one-dimensional and two-dimensional tight-binding lattices, emulated by a fully controllable 3 × 3 array of superconducting qubits. We probe the propagation of entanglement throughout the lattice and extract the degree of localization in the Anderson and Wannier-Stark regimes in the presence of site-tunable disorder strengths and gradients. Our results are in quantitative agreement with numerical simulations and match theoretical predictions based on the tight-binding model. The demonstrated level of experimental control and accuracy in extracting the system observables of interest will enable the exploration of larger, interacting lattices where numerical simulations become intractable.


INTRODUCTION
Single-particle dynamics and transport in periodic solids is well described by the tight-binding model 1 . This model is widely used to calculate the electronic band structure of condensed matter systems and to predict their transport properties, such as the conductance 2,3 . In the presence of a periodic lattice potential, the wavefunction of a given quantum particle overlaps neighboring lattice sites, leading to extended Bloch wavefunctions 1 . In the absence of lattice disorder, the particle propagation is ballistic and is described by a continuous-time quantum random walk 4 . This is in contrast to classical diffusive transport, where the propagation is quadratically slower in time 5 .
The quantum nature of transport in lattices leads to the emergence of non-local correlations and entanglement between lattice sites. However, lattice inhomogeneity causes scattering and leads to quantum interference that tends to inhibit particle propagation, a signature of localization 6,7 . The wavefunction of a localized particle rapidly decays away from its initial position, effectively confining the particle to a small region of the lattice.
Here, we study Anderson localization and Wannier-Stark localization in one-dimensional (1d) and two-dimensional (2d) tight-binding lattices. In the presence of random disorder of the lattice site energies, a particle wavefunction becomes spatially localized, known as Anderson localization 6 . In the context of electrical transport, this phenomenon alters the transport properties, e.g., decreasing the conductance by increasing the degree of localization, ultimately leading to an insulating state. Alternatively, the particle can be localized in the presence of a potential gradient across the lattice, e.g., as created in the transport case by an external, static electric field, referred to as Wannier-Stark localization 7 . The potential gradient induces a positiondependent phase shift in the particle wavefunction and causes the particle to undergo periodic Bloch oscillations in a confined region 8 . While the dynamics of these oscillations are challenging to observe in naturally occurring solid-state materials 9 , they can be directly emulated and experimentally studied using engineered quantum systems.
Anderson localization has been experimentally realized in Bose-Einstein condensates 10,11 , atomic Fermi gases 12 , and photonic lattices 13,14 . Similarly, Bloch oscillations and Wannier-Stark localization have been observed in optical lattices 9,15 , semiconductor superlattices 16 , and photonic lattices 17,18 . However, these demonstrations were limited to varying degrees by a lack of control and simultaneous readout of each individual site. Therefore, these experiments could not fully explore the different localization regimes and spatial dimensions of the tight-binding model in a single experimental instantiation.
Superconducting quantum circuits are highly-controllable quantum systems that enable us to probe the properties of lattice models. By engineering the Hamiltonian of the superconducting quantum circuit, we can emulate the dynamics of the lattice Hamiltonian, and we can use single-and two-qubit gate operations for state preparation and tomographic state readout [19][20][21][22][23] . Site-selective qubit control, strong qubit-qubit interactions, and long coherence times relative to typical gate times combined with the capability of simultaneously applying gates and performing tomographic state readout on all sites make this a promising platform for studying models of information propagation, many-body entanglement, and quantum transport.
In this work, we use a 9-qubit superconducting circuit to emulate the dynamics of single-particle quantum transport and localization in 1d and 2d tight-binding lattices. We probe the entanglement formed in the lattice as a result of the particle propagation. We then realize Anderson localization by implementing random disorder of the on-site energies with tunable strength, and extract the dependence of the propagation mean free path on the disorder strength in a regime with no analytical solution. Additionally, we study Wannier-Stark localization in the presence of isotropic and anisotropic potential gradients, and find close agreement between our experiments and theoretical predictions.
The close agreement between our experimental results, simulations, and theoretical predictions results from high-fidelity, simultaneous qubit control and readout, an accurate calibration, and an accounting for the relevant decoherence mechanisms in this system. Although our experiments are performed on a small lattice that can still be simulated on a classical computer, they demonstrate a platform for exploring larger, interacting systems where numerical simulations become intractable.

Experimental setup
Our device consists of a 3 × 3 array of capacitively coupled superconducting transmon qubits 24 (Fig. 1b), where qubit excitations correspond to particles in the lattice model. In the single-particle regime, the effective system Hamiltonian in the rotating frame with reference frequency ω ref is described by the tight-binding Hamiltonian: whereσ y i (σ i ) are the two-level creation (annihilation) operators for each lattice site i. The first term in the Hamiltonian represents particle tunneling between neighboring sites with rate J i,j , realized here with an average strength of J/2π = 8.1 MHz at ω ref = 5.5 GHz. The second term represents the particle occupation energy of sites i with the corresponding site energies ϵ i = ω i − ω ref , where ω i is the fundamental transition frequency of the transmon at site i. By using flux-tunable transmons, the energy of each site can be individually set over the range 3.0 GHz ≲ ω i /2π ≲ 5.5 GHz. This tunability enables us to realize arbitrary energy landscapes and to isolate sub-lattices (e.g., one-dimensional chains) by detuning certain qubits from their neighbors so that they do not interact with the rest of the lattice. The transmons have an average anharmonicity of U/2π = − 244 MHz, which corresponds to the on-site interaction energy of two particles occupying the same site. The system operates in the J ≪ |U| limit such that b Optical image of the 3 × 3 superconducting qubit lattice used in the experiment. The capacitor pads of the transmons qubits are false colored. c Operation sequence of the experiments. Each qubit is initially detuned from ω ref and can be addressed individually, and simultaneously. They are brought on resonance to emulate evolution underĤ for time t. The system dynamics are frozen by once again detuning the qubits, and all nine qubits are measured simultaneously. d A 1d quantum random walk with a single excitation (particle) initialized at the edge of the lattice (site 7). e A 2d quantum random walk with a single excitation (particle) initialized at the corner of the lattice (site 7). f The spatial propagation of the particle is measured by the root of the second moment of position ( ffiffiffiffiffiffiffiffiffi ffi hM 2 i p ) in the 1d and 2d lattices. The excitation initialized at the edge of the 1d chain exhibits ballistic propagation in time and with propagation speed v g ¼ ffiffi ffi 3 p J (red line), whereas an excitation initialized in the corner of a 2d lattice exhibits an increased propagation speed v g ¼ ð1 þ ffiffiffiffiffiffiffi ffi 3=2 p ÞJ (blue line).
A.H. Karamlou et al. each lattice site can be occupied by at most a single particle, generally realizing the hard-core Bose-Hubbard model 25,26 . Our experiments feature simultaneous, site-resolved, single-shot dispersive qubit readout with an average qubit state assignment fidelity of 95%. In addition, we are able to control the energy of each site with an average precision 〈Δϵ〉 < (2π) 200 kHz (≈J/40). These features are crucial to accurately emulate the 1d and 2d tight-binding model with different lattice potential configurations, and allow us to achieve close agreement between numerical simulations and experiment.
Quantum random walk Quantum random walks (QRWs) are the quantum mechanical analog of classical random walks. The spatial propagation of the particle throughout the lattice, relative to its initial location, is quantified by the mean square displacement where p i is the probability of finding the particle on site i, and M i is the Manhattan (1-norm) distance between site i and the initialization site. Quantum properties such as single-particle superposition and interference result in a qualitative difference between classical and quantum random walks: a classical random walk propagates diffusively in time with ffiffiffiffiffiffiffiffiffi ffi , whereas a QRW exhibits ballistic propagation with a mean-square displacement ffiffiffiffiffiffiffiffiffi ffi hM 2 i p / t 27 . In order to experimentally observe particle propagation through QRWs, we first implement a tight-binding lattice with uniform site energies (ϵ i = 0). We compare the respective propagation speeds for QRWs in a seven-site 1d chain (Fig. 1d) and a 2d 3 × 3 lattice (Fig. 1e). We inject a particle at the end (corner) of the 1d (2d) lattice and observe its propagation by tracking the excitation numbers hn i i ¼ hσ y iσ i i on each lattice site versus evolution time t. In the 1d case, the particle traverses the lattice and eventually reflects off the opposite end of the chain 28 . All intermediate sites are coupled to their two nearest neighbors, with the end sites being coupled to only one. Hence, the QRW propagates in both directions: an excitation is reemitted in both the forward and reverse directions (interim sites), or the direction is reversed (end sites). Quantum interference between these multiple trajectories alters the particle wavepacket as it evolves in time. The resulting QRW pattern for a seven-site chain features a non-trivial evolution that agrees well with numerical simulation (Fig. 1d).
In the 2d QRW, the particle propagates along its diagonal symmetry axis, analogous to a QRW on a binary tree 5 . Here, the particle position is defined by the sum of the site populations at a given Manhattan distance M from the injection site. The quantum interference leading to the simple back-and-forth pattern as a result of reflection is exceptional and arises from the high-degree of symmetry present in a 3 × 3 lattice (Fig. 1e). Such a high-degree of symmetry is similarly observed in 1d for a five-site chain (see Supplementary Material).
To verify these experimental results, we perform Lindblad master equation simulations (see Supplementary Material) based on the Hamiltonian in Eq. (1). We find excellent agreement between experimental data and numerical simulations by using the measured coupling strengths J ij between neighboring lattice sites i, j, and taking into account qubit relaxation and decoherence.
Both 1d and 2d QRWs exhibit ballistic propagation, and the propagation speed in 2d is faster than in 1d due to the constructive interference between multiple propagation paths leading to each site in 2d (Fig. 1f). The average group velocity v g of transport depends on the dimensionality of the lattice and the starting position of the particle with respect to the lattice boundaries. For a 1d QRW, a particle prepared at one end of the lattice initially propagates with v g = J, because it interacts with only one other lattice site. At later times, prior to reflection from the other end, it reaches a steady state value of v g ¼ ffiffi ffi 3 p J, in agreement with theory for infinitely long chains 27 . In contrast, we observe that a particle prepared in the corner of a 3 × 3 lattice propagates with an initial group velocity v g ¼ ffiffi ffi 2 p J (due to the coupling to its two nearest neighbors) and a steady state average velocity of v g ¼ ð1 þ ffiffiffiffiffiffiffi ffi 3=2 p ÞJ (see Supplementary Material). During the QRW, certain sites become entangled as the particle propagates through the lattice, a phenomenon that underpins its intrinsic quantum character 28,29 . We observe the emergence of entanglement as a result of the system time evolution via nonlocal spatial correlations. We quantify the entanglement formed among the qubits at the same distance from the particle initialization site, and study how the QRW entangles this site with the rest of the lattice in a coherent manner. The amount of entanglement within a two-qubit sub-system, described by ρ i,j , can be quantified using the pairwise concurrence 30 : where {λ μ } (in decreasing order) are the eigenvalues of the Here,ρ ¼ ðσ y σ y Þρ Ã ðσ y σ y Þ is the spin-flipped density matrix obtained through complex conjugation and applying Pauli-Y operators (σ y ) to each qubit. Concurrence is a monotone entanglement metric for two-qubit states that takes values between zero (product state) and one (maximally entangled). To reconstruct the two-qubit density matrix ρ i,j (t) and calculate the pairwise concurrence C i;j ðtÞ, we perform two-qubit state tomography for various evolution times t of the QRW.
We measure the concurrence formed during a 2d QRW between site pairs at the same Manhattan distance from the particle initialization site (Fig. 2a). Lattice sites at the same Manhattan distance become partially entangled as the particle traverses them, as indicated by an increase in the concurrence. This approach faithfully describes the quantum correlations in the subsystems comprising of two sites M = 1, 3. However, the subsystem of sites at M = 2 contains three qubits, and solely considering pairwise concurrence values does not fully capture the entanglement within this sub-system.
Using the pairwise concurrence values between the sites at M = 2, we calculate the lower-bound of the distributed concurrence C min iðj;kÞ for each site in this sub-system (Fig. 2b). For a state consisting of three qubits (i, j, k), the concurrence C iðj;kÞ between site i and the system comprised of the two remaining sites (j, k) can be lower-bounded using the pairwise concurrences 31 : Using the extracted pairwise and distributed concurrence values, we compute the average of the lower-bound concurrence C for sub-systems consisting of qubits that are at the same distance from the initial particle position (Fig. 2c). C is equal to the pairwise concurrence for sub-systems containing two qubits (M = 1, 3), and an average of C min iðj;kÞ for the sub-system containing three qubits (M = 2). The emergence of entanglement reflects the particle trajectory during the QRW (see Supplemental Material) and is in agreement with the light-cone of the measured quantum information propagation using out-of-time-ordered correlators 26 within the Lieb-Robinson bounds 32 (see Supplementary Material).
While the initially prepared single-particle state is separable from the rest of the lattice, the QRW, mediated by the nearestneighbor interactions, creates and distributes entanglement to varying degrees throughout the lattice. For a system containing exactly a single particle, Eq. (3) becomes an equality and can be generalized to a system with more than three sites 31 . In our 3 × 3 lattice, the concurrence between the initialization site (site 7) and the rest of the system C 7;lattice can be extracted exactly using pairwise concurrence values (see Supplementary Material): We observe that the quantum state of the initialization site becomes fully entangled with the larger system at time t ≃ (2J) −1 (Fig. 2d). The coherent propagation leads to the fall and revival of the concurrence as the particle wavefunction evolves with time, with C 7;lattice ¼ 0 when site 7 is in a pure ground or excited state. This behavior is in agreement with the evolution of the von Neumann entropy Sðρ 7 Þ ¼ Àtrðρ 7 lnρ 7 Þ (Fig. 2e), where ρ 7 is the single-qubit density matrix of site 7. The von Neumann entropy is a measure for the entanglement of the lattice site with its environment, and takes a value of lnð2Þ for a maximally entangled state. The dynamical revival of both C 7;lattice and S(ρ 7 ) indicates that site 7 becomes predominantly entangled with the rest of the lattice, rather than with the uncontrolled environmental degrees of freedom related to decoherence. Time evolution under an interacting Hamiltonian creates entanglement among the nine sites in our system. We quantified the entanglement between different sub-systems of the lattice using concurrence and von Neumann entropy, however, the overall entanglement formed in the lattice is also of interest for studying the quantum properties of the system. Without the need for full quantum state tomography, we probe the global entanglement E gl in the many-body system, beyond the nonlocal correlations formed between sub-systems, by extracting the average purity of all N sites 33,34 where ρ j is the reduced density matrix describing site j. In Fig. 2f, we observe that the global entanglement reaches a maximum value when the particle is not concentrated on either corner of the lattice, but rather extends across all lattice sites. In our 2d lattice, E gl is upper-bounded by the global entanglement in the 9-qubit W-state W 9 j i 35 , where W N j i ¼ 1 ffiffi ffi N p P π πð0:::01Þ j i is the superposition of all single-particle state permutations. At later times, the maximum global entanglement in the system decreases as a result of system decoherence. Using this metric, we study the degree of entanglement in the lattice as a result of the particle evolution.

Anderson localization
We have so far explored quantum transport in lattices with uniform site energies. However, particle propagation is significantly impacted by the introduction of random disorder in the site energies, leading to Anderson localization 6 . In the presence of disorder in 1d and 2d lattices, the particle wavefunction becomes spatially localized, regardless of the disorder strength in the thermodynamic limit 36 . In our experiments, we emulate the Anderson localization regime by introducing disorder with strength δ, where site energies are randomly sampled from a uniform distribution ϵ i ∈[−δ/2, δ/2]. As a result, wavefunction scattering causes the particle to lose phase coherence on the length scale of the mean-free path l. While an analytical form for the mean free path can be derived in the limit of weak (δ ≪ J) and strong (δ ≫ J) disorder 37,38 , there are no known analytical expressions in the intermediate disorder regime 39 .
We examine the degree of localization by measuring the increase in average population hn s i of the initially prepared lattice site and the decrease in particle spread ffiffiffiffiffiffiffiffiffi ffi hM 2 i q for increasing disorder strengths δ, each averaged over 50 random lattice realizations. As the disorder strength increases, a higher steadystate population remains on the source site as a result of the tightbinding interaction (Fig. 3a, c), while the average particle spread away from the source decreases (Fig. 3b, d) for both 1d and 2d lattices. For a given disorder strength, the particle is more confined in a 1d lattice compared to a 2d lattice; in 2d, the probability of transport is relatively greater in the presence of represents the minimum concurrence between qubit i and the sub-system consisting of qubits j and k. c Average minimum concurrence formed between the sites along the diagonal symmetry axis of the lattice. We display the population on the sites at distance 0 and 4 as a reference (shaded regions). d Concurrence of the propagation source and the rest of the lattice. The source becomes maximally entangled with the remaining of the lattice at t ≃ (2J) −1 . e The von Neumann entropy of the propagation source (site 7), reaching a maximum value of ð2Þ at time t ≃ (2J) −1 . f Global entanglement E gl among the lattice sites. The upper-bound of E gl for a single-particle tight-binding model is set by the global entanglement of the 9-qubit W-state W 9 j i.
A.H. Karamlou et al. disorder as propagation occurs along multiple paths between two sites. As a result of localization caused by disorder, the emergence of multipartite entanglement in the tight-binding lattice is inhibited: in Fig. 3e, f, we report a decrease in the steady state value for the average global entanglement of the system, E gl , with increasing disorder strength, due to spatial localization of the particle wavefunction. With one excitation, the confinement effects for seven-site chain causes a relatively greater deviation from the maximum value (E gl of W 7 j i) compared to the 2d lattice. We quantify the localization of the wavefunction in a lattice with N sites using the participation ratio PR, defined as 40 If the particle wavefunction is completely delocalized, then PR = N, whereas PR = 1 for a wavefunction fully localized to a single site (Fig. 3g). In our experiments, PR for each lattice realization is calculated using the time-averaged population on each site after the tight-binding model dynamics reach steady state 40 . In order to reduce the impact of qubit relaxation on these measurements, we use the averaged site populations between 100 ns and 400 ns (5/J ≲ t ≲ 20/J) to calculate PR (see Supplementary Material for an example). We experimentally extract the participation ratio PR averaged over different random lattice realizations, at different disorder strengths for both a 1d and 2d lattice. As δ increases, the particle wavefunction becomes more spatially confined and PR decreases. For weak disorder strengths, we find the localization length to be larger than the lattice size (gray region in Fig. 3h, i), imposing a limitation on our experiments caused by boundary effects. In a finite 1d lattice, the participation ratio is related to the localization length ξ 1d through the expression PRðξ 1d Þ ¼ cothð1=ξ 1d Þ tanhðN=ξ 1d Þ (see Supplementary Material). The 1d localization length ξ 1d scales directly with the mean free path l as ξ 1d = l 36 . We experimentally extract the dependence of the mean free path on the disorder strength, taking the form l = a(J/δ) γ , in our seven-site chain by realizing 60 random lattice disorders for different disorder strengths and calculating PR (Fig. 3h). In the 1d case, we find γ = 1.0 ± 0.03 and a = 17.0 ± 1.05 (in units of the lattice spacing) through fitting.
In a 2d n × n lattice, the participation ratio is related to the localization length through the expression PRðξ 2d Þ ¼ coth 2 ð1=ξ 2d Þ tanh 2 ðn=ξ 2d Þ (see Supplementary Material), where the localization length takes the form ξ 2d ¼ l e π 2 kl with a latticedependent factor k 36 . For weak disorder, the 2d PR scales exponentially with the mean free path. Consequently, the scaling for small values of δ is difficult to observe in our small lattice, and hence the value of k cannot be obtained with high accuracy. We extract the 2d PR in our 3 × 3 lattice (Fig. 3i) over 180 random lattice realizations for different disorder strengths, and for large values of δ we find the parameters for the mean free path γ = 0.80 ± 0.02 and a = 12.98 ± 0.60 (in units of the lattice spacing). In a 2d lattice, we notice that the extracted mean free path exhibits a Fig. 3 Anderson localization. Population at the propagation source (a) and particle root mean-squared distance from the source (b) for different 1d lattice disorders. Each point is averaged over 50 disorder realizations. Population at the propagation source (c) and particle meansquared distance from the source (d) for different 2d lattice disorders. Each point is averaged over 50 disorder realizations. The global entanglement among the sites in the 1d lattice (e) and 2d lattice (f) measured at different disorder strengths averaged over 50 realizations. g The wavefunction for a completely delocalized state takes the value of PR = N, whereas PR = 1 for a wavefunction fully localized to a single site. h 1d Anderson localization participation ratio PR scaling with disorder strength, obtained using 60 random lattice realizations for each value of disorder strength δ. The blue shaded region indicates the simulated statistical uncertainty interval (±1 standard deviation) as a result of lattice randomization. i 2d Anderson localization PR scaling with disorder strength, obtained using 180 random lattice realization for each value of δ.
A.H. Karamlou et al. weaker dependence on the disorder strength, consistent with our earlier observations. Wannier-Stark localization A particle in the tight-binding lattice is localized also in the presence of a static electric field, which creates a potential gradient across lattice site energies 7 . We emulate the effect of an electric field by creating a gradient in the lattice site energies ϵ i . The resulting potential gradient causes the particle to become spatially confined due to the emergence of a band-gap in the lattice band structure 9 . This phenomenon is referred to as Wannier-Stark localization 7,15,41 . The effective field creates a linear gradient F ! along the lattice axis, where F ! ¼ ∇ϵ i , and causes the particle to undergo periodic Bloch oscillations. We demonstrate Bloch oscillations in a 1d chain for the field strength F/J = 1.5 (F ¼ j F ! j) observing a spatially periodic breathing motion with a probability up to 99% of the particle reviving at the initialization site (Fig. 4a). By repeating the experiment for different values of F, we observe that a particle initialized in the center of the lattice oscillates with a Bloch period T B = 2π/F (Fig. 4b) and a maximum particle spread of d max (Fig. 4c) from the initialization site, in agreement with theory 8 . The finite lattice size causes the boundaries to have a notable effect on the oscillation for small values of F (gray region of Fig. 4b, c), namely the decrease in the oscillation period and limiting the particle spread. We observe a similar periodic motion along the diagonal symmetry axis of the lattice in 2d (Fig. 4d) for an isotropic field (F = F x = F y ). The Bloch period of a particle initialized in the corner of a 2d lattice with this symmetric gradient exhibits the same scaling with F as in 1d. We find that the maximum particle spread in a 2d lattice scales as d max B ¼ maxð ffiffiffiffiffiffiffiffiffi ffi hM 2 i p Þ ¼ 3:01 ± 0:01 J=F, where we obtained the value 3.01 from a data fitting procedure. Wannier-Stark localization is less pronounced in 2d compared to 1d for the same potential gradient, similar to the trend we observe in the Anderson localization regime and in agreement with theoretical predictions 42 .
We also explore the effect of non-isotropic fields in a 2d lattice by independently controlling the potential gradients along the horizontal (F x ) and the vertical (F y ) axes, respectively. In Fig. 4g, we show the time-averaged population of each site during the first 500 ns ( ≈ 25/J) of evolution as a result of different field ratios r = F y /F x . In the presence of the same field along both lattice axes (r = 1), the particle is localized along the diagonal symmetry axis of the lattice (Fig. 4g top left). As r increases by decreasing F x while keeping F y fixed, the particle becomes less localized along the F x -direction, skewing the propagation direction.
In the absence of boundary effects, the relation between Bloch oscillation periods T Bx and T By along each lattice axis depends on the net field direction, namely T Bx ¼ r T By (Fig. 4h).
For rational values of r, the total 2d Bloch period T B is the least Fig. 4 Wannier-Stark localization. a Exemplary 1d Bloch oscillation with a particle initialized at the center of the lattice. The particle undergoes a spatially periodic breathing motion before re-converging at the initialization site. b, c 1d Bloch oscillation period T B and maximum particle spread d max B scaling with the potential gradient strength F. d Exemplary 2d Bloch oscillation with an isotropic potential gradient with a particle initialized at the corner of the lattice. The particle undergoes periodic Bloch oscillations along the diagonal symmetry axis of the lattice. e, f Isotropic 2d Bloch oscillation period T B and maximum particle spread d max B depending on F. g Time-averaged population hn i i on each site for different anisotropy ratio r = F y /F x . As r increases, the propagation gets skewed towards the horizontal axis. h, i T B and d max common multiple of the one-dimensional Bloch periods along each axis 42 . A similar relationship holds for the maximum particle spread in each direction: d max Bx ¼ r d max By (Fig. 4i). In our experiments, we investigate these relations by measuring T Bx (see Fig. 4h) and d max Bx (see Fig. 4i) for r = 1, 2, 3 (solid circles) while varying F y . The data match simulation (solid lines) of the 3 × 3 system very well. We additionally measured T By and d max By , which are nominally unchanged for r = 1, 2, 3, and present the values measured for isotropic potential gradients (open circles). This enables us to obtain experimental estimates for r by taking the ratio between the experimentally measured Bloch periods T Bx =T By (Fig. 4h, right panel) and the maximum spread d max Bx =d max By (Fig. 4i, right panel). The ratios for each r are estimated from highlighted points in Fig. 4h, i.

DISCUSSION
In this work, we have emulated quantum particle propagation and localization in 1d and 2d tight-binding lattices using a 3 × 3 array of superconducting qubits. We have measured the coherent dynamics of different entanglement metrics, such as concurrence and the von Neumann entropy, during quantum transport using simultaneous control and readout. We have further investigated localization in different regimes of the tight-binding model, with random disorder resulting in Anderson localization and with a potential gradient causing Wannier-Stark localization and Bloch oscillations. We have studied Anderson localization in a disorder regime that lacks an analytical form, and used our data to extract the dependence of the propagation mean free path in finite 1d and 2d lattices on the disorder strength. We have measured the degree of Wannier-Stark localization as a result of different potential gradients for both isotropic and anisotropic potentials, and have observed quantitative agreement in the properties of Bloch oscillations with theoretical predictions. We can control all of the on-site energies with high precision compared to the coupling rate between neighboring sites. Furthermore, we can manipulate the lattice (qubit) energy landscape on timescales much shorter than the particle exchange (qubit-qubit interaction) and system coherence times by two orders of magnitude. This ability to accurately initialize and manipulate our system Hamiltonian combined with high-fidelity tomographic readout enables us to study tight-binding lattices under various configurations. For superconducting qubits, models based on calibrated qubit parameters and decoherence 43,44 have been demonstrated to accurately predict the behavior of multiqubit processors 45,46 . These features form the basis for agreement between our experimental results and our simulations. The demonstrated degree of control in realizing different regimes of the tight-binding model in 1d and 2d serves as a blueprint for exploring larger and strongly interacting quantum lattices in the few-body 47 and the many-body localization regimes 40,48,49 or that contain a time-dependent Hamiltonian leading to interesting phenomena such as quantum scars 50 topological order 51 and the breakdown of the Magnus expansion 52 .

DATA AVAILABILITY
The raw data used for generating the plots in this paper are available upon request.