CyberSec Research

We present a MATLAB script which can use GPU acceleration to simulate a trapped ion interacting with a low-density cloud of atoms. This script, called atomiongpu.m, can massively parallelize MD simulations of trajectories of a trapped ion and an atom starting far away. The script uses ode45gpu, which is our optimized and specialized implementation of the Runge-Kutta algorithm used in MATLAB's ODE solver ode45. We first discuss the physical system and show how ode45gpu can solve it up to 22x faster than MATLAB's ode45. Then, we show how to easily modify the inputs to atomiongpu.m to account for different kinds of atoms, ions, atom-ion interactions, trap potentials, simulation parameters, initial conditions, and computational hardware, so that atomiongpu.m automatically finds the probability of complex formation, the distribution of observables such as the scattering angle and complex lifetime, and plots of specific trajectories.

Landau levels (LLs) are the massively-degenerate discrete energy spectrum of a charged particle in a transverse magnetic field and lie at the heart of many intriguing phenomena such as the integer and fractional quantum Hall effects as well as quantized vortices. In this Letter, we consider coupling of LLs of a transversely driven charge neutral particle in a synthetic gauge potential to a quantized field of an optical cavity -- a setting reminiscent of superradiant self-ordering setups in quantum gases. We uncover that this complex system can be surprisingly described in terms of two highly nonlinearly-coupled quantum harmonic oscillators, thus enabling a full quantum mechanical treatment. Light-matter coupling mixes the LLs and the superradiant photonic mode, leading to the formation of hybrid states referred to as ``Landau polaritons''. They inherit partially the degeneracy of the LLs and possess intriguing features such as non-zero light-matter entanglement and quadrature squeezing. Depending on the system parameters and the choice of initial state, the system exhibits diverse nonequilibrium quantum dynamics and multiple steady states, with distinct physical properties. This work lays the foundation for further investigating the novel, driven-dissipative Landau-polariton physics in quantum-gas--cavity-QED settings.

The Hopfield model describes a neural network that stores memories using all-to-all-coupled spins. Memory patterns are recalled under equilibrium dynamics. Storing too many patterns breaks the associative recall process because frustration causes an exponential number of spurious patterns to arise as the network becomes a spin glass. Despite this, memory recall in a spin glass can be restored, and even enhanced, under quantum-optical nonequilibrium dynamics because spurious patterns can now serve as reliable memories. We experimentally observe associative memory with high storage capacity in a driven-dissipative spin glass made of atoms and photons. The capacity surpasses the Hopfield limit by up to seven-fold in a sixteen-spin network. Atomic motion boosts capacity by dynamically modifying connectivity akin to short-term synaptic plasticity in neural networks, realizing a precursor to learning in a quantum-optical system.

The interplay between topology and strong interactions gives rise to a variety of exotic quantum phases, including fractional quantum Hall (FQH) states and their lattice analogs - fractional Chern insulators (FCIs). Such topologically ordered states host fractionalized excitations, which for spinful systems are often accompanied by ferromagnetism and skyrmions. Here, we study a Hofstadter-Hubbard model of spinful fermions on a square lattice, extended by nearest-neighbor interactions. Using large-scale density matrix renormalization group (DMRG) simulations, we demonstrate the emergence of a spin-polarized $\frac{1}{3}$-Laughlin-like FCI phase, characterized by a quantized many-body Chern number, a finite charge gap, and hidden off-diagonal long-range order. We further investigate the quantum Hall ferromagnet at $\nu=1$ and its skyrmionic excitations upon doping. In particular, we find that nearest-neighbor repulsion is sufficient to stabilize both particle- and hole-skyrmions in the ground state around $\nu=1$, whereas we do not find such textures around $\nu=\frac{1}{3}$. The diagnostic toolbox presented in this work, based on local densities, correlation functions, and spin-resolved observables, is directly applicable in quantum gas microscopy experiments. Our results open new pathways for experimental exploration of FCIs with spin textures in both ultracold atom and electronic systems.

