Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Proceedings of SPIE, the International Society for Optical Engineering/Proceedings of SPIE

Quantum lattice gas algorithm for quantum turbulence and vortex reconnection in the Gross-Pitaevskii equation

George Vahala, Jeffrey Yepez, Linda L. Vahala

Published Mar 5, 2008• 3 citations
Publisher

Abstract

The ground state wave function for a Bose Einstein condensate is well described by the Gross-Pitaevskii equation. A Type-II quantum algorithm is devised that is ideally parallelized even on a classical computer. Only 2 qubits are required per spatial node. With unitary local collisions, streaming of entangled states and a spatially inhomogeneous unitary gauge rotation one recovers the Gross-Pitaevskii equation. Quantum vortex reconnection is simulated - even without any viscosity or resistivity (which are needed in classical vortex reconnection).

DOI
10.1117/12.777722

Related papers

Unitary Quantum Lattice Gas Algorithms for Quantum to Classical Turbulence

George Vahala, Jeffrey Yepez, Min Soe +2 more

Using a set of interleaved unitary collision-stream operators, a three-dimensional (3D) quantum lattice gas algorithm is devised which, on taking moments, recovers the Gross-Pitaevskii (GP) equation. If a zero-temperature Bose-Einstein condensate (BEC) is trapped in an a magnetic well, the evolution of the ground-state wave function satisfies the scalar GP equation, while if the BEC is trapped in an optical trap the ground-state wave function satisfies spin or GP equations. Quantum turbulence is studied in a scalar GP system on 5,7603 grid yielding not only the classical Kolmogorov k-5/3 cascade but also the quantum vortex k-3 spectrum. For a certain class of initial conditions, one finds an intermittent loss of tangled quantum vortices as the vortex cores attain minimal size, and thus prevent the Kelvin wave cascade (due to helical wave-wave coupling on the vortex). A coupled set of GP equations are solved for spin or BEC. Skrymions, which describe topologically-linked quantum vortices, are examined. One finds, for certain initial conditions that the incompressible kinetic energy spectrum for the condensate component of a vortex ring core rapidly departs from the k-3 linear quantum vortex spectrum.

DoD High Performance Computing Modernization Program Users Group Conference (HPCMP-UGC)•Jun 1, 2010•0 citations

Qubit unitary lattice algorithm for spin-2 Bose-Einstein condensates: II - vortex reconnection simulations and non-Abelain vortices

George Vahala, Min Soe, Linda L. Vahala

Using the qubit unitary algorithm developed in the previous paper, quantum reconnection of non-Abelian quantum vortices is simulated and rung-like vortices are seen during the reconnection. Vortices in two of the four tetrahedral symmetry group are examined. Moreover, the algorithm performs exceptionally on classical supercomputers, as shown on timings on the IBM-Q supercomputer MIRA. It is expected that one should see the exponential speedup on a quantum computer since quantum entanglement is crucial for the qubit algorithm to recover the Gross-Pitaevskii equations for spin-2 systems.Keywords : spin-2 BECs, non-Abeliam vortices

Radiation effects and defects in solids•Jan 2, 2020•8 citations

Unitary Quantum Lattice Algorithms for Turbulence

George Vahala, Linda L. Vahala

Abstract : A qubit unitary lattice algorithm is developed for the mean field evolution of the ground state of a Bose-Einstein condensate whether in a magnetic or optical trap. The algorithm consists of interleaved collide-stream operators which in diffusionordering recovers the spinor Gross Pitaevskii equations and is ideally parallelized. The algorithm is benchmarked againstexact one dimensional vector inelastic soliton collision solutions. Three dimensional quantum turbulence is examined forscalar Bose Einstein condensates and the compressible kinetic energy spectrum exhibits three different cascades

•May 23, 2016•0 citations

Modeling The Gross-Pitaevskii Equation Using The Quantum Lattice Gas Method

Armen Oganesov

We present an improved Quantum Lattice Gas (QLG) algorithm as a mesoscopic unitary perturbative representation of the mean field Gross Pitaevskii (GP) equation for Bose–Einstein Condensates (BECs). The method employs an interleaved sequence of unitary collide and stream operators. QLG is applicable to many different scalar potentials in the weak interaction regime and has been used to model the Korteweg–de Vries (KdV), Burgers and GP equations. It can be implemented on both quantum and classical computers and is extremely scalable. We present results for 1D soliton solutions with positive and negative internal interactions, as well as vector solitons with inelastic scattering. In higher dimensions we look at the behavior of vortex ring reconnection. A further improvement is considered with a proper operator splitting technique via a Fourier transformation. This is great for quantum computers since the quantum FFT is exponentially faster than its classical counterpart which involves non-local data on the entire lattice (Quantum FFT is the backbone of the Shor algorithm for quantum factorization). We also present an imaginary time method in which we transform the Schrödinger equation into a diffusion equation for recovering ground state initial conditions of a quantum system suitable for the QLG algorithm.

