Papers

Entanglement dynamics depend not only on how a quantum system is partitioned, but critically on how interactions across that partition are structured. For a spatial bipartition of a locally interacting system, entanglement is generated near the boundary and then propagates into the bulk. By contrast, when two extended quantum fields are coupled locally along their entire length, the interaction crosses the field-space partition everywhere, and this generates correlations throughout the system. Here, we study the entanglement dynamics between two gapless one-dimensional quantum many-body systems described by Luttinger liquid theory. The systems are initially decoupled and prepared at zero or finite temperature, after which a time-dependent tunneling interaction is activated uniformly along their length. Within a Gaussian approximation, we derive general analytical expressions for the logarithmic negativity, mutual information, and Rényi entropies under arbitrary coupling protocols. At zero temperature, entanglement displays an early-time power-law growth whose exponent is fixed solely by the first non-null derivative of the tunneling protocol. Once the coupling saturates, we obtain exact long-time averages of the information-theoretic quantities and characterise how the correlations scale with temperature and the final coupling strength. We also analyse how mutual information and logarithmic negativity approach the adiabatic limit for a very slow protocol with respect to intrinsic system timescale. This work extends the study of entanglement dynamics in nonequilibrium field theory to field-space partitions and mixed initial states.

We investigate the entanglement entropy of a class of $N$-particle excited state of fermions confined in a two-dimensional harmonic trap rotating at an angular frequency $\Omega$. The excited state is constructed by filling a particular set of $N$ single-particle energy levels. We analytically compute the Rényi entropies of order $q$ in a disc of radius $r$ around the centre of the trap, and the cumulants corresponding to number fluctuations of fermions within the disc. We found that the area law scaling of entanglement entropy holds even for a class of excited states. We also verified the well-known series expansion of entanglement entropy in terms of the particle number cumulants for non-interacting fermions. We further derive the centered cumulant generating function, demonstrating that the associated probability distribution function in the disc has identical scaling properties, up to a variable shift, to the known ground-state result. Finally, we extend our result to an annular region, showing that both the Rényi entropy and particle number cumulants of the annulus decompose into sums of the corresponding quantities of the two bounding discs. These additive relations hold as long as the width of the annulus is sufficiently large.

We report the observation of large anomalous shifts, up to +7.5 G, of the positions of the 33.6 G and 39.9 G Feshbach resonances in potassium-39 atoms confined in a 1063.9 nm optical dipole trap (ODT) at temperatures up to around 35 {\mu}K and trap depths up to about 136 {\mu}K. When the atom cloud is cooled to lower temperatures, by reducing the trap depth of the ODT, the shifts decrease proportionally with trap depth and approach zero at zero depth. We show that the large observed shifts originate from a large differential ac Stark shift between the incoming pair of potassium-39 atoms and the weakly bound Feshbach molecule, which in turn originates from an unexpectedly large dynamic polarizability of the Feshbach molecule. The polarizabilities of the Feshbach molecules extracted from the measured shifts of the 33.6 G and 39.9 G resonances are about four to seven times the sum of the polarizabilities of the two incoming potassium-39 atoms, that is, about four to seven times larger than the usual polarizability of weakly bound Feshbach molecules. The large polarizabilities of the Feshbach molecules are attributed to a near-coincidence between the frequency of the 1063.9 nm ODT laser and the frequency of a molecular transition from the last vibrational level of the lowest triplet a3{\Sigma}+u potential to a vibrational level of the excited b3{\Sigma}+g potential. Other potassium-39 Feshbach resonances we have studied exhibit a zero or very small shift, corresponding to molecular polarizabilities close to the sum of the polarizabilities of the two incoming potassium-39 atoms.

