Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Quantum-enhanced Markov chain Monte Carlo sampling to model Lagrangian tracer dispersion in turbulent boundary layer

Fabian Schindler and Jörg Schumacher*

  • *Contact author: joerg.schumacher@tu-ilmenau.de

Phys. Rev. Research 8, 033031 – Published 7 July, 2026

DOI: https://doi.org/10.1103/5qcg-fhk1

Abstract

We present a quantum-enhanced Markov chain Monte Carlo (QE-MCMC) method to sample turbulent acceleration vectors from a joint target distribution that depends on all three components and height to model the transport and dispersion of massless Lagrangian tracer particles in two different turbulent shear flows. A homogeneous shear flow, which is characterized by a uniform shear rate S, is considered as the starting point. Second, a turbulent boundary layer, which forms in both halves of a plane turbulent channel flow (TCF) at a friction Reynolds number of Reτ=1000, is considered. Then, the mean shear rate S(y) varies with distance from the wall y. In this hybrid quantum-classical method, the proposal distribution Q for the first of the two Metropolis-Hastings sampling substeps is constructed by a parametric quantum circuit. The algorithm generates synthetic tracer particle tracks. The resulting scaling laws for the tracer particle pair dispersion, a central quantity to probe the turbulent mixing properties of the shear flow from a Lagrangian perspective, are found to agree with a stochastic transport model consisting of a set of coupled Langevin equations and the classical MCMC counterpart. Differently to the classical sampling method, QE-MCMC takes to this end a tempered target distribution. Due to the height dependence of the tracer dynamics in TCF, an effective height-weighted spectral gap between the first and second eigenvalues of the Markov chain transition matrix is introduced. The latter is found to significantly exceed the one of classical MCMC when the sampling is done from a multivariate distribution with cross-correlations at the highest qubit numbers and thus resolutions. Consequently, our results support the applicability of this one-shot algorithm as a generative Lagrangian quantum computing module, possibly embedded in a complex fluid flow problem. Our module is found to work reliably for a relatively small number of qubits per spatial dimension of Nq≤6.

View figure in article

Physics Subject Headings (PhySH)

Article Text