We investigate the ground-state phases and spin-scissors dynamics of binary Bose-Einstein condensates confined in a twisted two-dimensional harmonic trap. The ground state hosts three distinct phases-phase-separated, polarized, and phase-mixed-determined by the Rabi coupling, interaction ratio G (between inter-component and intra-component interactions), and spin-twisting which induces edge-localized polarization through position-dependent detuning. In the phase-mixed regime, the ground state is characterized by a finite spin-scissors susceptibility and can be accurately described using local density approximation. In the dynamics, the system exhibits stable periodic beating in the phase-mixed state for G<1. For G>1, its evolution progresses from beat damping (phase-separated state) to polarized relaxation (polarized state), finally reaching stable periodic beating (phase-mixed state) after a finite waiting time. The dependence of the waiting time contrasts sharply with the monotonic behavior of one-dimensional spin-dipole dynamics, revealing qualitatively distinct mechanisms governed by geometry and interactions. In summary, these results establish a unified link between ground-state properties and nonequilibrium responses in twisted spinor condensates, offering a versatile platform for exploring spin-related quantum many-body phenomena.

Simulating out-of-equilibrium dynamics of quantum field theories in nature is challenging with classical methods, but is a promising application for quantum computers. Unfortunately, simulating interacting bosonic fields involves a high boson-to-qubit encoding overhead. Furthermore, when mapping to qubits, the infinite-dimensional Hilbert space of bosons is necessarily truncated, with truncation errors that grow with energy and time. A qubit-based quantum computer, augmented with an active bosonic register, and with qubit, bosonic, and mixed qubit-boson quantum gates, offers a more powerful platform for simulating bosonic theories. We demonstrate this capability experimentally in a hybrid analog-digital trapped-ion quantum computer, where qubits are encoded in the internal states of the ions, and the bosons in the ions' motional states. Specifically, we simulate nonequilibrium dynamics of a (1+1)-dimensional Yukawa model, a simplified model of interacting nucleons and pions, and measure fermion- and boson-occupation-state probabilities. These dynamics populate high bosonic-field excitations starting from an empty state, and the experimental results capture well such high-occupation states. This simulation approaches the regime where classical methods become challenging, bypasses the need for a large qubit overhead, and removes truncation errors. Our results, therefore, open the way to achieving demonstrable quantum advantage in qubit-boson quantum computing.

The paper presents a survey of some dynamical transitions in nonequilibrium trapped Bose-condensed systems subject to the action of alternating fields. Nonequilibrium states of trapped systems can be realized in two ways, resonant and nonresonant. Under resonant excitation, several coherent modes are generated by external alternating fields, whose frequencies are tuned to resonance with some transition frequencies of the trapped system. A Bose system of trapped atoms with Bose-Einstein condensate can display two types of the Josephson effect, the standard one, when the system is separated into two or more parts in different locations or when there are no any separation barriers, but the system becomes nonuniform due to the coexistence of several coherent modes interacting with each other, which is termed internal Josephson effect. The mathematics in both these cases is similar. We concentrate on the internal Josephson effect. Systems with nonlinear coherent modes demonstrate rich dynamics, including Rabi oscillations, Josephson effect, and chaotic motion. Under Josephson effect, there exist dynamic transitions that are similar to phase transitions in equilibrium systems. The bosonic Josephson effect is shown to be realizable not only for weakly interacting systems, but also in superfluids, with not necessarily weak interactions. Sufficiently strong nonresonant excitation can generate several types of nonequilibrium states comprising vortex germs, vortex rings, vortex lines, vortex turbulence, droplet turbulence, and wave turbulence. Nonequilibrium states can be characterized and distinguished by effective temperature, effective Fresnel number, and dynamic scaling laws.

It is challenging to build scalable quantum processors capable of both parallel control and local operation. As a promising platform to overcome this challenge, optical lattices offer exceptional parallelism. However, it has been struggling with precise local operations due to relatively narrow lattice spacings. Here, we introduce a new quantum processor incorporating orbit-qubit encoding and internal states (as auxiliary degrees of freedom) to achieve spatially selective operations together with parallel control. With this processor, we generate one-dimensional and two-dimensional cluster states using minimal layers of controlled-Z gates. We experimentally detect the multipartite entanglement of a two-dimensional cluster state involving 123 orbit qubits through direct stabilizer measurements, verifying the full bipartite non-separability. Furthermore, we demonstrate measurement-based quantum computation by implementing single-qubit and two-qubit logical gates, highlighting the flexibility of orbit-qubit operations. Our results establish orbit-qubit optical lattices as a scalable quantum processing architecture, opening new pathways for quantum computation applications.