We adapt tools from the theory of quantum random walks to investigate slow relaxation dynamics through the many-body state graph. Specifically, we construct a probe of heterogeneity between basis states defined using hitting times derived from the unitary time-evolution operator. We find that the state-graph geometry, encoded by the pairwise hitting time of basis states, is a highly sensitive indicator of slow relaxation dynamics in a variety of systems. We study three paradigmatic models: the Rosenzweig-Porter model, the quantum East model, and the triangular lattice gas model, exhibiting a sudden onset of slow dynamics upon tuning of a control parameter. As a global characterization of the graph geometry, we analyze the spectral radius of the hitting matrix. We find that it increases sharply at the onset of slow dynamics, spanning many orders of magnitude, with a characteristic crossing point at the transition. Our work provides a geometric framework for describing and identifying phases with slow relaxation using a unified graph-theoretic formalism.

One-dimensional strongly repulsive quantum mixtures exhibit non-trivial exchange statistics arising from the interplay between orbital and spin degrees of freedom. Using an exact solution, we demonstrate that in the single-impurity limit the wavefunction of Fermi-Fermi and Bose-Bose mixtures on a ring threaded by an artificial gauge field display exact anyonic statistics under exchange of the impurity with the majority particles. The resulting fractional exchange phase is fixed by the angular momentum sector selected through the applied flux. We find that the impurity momentum distribution coincide exactly with the anyonic one, independently of the bosonic or fermionic nature of the mixture, with the large-momentum tails encoding a direct signature of the anyonic statistical angle. Finally, we devise a quench protocol for reversible dynamical anyonization. Our results provide a path for realizing and manipulating anyonized states with ultracold atoms.

Classical waves in dispersive media naturally exhibit distinct phase and group velocities. Whether an analogous separation of velocities can emerge in matter waves under strong many-body interactions has remained experimentally unexplored. Here, we demonstrate the emergence of a velocity hierarchy in a strongly interacting lattice gas. Using quench spectroscopy together with time-resolved correlation measurements, we independently determine the sound, group, and phase velocities across the superfluid-to-Mott-insulator transition. These velocities are nearly degenerate close to the transition, but progressively separate as the Mott gap opens and the quasiparticle dispersion acquires a massive relativistic-like form. Strikingly, phase-coherence fronts propagate faster than the Lieb-Robinson velocity scale while remaining fully consistent with locality. The measured velocities satisfy a relativistic-like invariance relation in the insulating regime. Our results establish propagation-velocity hierarchies as emergent signatures of strongly correlated quantum dynamics.

We study the ultimate quantum limit of rotation sensing with a few strongly interacting bosons confined in a quasi-one-dimensional ring trap with two weak links. It is demonstrated that a self-consistent many-body solution of the problem is required to correctly predict the ultimate sensitivity of this strongly correlated many-body gyroscope to rotation. For both small rotation velocities and small particle numbers, the many-body quantum Fisher information becomes maximal for large interaction couplings, showing the potential of strongly interacting miniaturized many-body sensors to precisely estimate slow rotations with high spatial resolution.

The recent discovery of several families of kagome materials and experimental realization of optical kagome lattices have stimulated growing numerical studies of interaction-driven correlated states on the kagome lattice. Among the available numerical approaches, determinant quantum Monte Carlo (DQMC) is a powerful method for investigating such strongly correlated states. However, the accessible system sizes of existing DQMC simulations remain limited, preventing reliable finite-size scaling analyses. Here we develop a general acceleration scheme based on fast Fourier transform (FFT) for propagator multiplications on composite lattices and combine it with the delay-update algorithm, enabling simulations on system sizes twice as large as those of previous DQMC studies, allowing reliable finite-size scaling analyses of the attractive kagome-lattice Hubbard model. Our large-scale simulations reveal the interaction-driven zero-temperature superfluid quantum criticality at the Dirac filling and provide reliable estimates of the associated critical exponents. Besides, we find no evidence that the previously proposed triangle-rule charge-density-wave order survives in the thermodynamic limit, suggesting that it is likely a finite-size effect. Moreover, for system sizes accessible in current two-dimensional optical lattice experiments, the combined FFT and delay-update scheme exhibits an effective computational cost scaling as $N^{2.49}$, substantially below the $\mathcal{O}(N^3)$ computational cost of conventional DQMC simulations.

