SciPost Submission Page
Mott transition in a cavityboson system: A quantitative comparison between theory and experiment
by Rui Lin, Christoph Georges, Jens Klinder, Paolo Molignini, Miriam Büttner, Axel U. J. Lode, R. Chitra, Andreas Hemmerich, Hans Keßler
 Published as SciPost Phys. 11, 030 (2021)
Submission summary
As Contributors:  Hans Keßler · Rui Lin · Axel U. J. Lode · Paolo Molignini 
Preprint link:  scipost_202104_00032v2 
Code repository:  http://ultracold.org/ 
Date accepted:  20210720 
Date submitted:  20210630 19:37 
Submitted by:  Lin, Rui 
Submitted to:  SciPost Physics 
Academic field:  Physics 
Specialties: 

Approaches:  Experimental, Computational 
Abstract
The competition between shortrange and cavitymediated infiniterange interactions in a cavityboson system leads to the existence of a superfluid phase and a Mottinsulator phase within the selforganized regime. In this work, we quantitatively compare the steadystate phase boundaries of this transition measured in experiments and simulated using the Multiconfigurational TimeDependent Hartree Method for Indistinguishable Particles. To make the problem computationally feasible, we represent the full system by the exact manybody wave function of a twodimensional fourwell potential. We argue that the validity of this representation comes from the nature of both the cavityatomic system and the BoseHubbard physics. Additionally we show that the chosen representation only induces small systematic errors, and that the experimentally measured and theoretically predicted phase boundaries agree reasonably. We thus demonstrate a new approach for the quantitative numerical determination of the superfluidMottinsulator phase boundary.
Published as SciPost Phys. 11, 030 (2021)
Author comments upon resubmission
We would also like to thank the three Referees for their careful examination of our work. Referee 1 and 2 recommended our manuscript for publication in SciPost Physics. We have revised the manuscript and implemented their comments to improve the manuscript. Referee 3 wrote “It could there be published if the authors address my concerns”. We hope that we address in the following his/her concerns satisfactorily.
Particularly, the three Referees have pointed out two main weaknesses of our manuscript: the systematic errors from the experiments are not discussed in a detailed manner; and a more substantial comparison between our numerical methods with existing ones is missing. These two points are both more properly addressed in the new version.
In the new manuscript, the major changes related to the Referee’s comments are marked in blue. Please note that Figs. 4 and 5 in the previous version now become Figs. 5 and 6, respectively.

Reply to comments from Referee 1.

(Retardation effects) Indeed the retardation effect of the cavity field would also affect the SSFSMI boundary. We now discuss it in more detail in Sec. 4.1. As we further discuss in the same section, this retardation effect is inevitable due to the dilemma between the retardation effect and the atom loss.
2 & 3. (Monochromaticity) In experiments, the main loss of the laser pump strength comes from the two side peaks due to servo electronics, which are in the order of MHz away from the central peak. This shift is much larger than the cavity dissipation rate, and thus these two side peaks cannot contribute to the scattering of the cavity field. On the other hand, the linewidth of the central peak cannot be accurately measured because we do not have another stable laser to beat with. However, we estimate that the linewidth is roughly 100Hz, which is much smaller than the cavity dissipation rate. It can thus fully scatter into the cavity. The calibration factor 1.36 seems reasonable to us. We add a few lines discussing this in Sec. 2.2.

(Background signal of the cavity field). This large background counts originates from diffuse scattering on the optics implemented in the system, and not from the cavity. It can thus be safely subtracted from the data for the determination of the phase boundary. We add a sentence in the caption of Fig. 6 discussing that. Moreover, as one can see in our recent preprint Ref. 54, the symmetry in our system is very well established and there is no significant breaking of the symmetry due to diffusively scattered light. We clarify this in Sec. 2.1.

(Comparison to existing approaches). We thank the Referee for his/her appreciation on this discussion. Inspired by this comment, we now further perform a more quantitative comparison between the MCTDHX approach and the existing approach based on Wannier function analysis. This is now discussed in the new Appendix F. In the main text, we also move the relevant paragraph to Sec. 4.2. We decide to keep the discussion after the main result instead of moving it to an earlier section, because it requires understanding of the SSFSMI transition, which is mainly discussed in Sec. 4.1.

