ExplorerQuantum ComputingQuantum Physics
Research PaperResearchia:202605.30013

Koopman--von Neumann Molecular Dynamics for Green--Kubo Transport Coefficients

Masari Watanabe

Abstract

We formulate the Green--Kubo transport coefficients of classical molecular dynamics as a readout problem for quantum algorithms using the Koopman--von Neumann (KvN) representation. Both NVE and Nosé--Hoover-type NVT dynamics are derived as unitary evolutions on Hilbert spaces associated with the corresponding classical phase spaces. Numerical benchmarks on finite grids show that the discretization error in the correlation function decreases as a power law in the number of grid points $N_z$. Equi...

Submitted: May 30, 2026Subjects: Quantum Physics; Quantum Computing

Description / Details

We formulate the Green--Kubo transport coefficients of classical molecular dynamics as a readout problem for quantum algorithms using the Koopman--von Neumann (KvN) representation. Both NVE and Nosé--Hoover-type NVT dynamics are derived as unitary evolutions on Hilbert spaces associated with the corresponding classical phase spaces. Numerical benchmarks on finite grids show that the discretization error in the correlation function decreases as a power law in the number of grid points NzN_z. Equivalently, with Nz=2nzN_z=2^{n_z}, the error decreases exponentially in the register size nzn_z, so a target accuracy εε requires nz=O(log(1/ε))n_z=\mathcal{O}(\log(1/ε)) qubits. To read out a transport coefficient, we input a flux-excited state to quantum phase estimation (QPE). The probability P0P_0 of measuring the QPE ancilla register in the all-zero state corresponds to a Bartlett-windowed Green--Kubo integral. With maximum-likelihood amplitude estimation, the statistical estimation of P0P_0 defined by this QPE oracle improves from the Nqueries1/2N_{\rm queries}^{-1/2} scaling of direct shot sampling to scaling close to Nqueries1N_{\rm queries}^{-1}. Our circuit-resource analysis shows that one step of the NVE propagator can be built with O(n2)\mathcal{O}(n^2) CX gates, where n=nx+npn=n_x+n_p is the total number of position and momentum qubits. For the NVT propagator, the centered-difference Pauli-decomposition implementation of the Nosé--Hoover friction term scales as O(nξnp2np)\mathcal{O}(n_ξn_p\,2^{n_p}), where npn_p and nξn_ξ are the numbers of momentum and thermostat qubits, respectively. The proposed framework is a concrete step toward translating the principles of quantum algorithms into the transport-coefficient calculations required in practical molecular simulation.


Source: arXiv:2605.30142v1 - http://arxiv.org/abs/2605.30142v1 PDF: https://arxiv.org/pdf/2605.30142v1 Original Link: http://arxiv.org/abs/2605.30142v1

Please sign in to join the discussion.

No comments yet. Be the first to share your thoughts!

Access Paper
View Source PDF
Submission Info
Date:
May 30, 2026
Topic:
Quantum Computing
Area:
Quantum Physics
Comments:
0
Bookmark