Quantum master equation for nanoelectromechanical systems beyond the wide-band limit

Federico Cerisola; Janet Anders; Sofia Sevitz

REVIEW 3 major objections 5 minor 2 cited by

Quantum master equation for nanoelectromechanical systems beyond the wide-band limit

T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper derives a quantum master equation for a quantum dot coupled to a mechanical oscillator in the slow-tunneling regime, beyond the wide-band limit, and shows its non-equilibrium steady state matches numerically exact calculations.

desk verdict Useful new NEMS master equation with honest benchmarking, but the Lamb-shift neglect is not justified by the appendix and the benchmark is too narrow to cover the gap. read the letter →

arxiv 2506.20593 v1 pith:W5A3GK5Q submitted 2025-06-25 quant-ph cond-mat.mes-hall

Sofia Sevitz , Federico Cerisola , Janet Anders This is my paper · ORCID

classification quant-phcond-mat.mes-hall

keywords quantummasterequationnanoelectromechanicalsystemsRedfieldpolarontransformationslowtunnelingregimebeyondwide-bandlimitnon-equilibriumsteadystateenergy-dependentrates

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper derives a quantum master equation for a quantum dot coupled to a mechanical oscillator and two fermionic reservoirs, aimed at the regime where electrons tunnel more slowly than the oscillator vibrates ($\Gamma_\nu \ll \omega$). In that regime the usual semiclassical treatment is invalid, and the authors supply a fully quantum description that also keeps the tunneling rates energy-dependent instead of imposing the wide-band approximation. The central result, a projected second-order Redfield equation in a polaron frame, is benchmarked against numerically exact hierarchical equations-of-motion calculations and reproduces the steady state within about ten percent trace distance across coupling strengths. The paper also derives a closed particle-current expression that recovers thermal currents, Franck-Condon blockade, and the experimentally observed asymmetric suppression of conductance at Coulomb-diamond edges.

What carries the argument

