Quantum computing for classical problems: Variational Quantum Eigensolver for activated processes

The theory of stochastic processes impacts both physical and social sciences. At the molecular scale, stochastic dynamics is ubiquitous because of thermal fluctuations. The Fokker-Plank-Smoluchowski equation models the time evolution of the probability density of selected degrees of freedom in the diffusive regime and it is therefore a workhorse of physical chemistry. In this paper we report the development and implementation of a Variational Quantum Eigensolver procedure to solve the Fokker-Planck-Smoluchowski eigenvalue problem. We show that such an algorithm, typically adopted to address quantum chemistry problems, can be applied effectively to classical systems paving the way to new applications of quantum computers. We compute the conformational transition rate in a linear chain of rotors experiencing nearest-neighbour interaction. We provide a method to encode on the quantum computer the probability distribution for a given conformation of the chain and assess its scalability in terms of operations. Performance analysis on no isy quantum emulators and quantum devices (IBMQ Santiago) is provided for a small chain showing results in good agreement with the classical benchmark without further addition of any error mitigation technique.