Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2511.22898 · Physical Review Letters

Algorithmic Quantum Simulations of Quantum Thermodynamics

Yangsen Ye, Jue Nan, Dong Chen, Torsten V. Zache, Qingling Zhu, Yiming Zhang, Yuan Li, Xiawei Chen, Chong Ying, Chen Zha, Sirui Cao, Shaowei Li, Shaojun Guo, Haoran Qian, Hao Rong, Yulin Wu, Kai Yan, Feifan Su, Hui Deng, Yu Xu, Jin Lin, Ming Gong, Fusheng Chen, Gang Wu, Yong-Heng Huo, Chao-Yang Lu, Cheng-Zhi Peng, Xiaobo Zhu, Xiaopeng Li, Jian-Wei Pan

quant-phSubmitted Nov 28, 2025• Announced Dec 1, 2025• 0 citations
arXivHTMLPDF

Abstract

Characterizing quantum phases-of-matter at finite-temperature is essential for understanding complex materials and large-scale thermodynamic phenomena. Here, we develop algorithmic protocols for simulating quantum thermodynamics on quantum hardware through quantum kernel function expansion (QKFE), producing the free energy as an analytic function of temperature with uniform convergence. These protocols are demonstrated by simulating transverse field Ising and XY models with superconducting qubits. In both analogue and digital implementations of the QKFE algorithms, we exhibit quantitative agreement of our quantum simulation experiments with the exact results. Our approach provides a general framework for computing thermodynamic potentials on programmable quantum devices, granting access to key thermodynamic properties such as entropy, heat capacity and criticality, with far-reaching implications for material design and drug development.

Comments
31 pages, 19 figures, 7 tables

Related papers

Kernel function based quantum algorithms for finite temperature quantum simulation

Hai Wang, Jue Nan, Tao Zhang +3 more

Computing finite temperature properties of a quantum many-body system is key to describing a broad range of correlated quantum many-body physics from quantum chemistry and condensed matter to thermal quantum field theories. Quantum computing, which has seen rapid developments in recent years, has a huge potential to impact the computation of quantum thermodynamics. To fulfill these potential impacts, it is crucial to design quantum algorithms that utilize the computation power of quantum computing devices. Here, we present a quantum kernel function expansion (QKFE) algorithm for solving thermodynamic properties of quantum many-body systems. In this quantum algorithm, the many-body density of states is approximated by a kernel Fourier expansion, whose expansion moments are obtained by random state sampling and quantum interferometric measurements. Compared to its classical counterpart, namely, the kernel polynomial method (KPM), QKFE has an exponential advantage in the cost of both time and memory. In computing low temperature properties, QKFE becomes inefficient, similar to classical KPM. To resolve this difficulty, we further construct a thermal ensemble iteration (THEI) protocol which starts from the trivial limit of an infinite temperature ensemble and approaches the low temperature regime step by step. For quantum Hamiltonians, whose ground states are preparable with polynomial quantum circuits, THEI has an overall polynomial complexity. We demonstrate its efficiency with applications to one- and two-dimensional quantum spin models and a fermionic lattice. With our analysis of the realization with digital and analog quantum devices, we believe the quantum algorithm is accessible to current quantum technology.

arXiv:2202.01170•Physical review. B./Physical review. B•Aug 1, 2023•4 citations

Inverse Quantum Simulation for Quantum Material Design

Christian Kokail, Pavel E. Dolgirev, Rick van Bijnen +3 more

Quantum simulation provides a powerful route for exploring many-body phenomena beyond the capabilities of classical computation. Existing approaches typically proceed in the forward direction: a model Hamiltonian is specified, implemented on a programmable quantum platform, and its phase diagram and properties are explored. Here we present a quantum algorithmic framework for inverse quantum simulation, enabling quantum material design with desired properties. Target material characteristics are encoded as a cost function, which is minimized on quantum hardware to prepare a many-body state with the desired properties in quantum memory. Hamiltonian learning is then used to reconstruct a low-energy Hamiltonian for which this state is an approximate ground state, yielding a physically interpretable model that can guide experimental synthesis. As illustrative applications, we outline how the method can be used to search for high-temperature superconductors within the fermionic Hubbard model, enhancing $d$-wave correlations over a broad range of dopings and temperatures, design quantum phases by stabilizing a topological order through continuous Hamiltonian modifications, and optimize dynamical properties relevant for photochemistry and frequency- and momentum-resolved condensed-matter data. These results extend the scope of quantum simulators from exploring quantum many-body systems to designing and discovering new quantum materials.