The load-bearing object is the polaron transformation $U=e^{\lambda(\hat{b}^\dagger-\hat{b})\hat{d}^\dagger\hat{d}}$, which diagonalizes the dot-oscillator coupling by shifting the dot level to $\tilde{\mu}=\mu-\omega\lambda^2$ and dresses each electron-tunneling event with a displacement operator $\hat{D}(\lambda)=e^{\lambda(\hat{b}^\dagger-\hat{b})}$. From the von Neumann equation with an initially uncorrelated system-reservoir product state, a second-order Born-Markov reduction (the kernel $K_I(t,t')$ replaced by a Heaviside step) produces the projected Redfield equation (15). The resulting Redfield tensors contain rates evaluated at oscillator-shifted frequencies $\omega_{k,l}=\tilde{\mu}-\omega(k-l)$, so mechanical transitions enter the electronic rates directly. The same machinery yields a GKLS-type equation under a secular approximation valid when $\Gamma_\nu \ll 2\max\{\tilde{\mu},\omega\}$, and a particle-current expression obtained from the correlated part of the density matrix.

What would settle it

Take the same dot-oscillator-reservoir model with Lorentzian rates and compute the steady state without the Heaviside replacement, retaining the time-nonlocal kernel, at parameters near the boundary of the claimed regime, for example $\Gamma_\nu/\omega = 0.1$ with $\Gamma_\nu/T_\nu$ not very small; if the trace distance to Eq. (15) grows well beyond the roughly ten percent reported in Fig. 3, the Born-Markov step is the failing assumption.

Watch

Extended reading notes

Core claim

The central claim is that Eq. (15) of the paper, the projected second-order Redfield master equation in the polaron frame, correctly describes the non-equilibrium steady state of the quantum-dot-oscillator system when electronic tunneling is slower than the oscillator frequency. The Redfield tensors are built from energy-dependent transition rates, $R_\nu^{0\to1}(\epsilon)=\Upsilon_\nu(\epsilon)f_\nu(\epsilon)$ and $R_\nu^{1\to0}(\epsilon)=\Upsilon_\nu(\epsilon)[1-f_\nu(\epsilon)]$, with a Lorentzian spectral density $\Upsilon_\nu(\epsilon)$ that reduces to the wide-band limit as its width goes to infinity. Unlike earlier treatments, the oscillator density matrix is allowed to carry coherences, and the paper shows those coherences are necessary to match the exact steady state. The same derivation yields a particle-current formula, Eq. (22), that reduces to the standard quantum-dot current for zero coupling and reproduces known transport features.

Load-bearing premise

The calculation assumes the reservoirs forget their past instantly: the time-nonlocal kernel in Eq. (A12) is replaced by a Heaviside step, a move the paper describes as second-order perturbation while conceding that the precise link between the transport condition $\Gamma_\nu \ll T_\nu$ and open-systems weak coupling is not yet derived; if this memoryless replacement fails, the steady state could miss non-Markovian or higher-order corrections.

Editorial extensions

If this is right

Eq. (15) gives the missing fully quantum description of the non-equilibrium steady state in the slow-tunneling regime, where semiclassical Fokker-Planck or Langevin models are no longer valid.
Retaining energy-dependent tunneling rates changes the transport window and produces asymmetric suppression of the differential conductance at the edges of Coulomb diamonds, a feature reported in experiments.
The particle-current expression (Eq. (22)) reduces to the standard quantum-dot rate-equation current at zero dot-oscillator coupling and reproduces thermally driven currents and Franck-Condon blockade.
The polaron-frame Redfield steady state stays within roughly ten percent trace distance of the numerically exact solution for couplings up to $\lambda=1.5$, whereas ignoring oscillator coherences gives a noticeably worse match.
The secular GKLS version of the equation provides completely positive evolution under the stated parameter bound, making the model usable for quantum-thermodynamics applications.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A concrete derivation connecting the transport condition $\Gamma_\nu\ll T_\nu$ to the open-systems weak-coupling condition would sharpen the regime of validity; until then, the slow-tunneling boundary rests on the memoryless kernel replacement.
The model opens a route to study quantum self-oscillations, work extraction, and mechanical batteries in the slow-transport regime, applications the authors flag but do not develop.
The Lamb-shift analysis suggests a restriction $\beta_\nu\delta_\nu<\pi$ for Lorentzian spectral densities; extensions to multi-Lorentzian or non-Lorentzian bands should re-check such singularities before being used.
A natural experimental test is to measure the asymmetry of the Coulomb-diamond conductance suppression as a function of the reservoir spectral widths; the predicted dependence on $\delta_L$ and $\delta_R$ distinguishes the model from a wide-band treatment.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Dear colleague,

This paper derives a Redfield master equation for a quantum dot coupled to a mechanical oscillator in the slow-tunneling regime (Γ ≪ ω) with energy-dependent hopping rates, and gives an explicit particle current expression. The combination is genuinely new as far as I can tell: earlier work either assumed wide-band or fast tunneling. The derivation is careful, the polaron frame is appropriate, and the λ=0 limit recovers the standard QD rate equations. The HEOM benchmark in Figs. 3 and 4 is a real strength; for the single parameter set tested, the trace distance stays below 10%, and the oscillator coherences are captured, which a semiclassical approximation misses. The current expression is a useful byproduct.

The soft spots are real but not hidden. The authors state up front that the Lamb shift is neglected, and Appendix B tries to justify it. That justification does not work as written. The triangle inequality gives |Im G| ≤ |Re G| plus two positive terms; the claim that those terms scale as 1/N in the phonon number is not valid because the relevant E = μ̃ − ω(k−l) is set by populated levels and can be small or zero. And even if the positive terms were small, |Im| ≤ |Re| does not imply |Im| ≪ |Re|; it is a trivial bound. So Eq. (15) rests on an unjustified approximation. The HEOM comparison is reassuring for the one parameter set, but it is a single point in parameter space, and the neglected Lamb shift is the same order as Γ (0.1 GHz in the benchmark). The current expression (22) is not benchmarked at all. The Born-Markov / “high temperature limit” step in Appendix A is also admitted to lack a concrete derivation; that is another gap, though the authors flag it clearly.

I do not think these flaws are fatal to the paper's value. The approach is plausible, the numerics suggest it works where tested, and the community needs a model for this regime. But the paper currently overclaims: the appendix does not establish the bound it says it establishes, and the benchmark is too narrow to carry that weight alone.

I would send it to peer review, with the expectation that the authors either provide a correct estimate of the Lamb shift or present a broader HEOM comparison (varying temperatures, biases, and spectral densities) before publication. As it stands, I would trust the equations as a starting point, not as a settled result.

Yours,

[Name]

Referee Report

3 major / 5 minor

Summary. The manuscript derives a projected second-order Redfield master equation, Eq. (15), for a quantum dot coupled to a mechanical oscillator and fermionic reservoirs in the polaron frame, targeting the slow-tunneling regime where the semiclassical approximation fails. It goes beyond the wide-band limit by using Lorentzian energy-dependent tunneling rates, and it presents a closed particle-current expression, Eq. (22). The steady state is benchmarked against hierarchical equations of motion (HEOM) results with reported trace distances below 10%, and the current is compared with known experimental features such as thermal currents, Franck-Condon blockade, and Coulomb-diamond edge suppression. The authors also show that the master equation reduces to standard quantum-dot rate equations for vanishing dot-oscillator coupling.

Significance. If the result is valid, the paper fills a genuine gap: a fully quantum, non-wide-band treatment of NEMS in the slow-tunneling regime, including oscillator coherences that semiclassical models discard. The work has concrete strengths: the derivation is carried out in detail, there are no fitted parameters, the λ=0 limit correctly recovers standard rate equations, and the HEOM comparison provides an independent numerical check. The current expression is simple enough to be used by experimental groups. However, the central approximation of neglecting the Lamb shift is not justified by the analytic argument in Appendix B.1, and the numerical benchmark covers only a narrow set of parameters and does not validate the current expression Eq. (22). These issues are fixable in revision, but they are load-bearing for the claim that Eq. (15) correctly describes the non-equilibrium steady state.

major comments (3)

[Appendix B.1, Eqs. (B30)-(B33)] The argument for neglecting the Lamb shift is invalid. The triangle inequalities upper-bound |Im(G)| by |Re(G)| plus two positive remainder terms; because |Re(G)| appears additively, the bound cannot imply |Im(G)| is much smaller than |Re(G)|, and even if the remainder terms vanished the inequality would only give the trivial |Im(G)| ≤ |Re(G)|. The further claim that the remainder terms scale as 1/N because ω_{k,l}=μ̃−ω(k−l) is '∼ N' is not justified: the relevant frequency arguments are differences between oscillator levels that are actually populated, not the Hilbert-space truncation size N. Since the imaginary parts of the integrals in Eqs. (A32)-(A35) are discarded to arrive at Eq. (15), the central approximation is not backed by the presented analytic estimate. Please either include the Lamb-shift terms (the residue expressions in Eqs. (B25)-(B26) give them explicitly) or provide a numerical estimate of the imaginary part over the benchmark parameter range, including its effect on the steady state and on the current Eq. (22).
[Appendix A, Eq. (A12)] The replacement of the memory kernel K_I(t,t′) by the Heaviside function is a load-bearing step in the derivation of the Born-Markov Redfield equation. The manuscript itself states that this is 'understood as a second order perturbation' and that a concrete derivation connecting the transport-community high-temperature condition Γν ≪ Tν to the open-quantum-system weak-coupling condition is left to future work. Please make the small parameter explicit, for example a dimensionless ratio Γν divided by the relevant dot or oscillator energy scale, and state the resulting validity window. The single-parameter HEOM test in Fig. 3 does not yet map out this window, so the reader cannot currently assess where non-Markovian or higher-order corrections might become relevant.
[Section IV, Eq. (22) and Fig. 5] The particle-current expression is a central result of the paper, but it is not benchmarked against the numerically exact HEOM solution. The comparisons in Fig. 5 are qualitative, and the current inherits the Lamb-shift omission of Eq. (15), so the 'ready to use' claim needs quantitative support. Please benchmark I_R against HEOM for the parameter set of Figs. 3-4, and if feasible for a line of μ̃ and Δμ values, to show that the current expression is accurate and not accidentally reproducing only qualitative features.

minor comments (5)

[Section III.B, after Eq. (15)] The text 'we benchmark the sates obtained' contains a typo and should read 'we benchmark the states obtained'.
[Figure 2 caption] The caption lists 'T_L = 6.546 ... T_L = 5.237'; the second temperature should presumably be T_R, not T_L.
[Equations (16)-(19)] The superscript notation A^{n,n}, A^{n+1,n}, A^{n-1,n} is used before the convention for the lower index n is explained; please define the notation at first use in the main text.
[Figure 5(c) caption] The notation 'γ_L = −γ_R = −10 GHz' is ambiguous about the sign of γ_R; please write the two values explicitly.
[References] References [18] and [47] are listed as 'In preperation' (also misspelled); if these works are not yet available, mark them clearly as unpublished or remove them from the reference list.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the Redfield equation and current expression are derived from the microscopic Hamiltonian with stated assumptions, benchmarked against independent HEOM, and self-citations are not load-bearing.

full rationale

The paper's central derivation is self-contained. Equation (15) is obtained from the Hamiltonian in Eqs. (1)-(5) via the polaron transformation (10), the second-order perturbation step (A12), and the Markov/Redfield projection in Appendix A, with no parameter fitted to the target steady state or current. The benchmark in Sec. III C uses HEOM (Ref. [57]), an independent numerically exact method, so the comparison genuinely tests the Born-Markov and Lamb-shift approximations rather than re-inserting the result. The particle current expression (22) is derived from the same Redfield equation in Appendix E and is only compared qualitatively with external experimental features (e.g., Ref. [5]), so there is no fitted-input-called-prediction structure. The self-citations present (Refs. [18] and [47] are 'In preparation' works by overlapping authors; Refs. [3] and [11] include overlapping experimental authors) are motivational or experimental and are not load-bearing premises of the derivation. The paper itself flags its main limitations: Appendix A postpones a concrete derivation connecting the transport-community 'high temperature limit' to open-quantum-system weak coupling ('Future work will focus on a concrete derivation'), and Appendix B1 gives a loose, not rigorous, bound for neglecting the Lamb shift. These are validity/soundness concerns, not circular reductions, and the HEOM comparison partially supports the neglect in the benchmarked regime. Overall, no equation in the claimed derivation chain reduces by construction to its own inputs.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The derivation assumes standard open-quantum-system machinery: weak coupling, Markovianity, uncorrelated initial state, and polaron transformation. The only model-specific input is the Lorentzian spectral density. No new physical entities are posited.

free parameters (3)

Γ_L, Γ_R (dot-reservoir tunneling rates) = 0.1 GHz (benchmark)
Model inputs chosen to satisfy the slow-tunneling and weak-coupling conditions; not fitted to target results.
γ_L, γ_R (Lorentzian spectral density centers) = γ_L=2.5 GHz, γ_R=-2.5 GHz (Figs. 3-4); γ_L=-10 GHz, γ_R=10 GHz (Fig. 5)
Hand-chosen centers for the Lorentzian model used in the demonstrations; they are physical model parameters, not fitted to the benchmark.
δ_L, δ_R (Lorentzian spectral density widths) = δ_L=δ_R=2 GHz (Figs. 3-4); δ_L=20 GHz, δ_R=10 GHz (Fig. 5)
Hand-chosen widths; must satisfy δ >> Γ and βδ < π for validity.

assumptions (6)

domain assumption Born-Markov approximation: kernel K_I(t,t') is replaced by the Heaviside function in Eq. (A12), making the master equation time-local and second order in the dot-reservoir coupling.
Invoked in Appendix A between Eqs. (A11) and (A12); the authors acknowledge the justification is incomplete and defer it to future work.
standard math Initial uncorrelated system-reservoir state at t0=-∞.
Standard open-quantum-systems assumption, stated in Sec. III B and Eq. (A9).
domain assumption Particle superselection rule: only diagonal elements of the dot density matrix are kept.
Used to project the Redfield equation onto dot occupation states, see App. A 1.
domain assumption Lorentzian spectral density for tunneling rates (Eq. 5).
The master equation is derived for this specific energy dependence; extension to superpositions of Lorentzians is sketched in App. B.
domain assumption Neglect of the Lamb shift.
The imaginary parts of the bath correlation integrals are dropped; App. B gives a loose bound (terms scaling as 1/N) rather than a rigorous justification.
domain assumption Regime conditions Γν << Tν, δν >> Γν, βν δν < π.
Needed for the second-order and Markov approximations and to avoid divergences, as stated in Sec. III B and App. B.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Quantum master equation for nanoelectromechanical systems beyond the wide-band limit." pith.science (2026). https://pith.science/paper/W5A3GK5Q

@misc{pith2026250620593,
  author       = {Pith},
  title        = {Pith review of: Quantum master equation for nanoelectromechanical systems beyond the wide-band limit},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/W5A3GK5Q}},
  note         = {Machine review of arXiv:2506.20593}
}

read the original abstract

Coupling the vibrations of an oscillator to electronic transport is a key building block for nanoelectromechanical systems. They describe many nanoscale electrical components such as molecular junctions. Inspired by recent experimental developments, we derive a quantum master equation that describes nanoelectromechanical systems in a generally overlooked situation: when the electronic transport is slower than the natural frequency of the oscillator. Here, a semi-classical model is no longer valid and we develop the missing fully quantum approach. Moreover, we go beyond the wide-band limit and study the consequence of maintaining energy dependent tunneling rates, which are required to describe effects found in real devices. To benchmark our results, we compare with numerically exact results obtained with the hierarchical equations of motion method, and find overall good agreements in the experimentally accessible steady state regime. Furthermore, we derive from the microscopic model a ready to use particle current expression that replicates features already observed experimentally.

Figures

Figures reproduced from arXiv: 2506.20593 by the authors.

**Figure 1.** Illustration of a nanoelectromechanical system. A quantum dot with chemical potential µ is coupled (with arbitrary strength g) to a quantum harmonic oscillator with natural frequency ω. Moreover, the quantum dot is coupled to two fermionic reservoirs to the left (L) and right (R) with temperatures TL and TR, and chemical potentials µL and µR, respectively. A temperature and/or chemical imbalance between the reser… view at source ↗

**Figure 2.** Illustrative example to showcase the physical consequence of going beyond the wide-band limit. (a) Energy spectrum of the dot and the two reservoirs under the wide-band limit (δL, δR → ∞ in Eq. (5)). The shaded area denotes the fermionic occupation corresponding to the left reservoir (red) and right reservoir (blue). The shaded orange area denotes the transport window of the dot, that is the energy range where elec… view at source ↗

**Figure 3.** Benchmark of non-equilibrium steady state ρs compared to σs. Trace distance between the non-equilibrium steady state using the Redfield equation derived in this work ρs and the exact solution computed using HEOM σs as a function of QD+QHO coupling λ in blue. For comparison, in gray we exhibit the trace distance between a classical analogous of our solution ρ diag s with σs. Parameters: ΓL/2π = ΓR/2π = 0.1 GHz, ω/2π… view at source ↗

Figures from the paper (4 more)

**Figure 4.** Figure 4: Non-equilibrium steady state elements of the QHO density matrix ρQHO. The upper rows are solved using HEOM for increasing QD+QHO couplings λ = 0.5 (a-I), λ = 1 (a-II), and λ = 1.5 (a-III). The lower rows are the corresponding elements obtained using the Redfield equat…

**Figure 5.** Figure 5: Particle current obtained from the microscopic model. (a) Particle current corresponding to the right reservoir, IR as function of the dot’s normalized energy ˜µ for increasing dot-resonator coupling strength λ. Here there is no chemical potential imbalance ∆µ = µL − µ…

**Figure 6.** Figure 6: Integration curves used to compute the Lamb-shift. Panel (a) corresponds to the curve used for Eq. (B9) while panel (b) was used for Eq. (B10). For simplicity, in the graphics we consider E ≥ 0. (Eq. (B4)). We recover these expressions with the integration over the rea…

**Figure 7.** Figure 7: (a)-(b) Integration functions of the bath correlation functions as a function of frequency ω ′ (c)-(d) Bath correlation functions detailed in Eq. (A21) and (A22) respectively as a function of time s. Parameters: ΓL/2π = ΓR/2π = 0.1 GHz, ω/2π= 1 GHz, TL = TR = 2 GHz (co…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Autonomous oscillations in quantum electromechanics: tensor network treatment
quant-ph 2026-05 unverdicted novelty 7.0 of 10

A tensor-network method maps out when and how electron transport through a vibrating quantum dot produces stable self-sustained mechanical oscillations.
Autonomous conversion of particle-exchange to quantum self-oscillations
quant-ph 2025-08 unverdicted novelty 6.0 of 10

A quantum-dot particle-exchange machine coupled to a mechanical resonator produces sustained self-oscillations in the slow-transport regime, detectable in the electrical current, and only in the heater regime.

Reference graph

Works this paper leans on

73 extracted references · 61 canonical work pages · cited by 2 Pith papers

[1]

H. G. Craighead, Nanoelectromechanical Systems, Sci- ence290, 1532 (2000)

work page 2000
[2]

Blencowe, Nanoelectromechanical systems, Contem- porary Physics46, 249 (2005)

M. Blencowe, Nanoelectromechanical systems, Contem- porary Physics46, 249 (2005)

work page 2005
[3]

Vigneau, J

F. Vigneau, J. Monsel, J. Tabanera, K. Aggarwal, L. Bresque, F. Fedele, F. Cerisola, G. A. D. Briggs, J. An- ders, J. M. R. Parrondo, A. Auff` eves, and N. Ares, Ultra- strong coupling between electron tunneling and mechan- ical motion, Phys. Rev. R4, 043168 (2022)

work page 2022
[4]

Luo, Z.-Z

G. Luo, Z.-Z. Zhang, G.-W. Deng, H.-O. Li, G. Cao, M. Xiao, G.-C. Guo, and G.-P. Guo, Coupling graphene nanomechanical motion to a single-electron transistor, Nanoscale9, 5608 (2017)

work page 2017
[5]

H. B. Meerwaldt, G. Labadze, B. H. Schneider, A. Taspinar, Y. M. Blanter, H. S. J. van der Zant, and G. A. Steele, Probing the charge of a quantum dot with a nanomechanical resonator, Phys. Rev. B86, 115454 (2012)

work page 2012
[6]

G. A. Steele, A. K. H¨ uttel, B. Witkamp, M. Poot, H. B. Meerwaldt, L. P. Kouwenhoven, and H. S. J. van der Zant, Strong Coupling Between Single-Electron Tun- neling and Nanomechanical Motion, Science325, 1103 (2009)

work page 2009
[7]

Samanta, S

C. Samanta, S. L. De Bonis, C. B. Møller, R. Tormo- Queralt, W. Yang, C. Urgell, B. Stamenic, B. Thibeault, Y. Jin, D. A. Czaplewski, F. Pistolesi, and A. Bach- told, Nonlinear nanomechanical resonators approaching the quantum ground state, Nature Physics19, 1340–1344 (2023)

work page 2023
[8]

E. A. Laird, F. Pei, W. Tang, G. A. Steele, and L. P. Kouwenhoven, A high quality factor carbon nanotube mechanical resonator at 39 ghz, Nano Letters12, 193 (2012)

work page 2012

Show all 73 references

[9]

Moser, A

J. Moser, A. Eichler, J. G¨ uttinger, M. I. Dykman, and A. Bachtold, Nanotube mechanical resonators with quality factors of up to 5 million, Nature Nanotech9, 1007–1011 (2014)

2014
[10]

D. R. Schmid, P. L. Stiller, C. Strunk, and A. K. H¨ uttel, Liquid-induced damping of mechanical feedback effects in single electron tunneling through a suspended carbon nanotube, Applied Physics Letters107, 123110 (2015)

2015
[11]

Tabanera-Bravo, F

J. Tabanera-Bravo, F. Vigneau, J. Monsel, K. Aggar- wal, L. Bresque, F. Fedele, F. Cerisola, G. A. D. Briggs, J. Anders, A. Auff` eves, J. M. R. Parrondo, and N. Ares, Stability of long-sustained oscillations induced by elec- tron tunneling, Phys. Rev. R6, 013291 (2024)

2024
[12]

Leijnse, M

M. Leijnse, M. R. Wegewijs, and K. Flensberg, Nonlin- ear thermoelectric properties of molecular junctions with vibrational coupling, Phys. Rev. B82, 045412 (2010)

2010
[13]

Mitra, I

A. Mitra, I. Aleiner, and A. J. Millis, Phonon effects in molecular transistors: Quantal and classical treatment, Phys. Rev. B69, 245302 (2004)

2004
[14]

Erpenbeck, C

A. Erpenbeck, C. Schinabeck, U. Peskin, and M. Thoss, Current-induced bond rupture in single-molecule junc- tions, Phys. Rev. B97, 235452 (2018)

2018
[15]

A. K. H¨ uttel, G. A. Steele, B. Witkamp, M. Poot, L. P. Kouwenhoven, and H. S. J. van der Zant, Carbon Nanotubes as Ultrahigh Quality Factor Mechanical Res- onators, Nano Letters9, 2547 (2009)

2009
[16]

Chaste, A

J. Chaste, A. Eichler, J. Moser, G. Ceballos, R. Rurali, and A. Bachtold, A nanomechanical mass sensor with yoctogram resolution, Nature Nanotech7, 301 (2012)

2012
[17]

De and B

B. De and B. Muralidharan, Thermoelectric study of dissipative quantum-dot heat engines, Phys. Rev. B94, 165416 (2016)

2016
[18]

Sevitz, K

S. Sevitz, K. Aggarwal, F. Cerisola, J. Tabanera-Bravo, J. Monsel, F. Vigneau, F. Fedele, J. Dunlop, J. M. Parrondo, G. J. Milburn, J. Anders, and N. Ares, Sources of nonlinearity in the response of a driven nano- electromechanical resonator, In preperation (2025)

2025
[19]

Y. Wen, N. Ares, F. J. Schupp, T. Pei, G. A. D. Briggs, and E. A. Laird, A coherent nanomechanical oscillator driven by single-electron tunnelling, Nature Physics16, 75 (2020)

2020
[20]

C. W. W¨ achtler, P. Strasberg, S. H. L. Klapp, G. Schaller, and C. Jarzynski, Stochastic thermodynam- ics of self-oscillations: the electron shuttle, New Journal of Physics21, 073009 (2019)

2019
[21]

Culhane, M

O. Culhane, M. T. Mitchison, and J. Goold, Extractable work in quantum electromechanics, Phys. Rev. E106, L032104 (2022)

2022
[22]

Schaller,Open Quantum Systems Far from Equilib- rium, Lecture Notes in Physics, Vol

G. Schaller,Open Quantum Systems Far from Equilib- rium, Lecture Notes in Physics, Vol. 881 (Springer Inter- national Publishing, Cham, 2014)

2014
[23]

A. A. Clerk and S. Bennett, Quantum nanoelectrome- chanics with electrons, quasi-particles and Cooper pairs: effective bath descriptions and strong feedback effects, New Journal of Physics7, 238 (2005). 8

2005
[24]

Pistolesi and S

F. Pistolesi and S. Labarthe, Current blockade in classical single-electron nanomechanical resonator, Phys. Rev. B 76, 165317 (2007)

2007
[25]

Micchi, R

G. Micchi, R. Avriller, and F. Pistolesi, Mechanical Sig- natures of the Current Blockade Instability in Suspended Carbon Nanotubes, Phys. Rev. Lett.115, 206802 (2015)

2015
[26]

A. W. Barnard, M. Zhang, G. S. Wiederhecker, M. Lip- son, and P. L. McEuen, Real-time vibrations of a carbon nanotube, Nature566, 89 (2019)

2019
[27]

X. Wang, L. Cong, D. Zhu, Z. Yuan, X. Lin, W. Zhao, Z. Bai, W. Liang, X. Sun, G.-W. Deng, and K. Jiang, Visualizing nonlinear resonance in nanomechanical sys- tems via single-electron tunneling, Nano Research14, 1156 (2021)

2021
[28]

Timm, Tunneling through molecules and quantum dots: Master-equation approaches, Phys

C. Timm, Tunneling through molecules and quantum dots: Master-equation approaches, Phys. Rev. B77, 195416 (2008)

2008
[29]

C. Schinabeck,Hierarchical quantum master equa- tion approaches to nonequilibrium charge transport through single-molecule junctions, Phd thesis, Friedrich- Alexander-Universit¨ at Erlangen-N¨ urnberg (F AU) (2019)

2019
[30]

Piovano, F

G. Piovano, F. Cavaliere, E. Paladino, and M. Sassetti, Coherent properties of nanoelectromechanical systems, Phys. Rev. B83, 245311 (2011)

2011
[31]

P. P. Potts, Quantum thermodynamics (2024), arXiv:2406.19206 [quant-ph]

2024 arXiv
[32]

A similar argument can also be made using the Jordan- Wigner transformation
[33]

Kirˇ sanskas, J

G. Kirˇ sanskas, J. N. Pedersen, O. Karlstr¨ om, M. Lei- jnse, and A. Wacker, QmeQ 1.0: An open-source Python package for calculations of transport through quantum dot devices, Computer Physics Communications221, 317 (2017)

2017
[34]

G. D. Mahan and J. O. Sofo, The best thermoelectric., Proceedings of the National Academy of Sciences93, 7436 (1996)

1996
[35]

R. D. Mayrhofer, C. Elouard, J. Splettstoesser, and A. N. Jordan, Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine, Phys. Rev. B103, 075404 (2021)

2021
[36]

Esposito, N

M. Esposito, N. Kumar, K. Lindenberg, and C. Van Den Broeck, Stochastically driven single-level quantum dot: A nanoscale finite-time thermodynamic machine and its various operational modes, Phys. Rev. E85, 031117 (2012)

2012
[37]

J. M. Thijssen and H. S. J. Van Der Zant, Charge trans- port and single-electron effects in nanoscale systems, physica status solidi (b)245, 1455 (2008)

2008
[38]

A. Usui, K. Ptaszy´ nski, M. Esposito, and P. Stras- berg, Microscopic contributions to the entropy produc- tion at all times: from nonequilibrium steady states to global thermalization, New Journal of Physics26, 023049 (2024)

2024
[39]

Moulhim, B

A. Moulhim, B. Tripathi, and M. Kumar, Nonequilibrium green function technique for analyzing electron transport through single and two levels of interacting quantum dot, Physica Scripta96, 125802 (2021)

2021
[40]

P. P. Potts, A. A. S. Kalaee, and A. Wacker, A thermody- namically consistent markovian master equation beyond the secular approximation, New Journal of Physics23, 123013 (2021)

2021
[41]

B. Wang, Y. Xing, L. Zhang, and J. Wang, Transient dynamics of molecular devices under a steplike pulse bias, Phys. Rev. B81, 121103 (2010)

2010
[42]

Cheng, X

J. Cheng, X. Zuo, J. Wang, and Y. Xing, Quasi- particle trapping in quench dynamics of superconduc- tor/quantum dot/ superconductor josephson junctions, Phys. Rev. B110, 125417 (2024)

2024
[43]

Y. M. Blanter, O. Usmani, and Y. V. Nazarov, Single- electron tunneling with strong mechanical feedback, Phys. Rev. Lett.93, 136802 (2004)

2004
[44]

Erpenbeck, C

A. Erpenbeck, C. Hertlein, C. Schinabeck, and M. Thoss, Extending the hierarchical quantum master equation ap- proach to low temperatures and realistic band structures, The Journal of Chemical Physics149, 064106 (2018)

2018
[45]

Schaller, T

G. Schaller, T. Krause, T. Brandes, and M. Esposito, Single-electron transistor strongly coupled to vibrations: counting statistics and fluctuation theorem, New Journal of Physics15, 033032 (2013)

2013
[46]

Dorsch,Transport in nanowire-based quantum dot sys- tems: Heating electrons and confining holes, Ph.D

S. Dorsch,Transport in nanowire-based quantum dot sys- tems: Heating electrons and confining holes, Ph.D. thesis, Department of Physics, Lund University (2022)

2022
[47]

Sevitz, F

S. Sevitz, F. Cerisola, K. V. Hovhannisyan, and J. An- ders, Autonomous conversion of particle exchange to quantum self-oscillations, In preperation (2025)

2025
[48]

C. Wu, Y. D. Y. Yan, Y. Su, E. O. Ayieta, S. Radoˇ sevi´ c, G. Engelhardt, G. Schaller, and J. Luo, Anomalous current-electric field characteristics in trans- port through a nanoelectromechanical systems (2025), arXiv:2503.12106 [cond-mat]

2025 arXiv
[49]

D. W. Utami, H.-S. Goan, C. A. Holmes, and G. J. Mil- burn, Quantum noise in the electromechanical shuttle: Quantum master equation treatment, Phys. Rev. B74, 014303 (2006)

2006
[50]

Pistolesi and R

F. Pistolesi and R. Fazio, Charge Shuttle as a Nanome- chanical Rectifier, Phys. Rev. Lett.94, 036806 (2005)

2005
[51]

Nazir and D

A. Nazir and D. P. S. McCutcheon, Modelling exci- ton–phonon interactions in optically driven quantum dots, J. Phys.: Condens. Matter28, 103002 (2016)

2016
[52]

W¨ urger, Strong-coupling theory for the spin-phonon model, Phys

A. W¨ urger, Strong-coupling theory for the spin-phonon model, Phys. Rev. B57, 347 (1998)

1998
[53]

Kolli, A

A. Kolli, A. Nazir, and A. Olaya-Castro, Electronic exci- tation dynamics in multichromophoric systems described via a polaron-representation master equation, J. Chem. Phys.135, 154112 (2011)

2011
[54]

Siddiqui, A

L. Siddiqui, A. W. Ghosh, and S. Datta, Phonon run- away in carbon nanotube quantum dots, Phys. Rev. B 76, 085433 (2007)

2007
[55]

Zedler, G

P. Zedler, G. Schaller, G. Kiesslich, C. Emary, and T. Brandes, Weak-coupling approximations in non- Markovian transport, Phys. Rev. B80, 045309 (2009)

2009
[56]

S. D. Bartlett, T. Rudolph, and R. W. Spekkens, Refer- ence frames, superselection rules, and quantum informa- tion, Rev. Mod. Phys.79, 555 (2007)

2007
[57]

Huang, P.-C

Y.-T. Huang, P.-C. Kuo, N. Lambert, M. Cirio, S. Cross, S.-L. Yang, F. Nori, and Y.-N. Chen, An efficient Ju- lia framework for hierarchical equations of motion in open quantum systems, Communications Physics6, 313 (2023)

2023
[58]

M. A. Nielsen and I. L. Chuang,Quantum computation and quantum information(Cambridge university press, 2010)

2010
[59]

This was not a problem previously when we studied the non-equilibrium steady state given that this condition is reduced to making sure that the bath correlation function decay at some time scale [66], see Appendix F for details
[60]

Breuer and F

H.-P. Breuer and F. Petruccione,The theory of open quantum systems(OUP Oxford, 2002)

2002
[61]

Hartmann and W

R. Hartmann and W. T. Strunz, Accuracy assessment of 9 perturbative master equations: Embracing nonpositivity, Phys. Rev. A101, 012103 (2020)

2020
[62]

Montoya-Castillo, T

A. Montoya-Castillo, T. C. Berkelbach, and D. R. Reich- man, Extending the applicability of Redfield theories into highly non-Markovian regimes, The Journal of Chemical Physics143, 194108 (2015)

2015
[63]

Josefsson, A

M. Josefsson, A. Svilans, A. M. Burke, E. A. Hoffmann, S. Fahlvik, C. Thelander, M. Leijnse, and H. Linke, A quantum-dot heat engine operating close to the ther- modynamic efficiency limits, Nature Nanotechnology13, 920 (2018)

2018
[64]

Dorsch, S

S. Dorsch, S. Fahlvik, and A. Burke, Characterization of electrostatically defined bottom-heated InAs nanowire quantum dot systems, New Journal of Physics23, 125007 (2021)

2021
[65]

Koch and F

J. Koch and F. von Oppen, Franck-condon blockade and giant fano factors in transport through single molecules, Phys. Rev. Lett.94, 206804 (2005)

2005
[66]

Strasberg,Quantum Stochastic Thermodynamics: Foundations and Selected Applications(Oxford Univer- sity Press, 2022)

P. Strasberg,Quantum Stochastic Thermodynamics: Foundations and Selected Applications(Oxford Univer- sity Press, 2022)

2022
[67]

Lambert, E

N. Lambert, E. Gigu` ere, P. Menczel, B. Li, P. Hopf, G. Su´ arez, M. Gali, J. Lishman, R. Gadhvi, R. Agar- wal, A. Galicia, N. Shammah, P. Nation, J. R. Johans- son, S. Ahmed, S. Cross, A. Pitchford, and F. Nori, Qutip 5: The quantum toolbox in python (2024), arXiv:2412.04705...

2024 arXiv
[68]

(8) of the main text it is convenient to project the density matrix into the energy eigenbasis [see Eq

Projection of Redfield To compute the non-equilibrium steady state defined in Eq. (8) of the main text it is convenient to project the density matrix into the energy eigenbasis [see Eq. (13) of the main text]. First we return to the Schr¨ odinger picture as ∂t ˜ρs(t) =−i[ ˜Hs,...
[69]

(A24) in order to bring it to GKLS form

GKLS form In this section we make an additional approximation (the secular approximation) over Eq. (A24) in order to bring it to GKLS form. Before this, we write the system operators in the interaction picture in a Fourier series [31] ( ˜S0)I (t) = X α e−iωαt ˜Sα 0 ( ˜S1)I (t)...
[70]

g0(z∓E) z = iΓνδν 2z2 0 [1−f ν(z2 0 ∓E)] (B17) Res g1 z , z1 1 = lim z→z 1 1 (z−z 1
[71]

g1(z∓E) z = −iΓνδν 2z1 1 fν(z1 1 ∓E) (B18) Res g0 z , zk 0 = lim z→z k 0 (z−z k 0 ) g0(∓E) z = +1 βν Υν(zk 0 ∓E) zk 0 (B19) Res g1 z , zk 1 = lim z→z k 1 (z−z k 1 ) g1(z∓E) z = −1 βν Υν(zk 1 ∓E) zk 1 .(B20) Where in the last two lines we implemented L’Hˆ opital’s rule. Next, w...
[72]

(B25) and Eq

Bound on the Lamb shift In this section we study what are the conditions over the physical parameters such that the Lamb shift computed in Eq. (B25) and Eq. (B26) can be neglected. In other words we want to know when|Im(G ν 0/1(E))| ≪ |Re(Gν 0/1(E))| 19 is satisfied. For this,...
[73]

(A49) satisfied

V alidity of secular approximation In this section we study in detail under what conditions the secular proximation carried-out in Appendix A 2 is valid, i.e., under what conditions is Eq. (A49) satisfied. For clarity we re-write Eq. (A49) below |γν k (ωα)| ≪ |ωα −ω α′| ∀ω α ̸...

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.

[1] [1]

H. G. Craighead, Nanoelectromechanical Systems, Sci- ence290, 1532 (2000)

work page 2000

[2] [2]

Blencowe, Nanoelectromechanical systems, Contem- porary Physics46, 249 (2005)

M. Blencowe, Nanoelectromechanical systems, Contem- porary Physics46, 249 (2005)

work page 2005

[3] [3]

Vigneau, J

F. Vigneau, J. Monsel, J. Tabanera, K. Aggarwal, L. Bresque, F. Fedele, F. Cerisola, G. A. D. Briggs, J. An- ders, J. M. R. Parrondo, A. Auff` eves, and N. Ares, Ultra- strong coupling between electron tunneling and mechan- ical motion, Phys. Rev. R4, 043168 (2022)

work page 2022

[4] [4]

Luo, Z.-Z

G. Luo, Z.-Z. Zhang, G.-W. Deng, H.-O. Li, G. Cao, M. Xiao, G.-C. Guo, and G.-P. Guo, Coupling graphene nanomechanical motion to a single-electron transistor, Nanoscale9, 5608 (2017)

work page 2017

[5] [5]

H. B. Meerwaldt, G. Labadze, B. H. Schneider, A. Taspinar, Y. M. Blanter, H. S. J. van der Zant, and G. A. Steele, Probing the charge of a quantum dot with a nanomechanical resonator, Phys. Rev. B86, 115454 (2012)

work page 2012

[6] [6]

G. A. Steele, A. K. H¨ uttel, B. Witkamp, M. Poot, H. B. Meerwaldt, L. P. Kouwenhoven, and H. S. J. van der Zant, Strong Coupling Between Single-Electron Tun- neling and Nanomechanical Motion, Science325, 1103 (2009)

work page 2009

[7] [7]

Samanta, S

C. Samanta, S. L. De Bonis, C. B. Møller, R. Tormo- Queralt, W. Yang, C. Urgell, B. Stamenic, B. Thibeault, Y. Jin, D. A. Czaplewski, F. Pistolesi, and A. Bach- told, Nonlinear nanomechanical resonators approaching the quantum ground state, Nature Physics19, 1340–1344 (2023)

work page 2023

[8] [8]

E. A. Laird, F. Pei, W. Tang, G. A. Steele, and L. P. Kouwenhoven, A high quality factor carbon nanotube mechanical resonator at 39 ghz, Nano Letters12, 193 (2012)

work page 2012