Simulating dense QCD matter with ultracold atomic boson-fermion mixtures

We delineate, as an analog of two-flavor dense quark matter, the phase structure of a many-body mixture of atomic bosons and fermions in two internal states with a tunable boson-fermion attraction. The bosons b correspond to diquarks, and the fermions f to unpaired quarks. For weak b-f attraction, the system is a mixture of a Bose-Einstein condensate and degenerate fermions, while for strong attraction composite b-f fermions N, analogs of the nucleon, are formed, which are superfluid due to the N-N attraction in the spin-singlet channel. We determine the symmetry breaking patterns at finite temperature as a function of the b-f coupling strength, and relate the phase diagram to that of dense QCD.

💡 Research Summary

The paper proposes a concrete quantum‑simulation scheme for dense two‑flavor quark matter using an ultracold atomic mixture of bosons and fermions with two internal spin states. In the mapping, the bosonic atoms (b) play the role of diquarks, while the fermionic atoms (f) represent unpaired quarks. By exploiting a magnetic‑field‑tuned Feshbach resonance, the boson‑fermion attraction g_{bf} can be varied continuously, allowing the system to explore regimes that correspond to weakly coupled quark matter, a diquark‑condensed phase, and a nucleon‑like phase with superfluidity.

The authors first construct a low‑energy effective Hamiltonian that contains (i) a Gross‑Pitaevskii term for the bosons, (ii) a kinetic term for the two‑component fermions, and (iii) a contact interaction between b and f. Using a Hubbard‑Stratonovich transformation they introduce a composite fermion N ≡ bf, which mimics a nucleon (a bound state of a diquark and a quark). In the weak‑coupling limit (small |g_{bf}|) the bosons form a Bose‑Einstein condensate (BEC) while the fermions remain a degenerate Fermi gas; the symmetry breaking pattern is U(1)_b × SU(2)_f → 1, i.e. both the global phase symmetry of the condensate and the spin‑rotation symmetry of the fermions are spontaneously broken. Raising the temperature destroys the BEC, restoring the U(1) symmetry, while the fermionic sector stays normal.

When the attraction is increased beyond a critical value, the b‑f binding energy exceeds the Fermi energy and a stable N fermion emerges. The N particles inherit spin‑½ from the fermion and carry a residual internal “isospin” corresponding to the two fermionic spin states. Crucially, the N‑N interaction in the spin‑singlet channel is taken to be attractive, which leads to a BCS‑type pairing of N’s and a second superfluid transition. The resulting phase diagram in the (temperature, g_{bf}) plane therefore contains two continuous transition lines: (1) a BEC‑to‑N crossover where the bosonic condensate dissolves into bound composites, and (2) an N‑pairing line where the composite fermions undergo Cooper pairing. The critical temperature for the second transition follows the usual exponential dependence T_c ∝ n^{2/3} exp(−1/|g_{NN}|), with n the total density.

The symmetry‑breaking sequence in the strong‑coupling regime is thus U(1)_b → 1 (BEC disappears) followed by U(1)_N → 1 (N‑pairing), mirroring the hypothesized cascade in dense QCD: first color‑antitriplet diquark condensation (breaking of color gauge symmetry) and then color‑singlet nucleon superfluidity (analogous to nuclear superfluidity). The authors argue that the composite N plays the role of a nucleon, while the original boson b represents a diquark condensate; the two‑step symmetry breaking reproduces the essential features of the color‑superconducting phases predicted for quark matter at several times nuclear saturation density.

To make the proposal experimentally realistic, the paper discusses concrete atomic species such as ^87Rb (boson) and ^40K (fermion), for which wide Feshbach resonances have been characterized. Using typical densities (10^12–10^13 cm⁻³) and temperatures below 50 nK, the calculated binding energy of the bf molecule is on the order of a few kHz, sufficient to form a stable N fermion. The authors suggest detection strategies: (i) time‑of‑flight imaging to observe the disappearance of the BEC peak and the emergence of a fermionic momentum distribution characteristic of bound composites, (ii) radio‑frequency spectroscopy to measure the binding energy and the pairing gap of N‑N superfluidity, and (iii) Bragg scattering to probe collective modes associated with the two broken U(1) symmetries.

In the concluding section, the authors emphasize the broader significance of their work. By providing a tabletop platform that reproduces the hierarchy of symmetry breakings, pairing mechanisms, and phase transitions expected in dense QCD, the ultracold boson‑fermion mixture becomes a powerful analog quantum simulator. It opens the door to systematic studies of color‑superconductivity, nucleon superfluidity, and the equation of state of strongly interacting matter under controlled conditions. Future extensions could include adding a second fermionic species to simulate three‑flavor QCD, imposing optical lattices to explore crystalline color‑superconducting phases, or rotating the trap to generate vortices and study topological defects in the analog of quark‑matter vortices. The paper thus bridges high‑energy nuclear physics and low‑temperature atomic physics, offering a concrete pathway to experimentally probe phenomena that are otherwise inaccessible.