---
title: Two-Component Rydberg-Dressed BECs
url: https://www.emergentmind.com/topics/two-component-rydberg-dressed-bose-einstein-condensates
type: topic
---

# Two-Component Rydberg-Dressed BECs

Two-component Rydberg-dressed Bose-Einstein condensates are binary condensates in which each component acquires an effective isotropic soft-core van-der-Waals interaction through weak off-resonant coupling to a high-lying Rydberg state. In mean-field descriptions, these systems combine contact interactions with nonlocal intra- and inter-component soft-core repulsion, and in quasi-two-dimensional settings they may additionally include Raman-induced spin-orbit coupling (SOC), detuning, and harmonic confinement. Recent work has identified two complementary regimes: a quasi-2D balanced mixture with Raman-induced SOC, where the interplay of SOC, contact interactions, and soft-core Rydberg forces yields topological supersolids, chiral symmetry breaking, finite-window superfluidity, and a dissipative continuous time crystal [2501.06555]; and a three-dimensional bosonic mixture, where the competition between inter- and intra-component blockade radii yields ionic and segregated supersolid structures with face-centered cubic, simple cubic, tubular, and planar order [2507.16277].

## 1. Microscopic formulation and interaction structure

In the quasi-2D setting, the zero-temperature energy functional for a two-component condensate with Raman-induced SOC and Rydberg dressing is
$$
E[\psi_1,\psi_2]
=
\int d^2r\, \psi^\dagger(\mathbf r)\bigl(\hat H_0+V(\mathbf r)+V_{\rm so}\bigr)\psi(\mathbf r)
+\frac12\sum_{j,l=1,2}\int d^2r\, g_{jl}|\psi_j(\mathbf r)|^2|\psi_l(\mathbf r)|^2
+\frac12\sum_{j,l=1,2}\int d^2r\, d^2r'\,
\psi_j^*(\mathbf r)\psi_l^*(\mathbf r')U_{jl}(\mathbf r-\mathbf r')\psi_l(\mathbf r')\psi_j(\mathbf r),
$$
with $\psi=(\psi_1,\psi_2)^T$, $\hat H_0=-(\hbar^2/2m)\nabla^2$, and $V(\mathbf r)=(m\omega_\perp^2/2)(x^2+y^2)$. The Raman-induced SOC term is
$$
V_{\rm so}=-i\hbar\kappa(\sigma_x\partial_x+\sigma_y\partial_y)+\frac{\Omega_R}{2}\sigma_z-\frac{\delta}{2}\sigma_x,
$$
where $\kappa$ is the SOC strength, $\Omega_R$ the Raman coupling, and $\delta$ the detuning. The soft-core interaction is
$$
U_{jl}(\mathbf r)=\frac{\tilde V_6}{R_c^6+|\mathbf r|^6},
$$
with $\tilde V_6=(\Omega/2\Delta)^4 C_6/(\sqrt{2\pi}a_z)$ and $R_c=(C_6/2\hbar\Delta)^{1/6}$ [2501.06555].

In three spatial dimensions, the microscopic model is written in second-quantized form as
$$
H=
\sum_{\alpha=1,2}\int d^3r\, \hat\psi_\alpha^\dagger(\mathbf r)\left[-\frac{\hbar^2\nabla^2}{2M}\right]\hat\psi_\alpha(\mathbf r)
+\frac12\sum_{\alpha,\beta}\int d^3r\, d^3r'\,
\hat\psi_\alpha^\dagger(\mathbf r)\hat\psi_\beta^\dagger(\mathbf r')V_{\alpha\beta}(|\mathbf r-\mathbf r'|)\hat\psi_\beta(\mathbf r')\hat\psi_\alpha(\mathbf r)
+\frac12\sum_{\alpha,\beta}g_{\alpha\beta}\int d^3r\, \hat\psi_\alpha^\dagger(\mathbf r)\hat\psi_\beta^\dagger(\mathbf r)\hat\psi_\beta(\mathbf r)\hat\psi_\alpha(\mathbf r),
$$
where $g_{\alpha\beta}=4\pi\hbar^2 a_s^{(\alpha\beta)}/M$ and
$$
V_{\alpha\beta}(r)=\frac{\tilde C_6^{(\alpha\beta)}}{1+\bigl(r/R_c^{(\alpha\beta)}\bigr)^6}.
$$
For most of that analysis, the couplings are chosen fully symmetric in strength, $g\equiv g_{11}=g_{22}=g_{12}$ and $\tilde C_6\equiv \tilde C_6^{(11)}=\tilde C_6^{(22)}=\tilde C_6^{(12)}$, while allowing $R_c^{(12)}\neq R_c$ to expose specifically binary effects [2507.16277].

The common structural feature is the soft-core nonlocality. In both dimensions, the interaction saturates within a blockade radius and crosses over to a van-der-Waals tail outside it. This makes the relevant competition not merely one of interaction strength, but also of interaction range and of how that range differs between intra- and inter-component channels.

## 2. Mean-field reduction and control parameters

In the quasi-2D problem, the natural units are the blockade radius $R_c$ and the time scale $\tau=mR_c^2/\hbar$. After the rescalings $t\to t/\tau$, $\mathbf r\to \mathbf r/R_c$, and $\psi_j\to \sqrt{R_c^2/N}\,\psi_j$, the coupled Gross–Pitaevskii equations become
$$
i\partial_t\psi_1=
\left[
-\frac12\nabla^2+V(\mathbf r)+\beta_{11}|\psi_1|^2+\beta_{12}|\psi_2|^2+\frac{\Omega_R}{2}
+\int U_{11}(\mathbf r-\mathbf r')|\psi_1(\mathbf r')|^2\, d^2r'
+\int U_{12}(\mathbf r-\mathbf r')|\psi_2(\mathbf r')|^2\, d^2r'
\right]\psi_1
-\kappa(i\partial_x+\partial_y)\psi_2-\frac{\delta}{2}\psi_2,
$$
$$
i\partial_t\psi_2=
\left[
-\frac12\nabla^2+V(\mathbf r)+\beta_{22}|\psi_2|^2+\beta_{12}|\psi_1|^2-\frac{\Omega_R}{2}
+\int U_{22}(\mathbf r-\mathbf r')|\psi_2(\mathbf r')|^2\, d^2r'
+\int U_{12}(\mathbf r-\mathbf r')|\psi_1(\mathbf r')|^2\, d^2r'
\right]\psi_2
-\kappa(i\partial_x-\partial_y)\psi_1-\frac{\delta}{2}\psi_1,
$$
with $V(\mathbf r)=\frac12\omega_c^2(x^2+y^2)$, $\omega_c=\omega_\perp\tau$, $\beta_{jj}=2\sqrt{2\pi}(a_jN)/(a_zR_c)$, $\beta_{12}=2\sqrt{2\pi}(a_{12}N)/(a_zR_c)$, and $U_{jl}(\mathbf r)=C_6/(1+|\mathbf r|^6)$ with $C_6\equiv mN\tilde V_6/(\hbar^2R_c^6)$ [2501.06555].

In the three-dimensional mixture, the Gross–Pitaevskii approximation replaces $\hat\psi_\alpha(\mathbf r)\to \Psi_\alpha(\mathbf r)$, normalized by $\int |\Psi_\alpha|^2\,d^3r=N_\alpha$. The energy functional is
$$
E[\Psi_1,\Psi_2]
=
\sum_{\alpha=1,2}\int d^3r
\left[
\frac{\hbar^2}{2M}|\nabla\Psi_\alpha(\mathbf r)|^2+\frac12 g_{\alpha\alpha}|\Psi_\alpha(\mathbf r)|^4
\right]
+g_{12}\int d^3r\, |\Psi_1(\mathbf r)|^2|\Psi_2(\mathbf r)|^2
+\frac12\sum_{\alpha,\beta=1}^2\int d^3r\, d^3r'\,
|\Psi_\alpha(\mathbf r)|^2V_{\alpha\beta}(|\mathbf r-\mathbf r'|)|\Psi_\beta(\mathbf r')|^2.
$$
That work studies the zero-temperature phase diagram in the plane $(a_s/R_c,\, R_c^{(12)}/R_c)$ at fixed $\rho R_c^3=1.5\times 10^3$ and $\tilde C_6=0.15\,\hbar^2/(MR_c^2)$ [2507.16277].

A central organizing principle in the 3D mixture is the competition of ranges. If $R_c^{(12)}<R_c$, the intra-component repulsion is longer ranged than the cross-component repulsion. At weak $g$, each species prefers to crystallize on its own cluster lattice and to interleave, yielding “ionic” crystals. At larger $g$, the contact-energy cost of large local density suppresses strongly localized clusters and lowers the effective dimensionality of symmetry breaking, producing tubular-segregated or planar-segregated density profiles. This suggests that blockade-radius mismatch functions as a structural control parameter, rather than merely as a perturbation.

## 3. Quasi-two-dimensional phase structure with Raman-induced SOC

Imaginary-time evolution of the quasi-2D equations yields two ground-state phase diagrams: one in the $(\tilde V_6,\kappa)$ plane at fixed $\Omega_R=10$, and one in the $(\tilde V_6,\Omega_R)$ plane at fixed $\kappa=4$. The identified phases include half-quantum vortex, stripe supersolid, toroidal stripe with a central Anderson–Toulouse coreless vortex, checkerboard supersolid, mirror-symmetric supersolid, chiral supersolid, and standing-wave supersolid [2501.06555].

| Phase | Real-space characterization | Momentum / topology |
|---|---|---|
| Half-quantum vortex (HQV) | One component hosts a vortex, the other is uniform | Single peak at $k=0$; uniform spin texture |
| Stripe supersolid (SSS) | One-dimensional stripe modulation | Two peaks at $\pm k_y$; stripe $S_x,S_y$ pattern |
| Toroidal stripe (TS) | Concentric rings; component 2 fills the core of rotating component 1 | Ring plus side rings; Anderson–Toulouse vortex |
| Checkerboard supersolid (CBSS) | Two-dimensional checkerboard of droplets | Four-fold lattice of peaks on a ring; hidden vortex-antivortex pairs in each droplet |
| Mirror-symmetric supersolid (MSSS) | Phase-separated supersolid with mirror symmetry at $\Omega_R=0$ and strong $\tilde V_6$ | Irregular discrete peaks; skyrmion lattice for $y<0$ and antiskyrmion lattice for $y>0$, with $Q=\pm1$ |
| Chiral supersolid (CSS) | Concentric phase-separated rings with net circulation | Seven-point triangular lattice; helical antiskyrmion lattice with $Q=-1$ |
| Standing-wave supersolid (SWSS) | Ring of stripe standing waves | Continuous ring structure |

The chiral supersolid is singled out by the parameter regime $\Omega_R\gtrsim 10$, $\tilde V_6\gtrsim 10^3$, and $\kappa\approx 4$. Its mechanism is stated explicitly: Raman SOC breaks the mirror symmetry of the mirror-symmetric supersolid and drives all vortices to co-circulate, thereby producing net chirality, while long-range soft-core repulsion stabilizes concentric droplet rings. In the spin sector, the signature is a helical antiskyrmion lattice with $Q=-1$ and no Neél walls. In momentum space, the signature is a discrete triangular structure factor. The same study emphasizes a strong contrast with Rashba SOC, for which the supersolid phase is mirror symmetric and contains a skyrmion-antiskyrmion lattice pair rather than a chiral helical antiskyrmion lattice.

The phase characterization is not based solely on density profiles. The Fourier structure factor $S(\mathbf k)$ distinguishes stripe, checkerboard, and triangular ordering; the topological charge density
$$
q(\mathbf r)=\frac{1}{4\pi}\,\mathbf S\cdot(\partial_x\mathbf S\times \partial_y\mathbf S),
$$
with $Q=\int q(\mathbf r)\, d^2r$, distinguishes skyrmion from antiskyrmion textures; and the averaged angular momentum $\langle L_z\rangle$ exhibits jumps across the TS$\to$CBSS and CSS$\to$SWSS transitions. A plausible implication is that topology, crystalline order, and circulation are tightly locked in the SOC-dressed regime rather than constituting separable diagnostics.

## 4. Three-dimensional ionic and segregated supersolids

For $R_c^{(12)}/R_c\lesssim 1$, four principal three-dimensional phases appear in the $(a_s/R_c,\, R_c^{(12)}/R_c)$ phase diagram: ionic FCC crystal, ionic SC crystal, tubular-segregated phase, and planar-segregated phase [2507.16277].

| Phase | Density organization | Structural interpretation |
|---|---|---|
| Ionic FCC crystal | Each species forms an FCC lattice displaced by $\delta=(a/2,a/2,a/2)$ | Total density is a simple-cubic arrangement with a two-site basis |
| Ionic SC | Each species forms a simple-cubic lattice displaced by $\delta'=(a/2,a/2,a/2)$ | Total density becomes a body-centered cubic crystal |
| Tubular-segregated (ST) | Columns of one component arranged on a 2D square lattice; the other occupies interstitial tubes | Symmetry breaking in two directions, with $\sigma_z\to\infty$ in the Gaussian ansatz |
| Planar-segregated (SP) | Alternating planes of the two components | Symmetry breaking in one direction, with $\sigma_{y,z}\to\infty$ |

The ionic FCC phase occupies the lower-left region of the phase diagram. The component densities are represented by Gaussian lattice sums
$$
n_\alpha(\mathbf r)=\sum_{l,m,n}\exp\!\left[-\frac{(\mathbf r-\mathbf R^{(\alpha)}_{lmn})^2}{\sigma^2}\right],
$$
with $\mathbf R^{(1)}_{lmn}=(l,m,n)a$ and $\mathbf R^{(2)}_{lmn}=\mathbf R^{(1)}_{lmn}+\delta$. The ionic SC phase appears at intermediate $R_c^{(12)}/R_c$ and weak $g$; each component sits on its own simple-cubic lattice with the same spacing $a$ and relative displacement $(a/2,a/2,a/2)$, so that the combined density is body-centered cubic. As $g$ increases, the system transitions first to the ST phase, where one finds columns of one component running along, for example, $z$ and arranged on a square lattice in the transverse plane, and then to the SP phase, where alternating planes of the two components appear perpendicular to, for example, the $x$ axis.

The phase boundaries are obtained from full numerical imaginary-time evolution of coupled Gross–Pitaevskii equations and confirmed by a Gaussian variational calculation. Along the horizontal cut $R_c^{(12)}/R_c=0.95$, the typical critical parameters are FCC$\to$ST/SC at $a_s/R_c\simeq 5\times 10^{-4}$ and ST$\to$SP at $a_s/R_c\simeq 2.3\times 10^{-3}$. At $R_c^{(12)}/R_c=1$, the binary problem reduces to the single-component Rydberg-dressed condensate, recovering the sequence FCC for $a_s/R_c<1.2\times 10^{-4}$, then BCC through a first-order transition, and then an unmodulated plane-wave state through a first-order transition at $a_s/R_c\simeq 1.51\times 10^{-3}$.

These results make clear that the 3D mixture supports not only conventional crystalline cluster arrangements but also lower-dimensional segregation patterns. This suggests that increasing contact repulsion does not simply melt crystalline order; it can redirect it into filamentary or lamellar supersolid organization.

## 5. Excitations, roton physics, and superfluid response

In free space for the quasi-2D SOC system, the soft-core integral is approximated by a local term
$$
\gamma'=\int d^2r\, \frac{\tilde V_6}{1+|\mathbf r|^6}=\frac{2\pi^2}{3^{3/2}}\tilde V_6,
$$
so that the effective intra- and inter-component couplings become $\alpha_1\equiv \beta_{jj}+\gamma'$ and $\alpha_2\equiv \beta_{12}+\gamma'$. For a plane wave with momentum $k\hat x$, one has
$$
\psi_j=\frac{1}{\sqrt2}e^{i(kx-\mu t)}(e^{i\theta_k},\pm1)^T,\qquad \tan\theta_k=\frac{k_x}{k_y},\qquad \mu=\frac{k^2}{2}-\kappa k+\frac{\alpha_1+\alpha_2}{2},
$$
and linearization with $\delta\psi_j=e^{i(kx-\mu t)}[u_j e^{-i\omega t}-v_j^*e^{i\omega t}]$ leads to a Bogoliubov–de Gennes eigenvalue problem. The stability criteria are explicit: dynamical instability occurs for $|k|<\kappa$, where $\omega$ becomes complex; Landau instability occurs for $|k|>k_c$, where the lowest $\omega<0$; and superfluidity survives only in the finite window $\kappa\le |k|\le k_c$. For strong $\tilde V_6$, the typical spectra have two branches that contact at one or two points [2501.06555].

In the three-dimensional mixture, Bogoliubov analysis of the uniform miscible state yields two branches,
$$
\omega_\pm(k)=k\sqrt{\frac{k^2}{4}+\rho\bigl(g+\mathcal U(k)\bigr)\pm \rho\bigl(g+\mathcal U^{(12)}(k)\bigr)},
$$
with $\mathcal U(k)$ and $\mathcal U^{(12)}(k)$ the Fourier transforms of the soft-core potentials. The roton minimum of the out-of-phase branch $\omega_-(k)$ fixes the emergent length scale $a_{\rm rot}=2\pi/k_{\rm rot}$ and correctly predicts the lattice constants of the SP, ST, and ionic crystals. Tracking the softening of $\omega_-(k)$ locates the uniform-to-crystal transition boundaries [2507.16277].

Several order parameters clarify the relation between crystalline order and superfluidity. In 2D, the structure factor and topological charge distinguish density ordering from spin topology; in 3D, Bragg peaks in $S_\alpha(\mathbf k)=|n_\alpha(\mathbf k)|^2$ identify FCC, SC, ST, and SP order. The zero-temperature superfluid fraction is defined via the Leggett criterion,
$$
f_s=\frac{M}{N\hbar^2}\frac{\partial^2E(\theta)}{\partial\theta^2}\Big|_{\theta=0},
$$
and all four 3D phases maintain a nonzero $f_s$. The first-order coherence function $g^{(1)}_\alpha(\mathbf r,\mathbf r')=\langle \Psi_\alpha^\dagger(\mathbf r)\Psi_\alpha(\mathbf r')\rangle$ decays slowly, described as power-law or algebraically, off the lattice sites. A common misconception is that pronounced density modulation alone implies loss of superfluid transport; the reported nonzero $f_s$ and slow decay of $g^{(1)}$ show that, within these mean-field analyses, the crystalline phases remain supersolid rather than purely solid.

## 6. Nonequilibrium formation, rotation, and time-domain phenomena

The three-dimensional study addresses dynamical accessibility by real-time Gross–Pitaevskii simulations starting from uniform $\Psi_1$ and $\Psi_2$ of equal densities with small random white noise. At $t=0$, the Rydberg-dressed interactions are quenched to final values of $R_c$, $\tilde C_6$, $R_c^{(12)}/R_c$, and $a_s$, and the dynamics are propagated to approximately $t\approx 17\,{\rm ms}$ for Cesium mass and blockade radii $R_c\approx 4.5\times 10^4\,a_0$. For $a_s=40\,a_0$ with $R_c^{(12)}/R_c=0.95$, the SC ionic crystal nucleates around $t\sim 4\,{\rm ms}$ and stabilizes by $t\sim 16\,{\rm ms}$; for $a_s=80\,a_0$, ST filaments appear and coarsen; and for $a_s=120\,a_0$, SP layers form. The instantaneous density-correlation function $g(r)=\langle n(0)n(r)\rangle/\langle n\rangle^2$ develops sharp peaks on the emergent lattice vectors, confirming long-range order [2507.16277].

In the quasi-2D SOC system, rotation and dissipation generate a different nonequilibrium phenomenon. The time-dependent equations are modified to
$$
(i-\lambda)\partial_t\psi_1=\cdots-\Omega_rL_z\psi_1-\kappa(i\partial_x+\partial_y)\psi_2-\frac{\delta}{2}\psi_2,
$$
$$
(i-\lambda)\partial_t\psi_2=\cdots-\Omega_rL_z\psi_2-\kappa(i\partial_x-\partial_y)\psi_1-\frac{\delta}{2}\psi_1,
$$
with $\lambda=0.005$ and $\Omega_r\sim 0.6$. Starting from the chiral supersolid, one finds quasi-periodic, self-sustained oscillations in $\langle L_z\rangle(t)$ and in the real-space density radius. These persistent oscillations at fixed drive are identified as a dissipative continuous time crystal, with symmetry broken in time without external modulation. The same work also reports that, when rotation or an in-plane quadrupole magnetic field is included in the ground-state analysis, the chiral supersolid is broken and the ground state tends toward a miscible phase [2501.06555].

Taken together, the dynamical studies indicate that the ordered states are not only variational minima but also accessible outcomes of time evolution under experimentally motivated protocols. This suggests two distinct nonequilibrium roles for two-component Rydberg-dressed condensates: as pattern-forming media that self-organize into binary supersolids, and as driven-dissipative media that support persistent temporal order when initialized in a chiral supersolid configuration.

Source: https://www.emergentmind.com/topics/two-component-rydberg-dressed-bose-einstein-condensates