Papers

We study symmetric gapped ground states of one-dimensional spin $S$ chains with a continuous exponential symmetry $U(1)\rtimes\mathbb Z_2$ with integer base $q$, where the U(1) charge has a modulated unit of charge $q^{L-j}$ decaying exponentially with site $j$. The simplest model is a quantum breakdown model that can be viewed as a $q$-deformation of the XXZ Haldane spin chain. We show that the topology of their symmetric ground states crucially depends on $S$ and $q$. For integer $S$, odd $q$ allows a nontrivial $\mathbb Z_2$ Haldane-type symmetry-protected topological (SPT) phase with protected edge modes, while even $q$ allows no nontrivial SPT phase. For half-odd integer $S$, odd $q$ prohibits any unique symmetric gapped ground state due to a Lieb-Schultz-Mattis obstruction, whereas even $q$ enforces a novel symmetric gapped phase which we name as the “half-SPT" phase. Unconventionally, the half-SPT phase only has a symmetry protected edge mode on one edge, and has half of the string orders of an SPT phase. Our results uncover a new type of symmetry-protected gapped phases beyond existing classification schemes.

We have presented an analytical and numerical study of efficiency fluctuations in quantum Carnot and Otto engines with pistons operating under non-quasi-static adiabatic expansion and compression processes at constant speed, along with quasi-static isothermal (for the Carnot engine) or isochoric processes (for the Otto engine). We have considered the two-point measurement scheme during the adiabatic strokes in this regard. By considering each of the pistons as a single quantum particle in a 1-D box with a non-quasi-statically movable boundary, we capture the stochastic nature of the work output resulting from the unitary evolution of the quantum mechanical state of the system during non-equilibrium processes. Furthermore, we have explicitly shown how the piston speed during quantum non-equilibrium strokes introduces additional quantum fluctuations that modify the average efficiency of both the engines.

A major application of analog quantum hardware is quantum simulation, where a system evolves over time according to a given Hamiltonian. One challenge when implementing the Hamiltonian on a programmable quantum simulator is the availability of distinct coupling types, many of which go beyond the conventional ZZ-couplings found in Ising Hamiltonians that are naturally implemented by various platforms. In prior work, we have developed a method based on domain wall encoding to simulate one-dimensional spin chains using only Ising Hamiltonians, which would otherwise require XX+YY-couplings as well. Here, we implement the method on a commercially accessible Rydberg atom quantum device, and simulate the time evolution of the Heisenberg XXZ model with various anisotropies. We successfully probe the XY and Ising phases of the model and reproduce the qualitative behavior of several one- and two-body observables in both phases, as well as at the critical point. Our work augments the toolbox to simulate canonical quantum many-body systems on analog quantum simulators.

The stability of the Berezinskii–Kosterlitz–Thouless phase against long-range interactions remains an open problem. We address this question in the two-dimensional XY model with couplings decaying as $r^{-2-\sigma}$, by developing a real-space renormalization-group treatment near the Gaussian manifold of quasi-long-range ordered (QLRO) states. We retain the full many-body vortex interaction and show that, at leading order in the fugacity, screening by small vortex–antivortex pairs preserves its Gaussian-average representation without an expansion in the long-range coupling. Including spin waves perturbatively then yields a local stability criterion: a stable QLRO interval exists only for $\sigma>7/4$. We also investigate a Wilson–Fisher-like fixed point at finite coupling for $\sigma<7/4$, for which the retained flow predicts $\nu^{-1}\sim\sqrt{7/4-\sigma}$ as it merges with the Gaussian line.

