Abstract
Geometrical frustration and long-range couplings are key contributors to create quantum phases with different properties throughout physics. We propose a scheme where both ingredients naturally emerge in a Raman induced subwavelength lattice. We first demonstrate that Raman-coupled multicomponent quantum gases can realize a highly versatile frustrated Hubbard Hamiltonian with long-range interactions. The deeply subwavelength lattice period leads to strong long-range interparticle repulsion with tunable range and decay. We numerically demonstrate that the combination of frustration and long-range couplings generates many-body phases of bosons, including a range of density-wave and superfluid phases with broken translational and time reversal symmetries, respectively. Our results thus represent a powerful approach for efficiently combining long-range interactions and frustration in quantum simulations.
Similar content being viewed by others
Introduction
In a diverse range of systems from neutron stars1 to nuclei2,3 and electrons4,5, intrinsic strong long-range interactions are essential for stabilizing strongly correlated states6. Large scale quantum correlations are an essential element in a wide variety of phenomena7. These couplings between far-spaced elementary constituents can lead to interesting properties such as static entanglement8,9 and symmetry breaking10,11,12, as well as efficient spreading of highly non-local correlations in out-of-equilibrium configurations13,14. In condensed matter systems this situation is further enriched by geometric frustration which can compete with the electrons’ native screened Coulomb interaction15. Under these conditions, a large ground state degeneracy can occur, topological16,17,18 and spontaneously-symmetry-broken phases19,20 can take place. Here we describe a 1D lattice for ultracold atoms with effective geometric frustration, and interactions extending over several lattice sites.
Our analysis is particularly relevant because the simultaneous presence of long-range couplings, geometrical frustration and quantum fluctuations challenge all current theoretical treatments21,22; it is, therefore, crucial to back theoretical predictions with accurate experiments. While very recent solid state experiments tackle configurations where long-range couplings coexist with kinetic frustration23, the injection of geometrical frustration remains a distant goal. In this respect, promising initial results in specific geometries of tweezer-arrays of Rydberg atoms have been obtained24,25, however, complimentary experimental realizations for itinerant systems such as neutral atoms in optical lattices26,27 are lacking. Without frustration, the role of long-range interactions have been explored for magnetic-atom28,29, polar molecules30, and cavity QED31 systems. Without long-range interactions, geometrical frustration has also been experimentally investigated only in weakly interacting regime32,33,34,35,36,37, while strong interactions have never been explored. Even theoretical proposals to engineer geometrically frustrated strongly correlated phases mainly concentrate on systems with contact interactions38,39,40,41,42,43,44 or, very recently, nearest-neighbor repulsion45. This work provides a significant step forward by describing quantum gases in strongly interacting regimes where quantum fluctuations, geometrical frustration and long-range couplings strongly compete.
We consider the many-body physics of a recently realized class of subwavelength 1D optical lattices46,47,48,49,50,51,52,53,54,55 and show that they are a suitable platform for combining geometric frustration and finite-ranged interactions. These lattices use Raman transitions to couple N internal atomic states with lasers of wavelength λ slightly detuned from the Raman resonance condition [Fig. 1a]. This setup is described by an extended Hubbard Hamiltonian (Bose or Fermi) where the lattice period is reduced from λ/2 to λ/(2N). In contrast to existing optical lattice systems—such as magnetic atoms56,57, weakly dressed Rydberg atoms58, and polar molecules59,60—where the spatial decay of the interaction is fixed and tunneling processes occur between neighboring lattice sites, our lattice allows for interactions and tunneling with a tunable range. In particular, we show that: (1) the range and sign of tunneling processes can be controlled giving rise to effective geometric frustration [see Fig. 1b]; and (2) the interactions can be approximated by a power law whose exponent is a function of the Raman coupling strength.
a Lasers induce two-photon Raman transitions of intensity Ω that cyclically couple N = 3 consecutive internal states (labeled by m) with energy difference ϵm. b Subwavelength lattice with long-range tunneling; all links emanating from the j = 0 site (with dressed state index n = 0 and unit cell ℓ = 0) have their tunneling strength labeled. Top: synthetic dimension picture with triangular plaquettes and the potential for geometric frustration, with representative tunneling strengths labeled. Bottom: corresponding 1D lattice with explicit long-range links.
We then turn to a specific implementation based on bosonic 87Rb and identify a range of strongly correlated regimes through a matrix-product-states (MPS) analysis61. When the N = 3 states of the F = 1 hyperfine ground state are considered, we find a normal superfluid (SF) phase in the regime of weak interaction and short range tunneling. Detuning from Raman resonance introduces geometrical frustration leading to a chiral superfluid (CSF) phase with broken time-reversal (TR) symmetry; similar CSFs have been predicted in far ranging systems from cold atoms in p-orbitals62,63,64 to hadrons65. Extending the interaction range destabilizes the SF phases in favor of a spontaneous symmetry broken (SSB) density wave (DW1/2) insulator consisting of alternating occupied and empty sites. These phases persist for N = 5 internal states (modeling the F = 2 hyperfine manifold of 87Rb). The effective lattice period decreases as N increases, making long-range repulsion more significant. At filling factor 1/3 this stabilizes a period-3 density wave (DW1/3) never achieved in cold atoms setups with individual bosons always separated by two empty lattice sites. Finally, we provide a detailed protocol for state preparation and detection, providing complete experimental access to all the interesting many-body regimes. Our results provide a solid and alternative route to explore geometrically frustrated quantum matter in presence of strong long-range correlations.
Results
Physical setup
We consider the one dimensional sample of ultracold bosons of mass ma illuminated by a pair of counterpropagating lasers of wavelength λ. This geometry serves to define the two photon recoil wavenumber kR = 2π/λ and energy \({E}_{{{{\rm{R}}}}}={\hslash }^{2}{k}_{{{{\rm{R}}}}}^{2}/(2{m}_{a})\). The optical fields induce two photon Raman transitions that cyclically couple N internal atomic states (labeled by the index m = 0, 1, …, N − 1) with strength Ω. The cyclic coupling condition is fulfilled by Raman-coupling the m = 0 and the m = N − 1 states [show for N = 3 in Fig. 1a]. We explore the configuration of nearly-resonant Raman coupling, where each transition is detuned by a small amount δm = ϵm−1 − ℏδωm−1 from resonance; as shown in Fig. 1a, ϵm is the energy difference between consecutive states m and m + 1, and δωm is the corresponding Raman frequency-difference.
This scheme can be described by the light-matter Hamiltonian density \({\hat{H}}_{{{{\rm{LM}}}}}(x)={\hat{H}}_{{{{\rm{R}}}}}(x)+{\hat{H}}_{{{{\rm{d}}}}}(x)\), with contributions from Raman coupling
and detunings
both expressed in terms of bosonic field operators \({\hat{\phi }}_{m}^{{{\dagger}} }(x)\) and \({\hat{\phi }}_{m}(x)\). These describe the creation and annihilation of a particle in internal state m = 0, 1, . . . , N − 1 at position x. Owing to the cyclic coupling, we label the internal states periodically so that \({\hat{\phi }}_{m}^{{{\dagger}} }(x)={\hat{\phi }}_{m+N}^{{{\dagger}} }(x)\), and we adopt an energy zero such that the detunings sum to zero, \({\sum}_{m}{\delta }_{m}=0\). The operator \({\hat{{{{\mathcal{H}}}}}}_{{{{\rm{R}}}}}\) effects a tight-binding lattice in a synthetic dimensional space where each internal state m corresponds to a synthetic lattice site66. In this synthetic dimension picture, the Raman coupling in \({\hat{{{{\mathcal{H}}}}}}_{{{{\rm{R}}}}}\) includes a Peierls phase 2kRx on each hopping term, while the detuning term \({\hat{{{{\mathcal{H}}}}}}_{{{{\rm{d}}}}}\) captures on-site energies. In analogy with conventional tight-binding lattices in real space, we rewrite \({\hat{{{{\mathcal{H}}}}}}_{{{{\rm{LM}}}}}\) in a dressed state representation using the synthetic-dimension momentum states basis
with n ∈ {0, 1, ⋯ , N − 1}. This transformation diagonalizes the Raman coupling operator
with energies
describing the usual cosinusoidal tight-binding dispersion with minima shifted from zero “crystal momentum” by the Peierls phase 2kRx. In terms of the real-space coordinate x, \({\hat{H}}_{{{{\rm{R}}}}}(x)\) defines a set of N cosinusoidal adiabatic potentials with period λ/2. The potentials corresponding to neighboring dressed states are separated from each other by a subwavelength spacing a = λ/(2N), as illustrated in Fig. 2. The potential minima are located at spatial positions xj = aj given by the subwavelength lattice site index j = n + Nℓ, itself defined by both the unit cell ℓ of the underlying λ/2 lattice as well as the dressed state n.
In the synthetic-dimension momentum representation the detuning Hamiltonian density
has off-diagonal terms that induce long-range tunneling. For odd N this can be expressed in a conventional tunneling form
with matrix elements
given by a discrete Fourier transform of the detunings. For even N the Δn sum runs from 1 to N/2; for Δn = N/2 the tunneling matrix element must be divided by 2 to avoid double counting. In any case, these complex valued tunneling matrix elements can be expressed as the product \({\gamma }_{\Delta n}=| {\gamma }_{\Delta n}| \exp (i{\phi }_{\Delta n})\) of a strength \(| {\tilde{\gamma }}_{\Delta n}|\) with γ0 = 0 (since ∑mδm = 0) and a Peierls phase ϕΔn with ϕΔn = −ϕ−Δn (since δm is real valued). We focus on symmetric patterns of detuning (i.e., δm = δ−m), in which case γΔn is additionally real-valued, but still can be long-ranged with a combination of positive and negative contributions.
Band structure and tight binding description
The preceding discussion concluded with a continuum description of our sub-wavelength lattice; to describe the many-body physics of this configuration we now construct the 1D lattice model for atoms in the lowest Bloch band. Without interactions this provides an exact description of the low-energy physics, and for large enough Raman coupling interaction-mixing of higher bands can be neglected.
In the absence of detuning, the lattice consists of N independent sinusoidal potentials each with the ground-band Wannier states shown in Fig. 2. In general the operator
describes the creation of an atom in the r-th Bloch band of the n-th sublattice with Wannier amplitudes wr(x) computed for a sinusoidal-lattice67, and as above j = n + Nℓ. In what follows we focus on the lowest (r = 0) band and succinctly label these operators via \({\hat{b}}_{j}^{{{\dagger}} }\), using the subwavelength lattice index j alone.
The native tunneling strength within these independent lattices
couples states separated by ∣Δj∣ = N sublattice sites, depends only on the Raman coupling strength Ω, and is defined to be zero for ∣Δj∣ ≠ N.
Additional couplings, which are significant for distances Δj < N, are provided by detuning induced tunneling
that is proportional to the overlap integral between Wannier functions associated with different atomic dressed states.
Figure 3 demonstrates that the combined tunneling \({J}_{\Delta j}={J}_{N}^{\Omega }+{J}_{\Delta j}^{\delta }\) has a variable sign and significant long-ranged contributions; markers are computed directly from Wannier functions and curves use the Gaussian approximation and γΔj (see “Methods” for details). Panel (a) illustrates the simple N = 3 case for three different values of γ1 (as we note in Methods, this suffices to fully quantify the detuning induced tunneling in this case), with fixed native tunneling. This case illustrates both the long-range character of JΔj as well as the sign-inversion between short and long range. Panel (b) turns to the case of N = 5 internal states—specified by both γ1 and γ2—enabling more complicated tunneling configurations, such as shown where γ1 and γ2 have opposite sign.
Computed directly from Wannier functions (markers, including native tunneling as marked), or from a Gaussian variational ansatz presented in Methods (curves, excluding native tunneling). a N = 3 internal state case, computed for Raman coupling Ω = 3.0ER and first detuning Fourier component γ1/ER ∈ {0.02, 0.1, 0.18}. b N = 5 internal state case, computed for Raman coupling Ω = 3.5ER, first detuning Fourier component γ1 = −0.06ER and second detuning Fourier component γ2/ER ∈ {0.03, 0.06, 0.09}.
Contrary to atoms in the ground-band of an optical lattice, where the tunneling strength—fixed by the shape of the Wannier function w(x)—is strictly positive and the nearest-neighbor contribution dominates, Eq. (11) along with Eq. (7) show that proper selection of detuning parameters δm allows for long-range tunneling with a combination of positive and negative contributions. This enables effective geometric frustration even in 1D, and in the following sections we explore the many-body interplay between effective geometric frustration and interaction processes.
Interacting processes
The preceding section defined the single-particle tunneling contribution to our system’s Hubbard model description. Here we continue by computing the two-body bosonic interactions with strength given by the overlap integral of atomic densities
where the pre-factor g1D describes the strength of the contact interaction that does not depend on the atomic internal state (a good approximation for ultracold 87Rb atoms in their ground electronic state68). Even for such simple underlying interactions, effective interactions between laser-dressed atoms generally include additional non-local contributions such as density assisted tunneling and pair tunneling (see refs. 69,70 for examples in the continuum). However, in the present case such terms are absent because the transformation in Eq. (3) is independent of the spatial coordinate and therefore leaves density-density interactions in Eq. (12) unchanged. As suggested by the Wannier orbitals in Fig. 2, the overall strength and range of VΔj can be tuned by modifying the effective lattice depth 2Ω which predominately affects the width of Wannier functions.
The detailed properties of VΔj are summarized in Fig. 4 (with numerical values suitable for 87Rb; see “Methods” for details), with markers denoting explicit numerical evaluation of Eq. (12) for N = 3 (red circles) and N = 5 (blue stars) internal states. The solid curves plot the result of a variational calculation using a Gaussian ansatz for the Wannier functions (see “Methods” for details). Panel (a) confirms that VΔj can be a significant fraction of V0 over the range of a few subwavelength lattice sites, and the inset shows that, as expected for the underlying sinusoidal lattice, the overall strength depends weakly on the Raman coupling strength Ω with a nominal V0 ~ Ω1/4 scaling.
Computed directly from Wannier functions (markers), or from a Gaussian variational ansatz (curves). In both cases, red and blue mark N = 3 and N = 5 internal states, respectively. a Dependence of interaction strength VΔj on distance d for Raman intensity Ω/ER = 2.5. (Inset) On-site interaction strength. b Top: Relative strength of long-range interactions V1/V0 as a function of Raman coupling Ω. Bottom: Effective power law exponent α versus Raman coupling Ω.
Comparison to other long-range interactions
For the usual case of ultracold atoms in the ground band of an optical lattice, the on-site interaction V0 greatly exceeds longer ranged contributions because Wannier functions are highly localized to individual lattice sites. As a result, additional contributions such as dipolar interactions are required to induce long-range interactions in cold-atom systems.
In conjunction with local interactions, long-ranged interactions can be modeled by the power-law interaction potential
In the dipolar case an applied electric or magnetic field induces interactions with α = 328; this limits the range of many-body phenomena that can be realized. Our scheme is not subject to this limitation and α is not fixed a priori.
For many-body physics, the very long-ranged tail of this interaction is often unimportant, making the interaction strengths at Δj = 0, 1 and 2 the only relevant contributions71. These can be quantified by the relative strength Vα/V0 of the power-law to local potentials, as well as the power law exponent \(\alpha ={\log }_{2}({V}_{1}/{V}_{2})\). The relative strength Vα/V0, shown in Fig. 4(b-top), confirms that for both N = 3 and N = 5 the long-range contribution can be significant; because the interactions ultimately derive from overlap integrals, we have 1 > Vα/V0 ≥ 0. Owing to the reduced spacing between subwavelength lattice sites for increasing N, Vα/V0 is larger for N = 5 than N = 3. Figure 4(b-bottom) shows that the power law exponent is not fixed (as it would be for dipolar or Van der Waals systems), but crucially it can be tuned simply by varying Ω; for our parameters α approximately resides in 5 ≲ α ≲ 8 for N = 3 and 2 ≲ α ≲ 3 for N = 5. This highlights the flexibility of our setup compared to conventional optical lattice realizations, and provides an avenue for realizing interesting many-body phases.
More broadly speaking, tunable power-law like scaling of the interaction’s range in two-level systems has been realized for spin-spin couplings in trapped ion systems72 and predicted for transversely confined hard-core dipolar bosons71. In contrast, our construction is applicable to itinerant gases of both bosonic and fermionic ultracold atoms, and, as we focus on below, it combines effective geometrical frustration in the single-particle degrees of freedom with power-law like scaling of the interactions.
Table 1 summarizes the interaction strengths V0,1,2 for the range of Raman coupling Ω that we focus on. As compared to the typical ≲kB × 5 nK = h × 100 Hz thermal energy scales for ultracold atoms in optical lattices, this shows that for N = 3 interactions are relevant for Δj = 0 and 1; for N = 5 the Δj = 2 interaction is also appreciable. The Δj = 1 nearest neighbor interaction strengths largely exceed those of magnetic lanthanide atoms56,57, and instead are comparable those predicted for polar molecules59.
Many-body phase diagram
The combination of terms derived in the previous section can be assembled into an extended 1D Bose-Hubbard (BH) Hamiltonian
where \({\hat{n}}_{j}={\hat{b}}_{j}^{{{\dagger}} }{b}_{j}\), off resonant couplings to higher bands can be neglected in the regime of large Raman coupling (Ω ≳ 3 for 87Rb, see “Methods” for details. For smaller Raman couplings, one would need to calculate renormalized Hamiltonian matrix elements, which arise due to higher bands73. While extended BH models including either geometric frustration or long ranged interactions have been widely studied74,75,76,77,78,79,80,81,82,83, our realization embodied by Eq. (14) is the first proposal to include both long-range interactions and effective geometrical frustration. In what follows, we use MPS calculations61 to obtain the resulting ground state phases for N = 3 and N = 5, and then we quantify the resulting phases using three quantities.
First, the staggered density
where \(\bar{n}={L}^{-1}{\sum}_{j}\langle {\hat{n}}_{j}\rangle\), signals period-2 density modulations, and serves as an indicator of spontaneously broken translational symmetry. In the thermodynamic limit a true period-2 SSB ground state would generally be an equally weighted superposition of these two symmetry broken configurations, making δNj = 0; in this case a higher order correlation function would be required to extract this order. When studying this order parameter we use an odd number of lattice sites which serves to explicitly break the degeneracy between the two configurations. This order parameter would be non-zero for either density-wave solids.
Second, the single particle Green function
quantifies the degree of spatial phase coherence; in 1D, an algebraic decay of this quantity at long-range (i.e., quasi long-range off-diagonal order) reveals the presence of gapless phases with SF properties.
Lastly we consider the correlation function
where \({\hat{\kappa }}_{j}=i({\hat{b}}_{j}{\hat{b}}_{j+1}^{{{\dagger}} }-{\hat{b}}_{j}^{{{\dagger}} }{\hat{b}}_{j+1})\) is the local current operator for the link between j and j + 1. The long-range order of \({\kappa }_{j}^{2}(\Delta j)\) indicates correlations between currents on links a distance Δj apart, and is associated with spontaneous breaking of TR symmetry. The same conclusions can be derived by calculating directly the order parameter \({\hat{\kappa }}_{j}\). Importantly, this strategy requires the addition of a weak term \({\hat{\kappa }}_{j}\) in Eq. (14) which allows breaking the ground state degeneracy associated to the currents going from left to right and vice versa with the same amplitude.
Our MPS calculations were performed on large systems with L ≈ 180 sites (figure captions provide exact details), and all reported quantities were evaluated in the Lcen = 100 central sites to minimize boundary effects (from open boundary conditions). In all cases truncation errors <10−7 were achieved by using bond-dimensions up to 1000. We check for quasi-long range order (LRO) by evaluating correlation functions at this maximum possible range with \({j}_{\min }=(L-{L}_{{{{\rm{cen}}}}})/2\) and Δj = Lcen, for example, one would quantify long range phase coherence and current correlations in terms of
N = 3 internal states
We begin our MPS analysis with N = 3 internal states and at a fixed particle density of \(\bar{n}=1/2\) atoms per subwavelength lattice site. Geometric frustration is induced by selecting the Fourier transformed detuning γ1 of Eq. (8) to be positive, which makes both \({J}_{1}^{\delta },{J}_{2}^{\delta } \, < \, 0\) while the bare tunneling \({J}_{3}^{\Omega }\) remains, as always, positive [see Fig. 3a]. In this case, the extended BH in Eq. (14) models a triangular ladder with both ferromagnetic and antiferromagnetic tunnel couplings [see Fig. 1b, top panel] and long-range interactions. Figure 5a shows that staggered density order (with δN > 0) is present for small γ1 and range of coupling strengths Ω (different curves). This indicates the presence of a SSB phase, but does not yet distinguish between supersolid and density wave (DW) insulating phases. Next, Fig. 5b shows that, by quantifying off-diagonal order, the one-body Green’s function \({g}_{{{{\rm{cen}}}}}^{(1)}\) delineates these cases. At small γ1 there is no phase coherence, implying that the system is a DW1/2 insulator. Changes to γ1 proportionately change the detuning induced tunneling amplitudes (here J1 and J2), but have no impact on the native tunneling (given by J3) nor the interaction strengths VΔj. Therefore increasing γ1 increases the kinetic contribution to the Hamiltonian, ultimately melting the DW insulator. Comparing (a) and (b) shows that \({g}_{{{{\rm{cen}}}}}^{(1)}\) becomes non-negligible concurrently with the vanishing of δN; as a result we conclude that no supersolid is present and the transition is from DW1/2 to conventional SF.
a Staggered density δN(1); b One-body Green’s function \({g}_{{{{\rm{cen}}}}}^{(1)}\); and c vector order parameter \({\kappa }_{{{{\rm{cen}}}}}^{2}\). \({g}_{{{{\rm{cen}}}}}^{(1)}\) and \({\kappa }_{{{{\rm{cen}}}}}^{2}\) are defined in Eq. (18). All three cases include Raman couplings Ω/ER = 3.0 (empty markers) and Ω/ER = 3.5 (filled markers). These were obtained in a L = 181 lattice site chain with 91 particles.
Lastly, Fig. 5c shows that the current-current correlation function \({\kappa }_{{{{\rm{cen}}}}}^{2}\) becomes non-zero for a range of γ1 and smaller Ω (empty markers). Comparing (b) and (c) shows that while this order parameter is only non-zero when \({g}_{{{{\rm{cen}}}}}^{(1)} \, > \, 0\), the reverse is not true. This allows us to disambiguate a conventional SF phase [\({g}_{{{{\rm{cen}}}}}^{(1)} \, > \, 0\) and \({\kappa }_{{{{\rm{cen}}}}}^{2}=0\)] from a TR broken CSF [\({g}_{{{{\rm{cen}}}}}^{(1)} \, > \, 0\) and \({\kappa }_{{{{\rm{cen}}}}}^{2} \, > \, 0\)].
In the next section, we show that by increasing the number of possible internal states from N = 3 to N = 5, CSF and more intriguing SSB phases can be engineered.
N = 5 internal states
Increasing to N = 5 internal states offers more control over the extended BH owing to the independent tunneling parameters γ1 and γ2 (see “Methods” for details). Consequently, more complex configurations of tunneling amplitudes are possible. Here, we consider γ1 < 0 and γ2 > 0 so that \({J}_{1}^{\delta },{J}_{4}^{\delta } \, > \, 0\) and \({J}_{2}^{\delta },{J}_{3}^{\delta } \, < \, 0\) [see Fig. 3b]; this directly yields a geometrically frustrated lattice structure. In what follows we fix Ω/ER = 3.5, however, we verified that for 3.0 < Ω/ER < 4.5 the simulation results change only quantitatively.
We first connect to our N = 3 results by obtaining the phase diagram in the γ1-γ2 plane at half filling (\(\bar{n}=1/2\)) shown in Fig. 6a, and representative values of the correlation functions are shown in (b) evaluated at γ2 = −0.06. Individual phases are identified using the logic employed in the preceding section.
Both subplots were computed for N = 5 internal states at \(\bar{n}=1/2\) filling, Raman coupling Ω/ER = 3.5, a system size of L = 181 sites, and 91 particles. a Phase diagram as a function of detuning Fourier components γ1 and γ2. b Asymptotic correlation functions' (one-body Green’s function \({g}_{{{{\rm{cen}}}}}^{(1)}\), staggered density δN and asymptotic vector order parameter \({\kappa }_{{{{\rm{cen}}}}}^{2}\)) dependence on γ2 for fixed γ1/ER = −0.06, corresponding to the dashed black line in (a).
For small γ1 and γ2 all tunneling coefficients JΔj are small compared to the long-range repulsion VΔj, at half filling this stabilizes a DW1/2 phase. Making γ1 increasingly negative has the predominant effect of introducing a proportionally negative J1. As with the N = 3 case, this simply reduces the ratio \({V}_{\Delta j}/{J}_{1}^{\delta }\) of interaction to kinetic energy and melts the DW1/2 insulator. The resulting conventional SF phase restores the bulk translational symmetry and has quasi-LRO only in \({g}_{{{{\rm{cen}}}}}^{(1)}\) [see Fig. 6b]. In contrast, increasing γ2 introduces significant effective geometrical frustration [see Fig. 3b]. The associated increased kinetic energy still destabilizes the DW1/2 insulator, but favors a CSF where both \({g}_{{{{\rm{cen}}}}}^{(1)}\) and \({\kappa }_{{{{\rm{cen}}}}}^{2}\) are non-zero. As a result this lattice provides a unique opportunity for controlled studies of CSFs.
Period-3 order at 1/3 filling
As is visible in Fig. 4b, the N = 5 long-range interactions are significant even beyond the Δj = 1 nearest neighbor scale. Repulsive interactions that are significant up to VΔj tend to favor ordered phases with a period of Δj + 1. We search for the impact of Δj = 2 next nearest neighbor interactions by reducing the particle density to \(\bar{n}=1/3\), where these interactions would favor a DW1/3 insulator in the limit of zero tunneling. This expectation is confirmed in Fig. 7a where, at small γ2 and for three values of γ1, the expectation value of the on-site number operator has a period-3 oscillatory contribution.
Computed for N = 5 internal states at \(\bar{n}=1/3\) filling, Raman coupling Ω/ER = 3.5, detuning Fourier component γ2/ER = 0.01, a system size of L = 175, and 59 particles. a Local density \(\langle {\hat{n}}_{j}\rangle\); b structure factor S(k); and c asymptotic single particle Green function \({g}_{{{{\rm{cen}}}}}^{(1)}\).
Rather than quantify this structure in terms of a specific correlation function suited only to period-3 density order, we turn to the density-density correlation function \({C}_{j}(\Delta j)=\langle {\hat{n}}_{j+\Delta j}{\hat{n}}_{j}\rangle -\langle {\hat{n}}_{j+\Delta j}\rangle \langle {\hat{n}}_{j}\rangle\) that is sensitive to density fluctuations at a range of Δj and its Fourier transform
the static structure factor. Figure 7b shows that in this parameter regime the structure factor has peaks at k = ±2π/3 indicative of local order associated with spontaneously broken translational symmetry. The sharp peaks present for γ1 = −0.01 are indicative of quasi-LRO, while the broad Lorentzian-like peaks for more negative γ1 suggest an exponential decay of density-density correlations and a lack of SSB. Finally, Fig. 7(c) shows that in the small-negative γ1 SSB case \({g}_{j}^{(1)}(\Delta j)\) vanishes exponentially, thereby confirming the presence of a DW1/3 phase. For more negative γ1 long-range off-diagonal order is established suggesting a normal SF. Thus this transition (as a function of γ1 and for small γ2) from DW1/3 to SF at \(\bar{n}=1/3\) filling is analogous to the DW1/2 to SF transition at \(\bar{n}=1/2\) filling in Fig. 6a. This analysis rigorously proves that the strong long-range repulsion present in our model enables the realization of period three DW insulators, resembling \({{\mathbb{Z}}}_{3}\) Mott insulators predicted in chiral clock models84.
State preparation and detection
The previous sections identified a wide array of quantum states of matter that our setup can access. Although these phases are described by a spinless single band extended BH model, the constituent atoms exist in a dressed state representation, therefore conventional detection and measurement techniques are not effective here. In the following sections we therefore present alternative approaches to detect and prepare low energy states of our model.
Measurement opportunities
This work has focused on distinguishing between conventional SFs, CSFs, and DW solids (where the unit-filled Mott insulator would be DW1) using a range of correlation functions: the spatial density \(\langle {\hat{n}}_{i}\rangle\), the static structure factor S(k), the single particle Green’s function \({g}_{j}^{(1)}(\Delta j)\), and the vector order parameter \({\kappa }_{j}^{2}(\Delta j)\). Given the deeply subwavelength nature of these lattices—λ/6 for N = 3 and λ/10 for N = 5—even today’s highest resolution quantum gas microscopes85,86,87 are unable to directly resolve DW order in \(\langle {\hat{n}}_{i}\rangle\). Fortunately, the underlying physical structure of our system offers a unique opportunity to experimentally access our target observables.
Standard time-of-flight
In this section, we relate momentum distributions observed in time-of-flight (ToF) images with the crystal momentum distribution n(k) characterizing states in the subwavelength lattice. The crystal momentum distribution is directly related to the Fourier transformed total one-body Green’s function \({g}^{(1)}(\Delta j)={\sum}_{j}{g}_{j}^{(1)}(\Delta j)\) via
where \({\hat{b}}_{r}^{{{\dagger}} }(k)\) is the creation operator for a particle with crystal momentum k occupying the r-th band of the subwavelength lattice. Notice that k is dimensionless and the Brillouin zone (BZ) extends over a range of 2π.
By employing the Stern-Gerlach effect during ToF the final observed quantities are the internal state resolved momentum distributions
where \({\hat{\phi }}_{m}^{{{\dagger}} }(q)\) describes the creation of a boson with wavevector q in internal state m. In this expression q has dimensions of inverse length and the subwavelength lattice BZ has an extent of 2NkR; therefore we introduce a factor c = (2π)/(2NkR) to convert from these physical units to the dimensionless units of the discrete lattice.
As we show in Methods, the state resolved momentum distributions are
where \(\tilde{w(q)}\) is the Fourier transformed Wannier function. As such, the crystal momentum distribution, and therefore the single body Green’s function g(1)(Δj), can be obtained from the internal state resolved momentum distributions, but not from the total momentum density \(\hat{\rho (q)}={\sum }_{m}{\rho }_{m}(q)\).
Accessing n(k) provides a powerful tool for distinguishing the many-body phases described in the previous section. Figure 8 shows that, owing to quasi LRO in g(1)(Δj), SF phases (light and dark blue) give rise to sharp peaks that vanish in insulating DW phases (red) where g(1)(Δj) decays exponentially. More specifically the normal SF exhibits a single sharp peak at k = 0, while the CSF has two peaks. The two peaks at incommensurate k in the momentum distribution signal two minima in the dispersion relation. The interaction favors the predominant population of one of the minima and, as a consequence, the system enters a CSF phase with a non-zero local boson current characterized by a finite chirality \(\langle {\hat{\kappa }}_{j}\rangle\).
These crystal momentum distributions provide little information regarding the structure of DW solids, however, higher order correlation functions do. For example, the second order function
provides direct information regarding density order in gapped solids88,89,90.
Staggered readout
An alternate measurement protocol that is unique to this specific type of subwavelength lattice transforms each dressed state into a specific internal atomic state just prior to ToF; as above, this yields N independent momentum distributions, each of which samples every N-th site of the subwavelength lattice (recall that j = n + ℓN, where n is the dressed state index and ℓ is the λ/2 unit cell).
In order to implement this mapping we introduce two new degrees of freedom: (1) a new coupling Ωrf nearly identical with the Raman coupling in Eq. (1) except that it lacks any spatial dependence [this might be implemented with a radio frequency (rf) magnetic field, or with Raman transitions in a co-propagating geometry]; and (2) a detuning proportional to the internal state index, i.e., δm = ΔFm, as would be given by the usual linear Zeeman effect.
Our protocol adiabatically transforms dressed states into internal atomic states, where the adiabatic timescale T is selected to be rapid as compared to the ground-band atomic dynamics, but slow compared to the band splitting. In the following, each step is correspondingly marked in Fig. 9a:
-
(i)
In the first step we quench to zero the detunings δm (used to generate long-range tunneling) in a timescale τ, which is rapid compared to the adiabatic timescale T used for the following steps. An adiabatic timescale would cause detuning-induced Rabi oscillations, which would ruin the one-to-one mapping between dressed and bare states. Since the detuning quench is fast, it is not shown in Fig. 9a. This step returns the system to N interpenetrating but decoupled lattices.
-
(ii)
Next, the spatially uniform coupling Ωrf is ramped on while the the Raman coupling Ω is simultaneously ramped off. This transforms each independent Raman lattice with energy \(-2\Omega \cos (2\pi n/N-2{k}_{{{{\rm{R}}}}}x)\) into a spatially uniform dressed state with energy \(-2{\Omega }_{{{{\rm{rf}}}}}\cos (2\pi n/N)\).
-
(iii)
Lastly, Ωrf is ramped off while the conventional detuning is ramped to a final value of Δ.
We therefore conclude that this process transforms states in the \(\left\vert n\right\rangle\) dressed state into the \(\left\vert j=n\right\rangle\) internal atomic state.
Furthermore, reversing this protocol allows us to transform a 1D quasi-condensate prepared in a single internal atomic state into a corresponding SF state in the subwavelength lattice. Finally, once the SF state is prepared, the CSF and density wave (DW1/2, DW1/3) states can be accessed through adiabatic ramps that modify the system parameters to the required regimes.
Conclusions
In this manuscript, we showed that a recently realized class of 1D subwavelength optical lattices46,47,48,49,50 lead to extended BH models with effective geometric frustration and long-ranged interactions. These lattices Raman couple internal atomic states, and a lattice potential emerges in a dressed state basis whose spacing is reduced by a factor equal to the number of coupled internal states. This configuration features significant interparticle repulsion over the scale of several lattice sites, in the absence of dipolar or Coulomb couplings. On short scales the functional form can be approximated as local repulsion combined with a power-law tail, the exponent of which can range from −2 to −8 depending on the Raman coupling strength and the number of internal states. Tunneling on the subwavelength scale is induced by small deviations from the Raman resonance condition; specifically, tuning the detuning parameters modifies the sign, Peierls phase, and strength of tunneling connecting lattice sites spaced by distances equal to number of coupled internal states.
Controlling the relative sign of the tunneling amplitude at different ranges leads to regimes of effective geometrical frustration. This is in contrast with conventional optical lattice platforms where the range and sign of hopping processes are fixed. Other approaches inducing geometric frustration by periodically modulating one or more parameter in the single particle Hamiltonian91 with interactions, this process induces many-body dephasing which manifests as heating and then atom loss. In the present work, spontaneous emission from the Raman lasers is the only intrinsic heating process; for the example presented here, the associated 1/e lifetimes would be about 500 ms or \(\approx 400\times {(h{V}_{0})}^{-1}\)92, enabling quantum fluctuations to remain stable for long times.
We explored the many-body potentialities offered by this setup with a detailed numerical analysis using matrix product states across a range of interactions and frustration. We found that geometrical frustration favors CSFs with spontaneously broken TR symmetry; these are of great interest both in condensed matter62 and high energy65 physics. The competing presence of strong long-range repulsion favored insulating DWs1,93,94 with spontaneous spatial symmetry breaking.
Although the parameters used in matrix product states analysis are specific to 87Rb, our laser-coupling scheme can be applied to a large variety of atomic species including both bosons and fermions. In the case of fermions, the extended Hubbard model analogous to Eq. (14) still describes spinless particles; owing to Pauli repulsion, only long-range interactions with p-wave (and higher order) character are present70. Embedding such a fermionic sub-wavelength system in a Bose-Einstein condensate would be an intriguing next step. In analogy with materials with phonon mediated interactions between electrons, we expect the emergence of an oscillating bosonic-mediated interaction between fermions of Ruderman-Kittel-Kasuya-Yosida (RKKY)-type95,96,97. Noticeably, cold atomic systems have only been able to present preliminary results on the possible appearance of fermionic-mediated RKKY-type interactions between bosons98,99. Our proposed setup thus poses itself as a relevant source to engineer complex interacting processes analogous of real materials.
Figure 4 modeled interactions in our lattice as an effective power-law, valid for the first three subwavelength lattice sites. A similar procedure can be followed to instead frame these interaction terms as a screened Coulomb interaction, for example, of the repulsive Yukawa form, \({V}_{j}={V}_{{{{\rm{y}}}}}{e}^{-{\gamma }_{{{{\rm{y}}}}}j}/j\) with \({\gamma }_{{{{\rm{y}}}}}=\ln [{V}_{1}/(2{V}_{2})]\). Therefore, our lattice might be applicable for quantum simulation of plasma physics100 and screened electronic systems101. We considered interactions of the SU(N) that were the same for all atomic internal states. This can be a good approximation in some cases (such as 87Rb atoms) and nearly exact in others (such as 87Sr and 173Yb)102. In other cases the interaction strengths can differ greatly; for example, Feshbach resonances can induce significant differences103. This would result in non-local interaction driven tunneling processes like pair hopping that give rise to pair superfluids (PSFs) in bosonic models104,105,106.
In conclusion, our results represent an alternative proposal which can finally shed light on the investigation of long-range frustrated quantum systems.
Methods
Physical parameters
Here we summarize the explicit numerical values of the physical parameters used in our many-body calculations (all taken for 87Rb).
The one-dimensional interaction strength107 is
having used the constant C ≈ 1.4603107, the reduced Planck’s constant ℏ = 1.05 × 10−34 m2kg/s, the atomic mass m = 86.9AMU = 1.44 × 10−25 kg, the F = 2 manifold s-wave scattering length a22 = 95aBohr = 5.02 × 10−9 m108 and transverse confinement length a⊥ = 1.25 × 10−7 m. We used the aforementioned value of g1D for both N = 3 and N = 5. This is reasonable since a11 = 100aBohr and a22 = 95aBohr are fairly close to one another and thus the presented results will not qualitatively change.
The single photon recoil energy
and wavevector kR = 2π/λ both require additional knowledge of the optical wavelength (here λ = 790 nm) used to create the lattice potential.
Interaction Hamiltonian
Here we consider properties of the interaction Hamiltonian for the case of state independent [sometimes called SU(N)] interactions. This makes the interaction Hamiltonian a function of the total local density \({\hat{n}}_{{{{\rm{tot}}}}}(x)={\sum }_{m}{\hat{\phi }}_{m}^{{{\dagger}} }(x){\hat{\phi }}_{m}(x)={\sum }_{n}{\hat{\psi }}_{n}^{{{\dagger}} }(x){\hat{\psi }}_{n}(x)\), which takes the same form in the bare atomic basis (with creation operators \({\psi }_{n}^{{{\dagger}} }(x)\)) and the dressed basis (with creation operators \({\phi }_{m}^{{{\dagger}} }(x)\)). As a result, the interaction Hamiltonian is also unchanged with
where : ⋯ : denotes the normal ordering operation.
Expanding the dressed field operators in terms of lowest band Wannier functions
leads directly to the density-density interaction
that appeared in the Hubbard model [Eq. (14)], where j = Nl + n denotes lattice site index and Δj denotes distance between lattice sites.
Additional terms, such as density-induced tunneling (DIT), can appear when interaction strength \({g}_{j,{j}^{{\prime} }}\) becomes a function j and \({j}^{{\prime} }\). In this case the interaction Hamiltonian
is no longer invariant with respect to a change of basis, and the dressed state Hamiltonian
contains every possible combination of field operators, where
Once again expanding the field operators in the Wannier basis, one obtains
where the sum is over all subwavelength lattice sites and interaction strengths
Coefficients such as Vijki lead to DIT.
In the considered experimental situation, i.e., National Institute of Standards and Technology (NIST) 87Rb cyclic coupling experiment47, interactions are homogeneous at the 0.995 fractional level making DIT terms negligible.
Detuning induced tunneling
As mentioned previously, the synthetic dimension tunneling parameters \({\gamma }_{\Delta n}=| {\gamma }_{\Delta n}| \exp (i{\phi }_{\Delta n})\) are significantly constrained by the properties of the discrete Fourier transform, as well as our restrictions on the allowed detunings:
-
1.
γΔn = γΔn+N owing to the periodicity of Fourier transforms.
-
2.
γ0 = 0, because \({\sum}_{m}{\delta }_{m}=0\).
-
3.
ϕΔn = −ϕ−Δn, because the detunings δm are real valued.
-
4.
For our current subwavelength lattice we focus on the simplification δm = δ−m, making γΔn real-valued and symmetric. Note that this condition is violated for our staggered readout procedure for mapping dressed states to bare states.
For example, for the N = 3 case, the two detuning constraints reduce the number of free degrees of freedom to one, implying that γ1 alone quantifies the detuning induced tunneling. Similarly the N = 5 configuration has two independent degrees of freedom, γ1 and γ2.
For odd N these constraints allow Eq. (8) to be simplified as
where we reindexed the sum to run from − (N − 1)/2 to (N − 1)/2 and combined exponentials at positive and negative m into cosine terms. This shows that δ1, ⋯ , δ(N−1)/2 are the independent degrees of freedom. This expression can be inverted to provide a relation between a desired set of γΔn and the experimental parameters δm.
Variational Gaussian Wannier approximation
Here we derive the approximate Wannier functions yielding the continuous curves in Fig. 4. In brief, these begin with a simple Gaussian approximation for a wavepacket centered on a single lattice site, and then we use a variational ansatz to optimize the width.
We begin with the dimensionless (with energy in units of ER and length in units of \({k}_{{{{\rm{R}}}}}^{-1}\)) Hamiltonian
for a particle moving in a lattice potential of depth s = 4Ω/ER. The second order series expansion around x = 0 yields the harmonic oscillator Hamiltonian
with oscillator frequency \(\omega =2\sqrt{s}\), length ℓ = s−1/4, and ground state wavefunction
Changing Eqs. (10), (11) and (12) to dimensionless quantities, this gives explicit relations for: native tunneling
detuning induced tunneling
and interactions
where we have introduced the dimensionless displacement d = πΔj/N between the center of the Wannier orbitals. We include an expression for the native tunneling \({J}_{N}^{\Omega }\), but note that the Gaussian ansatz leads a nonphysical dependency on the zero of energy (because Gaussian wavepackets at neighboring lattices sites are not orthogonal).
The dashed curves in Fig. 10b plot the on-site interaction computed using this expression along with the numerically computed points. The agreement is poor.
Plotted quantities are computed directly from Wannier functions (markers), from Gaussian variational ansatz (solid curves), or from standard Gaussian ansatz (dotted curves). a Normalized detuning-induced nearest neighbor tunneling \(-{J}_{1}^{\delta }/{\gamma }_{1}\). b Interaction strength VΔj in units of recoil energy ER. c First energy gap ΔE in units of recoil energy ER. d Variational parameter β. Dashed vertical lines mark the critical lattice depth sc = (e/2)2; for arguments below sc the variational parameter β becomes complex.
We improved the accuracy of wavefunctions of this form using the variational principle where we replaced s → β2s and minimized the energy functional
with respect to β. This yields the condition
which is solved by
where W(x) is the notorious Lambert W function109,110; Fig. 10d plots β as a function of laser coupling strength Ω. Because W(x) becomes imaginary for arguments below −1/e, the expression for β is only defined for s > (e/2)2, and for large s, β approaches unity confirming that the standard Gaussian approximation is accurate for very deep lattices. The solid curves in Fig. 10b show the improved on-site interaction energy computed using this correction factor.
First excited band
One can also accurately compute the first energy gap. To do this, we approximate the excited state Wannier function with the first excited harmonic oscillator wavefunction
and in analogy with Eq. (42) we obtain the excited state energy
Subtracting the ground state energy yields the energy gap
Our DMRG computations were performed for Ω/ER > 3, where the impact of higher Bloch bands of the sinusoidal adiabatic potentials are negligible. Using Eq. (46), the first band gap can be approximated as ΔE/ER = 5.8 (direct numerics give ΔE/ER = 5.3ER) which is large compared to the interaction scales (with Vj/ER ≲ 0.25, see Table 1), and the tunneling scales [JΔj/ER ≲ 0.5, see Fig. 10a]. Because the many-body physics under study require temperatures kBT ≲ (J1, V0), thermal excitations are also negligible.
Relation to continuum degrees of freedom
Here we relate the momentum distribution observed in ToF images with the crystal momentum distribution of states in the subwavelength lattice. By employing the Stern-Gerlach effect during ToF the final observed quantities are the internal state resolved momentum density operators
where \({\hat{\phi }}_{m}^{{{\dagger}} }(q)\) describes the creation of a boson at momentum q in internal state m. In terms of the continuum dressed state field operators in Eq. (3) this becomes
leading to the final expression in the discrete Wannier basis
having made use of
where
Here we again assume that only the ground band (r = 0) is relevant and thus we write
where \({\phi }_{\ell ,n,m}(q)=q\left(n+N\ell \right)-2{k}_{{{{\rm{R}}}}}nm\) and we have made use of the periodicity of the exponential function \(\exp (2\pi inm/N)=\exp (2\pi i(n+N\ell )m/N)\). Implementing the dressed state transformation gives
Because we defined crystal momentum to be dimensionless with a BZ 2π in extent. We see that this expression links momentum q in state m with a crystal momentum \(2\pi \left(2/(2N{k}_{{{{\rm{R}}}}})-m/N\right)\) where 2NkR is the extent of the BZ in physical units, shifted by 2πm/N which can be interpreted as a result of the recoil kick imparted by each two-photon Raman transition.
This result connects the internal-state resolved momentum density operator and the crystal momentum density
and shows that the probability is split equally between the N internal states.
Data availability
Data of the study can be obtained from the corresponding authors upon request.
Code availability
Code used in the study can be obtained from the corresponding authors upon request.
References
Haskell, B. & Sedrakian, A. Superfluidity and Superconductivity in Neutron Stars 401–454 (Springer International Publishing, 2018).
McIntosh, J. S., Park, S. C. & Rawitscher, G. H. Nucleus-nucleus long-range interaction potential in elastic scattering. Phys. Rev. 134, B1010–B1021 (1964).
Feinberg, G. & Sucher, J. Long-range forces from neutrino-pair exchange. Phys. Rev. 166, 1638–1644 (1968).
Yukawa, H. On the interaction of elementary particles. I. Proc. Phys. Math. Soc. Jpn. 3rd Ser. 17, 48–57 (1935).
Bohm, D. & Pines, D. A collective description of electron interactions: III. coulomb interactions in a degenerate electron gas. Phys. Rev. 92, 609–625 (1953).
Defenu, N. et al. Long-range interacting quantum systems. Rev. Mod. Phys. 95, 035002 (2023).
Chiara, G. D. & Sanpera, A. Genuine quantum correlations in quantum many-body systems: a review of recent progress. Rep. Prog. Phys. 81, 074002 (2018).
Eisert, J., van den Worm, M., Manmana, S. R. & Kastner, M. Breakdown of quasilocality in long-range quantum lattice models. Phys. Rev. Lett. 111, 260401 (2013).
Gong, Z.-X., Foss-Feig, M., Brandão, F. G. S. L. & Gorshkov, A. V. Entanglement area laws for long-range interacting systems. Phys. Rev. Lett. 119, 050501 (2017).
Bruno, P. Absence of spontaneous magnetic order at nonzero temperature in one- and two-dimensional Heisenberg and XY systems with long-range interactions. Phys. Rev. Lett. 87, 137203 (2001).
Peter, D., Müller, S., Wessel, S. & Büchler, H. P. Anomalous behavior of spin systems with dipolar interactions. Phys. Rev. Lett. 109, 025303 (2012).
Maghrebi, M. F., Gong, Z.-X. & Gorshkov, A. V. Continuous symmetry breaking in 1d long-range interacting quantum systems. Phys. Rev. Lett. 119, 023001 (2017).
Hauke, P. & Tagliacozzo, L. Spread of correlations in long-range interacting quantum systems. Phys. Rev. Lett. 111, 207202 (2013).
Richerme, P. et al. Non-local propagation of correlations in quantum systems with long-range interactions. Nature 511, 198–201 (2014).
Lacroix, C., Mendels, P. & Mila, F. Introduction to frustrated magnetism. Springer Series in Solid-State Sciences (Springer, 2011).
Kane, C. L. & Mele, E. J. Quantum spin hall effect in graphene. Phys. Rev. Lett. 95, 226801 (2005).
Qi, X.-L. & Zhang, S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057–1110 (2011).
Balents, L. Spin liquids in frustrated magnets. Nature 464, 199–208 (2010).
Majumdar, C. K. & Ghosh, D. K. On next-nearest-neighbor interaction in linear chain. I. J. Math. Phys. 10, 1388–1398 (1969).
Haldane, F. D. M. Spontaneous dimerization in the \(s=\frac{1}{2}\) Heisenberg antiferromagnetic chain with competing interactions. Phys. Rev. B 25, 4925–4928 (1982).
Läuchli, A. M. Numerical Simulations of Frustrated Systems 481–511 (Springer Berlin Heidelberg, 2011).
Crosswhite, G. M., Doherty, A. C. & Vidal, G. Applying matrix product operators to model systems with long-range interactions. Phys. Rev. B 78, 035116 (2008).
Ciorciaro, L. et al. Kinetic magnetism in triangular moirématerials. Nature 623, 509–513 (2023).
Scholl, P. et al. Quantum simulation of 2D antiferromagnets with hundreds of Rydberg atoms. Nature 595, 233–238 (2021).
Semeghini, G. et al. Probing topological spin liquids on a programmable quantum simulator. Science 374, 1242–1247 (2021).
Gross, C. & Bloch, I. Quantum simulations with ultracold atoms in optical lattices. Science 357, 995–1001 (2017).
Lewenstein, M. et al. Ultracold atomic gases in optical lattices: mimicking condensed matter physics and beyond. Adv. Phys. 56, 243–379 (2007).
Lahaye, T., Menotti, C., Santos, L., Lewenstein, M. & Pfau, T. The physics of dipolar bosonic quantum gases. Rep. Prog. Phys. 72, 126401 (2009).
Chomaz, L. et al. Dipolar physics: a review of experiments with magnetic quantum gases. Rep. Prog. Phys. 86, 026401 (2022).
Langen, T., Valtolina, G., Wang, D. & Ye, J. Quantum state manipulation and cooling of ultracold molecules. Nat. Phys. 20, 702–712 (2024).
Farokh Mivehvar, T. D., Francesco, P. & Ritsch, H. Cavity qed with quantum gases: new paradigms in many-body physics. Adv. Phys. 70, 1–153 (2021).
Struck, J. et al. Quantum simulation of frustrated classical magnetism in triangular optical lattices. Science 333, 996–999 (2011).
Struck, J. et al. Engineering ising-xy spin-models in a triangular lattice using tunable artificial gauge fields. Nat. Phys. 9, 738–743 (2013).
Leung, T.-H. et al. Interaction-enhanced group velocity of bosons in the flat band of an optical kagome lattice. Phys. Rev. Lett. 125, 133001 (2020).
Wang, X.-Q. et al. Evidence for an atomic chiral superfluid with topological excitations. Nature 596, 227–231 (2011).
Brown, C. D. et al. Direct geometric probe of singularities in band structure. Science 377, 1319–1322 (2022).
Saugmann, P. et al. Route toward classical frustration and band flattening via optical lattice distortion. Phys. Rev. A 106, L041302 (2022).
Damski, B. et al. Quantum gases in trimerized kagomé lattices. Phys. Rev. A 72, 053612 (2005).
Eckardt, A. et al. Frustrated quantum antiferromagnetism with ultracold bosons in a triangular lattice. EPL 89, 10010 (2010).
Zhang, T. & Jo, G.-B. One-dimensional sawtooth and zigzag lattices for ultracold atoms. Sci. Rep. 5, 16044 (2015).
Yamamoto, D., Fukuhara, T. & Danshita, I. Frustrated quantum magnetism with Bose gases in triangular optical lattices at negative absolute temperatures. Commun. Phys. 3, 56 (2020).
Anisimovas, E. et al. Semisynthetic zigzag optical lattice for ultracold bosons. Phys. Rev. A 94, 063632 (2016).
Cabedo, J., Claramunt, J., Mompart, J., Ahufinger, V. & Celi, A. Effective triangular ladders with staggered flux from spin-orbit coupling in 1d optical lattices. Eur. Phys. J. D 74, 123 (2020).
Barbiero, L., Cabedo, J., Lewenstein, M., Tarruell, L. & Celi, A. Frustrated magnets without geometrical frustration in bosonic flux ladders. Phys. Rev. Res. 5, L042008 (2023).
Baldelli, N., Cabrera, C. R., Julià-Farré, S., Aidelsburger, M. & Barbiero, L. Frustrated extended Bose-Hubbard model and deconfined quantum critical points with optical lattices at the antimagic wavelength. Phys. Rev. Lett. 132, 153401 (2024).
Wang, Y. et al. Dark state optical lattice with a subwavelength spatial structure. Phys. Rev. Lett. 120, 083601 (2018).
Anderson, R. P. et al. Realization of a deeply subwavelength adiabatic optical lattice. Phys. Rev. Res. 2, 013149 (2020).
Lauria, P., Kuo, W.-T., Cooper, N. R. & Barreiro, J. T. Experimental realization of a fermionic spin-momentum lattice. Phys. Rev. Lett. 128, 245301 (2022).
Li, C.-H. et al. Bose-einstein condensate on a synthetic topological hall cylinder. PRX Quantum 3, 010316 (2022).
Burba, D., Račiūnas, M., Spielman, I. B. & Juzeliūnas, G. Topological charge pumping with subwavelength Raman lattices. Phys. Rev. A 107, 023309 (2023).
Han, J. H., Kang, J. H. & Shin, Y. Band gap closing in a synthetic hall tube of neutral fermions. Phys. Rev. Lett. 122, 065303 (2019).
Yan, Y., Zhang, S.-L., Choudhury, S. & Zhou, Q. Emergent periodic and quasiperiodic lattices on surfaces of synthetic hall tori and synthetic hall cylinders. Phys. Rev. Lett. 123, 260405 (2019).
Zhang, R., Yan, Y. & Zhou, Q. Localization on a synthetic hall cylinder. Phys. Rev. Lett. 126, 193001 (2021).
Liang, Q.-Y. et al. Coherence and decoherence in the Harper-Hofstadter model. Phys. Rev. Res. 3, 023058 (2021).
Fabre, A., Bouhiron, J.-B., Satoor, T., Lopes, R. & Nascimbene, S. Laughlin’s topological charge pump in an atomic hall cylinder. Phys. Rev. Lett. 128, 173202 (2022).
Baier, S. et al. Extended Bose-Hubbard models with ultracold magnetic atoms. Science 352, 201–205 (2016).
Su, L. et al. Dipolar quantum solids emerging in a Hubbard quantum simulator. Nature 622, 724–729 (2023).
Guardado-Sanchez, E. et al. Quench dynamics of a fermi gas with strong nonlocal interactions. Phys. Rev. X 11, 021036 (2021).
Christakis, L. et al. Probing site-resolved correlations in a spin system of ultracold molecules. Nature 614, 64–69 (2023).
Carroll, A. N. et al. Observation of generalized t-j spin dynamics with tunable dipolar interactions (2024).
Schollwöck, U. The density-matrix renormalization group in the age of matrix product states. Ann. Phys. 326, 96–192 (2011).
Li, X. & Liu, W. V. Physics of higher orbital bands in optical lattices: a review. Rep. Prog. Phys. 79, 116401 (2016).
Di Liberto, M. & Goldman, N. Chiral orbital order of interacting bosons without higher bands. Phys. Rev. Res. 5, 023064 (2023).
Goldman, N. et al. Floquet-engineered nonlinearities and controllable pair-hopping processes: from optical Kerr cavities to correlated quantum matter. PRX Quantum 4, 040327 (2023).
Son, D. T. & Surówka, P. Hydrodynamics with triangle anomalies. Phys. Rev. Lett. 103, 191601 (2009).
Ozawa, T. & Price, H. M. Topological quantum matter in synthetic dimensions. Nat. Rev. Phys. 1, 349–357 (2019).
Jiménez-García, K. & Spielman, I. B. Annual Review of Cold Atoms and Molecules Vol. 2, 145–191 (World Scientific, 2013).
Widera, A. et al. Precision measurement of spin-dependent interaction strengths for spin-1 and spin-2 87 rb atoms. N. J. Phys. 8, 152 (2006).
Spielman, I. B. Raman processes and effective gauge potentials. Phys. Rev. A 79, 063613 (2009).
Williams, R. A. et al. Synthetic partial waves in ultracold atomic collisions. Science 335, 314–317 (2012).
Korbmacher, H., Domínguez-Castro, G. A., Li, W.-H., Zakrzewski, J. & Santos, L. Transversal effects on the ground state of hard-core dipolar bosons in one-dimensional optical lattices. Phys. Rev. A 107, 063307 (2023).
Monroe, C. et al. Programmable quantum simulations of spin systems with trapped ions. Rev. Mod. Phys. 93, 025001 (2021).
Lühmann, D.-S., Jürgensen, O. & Sengstock, K. Multi-orbital and density-induced tunneling of bosons in optical lattices. N. J. Phys. 14, 033021 (2012).
Dutta, O. et al. Non-standard hubbard models in optical lattices: a review. Rep. Prog. Phys. 78, 066001 (2015).
Dalla Torre, E. G., Berg, E. & Altman, E. Hidden order in 1d Bose insulators. Phys. Rev. Lett. 97, 260401 (2006).
Dhar, A. et al. Bose-Hubbard model in a strong effective magnetic field: emergence of a chiral mott insulator ground state. Phys. Rev. A 85, 041602(R) (2012).
Greschner, S., Santos, L. & Vekua, T. Ultracold bosons in zig-zag optical lattices. Phys. Rev. A 87, 033609 (2013).
Zaletel, M. P., Parameswaran, S. A., Rüegg, A. & Altman, E. Chiral bosonic mott insulator on the frustrated triangular lattice. Phys. Rev. B 89, 155142 (2014).
Romen, C. & Läuchli, A. M. Chiral mott insulators in frustrated Bose-Hubbard models on ladders and two-dimensional lattices: a combined perturbative and density matrix renormalization group study. Phys. Rev. B 98, 054519 (2018).
Fraxanet, J. et al. Topological quantum critical points in the extended Bose-Hubbard model. Phys. Rev. Lett. 128, 043402 (2022).
Singha Roy, S., Carl, L. & Hauke, P. Genuine multipartite entanglement in a one-dimensional Bose-Hubbard model with frustrated hopping. Phys. Rev. B 106, 195158 (2022).
Halati, C.-M. & Giamarchi, T. Bose-Hubbard triangular ladder in an artificial gauge field. Phys. Rev. Res. 5, 013126 (2023).
Chanda, T., Barbiero, L., Lewenstein, M., Mark, M. J. & Zakrzewski, J. Recent progress on quantum simulations of non-standard Bose-Hubbard models (2024).
Samajdar, R., Choi, S., Pichler, H., Lukin, M. D. & Sachdev, S. Numerical study of the chiral \({{\mathbb{z}}}_{3}\) quantum phase transition in one spatial dimension. Phys. Rev. A 98, 023614 (2018).
Bakr, W. S., Gillen, J. I., Peng, A., Fölling, S. & Greiner, M. A quantum gas microscope for detecting single atoms in a Hubbard-regime optical lattice. Nature 462, 74 EP – (2009).
Sherson, J. F. et al. Single-atom-resolved fluorescence imaging of an atomic Mott insulator. Nature 467, 68–72 (2010).
Asteria, L., Zahn, H. P., Kosch, M. N., Sengstock, K. & Weitenberg, C. Quantum gas magnifier for sub-lattice-resolved imaging of 3d quantum systems. Nature 599, 571–575 (2021).
Altman, E., Demler, E. & Lukin, M. D. Probing many-body states of ultracold atoms via noise correlations. Phys. Rev. A 70, 013603 (2004).
Fölling, S. et al. Spatial quantum noise interferometry in expanding ultracold atom clouds. Nature 434, 481–484 (2005).
Spielman, I. B., Phillips, W. D. & Porto, J. V. Mott-insulator transition in a two-dimensional atomic Bose gas. Phys. Rev. Lett. 98, 080404 (2007).
Eckardt, A. Colloquium: atomic quantum gases in periodically driven optical lattices. Rev. Mod. Phys. 89, 011004 (2017).
Lin, Y. J., Compton, R. L., Jiménez-García, K., Porto, J. V. & Spielman, I. B. Synthetic magnetic fields for ultracold neutral atoms. Nature 462, 628–632 (2009).
Balibar, S. The enigma of supersolidity. Nature 464, 176–182 (2010).
Grüner, G. The dynamics of charge-density waves. Rev. Mod. Phys. 60, 1129–1181 (1988).
Ruderman, M. A. & Kittel, C. Indirect exchange coupling of nuclear magnetic moments by conduction electrons. Phys. Rev. 96, 99–102 (1954).
Kasuya, T. A theory of metallic ferro- and antiferromagnetism on Zener’s model. Prog. Theor. Phys. 16, 45–57 (1956).
Yosida, K. Magnetic properties of cu-mn alloys. Phys. Rev. 106, 893–898 (1957).
DeSalvo, B. J., Patel, K., Cai, G. & Chin, C. Observation of fermion-mediated interactions between bosonic atoms. Nature 568, 61–64 (2019).
Edri, H., Raz, B., Matzliah, N., Davidson, N. & Ozeri, R. Observation of spin-spin fermion-mediated interactions between ultracold bosons. Phys. Rev. Lett. 124, 163401 (2020).
Zhang, C., Capogrosso-Sansone, B., Boninsegni, M., Prokof’ev, N. V. & Svistunov, B. V. Superconducting transition temperature of the Bose one-component plasma. Phys. Rev. Lett. 130, 236001 (2023).
Batista, C. D., Lin, S.-Z., Hayami, S. & Kamiya, Y. Frustration and chiral orderings in correlated electron systems. Rep. Prog. Phys. 79, 084504 (2016).
Fallani, L. Multicomponent spin mixtures of two-electron fermions (2023).
Chin, C., Grimm, R., Julienne, P. & Tiesinga, E. Feshbach resonances in ultracold gases. Rev. Mod. Phys. 82, 1225–1286 (2010).
Zhou, X.-F., Zhang, Y.-S. & Guo, G.-C. Pair tunneling of bosonic atoms in an optical lattice. Phys. Rev. A 80, 013605 (2009).
Takayoshi, S., Katsura, H., Watanabe, N. & Aoki, H. Phase diagram and pair Tomonaga-Luttinger liquid in a Bose-Hubbard model with flat bands. Phys. Rev. A 88, 063613 (2013).
Tovmasyan, M., van Nieuwenburg, E. P. L. & Huber, S. D. Geometry-induced pair condensation. Phys. Rev. B 88, 220510(R) (2013).
Olshanii, M. Atomic scattering in the presence of an external confinement and a gas of impenetrable bosons. Phys. Rev. Lett. 81, 938–941 (1998).
Egorov, M. et al. Measurement of s-wave scattering lengths in a two-component bose-einstein condensate. Phys. Rev. A 87, 053614 (2013).
Lambert, J. H. Observationes Variae in Mathesin Puram (Acta Helveticae physico-mathematico-anatomico-botanico-medica, 1758).
NIST Digital Library of Mathematical Functions. https://dlmf.nist.gov/, Release 1.2.1. (eds Olver, W. J., Olde Daalhuis, A. B., Lozier, D. W., Schneider, B. I., Boisvert, R. F., Clark, C. W., Miller, B. R., Saunders, B. V., Cohl, H. S. & McClain, M. A) (2024).
Acknowledgements
We thank F. Schreck and P. Thekkeppat for discussions. L.B. acknowledges funding from Politecnico di Torino, starting package Grant No. 54 RSG21BL01 and from the Italian MUR (PRIN DiQut Grant No. 2022523NA7). I.B.S. was partially supported by the National Institute of Standards and Technology; the National Science Foundation through the Quantum Leap Challenge Institute for Robust Quantum Simulation (grant OMA-2120757); and the Air Force Office of Scientific Research Multidisciplinary University Research Initiative “RAPSYDY in Q” (FA9550-22-1-0339). D.B. and G.J. acknowledge support from the Lithuanian Research Council (Grant No. S-MIP-24-97). D.B. used resources at the High Performance Computing Center (HPCC), “HPC Sauletekis” in Vilnius University, Faculty of Physics.
Author information
Authors and Affiliations
Contributions
L.B. initiated the work. D.B. and G.J. performed theoretical analysis on the subwavelength Raman lattice. D.B. numerically calculated Wannier functions and Hamiltonian matrix elements. D.B. and G.J. explicitly derive the extended BH Hamiltonian. L.B. performed all of the quantum many-body analysis with DMRG calculations (asymptotic correlation functions, phase diagram, structure factor, etc.) on the extended BH Hamiltonian. I.B.S. considered experimental feasibility and implementation questions. I.B.S. also studied experimental state preparation and measurement details. D.B. numerically calculated dynamics during staggered readout. I.B.S. proposed variational Gaussian Wannier approximation.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Communications Physics thanks Thomas Bilitewski and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Burba, D., Juzeliūnas, G., Spielman, I.B. et al. Many-body phases from effective geometrical frustration and long-range interactions in a subwavelength lattice. Commun Phys 8, 141 (2025). https://doi.org/10.1038/s42005-025-02043-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s42005-025-02043-y