Two-dimensional quantum gases provide a distinctive setting in which enhanced thermal fluctuations, finite-size effects, and two-body bound-state formation are intrinsically intertwined. In this work, we study a single attractive impurity immersed in a finite, weakly interacting two-dimensional Bose gas, where finite size stabilizes a nonzero condensate fraction by introducing an infrared momentum scale, thereby enabling a Bogoliubov description of the bath. Using a hybrid approach that combines finite-temperature many-body scattering theory with input from path-integral Monte Carlo, we analyze the impurity quasiparticle energy across the condensate and normal regimes. The infrared scale generates a phonon-activation temperature below which the impurity energy remains nearly temperature independent. Once the resolved phonon modes become thermally populated, their contribution competes with condensate depletion, producing a nonmonotonic temperature dependence of the polaron energy. These results suggest that attractive Bose polarons may serve as sensitive probes of finite-size thermal fluctuations, phonon dressing, and bound-state physics in low-dimensional Bose gases.

We investigate Faraday instabilities in three-dimensional parametrically driven trapped dipolar quantum gases within the framework of the extended Gross-Pitaevskii equation. In the superfluid regime, periodic modulation of the short-range interactions induces the resonant excitation of discrete polygonal surface modes displaying sub-harmonic response. The resonance frequencies and the parametric windows of instability are accurately captured by an appropriate Mathieu equation, further corroborating our simulations. It is also shown that dipolar interactions in superfluids shift the resonance frequencies to lower values compared to their non-dipolar counterparts. Near the superfluid-to-supersolid transition, intrinsic density undulations associated with the softened roton mode accelerate pattern formation, yielding hybrid surface and bulk excitations. The bulk patterns prevail deeper in the supersolid regime, emerging from droplet collisions with the superfluid background. Our results reveal a crossover from surface collective modes to hybrid surface-bulk excitations and demonstrate how parametric driving dictates pattern formation in long-range interacting quantum fluids.

High-power laser sources that exhibit low relative intensity noise and allow simultaneous dynamic control of their light level are required for a broad range of applications in various fields of physics. Acousto-optic modulators (AOMs) are widely used for active power stabilization and regulation due to their simple drive electronics requirements and high optical power handling capability in free-space. However, the rather slow propagation speed of the sound wave within the AOM crystal typically limits their control bandwidth to a few hundred kHz. In this work, we present a novel AOM-based control system that is capable of significantly suppressing intensity noise of high-power lasers up to the MHz range. By combining two standard feedback loops with one feedforward control branch and optimizing the beam path in the AOM crystal, ultra-low relative intensity noise levels down to $-155\, \text{dB}\,\text{Hz}^{-1}$ even at several hundred kHz are achieved. Our results are relevant for applications that require ultra-low intensity noise at Fourier frequencies up to the MHz range, such as optical lattice experiments with light ultracold atoms.

Gauge invariance is central to modern physics and underpins quantum simulations of lattice gauge theories (LGTs). Existing quantum simulation approaches employ Gauss's law either to energetically suppress gauge-violating processes in analog platforms or to detect and discard gauge-violating outcomes in digital devices. Here we introduce a third paradigm, in which Gauss's law is used to dynamically generate the gauge theory itself from a substantially simpler Hamiltonian. Starting from a readily programmable three-body XXX model, we employ experimentally efficient single-qubit U(1) gauge symmetry-generator terms that induce the dynamics of a U(1) LGT. We implement this approach using 101 qubits on a 156-qubit IBM quantum processor and observe real-time dynamics in quantitative agreement with the target LGT while reducing the entangling-gate depth per Trotter step by a factor of five compared with a direct implementation. Our results establish gauge protection as a resource for Hamiltonian engineering rather than merely symmetry preservation, opening a scalable resource-efficient route towards digital quantum simulations of increasingly complex gauge theories in higher spatial dimensions.