Microwave shielding has recently emerged as a powerful tool for engineering interactions in ultracold polar molecules, yet its potential for controlling cavity quantum electrodynamics remains largely unexplored. Here, we investigate a molecular cavity quantum electrodynamics platform in which two microwave-shielded polar molecules are coupled to a single optical cavity mode and demonstrate that shielding-induced interactions provide a unified mechanism for both photon blockade and Bell-state protection. The anisotropic interaction reshapes the few-excitation spectrum by enhancing its anharmonicity, thereby suppressing multiphoton transitions and improving the single-photon purity by more than three orders of magnitude. The enhanced blockade is accompanied by the emergence of negative longitudinal spin correlations, revealing the interaction-induced suppression of simultaneous molecular excitations. We further show that the photon statistics are highly sensitive to the relative molecular configuration, with positional variations on the scale of the relative zero-point fluctuation substantially modifying the blockade performance. Beyond few-photon nonlinear optics, the same interaction protects an initially prepared molecular Bell state by dispersively decoupling molecular excitations from the lossy cavity mode, thereby suppressing cavity-mediated dissipation and slowing the fidelity decay. Our results establish microwave-shielded interactions as a unified interaction resource for engineering few-photon nonlinearities and protecting quantum states in molecular cavity-QED systems.

Characterizing intermediate quantum chaos in systems with mixed phase spaces requires probing spectral correlations across multiple scales. Here, we present a multi-scale statistical analysis of the resonance spectra of ultracold 166Er and 168Er atoms using higher-order spacing ratios, spacing increments, and the accumulated power spectral density of decimated level sequences. Our analysis reveals that while the spectra exhibit clear level repulsion, they are consistently better described by semi-Poisson than by standard Wigner-Dyson ensembles, with this intermediate character persisting across higher-order correlations. The accumulated power spectrum reveals an intermediate 1/f^alpha regime, with alpha approximately 1.5 at full spectral resolution and a tendency toward the integrable random-walk limit alpha = 2 at intermediate decimation scales, indicating a scale-dependent loss of spectral rigidity. Furthermore, the spacing-increment distributions independently reveal statistics intermediate between Laplace and Gaussian limits. Together, these complementary diagnostics show that the semi-Poisson character reflects persistent non-Wigner-Dyson correlations beyond the nearest-neighbor scale, providing a consistent picture of intermediate quantum chaos and a possible connection to hierarchical trapping in the underlying collision dynamics.

Hydrodynamic relaxation in open quantum many-body systems can be understood in terms of diffusive Nambu–Goldstone (NG) modes associated with strong-to-weak spontaneous symmetry breaking (SWSSB). Here, we find a common structure underlying weak ergodicity breaking in closed quantum systems and diffusive hydrodynamics in open quantum systems. We establish a spectral correspondence between the parent Hamiltonian of rainbow quantum many-body scar (RQMBS) states in closed systems and the effective Lindbladian of strongly dissipative open systems with strong $U(1)$ symmetry and local charge-dephasing Lindblad operators. For broad classes of spin, fermionic, and bosonic systems, we construct asymptotic RQMBS (ARQMBS) states. Within the enlarged scar subspace, the parent Hamiltonian admits a local Rokhsar–Kivelson-type frustration-free decomposition. A single-mode variational construction gives an $O(k^2)$ upper bound on its lowest excitation energy and yields states satisfying the defining ARQMBS criteria: orthogonality, vanishing energy variance, and characteristic entanglement scaling. Under vectorization, RQMBS states are mapped to maximally mixed infinite-temperature states in fixed-charge sectors, which exhibit long-range SWSSB order under finite-density conditions. Assuming uniformly bounded positive transition rates, the effective Lindbladian has the same local projectors as the parent Hamiltonian but with positive transition-dependent weights, yielding two-sided spectral bounds and the same system-size scaling of their gaps. For uniform rates, the two operators are proportional to each other and their eigenmodes coincide. Therefore, an exact quadratic parent-Hamiltonian branch maps directly to a diffusive NG mode. Our results connect different mechanisms of slow relaxation in closed and open quantum many-body systems.

