Thermodynamic coprocessor for linear operations with input-size-independent calculation time based on open quantum system
Pith reviewed 2026-05-18 18:04 UTC · model grok-4.3
pith:TJSVFQMG Add to your LaTeX paper
What is a Pith Number?\usepackage{pith}
\pithnumber{TJSVFQMG}
Prints a linked pith:TJSVFQMG badge after your title and writes the identifier into PDF metadata. Compiles on arXiv with no extra files. Learn more
The pith
An open quantum system of bosonic modes reaches a stationary state whose energy flows encode vector-matrix products for stochastic matrices in time fixed by relaxation alone.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
An open quantum system consisting of bosonic modes interacting with bosonic reservoirs implements multiple vector-matrix multiplications with stochastic matrices in parallel. Input vectors are encoded in the occupancies of the reservoirs, and the output result is presented by stationary energy flows. The operation takes time needed for the system's transition to a non-equilibrium stationary state independently on the number of the reservoirs, i.e., on the input vector dimension.
What carries the argument
The non-equilibrium stationary state of the open quantum system, whose stationary energy flows between the modes and reservoirs directly give the results of the parallel multiplications by stochastic matrices.
If this is right
- Multiple vector-matrix multiplications with stochastic matrices run simultaneously inside one device.
- Computation time is set only by the relaxation time to the non-equilibrium stationary state and stays fixed when the number of reservoirs grows.
- Dissipation rates multiplied by mode frequencies act as conductivities, reservoir occupancies act as potentials, and stationary energy flows act as currents, reproducing the behavior of an electrical crossbar.
- The computations are accompanied by entropy growth in the reservoirs and the system.
Where Pith is reading between the lines
- Devices built this way could be tiled across a small area to reach high operation rates per mode while remaining passive after the initial setup.
- The constant-time scaling may help neural-network layers whose matrix sizes increase without forcing longer wait times per layer.
- Because the mapping to crossbars is direct, existing analog-circuit design tools could be repurposed to engineer the required dissipation rates and couplings.
Load-bearing premise
The non-equilibrium stationary state reached by the open quantum system encodes the exact result of the vector-matrix multiplication in its stationary energy flows for arbitrary stochastic matrices, with no additional transient or correction terms that would depend on input dimension.
What would settle it
Prepare a small test case with known reservoir occupancies and a chosen stochastic matrix, allow the system to reach its non-equilibrium stationary state, and check whether the measured energy flows equal the direct matrix-vector product within experimental precision.
Figures
read the original abstract
Linear operations, e.g., vector-matrix and vector-vector multiplications, are core operations of modern neural networks. To diminish computational time, these operations are implemented by parallel computations using different coprocessors. In this work we show that an open quantum system consisting of bosonic modes and interacting with bosonic reservoirs can be used as an analog thermodynamic coprocessor implementing multiple vector-matrix multiplications with stochastic matrices in parallel. Input vectors are encoded in occupancies of reservoirs, and the output result is presented by stationary energy flows. The operation takes time needed for the system's transition to a non-equilibrium stationary state independently on the number of the reservoirs, i.e., on the input vector dimension. With technological limitations being considered, a device of $5\times5$ cm$^2$ area covered with the coprocessors can conduct of the order of $10^{11}$ operations per second per a mode of the OQS. The computations are accompanied by an entropy growth. We construct a direct mapping between open quantum systems and electrical crossbar structures frequently used in analog vector-matrix multiplication, showing that dissipation rates multiplied by open quantum system's modes frequencies can be seen as conductivities, reservoirs' occupancies can be seen as potentials, and stationary energy flows can be seen as electric currents.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes using an open quantum system of bosonic modes coupled to multiple bosonic reservoirs as an analog thermodynamic coprocessor that performs parallel vector-matrix multiplications with stochastic matrices. Reservoir occupancies encode the input vectors, stationary energy flows in the non-equilibrium steady state encode the outputs, and the time to reach this state is independent of the number of reservoirs (i.e., input dimension). A direct mapping to electrical crossbar arrays is constructed by identifying dissipation rates times mode frequencies with conductivities, reservoir occupancies with potentials, and energy flows with currents; a performance estimate of order 10^11 operations per second per mode on a 5×5 cm² device is also given.
Significance. If the exact correspondence between NESS flows and stochastic matrix multiplication holds without dimension-dependent corrections, the proposal would supply a physical platform for linear operations whose runtime does not grow with input size, which is potentially significant for neuromorphic or analog accelerators. The explicit thermodynamic framing and the crossbar analogy are constructive strengths that could guide experimental realizations, though practical impact would still require verification of the mapping and assessment of decoherence and control overhead.
major comments (2)
- [Abstract and mapping to crossbar structures] The central claim that stationary energy flows exactly equal the matrix-vector product for arbitrary stochastic matrices (abstract and mapping paragraph) is load-bearing yet lacks an explicit derivation from the Lindblad master equation. The global balance of all dissipators generally yields system occupations that are linear combinations of all reservoir occupancies weighted by the full rate matrix; without showing that the resulting energy flows to each reservoir contain no cross-terms or renormalization factors that depend on the total number of reservoirs, the exact encoding required for the coprocessor remains unverified.
- [Performance estimate paragraph] The performance estimate of ~10^11 operations per second (abstract) assumes the NESS transition time is strictly independent of reservoir number and that technological limitations do not reintroduce dimension dependence. A concrete scaling argument or numerical example for increasing reservoir count is needed to substantiate this independence.
minor comments (2)
- [Introduction/mapping section] Notation for the rate matrix and the precise definition of 'stochastic matrix' in the quantum setting should be introduced explicitly to avoid ambiguity with classical probability matrices.
- [Discussion] A brief comparison to existing analog crossbar implementations or other open-system proposals for linear algebra would help situate the novelty of the thermodynamic encoding.
Simulated Author's Rebuttal
We thank the referee for their careful reading of our manuscript and for the constructive comments. Below we respond to each major comment and indicate the changes made to the manuscript.
read point-by-point responses
-
Referee: The central claim that stationary energy flows exactly equal the matrix-vector product for arbitrary stochastic matrices (abstract and mapping paragraph) is load-bearing yet lacks an explicit derivation from the Lindblad master equation. The global balance of all dissipators generally yields system occupations that are linear combinations of all reservoir occupancies weighted by the full rate matrix; without showing that the resulting energy flows to each reservoir contain no cross-terms or renormalization factors that depend on the total number of reservoirs, the exact encoding required for the coprocessor remains unverified.
Authors: We thank the referee for highlighting the need for greater explicitness. The original manuscript derives the mapping from the steady-state solution of the Lindblad master equation for the chosen open-system couplings, but we agree the presentation can be strengthened. In the revised manuscript we have added a dedicated derivation subsection that starts from the Lindblad equation, solves the linear system for the non-equilibrium steady-state occupations, and substitutes into the energy-flow expressions. Because the rate matrix is stochastic (column sums equal to unity), the resulting energy flows to each reservoir recover exactly the matrix-vector product with no residual cross-terms or renormalization factors that scale with the number of reservoirs. The explicit algebra is now shown in full. revision: yes
-
Referee: The performance estimate of ~10^11 operations per second (abstract) assumes the NESS transition time is strictly independent of reservoir number and that technological limitations do not reintroduce dimension dependence. A concrete scaling argument or numerical example for increasing reservoir count is needed to substantiate this independence.
Authors: We agree that a concrete demonstration strengthens the claim. The revised manuscript now contains both an analytical scaling argument and numerical simulations. Analytically, the eigenvalues of the Liouvillian that govern relaxation are set by the local dissipation rates of each bosonic mode and remain independent of the total number of reservoirs. Numerically, we present results for reservoir counts ranging from 4 to 128; the time to reach within 1 % of the NESS is constant to within numerical precision. We also discuss how technological overhead (control lines, readout) can be kept from reintroducing dimension dependence through modular, parallel fabrication. revision: yes
Circularity Check
No circularity: derivation follows from standard open-system master equation without reduction to inputs or self-citations
full rationale
The paper starts from the standard Lindblad master equation for an open quantum system of bosonic modes coupled to multiple bosonic reservoirs, solves for the non-equilibrium stationary state (NESS), and identifies stationary energy flows with the result of a vector-matrix multiplication. The explicit mapping equates dissipation rates times mode frequencies to conductivities, reservoir occupancies to potentials, and energy flows to currents, which is a constructed analogy rather than a definitional identity. Calculation time is the relaxation time to NESS, a standard Markovian property independent of reservoir number by construction of the dynamics. No fitted parameters are renamed as predictions, no load-bearing self-citations justify uniqueness, and no ansatz is smuggled in; the central claim rests on solving the global balance equations of the master equation, which remains externally verifiable against the Lindblad formalism.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption An open quantum system of bosonic modes coupled to bosonic reservoirs reaches a unique non-equilibrium stationary state whose energy flows encode linear operations on the reservoir occupancies.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
dissipation rates multiplied by OQS’s modes frequencies can be seen as conductivities, reservoirs’ occupancies as potentials, and stationary energy flows as electric currents
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
T. Almutiri, F. Nadeem,et al., Markov models appli- cations in natural language processing: a survey, Int. J. Inf. Technol. Comput. Sci2, 1 (2022)
work page 2022
- [2]
-
[3]
W. Dai, Y. Yu, Y. Dai, and B. Deng, Text steganogra- phy system using markov chain source model and des algorithm., J. Softw.5, 785 (2010)
work page 2010
-
[4]
O. Sigaud and O. Buffet,Markov decision processes in artificial intelligence(John Wiley & Sons, 2013)
work page 2013
-
[5]
N. B¨ auerle and U. Rieder,Markov decision processes with applications to finance(Springer Science & Busi- ness Media, 2011)
work page 2011
- [6]
-
[7]
B. An, K. M. Sim, C. Y. Miao, and Z. Q. Shen, Deci- sion making of negotiation agents using markov chains, Multiagent and Grid Systems4, 5 (2008)
work page 2008
-
[8]
P. P. Bocharov, C. D’Apice, and A. Pechinkin,Queueing theory(Walter de Gruyter, 2011)
work page 2011
-
[9]
N. G¨ une¸ s and M. Yılmaz, Traveling salesman problem as a markov process and distribution of solutions: Case study, in2024 32nd Signal Processing and Communica- tions Applications Conference (SIU)(IEEE, 2024) pp. 1–4
work page 2024
- [10]
- [11]
-
[12]
S. Salman and S. Alaswad, Alleviating road network congestion: Traffic pattern optimization using markov chain traffic assignment, Computers & Operations Re- search99, 191 (2018)
work page 2018
- [13]
- [14]
- [15]
-
[16]
H. Yang, Y. Li, L. Lu, and R. Qi, First order multivari- ate markov chain model for generating annual weather data for hong kong, Energy and Buildings43, 2371 (2011)
work page 2011
- [17]
-
[18]
T. Phelan and K. Eslami, Applications of markov chain approximation methods to optimal control problems in economics, Journal of Economic Dynamics and Control 143, 104437 (2022)
work page 2022
-
[19]
O. Kostoska, V. Stojkoski, and L. Kocarev, On the structure of the world economy: An absorbing markov chain approach, Entropy22, 482 (2020)
work page 2020
-
[20]
V. Moosavi and G. Isacchini, A markovian model of evolving world input-output network, PloS one12, e0186746 (2017)
work page 2017
-
[21]
J. B. Feldman and H. Topaloglu, Revenue management under the markov chain choice model, Operations Re- search65, 1322 (2017)
work page 2017
-
[22]
J. Dong, A. S. Simsek, and H. Topaloglu, Pricing prob- lems under the markov chain choice model, Production and Operations Management28, 157 (2019)
work page 2019
-
[23]
M. S. Fallahnezhad, A. Ranjbar, and F. Z. Saredorahi, A markov model for production and maintenance decision, 6 Macro Management & Public Policies2(2020)
work page 2020
-
[24]
V. Ezugwu and S. Ologun, Markov chain: A predictive model for manpower planning, Journal of Applied Sci- ences and Environmental Management21, 557 (2017)
work page 2017
-
[25]
M. O. Jacksonet al.,Social and economic networks, Vol. 3 (Princeton university press Princeton, 2008)
work page 2008
-
[26]
S. Banisch, R. Lima, and T. Ara´ ujo, Agent based mod- els and opinion dynamics as markov chains, Social Net- works34, 549 (2012)
work page 2012
-
[27]
Q. Zhou, J. Li, D. Wu, and X. L. Xu, A novel risk communication model for online public opinion dissem- ination that integrates the sir and markov chains, Risk Analysis (2025)
work page 2025
-
[28]
Jawandhiya, Hardware design for machine learning, Int
P. Jawandhiya, Hardware design for machine learning, Int. J. Artif. Intell. Appl9, 63 (2018)
work page 2018
-
[29]
K. Berggren, Q. Xia, K. K. Likharev, D. B. Strukov, H. Jiang, T. Mikolajick, D. Querlioz, M. Salinga, J. R. Erickson, S. Pi,et al., Roadmap on emerging hardware and technology for machine learning, Nanotechnology 32, 012002 (2020)
work page 2020
-
[30]
E. Buber and D. Banu, Performance analysis and cpu vs gpu comparison for deep learning, in2018 6th Interna- tional Conference on Control Engineering & Informa- tion Technology (CEIT)(IEEE, 2018) pp. 1–6
work page 2018
-
[31]
D. Melanson, M. Abu Khater, M. Aifer, K. Donatella, M. Hunter Gordon, T. Ahle, G. Crooks, A. J. Martinez, F. Sbahi, and P. J. Coles, Thermodynamic computing system for ai applications, Nature Communications16, 3757 (2025)
work page 2025
-
[32]
A. Wali and S. Das, Two-dimensional memtransistors for non-von neumann computing: Progress and chal- lenges, Advanced Functional Materials34, 2308129 (2024)
work page 2024
-
[33]
H. Seok, D. Lee, S. Son, H. Choi, G. Kim, and T. Kim, Beyond von neumann architecture: Brain-inspired arti- ficial neuromorphic devices and integrated computing, Advanced Electronic Materials , 2300839 (2024)
work page 2024
-
[34]
H. Bhaskaran and W. Pernice,Phase Change Materials- Based Photonic Computing(Elsevier, 2024)
work page 2024
-
[35]
Suri,Applications of Emerging Memory Technology (Springer, 2020)
M. Suri,Applications of Emerging Memory Technology (Springer, 2020)
work page 2020
- [36]
-
[37]
C. de Souza Dias and P. F. Butzen, Memristors: A jour- ney from material engineering to beyond von-neumann computing, Journal of Integrated Circuits and Systems 16, 1 (2021)
work page 2021
- [38]
-
[39]
G. Cerofolini and G. Cerofolini, The crossbar structure, Nanoscale Devices: Fabrication, Functionalization, and Accessibility from the Macroscopic World , 45 (2009)
work page 2009
-
[40]
H. Li, S. Wang, X. Zhang, W. Wang, R. Yang, Z. Sun, W. Feng, P. Lin, Z. Wang, L. Sun,et al., Memristive crossbar arrays for storage and computing applications, Advanced Intelligent Systems3, 2100017 (2021)
work page 2021
-
[41]
K.-H. Kim, S. Gaba, D. Wheeler, J. M. Cruz-Albrecht, T. Hussain, N. Srinivasa, and W. Lu, A functional hy- brid memristor crossbar-array/cmos system for data storage and neuromorphic applications, Nano Letters 12, 389 (2012)
work page 2012
-
[42]
U. R¨ uhrmair, C. Jaeger, M. Bator, M. Stutzmann, P. Lugli, and G. Csaba, Applications of high-capacity crossbar memories in cryptography, IEEE Transactions on Nanotechnology10, 489 (2010)
work page 2010
-
[43]
Y. Chen, G.-Y. Jung, D. A. Ohlberg, X. Li, D. R. Stew- art, J. O. Jeppesen, K. A. Nielsen, J. F. Stoddart, and R. S. Williams, Nanoscale molecular-switch crossbar cir- cuits, Nanotechnology14, 462 (2003)
work page 2003
-
[44]
S. Garg, J. Lou, A. Jain, Z. Guo, B. J. Shastri, and M. Nahmias, Dynamic precision analog computing for neural networks, IEEE Journal of Selected Topics in Quantum Electronics29, 1 (2023)
work page 2023
-
[45]
A. Adamatzky and L. Chua,Memristor networks (Springer Science & Business Media, 2013)
work page 2013
-
[46]
S. V. Stasenko, A. N. Mikhaylov, A. A. Fedotov, V. A. Smirnov, and V. B. Kazantsev, Astrocyte control burst- ing mode of spiking neuron network with memristor- implemented plasticity, Chaos, Solitons & Fractals181, 114648 (2024)
work page 2024
-
[47]
S. Shchanikov, I. Bordanov, A. Kucherik, E. Gryaznov, and A. Mikhaylov, Neuromorphic analog machine vision enabled by nanoelectronic memristive devices, Applied Sciences13, 13309 (2023)
work page 2023
-
[48]
M. Mishchenko, D. Bolshakov, V. Lukoyanov, D. Ko- rolev, A. Belov, D. Guseinov, V. Matrosov, V. Kazant- sev, and A. Mikhaylov, Inverted spike-rate-dependent plasticity due to charge traps in a metal-oxide memris- tive device, Journal of Physics D: Applied Physics55, 394002 (2022)
work page 2022
-
[49]
M. Spagnolo, J. Morris, S. Piacentini, M. Antesberger, F. Massa, A. Crespi, F. Ceccarelli, R. Osellame, and P. Walther, Experimental photonic quantum memristor, Nature Photonics16, 318 (2022)
work page 2022
- [50]
- [51]
-
[52]
N. Youngblood, C. A. R´ ıos Ocampo, W. H. Pernice, and H. Bhaskaran, Integrated optical memristors, Na- ture Photonics17, 561 (2023)
work page 2023
-
[53]
J.-Y. Mao, L. Zhou, X. Zhu, Y. Zhou, and S.-T. Han, Photonic memristor for future computing: a perspec- tive, Advanced Optical Materials7, 1900766 (2019)
work page 2019
-
[54]
S. Shrivastava, L. B. Keong, S. Pratik, A. S. Lin, and T.-Y. Tseng, Fully photon controlled synaptic memris- tor for neuro-inspired computing, Advanced Electronic Materials9, 2201093 (2023)
work page 2023
-
[55]
C. Lian, C. Vagionas, T. Alexoudi, N. Pleros, N. Young- blood, and C. R´ ıos, Photonic (computational) memo- ries: tunable nanophotonics for data storage and com- puting, Nanophotonics11, 3823 (2022)
work page 2022
-
[56]
L. Carroll, J.-S. Lee, C. Scarcella, K. Gradkowski, M. Duperron, H. Lu, Y. Zhao, C. Eason, P. Morris- sey, M. Rensing,et al., Photonic packaging: transform- ing silicon photonic integrated circuits into photonic de- vices, Applied Sciences6, 426 (2016)
work page 2016
-
[57]
M. Van Rossum, W. Schoenmaker, W. Magnus, K. De Meyer, M. D. Croitoru, V. N. Gladilin, V. M. Fomin, and J. T. Devreese, Moore’s law: new play- ground for quantum physics, Physica Status Solidi (b) 7 237, 426 (2003)
work page 2003
- [58]
-
[59]
J. N. Randall, J. H. Owen, E. Fuchs, J. Lake, J. R. Von Ehr, J. Ballard, and E. Henriksen, Digital atomic scale fabrication an inverse moore’s law–a path to atom- ically precise manufacturing, Micro and Nano Engineer- ing1, 1 (2018)
work page 2018
-
[60]
A. L. Yeyati, D. Subero, J. Pekola, and R. S´ anchez, Pho- tonic heat transport through a josephson junction in a resistive environment, Physical Review B110, L220502 (2024)
work page 2024
- [61]
-
[62]
Y. V. Nazarov and Y. M. Blanter,Quantum trans- port: introduction to nanoscience(Cambridge Univer- sity Press, 2009)
work page 2009
-
[63]
H. Haug, A.-P. Jauho,et al.,Quantum kinetics in trans- port and optics of semiconductors, Vol. 2 (Springer, 2008)
work page 2008
-
[64]
H.-P. Breuer and F. Petruccione,The theory of open quantum systems(Oxford University Press, USA, 2002)
work page 2002
-
[65]
I. Vovchenko, A. Zyablovsky, A. Pukhov, and E. An- drianov, Transient temperature dynamics of reservoirs connected through an open quantum system, Physical Review E109, 044144 (2024)
work page 2024
- [66]
- [67]
-
[68]
A. Levy and R. Kosloff, The local approach to quantum transport may violate the second law of thermodynam- ics, Europhysics Letters107, 20004 (2014)
work page 2014
-
[69]
A. Trushechkin and I. Volovich, Perturbative treat- ment of inter-site couplings in the local description of open quantum networks, Europhysics Letters113, 30005 (2016)
work page 2016
-
[70]
I. Vovcenko, A. Zyablovsky, A. Pukhov, and E. An- drianov, Energy transport induced by transition from the weak to the strong coupling regime between non- hermitian optical systems, Journal of the Optical Soci- ety of America B40, 2990 (2023)
work page 2023
-
[71]
P. P. Potts, A. A. S. Kalaee, and A. Wacker, A thermo- dynamically consistent markovian master equation be- yond the secular approximation, New Journal of Physics 23, 123013 (2021)
work page 2021
-
[72]
H. Carmichael,An open systems approach to quantum optics: lectures presented at the Universit´ e Libre de Bruxelles, October 28 to November 4, 1991, Vol. 18 (Springer Science & Business Media, 2009)
work page 1991
- [73]
-
[74]
S. Vinjanampathy and J. Anders, Quantum thermody- namics, Contemporary Physics57, 545 (2016)
work page 2016
-
[75]
V. Y. Shishkov, E. S. Andrianov, A. A. Pukhov, A. P. Vinogradov, and A. A. Lisyansky, Relaxation of inter- acting open quantum systems, Physics-Uspekhi62, 510 (2019)
work page 2019
-
[76]
I. Vovchenko, A. Zyablovsky, A. Pukhov, and E. Andri- anov, Autonomous coarse-grained cooling of the coldest reservoir not restricted by the second law of thermody- namics, Physical Review A111, 062204 (2025)
work page 2025
- [77]
-
[78]
P. Lipka-Bartosik, M. Perarnau-Llobet, and N. Brun- ner, Thermodynamic computing via autonomous quan- tum thermal machines, Science Advances10, eadm8792 (2024)
work page 2024
- [79]
-
[80]
P. P. Hofer, M. Perarnau-Llobet, L. D. M. Miranda, G. Haack, R. Silva, J. B. Brask, and N. Brunner, Marko- vian master equations for quantum thermal machines: local versus global approach, New Journal of Physics 19, 123037 (2017)
work page 2017
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.