quant-pharXiv:2601.12239•Jan 20, 2026•0 citations

Quantum simulation of thermodynamics: Maxwell relations for pair correlations

Rist, F., Raymon S. Watson, H. L. Nourse +2 more

Quantum simulators hold enormous promise for advancing the modelling of materials and understanding emergent physics, such as high temperature superconductivity and topological order. While correlation functions are, typically, straightforward to measure in quantum simulators, thermodynamic properties are not. This limits our ability to directly compare the results of quantum simulations to experiments on the materials being modelled. Maxwell relations are an extremely powerful tool for characterising complex materials, as they enable the determination of challenging-to-measure thermodynamic properties from more accessible ones. Here, we introduce generalised Maxwell relations that relate every thermodynamic quantity to a single local correlation function. We illustrate their utility by deducing the thermodynamic properties of several iconic quantum many-body models from pair correlation functions using the generalised Maxwell relations. We show that this {universal} approach is readily accessible in quantum simulators and suggest applications to condensed matter systems where thermodynamic measurements are challenging, such as atomically thin materials.

arXiv:2506.19407•Jun 24, 2025•0 citations

Quantum Phase Transitions in the Transverse-Field Ising Model: A Comparative Study of Exact, Variational, and Hardware-Based Approaches

Rudraksh Sharma

The quantum phase transitions provide a paradigm for studying collective quantum phenomena that are a result of competing non-commuting interactions. This paper will study the ground state properties and quantum critical dynamics of the one-dimensional transverse field Ising model through a combined perspective that includes exact diagonalisation, variational quantum eigensolver (VQE) simulations, and simulations on realistic physical quantum devices. We focus on a lattice of four spins, where we calculate the ground-state energies, magnetic order parameters and correlation functions at uniformly applied conditions, which is repeated by all systems. Precise diagonalisation provides both a benchmark, which is symmetry-conserving, and a depth-two, physics inspired variational approximation, which provides simulations accessible to hardware. The circuits that have been optimised identically are then placed on the IQM Garnet quantum processor, using a resource-efficient batched protocol. We find that the ground-state energies of shallow variational circuits are reliably captured by the circuit over the entire parameter space; the magnetic arrangement parameters and observables sensitive to correlation signal significantly more noise. The error analysis of quantitative analysis reveals a strong broadening of critical crossover on hardware, which is consistent with the noise attenuation of long-range correlations. These findings highlight the current capabilities as well as the fundamental limitations of noisy intermediate-scale quantum systems in modelling quantum critical phenomena as a benchmark to future enhancements in obtaining quantum hardware and quantum algorithms development.

quant-pharXiv:2601.17515•Jan 27, 2026•0 citations

Calculating the many-body density of states on a digital quantum computer

Alessandro Summer, Cecilia Chiaracane, Mark T. Mitchison +1 more

Quantum statistical mechanics allows us to extract thermodynamic information from a microscopic description of a many-body system. A key step is the calculation of the density of states, from which the partition function and all finite-temperature equilibrium thermodynamic quantities can be calculated. In this work, we devise and implement a quantum algorithm to perform an estimation of the density of states on a digital quantum computer which is inspired by the kernel polynomial method. Classically, the kernel polynomial method allows us to sample spectral functions via a Chebyshev polynomial expansion. Our algorithm computes moments of the expansion on quantum hardware using a combination of random-state preparation for stochastic trace evaluation and a controlled unitary operator. We use our algorithm to estimate the density of states of a nonintegrable Hamiltonian on the Quantinuum H1-1 trapped ion chip for a controlled register of 18 qubits. This not only represents a state-of-the-art calculation of thermal properties of a many-body system on quantum hardware, but also exploits the controlled unitary evolution of a many-qubit register on an unprecedented scale. Published by the American Physical Society 2024

arXiv:2303.13476•Physical Review Research•Jan 26, 2024•11 citations

Finite-Temperature Dynamical Phase Diagram of the $2+1$D Quantum Ising Model

Lucas Katschke, Roland C. Farrell, Umberto Borla +2 more