Impurity nonequilibrium dynamics in quantum many-body systems is a frontier subject in ultracold-atom physics, which helps uncover microscopic mechanisms of polaron and collective excitation phenomena. Although repulsive interacting systems have been extensively investigated, attractive-interaction dynamics lacks systematic studies owing to intricate couplings among different states. We study a one-dimensional ideal Fermi gas containing an attractive spin-down impurity with initial momentum and analyze time evolution of interspin two-body correlations and impurity momentum. Using exact Bethe ansatz solutions, we simplify correlation matrix elements to finite sums for efficient calculations of eigenstate occupations and long-time dynamical evolution. In the weakly attractive regime, bound-state features appear when the total momentum is below or equal to the Fermi momentum, while mixed oscillatory behaviors arise at larger total momentum. The strongly attractive regime shows locally bound-state dominated dynamics with localized correlation peaks and scattering-induced Friedel-like oscillations. We characterize quantum flutter, the periodic oscillation of impurity momentum, and obtain consistent critical chemical potential values of the impurity via four independent Bethe ansatz approaches. This work clarifies quench dynamics and quantum flutter under attractive interactions, improves the understanding of nonequilibrium quantum many-body properties, and provides theoretical support for relevant ultracold-atom impurity experiments.

While a Fabry-Perot (FP) cavity interferometer is a standard tool in modern optics, analogous physics can also be observed with matter waves such as Bose-Einstein condensates (BECs). Interatomic collisions in BECs can induce nonlinear features such as multi-valued transmission spectra for interferometry experiments. Observation of interaction effects in the matter-wave analogue of FP interferometry has not yet been experimentally achieved, mostly due to broadening of transmission peaks by interaction-induced shifts of cavity resonances. In this work, we theoretically show that this limitation can be overcome by engineering a momentum chirp in a BEC wavepacket. By balancing this chirp against the dynamical mean-field energy shift of the cavity during transmission, we predict enhanced resonant spectra and bistability features due to interaction-assisted tunneling that would be otherwise inaccessible. We show that this physics is readily observable with quasi-one-dimensional BECs of cold atoms colliding with a pair of light-induced Gaussian potential barriers forming the cavity. We also benchmark the requisite optimum chirp needed and characterize the contrast of near-discontinuous FP transmission spectra at an interacting resonance. Our findings enable the realization of technologies like matter-wave switching protocols, and high-precision velocity measurements; estimates of achievable sensitivity for the latter are provided.

Beyond the frozen-pairing contribution, the superfluid weight of a multiband superconductor contains a correction from the self-consistent relaxation of the pairing amplitudes under a phase twist. We show that this correction has a quantum-geometric origin, not in the space of Bloch states but in the space of collective pairing fields. For an attractive multiband Hubbard model at the mean-field level, the thermodynamic Hessian governing the pairing response coincides with the static Gaussian pair-fluctuation kernel at zero momentum at every temperature below $T_c$. Since gauge invariance relates a phase twist to a finite center-of-mass momentum of the pair field, the superfluid weight is set by the small-momentum curvature of the eigenvalue branch that evolves from the Goldstone mode, and the relaxation correction is a stiffness-weighted quantum metric of the corresponding soft eigenvector on the manifold of pairing configurations. The same construction yields the Ginzburg-Landau gradient coefficient with the identical geometric correction, hence the coherence-length tensor, and, with the low-frequency dynamics of the kernel, the pair-mass tensor. Lieb-lattice calculations verify the equivalence of the thermodynamic, fluctuation-kernel, and soft-mode formulations and the independence of the superfluid weight from the orbital embedding.

How a mobile impurity exchanges momentum and energy with a many-body environment is a central question in nonequilibrium quantum physics. This exchange becomes particularly complex when the impurity moves rapidly and is strongly coupled to its bath. Here, we realize an interaction-tunable cold-atom collider with fermionic $^{40}$K impurities immersed in a $^{23}$Na Bose–Einstein condensate. A species-selective Raman pulse simultaneously launches the condensate and quenches the interspecies scattering length, producing an initial relative speed 36 times the condensate speed of sound. We track the ensuing relative motion as the interspecies interaction is tuned from weak coupling to resonance. At weak and intermediate coupling, the impurity dynamics are well described by a finite-energy two-body collision model. Near resonance, however, the early-time impurity acceleration is markedly suppressed relative to the two-body prediction, corresponding to a strong enhancement of the apparent dynamical inertia. The enhancement is observed near two distinct Feshbach resonances. These observations provide evidence for a polaronic response in the strongly coupled, supersonic regime of a degenerate Bose–Fermi mixture.

