
Thermodynamic Phase Diagram of Two-Dimensional Bosons in a Quasicrystal
Potential
Zhaoxuan Zhu,1Hepeng Yao,2and Laurent Sanchez-Palencia1
1CPHT, CNRS, Ecole Polytechnique, IP Paris, F-91128 Palaiseau, France
2Department of Quantum Matter Physics, University of Geneva,
24 Quai Ernest-Ansermet, CH-1211 Geneva, Switzerland
(Dated: July 13, 2023)
Quantum simulation of quasicrystals in synthetic bosonic matter now paves the way to the explo-
ration of these intriguing systems in wide parameter ranges. Yet thermal fluctuations in such systems
compete with quantum coherence, and significantly affect the zero-temperature quantum phases.
Here we determine the thermodynamic phase diagram of interacting bosons in a two-dimensional,
homogeneous quasicrystal potential. Our results are found using quantum Monte Carlo simulations.
Finite-size effects are carefully taken into account and the quantum phases are systematically dis-
tinguished from thermal phases. In particular, we demonstrate stabilization of a genuine Bose glass
phase against the normal fluid in sizable parameter ranges. Our results for strong interactions are
interpreted using a fermionization picture and experimental relevance is discussed.
The discovery of quasiperiodic structures in plane
tilings [1] and material science [2,3] has profoundly al-
tered our dichotomous perception of order and disorder.
Lying at the interface of the two realms, quasicrystals
display a number of intriguing properties, including un-
usual localization and fractal properties, anomalous crit-
ical scalings, and phasonic degrees of freedom [4–9]. So
far, quasicrystals have been observed in their natural
state in meteorites [10,11] and nuclear blast residues [12]
or in the laboratory after fast solidification of certain al-
loys [2,13], and have been extensively studied in solid-
state physics [2,5,6,14–16]. Moreover, artificial qua-
sicrystals can now be engineered in synthetic quantum
matter with unique control knobs, using photonic crys-
tals [8,17–19], quantum fluids of light [20–22], and ultra-
cold quantum gases [23–25]. In the latter, defectless and
phononfree quasicrystal potentials can be emulated in
a variety of configurations using appropriately-arranged
sets of laser beams [26–31]. Furthermore, two-body in-
teractions can be tuned using magnetic control [32–35],
hence paving the way to the exploration of quantum
phase diagrams in wide parameter ranges.
In past years, one-dimensional (1D) quasiperiodic
models of ultracold atoms have been discussed quite ex-
haustively [36–54] but exploration of their 2D counter-
parts has only recently gained momentum, mostly in
tight-binding models [55–57]. So far, theoretical and
experimental work has demonstrated the emergence of
quasicrystalline order through matterwave interferome-
try [28,30], Anderson-like localization [28,31,58], and
Bose glass (BG) physics [31,55,56,59]. The BG is an
emblematic compressible insulator, characteristic of dis-
ordered or quasi-disordered systems and distinct from the
superfluid (SF) and Mott insulator (MI) phases, which
also appear in periodic systems [60–62]. In bosonic mod-
els, however, thermal fluctuations compete with (quasi-
)disorder, which has so far hindered the observation of
the BG phase [43,44]. It has been recently proposed that
this issue may be overcome by scaling up characteristic
energies using shallow quasiperiodic potentials [47]. Up
to now, this has been investigated only in 1D [48] and 2D
harmonically trapped [63] systems. In contrast, the case
of a 2D Bose gas with genuine long-range quasicrystal
order remains unexplored. Moreover, the central issue of
discriminating the BG phase from trivial thermal phases
has been hardly addressed. As argued below, this cannot
be achieved as in 1D and requires specific analysis in 2D.
In this Letter, we determine the first thermodynamic
phase diagrams of weakly to strongly interacting 2D Bose
gases in a shallow quasicrystal potential at finite temper-
atures. Quantum Monte Carlo simulations are performed
in quasicrystal, homogeneous potentials and finite-size ef-
fects are carefully taken into account. The SF, MI, and
BG quantum phases, induced by the competition of in-
teractions and quasicrystal potential, are systematically
discriminated from the normal fluid (NF), which is in-
stead dominated by thermal fluctuations. Most impor-
tantly, we find that the BG phase survives up to signif-
icantly high temperatures. Our results in the strongly-
interacting regime are interpreted using a fermionization
picture and implications to experiments in ultracold atom
systems are discussed.
Model.— The dynamics of the 2D Bose gas is gov-
erned by the Hamiltonian
ˆ
H=ˆdrΨ(r)†−ℏ2∇2
2m+V(r)Ψ(r),(1)
+1
2ˆdrdr′Ψ(r)†Ψ(r′)†U(r−r′)Ψ(r′)Ψ(r),
where Ψ(r)is the bosonic field operator at position rand
mis the particle mass. The quasicrystal potential,
V(r) = V0
4
X
k=1
cos2(Gk·r),(2)
arXiv:2210.15526v2 [cond-mat.quant-gas] 11 Jul 2023