Arxiv:2106.03883V2 [Cond-Mat.Quant-Gas] 21 Jun 2021
Total Page:16
File Type:pdf, Size:1020Kb
Cooling and state preparation in an optical lattice via Markovian feedback control Ling-Na Wu∗ and Andr´eEckardty Institut f¨urTheoretische Physik, Technische Universit¨atBerlin, Hardenbergstraße 36, Berlin 10623, Germany (Dated: 22nd June 2021) We propose and investigate a scheme for in-situ cooling a system of interacting bosonic atoms in a one-dimensional optical lattice without particle loss. It is based on Markovian feedback control. The system is assumed to be probed weakly via the homodyne detection of photons that are scattered off-resonantly by the atoms from a structured probe beam into a cavity mode. By applying an inertial force to the system that is proportional to the measured signal, the system can be guided into a pure target state. Cooling is achieved by preparing the system's ground state, but also excited states can be prepared. While the approach always works for weak interactions, for integer filling the ground state can be prepared with high fidelity for arbitrary interaction strengths. The scheme is found to be robust against reduced measurement efficiencies. Introduction.| Atomic quantum gases in optical lat- kHz scale is sufficient, which can be achieved easily us- tices constitute a unique experimental platform for ing digital signal processors. The Markovian feedback studying quantum many-body systems under extremely method has been applied to various control problems, in- clean and flexible conditions [1]. An important property cluding the stabilization of arbitrary one-qubit quantum of atomic quantum gases is their extremely high degree states [43, 44], the manipulation of quantum entangle- of isolation from the environment, which is provided by ment between two qubits [45{48] as well as optical and the optical trapping of atoms inside ultrahigh vacuum. spin squeezing [41, 49, 50]. It offers the rather unique opportunity to experiment- We consider a system of interacting bosonic atoms in ally study coherent nonequilibrium dynamics of many- a one-dimensional optical lattice. The system shall be body systems [2] over hundreds of tunneling times, before weakly measured via off-resonant scattering of photons eventually residual heating processes make themselves from an incoming structured probe beam, as recently felt. However, this isolation from the environment also proposed by Mehkov et al. [51, 52]. In order to enhance comes with a cost. Namely, whenever excitations are cre- the detection efficiency the system shall be placed inside ated, e.g. by time-dependent parameter variations, the an optical cavity [36, 53{56], so that (via the Purcell ef- induced energy remains in the system and causes heat- fect [57]) the photons are scattered predominantly into ing. Therefore, the ability to cool the system (without one cavity mode. We will show that for appropriately losing particles) in situ or even to mimic its coupling to a shaped probe beams, the system can then be cooled in low-temperature bath with tunable properties would be situ, by forcing it via a lattice acceleration that is pro- highly desirable. portional to the signal obtained from a homodyne detec- tion of the scattered photons leaking out of the cavity. Here we propose and investigate a scheme for in-situ We first investigate the non-interacting case and propose cooling optical lattice systems by means of measurement- a feedback scheme to stabilize its ground state. This based feedback control [3, 4], where depending on scheme is found to perform also well in the presence the outcome of a continuous measurement, the system of weak interactions. Moreover, for strong interactions, Hamiltonian is manipulated in a coherent fashion. As a we propose a method for cooling in the Mott phase (at valuable resource, measurement-based feedback control integer fillings) [58]. Numerical results show that the has been proposed and used for quantum state prepara- scheme performs very well even in the intermediate in- tion [5{17], cooling [18{27], simulating nonlinear dynam- teractions regimes and continuously connects to the non- ics [28{30], as well as the control of dynamics of a sys- interacting scheme so that for systems at integer filling tem [31{36] and phase transition [37{41], just to name efficient cooling is possible for the whole range of interac- a few. There are various feedback control strategies. In tion strengths. Moreover our approach is robust against this paper, we consider Markovian feedback control as a reduction of the measurement efficiency to about 50%. proposed by Wiseman [42], where a feedback Hamilto- Model.| We describe a system of N bosons with open arXiv:2106.03883v2 [cond-mat.quant-gas] 21 Jun 2021 nian proportional to the measurement signal is continu- boundary conditions using the one-dimensional Bose- ously added to the system. Such a feedback is assumed Hubbard model, to be instantaneous, namely, the delay time between the measurement and the application of the control field M−1 M X U X is small compared to the typical timescales of the sys- H = −J (aya + ay a ) + n (n − 1); (1) l l+1 l+1 l 2 l l tem. This assumption ensures the Markovianity of the l=1 l=1 dynamical description and thus makes it possible to y write down a feedback-modified Lindblad master equa- where al annihilates a particle on site l and nl = a al P l tion (ME) [42]. For our scheme, a control on the higher counts the particle number on site l, with l nl = N. 2 The first term in (1) describes tunneling between neigh- be neglected, since it is the only dissipative term. This boring sites with rate J, and the second term denotes strategy is confirmed in our simulations below, where Hfb on-site interactions with strength U. is fully taken into account. We are interested in cooling the system to its ground Feedback scheme for non-interacting case.| In order state. For this purpose, let us consider a homodyne to get an insight of how to choose the measurement and measurement of an operator c. The dynamical evolu- feedback operator, let us first take a look at the double- tion of the system is then described by the stochastic well system (M = 2). It can be mapped to a spin-N=2 master equation (SME) [3] (¯h = 1 hereafter), dρc = system [60] with the collective spin operators defined as −i[H; ρ ]dt + D[c]ρ dt + H[c]ρ dW , with H[c]ρ := cρ + 1 y y i y y c c c Jx = 2 (a1a2 + a2a1), Jy = − 2 (a1a2 − a2a1), and Jz = ρcy − Tr[(c + cy)ρ]ρ and D[c]ρ := cρcy − 1 (cycρ + ρcyc). 1 y y 2 (a1a1 − a2a2), which are angular momentum operators Here ρ denotes the quantum state conditioned on the 2 c satisfying the commutation relations [Jµ;Jν ] = iµνσJσ measurement result, I = Tr[(c + cy)ρ] + ξ(t), with hom for µ, ν; σ = x; y; z and the Levi-Civita symbol µνσ. The ξ(t) = dW=dt and dW being the standard Wiener in- ground state of the non-interacting Hamiltonian H2 = crement with mean zero and variance dt. The quantum y y −J(a1a2 + a2a1) = −2JJx is the eigenstate of Jx, jj = backaction of a weak measurement can be used for tailor- N=2; m = N=2i, with the largest eigenvalue m = N=2. ing the system's dynamics and to prepare target states. From the theory of angular momentum, we know that For instance, a quantum nondemolition measurement this state is also an eigenstate of operator J− = Jz − iJy of light allows for the preparation of different types of with eigenvalue 0. Hence, by choosing the measurement atom-number squeezed and macroscopic superposition p p operator c as γJz and feedback operator F as γJy, states [51, 52]. While the state generated in this way is p so that A = γ(Jz − iJy), we will reach a steady state conditional due to the nondeterministic nature of meas- very close to the ground state of the Hamiltonian H2 for urement, the introduction of feedback using the informa- weak measurement with strength γ J. tion acquired from the measurements allows to steer the Now let us turn to the multi-site case. As a generaliza- system's dynamics into a desired state. tion of the double-well system, we consider the following Here we consider a direct feedback strategy, where a measurement operator signal-dependent, i.e. conditional, feedback term IhomF M is added to the Hamiltonian. According to Markovian p X feedback control theory [42], the system is then governed c = γ zlnl: (5) by the feedback-modified SME l=1 It is a weighted sum of site occupations, which can be dρc = −i[H + Hfb; ρc]dt + D[A]ρcdt + H[A]ρcdW; (2) implemented via homodyne detection of the off-resonant with operators scattering of probe light from the atoms [52, 55, 61]. The coefficients z will be fixed later. For the feedback oper- 1 l A = c − iF; H = (cyF + F c): (3) ator, we take fb 2 M−1 p X By taking the ensemble average of the possible meas- F = γ(−i) (aya − ay a ); (6) urement outcomes, we arrive at the feedback-modified l l+1 l+1 l l=l ME [42] which describes tunneling between neighboring sites with ρ_ = −i[H + Hfb; ρ] + D[A]ρ. (4) complex rate and can be realized simply by an iner- tial force induced by accelerating the lattice, as will The effect induced by the feedback loop is seen to replace be discussed in detail below. For the non-interacting the collapse operator c by A and to add an extra term case (U = 0), the ground state jGi of (1) is given H to the Hamiltonian.