We show that Hawking radiation spontaneously generated by dynamical black-hole formation can be realized in a local quantum spin model. Starting from the ground state corresponding to the early-time vacuum, with no initially prepared excitations, we follow the dynamics induced by the evolving geometry. In the corresponding continuum quantum field theory, Hawking particle creation is associated with nontrivial Bogoliubov mixing between the in- and out-mode decompositions. We derive the resulting renormalized local energy density in the continuum quantum field theory, which corresponds to an energy-weighted integral of the Hawking spectrum, and identify its counterpart as a local spin observable whose expectation value quantitatively reproduces the continuum prediction, including its dependence on the surface gravity. The agreement persists over a finite time window whose duration grows logarithmically with system size, revealing the crossover from universal continuum Hawking physics to lattice-cutoff dynamics. These results establish a microscopic route to quantum simulation of quantum field theory in curved spacetime, connecting field-theoretic particle creation directly to experimentally accessible many-body observables.

The Jx-Jy Blume-Capel model, with and without an external magnetic field, is investigated using the Monte Carlo method. Our aim is to verify the location of the tricritical point, examine the existence of reentrant behavior, and study the magnetocaloric properties of this specifically designed system. In the non-field case, the critical temperature given by the Monte Carlo simulations is lower compared to the conventional antiferromagnetic square lattice. When an external magnetic field is applied, we observed signs of a one-dimensional spontaneous symmetry-breaking phase in the mixed-states regime, which is counterintuitive with secondary specific-heat peak. We take the snapshot of this state and elucidated the underlying mechanisms using mean field theory in the field of spin-1 spinor Bose-Einstein condensate . Regarding reentrant behavior, we confirmed its occurrence as the magnetic field varies, but only when the temperature is kept constant. Temperature-dependent reentrant behavior could not be identified within this system. Finally, we examined the magnetocaloric properties and found no evidence of first-order phase transitions within the region of the considered phase-transition point. Although both inverse and direct magnetocaloric effects are present and related to the mixed state, we identified that the relevant phase transitions are all of the second order.

We propose a four-level scheme for generating a uniform imaginary magnetic field for ultracold atoms using Raman coupling and state-dependent particle losses. Eliminating two lossy auxiliary states and subsequently projecting onto one of the two remaining Raman-dressed atomic states, we find that the atomic motion is affected by an imaginary vector potential which is engineered to have a constant curl corresponding to a spatially uniform imaginary synthetic magnetic field. We analytically derive the Gaussian wavepacket dynamics for atoms in such a synthetic magnetic field and study the width-dependent transport, unlike motion in a uniform real magnetic field. We find numerical agreement between the two-level non-Hermitian Hamiltonian and the adiabatic projection description, confirming the applicability of the latter adiabatic approach.