Recent experiments in moiré materials have enabled the realization of a variety of exotic quantum phases. In this context, the Hofstadter-Hubbard model has been proposed as a possible setting for hosting chiral spin liquid. Concurrently, significant progress has been recently made in the computational methods for two-dimensional many-body fermion systems, which makes numerically studying this challenging model a real possibility in genuine 2D geometry. Motivated by these advances, we investigate the putative chiral spin liquid phase in the triangular-lattice Hofstadter-Hubbard model using variational Monte Carlo with neural quantum states (NQS) and projected entangled pair states (PEPS). We observe spin-charge separation directly in real space through numerical spin-pumping simulation and real-time spin and charge motion. In addition, in the context of anyonic superconductivity conjectured in this model, we find a positive two-electron binding energy on small systems, but it decreases below our numerical resolution as the system size increases. Our work demonstrates NQS and PEPS as powerful tools, capable of cross-checking each other, for diagnosing topological order and fractionalized excitations in strongly correlated electronic systems.

We investigate the validity of effective one-dimensional models for quasi-1D fermionic systems by benchmarking a coupled-channel approach against the exact low-energy theory derived from the underlying three-dimensional problem in the weakly and strongly attractive limits. We show that reproducing the exact two-body scattering amplitude is insufficient to construct the correct effective low-energy theory of quasi-1D fermions. In the weakly attractive regime, it does not capture the emergent three-body interaction induced by transverse excitations. In the strongly attractive regime, it yields an atom-dimer scattering length with an incorrect dependence on the three-dimensional scattering length. These results demonstrate the limitations of effective one-dimensional descriptions based solely on two-body physics and highlight the need to explicitly include few-body correlations in quasi-1D systems.

We study the spin-partitioned entanglement Hamiltonian of the one-dimensional repulsive Hubbard model. By combining bosonization with exact diagonalization, we show that tracing out one spin species produces a modular Luttinger liquid, whose properties fundamentally differ from those of the physical system. While the modular spectrum is fully dispersionless and possesses a momentum-independent entanglement gap, its eigenstates exhibit algebraic correlations governed by a single effective Luttinger parameter equal to the geometric mean of the charge and spin Luttinger parameters. The resulting entanglement spectrum displays a universal branching hierarchy in excellent agreement with exact diagonalization. We further demonstrate that the modular ground state is nearly identical to that of the spinless Luttinger liquid. These results uncover a universal modular structure in interacting one-dimensional fermionic systems.

One-dimensional (1D) superlattices provide one simplified platform for exploring moiré physics from a low-dimensional perspective, with the ratio of lattice constants playing a role analogous to the twist angle in two-dimensional bilayers. Here, we propose a 1D $\mathcal{PT}$-symmetric bichromatic optical lattice for a weakly repulsive Bose-Einstein condensate and investigate how the interplay of dissipation and interaction impacts the lowest moiré flat band. Without interaction, we find that the lowest-band flatness induced by commensurate ratios exhibits a parity-dependent response to the $\mathcal{PT}$-symmetric imaginary potential due to the distinct $\mathcal{PT}$ pairing mechanism for the energy spectrum. For ratios with even denominators (i.e., even parities), the level attraction and thus the $\mathcal{PT}$-symmetry breaking occur within the lowest two bands, leading to a monotonic broadening of the lowest flat band, whereas odd denominators (i.e., odd parities) yield a nonmonotonic response due to the $\mathcal{PT}$-symmetry breaking within the second and the third lowest bands instead while the lowest band remains purely real. This parity-dependent phenomenon can be understood by the perturbation theory. Furthermore, by solving the Gross-Pitaevskii equation, we also find that although the weak repulsive interaction can broaden the moiré bands alone, the combined effects of interaction and imaginary potential also lead to parity-dependent behaviors. For even parities, band flattening is consistently diminished, whereas for odd parities, the imaginary potential can either enhance or reduce the degree of flattening. These results pave the way for experimental studies of dissipation and interaction effects on band flatness in moiré systems.