2026-07-30 | | Total: 7
Ultracold polar molecules provide strong, long-range interactions that microwave shielding makes tunable and nearly nuclear-spin independent, giving an emergent SU($N$) symmetry. However, extended Hubbard models of polar molecules in optical lattices lack, so far, controllable finite on-site interactions, a key ingredient of strong correlated physics. We show that tuning the Rabi frequency of the microwave coupling can bring two individual molecules (monomers) on neighboring lattice sites into resonance with a field-linked dimer (doublon) on one of the sites, enabling coherent doublon--monomer-pair conversion. In this model, we characterize the key Hubbard parameters and the dimer lifetime, demonstrating that the on-site and off-site interactions can be tuned nearly independently through the microwave amplitude and orientation, respectively. Our results provide a roadmap for implementing SU($N$)-symmetric extended Hubbard models with controllable doublon fluctuations, providing access to quantum-simulation in the strongly dipolar regime.
Using a non-perturbative quantum kinetic framework, we develop a unified description of the far-from-equilibrium dynamics of three-dimensional Bose gases following cooling quenches across the Bose-Einstein condensation transition. By tracking the spatio-temporal evolution of the momentum distribution, we show that the equilibrium condensation threshold simultaneously acts as a dynamical critical point, separating distinct far-from-equilibrium universality classes governed by different non-equilibrium attractors. While quenches above the transition exhibit a single-timescale relaxation toward a thermal fixed point, quenches below the transition display a crossover from a transient weak-turbulence regime to a coarsening fixed point governed by the diffusive recombination of vortex lines. Quenches directly to the condensation threshold, finally, are controlled by a previously unidentified critical fixed point characterized by the superdiffusive spreading of critical fluctuations and a distinct set of dynamical exponents. Together, these dynamical scaling laws establish a far-from-equilibrium counterpart of the condensation phase transition, in which the equilibrium critical point also organizes the long-time non-equilibrium dynamics.
We characterize the Hall response of non-interacting fermionic $M$-leg ladders in the presence of an artificial magnetic flux, that can be realized in one-dimensional optical lattices supplemented with a synthetic dimension. We focus on the Hall imbalance, which can be directly measured in experiments with ultracold fermionic atoms. At relatively large synthetic flux we find a dependence of the density that contrasts with previously reported density-independent behavior. In particular, the Hall imbalance can be significantly enhanced in the limit of small density of particles (or holes), or it can vanish and change sign at specific fluxes, which we obtain analytically in the large inter-leg coupling regime. This behavior is explained in terms of Lifshitz transitions of the band structure, where the number of Fermi points changes as a function of the flux and density. Finally, we explore the connection between these transitions and the so-called Meissner-vortex transition for fermionic ladders by computing the site-resolved leg and rung currents and discussing the similarities and differences with the bosonic counterpart.
A recent \texttt{Google Quantum AI} experiment [\href{https://www.science.org/doi/10.1126/science.adr9680}{Gyawali \textit{et al.}, Science \textbf{393}, 71 (2026)}] has exploited quantum parallelism to emulate disorder-averaged many-body dynamics, with conserved local degrees of freedom generating an effective disorder potential. We investigate how the local spectrum of these static variables controls localization in a flavor-extended ${\mathbb Z}_2$ lattice gauge theory, which maps onto a mixed-field Ising chain with $n$-level bond disorder. Combining finite-size spectral and entanglement diagnostics with infinite matrix-product state dynamics, we find a qualitative distinction between binary and multilevel disorder. For $n=2$, apparent localization ultimately gives way to thermalization; the long-lived transient arises from energy-scale separation, degenerate spectral towers, and approximate Hilbert-space fragmentation. By contrast, $n=4$ displays consistent localization signatures, including Poissonian level statistics, area-law eigenstate entanglement, nonthermal entanglement spectra, and persistent local memory over accessible times in the thermodynamic limit. Our results show that, despite its larger variance, binary disorder lacks the local amplitude diversity needed to suppress resonances. Thus, localization is governed not simply by disorder strength, but by the local disorder spectrum and the resulting resonant connectivity of the many-body Hilbert space.
Simulating the real-time dynamics of particle collisions is a promising application of quantum simulators, because classical methods such as tensor networks struggle to capture the highly entangled states generated in high-energy scattering. Realizing such simulations requires both the preparation of incoming wave packets and the detection of the outgoing scattering products. In this work, we propose protocols that address both challenges on programmable analog and digital quantum simulation platforms. Our state-preparation scheme uses a weakly coupled auxiliary qubit - or, more generally, a customized local quench - to inject a single quasiparticle with well-defined momentum. Because it relies only on conservation of energy, this scheme requires no fine-tuning or prior knowledge about particle eigenstates, making it robust against errors in calibration and implementation. The momenta of scattering products are then extracted, using only local measurements, from the interference pattern that arises when particles are reflected at the system's boundary. We validate our protocols through numerical simulations, first in a simple single-particle model and subsequently in two interacting many-body systems: a Rydberg atom chain and an Ising chain in a mixed field. We demonstrate how high-energy regimes, necessary to access inelastic scattering processes, can be reached through an adiabatic ramp, and how the wave packet shape can be optimized by spatially modulating the Hamiltonian. Finally, we show how the protocol can be generalized to systems with more than one spatial dimension. Our proposal provides a versatile approach to the quantum simulation of scattering phenomena, and is compatible with several quantum simulation platforms that are already experimentally available.
Non-Abelian gauge theories with dynamical matter govern the strong interaction and a broad class of strongly correlated quantum systems, yet their ground-state properties remain difficult to obtain from first principles. Using a continuous-group variational Monte Carlo approach that retains the full SU$(2)$ gauge field without truncation, we determine the ground-state behavior of the SU$(2)$ lattice gauge theory with staggered fermions on an $L\times L$ square lattice. Treating the magnetic and electric couplings $λ$ and $g^2$ independently, we find a magnetic-flux transition at $λ^\ast=-0.040\pm 0.005$, with no resolvable drift of the transition point as the electric coupling is varied. Along the physical coupling line $λ=4/g^2$, for $L=4,6,8$, we uncover a gauge-matter delocalization crossover from a flux-disordered regime at strong electric coupling to an ordered unity-flux regime at weak coupling. The chiral condensate, a gauge-invariant Wilson-line meson correlator, and the local color density consistently reveal the emergence of coherent gauge-assisted matter dynamics. Together, these results provide a unified physical picture of how magnetic-flux ordering and fermionic coherence develop in an untruncated non-Abelian lattice gauge theory.
Determining the ground state of non-Abelian lattice gauge theories coupled to dynamical fermions is key to understanding confinement and the phase structure of gauge--matter systems. We present a variational Monte Carlo framework for the ground state of the untruncated fully-continuous SU$(2)$ lattice gauge theory coupled to dynamical staggered fermions on an $L\times L$ square lattice. We work in the magnetic basis with a neural-network representation of the gauge wavefunction. The fermions are described by a gauge-covariant Gaussian fermionic correction built on a fixed Néel reference state where, for each sampled gauge configuration $\mathbf{U}$, the correction is generated by a Hermitian operator. This operator is constructed from short Wilson lines and the eigenvectors of the mass--hopping Hamiltonian, with number of variational parameters polynomial in the system size. This Gaussian structure also gives analytical expressions for all fermionic contributions to the energy and related observables in terms of the fermion occupation matrix. The results are validated against strong-coupling perturbation theory, where they recover the expected effective antiferromagnetic spin Hamiltonian. Using this framework, we map a coarse ground state phase diagram in the plane of independent electric and magnetic couplings $(g^2, λ)$ and show that a hysteresis analysis can identify the existence of phase transitions. Restoring the physical relation $λ=4/g^2$, we characterize how increasing the system size and changing the electric coupling $g^2$ move the state away from the reference Néel state, for lattice sizes $L=4,6,8$. More broadly, the method offers a sign-problem-free variational framework for continuous non-Abelian gauge groups with dynamical matter that should extend to other matter content and higher-dimensional lattices.