We study higher moments of quantum fluctuations and microscopic correlation energy in the three-dimensional Gross–Pitaevskii regime. For a smooth, compactly supported, nonnegative radial interaction and an $N$-independent normalized initial condensate wave function in $H^6$, we propagate the third spectral moment of a positive correlated energy along the exact many-body evolution. A uniform initial bound controls the third moment of the depletion number and sixth moments of collective one-particle fluctuations on every compact time interval, independently of $N$. For a nonzero interaction and a smooth, compactly supported initial condensate, we also identify the excitation-energy measure at scale $N^2$, weighted by the first energy. This limit holds for norm-convergent pure-state families with a uniform bound on the initial third correlated-energy moment. Its positive-energy part is a spectral measure of the relative scattering operator $-2\Delta+v$, multiplied by a time-dependent dressed contact. We compute the radial density and prove strict positivity of the contact for nonzero references with finite kinetic energy and finite mean excitation number. The same evolving states have bounded physical fluctuation moments and excitation-energy moments of order $p>1$ bounded below by a positive multiple of $N^{2p-2}$, in the extended spectral sense. Inverse powers of the positive initial energy construct admissible pure and mixed families. For the norm-convergent pure-state class, we determine the scaled interaction-energy moments for $1\le p<2$ and bound the second moment by $O(N^2)$. A finite cubic approximation has norm error $O(N^{-3/16})$ on the energy form domain and $O(N^{-5/4})$ on a stronger initial domain of its generator.

We show that trap rotation overrides Kibble-Zurek (KZ) universality when a 2D Bose-Einstein condensate forms via a finite-time quench. In stochastic Gross-Pitaevskii simulations, the rotation-induced chiral bias drives the vortex number close to the Feynman-Onsager (FO) value, independent of the quench time. Ramping up rotation with the quench restores the KZ scaling at low frequencies, with vortex numbers below the FO value. The vortices form a disordered lattice that a pinning potential can order into an Abrikosov lattice, even at low rotation. Trap rotation thus offers a control knob for critical dynamics.

We show that quasi-paraxial fluids of light propagating in nonlinear Kerr media naturally experience an effective gauge field arising from the intrinsic spin-orbit interaction of light encoded in Maxwell equations. By deriving the leading nonparaxial corrections to the nonlinear Schrödinger equation governing the fluid evolution, we find that, for circularly polarized beams, these corrections generate an effective Lorentz force that modifies the fluid dynamics. Remarkably, the corresponding gauge field is not externally imposed but self-consistently generated by the fluid itself, with its sign controlled by the light helicity and its magnitude governed by the spatial variations of the nonlinear refractive index. We illustrate the physical consequences of this gauge field in two configurations: the chiral deformation of normal and superfluid flows past an obstacle, accompanied by a transverse spin Hall force, and the cumulative precession of an optical wave packet trapped in a focusing fluid.

We study $N$ bosons on a three-site ring with nearest-neighbour hopping, density-assisted hopping of the opposite sign, and an on-site interaction, written equivalently as three momentum modes $q=0,\pm1$ with a channel-dependent quartic vertex. As the density-assisted term grows it inverts the band, and exact diagonalisation up to $N=400$ shows a first-order ground-state transition from the $q=0$ condensate to a state with $\avg{n_{+1}}=\avg{n_{-1}}\simeq N/2$; the avoided-crossing gap closes faster than exponentially with $N$. In its symmetric eigenstates this state is fragmented in the Penrose-Onsager sense, and at leading order in $N$ it has the energy of a Gross-Pitaevskii (GP) standing wave, as expected for a few-mode mean-field system. We distinguish two cases by switching off umklapp scattering, which is available because $2\cdot(+1)\equiv-1$ on three sites, and by varying the on-site interaction. Without umklapp, $n_{+1}-n_{-1}$ is conserved, the eigenstates form a rotor tower whose lowest member is a fragmented Fock state, and the transition leads to fragmentation exactly when the $\pm1$ interaction is attractive, as predicted by Nozières' exchange argument with the sign of the interaction reversed; otherwise the system forms a single current-carrying condensate. With umklapp, the standing wave is pinned to three positions, the three momentum sectors become degenerate exponentially in $N$, and coherent admixture of the $q=0$ mode stabilises the standing wave over a wider range; the exact ground state is then a number-squeezed, momentum-projected standing-wave condensate rather than a robustly fragmented state.