We propose a field-theoretic framework for ideal hydrodynamics of charged relativistic fluids formulated in terms of a complex scalar field defined on a spacelike hypersurface comoving with the fluid. In the normal phase, the dynamics of charge-carrying fluids is constrained by the restrictive chemical shift symmetry, which locks charges to fixed positions in the comoving plane as they are transported through space by the fluid's motion. On the other hand, in the superfluid phase, the chemical shift symmetry is relaxed to a constant shift, allowing charges to redistribute freely across the comoving hypersurface. We demonstrate that both models recover the respective nonlinear hydrodynamic equations and provide explicit expressions for the collective variables of hydrodynamics in terms of the theory's fields. Introduced models provide a UV completion to the effective field theories of hydrodynamics constructed in terms of the Goldstone fields. Finally, we propose a relativistic fracton fluid phase as a natural interpolation between the normal and superfluid phases, in which the mobility of elementary charges is constrained by a linear shift symmetry in the comoving space.

As in the gravitational context, one of the most challenging open question in analogue black holes formed in Bose-Einstein condensates concerns the backreaction of Hawking-like radiation on the condensate and its subsequent evolution. In this work we derive the basic equations describing this backreaction within the density-phase formalism, which avoids infrared divergences and is particularly well suited to one-dimensional configurations.

We study the impact of weak quenched disorder on the BCS-BEC crossover in the Hubbard model within a functional-integral framework. By deriving the thermodynamic potential up to second order in both the disorder potential and pairing fluctuations, we obtain self-consistent expressions for the number equation, condensate fraction, superfluid fraction and sound speed at zero temperature. In the dilute BEC limit, our results analytically reproduce the known continuum limits of weakly interacting bosons, where weak disorder depletes the superfluid more strongly than the condensate due to broken translational symmetry, and enhances the sound speed through the overcompensation of the static compressibility. These findings establish a unified and controlled framework for describing the BCS-BEC crossover in disordered lattice models, and they provide a foundation for future extensions to finite temperatures and multiband Hubbard models.

The synthetic frequency dimension, which realizes fictitious spatial dimensions from the spectral degree of freedom, has emerged as a promising platform for engineering artificial gauge fields in studying quantum simulations and topological physics with photons. A current central task for frequency-domain photons is the creation and manipulation of nontrivial non-Abelian field strength tensors and observing their governing dynamics. Here, we experimentally demonstrate a miniaturized scheme for creating non-Abelian electric fields in a photonic frequency chain using a polarization-multiplexed, time-modulated ring resonator. By engineering spin-orbit coupling via modulation dephasing, polarization rotation, and polarization retardation, we achieve programmable control over synthetic Floquet bands and their quasimomentum spin-resolved textures. Leveraging self-heterodyne coherent detection, we demonstrate Zitterbewegung -- a trembling motion of photons -- induced by non-Abelian electric fields on the frequency chain. We further observe the interference between Zitterbewegung and Bloch oscillations arising from the coexistence of non-Abelian and Abelian electric fields. Our work bridges synthetic dimensions with non-Abelian gauge theory for versatile photonic emulation of relativistic quantum mechanics and spinor dynamics, and can be instrumental in applications like frequency-domain optical computation and multimodal frequency comb control.

We give a bosonic realization of the SU(3) antiferromagnetic Heisenberg (AFH) chain in the alternating conjugate representation, and study its phase diagram as a function of staggered interactions and anisotropy along the $T^3$ and $T^8$ directions. Unlike the SU(2) case, we observe a chiral-reversed quantum phase transition, where each competing phase is adiabatically connected to one of the chiral Haldane phases predicted in the SU(3) AFH chain with local adjoint representation. In the vicinity of the Heisenberg point, we identify a symmetry-protected topological state that appears at the first excited energy level. We also study the spontaneous $\mathbb{Z}_3$ symmetry breaking of the system, and provide a variational wavefunction that captures the transition from the topological phase to the trivial phase. Finally, we propose an experimental realization of our bosonic model by two spin-1/2 bosons in an optical lattice.