•Jan 1, 2017•1 citation

Quantum Lattice-Gas Algorithm for Quantum Turbulence - CAP Simulations on 12,288 Cores of Cray XT-5 Einstein at NAVO

George Vahala, Jeffrey Yepez, Min Soe +2 more

A novel unitary quantum lattice algorithm is developed to explore quantum turbulence. Because of its low memory requirements and its near perfect parallelization to the full 12,288 cores on the Cray XT5, simulations were run up to spatial grids of 5,7603. The Gross-Pitaevskii equation, which describes the ground state of a Bose Einstein condensate (BEC), is solved and it is found that the incompressible kinetic energy spectrum exhibits 3 distinct power laws: classical Kolmogorov k-5/3spectrum at scales much larger than the individual quantum vortex cores, and a quantum Kelvin wave cascade spectrum of k-3at scales of the order of the quantum cores. In the adjoining semiclassical regime, there is a steeper spectral decay transitioning between the classical and quantum regimes. However, its spectral exponent does not seem to be universal. This is the first, first-principle simulation yielding the universal quantum Kelvin cascade exponent.

DoD High Performance Computing Modernization Program Users Group Conference (HPCMP-UGC)•Jun 1, 2009•1 citation

Quantum Algorithm for Bose-Einstein Condensate Quantum Fluid Dynamics: Twisting of Filamentary Vortex Solitons Demarcated by Fast Poincare Recursion

Jeffrey Yepez, George Vahala, Linda L. Vahala

The dynamics of vortex solitons is studied in a BEC superfluid. A quantum lattice-gas algorithm (measurementbased quantum computation) is employed to examine the dynamical behavior vortex soliton solutions of the Gross-Pitaevskii equation (ø4 interaction nonlinear Schroedinger equation). Quantum turbulence is studied in large grid numerical simulations: Kolmogorov spectrum associated with a Richardson energy cascade occurs on large flow scales. At intermediate scales, a new k-6 power law emerges, due to vortex filamentary reconnections associated with Kelvin wave instabilities (vortex twisting) coupling to sound modes and the exchange of intermediate vortex rings. Finally, at very small spatial scales a k-3power law emerges, characterizing fluid dynamics occurring within the scale size of the vortex cores themselves. Poincaré recurrence is studied: in the free non-interacting system, a fast Poincaré recurrence occurs for regular arrays of line vortices. The recurrence period is used to demarcate dynamics driving a nonlinear quantum fluid towards turbulence, since fast recurrence is an approximate symmetry of the nonlinear quantum fluid at early times. This class of quantum algorithms is useful for studying BEC superfluid dynamics and, without modification, should allow for higher resolution simulations (with many components) on future quantum computers.