We show that a collisionless point-vortex gas behaves as a genuine Vlasov system. Starting from the non-canonical Hamiltonian structure of point-vortex dynamics, we derive a reduced Klimontovich equation and the corresponding Vlasov description directly on the physical plane, without introducing auxiliary momentum variables. The theory predicts a rich spectrum of collective phenomena, including vortex Rossby-like waves, global oscillations, resonant critical layers, and a dielectric response operator analogous to that of plasmas and self-gravitating systems. For smooth vortex distributions, differential rotation gives rise to a vortex analogue of Landau damping, whereby collective modes transfer energy and angular momentum to resonant vortex annuli in the absence of collisions or dissipation. Remarkably, differential rotation acts as a stringent mode-selection mechanism: for a Gaussian vortex cloud, only a single quadrupolar quasi-mode survives, while all higher multipoles are destroyed by phase mixing. We further predict a vortex counterpart of the beam–plasma instability, in which a thin vortex ring orbiting an Onsager patch selectively excites Kelvin modes with an azimuthal symmetry controlled by the ring population. These results establish kinetic theory as a natural framework for collective vortex dynamics and suggest new mechanisms for vortex transport, vortex avalanches, and rotational glitches in neutron-star superfluids.

Harnessing integrability at arbitrary interaction strength, we build long-lived and spatially localized soliton-like excitations in the 1D Bose gas with contact repulsion (Lieb-Liniger model). We show how these density profiles, built from selective superpositions of multi-hole eigenstates, display long lifetimes tunable through parameters like dispersion width and number of holes. We further illustrate a scaling protocol that in principle enables unbounded increase in lifetimes in the thermodynamic limit. Our results showcase classical-like excitations originating from unitary quantum dynamics, away from the semi-classical regime, and in regions where interactions cannot be treated perturbatively.

Spin-1 Bose-Einstein condensates provide a fundamental platform for exploring quantum magnetism and coherent spin dynamics. Their three Zeeman components are coupled by spin-exchange interactions under conserved particle number and magnetization, giving rise to rich magnetic phases and nonlinear spin structures; a rigorous classification of stationary states therefore clarifies the possible spin configurations and their physical origins. In this paper, as a continuation of our previous work \cite{LLWZ23,LLWZ24,LLWZ26}, we study the classification of solution for the spin-1 Bose-Einstein condensates in $\mathbb{R}$, by applying the non-standard Hirota's bilinear method, the determinant formula and the general left-uniqueness strategy, we present some novel explicit solution and give the complete classification of the solutions. In comparison with the classical Hirota's method, the non-standard Hirota's method needs to introduce auxiliary functions to construct Hirota's bilinear form, which essentially adds additional constraint on the wave functions, so that it possible to deal with more complex equations. As far as we know, this is the first systematic study the complete classification of solutions to the spin-1 Bose-Einstein condensates. Our classification results yields experimentally testable predictions for spin populations, relative phases, and magnetic textures

We study the effect of repulsive interactions on skin modes and entanglement in the unidirectional Bose–Hubbard model. Although the ground state of this integrable model undergoes a superfluid–Mott transition with suppressed skin accumulation, the corresponding change in entanglement depends on the selected many-body state. We compare the ground state with the steady state, defined as the largest-imaginary-energy eigenstate in the zero-momentum sector. For sixteen sites at unit filling, we find that the ground-state density approaches uniform filling and its entanglement decreases toward the Mott limit, whereas the steady states exhibit a nonmonotonic entanglement crossover. At intermediate interaction ($U/t=10$), the second Rényi entropy of the steady state develops four spatial maxima. We show that these maxima are accounted for by competing subsystem-number sectors within a dominant configuration manifold of four doublons and four holons, but not by the subsystem particle-number variance alone. The dominant global onsite-pair manifold retains the same entropy peak structure. At strong interaction, projection onto the one-doublon–one-holon manifold closely reproduces the full entropy profile and its two edge peaks. Furthermore, we demonstrate that the many-body winding about the origin remains nonzero at strong interaction; an exact two-site example shows how it becomes ill-defined in the singular projected $U/t\to\infty$ limit. Our results distinguish ground-state skin suppression from state-selective entanglement and point-gap topology.

