Scientists have developed a novel methodology for calculating Green-Kubo transport coefficients, crucial parameters in molecular simulations, by reformulating classical molecular dynamics as a problem amenable to quantum algorithms. Masari Watanabe at Quemix Inc, and colleagues at DENSO CORP, The University of Tokyo, Quantum Materials and Applications Research Centre, and National Institute, demonstrate that both NVE (Nominal Velocity and Energy conserving) and Nosé-Hoover dynamics, commonly used for simulating molecular systems, can be expressed as unitary evolutions within quantum Hilbert spaces. This innovative framework promises a substantial improvement in the efficiency of transport-coefficient calculations, achieving exponential error reduction with increasing qubit numbers and enhancing the statistical estimation of key probabilities from the typical inverse-square-root scaling to a near-linear scaling with the number of queries. Their detailed analysis, including estimations of the circuit resources required, constitutes a concrete step forward in the application of quantum algorithms to practical molecular simulation, addressing a long-standing need for more efficient computational methods in materials science.
Enhanced Green-Kubo calculations via quantum algorithms and optimised statistical estimation
The statistical estimation of the probability P_0, central to the accurate calculation of Green-Kubo transport coefficients, now improves from a scaling of $N$ queries to the power of -1/2, to a scaling approaching $N$ queries to the power of -1. This represents a significant advancement. Achieving this level of precision was previously computationally prohibitive for complex systems, exceeding the limitations of direct shot sampling, a conventional method which previously hindered accurate calculations of material transport properties such as thermal conductivity, viscosity, and diffusion coefficients. Quantum algorithms facilitate this enhanced efficiency by formulating classical molecular dynamics as a quantum readout problem, leveraging the Koopman-von Neumann (KvN) representation. The KvN representation maps classical observables to quantum operators, allowing for the application of quantum techniques to classical systems. Discretization error in the correlation function, a critical factor for accurate transport property calculations, diminishes as a power law dependent on the number of grid points used. Specifically, employing Nz=2nz grid points results in an exponential decrease in error with register size, nz. This implies that achieving a target accuracy of ε requires only $\mathcal{O}(\log(1/ε))$ qubits, representing a sharp reduction in computational demand compared to classical methods which often scale polynomially with the desired accuracy. The Green-Kubo framework relies on calculating time correlation functions of relevant fluxes, and minimising the error in these functions is paramount for obtaining reliable results.
A single step of the NVE propagator, which simulates the time evolution of particle positions and momenta, can be constructed using mathcalO(n2) CX gates, where ‘n’ represents the total number of position and momentum qubits. The CX gate, or controlled-NOT gate, is a fundamental building block in quantum circuits. The Nosé-Hoover propagator, incorporating a thermostat to control temperature and maintain a constant NVT (constant number of particles, volume, and temperature) ensemble, scales as mathcalO(n_ξnp,2np), dependent on the number of momentum (np) and thermostat (n_ξ) qubits. The Nosé-Hoover thermostat introduces additional degrees of freedom to the system, necessitating a larger quantum circuit. Accurate modelling of molecular behaviour is vital for calculating material properties, a task currently limited by the computational cost of simulating complex systems, particularly those with many interacting particles or long timescales. Classical simulations often struggle with the exponential growth of computational resources required to accurately represent the phase space of these systems.
These propagators demonstrate the potential for quantum speedups in molecular dynamics simulations, offering a pathway to overcome current computational bottlenecks. The ability to represent classical dynamics as quantum evolutions opens up possibilities for utilising quantum algorithms such as quantum phase estimation and variational quantum eigensolvers to accelerate the calculation of transport properties. Although the framework’s feasibility is demonstrated on finite grids, a simplification that doesn’t fully capture the intricacies of real-world molecular simulations with continuous variables, it represents a valuable advance in computational materials science. The use of finite grids introduces discretisation errors, but the analysis shows these errors decrease predictably with grid resolution. The authors acknowledge that current quantum hardware remains in its infancy, with limitations in qubit count, coherence times, and gate fidelity, raising questions regarding the practical limits of this approach on near-term quantum hardware. Further research is needed to optimise the quantum circuits and mitigate the effects of noise. Establishing a clear pathway for utilising future, more powerful quantum processors is crucial for realising the full potential of this methodology. Green-Kubo transport coefficients, essential for modelling material behaviours and predicting macroscopic properties from microscopic interactions, are now translated into a quantum computing framework. Classical molecular dynamics can be expressed as quantum evolutions by employing this mathematical technique to simplify complex systems, allowing for the potential application of quantum algorithms to enhance computational efficiency and unlock new possibilities in materials discovery and design. The development of efficient quantum algorithms for calculating transport properties could revolutionise fields such as materials science, chemistry, and engineering.
The researchers successfully reformulated Green, Kubo transport coefficients, used to model material behaviours, as a problem solvable by quantum algorithms. This representation allows classical molecular dynamics to be expressed as quantum evolutions, potentially enabling the use of quantum computation to improve efficiency. Numerical tests demonstrated that discretisation errors decrease exponentially with increasing grid resolution, requiring a number of qubits proportional to the logarithm of desired accuracy. The authors suggest future work will focus on optimising quantum circuits and addressing the limitations of current quantum hardware to fully realise this methodology.
👉 More information
🗞 Koopman–von Neumann Molecular Dynamics for Green–Kubo Transport Coefficients
🧠 ArXiv: https://arxiv.org/abs/2605.30142