•Jan 1, 2009•0 citations
The ground state wave function for a Bose Einstein condensate is well described by the Gross-Pitaevskii equation. A Type-II quantum algorithm is devised that is ideally parallelized even on a classical computer. Only 2 qubits are required per spatial node. With unitary local collisions, streaming of entangled states and a spatially inhomogeneous unitary gauge rotation one recovers the Gross-Pitaevskii equation. Quantum vortex reconnection is simulated - even without any viscosity or resistivity (which are needed in classical vortex reconnection).
Using a set of interleaved unitary collision-stream operators, a three-dimensional (3D) quantum lattice gas algorithm is devised which, on taking moments, recovers the Gross-Pitaevskii (GP) equation. If a zero-temperature Bose-Einstein condensate (BEC) is trapped in an a magnetic well, the evolution of the ground-state wave function satisfies the scalar GP equation, while if the BEC is trapped in an optical trap the ground-state wave function satisfies spin or GP equations. Quantum turbulence is studied in a scalar GP system on 5,7603 grid yielding not only the classical Kolmogorov k-5/3 cascade but also the quantum vortex k-3 spectrum. For a certain class of initial conditions, one finds an intermittent loss of tangled quantum vortices as the vortex cores attain minimal size, and thus prevent the Kelvin wave cascade (due to helical wave-wave coupling on the vortex). A coupled set of GP equations are solved for spin or BEC. Skrymions, which describe topologically-linked quantum vortices, are examined. One finds, for certain initial conditions that the incompressible kinetic energy spectrum for the condensate component of a vortex ring core rapidly departs from the k-3 linear quantum vortex spectrum.
Using the qubit unitary algorithm developed in the previous paper, quantum reconnection of non-Abelian quantum vortices is simulated and rung-like vortices are seen during the reconnection. Vortices in two of the four tetrahedral symmetry group are examined. Moreover, the algorithm performs exceptionally on classical supercomputers, as shown on timings on the IBM-Q supercomputer MIRA. It is expected that one should see the exponential speedup on a quantum computer since quantum entanglement is crucial for the qubit algorithm to recover the Gross-Pitaevskii equations for spin-2 systems.Keywords : spin-2 BECs, non-Abeliam vortices
Abstract : A qubit unitary lattice algorithm is developed for the mean field evolution of the ground state of a Bose-Einstein condensate whether in a magnetic or optical trap. The algorithm consists of interleaved collide-stream operators which in diffusionordering recovers the spinor Gross Pitaevskii equations and is ideally parallelized. The algorithm is benchmarked againstexact one dimensional vector inelastic soliton collision solutions. Three dimensional quantum turbulence is examined forscalar Bose Einstein condensates and the compressible kinetic energy spectrum exhibits three different cascades
We present an improved Quantum Lattice Gas (QLG) algorithm as a mesoscopic unitary perturbative representation of the mean field Gross Pitaevskii (GP) equation for Bose–Einstein Condensates (BECs). The method employs an interleaved sequence of unitary collide and stream operators. QLG is applicable to many different scalar potentials in the weak interaction regime and has been used to model the Korteweg–de Vries (KdV), Burgers and GP equations. It can be implemented on both quantum and classical computers and is extremely scalable. We present results for 1D soliton solutions with positive and negative internal interactions, as well as vector solitons with inelastic scattering. In higher dimensions we look at the behavior of vortex ring reconnection. A further improvement is considered with a proper operator splitting technique via a Fourier transformation. This is great for quantum computers since the quantum FFT is exponentially faster than its classical counterpart which involves non-local data on the entire lattice (Quantum FFT is the backbone of the Shor algorithm for quantum factorization). We also present an imaginary time method in which we transform the Schrödinger equation into a diffusion equation for recovering ground state initial conditions of a quantum system suitable for the QLG algorithm.
A novel unitary quantum lattice algorithm is developed to explore quantum turbulence. Because of its low memory requirements and its near perfect parallelization to the full 12,288 cores on the Cray XT5, simulations were run up to spatial grids of 5,7603. The Gross-Pitaevskii equation, which describes the ground state of a Bose Einstein condensate (BEC), is solved and it is found that the incompressible kinetic energy spectrum exhibits 3 distinct power laws: classical Kolmogorov k-5/3spectrum at scales much larger than the individual quantum vortex cores, and a quantum Kelvin wave cascade spectrum of k-3at scales of the order of the quantum cores. In the adjoining semiclassical regime, there is a steeper spectral decay transitioning between the classical and quantum regimes. However, its spectral exponent does not seem to be universal. This is the first, first-principle simulation yielding the universal quantum Kelvin cascade exponent.
The dynamics of vortex solitons is studied in a BEC superfluid. A quantum lattice-gas algorithm (measurementbased quantum computation) is employed to examine the dynamical behavior vortex soliton solutions of the Gross-Pitaevskii equation (ø4 interaction nonlinear Schroedinger equation). Quantum turbulence is studied in large grid numerical simulations: Kolmogorov spectrum associated with a Richardson energy cascade occurs on large flow scales. At intermediate scales, a new k-6 power law emerges, due to vortex filamentary reconnections associated with Kelvin wave instabilities (vortex twisting) coupling to sound modes and the exchange of intermediate vortex rings. Finally, at very small spatial scales a k-3power law emerges, characterizing fluid dynamics occurring within the scale size of the vortex cores themselves. Poincaré recurrence is studied: in the free non-interacting system, a fast Poincaré recurrence occurs for regular arrays of line vortices. The recurrence period is used to demarcate dynamics driving a nonlinear quantum fluid towards turbulence, since fast recurrence is an approximate symmetry of the nonlinear quantum fluid at early times. This class of quantum algorithms is useful for studying BEC superfluid dynamics and, without modification, should allow for higher resolution simulations (with many components) on future quantum computers.