We present release 3.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open-source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. As in previous releases, development is centered on common data formats, on libraries to simplify and speed up code development, and on full-featured simulation programs that let non-experts carry out serial or parallel numerical simulations using the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended-ensemble simulations, exact and full diagonalization (ED), the density matrix renormalization group (DMRG), and continuous-time QMC solvers for dynamical mean-field theory (DMFT). Major changes in release 3.0 include distribution of the pyalps binary through the Python Package Index (pip install pyalps) and through Spack for HPC systems; migration of development to GitHub with continuous integration and automated testing; relicensing of the package under the permissive MIT license; a rebuilt documentation and tutorial website, including a set of Jupyter-notebook tutorials and localized content; a broad modernization of the C++ codebase (C++17 compliance, Boost and NumPy 2.0 compatibility, and warning and dead-code cleanup) together with a major DMRG update and associated reliability and build-compatibility fixes; the removal of legacy components (the VisTrails provenance integration and the TEBD, MPS, and directed-worm-algorithm application codes); archival of release 3.0.0 with a Zenodo DOI; and a formal governance and sustainability model developed under the US National Science Foundation (NSF) POSE program. The software is available at https://alps.comp-phys.org.

Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconductivity, Cooper pairing at finite-momentum driven by spin imbalance, remains difficult to study numerically in two dimensions. The fermion sign problem in the spin-imbalanced attractive Hubbard model severely limits quantum Monte Carlo access to the low temperatures where pairing correlations are most pronounced. We show that the Holstein model, in which the effective attraction is retarded by phonon dynamics, strongly suppresses sign fluctuations under spin imbalance, enabling unbiased determinant quantum Monte Carlo simulations deep in this regime. With increasing magnetization, zero-momentum pairing evolves into finite-momentum FFLO pairing, evidenced by oscillating real-space pair correlations and a shift of the pairing susceptibility peak to finite $\vec{Q}$. The finite-momentum susceptibility substantially exceeds its zero-momentum value. At high density $n \sim 0.875$, near magnetization $m \approx 1-n$ where the majority-spin density approaches half filling, FFLO pairing coexists with checkerboard charge-density-wave order at short range, while staggered longitudinal spin correlations emerge as a magnetic halo at larger distances. These spatial textures provide experimentally accessible signatures of finite-momentum pairing and motivate cold-atom realizations of FFLO physics.

At the low densities encountered in the inner crust of neutron stars, neutron matter becomes superfluid, directly influencing macroscopic phenomena such as star cooling. However, a precise quantitative understanding of this system has been hindered by uncontrolled many-body approximations. In this letter, we exploit the nonperturbative Gorkov algebraic diagrammatic construction method to provide a unified ab initio description of the equation of state and the single-particle spectral functions of homogeneous neutron matter, spanning from the dilute regime to saturation density. We identify clear microscopic signatures of superfluidity in the excitation spectrum and quantify the pairing gaps associated with realistic chiral Hamiltonians, estimating for the first time their sensitivity to the many-body truncation. Our robust handling of both dynamical and pairing correlations paves the way for exploring electroweak processes across the regimes relevant to neutron stars.

We present a numerical study of a blockade Hamiltonian with a local $\mathbb{Z}_2$ symmetry. For such models, the existence of a topologically ordered ground state has been proven rigorously for weak Rabi driving [T. F. Maier et al., PRX Quantum 6, 030340 (2025)]. However, the phase diagram beyond weak driving is still an open question. We perform large-scale infinite-size density matrix renormalization group (iDMRG) simulations on a cylinder and map out the full phase diagram of the model. We find a single continuous phase transition between the topologically ordered phase and a trivial paramagnetic phase, and confirm numerically the stability of the topologically ordered phase up to Rabi drives comparable to the characteristic energy scale of the system. Furthermore, we compute the excitation gaps and show that the charge gap remains open across the phase transition and charges become confined in the trivial phase. On the other hand, the gap in the symmetric sector closes at the transition, consistent with the condensation of flux excitations.