Mapping finite-temperature dynamical phase diagrams of quantum many-body models is a necessary step towards establishing a framework of far-from-equilibrium quantum many-body universality. However, this is quite difficult due, in part, to the severe challenges in representing the volume-law entanglement that is generated under nonequilibrium dynamics at finite temperatures. Here, we address these challenges with an efficient equilibrium quantum Monte Carlo (QMC) framework for computing the finite-temperature dynamical phase diagram. Our method uses energy conservation and the self-thermalizing properties of ergodic quantum systems to determine observables at late times after a quantum quench. We use this technique to chart the dynamical phase diagram of the $2+1$D quantum Ising model generated by quenches of the transverse field in initial thermal states. Our approach allows us to track the evolution of dynamical phases as a function of both the initial temperature and transverse field. Surprisingly, we identify quenches in the ordered phase that cool the system as well as an interval of initial temperatures where it is possible to quench from the paramagnetic (PM) to ferromagnetic (FM) phases. Our method gives access to dynamical properties without explicitly simulating unitary time evolution, and is immediately applicable to other lattice geometries and interacting many-body systems. Finally, we propose a quantum simulation experiment on state-of-the-art digital quantum hardware to directly probe the predicted dynamical phases and their real-time formation.

quant-pharXiv:2602.16772•Feb 20, 2026•0 citations
Characterizing quantum phases-of-matter at finite-temperature is essential for understanding complex materials and large-scale thermodynamic phenomena. Here, we develop algorithmic protocols for simulating quantum thermodynamics on quantum hardware through quantum kernel function expansion (QKFE), producing the free energy as an analytic function of temperature with uniform convergence. These protocols are demonstrated by simulating transverse field Ising and XY models with superconducting qubits. In both analogue and digital implementations of the QKFE algorithms, we exhibit quantitative agreement of our quantum simulation experiments with the exact results. Our approach provides a general framework for computing thermodynamic potentials on programmable quantum devices, granting access to key thermodynamic properties such as entropy, heat capacity and criticality, with far-reaching implications for material design and drug development.
Computing finite temperature properties of a quantum many-body system is key to describing a broad range of correlated quantum many-body physics from quantum chemistry and condensed matter to thermal quantum field theories. Quantum computing, which has seen rapid developments in recent years, has a huge potential to impact the computation of quantum thermodynamics. To fulfill these potential impacts, it is crucial to design quantum algorithms that utilize the computation power of quantum computing devices. Here, we present a quantum kernel function expansion (QKFE) algorithm for solving thermodynamic properties of quantum many-body systems. In this quantum algorithm, the many-body density of states is approximated by a kernel Fourier expansion, whose expansion moments are obtained by random state sampling and quantum interferometric measurements. Compared to its classical counterpart, namely, the kernel polynomial method (KPM), QKFE has an exponential advantage in the cost of both time and memory. In computing low temperature properties, QKFE becomes inefficient, similar to classical KPM. To resolve this difficulty, we further construct a thermal ensemble iteration (THEI) protocol which starts from the trivial limit of an infinite temperature ensemble and approaches the low temperature regime step by step. For quantum Hamiltonians, whose ground states are preparable with polynomial quantum circuits, THEI has an overall polynomial complexity. We demonstrate its efficiency with applications to one- and two-dimensional quantum spin models and a fermionic lattice. With our analysis of the realization with digital and analog quantum devices, we believe the quantum algorithm is accessible to current quantum technology.
Quantum simulation provides a powerful route for exploring many-body phenomena beyond the capabilities of classical computation. Existing approaches typically proceed in the forward direction: a model Hamiltonian is specified, implemented on a programmable quantum platform, and its phase diagram and properties are explored. Here we present a quantum algorithmic framework for inverse quantum simulation, enabling quantum material design with desired properties. Target material characteristics are encoded as a cost function, which is minimized on quantum hardware to prepare a many-body state with the desired properties in quantum memory. Hamiltonian learning is then used to reconstruct a low-energy Hamiltonian for which this state is an approximate ground state, yielding a physically interpretable model that can guide experimental synthesis. As illustrative applications, we outline how the method can be used to search for high-temperature superconductors within the fermionic Hubbard model, enhancing ddd-wave correlations over a broad range of dopings and temperatures, design quantum phases by stabilizing a topological order through continuous Hamiltonian modifications, and optimize dynamical properties relevant for photochemistry and frequency- and momentum-resolved condensed-matter data. These results extend the scope of quantum simulators from exploring quantum many-body systems to designing and discovering new quantum materials.
Quantum simulators hold enormous promise for advancing the modelling of materials and understanding emergent physics, such as high temperature superconductivity and topological order. While correlation functions are, typically, straightforward to measure in quantum simulators, thermodynamic properties are not. This limits our ability to directly compare the results of quantum simulations to experiments on the materials being modelled. Maxwell relations are an extremely powerful tool for characterising complex materials, as they enable the determination of challenging-to-measure thermodynamic properties from more accessible ones. Here, we introduce generalised Maxwell relations that relate every thermodynamic quantity to a single local correlation function. We illustrate their utility by deducing the thermodynamic properties of several iconic quantum many-body models from pair correlation functions using the generalised Maxwell relations. We show that this {universal} approach is readily accessible in quantum simulators and suggest applications to condensed matter systems where thermodynamic measurements are challenging, such as atomically thin materials.
The quantum phase transitions provide a paradigm for studying collective quantum phenomena that are a result of competing non-commuting interactions. This paper will study the ground state properties and quantum critical dynamics of the one-dimensional transverse field Ising model through a combined perspective that includes exact diagonalisation, variational quantum eigensolver (VQE) simulations, and simulations on realistic physical quantum devices. We focus on a lattice of four spins, where we calculate the ground-state energies, magnetic order parameters and correlation functions at uniformly applied conditions, which is repeated by all systems. Precise diagonalisation provides both a benchmark, which is symmetry-conserving, and a depth-two, physics inspired variational approximation, which provides simulations accessible to hardware. The circuits that have been optimised identically are then placed on the IQM Garnet quantum processor, using a resource-efficient batched protocol. We find that the ground-state energies of shallow variational circuits are reliably captured by the circuit over the entire parameter space; the magnetic arrangement parameters and observables sensitive to correlation signal significantly more noise. The error analysis of quantitative analysis reveals a strong broadening of critical crossover on hardware, which is consistent with the noise attenuation of long-range correlations. These findings highlight the current capabilities as well as the fundamental limitations of noisy intermediate-scale quantum systems in modelling quantum critical phenomena as a benchmark to future enhancements in obtaining quantum hardware and quantum algorithms development.
Quantum statistical mechanics allows us to extract thermodynamic information from a microscopic description of a many-body system. A key step is the calculation of the density of states, from which the partition function and all finite-temperature equilibrium thermodynamic quantities can be calculated. In this work, we devise and implement a quantum algorithm to perform an estimation of the density of states on a digital quantum computer which is inspired by the kernel polynomial method. Classically, the kernel polynomial method allows us to sample spectral functions via a Chebyshev polynomial expansion. Our algorithm computes moments of the expansion on quantum hardware using a combination of random-state preparation for stochastic trace evaluation and a controlled unitary operator. We use our algorithm to estimate the density of states of a nonintegrable Hamiltonian on the Quantinuum H1-1 trapped ion chip for a controlled register of 18 qubits. This not only represents a state-of-the-art calculation of thermal properties of a many-body system on quantum hardware, but also exploits the controlled unitary evolution of a many-qubit register on an unprecedented scale. Published by the American Physical Society 2024
Finite-Temperature Dynamical Phase Diagram of the 2+12+12+1D Quantum Ising Model
Mapping finite-temperature dynamical phase diagrams of quantum many-body models is a necessary step towards establishing a framework of far-from-equilibrium quantum many-body universality. However, this is quite difficult due, in part, to the severe challenges in representing the volume-law entanglement that is generated under nonequilibrium dynamics at finite temperatures. Here, we address these challenges with an efficient equilibrium quantum Monte Carlo (QMC) framework for computing the finite-temperature dynamical phase diagram. Our method uses energy conservation and the self-thermalizing properties of ergodic quantum systems to determine observables at late times after a quantum quench. We use this technique to chart the dynamical phase diagram of the 2+12+12+1D quantum Ising model generated by quenches of the transverse field in initial thermal states. Our approach allows us to track the evolution of dynamical phases as a function of both the initial temperature and transverse field. Surprisingly, we identify quenches in the ordered phase that cool the system as well as an interval of initial temperatures where it is possible to quench from the paramagnetic (PM) to ferromagnetic (FM) phases. Our method gives access to dynamical properties without explicitly simulating unitary time evolution, and is immediately applicable to other lattice geometries and interacting many-body systems. Finally, we propose a quantum simulation experiment on state-of-the-art digital quantum hardware to directly probe the predicted dynamical phases and their real-time formation.