References (74)

  1. C. M. Bishop, Pattern Recognition and Machine Learning (Springer, New York, 2006).
  2. N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, Equation of state calculations by fast computing machines, J. Chem. Phys. 21, 1087 (1953).
  3. W. K. Hastings, Monte Carlo sampling methods using Markov chains and their applications, Biometrika 57, 97 (1970).
  4. J. A. Vrugt, C. J. F. ter Braak, C. G. H. Diks, and G. Schoups, Hydrologic data assimilation using particle Markov chain Monte Carlo simulation: Theory, concepts and applications, Adv. Water Resour. 51, 457 (2013).
  5. F. Zuccarello, G. Bilotta, F. Cannavò, A. Cappello, R. Guardo, and G. Ganci, A Markov chain Monte Carlo approach for complex lava flow simulations driven by satellite-derived data, Environ. Model. Software 188, 106426 (2025).
  6. M. Cerezo, A. Arrasmith, R. Babbush, S. C. Benjamin, S. Endo, K. Fujii, J. R. McClean, K. Mitarai, X. Yuan, L. Cincio, and P. J. Coles, Variational quantum algorithms, Nat. Rev. Phys. 3, 625 (2021).
  7. T. L. Patti, O. Shehab, K. Najafi, and S. F. Yelin, Markov chain Monte Carlo enhanced variational quantum algorithms, Quantum Sci. Technol. 8, 015019 (2023).
  8. M. Nielsen and I. Chuang, Quantum Computation and Quantum Information, 10th anniversary ed. (Cambridge University Press, Cambridge, UK, 2011).
  9. F. Arute, K. Arya, and R. Babbush, et al., Quantum supremacy using a programmable superconducting processor, Nature (London) 574, 505 (2019).
  10. D. S. Wild, D. Sels, H. Pichler, C. Zanoci, and M. D. Lukin, Quantum sampling algorithms for near-term devices, Phys. Rev. Lett. 127, 100504 (2021).
  11. D. Layden, G. Mazzola, R. V. Mishmash, M. Motta, P. Wocjan, J.-S. Kim, and S. Sheldon, Quantum-enhanced Markov-chain Monte Carlo, Nature (London) 619, 282 (2023).
  12. J. Graham, K. Kanov, X. I. A. Yang, M. K. Lee, N. Malaya, C. C. Lalescu, R. Burns, G. Eyink, A. Szalay, R. D. Moser, and C. Meneveau, A web services accessible database of turbulent channel flow and its use for testing a new integral wall model for LES, J. Turbul. 17, 181 (2016).
  13. D. A. Levin and Y. Peres, Markov Chains and Mixing Times, 2nd ed. (American Mathematical Society, Providence, RI, 2017).
  14. Y. Nakano, H. Hakoshima, K. Mitarai, and K. Fujii, Markov-chain Monte Carlo method enhanced by a quantum alternating operator ansatz, Phys. Rev. Res. 6, 033105 (2024).
  15. A. Orfi and D. Sels, Bounding speedup of quantum-enhanced Markov chain Monte Carlo, Phys. Rev. A 110, 052414 (2024).
  16. Y. Nakano, K. N. Okada, and K. Fujii, Neural-network-assisted Monte Carlo sampling trained by quantum approximate optimization algorithm, PRX Quantum 7, 010338 (2026).
  17. S. Ferguson and P. Wallden, Quantum-enhanced Markov chain Monte Carlo for systems larger than a quantum computer, Phys. Rev. Res. 7, 013231 (2025).
  18. J. Ingelmann, F. Schindler, and J. Schumacher, Quantum-assisted Lagrangian tracer dispersion in turbulent shear flow, Phys. Fluids 37, 085251 (2025).
  19. J. Kim, P. Moin, and R. D. Moser, Turbulence statistics in fully developed channel flow at low Reynolds number, J. Fluid Mech. 177, 133 (1987).
  20. P. K. Yeung, Lagrangian investigations of turbulence, Annu. Rev. Fluid Mech. 34, 115 (2002).
  21. F. Toschi and E. Bodenschatz, Lagrangian properties of particles in turbulence, Annu. Rev. Fluid Mech. 41, 375 (2009).
  22. K. R. Sreenivasan and J. Schumacher, Lagrangian views on turbulent mixing of passive scalars, Philos. Trans. R. Soc. A 368, 1561 (2010).
  23. I. M. Mazzitelli, F. Fornarelli, A. S. Lanotte, and P. Oresta, Pair and multi-particle dispersion in numerical simulations of convective boundary layer turbulence, Phys. Fluids 26, 055110 (2014).
  24. I. M. Mazzitelli, F. Toschi, and A. S. Lanotte, An accurate and efficient Lagrangian sub-grid model, Phys. Fluids 26, 095101 (2014).
  25. B. L. Sawford, Reynolds number effects in Lagrangian stochastic models of turbulent dispersion, Phys. Fluids A 3, 1577 (1991).
  26. B. Viggiano, J. Friedrich, R. Volk, M. Bourgoin, R. Bayoàn Cal, and L. Chevillard, Modelling Lagrangian velocity and acceleration in turbulent flows as infinitely differentiable stochastic processes, J. Fluid Mech. 900, A27 (2020).
  27. F. T. M. Nieuwstadt, A large-eddy simulation of a line source in a convective atmospheric boundary layer—I. Dispersion characteristics, Atmos. Environ. Part A 26, 485 (1992).
  28. A. Dosio and J. Vilà-Guerau de Arrelano, Statistics of absolute and relative dispersion in the atmospheric convective boundary layer: A large-eddy simulation study, J. Atmos. Sci. 63, 1253 (2006).
  29. A. F. Stein, R. R. Draxler, G. D. Rolph, B. J. B. Stunder, M. D. Cohen, and F. Ngan, NOAA's HYSPLIT atmospheric transport and dispersion modeling system, Bull. Am. Meteor. Soc. 96, 2059 (2015).
  30. L. Hoffmann, G. Günther, D. Li, O. Stein, X. Wu, S. Griessbach, Y. Heng, P. Konopka, R. Müller, B. Vogel, and J. S. Wright, From ERA-Interim to ERA5: The considerable impact of ECMWF's next-generation reanalysis on Lagrangian transport simulations, Atmos. Chem. Phys. 19, 3097 (2019).
  31. L. Hoffmann, T. Rößler, S. Griessbach, Y. Heng, and O. Stein, Lagrangian transport simulations of volcanic sulfur dioxide emissions: Impact of meteorological data products, J. Geophys. Res. Atmos. 121, 4651 (2016).
  32. M. Lubasch, J. Joo, P. Moinier, M. Kiffner, and D. Jaksch, Variational quantum algorithms for nonlinear problems, Phys. Rev. A 101, 010301(R) (2020).
  33. B. N. Todorova and R. de Steijl, Quantum algorithm for the collisionless Boltzmann equation, J. Comput. Phys. 409, 109347 (2020).
  34. F. Gaitan, Finding flows in a Navier-Stokes fluid through quantum computing, npj Quantum Inf. 6, 61 (2021).
  35. L. Budinski, Quantum algorithm for the advection–diffusion equation simulated with the lattice Boltzmann method, Quantum Inf. Process. 20, 57 (2021).
  36. P. Pfeffer, F. Heyder, and J. Schumacher, Hybrid quantum-classical reservoir computing of thermal convection flow, Phys. Rev. Res. 4, 033176 (2022).
  37. P. Pfeffer, F. Heyder, and J. Schumacher, Reduced-order modeling of two-dimensional turbulent Rayleigh-Bénard flow by hybrid quantum-classical reservoir computing, Phys. Rev. Res. 5, 043242 (2023).
  38. S. S. Bharadwaj and K. R. Sreenivasan, Hybrid quantum algorithms for flow problems, Proc. Natl. Acad. Sci. USA 120, e2311014120 (2023).
  39. D. Jaksch, P. Givi, A. J. Daley, and T. Rung, Variational quantum algorithms for computational fluid dynamics, AIAA J. 61, 1885 (2023).
  40. Z. Meng and Y. Yang, Quantum computing of fluid dynamics using the hydrodynamic Schrödinger equation, Phys. Rev. Res. 5, 033182 (2023).
  41. J. Ingelmann, S. S. Bharadwaj, P. Pfeffer, K. R. Sreenivasan, and J. Schumacher, Two quantum algorithms for solving the one-dimensional advection–diffusion equation, Comput. Fluids 281, 106369 (2024).
  42. O. Ahmed, F. Tennie, and L. Magri, Prediction of chaotic dynamics and extreme events: A recurrence-free quantum reservoir computing approach, Phys. Rev. Res. 6, 043082 (2024).
  43. A. J. Pool, A. D. Somoza, C. Mc Keever, M. Lubasch, and B. Horstmann, Nonlinear dynamics as a ground-state solution on quantum computers, Phys. Rev. Res. 6, 033257 (2024).
  44. P. Brearley and S. Laizet, Quantum algorithm for solving the advection equation using Hamiltonian simulation, Phys. Rev. A 110, 012430 (2024).
  45. Z. Meng, J. Zhong, S. Xu, et al., Simulating unsteady flows on a superconducting quantum processor, Commun. Phys. 7, 349 (2024).
  46. N. Köcher, H. Rose, S. S. Bharadwaj, J. Schumacher, and S. Schumacher, Numerical solution of nonlinear Schrödinger equation by a hybrid pseudospectral-variational quantum algorithm, Sci. Rep. 15, 23478 (2025).
  47. S. S. Bharadwaj and K. R. Sreenivasan, Compact quantum algorithms for time-dependent differential equations, Phys. Rev. Res. 7, 023262 (2025).
  48. P. Over, S. Bengoechea, T. Rung, F. C. L. Scandurra, E. de Villiers, and D. Jaksch, Boundary treatment for variational quantum simulations of partial differential equations on quantum computers, Comput. Fluids 288, 106508 (2025).
  49. P. Pfeffer, P. Brearley, S. Laizet, and J. Schumacher, Spectral quantum algorithm for passive scalar transport in shear flows, Sci. Rep. 15, 41172 (2025).
  50. S. S. Bharadwaj and K. R. Sreenivasan, Quantum computation of fluid dynamics, Indian Acad. Sci. Conf. Ser. 3, 77 (2020).
  51. S. Succi, W. Itani, K. R. Sreenivasan, and R. Steijl, Quantum computing for fluids: Where do we stand? Europhys. Lett. 144, 10001 (2023).
  52. F. Tennie, S. Laizet, S. Lloyd, and L. Magri, Quantum computing for nonlinear differential equations and turbulence, Nat. Rev. Phys. 7, 220 (2025).
  53. P. Shen and P. K. Yeung, Fluid particle dispersion in homogeneous turbulent shear flow, Phys. Fluids 9, 3472 (1997).
  54. S. B. Pope, Stochastic Lagrangian models of velocity in homogeneous turbulent shear flow, Phys. Fluids 14, 1696 (2002).
  55. J. Schumacher, Derivative moments in stationary homogeneous shear turbulence, J. Fluid Mech. 441, 109 (2001).
  56. J. Schumacher, Relation between shear parameter and Reynolds number in statistically stationary turbulent shear flows, Phys. Fluids 16, 3094 (2004).
  57. S. B. Pope, Turbulent Flows (Cambridge University Press, Cambridge, UK, 2000).
  58. J. I. Polanco, I. Vinkovic, N. Stelzenmuller, N. Mordant, and M. Bourgoin, Relative dispersion of particle pairs in turbulent channel flow, Int. J. Heat Fluid Flow 71, 231 (2018).
  59. Y. Li, E. Perlman, M. Wan, Y. Yang, C. Meneveau, R. Burns, S. Chen, A. Szalay, and G. Eyink, A public turbulence database cluster and applications to study Lagrangian evolution of velocity increments in turbulence, J. Turbul. 9, N31 (2008).
  60. J. Schumacher and B. Eckhardt, On statistically stationary homogeneous shear turbulence, Europhys. Lett. 52, 627 (2000).
  61. P. Gualtieri, C. M. Casciola, R. Benzi, G. Amati, and R. Piva, Scaling laws and intermittency in homogeneous shear flow, Phys. Fluids 14, 583 (2002).
  62. S. C. C. Bailey and A. J. Smits, Experimental investigation of the structure of large- and very-large-scale motions in turbulent pipe flow, J. Fluid Mech. 651, 339 (2010).
  63. A. J. Smits and I. Marusic, Wall-bounded turbulence, Phys. Today 66(9), 25 (2013).
  64. D. J. Higham and P. E. Kloeden, An Introduction to the Numerical Simulation of Stochastic Differential Equations (SIAM, Philadelphia, 2021).
  65. A. M. Obukhov, Distribution of energy in the spectrum of turbulent flow, Izv. Akad. Nauk Ser. Geogr. Geophys. 4, 453 (1941).
  66. G. Boffetta and I. M. Sokolov, Statistics of two-particle dispersion in two-dimensional turbulence, Phys. Fluids 14, 3224 (2002).
  67. B. Birnir and L. Angheluta, Scaling of Lagrangian structure functions, Phys. Rev. Res. 7, 023225 (2025).
  68. J. Schumacher, Lagrangian dispersion and heat transport in convective turbulence, Phys. Rev. Lett. 100, 134502 (2008).
  69. M. Ettel, R. J. Samuel, M. Chertkov, and J. Schumacher, Lagrangian single-particle, multi-particle and topological analyses in turbulent Rayleigh-Bénard convection, arXiv:2605.22326.
  70. T. Gubaev, P. Pfeffer, C. Dreßler, and J. Schumacher, Variational quantum algorithm for anion exchange across electrolyzer membrane, Phys. Rev. Appl. (2026).
  71. S. K. Lam, A. Pitrou, and S. Seibert, Numba: A LLVM-based Python JIT compiler, in Proceedings of the Second Workshop on the LLVM Compiler Infrastructure in HPC (ACM, New York, NY, 2015), pp. 1–6.
  72. J. Ingelmann, Quantum algorithms for simple fluid flows (version v3), Zenodo, 2025, https://zenodo.org/records/16035235.
  73. F. Schindler and J. Schumacher, Quantum-enhanced Markov chain Monte Carlo sampling to model Lagrangian tracer dispersion in turbulent boundary layer, Zenodo, 2026, https://zenodo.org/records/20672071.
  74. S. Särkkä and A. Solin, Applied Stochastic Differential Equations (Cambridge University Press, Cambridge, UK, 2019).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation