Skip to main content An official website of the United States government Here's how you know Here's how you know Official websites use .gov A .gov website belongs to an official government organization in the United States. Secure .gov websites use HTTPS A lock ( Lock Locked padlock icon ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites. Search Log in Dashboard Publications Account settings Log out Search… Search NCBI Primary site navigation Search Logged in as: Dashboard Publications Account settings Log in Search PMC Full-Text Archive Search in PMC Journal List User Guide PERMALINK Copy As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsement of, or agreement with, the contents by NLM or the National Institutes of Health. Learn more: PMC Disclaimer | PMC Copyright Notice Nat Commun . 2026 Feb 10;17:2602. doi: 10.1038/s41467-026-69168-8 Search in PMC Search in PubMed View in NLM Catalog Add to search Simulating fluid vortex interactions on a superconducting quantum processor Ziteng Wang Ziteng Wang 1 School of Aeronautics and Astronautics, Zhejiang University, Hangzhou, 310027 China Find articles by Ziteng Wang 1, # , Jiarun Zhong Jiarun Zhong 2 School of Physics, ZJU-Hangzhou Global Scientific and Technological Innovation Center, and Zhejiang Key Laboratory of Micro-nano Quantum Chips and Quantum Control, Zhejiang University, Hangzhou, 310027 China Find articles by Jiarun Zhong 2, # , Ke Wang Ke Wang 2 School of Physics, ZJU-Hangzhou Global Scientific and Technological Innovation Center, and Zhejiang Key Laboratory of Micro-nano Quantum Chips and Quantum Control, Zhejiang University, Hangzhou, 310027 China Find articles by Ke Wang 2 , Zitian Zhu Zitian Zhu 2 School of Physics, ZJU-Hangzhou Global Scientific and Technological Innovation Center, and Zhejiang Key Laboratory of Micro-nano Quantum Chips and Quantum Control, Zhejiang University, Hangzhou, 310027 China Find articles by Zitian Zhu 2 , Zehang Bao Zehang Bao 2 School of Physics, ZJU-Hangzhou Global Scientific and Technological Innovation Center, and Zhejiang Key Laboratory of Micro-nano Quantum Chips and Quantum Control, Zhejiang University, Hangzhou, 310027 China Find articles by Zehang Bao 2 , Chenjia Zhu Chenjia Zhu 1 School of Aeronautics and Astronautics, Zhejiang University, Hangzhou, 310027 China Find articles by Chenjia Zhu 1 , Wenwen Zhao Wenwen Zhao 1 School of Aeronautics and Astronautics, Zhejiang University, Hangzhou, 310027 China Find articles by Wenwen Zhao 1 , Yaomin Zhao Yaomin Zhao 3 State Key Laboratory for Turbulence and Complex Systems, School of Mechanics and Engineering Science, Peking University, and HEDPS-CAPT, Peking University, Beijing, 100871 China Find articles by Yaomin Zhao 3 , Yue Yang Yue Yang 3 State Key Laboratory for Turbulence and Complex Systems, School of Mechanics and Engineering Science, Peking University, and HEDPS-CAPT, Peking University, Beijing, 100871 China Find articles by Yue Yang 3 , Chao Song Chao Song 2 School of Physics, ZJU-Hangzhou Global Scientific and Technological Innovation Center, and Zhejiang Key Laboratory of Micro-nano Quantum Chips and Quantum Control, Zhejiang University, Hangzhou, 310027 China Find articles by Chao Song 2, ✉ , Shiying Xiong Shiying Xiong 1 School of Aeronautics and Astronautics, Zhejiang University, Hangzhou, 310027 China Find articles by Shiying Xiong 1, ✉ Author information Article notes Copyright and License information 1 School of Aeronautics and Astronautics, Zhejiang University, Hangzhou, 310027 China 2 School of Physics, ZJU-Hangzhou Global Scientific and Technological Innovation Center, and Zhejiang Key Laboratory of Micro-nano Quantum Chips and Quantum Control, Zhejiang University, Hangzhou, 310027 China 3 State Key Laboratory for Turbulence and Complex Systems, School of Mechanics and Engineering Science, Peking University, and HEDPS-CAPT, Peking University, Beijing, 100871 China ✉ Corresponding author. # Contributed equally. Received 2025 May 27; Accepted 2026 Jan 20; Collection date 2026. © The Author(s) 2026 Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/ . PMC Copyright notice PMCID: PMC13003102 PMID: 41667481 Abstract Vortex interactions are commonly observed in atmospheric turbulence, plasma dynamics, and collective behaviors in biological systems. However, accurately simulating these complex interactions is highly challenging due to the need to capture fine-scale details over extended timescales, which places computational burdens on traditional methods. In this study, we introduce a quantum vortex method, reformulating the Navier–Stokes (NS) equations within a quantum mechanical framework to enable the simulation of multi-vortex interactions on a quantum computer. We construct the effective Hamiltonian for the vortex system and implement a spatiotemporal evolution circuit to simulate its dynamics over prolonged periods. By leveraging eight qubits on a superconducting quantum processor with gate fidelities of 99.97% for single-qubit gates and 99.76% for two-qubit gates, we successfully reproduce natural vortex interactions. Overall, we establish a framework that reformulates vortex dynamics into a normalized wavefunction representation compatible with quantum system unitary evolution, combined with the designed spatiotemporal encoding scheme, providing a concrete pathway toward leveraging quantum resources in fluid systems. Subject terms: Fluid dynamics, Quantum simulation Despite recent progress, simulations of complex fluid dynamics phenomena on quantum computers remain challenging. Here, the authors present a quantum friendly reformulation of vortex dynamics with a spatiotemporal encoding scheme and apply it to simulations on a superconducting quantum processor. Introduction Vortices in fluids constitute a core component of complex flow behaviors, encompassing phenomena such as tropical cyclones 1 – 3 , ocean currents 4 – 6 , microfluidics 7 , 8 , as well as plasmas and magnetofluids 9 – 13 . Vortex interactions, which involve complex behaviors like vortex pairing and the leapfrogging effect (as shown in Fig. 1 a, b), affect energy transport, momentum exchange, and the scale cascade process in fluids, ultimately determining turbulence characteristics and their evolution 14 – 17 . However, simulating these critical and intricate structures using classical computation is highly challenging, as achieving the necessary spatial and temporal resolution to capture fine-scale details over extended timescales demands massive computational resources 18 – 21 , often exceeding practical limits. This complexity has spurred the development of advanced methods to address the computational bottlenecks while maintaining accuracy. Fig. 1. Overview for implementing vortex interactions using a superconducting quantum chip. Open in a new tab a Vortex pairs generated by paddling in natural fluid systems. b Laboratory-induced vortex interactions, leading to a leapfrogging configuration (Reprinted from Lim 69 , with the permission of AIP Publishing.). c Schematic of our superconducting quantum chip, where all qubits are arranged in a square lattice with nearest-neighbor couplings. Blue and green circles represent the spatial and temporal qubits used in our experiment. d The initial state of a four-vortex-particle system, which can be expressed as the tensor product of ∣ ψ 0 ⟩ that encodes the initial spatial information of the vortex particles and a uniform superposition state that encodes the temporal information. The quantum state is prepared with n p spatial qubits and n t temporal qubits. Here, we take n p = 2 and n t = 3 for illustration. e System’s final state, where each basis state of the temporal qubits is entangled with the corresponding spatial state, ∣ ψ i ⟩ , with i indexing the time step. f The quantum parallel evolution enabled by our circuit. F k represents the unitary that evolves the system for 2 k −1 time steps, with k ranging from 1 to n t . Recent progress in quantum computing presents a promising avenue to address these challenges, as emerging research on universal quantum partial differential equation (PDE)/ordinary differential equation (ODE) solvers 22 – 29 demonstrates potential for application in computational fluid dynamics (CFD) by leveraging quantum algorithms to replace key components of traditional solvers based on the fluid governing equations 30 – 37 . Additionally, alternative fluid dynamics descriptions optimized for quantum computing have been proposed, including quantum algorithms inspired by the lattice Boltzmann method 38 – 42 , quantum simulations based on Schrödingerization 43 – 45 , and the hydrodynamic Schrödinger equation 46 – 48 , which is inherently more suitable for quantum computing than the conventional Navier–Stokes (NS) equations 49 . Although quantum computing has demonstrated its potential in fluid mechanics, simulating fluid motion on actual quantum devices based on existing algorithms remains challenging. Current research into complex phenomena like vortex interactions remains partially constrained by the reliance on Eulerian methods, whose high spatial resolution requirements for accurate fluid behavior capture consequently increase quantum resource demands as qubit needs scale with grid resolution, making such implementations challenging on present noisy intermediate-scale quantum (NISQ) hardware 50 . We adopt Lagrangian vortex methods (LVM) that circumvent the limitations of Eulerian formulations, which are commonly used in the CFD community 51 . Furthermore, their intrinsic conservation laws governing vorticity evolution in high-Reynolds-number flows exhibit inherent compatibility with unitary quantum evolution 52 , 53 , thereby establishing a promising pathway for developing quantum simulations that preserve the conservation laws of fluid systems. Moreover, many quantum algorithms for simulating the time evolution of a system typically require a measurement at every time step to extract information necessary for studying dynamical behavior, computing physical quantities, and optimizing algorithms. However, since measurements collapse the quantum state, the quantum state at intermediate steps must be re-prepared to continue the computation. We design a novel spatiotemporal encoding scheme that embeds both spatial and temporal information directly into the quantum state. This approach prepares, in a single quantum-circuit structure, a quantum state that encodes information at multiple time steps in superposition, thereby eliminating the need for stepwise state preparation. We note that the measurement of a state for a specific time step requires repeated executions of the quantum circuit to collect samples, with the number of realizations for a given accuracy increasing with temporal qubits. In this work, we propose the quantum vortex method (QVM), which directly focuses on vortices themselves instead of relying on spatial discretization grids as in the Eulerian methods, thereby enabling the reformulation of complex vortex interactions in fluids within the framework of quantum computing. The QVM transforms the evolution of the vortex particle system into the evolution of a wavefunction. We adopt a data-driven strategy to train evolution modules that capture the dynamics of the wavefunction. Leveraging the trained modules, we then design an efficient spatiotemporal evolution circuit to implement the wavefunction propagation, where spatial qubits encode the spatial information of the vortex particle system, while auxiliary temporal qubits, initialized into a superposition state via Hadamard gates, act as placeholders for all time steps and later serve as control qubits to guide the evolution module in the spatial circuit. Building upon these theoretical developments, we implement the QVM on superconducting quantum processors to efficiently compute vortex interaction dynamics. This approach bridges classical fluid dynamics and quantum simulations, providing a new platform for exploring both quantum and classical vortex phenomena and offering a powerful tool for reinterpreting classical vortex dynamics from a quantum perspective. Results Quantum vortex method The fluid dynamics are governed by the NS equations for the velocity field u ( x , t ), which describe the evolution of the flow under the influence of pressure p , viscosity ν , and external forces f : D u D t = − 1 ρ ∇ p + ν ∇ 2 u + f , ∇ ⋅ u = 0 , 1 where D/D t = ∂/∂ t + u ⋅ ∇ is the material derivative, and ρ is the constant fluid density. To adapt the NS equations for quantum computing, we utilize the relationship between the vorticity field ω and the velocity field u ( ω = ∇ × u ), discretize the vorticity field into N p point vortices, and map their coordinates to complex variables, leading to the generalized Schrödinger equation: d ψ j ( t ) d t = i H j ( ψ 1 ( t ) , … , ψ N p ( t ) ) , j = 1 , … , N p . 2 Here H j denotes the j th component of a vortex-interaction-dependent Hamiltonian. The wave function ψ j is transformed from the j th vortex particle position ϕ j with ψ j = λ ϕ j − ∫ 0 t c ( τ ) d τ + C 0 , 3 where C 0 is an arbitrary constant, j indexes the vortex particles, and λ is a scaling factor that ensures the normalization condition: ∑ j = 1 N p ∣ ψ j ∣ 2 = 1 at t = 0 . The time-dependent function c ( t ) is defined as: c ( t ) = i λ 1 4 π ∑ j ≠ k , j = 1 , k = 1 N p Γ k ψ j ψ k * − ψ k ψ j * ∣ ψ j − ψ k ∣ 2 ∑ i N p ψ i * , 4 where (⋅) * denotes complex conjugation, and Γ k denotes the vortex strength of the k -th vortex. The evolution of the quantum vortex system is governed by equations that ensure the normalization of the quantum state, facilitating accurate simulations of fluid flows with reduced computational costs. Furthermore, we observe that when the vortex particle system exhibits collective motion in a certain direction, c ( t ) tends to remain relatively stable, exhibiting only minor fluctuations around a constant value. Therefore, we also provide a random sampling approximation method, in which we randomly select a subset of time instances of c ( t ) and average them to approximate their complete set. We remark that the normalization mapping the evolution of the vortex particle system to the evolution of a wavefunction is mathematically equivalent to the conventional vortex method. To enable experiments on present hardware, approximations and linearization are introduced to preserve leading-order vortex interactions in a linear representation. A detailed description of the QVM is in Supplementary Note 1 . Quantum encoding and evolution The fluid dynamics are governed by Eq. ( 1 ), with vortex particle positions represented by ϕ , which can be transformed into wave function ψ through Eq. ( 3 ), and ψ evolves according to Eq. ( 2 ). One may solve Eq. ( 1 ) on a grid, extract ϕ , and apply the transformation to obtain ψ , thus creating the data needed for training the nonlinear model described by Eq. ( 2 ). We investigate a system governed by Eq. ( 2 ) with N p vortices, discretizing time evolution into evenly spaced N t intervals. For the j th vortex particle, its position in configuration space is mapped to a complex variable ϕ j , with j ranging from 1 to N p . Subsequently, we introduce a transformation that shifts and scales the complex coordinates ϕ j to define new variables ψ j . This transformation ensures that each component ψ j of the wave function is properly normalized in terms of probability and remains conserved during the evolution governed by QVM. In the case of time discretization, the value of ψ j at the i -th time step is denoted by ψ j i . At each time step, the collection of ψ j i forms a state vector ∣ ψ i ⟩ = [ ψ 1 i , ψ 2 i , … , ψ N p i ] T . To achieve efficient evolution, the initial flow-field state ∣ ψ 0 ⟩ is first encoded into a larger quantum system. Specifically, the system’s initial state is prepared as a tensor product of the flow field’s initial state, which is encoded in n p = ⌈ log 2 N p ⌉ qubits, and a uniform superposition state encoded in n t = ⌈ log 2 N t ⌉ qubits. This procedure effectively generates multiple replicas of the flow field’s initial state, as shown in Fig. 1 d. These replicas explore different temporal evolutions simultaneously, and evolve in a branching manner as depicted in Fig. 1 f, eventually yielding a superposition of flow field states at 2 n t time steps ∣ ψ ⟩ = 1 N t [ ∣ ψ 0 ⟩ , ∣ ψ 1 ⟩ , … , ∣ ψ N t − 1 ⟩ ] T , as shown in Fig. 1 e. The evolution process is implemented with the quantum circuit illustrated in Fig. 2 . The core element of this circuit is the evolution unitary F k , with k ranging from 1 to n t , which evolves the state from i -th time step to ( i + 2 k −1 )-th time step as ∣ ψ i + 2 k − 1 ⟩ = F k ∣ ψ i ⟩ . At the beginning of the evolution, all qubits are initialized in the state 0 . The system then undergoes spatiotemporal evolution through a layered quantum circuit architecture, evolving the quantum state into ∣ ψ ⟩ . Specifically, the spatial qubits are initialized to the desired initial state through the “State Prep.” module, while the temporal qubits are prepared into a uniform superposition state via Hadamard gates. The temporal qubits then act as control qubits, which sequentially control the implementation of evolution unitaries F 1 , F 2 , …, F k on the spatial qubits. Fig. 2. The evolution circuit of the QVM method. Open in a new tab The circuit consists of n p spatial qubits encoding the spatial state of vortex particles and n t temporal qubits encoding the temporal information. The spatial qubits are initialized via the “State Prep." module, while the temporal qubits are prepared in a uniform superposition state using Hadamard gates. The quantum parallel evolution illustrated in Fig. 1 f is achieved through sequentially applying the controlled- F k operations to the system, with the temporal and spatial qubits being the control and target, respectively. With this controlled-unitary scheme, the quantum state undergoes a tree-like branching evolution from the initial state shown in Fig. 1 f, ultimately resulting in a superposition of all system states across 2 n t time steps. This design fully leverages the parallelism of quantum computing, significantly improving the efficiency of the simulation. See Supplementary Note 2 for more details. Experimental setup The algorithm is implemented with eight frequency-tunable transmon qubits on a flip-chip superconducting quantum processor, as shown in Fig. 1 c, where blue circles represent the spatial qubits and green circles represent the temporal qubits. Each qubit can be controlled and read out individually. The nearest-neighbouring qubits are connected with a tunable coupler for tuning on and off the effective coupling strength of the two qubits. The single-qubit gate with a length of 24 ns is realized by applying a Gaussian-shaped microwave pulse with DRAG correction 54 . The two-qubit CZ gate, with a duration of around 40 ns, is realized by tuning the frequencies and coupling strength of the two qubits to achieve a close-cycle diabatic transition, a global process involving both qubits as a combined system, between 11 (both qubits in their excited states) and 02 (the first qubit in the ground state and the second qubit in its second excited state), accumulating a π phase shift that transforms 11 into − 11 . The median parallel single-qubit gate and two-qubit gate fidelities are 99.97% and 99.76%, respectively. See Supplementary Note 6 for details. Nonlinear interactions in vortex systems The leapfrog vortex phenomenon 55 , described by the NS equation, refers to a mode of motion that occurs when two or more vortex rings interact, corresponding to four or more vortex particles in two dimensions. We consider the evolution of a four-vortex-particle system under Eq. ( 2 ), with the positions of the four particles initialized to be (0, 1), (0, 0.3), (0, − 1), and (0, − 0.3), respectively. We then apply the transformation defined in Eq. ( 3 ), with C 0 = −1.7903 and Γ = (1, 1, − 1, − 1) to obtain the corresponding quantum state. To learn the Hamiltonian of the vortex system, we select 100 vortex state pairs at time ( t i , t i + 1) to form the training set. Here, t i = 0.18( i − 1) with i = 1, …, 100 is equally sampled from a time range of [0, 18], which roughly corresponds to the period of a full leapfrogging cycle. In our experiment, we use two qubits to encode the spatial information of the four vortex particles. Additionally, six qubits are used to represent 64 time steps in the evolution. We then apply the QVM circuit based on the learned Hamiltonian to prepare the quantum state that encodes the spatiotemporal dynamics of the entire evolution process. To obtain the evolution trajectories of the vortex particles, we perform quantum state tomography (QST) on the two spatial qubits, while simultaneously conducting projective measurements on all temporal qubits. For each time step, we postselect the QST data for the respective temporal state to reconstruct the density matrix, from which the positions of the vortex particles can be extracted (see “Methods”). Note that additional global phases, which are experimentally unobservable, are numerically applied at each time step to preserve the symmetry of the system and constrain particle motion along the positive real axis. For comparison, we conduct ideal (noiseless) and noisy simulations using the same circuit as in the experiments. In the noisy simulation, we consider error models including depolarizing error of the single- and two-qubit gates, the qubit decoherence, and readout error, with the error rates obtained from experiments (see Supplementary Note 6 ). In Fig. 3 a, we plot the experimentally extracted trajectory of the four vortex particles for time steps outside of the training set, i.e., after t = 18. The results demonstrate a decent agreement between experimental data and noisy simulation. To characterize the experimental performance, we compare the reconstructed state of the spatial qubits and positions of the vortex particles with those obtained from noiseless numerical simulation (Fig. 3 b). The state fidelity values exceed 97% for all time steps (Fig. 3 b, upper panel). Besides, the position deviations of all four vortex particles from the noiseless simulated values remain below 0.2 throughout the evolution (Fig. 3 b, lower panel). Moreover, using the vortex particle position data from both experimental and noiseless simulation results, we reconstruct the velocity field for each time step based on the Biot-Savart formula and visualize it in Fig. 3 c, d, respectively. For all illustrated time steps, the velocity fields from experimental data are in close agreement with those of the noiseless simulation, demonstrating that two vortex rings, corresponding to four point vortices in 2D, alternately pass through and move forward in a real flow field while maintaining a degree of symmetry. An additional leapfrogging-type example in Supplementary Note 5 empirically probes how the learned evolution module F k responds to small perturbations of the initial condition, where we observe indications of generalization beyond the training trajectory. Fig. 3. Experimental results of nonlinear interactions in vortex systems. Open in a new tab a Trajectories of vortex particles obtained from ideal (noiseless) simulation, noisy simulation, and experimental data. b Fidelity and the position deviations as functions of time. The fidelity F at each time step t is defined as F = ∣ ⟨ ψ ideal t ∣ ψ exp t ⟩ ∣ 2 , where ∣ ψ exp t ⟩ and ∣ ψ ideal t ⟩ denote the state vector of the spatial qubits obtained through experiment and noiseless numerical simulation, respectively. The position deviation is defined as the Euclidean distance between the vortex particles in the experiment and that in noiseless simulation d = ∑ n = 1 N p ∣ r → exp n , t − r → ideal n , t ∣ , where r → exp n , t and r → ideal n , t denote the coordinates of the vortex particle n at time t obtained through experiment and noise-free simulation, respectively. c , d Velocity fields and streamlines induced by vortex particles at t = 24, 44, 62, and 81, obtained in the experiment ( c ) and noiseless simulation ( d ). Turbulent vortex particle system To demonstrate the robustness of our method, we further implement it to simulate the dynamics of an eight-vortex-particle system. The positions and vortex strengths of the eight vortex particles are initialized randomly, akin to a turbulent vortex particle system. We numerically perform the simulation using MindQuantum 56 , an open-source quantum computing framework for simulating and implementing quantum algorithms. We use three spatial qubits and nine temporal qubits to encode the positions and time steps, and then simulate the dynamics of the system under Eq. ( 2 ). Specifically, we select 64 equally spaced time steps, namely, N t ∈ {0, 4, 8, …, 252}, within the time range [0, 256] as training data to directly learn the implementation circuit using the variational quantum algorithm (VQA) described in the Methods. Applying the learned circuit to the first frame, we construct the wavefunctions for all time steps from 0 to 511. Figure 4 a–d visualizes the results at time steps 128 and 384, respectively. The progression from (a) to (c) illustrates how vortex dynamics evolve while maintaining coherent structures, whereas the corresponding velocity distributions from (b) to (d) quantitatively capture the spatial variations in flow velocity magnitude and direction. Fig. 4. The simulation results of the turbulent vortex particle system and the viscous vortex particle system. Open in a new tab a Flow field visualization rendered based on vortex particle position data from initial time to t = 128. b Velocity distribution of the flow field at t = 128. c Flow field visualization rendered based on vortex particle position data from the initial time to t = 384. d Velocity distribution of the flow field at t = 384. e , f Evolution of vortex particles in a viscous fluid computed using QVM ( e ) and LVM ( f ). Viscous vortex particle systems We now turn to a viscous system containing two vortex particles. For viscous vortex particle systems, our data-driven approach enables viscosity terms to be encoded within normalized quantum state vectors that preserve their physical properties during evolution. Vortex particle strengths are implicitly incorporated within the evolution module during training and rollout, and the velocity field can be recovered by applying empirical relations for the time evolution of strength. A detailed description can be found in Supplementary Note 1 . In this simple two-vortex system, both particle positions and viscous interactions can be represented through normalized wavefunctions, enabling our QVM to compute viscous evolution directly from the learned circuit using VQA on MindQuantum. In contrast, the classical Lagrangian vortex method faces limitations in directly incorporating viscosity terms into the vortex particle evolution. We employ high-precision grid-based Eulerian methods for solving the NS equation to compute the two-dimensional flow field and extract vortex particle positions as the dataset. The spatial information is encoded using n p = 1 qubits, while n t = 4 qubits are allocated for encoding the temporal steps, with the first four frames used to optimize the variational circuit parameters. A comparison of the computational results between the QVM and LVM in Fig. 4 e, f reveals that the former exhibits perfect agreement with ground truth data, whereas the latter demonstrates significant positional deviations indicative of strong viscous dissipation effects. Although 16 frames, corresponding to 4 qubits, are actually computed, only frames 0, 2, …, 14 are visualized due to space limitations. We remark that our method extends the idealized framework by learning viscous diffusion from data; the example of two co-rotating vortices in Fig. 4 shows vortex pre-merging dynamics. Although we implicitly handle vortex particle strengths, how to address topological changes induced by viscosity 57 , such as vortex splitting, merging, or reconnection, has not yet been investigated due to the fixed structure of the quantum representation. In conventional vortex methods, various particle reseeding and local re-meshing techniques, including vortex element splitting, merging, and deletion, have been developed to handle vortex reconnection 58 , 59 . From an implementation perspective, we could initialize partially empty vortex element data structures to accommodate variations in vortex element count, suggesting that incorporating auxiliary registers might provide a potential solution for handling vortex reconnection in future work. We provide a more detailed discussion in Supplementary Note 1 . Discussion This study introduces a quantum algorithmic framework designed to simulate intricate vortex interactions in fluid dynamics. By directly encoding vortex information into quantum states, the approach circumvents the inherent challenges associated with quantum encoding of fluid fields. The proposed algorithm is validated through numerical simulations of turbulent and viscous flows, as well as experimental simulations of the leapfrogging vortex phenomenon on a superconducting quantum processor. The present QVM is most effective when the flow admits a compact Lagrangian description with a small number of coherent vortices. In such cases, the qubit count scales with the number of vortices rather than a spatial grid; time-parallel encoding reduces repeated state preparations; and measured observables can be tailored to vortex-level quantities. By contrast, when smooth potential fields are of primary interest or full-field observables are essential, Schrödinger/Madelung-based formulations provide a natural Eulerian pathway and remain preferable. We therefore view the two families of approaches as complementary, and we outline extending our framework with data-driven reduced operators and improved measurement strategies as a path toward bridging more complex scenarios. We provide a detailed discussion in Supplementary Note 1 . The theoretical benefit in circuit evolution primarily depends on the complexity of the evolution module F k and its count O ( log N t ) , where N t denotes the number of discrete time steps. For different temporal qubit indices k , the complexity of F k shows no significant difference; as such, we characterize each module as O ( ḡ CZ ) . We conclude that the overall complexity of the quantum circuit execution is O ( ḡ CZ log N t ) , demonstrating a quantum benefit for predicting long-time flow evolution as the algorithmic complexity grows logarithmically with N t . We provide a detailed discussion in Supplementary Note 4 . Our spatiotemporal encoding approach leverages both spatial and temporal degrees of freedom to expand the Hilbert space available for information storage, enabling an exponential increase in capacity compared to classical systems of similar scale. Specifically, our method processes all time steps in a quantum superposition within a single workflow, requiring log N p + log N t qubits compared with the classical N p N t units of storage for the same functionality. For sequential time-stepping dynamics, classical methods require storage scaling linearly with system size N p , while our quantum approach exploits superposition and entanglement to encode spatial information efficiently. Despite logarithmic overhead from the time register, the quantum method still maintains a benefit as the problem size increases. This high-density encoding scheme is particularly well-suited for storing dynamic or high-dimensional data such as neural network parameters or physical trajectories. It also aligns naturally with quantum memory architectures, allowing efficient data retrieval via quantum algorithms, including Grover’s search and quantum random access memory 60 . Applications span a wide range of fields, including artificial intelligence 61 , where data and models can be encoded and processed in parallel, scientific simulations involving complex many-body or time-evolving systems, and quantum cryptography 62 , where secure and scalable storage of large keys or quantum states is essential. Phenomena governed by viscosity and small-scale physics, such as splitting, co-rotating merging, are not yet fully modeled. A complete quantum vortex particle solver for the incompressible NS equations remains a subject for further study. Future advances in measurement techniques and error mitigation strategies, such as quantum error correction 63 , noise filtering, and more efficient tomography 64 , along with the development of novel quantum algorithms that reduce the need for intensive measurements, could further alleviate the computational burden and enhance the efficiency of quantum simulations. Methods Implementing the evolution modules To implement the quantum circuit module, we leverage a data-driven approach. We approximate the N p -particle system as a linear system described by a parameterized effective Hamiltonian H eff ( θ ), expressed as an N p × N p complex Hermitian matrix. The training time range T train is uniformly divided into N train ( N train ≥ N p 2 ) segments, based on which we extract N train state pairs separated by a step size Δ t train as our training data. While Δ t train typically matches the evolution step size, it can be adjusted for specific purposes such as data interpolation. For the training process, we construct a temporary unitary evolution operator U ( θ ; Δ t train ) = e − i H eff ( θ ) Δ t train based on the chosen step size Δ t train . We optimize the parameterized unitary matrix to match the evolution for all state pairs in the training set, thereby obtaining a linear dynamics. The evolution matrices F k are then constructed similarly through e − i H eff ( θ ) 2 k − 1 Δ t predict , where the time step Δ t predict can be theoretically arbitrary. With the determined evolution matrices F k , the CPFlow method 65 is then employed to obtain the quantum circuit implementation with the optimal fidelity and number of CZ gates corresponding to F k . The complexity of the operator F k is heuristic and scales empirically with the system size, as in the quantum approximate optimization algorithm (QAOA) 66 , 67 . Further details on the practical implementation and complexity of F k are provided in Supplementary Note 3 . In addition to introducing a Hamiltonian as an intermediary to accommodate hardware limitations, the evolution module can alternatively be implemented by designing a parameterized quantum circuit ansatz and employing VQA to minimize the loss function, thereby learning vortex dynamics directly from data. We theoretically derive the complete quantum circuit for loss evaluation and gradient computation methods in Supplementary Note 3 . To validate the generalization and practical implementability of our QVM, we evaluate the transfer performance of its evolution module F k to nearby initial conditions against a configuration-space dynamic mode decomposition (DMD) surrogate trained on identical data. Results demonstrate that the solution error of our method remains stable across varying initial conditions, indicating that the trained operators can be effectively reused for systems with the same vortex count and similar configurations. Both methods show comparable performance within the training horizon under small deviations, while QVM exhibits slower error growth and tighter uncertainty bands under larger perturbations. We hypothesize that the normalization and unitary propagation act as an implicit regularizer under distribution shift, improving robustness at the cost of a small in-distribution bias. Future work could explore a hybrid approach combining DMD’s interpretability with QVM’s hardware-native generalization for enhanced performance and adaptability. A detailed configuration and representation are provided in Supplementary Note 5 . Extracting the spatial information In our experiment, we use QST to obtain the density matrix of the spatial qubits. We then extract the spatial information of the vortex particles from the eigenstate of the density matrix with the maximum eigenvalue. The method is valid for the depolarization error with a sufficiently small error rate p . Under the depolarizing error channel, the experimental density matrix can be modeled as ρ exp → (1 − p ) ρ + ( p /2 N ) I , where ρ = ∣ ψ ⟩ ⟨ ψ ∣ is the ideal density matrix, and N is the number of qubits. Thus, ∣ ψ ⟩ is also the eigenstate of the experimental density matrix with an eigenvalue of 1 − p + p /2 N , which remains the largest among all eigenvalues for a small p . However, in our experiment, there are also coherent error channels, which can change the eigenstates of ρ exp , introducing additional errors during the spatial information extraction procedure. To mitigate this error, we apply Pauli Twirling 68 , which can effectively transform all errors into depolarizing errors. Specifically, we average the QST data obtained from 50 equivalent circuits, which are generated by randomly replacing the CZ gates in the original circuit with 16 equivalent gates realized by applying additional Pauli gates before and after the CZ gate. Supplementary information Supplementary Information (829.6KB, pdf) Transparent Peer Review file (466.2KB, pdf) Acknowledgements We thank the Superconducting Quantum Computing Group at Zhejiang University for supporting the device and the experimental platform on which the experiment was carried out. The device was fabricated at the Micro-Nano Fabrication Center of Zhejiang University. We thank Prof. Bo Zhu for the photograph used in Fig. 1 a. The authors acknowledge the support of the National Key Research and Development Program of China (Grant No. 2023YFB4502600, received by C.S.) and the National Natural Science Foundation of China (Grants No. 12302294, received by S.X., Grants No. 12432010, received by Y.Y., and Grants No. 12525201, received by Y.Y.). Author contributions S.X. conceived the theoretical ideas. J.Z. and K.W. carried out the experiment under the supervision of C.S.; Z.W. and J.Z. designed the digital quantum circuits under the supervision of C.S.; Z.W., S.X., C.Z., W.Z., Y.Z. and Y.Y. conducted the theoretical analysis. J.Z., K.W., Z.Z., and Z.B. contributed to the experimental setup. All authors contributed to data analysis, discussion of the results, and writing of the paper. Peer review Peer review information Nature Communications thanks Ljubomir Budinski, and the other, anonymous, reviewer for their contribution to the peer review of this work. A peer review file is available. Data availability The data generated in this study have been deposited in the Figshare database (DOI: 10.6084/m9.figshare.30698342). Code availability The codes for numerical simulations are available in the Figshare database (DOI: 10.6084/m9.figshare.30698342). Competing interests The authors declare no competing interests. Footnotes Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. These authors contributed equally: Ziteng Wang, Jiarun Zhong. Contributor Information Chao Song, Email: [email protected]. Shiying Xiong, Email: [email protected]. Supplementary information The online version contains supplementary material available at 10.1038/s41467-026-69168-8. References 1. Fiorino, M. & Elsberry, R. L. Some aspects of vortex structure related to tropical cyclone motion. J. Atmos. Sci. 46 , 975 (1989). [ Google Scholar ] 2. Rios-Berrios, R. et al. A review of the interactions between tropical cyclones and environmental vertical wind shear. J. Atmos. Sci. 81 , 713 (2024). [ Google Scholar ] 3. Wang, Y. & Wu, C.-C. Current understanding of tropical cyclone structure and intensity changes—a review. Meteorol. Atmos. Phys. 87 , 257 (2004). [ Google Scholar ] 4. McWilliams, J. C. Submesoscale, coherent vortices in the ocean. Rev. Geophys. 23 , 165 (1985). [ Google Scholar ] 5. Zhang, Z., Wang, G., Wang, H. & Liu, H. Three-dimensional structure of oceanic mesoscale eddies. Ocean-Land-Atmos. Res. 3 , 0051 (2024). [ Google Scholar ] 6. Jebri, F. et al. Absence of the great whirl giant ocean vortex abates productivity in the Somali upwelling region. Commun. Earth Environ. 5 , 20 (2024). [ Google Scholar ] 7. Han, J. Y., La Fiandra, J. N. & DeVoe, D. L. Microfluidic vortex focusing for high throughput synthesis of size-tunable liposomes. Nat. Commun. 13 , 6997 (2022). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 8. Thurgood, P. et al. Dynamic vortex generation, pulsed injection, and rapid mixing of blood samples in microfluidics using the tube oscillation mechanism. Anal. Chem. 95 , 3089 (2023). [ DOI ] [ PubMed ] [ Google Scholar ] 9. Tanaka, M. Y. Vortex in plasma. Rev. Mod. Plasma Phys. 3 , 9 (2019). [ Google Scholar ] 10. Izhovkina, N. I., Artekha, S. N., Erokhin, N. S. & Mikhailovskaya, L. A. Interaction of atmospheric plasma vortices. Pure Appl. Geophys. 173 , 2945 (2016). [ Google Scholar ] 11. Xiong, S. & Yang, Y. Effects of twist on the evolution of knotted magnetic flux tubes. J. Fluid Mech. 895 , A28 (2020). [ Google Scholar ] 12. Aluie, H. Hydrodynamic and Magnetohydrodynamic Turbulence: Invariants, Cascades, and Locality , Ph.D. thesis (The Johns Hopkins University, 2009). 13. Belashov, V. Y. Interaction of n-vortex structures in a continuum, including atmosphere, hydrosphere and plasma. Adv. Space Res. 60 , 1878 (2017). [ Google Scholar ] 14. Shen, W., Yao, J. & Yang, Y. Designing turbulence with entangled vortices. Proc. Natl. Acad. Sci. USA 121 , e2405351121 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 15. Almoguera, A. P. & Kivotides, D. Vortex dynamics of turbulent energy cascades. Phys. Fluids 36 , 125156 (2024). [ Google Scholar ] 16. Majda, A. J. & Bertozzi, A. L. Vorticity and Incompressible Flow . (Cambridge University Press, 2002). 17. Morton, B. R. The generation and decay of vorticity. Geophys. Astrophys. Fluid Dyn. 28 , 277 (1984). [ Google Scholar ] 18. Ishihara, T., Gotoh, T. & Kaneda, Y. Study of high–reynolds number isotropic turbulence by direct numerical simulation. Annu. Rev. Fluid Mech. 41 , 165 (2009). [ Google Scholar ] 19. Ferziger, J. H., Perić, M. & Street, R. L. Computational Methods for Fluid Dynamics (Springer, 2019). 20. Kochkov, D. et al. Machine learning–accelerated computational fluid dynamics. Proc. Natl. Acad. Sci. 118 , e2101784118 (2021). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 21. Blocken, B. Computational Fluid Dynamics for urban physics: Importance, scales, possibilities, limitations and ten tips and tricks towards accurate and reliable simulations. Build. Environ. 91 , 219 (2015). [ Google Scholar ] 22. Liu, J.-P. et al. Efficient quantum algorithm for dissipative nonlinear differential equations. Proc. Natl. Acad. Sci. U.S.A. 118 , e2026805118 (2021). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 23. Tennie, F. & Magri, L. Integration of the Fokker--Planck equation on quantum computers: a new path to modelling nonlinear dynamics. Proc. R. Soc. A 481 , 2326 (2025). 24. Lloyd, S. et al. Quantum algorithm for nonlinear differential equations. Preprint at arXiv 10.48550/arXiv.2011.06571 (2020). 25. Pfeffer, P., Heyder, F. & Schumacher, J. Hybrid quantum-classical reservoir computing of thermal convection flow. Phys. Rev. Res. 4 , 033176 (2022). [ Google Scholar ] 26. Kyriienko, O., Paine, A. E. & Elfving, V. E. Solving nonlinear differential equations with differentiable quantum circuits. Phys. Rev. A 103 , 052416 (2021). [ Google Scholar ] 27. Lubasch, M., Joo, J., Moinier, P., Kiffner, M. & Jaksch, D. Variational quantum algorithms for nonlinear problems. Phys. Rev. A 101 , 010301 (2020). [ Google Scholar ] 28. Harrow, A. W., Hassidim, A. & Lloyd, S. Quantum algorithm for linear systems of equations. Phys. Rev. Lett. 103 , 150502 (2009). [ DOI ] [ PubMed ] [ Google Scholar ] 29. Childs, A. M., Kothari, R. & Somma, R. D. Quantum algorithm for systems of linear equations with exponentially improved dependence on precision. SIAM J. Comput. 46 , 1920 (2017). [ Google Scholar ] 30. Steijl, R. & Barakos, G. N. Parallel evaluation of quantum algorithms for computational fluid dynamics. Comput. Fluids 173 , 22 (2018). [ Google Scholar ] 31. Gaitan, F. Finding flows of a navier–stokes fluid through quantum computing. NPJ Quantum Inf. 6 , 61 (2020). [ Google Scholar ] 32. Chen, Z.-Y. et al. Quantum approach to accelerate finite volume method on steady computational fluid dynamics problems. Quantum Inf. Process. 21 , 137 (2022). [ Google Scholar ] 33. Lapworth, L. A hybrid quantum-classical CFD methodology with benchmark HHL solutions. Preprint at arXiv 10.48550/arXiv.2206.00419 (2022). 34. Demirdjian, R., Gunlycke, D., Reynolds, C. A., Doyle, J. D. & Tafur, S. Variational quantum solutions to the advection–diffusion equation for applications in fluid dynamics. Quantum Inf. Process. 21 , 322 (2022). [ Google Scholar ] 35. Succi, S., Itani, W., Sanavio, C., Sreenivasan, K. R. & Steijl, R. Ensemble fluid simulations on quantum computers. Comput. Fluids 270 , 106148 (2024). [ Google Scholar ] 36. Liu, B., Zhu, L., Yang, Z. & He, G. Quantum implementation of numerical methods for convection-diffusion equations: Toward computational fluid dynamics. Commun. Comput. Phys. 33 , 425 (2023). [ Google Scholar ] 37. Jaksch, D., Givi, P., Daley, A. J. & Rung, T. Variational quantum algorithms for computational fluid dynamics. AIAA J. 61 , 1885 (2023). [ Google Scholar ] 38. Budinski, L. Quantum algorithm for the Navier–Stokes equations by using the streamfunction-vorticity formulation and the lattice Boltzmann method. Int. J. Quantum Inf. 20 , 2150039 (2022). [ Google Scholar ] 39. Itani, W., Sreenivasan, K. R. & Succi, S. Quantum algorithm for lattice Boltzmann (QALB) simulation of incompressible fluids with a nonlinear collision term. Phys. Fluids 36 , 017112 (2024). [ Google Scholar ] 40. Sanavio, C. & Succi, S. Lattice Boltzmann–Carleman quantum algorithm and circuit for fluid flows at moderate Reynolds number. AVS Quantum Sci. 6 , 023802 (2024). [ Google Scholar ] 41. Budinski, L. Quantum algorithm for the advection–diffusion equation simulated with the lattice Boltzmann method. Quantum Inf. Process. 20 , 57 (2021). [ Google Scholar ] 42. Wang, B., Meng, Z., Zhao, Y. & Yang, Y. Quantum lattice Boltzmann method for simulating nonlinear fluid dynamics. NPJ Quantum Inf. 11 , 196 (2025). 43. Jin, S., Liu, N. & Yu, Y. Quantum simulation of partial differential equations via schrödingerization. Phys. Rev. Lett. 133 , 230602 (2024). [ DOI ] [ PubMed ] [ Google Scholar ] 44. Jin, S., Liu, N. & Yu, Y. Quantum simulation of partial differential equations: applications and detailed analysis. Phys. Rev. A 108 , 032603 (2023). [ Google Scholar ] 45. Jin, S. & Liu, N. Analog quantum simulation of partial differential equations. Quantum Sci. Technol. 9 , 035047 (2024). [ DOI ] [ PubMed ] [ Google Scholar ] 46. Meng, Z. & Yang, Y. Quantum computing of fluid dynamics using the hydrodynamic schrödinger equation. Phys. Rev. Res. 5 , 033182 (2023). [ Google Scholar ] 47. Meng, Z. et al. Simulating unsteady flows on a superconducting quantum processor. Commun. Phys. 7 , 349 (2024). [ Google Scholar ] 48. Salasnich, L., Succi, S. & Tiribocchi, A. Quantum wave representation of dissipative fluids. Int. J. Mod. Phys. C 35 , 2450100 (2024). [ Google Scholar ] 49. Succi, S., Itani, W., Sreenivasan, K. & Steijl, R. Quantum computing for fluids: Where do we stand?. Europhys. Lett. 144 , 10001 (2023). [ Google Scholar ] 50. Preskill, J. Quantum computing in the NISQ era and beyond. Quantum 2 , 79 (2018). [ Google Scholar ] 51. Mimeau, C. & Mortazavi, I. A review of vortex methods and their applications: from creation to recent advances. Fluids 6 , 68 (2021). [ Google Scholar ] 52. Helmholtz, H. Über Integrale der hydrodynamischen Gleichungen, welche den Wirbelbewegungen entsprechen. J. Reine Angew. Math. 55 , 25 (1858). [ Google Scholar ] 53. Ertel, H. Ein neuer hydrodynamischer Wirbelsatz. Meteorol. Z. 59 , 271 (1942). [ Google Scholar ] 54. Kelly, J. et al. Optimal quantum control using randomized benchmarking. Phys. Rev. Lett. 112 , 240504 (2014). [ DOI ] [ PubMed ] [ Google Scholar ] 55. Aref, H. Integrable, chaotic, and turbulent vortex motion in two-dimensional flows. Annu. Rev. Fluid Mech. 15 , 345 (1983). [ Google Scholar ] 56. Xu, X. et al. Mindspore quantum: a user-friendly, high-performance, and AI-compatible quantum computing framework. Preprint at arXiv 10.48550/arXiv.2406.17248 (2024). 57. Kida, S. & Takaoka, M. Vortex reconnection. Annu. Rev. Fluid Mech. 26 , 169 (1994). [ Google Scholar ] 58. Chorin, A. J. Hairpin removal in vortex interactions. J. Comput. Phys. 91 , 1 (1990). [ Google Scholar ] 59. Xiong, S., Tao, R., Zhang, Y., Feng, F. & Zhu, B. Incompressible flow simulation on vortex segment clouds. ACM Trans. Graph. 40 , 98 (2021). [ Google Scholar ] 60. Ladd, T. D. et al. Quantum computers. Nature 464 , 45 (2010). [ DOI ] [ PubMed ] [ Google Scholar ] 61. Dean, J. & Ghemawat, S. Mapreduce: Simplified data processing on large clusters. Commun. ACM 51 , 107 (2004). [ Google Scholar ] 62. Shor, P. W. Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer. SIAM J. Comput. 26 , 1484 (1997). [ Google Scholar ] 63. Google Quantum AI and Collaborators Quantum error correction below the surface code threshold. Nature 638 , 920 (2025). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 64. Aaronson, S. Shadow tomography of quantum states. SIAM J. Comput. 49 , 5 (2020). [ Google Scholar ] 65. Nemkov, N. A., Kiktenko, E. O., Luchnikov, I. A. & Fedorov, A. K. Efficient variational synthesis of quantum circuits with coherent multi-start optimization. Quantum 7 , 993 (2023). [ Google Scholar ] 66. Zhou, L., Wang, S.-T., Choi, S., Pichler, H. & Lukin, M. D. Quantum approximate optimization algorithm: Performance, mechanism, and implementation on near-term devices. Phys. Rev. X 10 , 021067 (2020). [ Google Scholar ] 67. He, Z. et al. Parameter setting heuristics make the quantum approximate optimization algorithm suitable for the early fault-tolerant era. in Proc. 43rd IEEE/ACM Int. Conf. Comput.-Aided Des. (ICCAD) pp 1–7 (2024). 68. Dür, W., Hein, M., Cirac, J. I. & Briegel, H.-J. Standard forms of noisy quantum operations via depolarization. Phys. Rev. A 72 , 052326 (2005). [ Google Scholar ] 69. Lim, T. T. A note on the leapfrogging between two coaxial vortex rings at low Reynolds numbers. Phys. Fluids 9 , 239 (1997). [ Google Scholar ] Associated Data This section collects any data citations, data availability statements, or supplementary materials included in this article. Supplementary Materials Supplementary Information (829.6KB, pdf) Transparent Peer Review file (466.2KB, pdf) Data Availability Statement The data generated in this study have been deposited in the Figshare database (DOI: 10.6084/m9.figshare.30698342). The codes for numerical simulations are available in the Figshare database (DOI: 10.6084/m9.figshare.30698342). Articles from Nature Communications are provided here courtesy of Nature Publishing Group ACTIONS View on publisher site PDF (8.2 MB) Cite Collections Permalink PERMALINK Copy RESOURCES Similar articles Cited by other articles Links to NCBI Databases Cite Copy Download .nbib .nbib Format: AMA APA MLA NLM Add to Collections Create a new collection Add to an existing collection Name your collection * Choose a collection Unable to load your collection due to an error Please try again Add Cancel Follow NCBI NCBI on X (formerly known as Twitter) NCBI on Facebook NCBI on LinkedIn NCBI on GitHub NCBI RSS feed Connect with NLM NLM on X (formerly known as Twitter) NLM on Facebook NLM on YouTube National Library of Medicine 8600 Rockville Pike Bethesda, MD 20894 Web Policies FOIA HHS Vulnerability Disclosure Help Accessibility Careers NLM NIH HHS USA.gov Back to Top