arxiv: 2604.06944 · v1 · submitted 2026-04-08 · ❄️ cond-mat.mes-hall · quant-ph

Recognition: unknown

Tensor-network simulation of quantum transport in many-quantum-dot systems

Maximilian Streitberger , Marko J. Ran\v{c}i\'c

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:58 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall quant-ph

keywords quantum transportquantum dotstensor networksopen quantum systemssteady-state currentsmany-body simulationnonequilibrium transportnumerical methods

0 comments

The pith

Tensor networks with a jump-counting estimator compute steady-state currents through arrays of up to fifty quantum dots.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a method to simulate electron transport in large arrays of interacting quantum dots by representing the system's state with tensor networks. Traditional solvers based on density matrices cannot handle the rapid growth in complexity beyond small numbers of dots. The new approach adds a jump-counting estimator that tallies tunneling events between leads and dots to obtain currents directly in the steady state. Benchmarks show the currents match those from a standard master-equation solver when both can be run. This reduction in memory and computation time makes systematic studies of size-dependent transport feasible for systems that were previously inaccessible.

Core claim

The central claim is that a tensor-network solver augmented by a jump-counting estimator for lead-induced tunneling events yields steady-state electron currents that agree quantitatively with established master-equation results in small systems while scaling to arrays of fifty quantum dots, thereby cutting memory requirements and runtime by orders of magnitude relative to density-matrix methods.

What carries the argument

The tensor-network compression of the open quantum system's density operator together with the jump-counting estimator that accumulates current contributions from individual tunneling events.

If this is right

Quantitative agreement holds with the master-equation solver across tested lead-dot and inter-dot coupling strengths in the regime where both methods apply.
Memory use and wall-clock time drop by orders of magnitude compared with density-matrix solvers.
Simulations become possible for interacting quantum-dot arrays far larger than those reachable by conventional transport codes.
Systematic exploration of how nonequilibrium currents depend on array length and interaction strength becomes practical.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same jump-counting idea could be combined with other tensor-network ansätze to treat two-dimensional or disordered dot lattices.
If the method remains accurate, it would let researchers map out the crossover from coherent to incoherent transport as dot number increases.
Extension to time-dependent driving or finite-bias transients would require only modest changes to the estimator.
The approach may help benchmark approximate analytic theories for large open quantum systems where exact references are unavailable.

Load-bearing premise

The tensor-network truncation and the jump-counting procedure continue to reproduce the exact steady-state currents without accumulating noticeable bias when the array size grows to fifty dots.

What would settle it

Compute the steady-state current for a fifteen-dot chain with moderate couplings using both the tensor-network method and an exact diagonalization of the Liouvillian on a sufficiently large classical computer, then check whether the two results differ by more than a few percent.

Figures

Figures reproduced from arXiv: 2604.06944 by Marko J. Ran\v{c}i\'c, Maximilian Streitberger.

**Figure 1.** Figure 1: Four-quantum-dot device schematic. Schematic of a linear array of four quantum dots tunnel-coupled in series and connected to source and drain reservoirs. 1 arXiv:2604.06944v1 [cond-mat.mes-hall] 8 Apr 2026 [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗

**Figure 2.** Figure 2: Current and method discrepancy as functions of lead–dot and inter-dot coupling. Top panels show the current as a function of the lead–dot coupling, Γ (left), and the inter-dot coupling, Ω (right), for single-quantum-dot (1-QD, blue) and four-quantum-dot (4-QD, green) systems. QmeQ results are shown as crosses, and TJM results as circles connected by solid lines. Bottom panels show the corresponding absolut… view at source ↗

**Figure 3.** Figure 3: Runtime and memory scaling of TJM and QmeQ with system size. (a): Wall-clock time as a function of the number of quantum dots for TJM and QmeQ. Although QmeQ is faster for small systems, its computational cost increases much more rapidly with system size. TJM becomes advantageous beyond the crossover regime and remains the more efficient approach for larger systems. (b): Memory requirement as a function of… view at source ↗

**Figure 4.** Figure 4: Source–Drain bias sweep. Here the left lead is acting as the source and the right lead as the drain. The potential difference corresponds to the bias Vbias = µL − µR. We can see that the current decreases systematically with increasing array length. number of trajectories, system size and maximal bond dimension. By contrast, the memory requirement of the Lindblad solver grows combinatorially with the numbe… view at source ↗

**Figure 5.** Figure 5: Current dynamics and runtime scaling beyond the direct benchmark regime. (a) Current dynamics in a 50-dot [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗

read the original abstract

Transport through correlated nanoscale systems underpins the operation of quantum-dot and molecular-scale devices, yet accurate simulations of large open quantum systems remain computationally challenging as system size increases. Tensor-network methods offer a promising route past this scaling barrier by efficiently compressing quantum states. Here we extend a tensor-based solver with a jump-counting estimator that enables direct computation of steady-state electron currents from lead-induced tunneling events. We benchmark the resulting currents against the state-of-the-art master-equation solver QmeQ across a range of lead-dot and inter-dot coupling parameters and find quantitative agreement in the tractable regime. Compared with classical approaches, TJM reduces memory requirements and wall-clock time by orders of magnitude, enabling simulations of interacting quantum-dot arrays far beyond the range accessible to density-matrix-based transport solvers and systematic studies of size-dependent nonequilibrium transport in larger arrays. Our approach allow us to model quantum transport in an array of up to fifty (50) quantum dots.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Tensor networks plus jump-counting reach 50 dots with good small-N benchmarks, but error scaling at large N is not shown.

read the letter

The main takeaway is that this work adds a jump-counting estimator to tensor-network solvers so currents can be tallied directly from tunneling events, letting them push simulations to arrays of 50 quantum dots. That is a concrete step past standard density-matrix methods for open systems. They benchmark the currents against QmeQ on smaller arrays across a range of lead-dot and inter-dot couplings and report quantitative agreement plus large drops in memory and runtime. Those parts look solid and useful on their own terms. The jump-counting trick avoids some of the overhead of full steady-state density matrices, which is a practical addition. The soft spot is the jump to N=50. The abstract states they can model up to fifty dots, but the quantitative checks are only in the tractable regime. There is no reported scaling of bond dimension, truncation error, or estimator variance with N, and no convergence tests that would confirm the errors stay controlled when the Hilbert space grows. If the required bond dimension rises faster than expected or if the jump estimator picks up a bias only visible at larger sizes, the 50-dot currents rest on an untested extrapolation. This is for people who need numerical transport results in mesoscopic systems larger than what master-equation codes can reach. A reader working on size-dependent nonequilibrium effects in dot arrays would get a workable method and the reported speedups, though they would still have to run their own error checks for the biggest cases. It deserves peer review because the core extension is well-motivated and the small-system evidence is clear; referees can sort out the large-N validation without much trouble.

Referee Report

2 major / 1 minor

Summary. The manuscript introduces a tensor-network method (TJM) that augments a tensor-based solver with a jump-counting estimator to compute steady-state electron currents in open quantum-dot arrays. It reports quantitative agreement with the QmeQ master-equation solver for small, tractable systems across various coupling parameters and claims orders-of-magnitude reductions in memory and wall-clock time, enabling simulations of interacting arrays with up to 50 quantum dots.

Significance. If the error control and accuracy claims hold at large N, the work would provide a practical route to nonequilibrium transport simulations in correlated quantum-dot systems far beyond the reach of density-matrix solvers, facilitating systematic studies of size-dependent phenomena in nanoscale devices.

major comments (2)

[Results and benchmarking sections] The central claim that the method enables reliable modeling of 50-dot arrays rests on an extrapolation from benchmarks performed only in the small-N regime where QmeQ is tractable. No data are presented on the scaling of required bond dimension, truncation error, or jump-counting estimator variance with N, which is load-bearing because the Hilbert space grows exponentially and uncontrolled growth in these quantities would invalidate the large-N results.
[Methods on the jump-counting estimator] The jump-counting estimator for lead-induced tunneling events is introduced to obtain currents directly from the tensor-network representation, but its unbiasedness and convergence properties at N=50 are not demonstrated; only agreement in the small-N limit is shown, leaving open the possibility of systematic bias that appears only at larger sizes.

minor comments (1)

[Abstract] The abstract contains a grammatical error in the final sentence ('Our approach allow us to model...').

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of the significance of our work and for the detailed comments. We address each major comment below and will update the manuscript to incorporate the suggested improvements.

read point-by-point responses

Referee: [Results and benchmarking sections] The central claim that the method enables reliable modeling of 50-dot arrays rests on an extrapolation from benchmarks performed only in the small-N regime where QmeQ is tractable. No data are presented on the scaling of required bond dimension, truncation error, or jump-counting estimator variance with N, which is load-bearing because the Hilbert space grows exponentially and uncontrolled growth in these quantities would invalidate the large-N results.

Authors: We acknowledge that explicit scaling data would strengthen the large-N claims. In the revised manuscript we will add figures in the results section showing bond dimension versus N, truncation error versus N for representative parameters, and jump-counting estimator variance versus N up to 50 dots. These will demonstrate that the tensor-network compression remains efficient and the monitored truncation error stays below a controlled threshold, supporting the reliability of the N=50 results even without direct current benchmarks at large N. revision: yes
Referee: [Methods on the jump-counting estimator] The jump-counting estimator for lead-induced tunneling events is introduced to obtain currents directly from the tensor-network representation, but its unbiasedness and convergence properties at N=50 are not demonstrated; only agreement in the small-N limit is shown, leaving open the possibility of systematic bias that appears only at larger sizes.

Authors: The jump-counting estimator follows directly from the quantum-jump unraveling of the Lindblad master equation and is therefore unbiased for any system size by construction. Its statistical convergence is governed by the number of trajectories, independent of N. In the revised methods section we will add a convergence analysis for the N=50 simulations, including plots of computed current versus number of trajectories that show the variance remains manageable and the estimator stabilizes to the same precision as in the small-N benchmarks. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation grounded by external benchmarks

full rationale

The paper presents a tensor-network method augmented by a jump-counting estimator for computing steady-state currents in open quantum-dot systems. It explicitly benchmarks the computed currents against the independent master-equation solver QmeQ, reporting quantitative agreement in the small-N regime where both are tractable. The extension to N=50 relies on the compression properties of the tensor representation and the estimator's direct counting of tunneling events, neither of which is defined in terms of the target results or fitted to the large-N data. No self-citations, ansatzes, or uniqueness theorems are invoked as load-bearing steps in the provided text, and no prediction is shown to reduce tautologically to an input parameter. The central claim therefore remains self-contained against external validation rather than circular by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Review limited to abstract; no specific free parameters, axioms, or new entities are described. The work builds on standard tensor-network techniques for open quantum systems and introduces the jump-counting estimator as the key addition.

pith-pipeline@v0.9.0 · 5467 in / 1270 out tokens · 83721 ms · 2026-05-10T16:58:48.398988+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

34 extracted references · 33 canonical work pages

[1]

P., Petta, J

Hanson, R., Kouwenhoven, L. P., Petta, J. R., Tarucha, S. & Vandersypen, L. M. K. Spins in few- electron quantum dots.Rev. Mod. Phys.79, 1217–1265 (2007).https://doi.org/10.1103/RevMod Phys.79.1217

work page doi:10.1103/revmod 2007
[2]

Zwanenburg, F. A. et al. Silicon quantum electronics.Rev. Mod. Phys.85, 961–1019 (2013).https: //doi.org/10.1103/RevModPhys.85.961

work page doi:10.1103/revmodphys.85.961 2013
[3]

Kouwenhoven, L. P. et al. Electron Transport in Quantum Dots. InMesoscopic Electron Transport, 105–214 (Springer, Dordrecht, 1997).https://doi.org/10.1007/978-94-015-8839-3_4 10

work page doi:10.1007/978-94-015-8839-3_4 1997
[4]

Kastner, M. A. The single-electron transistor.Rev. Mod. Phys.64, 849–858 (1992).https://doi. org/10.1103/RevModPhys.64.849

work page doi:10.1103/revmodphys.64.849 1992
[5]

Kastner, M. A. Artificial Atoms.Phys. Today46, 24–31 (1993).https://doi.org/10.1063/1.88 1393

work page doi:10.1063/1.88 1993
[6]

Computational Linguistics16(1), 22–29 (1990) https://doi.org/10

Breuer, H.-P. & Petruccione, F.The Theory of Open Quantum Systems(Oxford Univ. Press, Oxford, 2007).https://doi.org/10.1093/acprof:oso/9780199213900.001.0001

work page doi:10.1093/acprof:oso/9780199213900.001.0001 2007
[7]

Nazarov, Y. V. & Blanter, Y. M.Quantum Transport: Introduction to Nanoscience(Cambridge University Press, Cambridge, 2009).https://doi.org/10.1017/CBO9780511626906

work page doi:10.1017/cbo9780511626906 2009
[8]

Jacoboni, C.Theory of Electron Transport in Semiconductors: A Pathway from Elementary Physics to Nonequilibrium Green Functions(Springer, Berlin, Heidelberg, 2010).https://doi.org/10.100 7/978-3-642-10586-9

2010
[9]

Datta, S.Electronic Transport in Mesoscopic Systems(Cambridge University Press, Cambridge, 1995).https://doi.org/10.1017/CBO9780511805776

work page doi:10.1017/cbo9780511805776 1995
[10]

A short introduction to the Lindblad master equation.AIP Advances, 10(2):025106, 2020.doi:10.1063/1.5115323,arXiv:1906.04478

Manzano, D. A short introduction to the Lindblad Master Equation.AIP Adv.10, 025106 (2020). https://doi.org/10.1063/1.5115323

work page doi:10.1063/1.5115323 2020
[11]

& Fran¸ ca, V

Diniz, G., Quintino, S. & Fran¸ ca, V. V. Transport in Single Quantum Dots: A Review from Linear Response to Nonlinear Regimes.Braz. J. Phys.56, 27 (2026).https://doi.org/10.1007/s13538 -025-01953-0

work page doi:10.1007/s13538 2026
[12]

& Lee, E

Dorligjav, U., Seo, M. & Lee, E. Theoretical analysis and high-performance implementation of low- rank approximations in NEGF-based quantum transport.Commun. Nonlinear Sci. Numer. Simul. 152, 109265 (2026).https://doi.org/10.1016/j.cnsns.2025.109265

work page doi:10.1016/j.cnsns.2025.109265 2026
[13]

& L¨ auchli, A

Bonnes, L. & L¨ auchli, A. M. Superoperators vs. Trajectories for Matrix Product State Simulations of Open Quantum System: A Case Study. Preprint athttps://arxiv.org/abs/1411.4831(2014)

work page arXiv 2014
[14]

& Or´ us, R

Weimer, H., Kshetrimayum, A. & Or´ us, R. Simulation methods for open quantum many-body systems.Rev. Mod. Phys.93, 015008 (2021).https://doi.org/10.1103/RevModPhys.93.015008

work page doi:10.1103/revmodphys.93.015008 2021
[15]

A., Guti´ errez, R., Song, B

Ryndyk, D. A., Guti´ errez, R., Song, B. & Cuniberti, G. Green function techniques in the treatment of quantum transport at the molecular scale. InEnergy Transfer Dynamics in Biomaterial Systems (eds. Burghardt, I., May, V., Micha, D. A. & Bittner, E. R.) 213–335 (Springer, Berlin, Heidelberg, 2009).https://doi.org/10.1007/978-3-642-02306-4_9

work page doi:10.1007/978-3-642-02306-4_9 2009
[16]

Sander, A. et al. Large-scale stochastic simulation of open quantum systems.Nat. Commun.16, 11074 (2025).https://doi.org/10.1038/s41467-025-66846-x

work page doi:10.1038/s41467-025-66846-x 2025
[17]

& Dalibard, J

Mølmer, K., Castin, Y. & Dalibard, J. Monte Carlo wave-function method in quantum optics.J. Opt. Soc. Am. B10, 524–538 (1993).https://doi.org/10.1364/JOSAB.10.000524

work page doi:10.1364/josab.10.000524 1993
[18]

Daley, A. J. Quantum trajectories and open many-body quantum systems.Adv. Phys.63, 77–149 (2014).https://doi.org/10.1080/00018732.2014.933502

work page doi:10.1080/00018732.2014.933502 2014
[19]

Schollwöck,The density-matrix renormalization group in the age of matrix product states, Annals of Physics326(1), 96 (2011), doi:10.1016/j.aop.2010.09.012

Schollw¨ ock, U. The density-matrix renormalization group in the age of matrix product states.Ann. Phys.326, 96–192 (2011).https://doi.org/10.1016/j.aop.2010.09.012

work page doi:10.1016/j.aop.2010.09.012 2011
[20]

Saad,Numerical methods for large eigenvalue prob- lems: revised edition(SIAM, 2011), URLhttp://dx

Or´ us, R. Tensor networks for complex quantum systems.Nat. Rev. Phys.1, 538–550 (2019).https: //doi.org/10.1038/s42254-019-0086-7

work page doi:10.1038/s42254-019-0086-7 2019
[21]

Paeckel, S. et al. Time-evolution methods for matrix-product states.Ann. Phys.411, 167998 (2019). https://doi.org/10.1016/j.aop.2019.167998

work page doi:10.1016/j.aop.2019.167998 2019
[22]

N., Karlstr¨ om, O., Leijnse, M

Kirˇ sanskas, G., Pedersen, J. N., Karlstr¨ om, 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).https://doi.org/10.1016/j.cpc.2017.07.024

work page doi:10.1016/j.cpc.2017.07.024 2017
[23]

Prosen and M

Prosen, T. & ˇZnidariˇ c, M. Diffusive high-temperature transport in the one-dimensional Hubbard model.Phys. Rev. B86. 125118 (2012).https://doi.org/10.1103/PhysRevB.86.125118 11

work page doi:10.1103/physrevb.86.125118 2012
[24]

Rams, M. M. & Zwolak, M. Breaking the Entanglement Barrier: Tensor Network Simulation of Quantum Transport.Phys. Rev. Lett.124, 137701 (2020).https://doi.org/10.1103/PhysRevLet t.124.137701

work page doi:10.1103/physrevlet 2020
[25]

& von Delft, J

Schwarz, F., Goldstein, M., Dorda, A., Arrigoni, E., Weichselbaum, A. & von Delft, J. Lindblad- driven discretized leads for nonequilibrium steady-state transport in quantum impurity models: Recovering the continuum limit.Phys. Rev. B94, 155142 (2016).https://doi.org/10.1103/Ph ysRevB.94.155142

work page doi:10.1103/ph 2016
[26]

& Goldstein, M

Lotem, M., Weichselbaum, A., von Delft, J. & Goldstein, M. Renormalized Lindblad Driving: A Numerically-Exact Nonequilibrium Quantum Impurity Solver.Phys. Rev. Res.2. 043052 (2020). https://doi.org/10.1103/PhysRevResearch.2.043052

work page doi:10.1103/physrevresearch.2.043052 2020
[27]

& Guo, C

Chen, R., Xu, X. & Guo, C. Grassmann time-evolving matrix product operators for quantum impu- rity models.Phys. Rev. B109. 045140 (2024).https://doi.org/10.1103/PhysRevB.109.045140

work page doi:10.1103/physrevb.109.045140 2024
[28]

Tsupko, G.S

Thoenniss, J., Sonner, M., Lerose, A. & Abanin, D. A. Efficient method for quantum impurity problems out of equilibrium.Phys. Rev. B107. L201115 (2023).https://doi.org/10.1103/Phys RevB.107.L201115

work page doi:10.1103/phys 2023
[29]

L., Maniscalco, S

Cattaneo, M., Giorgi, G. L., Maniscalco, S. & Zambrini, R. Local vs global master equation with common and separate baths: superiority of the global approach in partial secular approximation. New J. Phys.21, 113045 (2019).https://doi.org/10.1088/1367-2630/ab54ac

work page doi:10.1088/1367-2630/ab54ac 2019
[30]

Haegeman, C

Haegeman, J., Lubich, C., Oseledets, I., Vandereycken, B. & Verstraete, F. Unifying time evolution and optimization with matrix product states.Phys. Rev. B94, 165116 (2016).https://doi.org/ 10.1103/PhysRevB.94.165116

work page doi:10.1103/physrevb.94.165116 2016
[31]

Haegeman, J. et al. Time-Dependent Variational Principle for Quantum Lattices.Phys. Rev. Lett. 107, 070601 (2011).https://doi.org/10.1103/PhysRevLett.107.070601

work page doi:10.1103/physrevlett.107.070601 2011
[32]

Campaioli, F., Cole, J. H. & Hapuarachchi, H. Quantum Master Equations: Tips and Tricks for Quantum Optics, Quantum Computing, and Beyond.PRX Quantum5, 020202 (2024).https: //doi.org/10.1103/PRXQuantum.5.020202

work page doi:10.1103/prxquantum.5.020202 2024
[33]

& Znidaric, M

Benenti, G., Casati, G., Prosen, T., Rossini, D. & Znidaric, M. Charge and spin transport in strongly correlated one-dimensional quantum systems driven far from equilibrium.Phys. Rev. B80, 035110 (2009).https://doi.org/10.1103/PhysRevB.80.035110

work page doi:10.1103/physrevb.80.035110 2009
[34]

T., Kewming, M

Landi, G. T., Kewming, M. J., Mitchison, M. T. & Potts, P. P. Current Fluctuations in Open Quan- tum Systems: Bridging the Gap Between Quantum Continuous Measurements and Full Counting Statistics.PRX Quantum5, 020201 (2024).https://doi.org/10.1103/PRXQuantum.5.020201 12 Supplementary Information Matrix-product states Matrix-product states (MPS) provide an...

work page doi:10.1103/prxquantum.5.020201 2024