Kinetic constraints are generally expected to slow down dynamics in many-body systems, obstructing or even completely suppressing transport of conserved charges. Here, we show how gauge theories can defy this wisdom by yielding constrained models with faster-than-diffusive dynamics. We first show how, upon integrating out the gauge fields, one-dimensional U(1) lattice gauge theories are exactly mapped onto XX models with non-local constraints. This new class of kinetically constrained models interpolates between free theories and highly constrained local fermionic models. We find that energy transport is superdiffusive over a broad parameter regime. Even more drastically, spin transport exhibits ballistic behavior, albeit with anomalous finite-volume properties as a consequence of gauge invariance. Our findings are relevant to current efforts in quantum simulations of gauge-theory dynamics and anomalous hydrodynamics in closed quantum many-body systems.

We extend the Gross-Pitaevskii equation to incorporate the effect of quantum fluctuations onto the flow of a weakly interacting Bose-Einstein condensate. Applying this framework to an analogue black hole in a quasi-one-dimensional, transonic flow, we investigate how acoustic Hawking radiation back-reacts on the background condensate. Our results point to the emergence of stationary density and velocity undulations in the supersonic region (analogous to the black hole interior) and enable to evaluate the change in upstream and downstream Mach numbers caused by Hawking radiation. These findings provide new insight into the interplay between quantum fluctuations and analogue gravity in Bose-Einstein condensates.

Understanding flux string dynamics can provide insight into quark confinement and hadronization. First-principles quantum and numerical simulations have mostly focused on toy-model Abelian lattice gauge theories (LGTs). With the advent of state-of-the-art quantum simulation experiments, it is important to bridge this gap and study string dynamics in non-Abelian LGTs beyond one spatial dimension. Using tensor network methods, we simulate the real-time string dynamics of a $2\!+\!1$D SU$(2)$ Yang--Mills LGT with dynamical matter. In the strong-coupling regime and at resonance, string breaking occurs through sharp Casimir reduction along with meson and baryon-antibaryon formation, a distinctively non-Abelian feature. At finite baryon density, we discover a \textit{baryon blockade} mechanism that delays string breaking. Away from resonance, the magnetic term drives purely non-Abelian fluctuations: glueball loops and self-crossed strings that resolve two SU$(2)$ intertwiners with distinct dynamics. For higher-energy strings, we uncover representation-dependent tension-reduction resonances. Our findings serve as a guide for upcoming quantum simulators of non-Abelian LGTs.

We implement an analogue Wheeler-DeWitt mini-universe constituted by a well-isolated atomic Bose-Einstein condensate in a time-independent conservative potential. In exact analogy with the Wheeler-DeWitt framework for the actual universe, our system has the fundamental problem of defining from within a meaningful time variable over which to order the events. Here, we partition the mini-universe into a bright and a dark sector, enabling entropy exchange between them through a potential barrier. We show that the Hamiltonian of the condensate in the bright sector is analogous to the one in canonical minisuperspace models. We define an entropic time and show with experimental data that it is robustly monotonic even when the bright sector undergoes several cycles that begin with a 'big bang' and end with a 'big crunch'. By tuning the barrier height, we control the rate of entropy production and thus the speed of the emergent entropic time and the dynamics of the bright universe. We finally derive an entropic time-dependent Schroedinger equation that could be considered as a generalization of the standard one, and use it to reproduce our data. This work experimentally validates the proposition that time in quantum cosmological models may not be fundamental, but instead emerges from thermodynamic gradients, while establishing a concrete experimental platform for evaluating several aspects of quantum gravity theories.

We demonstrate a continuous spectroscopy technique with frequency sensitivity well below the natural transition linewidth, while maintaining a locking range hundreds of times larger. The method exploits the position dependence of a continuous, broadband magneto-optical trap operating on the 7.5 kHz-wide intercombination line of strontium. We show that the frequency sensitivity is fundamentally insensitive to the effective MOT laser linewidth. By applying active feedback on the MOT position to a dispersion-optimized frequency comb, which serves as the reference for stabilizing the MOT laser, we achieve a frequency instability below $4.4\times10^{-13}$ after 400 s of averaging in both the optical and radio-frequency domains, surpassing the performance of conventional hot-vapor modulation transfer spectroscopy. Our method is a broadly applicable alternative route to long-term frequency references in the high $10^{-14}$ range around 100 s.