Reply to comments from Referee 2.
1, 3 & 4. (Heating and atom loss). Indeed there are multiple effects induced by heating and atom loss which we have not fully discussed in the previous version.
 (Comment 3): The NPSSF boundary is slightly shifted to higher pump strengths for the large detunings due to heating and atom loss because we reach the selfordered phase at larger pump strength and hence later times. Therefore, when we are fitting this boundary, we only consider the data points from low detunings (2pi5 kHz to 2pi20 kHz). This has already been implemented for the fitting in the last version, and is now explicitly clarified in Sec. 3.2.
 (Comments 1 & 3): Moreover, the heating and atom loss are also accountable for the plateau seen in $\alpha^2$ in Fig. 6(b). This will certainly affect the cavityinduced optical potential. However, for the fast ramp presented in Fig. 6(b), this plateau occurs later than the Mott transition, so it has little effect on the predicted boundary. This is discussed now in Sec. 4.1.
 (Comment 4): However, inside the SMI phase, indeed the atom loss will change the shape of the central Bragg peak. This indeed partly contributes to the difference between the behavior of W in experiments and simulations. We now also discuss this in Sec. 4.1.
 (Meanfield density distributions in momentum space). In the meanfield limit, if we blindly use the algorithm to calculate the “momentumspace density distribution”, we would obtain the Fourier transform of the realspace density distribution. Nevertheless, we would say that this result does not have physical significance, because it has not converged numerically. Therefore, we don’t think that it would be useful to present them.
5 & 6. (Terminology and typos) We thank the Referee for pointing out these issues. We have now replaced “superlattice” with “lattice”, and fixed the typos.
 Reply to comments from Referee 3.
In the report, the Referee posts two questions regarding to the applicability of the fourwell model and the relevance of the beyondGrossPitaevskii meanfield effects, respectively.

(Applicability of the fourwell model). Indeed, the symmetry of the checkerboard lattice makes us possible to choose a simple model with only four lattice sites, and this is not guaranteed for many other systems. As the number of lattice sites in a unit cell increases, the required computational resources scale up significantly. We now include an extra paragraph in Sec. 3.3 to discuss this.

(Beyond meanfield effects). If we constrain ourselves only in the GrossPitaevskii meanfield limit, we are able to reproduce the selforganization, but we cannot capture the Mott insulation. This is because the Mott insulation requires at least one singleparticle orbital for each lattice site. If we blindly do the simulations in the meanfield limit even in the fourwell model, we will consistently obtain “superfluid” states even for very large lattice depth. This has already been discussed in detail in Sec. 3.1, but is now also briefly discussed in the introduction as suggested by the Referee.
The Referee also requests us to clarify two statements.

(Quantumoptical BoseHubbard model). Here what we have in mind is a BoseHubbard model with longrange interaction between lattice sites. We have made the statement more clear now.

(Parameters). Indeed our sentence was a bit confusing in the previous version. We want to say that these are the last two parameters which we need to determine for the simulation. We have clarified this now. If we want to use other approaches instead, for example the Wannierbased BoseHubbard approach as presented in Sec. 4.2 now, the same parameters are also required.
List of changes
Sec. 1: We clarify the phrase “quantumoptic version of the BoseHubbard model”, and briefly discuss the beyondmeanfield effects in our system.
Sec. 2.2: We further discuss in detail the retardation and heating effects on both the NPSSF and SSFSMI boundaries. We also clarify the source of power loss of the pump laser.
Sec. 3.2: We now include a new figure showing the real space density distribution of the full 2D model. We also clarify that the fit for the experimental NPSSF boundary only takes the data points from small detunings into account, in order to minimize the heating effects.
Sec. 3.3: We discuss how the computational effort scales when the symmetry of the optical lattice becomes weaker.
Sec. 4.1: We now include new panels (sx) in Fig.5 showing the real space density distribution of the full 2D model. This is for comparison with the fourwell model in panels (mr). We also discuss now the retardation and heating effects observed in Fig. 6 in more detail. We have also fixed a typo in the experimental results in Fig. 6(b). The number of photons is in the order of 10^3.
Sec. 4.2 & Appendix F: We perform a comparison of our MCTDHX approach with an alternative approach based on Wannier functions and BoseHubbard model. The detail of the comparison is shown in the new Appendix F.
We also rewrite several paragraphs to improve the narrative.
Submission & Refereeing History
You are currently on this page
Reports on this Submission
Anonymous Report 1 on 2021712 (Invited Report)
Report
I would like to thank the authors for addressing my previous comments in detail and to perform related changes on the manuscript. In particular the in depth discussion of potential systematic errors, as well as, the comparison to methods based on Wannier states following the other referee's comments provide valuable additional information to the reader.
Based on this and the arguments mentioned in my previous report I recommend publication of the manuscript in the SciPost physics.