We investigate whether the geometry of a curved conductor–dielectric interface can modify the effective mass and propagation of scalar surface plasmons (SSPs). In the hydrodynamic limit, we describe SSPs as density perturbations in a semiclassical jellium electron gas. Using a thin–layer dimensional reduction procedure, we derive a massive Klein–Gordon equation for SSPs on the curved interface with a geometric potential. We find that under specific conditions, the geometric potential allows SSP propagation with a vanishing effective mass. We then apply the formalism to a well known analogue black hole configuration, determining the corresponding quasinormal mode frequencies. Our results show that interface curvature provides a geometric handle on plasmon propagation and localization, suggesting new possibilities for plasmon transport and analogue gravity experiments with surface plasmons.

We demonstrate a general mechanism to generate scalable, metrologically useful entanglement in the form of generalized spin squeezing in multi-level (qudit) spin systems with power-law interactions. We demonstrate this for spin-changing dynamics, extending the physics of spin-nematic squeezing known from spinor BECs to spatially extended lattice geometries. In the absence of all-to-all connectivity, the dynamics generally suffers from leakage into non-collective finite-momentum modes, destroying Heisenberg scaling. We show that introducing local $\mathfrak{su}(D)$ exchange interactions establishes an energy gap that protects the fully symmetric collective manifold, suppressing non-collective excitations, and achieving Heisenberg scaling $\xi^2\propto1/N$ in one to three dimensions for any power-law exponent. We provide an analytical understanding of the mechanism via the analysis of the Bogoliubov excitation spectrum, and demonstrate its effectiveness via numerical simulations of the quantum many-body dynamics. Our findings more broadly apply to other spin-squeezing Hamiltonians, and directly generalize to larger internal spin dimensions. Our results provide a practical pathway towards scalable Heisenberg-limited qudit entanglement generation in experimental platforms such as arrays of (polar) molecules, magnetic atoms, and Rydberg atoms realizing multi-level spin systems.

Counterdiabatic driving provides a universal route to fast quantum control, but its standard formulation presupposes a Hamiltonian generator. We develop a Hamiltonian-free theory of the adiabatic gauge potential (AGP) for parameter-dependent unitary operators $U_\lambda$. Exact and regularized versions of the AGP follow from a Hermitian, banded linear system at a cost comparable to the Lanczos construction in the Hermitian case, requiring neither diagonalization nor a matrix logarithm. We demonstrate rapid convergence of the truncated Krylov approximation for the kicked top and for an integrable brickwork circuit built from unitary $\check{R}$ matrices, a setting in which no local effective Hamiltonian is available. Our results extend shortcuts to adiabaticity, and AGP-based diagnostics of chaos, to Floquet systems and quantum circuits.

Measuring part of a many-body wave function generates a projected ensemble of pure quantum states on the unmeasured subsystem. Recent advances in deep thermalization have shown that, in chaotic systems, this ensemble universally converges to the maximally random ensemble compatible with its average density matrix, known as the Scrooge ensemble. By contrast, free-fermion systems are nonchaotic, and their quantum states are constrained to remain Gaussian. Motivated by this distinction, we introduce the fermionic Gaussian Scrooge ensemble to describe deep thermalization in generic free-fermion systems. This ensemble is defined as a distortion of the Gaussian Haar ensemble by the average density matrix, and its moments reveal an enlarged symmetry of Gaussian states in the replicated Hilbert space. We demonstrate the emergence of the fermionic Gaussian Scrooge ensemble in two complementary settings: (1) we analytically prove that the projected ensemble of the SYK$_2$ model follows the fermionic Gaussian Scrooge ensemble at arbitrary evolution times; (2) we provide numerical evidence that it emerges at sufficiently long times in random Gaussian circuits with charge conservation. Our results establish the fermionic Gaussian Scrooge ensemble as a universal description of deep thermalization in free-fermion systems.