Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Nonlinear atom-field dynamics in high-Q cavities: from a BEC to a thermal gas

Open Access Open Access

Abstract

A cold gas of polarizable particles moving in the optical potential of a standing wave high finesse optical resonator acts as a dynamic refractive index. For a sufficiently strong cavity pump the optical forces generated by the intra cavity field perturb the particles phase space distribution, which shifts the optical resonance frequency and induces a nonlinear optical response. By help of the corresponding Vlasov equation we predict that beyond the known phenomenon of optical bi-stability one finds regions in parameter space, where no stable stationary solution exists. The atom field dynamics then exhibits oscillatory solutions converging to stable limit cycles of the system. The linearized analytical predictions agree well with corresponding numerical solutions of the full time dependent equations and first experimental observation in both cases.

©2011 Optical Society of America

1. Introduction

Based on the surprisingly fast progress in the experimental technology of laser cooling and manipulation of dilute atomic gases as well as optical resonator and laser stabilization technology, it is now experimentally possible to confine larger and larger numbers of colder and colder atoms within the mode volume of stabilized very high finesse resonators [1,2]. For slow enough particles the optical dipole force induced by the intracavity field significantly influences the atomic motion even in the dispersive limit at very large detunings, where absorption and spontaneous emission only play a minor role. In this limit the dispersive scattering of the particles simply acts as a dynamic refractive index changing the intracavity field evolution [35]. The resulting complex coupled atom-field dynamics leads to a wealth of interesting physical phenomena and applications [6, 7], which can be analysed theoretically using a wide range of models of very different complexity.

In previous work we have developed a new approach based on a classical phase space density description of the particles involving a corresponding Vlasov equation together with classical equations for the field mode amplitudes [8]. This approach is particularly suitable to describe very large particle numbers dispersively coupled to a single or a few driven cavity modes, when spontaneous emission and direct interparticle interaction play only a minor role. Related approaches to cold atom dynamics have put forward also by other groups [9, 10].

It has been known theoretically and experimentally for some time now, that the particles in this case act as a nonlinear optical medium through their motional response to the field, even though they are only weakly excited in the regime of linear polarizability at this frequency and intensity [11, 12]. Even for the conceptually most simple case of a single mode standing wave resonator, the coupled atom-field system can exhibit optical bistability. For a BEC in such a cavity self pulsing solutions were experimentally found [13]. They can be qualitatively well described via a fairly simple two-mode expansion for the atomic mean field dynamics as a nonlinearly coupled oscillator model [14]. This model was shown to be mathematically equivalent to a typical optomechanical setup of a high finesse resonator with a single field mode and a movable mirror at one end. Here we study the appearance of a very similar behavior for a general thermal gas in resonator and show that we can treat both limits of very low and high temperatures in a unified mean field model.

The paper is organized as follows. After presenting our general model in section 2 we use solve the Vlasov equation to find stationary states in section 3, which can give several solutions for a single choice of parameters. Using a linear perturbation analysis we then check the stability regions of these solution in section 4. In section 5 we construct and analyze the long time limit cycle solutions in the unstable regime. These are confirmed by numerical solutions of the time dependent Vlasov equation [Eq. (64)].

2. Model and basic equations

 figure: Fig. 1

Fig. 1 Schematic drawing of the physical system. The particles are indicated as the shaded cloud in the center. The resonator mode decays at a rate κ and is driven by a laser characterised by the pump parameter η.

Download Full Size | PPT Slide | PDF

We consider a large number N → ∞ polarizable particles of mass m > 0 inside a standing wave cavity interacting with a single mode of the resonator. The cavity photons induce dipole forces on the particles and the distribution of particles in turn shifts the cavity resonance through their dipole moments, thus giving rise to a nonlinearly coupled system. The dynamical quantities of the mean field type model that will be examined in this work are the phase space (quasi-) distribution f (x,v,t), which is the Wigner transform of the reduced one particle density matrix ρp ,1(x,x′,t),

f(x,v,t)=m2πħeizmv/ħρp,1(x+z2,xz2,t)dz,
and the non-dimensional complex cavity field amplitude, the average of the annihilation operator: α(t) = 〈â(t)〉. As detailed in the appendix, in the limit of large particle numbers, the phase space density obeys Wigner’s equation:
ft+vfx+U0|α|22sin(2kx)(f(x,v+vR)f(x,vvR))=0,
where vR:=ħkm and the mode amplitude satisfies
α˙=[κ+i(ΔcNU0/2)]αiNU02αdv0Lf(x,v,t)cos(2kx)dx+η.
In these equations, Δc = ωLωc is the detuning between the pump laser frequency ωL and the cavity resonance frequency ωc. η > 0 is related to the amplitude of the pump laser and will be referred to as pump strength or pump parameter. Moreover, U 0 quantifies the polarisability of the particles, L > 0 denotes the cavity length and k > 0 the wave vector of the corresponding resonator mode. The latter decays at a rate κ. To simplify the problem, let us adopt periodic boundary conditions f (0,v,t) = f (L,v,t) and assume that the optical forces exerted on the particles are weak. Since in this case the spatial modulation of the particle distribution will be small, we may write
f(x,v,t)=f0(v)+δf(x,v,t),
where f 0 denotes the homogeneous bulk distribution and δ f designates the small deviation therefrom. Note however that due to there being many particles present, this modulation may nevertheless be accompanied by a considerable shift in the cavity resonance. These assumptions imply that δ f approximately satisfies the linearized Wigner equation
(t+vx)δf=U0|α|22sin(2kx)(f0(v+vR)f0(vvR)).
Defining ϕ(v,t):=0Lδf(x,v,t)e2ikxdx and F 0(v) := L f 0(v), we are left with
ϕt+2ikvϕ=iU0|α|24(F0(v+vR)F0(vvR)),
α˙=[κ+i(ΔcNU0/2)]αiNU02αRe{ϕ(v,t)}dv+η.
Note that for every fixed velocity v, ϕ (v,.) satisfies the differential equation of a driven harmonic oscillator with a natural frequency given by 2k|v|.

3. Steady states

Let’s turn to the determination of the steady states (ϕ 0,α 0) of Eqs. (5)(6). One easily finds

ϕ0(v)=vRU0|α0|241kvΔF0Δv,
where we have defined
ΔF0Δv:=F0(v+vR)F0(vvR)2vR.
Substituting Eq. (7) into Eq. (6) we get
0=(κ+iδ)α0ivRNU028k|α0|2α01vΔF0Δvdv+η,
where we have set δ := ΔcNU 0/2. Whenever the temperature is nonzero, we can define a function F(x), such that
F(x):=vTF0(xvT),
or equivalently F0(v)=vT1F(v/vT). To shorten the notation, let us introduce the factor
J(vR/vT):=121xF(x+vR/vT)F(xvR/vT)vR/vTdx,
and let us define the steady state photon number I 0 = |α 0|2. From Eq. (9) one gets a cubic equation for the steady state photon number
I0[κ2+(δI0ωRNU02J4ωT2)2]=η2.
Here, ωR := h̄k2/2m denotes the recoil frequency and ωT = kvT. Clearly, Eq. (12) can have up to three distinct real solutions. A typical response curve is depicted in Fig. (2). In the interval of effective detuning designated by the letter A the system allows for several possible steady states at once. Which one of these (if any) can be attained? A general investigation of the linear stability properties of steady states will shed some light on that question. Let us therefore establish the necessary and sufficient conditions for a given steady state solution (ϕ 0, I 0) to be stable or unstable.

 figure: Fig. 2

Fig. 2 Normalized steady state photon number versus effective detuning δ/κ for N = 105, U 0 = 0.04κ, η =18κ, ωT = κ ωR = 0.5 · 10−3 κ and a Gaussian velocity distribution.

Download Full Size | PPT Slide | PDF

4. Stability analysis

4.1. Small Perturbation Analysis

The starting point of our analysis will be Eqs. (5)-(6), linearized around a steady state solution. Hence let us write

ϕ(v,t)=ϕ0(v)+δϕ(v,t),
α(t)=α0+δα(t).
Substituting into the equations of motion and keeping only terms that are linear in the deviations δϕ and δα, we obtain
δϕt+2ikvδϕ=ivRU02(α0δα*+α0*δα)ΔF0Δv,
δα˙=(κ+iΔ)δαiNU02α0Re{δϕ(v,t)}dv,
where we have set Δ:=δI0ωRNU02J4ωT2. It is advantageous to introduce the Laplace transforms with respect to time of δϕ, δϕ *, δα and δα *. To keep the notation to a minimum, we will not introduce new symbols for the transformed quantities. Introducing the “susceptibility”
χ(s):=NU02vR8(ΔF0/Δvs+2ikvΔF0/Δvs2ikv)dv,
we obtain the linear system of equations:
(s+κiΔI0χ(s)α02χ(s)α0*2χ(s)s+κ+iΔ+I0χ(s))·(δαδα*)=(I1I2),
where (I 1, I 2) denotes an initial value contribution. Let D(s) designate the determinant of the matrix on the left hand side of this last equation, i.e.,
D(s)=(κ+s)2+Δ22iΔI0χ(s).
We will henceforth call the complex valued function D(s) dispersion relation. As explained e.g. in [15], instability is equivalent to the existence of zeros of the dispersion relation with positive real part.

4.2. Limit of zero temperature-BEC

Before we direct our attention to the limiting case of a classical gas, let us now look closer at the opposite limit of a condensate at zero temperature. In that case, the velocity distribution is given by

F0(v)=δ(v),
where δ (v) is the dirac delta function. The steady-state photon number can be calculated from Eq. (12) with the replacements: J → − 1 and ωT → 2ωR. In order to make contact with the method of analysis used in other works, we decompose the solution to Eq. (5) according to
ϕ(v,t)=ϕ+(t)δ(vvR)+ϕ(t)δ(v+vR)+ϕ0(v,t),
where ϕ 0 is a solution of the homogeneous equation. This last contribution can be neglected for sufficiently smooth initial conditions due to rapid “phase mixing”, i.e.,
ϕ0(v,0)e2ikvtdv0,ast.
Defining x(t) := Re {ϕ + (t) +ϕ (t) } one easily checks that Eq. (5) leads to
x¨+(4ωR)2x=ωRU0|α|2,
α˙=[k+i(ΔcNU0/2)]αiNU02αx+η,
which are precisely the equations used to study the dynamics of a BEC inside a standing wave cavity using a Gross-Pitaevskii equation approach [16]. The stability of steady states is governed by Eq. (19), in which the susceptibility takes on the form
χ(s)=iωRNU0221s2+(4ωR)2.
Without going into the detailed consequences of this susceptibility, we note that it can be demonstrated that no steady state with Δ=δ+NU02I032ωR>0 is stable. This also means that a positive effective detuning, δ > 0, implies instability.

4.3. Thermal gas - classical limit

Let us now investigate the dispersion relation for a nonzero temperature. Note that as it stands, Eq. (19) together with Eq. (38) are defined only for Re(s) > 0. To find the boundary that separates unstable from stable states, we need to find the limiting form of the dispersion relation as Re(s) → 0+. To this end we recall the Plemjel formula

limɛ0+1xy±iɛ=Pxyiπδ(xy).
Here and elsewhere, P denotes the cauchy principal value. Setting s = γ + and writing sloppily χ(ω) = limγ →0+ χ(s), we find
χ(ω)=NU02vR8k[1iPvΔF0/Δv(ω/2k2)v2dv+π2(ΔF0Δv|ω2kΔF0Δv|ω2k)].
It does not make much sense to try to extract general (i.e. true for arbitrary F 0) consequences of Eqs. (19) and (27), for even if we could analytically find the stability boundary in terms of a function I0crit=I0crit(δ,N,U 0,...), we would still have to find the intersection of this curve with the one specifying the actually possible combinations (I 0,δ,...) (i.e., Eq. (12)), and this is obviously a hopelessly complicated task. We will therefore concentrate on the physically most relevant velocity distribution, namely the thermal distribution. Well above the condensation threshold, it is given by the Maxwell-Boltzmann distribution given by
F0(v)=1vTπe(v/vT)2;J=2.
In this regime, we have vTvR, such that we may replace
ΔF0ΔvdF0dv
in all relevant equations, including Eq. (5). Finally, we should state the condition that the deviation of the particle distribution from the homogeneous be small, namely
ħU0|α|2mvT22,or4U0η2ωRκ2ωT21,
which simply means that the potential energy must be much smaller than the kinetic energy. We will not try to solve the stability problem fully analytically (which is impossible), but use graphical and numerical methods where necessary instead. Using Eq. (27) and writing limγ →0+ D(s) = Dr(ω) + iDi(ω), we find
Dr(ω)=κ2+Δ2ω2+NU02ΔI0ωRωT2(1ωωTDa(ω2ωT)),
Di(ω)=2κωπNU02ΔI0ωR2ωT2ωωTe(ω2ωT)2.
Here, Da(x) denotes Dawson’s integral which is defined as
Da(x)=ex20xet2dt.
Note that the dispersion relation belonging to a thermal velocity distribution can now be seen (upon substituting Δ → δ and U02I04η2) as being identical to the one considered in [8]. Therefore all the results exactly carry over to the present case: A steady state I 0 is unstable if and only if we have
Dr(ω0)>0
for all ω 0 such that
Di(ω0)=0.

4.4. Graphical solution of the stability problem

In order to extract the information included in the general criterion Eq. (34), we will resort to graphical means. To do so, we note that the two conditions Di = 0, Dr = 0 define a function

I0=I0crit(Δ)
(the dependence on the other parameters is understood implicitly), with the property that a steady state I 0 is unstable if and only if I0>I0crit. Hence (Δ, I0crit(Δ)) defines the stability boundary in the parameter-space spanned by (Δ,I 0). Using that
δ=ΔNU02ωR2ωT2I0,
we can parametrically plot the stability boundary, i.e. the mapping
Δ(ΔNU02ωR2ωT2I0crit(Δ),I0cirt(Δ))
in the space of interest, namely the one spanned by (δ,I 0). Together with the response curve defined by Eq. (12), this allows to analyze the stability properties of the system. Consider again the set of parameters used to plot Fig. (2). Figure (3) depicts the response curve and regions of instability defined by Eqs. (12) and (35) for two different values of the pump parameter. Clearly, there exists an interval of effective detunings, where two distinct steady photon numbers are stable (bistability). Perhaps more surprising however is the finding of an interval where there is no stable steady state at all (Fig. (3b))! To summarize the possible behavior of the system when varying the pump strength, we can state that below a certain value, the response curve is single-valued and lies entirely in the stable domain. Above that value, we find multi-valuedness and a region of bistability (Fig. (3a)). Only upon crossing a second threshold does the completely unstable interval appear (Fig. (3b)).

 figure: Fig. 3

Fig. 3 Response curve (red) and regions of instability (shaded) for η = 13κ (a) and η = 18κ (b). Those parts of the response curve that lie inside the unstable region (red dashed), correspond to linearly unstable steady states. The black dashed line separates the parameter space into positive and negative Δ. The intervalls designated A correspond to bistability, the intervall designated B supports no stable steady state at all. All other parameters as in Fig. (2).

Download Full Size | PPT Slide | PDF

5. Long-time behavior

5.1. Classification

Whenever a steady state (two steady states) is (are) available, we expect the system to relax to it (or one of them) starting from arbitrary initial conditions. How does the system behave for parameters that allow for no stable steady state? To answer this question and to see whether the first guess is correct, we discreetized and solved Eqs. (5)-(6) numerically. Figure (4)a depicts the time averaged photon number (divided by the maximally possible) for different effective detunings for the usual set of parameters. It is clearly visible that this averaged photon number deviates from the steady state photon number only for detunings within the unstable interval, confirming our calculations. To our surprise and for generic sets of parameters, the system exhibits limit cycle oscillations within the completely unstable interval. Unlike in the case of a collisionless BEC, where oscillatory solutions depend on the initial conditions, these are true limit cycles independent of the initial conditions. This difference can be attributed to the phenomenon of “phase mixing” of the infinite number of oscillators in the classical case. The “observed” frequencies of the latter versus detuning are shown in figure (4)b. For the chosen parameters, these frequencies are of the order of the thermal frequency. For weak instability, the limit cycle oscillation is almost monochromatic. Figure (5) depicts these oscillations of the photon number for two detunings inside the unstable interval. Associated to these oscillations is a standing density wave formed by the particles.

 figure: Fig. 4

Fig. 4 a: Response curve (stable parts red, unstable parts black dotted) and time averaged normalized photon number (blue circles) versus effective detuning. A deviation from the response curve occurs only inside the unstable intervall [−3.3,0], where the system exhibits limit cycle oscillations. b: Limit cycle frequency versus detuning. Parameters as usual.

Download Full Size | PPT Slide | PDF

 figure: Fig. 5

Fig. 5 (a) Photon number as a function of time for δ = −κ (upper row) and δ = −3κ (lower row) and (b) normalized power spectrum. (c) Shape of the limit cycle. All additional parameters as in Fig. (2).

Download Full Size | PPT Slide | PDF

In Fig. (6) we compare the numerical solution of the reduced model with the solution of the Vlasov equation [Eq. (64)] for the same parameters as in Fig. (5). The agreement between the two models is excellent, confirming the validity of the approximations we introduced.

 figure: Fig. 6

Fig. 6 Comparison between the reduced (blue) and the Vlasov model Eq. (64) (red). Apart from a phase shift, the results are almost identical. All parameters as in Fig. (2) and Fig. (5).

Download Full Size | PPT Slide | PDF

5.2. Limit cycles frequencies close to threshold

As mentioned, one finds numerically that for parameters such that the steady state with Δ > 0 is weakly unstable, in its place there exists a periodic oscillation (limit cycle), which is almost monochromatic. This enables us to construct these cycles in a perturbative manner. Without loss of generality, we may assume that the sought for solution satisfies

χ(t):=Re{ϕ(v,t)}dvχ0+χ1cos(ωt)
and
|α(t)|2=:I(t)=nIneinωt
with real coefficients χ 0 and χ 1 and In=In*. Substituting Eq. (38) into Eq. (6) and using eix cos( φ ) = Σn Jn(x)einφ where Jn are the Bessel functions of the first kind, we obtain the fourier coefficients of the photon number in terms of ω,χ 0, χ 1 as
In(χ0,χ1)=η2m=Jm(NU0χ12ω)Jmn(NU0χ12ω)[κ+iΔimω][κiΔi(nm)ω],
where Δ=δ+NU0χ02. In order to link χ 0,χ 1 to the photon number coefficients In self-consistently, we must evaluate the response of the particles to an applied intensity varying like ΣnIneinωt in the form ϕ = Σn ϕn(v)einωt, i.e. we have to solve
(ω+2kv)ϕn=vRU0In2F0v.
Note carefully that we cannot simply divide this equation by + 2kv to obtain the solution because of the singularity at v = –ω/2k. Instead, we have to regard ϕn(v) as a generalized function to obtain
ϕn(v)=vRU0In2(F0(v)vPnω+2kv+hn(v)δ(nω+2kv)),
where P denotes the cauchy principal value and hn(v) is arbitrary. As a remedy to this uniqueness problem, one may introduce collisions characterized by a collision frequency and let this frequency go to zero in the end. Alternatively, one may replace
einωtlimγ0+e(inω+γ)t,
which means to consider the intensity as turned on infinitely slowly from the infinitely remote past. In either case the arbitrariness is lifted and we get
hn(v)=iπF0(v)v.
From this solution we find
χ0=vRU0I02k1vF0(v)vdv=U0ωRJ2ωT2I0
and
χ1=U0vR4k(IaF0(v)vdvω/2k+vIbπF0v(ω/2k)),
where Ia=I1+I1* and Ib=i(I1I1*). As we have chosen 02π/ωχ(t)sin(ωt)dt=0, we have additionally
IbF0(v)vdvω/2k+v=πIaF0v(ω/2k).
By means of Eq. (40), Eqs. (45), (46), and (47) constitute three equations for the three unknowns ω,χ 0,χ 1 characterizing the limit cycle solution we wish to find. To simplify the problem at hand, let us assume that the cycle be such that
|NU0χ12ω|1.
Then we may use that for n0:Jn(x)1n!(x2)n, if |x|1+n and J n = (−1)n Jn. Then we have
χ0=χ0(0)+χ0(2)(ω)(NU0χ14ω)2+,
where χ0(0) solves
X0(0)=U0ωRJ2ωT2η2κ2+Δ02,
with Δ0=δ+NU0χ0(0)/2. This of course means that, to lowest order, the time averaged shift of the cavity resonance frequency in the cycle is the same as the one corresponding to the unstable steady state. With these results, one may write
I1=I1(1)(ω,χ0(0))NU0χ14ω+I1(3)(ω,χ0(0))(NU0χ14ω)3+,
where
I1(1)(ω,χ0(0))=2η2ωΔ0(κ2+Δ02)(κ2+Δ02ω2+2iκω).
Now, to get the cycle frequency ω, we solve Eq. (47) with Eq. (51) to lowest order in the small quantity. This yields the dispersion relation for limit cycles just above threshold as
2κωF0(v)vdvω/2k+v=(κ2+Δ02ω2)πF0v(ω/2k),
which, for a thermal distribution, is the equation determining the frequency of oscillations of perturbations of the unstable steady state. For ωωT we find
ω=κ2+Δ02+4κωTπ.
Note that as Δ0 (if positive) is monotonously increasing with increasing effective detuning δ, so is the limit cycle frequency.

 figure: Fig. 7

Fig. 7 Response curve and numerically determined limit cycle frequency (blue triangles) vs. theoretical frequency from Eq. (53) (green dashed). The blue circles show the time averaged photon number. Parameters: N = 105, kvT = κ, ωR = 10−3 κ and U 0 = 0.05κ.

Download Full Size | PPT Slide | PDF

6. Conclusions and outlook

The nonlinear dynamical response of a high-Q optical cavity filled with a cold gas of polarizable particles can be treated by an effective Vlasov type mean field approach over a large range of temperatures from close to zero (BEC) to a thermal gas. In both cases density waves of the medium are strongly coupled to the intra-cavity field dynamics. Besides the well known optical bi-stability effect concurrent with a hysteresis effect, we can clearly identify the regions in parameter space, where no stable stationary solutions exists and time evolution converges to a periodic limit cycle. Its frequency is characteristic for the temperature and gas properties and can be explicitly calculated in the BEC limit or close to threshold for a thermal gas. In both cases such oscillations have been reported experimentally and thoroughly studied [13, 17, 18]. The observations agree well with our model. Interestingly, the effect is predicted to persist at much higher temperatures creating a nonlinear response in a linear medium. To observe it in practice at high temperatures would require sufficiently high intensities and densities or very high Q cavities. Here collisions can no longer be ignored, which can obscure this effect or even lead to a qualitative similar behavior based on a quite different physical mechanism. The collisionless model is in general a valid approximation for laser cooled thermal gases but gets more doubtful in the BEC limit. In the current experiments on this system at ETHZ [7] or Tübingen their influence could largely be neglected so far. In the long run the system thus represents a well controllable and experimentally implementable toy system to study generic nonlinear dynamics in the transition range from classical to quantum mechanics. Adding extra modes or atomic species will allow to extend its scope and complexity in many directions.

7. Appendix: quantum and classical mean field limit

Here we show, how the mean field Wigner equation can be obtained from the quantum mechanical equations for N → ∞ polarizable particles coupled to one or more resonator modes. When spontaneous emission is neglected, the Heisenberg equation of motion for the particle field operator ψ^(x,t) reads

tψ^=iħ2m2x2ψ^+1iħV^ψ^,
wherein = (x, â 1,..., âM) = is the potential at x created by the photons of the M modes. Clearly we have [,ψ^] = 0. Let us introduce the phase space density operator
f^(x,p)=12πħeizp/ħψ^(xz2)ψ^(x+z2)dz,
where p = mv is the momentum. It may be noted that 〈 (x,p)〉 is the Wigner transform of the one body reduced density matrix ρp ,1(x,x′). It may be shown from Eq. (55) that the phase space density operator satisfies:
f^t+pmf^x+U^(x,q)f^(x,pħq)dq=0,
where
U^(x,q):=i2πħdz[V^(x+z/2)V(xz/2)]eiqz.
In case of a single mode of a standing wave resonator, we have (x,â) = h̄U 0 â â cos2(kx) and with this one readily finds from Eqs. (57) and (58) that the equation for the phase space density operator thus reads:
f^t+pmf^x+U02a^a^sin(2kx)[f^(x,p+ħk)f^(x,pħk)]=0.
This equation is to be supplemented by the Heisenberg-Langevin equation for the annihilation operator â of the mode:
a^˙=(κ+iδ)a^iU0a^0Ldxf^(x,p)cos(2kx)dp+η+a^in.
The mean field limit N ≫ 1 can be found by decomposing âj = αj + δâj and = f + δf̂ with 〈δâj〉 = 〈δf̂〉 = 0, taking the average of Eq. (57) and neglecting all correlations. This procedure leads to Wigner’s equation (sometimes called the Quantum Vlasov equation) for the expectation value of (x, p)
ft+pmfx+U02|α|2sin(2kx)[f(x,p+ħk)f(x,pħk)]=0,
α˙=(κ+iδ)αiU0α20Ldxf(x,p)cos(2kx)dp+η.
If we further assume that the typical momentum is much larger than h̄k, we can develop
f(x,p±ħk)f(x,p)±ħkfp(x,p)+ħ2k222fp2(x,p).
Then Eq. (61) becomes
ft+pmfx+ħkU0|α|2sin(2kx)fp=0,
which is Vlasov’s equation. Hence, the classical model employed in [8] and in this work is expected to be accurate if kBT 2 k 2/2m, i.e. if the thermal energy (or kinetic energy per particle) is much higher than the recoil energy.

Acknowledgments

This work was supported by the Austrian Science Fund FWF under grant Nr. I119-N16 and SFB F40 FOQUS. We thank Claus Zimmerman for stimulating discussions and Matteo Cristiani and Jürgen Eschner for helpful comments and early communication of their experimental results.

References and links

1. D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003). [CrossRef]  

2. M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008). [CrossRef]  

3. P. Domokos and H. Ritsch, “Mechanical effects of light in optical resonators,” J. Opt. Soc. Am. B 20, 1098–1130 (2003). [CrossRef]  

4. P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997). [CrossRef]  

5. P. Domokos, P. Horak, and H. Ritsch, “Semiclassical theory of cavity-assisted atom cooling,” J. Phys. B 34, 187–198 (2001). [CrossRef]  

6. J. M. Zhang, F. C. Cui, D. L. Zhou, and W. M. Liu, “Nonlinear dynamics of a cigar-shaped bose-einstein condensate in an optical cavity,” Phys. Rev. A 79, 033401 (2009). [CrossRef]  

7. S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009). [CrossRef]  

8. T. Grießer, H. Ritsch, M. Hemmerling, and G. Robb, “A Vlasov approach to bunching and selfordering of particles in optical resonators,” Eur. Phys. J. D 58, 349–368 (2010). [CrossRef]  

9. J. Javaloyes, M. Perrin, G. L. Lippi, and A. Politi, “Self-generated cooperative light emission induced by atomic recoil,” Phys. Rev. A 70, 023405 (2004). [CrossRef]  

10. R. Bach, K. Burnett, M. d’Arcy, and S. Gardiner, “Quantum-mechanical cumulant dynamics near stable periodic orbits in phase space: application to the classical-like dynamics of quantum accelerator modes,” Phys. Rev. A 71, 33417 (2005). [CrossRef]  

11. S. Gupta, K. L. Moore, K. W. Murch, and D. M. Stamper-Kurn, “Cavity nonlinear optics at low photon numbers from collective atomic motion,” Phys. Rev. Lett. 99, 213601 (2007). [CrossRef]  

12. A. Vukics, W. Niedenzu, and H. Ritsch, “Cavity nonlinear optics with few photons and ultracold quantum particles,” Phys. Rev. A 79, 013828 (2009). [CrossRef]  

13. F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose–Einstein condensate,” Science 322, 235–238 (2008). [CrossRef]   [PubMed]  

14. D. Nagy, P. Domokos, A. Vukics, and H. Ritsch, “Nonlinear quantum dynamics of two BEC modes dispersively coupled by an optical cavity,” Eur. Phys. J. D 55, 659–668 (2009). [CrossRef]  

15. D. C. Montgomery, Theory of the unmagnetized plasma (Gordon & Breach, 1971).

16. F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007). [CrossRef]   [PubMed]  

17. M. Cristiani and J. Eschner (personal communication, 2011).

18. T. Valenzuela, M. Cristiani, H. Gothe, and J. Eschner, “Cold Ytterbium atoms in high-finesse optical cavities: cavity cooling and collective interactions,” in “Lasers and Electro-Optics 2009 and the European Quantum Electronics Conference. CLEO Europe-EQEC 2009. European Conference on,” (IEEE, 2009), p. 1. [CrossRef]  

References

  • View by:

  1. D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
    [Crossref]
  2. M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
    [Crossref]
  3. P. Domokos and H. Ritsch, “Mechanical effects of light in optical resonators,” J. Opt. Soc. Am. B 20, 1098–1130 (2003).
    [Crossref]
  4. P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997).
    [Crossref]
  5. P. Domokos, P. Horak, and H. Ritsch, “Semiclassical theory of cavity-assisted atom cooling,” J. Phys. B 34, 187–198 (2001).
    [Crossref]
  6. J. M. Zhang, F. C. Cui, D. L. Zhou, and W. M. Liu, “Nonlinear dynamics of a cigar-shaped bose-einstein condensate in an optical cavity,” Phys. Rev. A 79, 033401 (2009).
    [Crossref]
  7. S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
    [Crossref]
  8. T. Grießer, H. Ritsch, M. Hemmerling, and G. Robb, “A Vlasov approach to bunching and selfordering of particles in optical resonators,” Eur. Phys. J. D 58, 349–368 (2010).
    [Crossref]
  9. J. Javaloyes, M. Perrin, G. L. Lippi, and A. Politi, “Self-generated cooperative light emission induced by atomic recoil,” Phys. Rev. A 70, 023405 (2004).
    [Crossref]
  10. R. Bach, K. Burnett, M. d’Arcy, and S. Gardiner, “Quantum-mechanical cumulant dynamics near stable periodic orbits in phase space: application to the classical-like dynamics of quantum accelerator modes,” Phys. Rev. A 71, 33417 (2005).
    [Crossref]
  11. S. Gupta, K. L. Moore, K. W. Murch, and D. M. Stamper-Kurn, “Cavity nonlinear optics at low photon numbers from collective atomic motion,” Phys. Rev. Lett. 99, 213601 (2007).
    [Crossref]
  12. A. Vukics, W. Niedenzu, and H. Ritsch, “Cavity nonlinear optics with few photons and ultracold quantum particles,” Phys. Rev. A 79, 013828 (2009).
    [Crossref]
  13. F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose–Einstein condensate,” Science 322, 235–238 (2008).
    [Crossref] [PubMed]
  14. D. Nagy, P. Domokos, A. Vukics, and H. Ritsch, “Nonlinear quantum dynamics of two BEC modes dispersively coupled by an optical cavity,” Eur. Phys. J. D 55, 659–668 (2009).
    [Crossref]
  15. D. C. Montgomery, Theory of the unmagnetized plasma (Gordon & Breach, 1971).
  16. F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
    [Crossref] [PubMed]
  17. M. Cristiani and J. Eschner (personal communication, 2011).
  18. T. Valenzuela, M. Cristiani, H. Gothe, and J. Eschner, “Cold Ytterbium atoms in high-finesse optical cavities: cavity cooling and collective interactions,” in “Lasers and Electro-Optics 2009 and the European Quantum Electronics Conference. CLEO Europe-EQEC 2009. European Conference on,” (IEEE, 2009), p. 1.
    [Crossref]

2010 (1)

T. Grießer, H. Ritsch, M. Hemmerling, and G. Robb, “A Vlasov approach to bunching and selfordering of particles in optical resonators,” Eur. Phys. J. D 58, 349–368 (2010).
[Crossref]

2009 (4)

J. M. Zhang, F. C. Cui, D. L. Zhou, and W. M. Liu, “Nonlinear dynamics of a cigar-shaped bose-einstein condensate in an optical cavity,” Phys. Rev. A 79, 033401 (2009).
[Crossref]

S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
[Crossref]

A. Vukics, W. Niedenzu, and H. Ritsch, “Cavity nonlinear optics with few photons and ultracold quantum particles,” Phys. Rev. A 79, 013828 (2009).
[Crossref]

D. Nagy, P. Domokos, A. Vukics, and H. Ritsch, “Nonlinear quantum dynamics of two BEC modes dispersively coupled by an optical cavity,” Eur. Phys. J. D 55, 659–668 (2009).
[Crossref]

2008 (2)

F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose–Einstein condensate,” Science 322, 235–238 (2008).
[Crossref] [PubMed]

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

2007 (2)

S. Gupta, K. L. Moore, K. W. Murch, and D. M. Stamper-Kurn, “Cavity nonlinear optics at low photon numbers from collective atomic motion,” Phys. Rev. Lett. 99, 213601 (2007).
[Crossref]

F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
[Crossref] [PubMed]

2005 (1)

R. Bach, K. Burnett, M. d’Arcy, and S. Gardiner, “Quantum-mechanical cumulant dynamics near stable periodic orbits in phase space: application to the classical-like dynamics of quantum accelerator modes,” Phys. Rev. A 71, 33417 (2005).
[Crossref]

2004 (1)

J. Javaloyes, M. Perrin, G. L. Lippi, and A. Politi, “Self-generated cooperative light emission induced by atomic recoil,” Phys. Rev. A 70, 023405 (2004).
[Crossref]

2003 (2)

P. Domokos and H. Ritsch, “Mechanical effects of light in optical resonators,” J. Opt. Soc. Am. B 20, 1098–1130 (2003).
[Crossref]

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

2001 (1)

P. Domokos, P. Horak, and H. Ritsch, “Semiclassical theory of cavity-assisted atom cooling,” J. Phys. B 34, 187–198 (2001).
[Crossref]

1997 (1)

P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997).
[Crossref]

Alt, W.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Bach, R.

R. Bach, K. Burnett, M. d’Arcy, and S. Gardiner, “Quantum-mechanical cumulant dynamics near stable periodic orbits in phase space: application to the classical-like dynamics of quantum accelerator modes,” Phys. Rev. A 71, 33417 (2005).
[Crossref]

Baumann, K.

S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
[Crossref]

Benhelm, J.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Bourdel, T.

F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
[Crossref] [PubMed]

Brennecke, F.

S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
[Crossref]

F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose–Einstein condensate,” Science 322, 235–238 (2008).
[Crossref] [PubMed]

F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
[Crossref] [PubMed]

Burnett, K.

R. Bach, K. Burnett, M. d’Arcy, and S. Gardiner, “Quantum-mechanical cumulant dynamics near stable periodic orbits in phase space: application to the classical-like dynamics of quantum accelerator modes,” Phys. Rev. A 71, 33417 (2005).
[Crossref]

Courteille, P. W.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Cristiani, M.

M. Cristiani and J. Eschner (personal communication, 2011).

T. Valenzuela, M. Cristiani, H. Gothe, and J. Eschner, “Cold Ytterbium atoms in high-finesse optical cavities: cavity cooling and collective interactions,” in “Lasers and Electro-Optics 2009 and the European Quantum Electronics Conference. CLEO Europe-EQEC 2009. European Conference on,” (IEEE, 2009), p. 1.
[Crossref]

Cui, F. C.

J. M. Zhang, F. C. Cui, D. L. Zhou, and W. M. Liu, “Nonlinear dynamics of a cigar-shaped bose-einstein condensate in an optical cavity,” Phys. Rev. A 79, 033401 (2009).
[Crossref]

d’Arcy, M.

R. Bach, K. Burnett, M. d’Arcy, and S. Gardiner, “Quantum-mechanical cumulant dynamics near stable periodic orbits in phase space: application to the classical-like dynamics of quantum accelerator modes,” Phys. Rev. A 71, 33417 (2005).
[Crossref]

Domokos, P.

D. Nagy, P. Domokos, A. Vukics, and H. Ritsch, “Nonlinear quantum dynamics of two BEC modes dispersively coupled by an optical cavity,” Eur. Phys. J. D 55, 659–668 (2009).
[Crossref]

P. Domokos and H. Ritsch, “Mechanical effects of light in optical resonators,” J. Opt. Soc. Am. B 20, 1098–1130 (2003).
[Crossref]

P. Domokos, P. Horak, and H. Ritsch, “Semiclassical theory of cavity-assisted atom cooling,” J. Phys. B 34, 187–198 (2001).
[Crossref]

Donner, T.

S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
[Crossref]

F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose–Einstein condensate,” Science 322, 235–238 (2008).
[Crossref] [PubMed]

F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
[Crossref] [PubMed]

Dotsenko, I.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Elsässer, T.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Eschner, J.

T. Valenzuela, M. Cristiani, H. Gothe, and J. Eschner, “Cold Ytterbium atoms in high-finesse optical cavities: cavity cooling and collective interactions,” in “Lasers and Electro-Optics 2009 and the European Quantum Electronics Conference. CLEO Europe-EQEC 2009. European Conference on,” (IEEE, 2009), p. 1.
[Crossref]

M. Cristiani and J. Eschner (personal communication, 2011).

Esslinger, T.

S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
[Crossref]

F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose–Einstein condensate,” Science 322, 235–238 (2008).
[Crossref] [PubMed]

F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
[Crossref] [PubMed]

Gardiner, S.

R. Bach, K. Burnett, M. d’Arcy, and S. Gardiner, “Quantum-mechanical cumulant dynamics near stable periodic orbits in phase space: application to the classical-like dynamics of quantum accelerator modes,” Phys. Rev. A 71, 33417 (2005).
[Crossref]

Gheri, K. M.

P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997).
[Crossref]

Gothe, H.

T. Valenzuela, M. Cristiani, H. Gothe, and J. Eschner, “Cold Ytterbium atoms in high-finesse optical cavities: cavity cooling and collective interactions,” in “Lasers and Electro-Optics 2009 and the European Quantum Electronics Conference. CLEO Europe-EQEC 2009. European Conference on,” (IEEE, 2009), p. 1.
[Crossref]

Grießer, T.

T. Grießer, H. Ritsch, M. Hemmerling, and G. Robb, “A Vlasov approach to bunching and selfordering of particles in optical resonators,” Eur. Phys. J. D 58, 349–368 (2010).
[Crossref]

Guerlin, C.

S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
[Crossref]

Gupta, S.

S. Gupta, K. L. Moore, K. W. Murch, and D. M. Stamper-Kurn, “Cavity nonlinear optics at low photon numbers from collective atomic motion,” Phys. Rev. Lett. 99, 213601 (2007).
[Crossref]

Hechenblaikner, G.

P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997).
[Crossref]

Hemmerich, A.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Hemmerling, M.

T. Grießer, H. Ritsch, M. Hemmerling, and G. Robb, “A Vlasov approach to bunching and selfordering of particles in optical resonators,” Eur. Phys. J. D 58, 349–368 (2010).
[Crossref]

Horak, P.

P. Domokos, P. Horak, and H. Ritsch, “Semiclassical theory of cavity-assisted atom cooling,” J. Phys. B 34, 187–198 (2001).
[Crossref]

P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997).
[Crossref]

Javaloyes, J.

J. Javaloyes, M. Perrin, G. L. Lippi, and A. Politi, “Self-generated cooperative light emission induced by atomic recoil,” Phys. Rev. A 70, 023405 (2004).
[Crossref]

Kampschulte, T.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Khudaverdyan, M.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Köhl, M.

F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
[Crossref] [PubMed]

Kruse, D.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Lenhard, K.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Lippi, G. L.

J. Javaloyes, M. Perrin, G. L. Lippi, and A. Politi, “Self-generated cooperative light emission induced by atomic recoil,” Phys. Rev. A 70, 023405 (2004).
[Crossref]

Liu, W. M.

J. M. Zhang, F. C. Cui, D. L. Zhou, and W. M. Liu, “Nonlinear dynamics of a cigar-shaped bose-einstein condensate in an optical cavity,” Phys. Rev. A 79, 033401 (2009).
[Crossref]

Meschede, D.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Montgomery, D. C.

D. C. Montgomery, Theory of the unmagnetized plasma (Gordon & Breach, 1971).

Moore, K. L.

S. Gupta, K. L. Moore, K. W. Murch, and D. M. Stamper-Kurn, “Cavity nonlinear optics at low photon numbers from collective atomic motion,” Phys. Rev. Lett. 99, 213601 (2007).
[Crossref]

Murch, K. W.

S. Gupta, K. L. Moore, K. W. Murch, and D. M. Stamper-Kurn, “Cavity nonlinear optics at low photon numbers from collective atomic motion,” Phys. Rev. Lett. 99, 213601 (2007).
[Crossref]

Nagorny, B.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Nagy, D.

D. Nagy, P. Domokos, A. Vukics, and H. Ritsch, “Nonlinear quantum dynamics of two BEC modes dispersively coupled by an optical cavity,” Eur. Phys. J. D 55, 659–668 (2009).
[Crossref]

Niedenzu, W.

A. Vukics, W. Niedenzu, and H. Ritsch, “Cavity nonlinear optics with few photons and ultracold quantum particles,” Phys. Rev. A 79, 013828 (2009).
[Crossref]

Perrin, M.

J. Javaloyes, M. Perrin, G. L. Lippi, and A. Politi, “Self-generated cooperative light emission induced by atomic recoil,” Phys. Rev. A 70, 023405 (2004).
[Crossref]

Politi, A.

J. Javaloyes, M. Perrin, G. L. Lippi, and A. Politi, “Self-generated cooperative light emission induced by atomic recoil,” Phys. Rev. A 70, 023405 (2004).
[Crossref]

Rauschenbeutel, A.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Reick, S.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Ritsch, H.

T. Grießer, H. Ritsch, M. Hemmerling, and G. Robb, “A Vlasov approach to bunching and selfordering of particles in optical resonators,” Eur. Phys. J. D 58, 349–368 (2010).
[Crossref]

A. Vukics, W. Niedenzu, and H. Ritsch, “Cavity nonlinear optics with few photons and ultracold quantum particles,” Phys. Rev. A 79, 013828 (2009).
[Crossref]

D. Nagy, P. Domokos, A. Vukics, and H. Ritsch, “Nonlinear quantum dynamics of two BEC modes dispersively coupled by an optical cavity,” Eur. Phys. J. D 55, 659–668 (2009).
[Crossref]

P. Domokos and H. Ritsch, “Mechanical effects of light in optical resonators,” J. Opt. Soc. Am. B 20, 1098–1130 (2003).
[Crossref]

P. Domokos, P. Horak, and H. Ritsch, “Semiclassical theory of cavity-assisted atom cooling,” J. Phys. B 34, 187–198 (2001).
[Crossref]

P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997).
[Crossref]

Ritter, S.

S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
[Crossref]

F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose–Einstein condensate,” Science 322, 235–238 (2008).
[Crossref] [PubMed]

F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
[Crossref] [PubMed]

Robb, G.

T. Grießer, H. Ritsch, M. Hemmerling, and G. Robb, “A Vlasov approach to bunching and selfordering of particles in optical resonators,” Eur. Phys. J. D 58, 349–368 (2010).
[Crossref]

Ruder, M.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Schörner, K.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Stamper-Kurn, D. M.

S. Gupta, K. L. Moore, K. W. Murch, and D. M. Stamper-Kurn, “Cavity nonlinear optics at low photon numbers from collective atomic motion,” Phys. Rev. Lett. 99, 213601 (2007).
[Crossref]

Stecher, H.

P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997).
[Crossref]

Valenzuela, T.

T. Valenzuela, M. Cristiani, H. Gothe, and J. Eschner, “Cold Ytterbium atoms in high-finesse optical cavities: cavity cooling and collective interactions,” in “Lasers and Electro-Optics 2009 and the European Quantum Electronics Conference. CLEO Europe-EQEC 2009. European Conference on,” (IEEE, 2009), p. 1.
[Crossref]

von Cube, C.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Vukics, A.

A. Vukics, W. Niedenzu, and H. Ritsch, “Cavity nonlinear optics with few photons and ultracold quantum particles,” Phys. Rev. A 79, 013828 (2009).
[Crossref]

D. Nagy, P. Domokos, A. Vukics, and H. Ritsch, “Nonlinear quantum dynamics of two BEC modes dispersively coupled by an optical cavity,” Eur. Phys. J. D 55, 659–668 (2009).
[Crossref]

Widera, A.

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Zhang, J. M.

J. M. Zhang, F. C. Cui, D. L. Zhou, and W. M. Liu, “Nonlinear dynamics of a cigar-shaped bose-einstein condensate in an optical cavity,” Phys. Rev. A 79, 033401 (2009).
[Crossref]

Zhou, D. L.

J. M. Zhang, F. C. Cui, D. L. Zhou, and W. M. Liu, “Nonlinear dynamics of a cigar-shaped bose-einstein condensate in an optical cavity,” Phys. Rev. A 79, 033401 (2009).
[Crossref]

Zimmermann, C.

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

Appl. Phys. B (1)

S. Ritter, F. Brennecke, K. Baumann, T. Donner, C. Guerlin, and T. Esslinger, “Dynamical coupling between a Bose–Einstein condensate and a cavity optical lattice,” Appl. Phys. B 95, 213–218 (2009).
[Crossref]

Eur. Phys. J. D (2)

T. Grießer, H. Ritsch, M. Hemmerling, and G. Robb, “A Vlasov approach to bunching and selfordering of particles in optical resonators,” Eur. Phys. J. D 58, 349–368 (2010).
[Crossref]

D. Nagy, P. Domokos, A. Vukics, and H. Ritsch, “Nonlinear quantum dynamics of two BEC modes dispersively coupled by an optical cavity,” Eur. Phys. J. D 55, 659–668 (2009).
[Crossref]

J. Opt. Soc. Am. B (1)

J. Phys. B (1)

P. Domokos, P. Horak, and H. Ritsch, “Semiclassical theory of cavity-assisted atom cooling,” J. Phys. B 34, 187–198 (2001).
[Crossref]

Nature (1)

F. Brennecke, T. Donner, S. Ritter, T. Bourdel, M. Köhl, and T. Esslinger, “Cavity QED with a Bose–Einstein condensate,” Nature 450, 268–271 (2007).
[Crossref] [PubMed]

New J. Phys. (1)

M. Khudaverdyan, W. Alt, I. Dotsenko, T. Kampschulte, K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schörner, A. Widera, and D. Meschede, “Controlled insertion and retrieval of atoms coupled to a high-finesse optical resonator,” New J. Phys. 10, 073023 (2008).
[Crossref]

Phys. Rev. A (5)

D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, P. W. Courteille, T. Elsässer, B. Nagorny, and A. Hemmerich, “Cold atoms in a high-Q ring cavity,” Phys. Rev. A 67, 051802 (2003).
[Crossref]

J. Javaloyes, M. Perrin, G. L. Lippi, and A. Politi, “Self-generated cooperative light emission induced by atomic recoil,” Phys. Rev. A 70, 023405 (2004).
[Crossref]

R. Bach, K. Burnett, M. d’Arcy, and S. Gardiner, “Quantum-mechanical cumulant dynamics near stable periodic orbits in phase space: application to the classical-like dynamics of quantum accelerator modes,” Phys. Rev. A 71, 33417 (2005).
[Crossref]

J. M. Zhang, F. C. Cui, D. L. Zhou, and W. M. Liu, “Nonlinear dynamics of a cigar-shaped bose-einstein condensate in an optical cavity,” Phys. Rev. A 79, 033401 (2009).
[Crossref]

A. Vukics, W. Niedenzu, and H. Ritsch, “Cavity nonlinear optics with few photons and ultracold quantum particles,” Phys. Rev. A 79, 013828 (2009).
[Crossref]

Phys. Rev. Lett. (2)

S. Gupta, K. L. Moore, K. W. Murch, and D. M. Stamper-Kurn, “Cavity nonlinear optics at low photon numbers from collective atomic motion,” Phys. Rev. Lett. 99, 213601 (2007).
[Crossref]

P. Horak, G. Hechenblaikner, K. M. Gheri, H. Stecher, and H. Ritsch, “Cavity-induced atom cooling in the strong coupling regime,” Phys. Rev. Lett. 79, 4974–4977 (1997).
[Crossref]

Science (1)

F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, “Cavity optomechanics with a Bose–Einstein condensate,” Science 322, 235–238 (2008).
[Crossref] [PubMed]

Other (3)

M. Cristiani and J. Eschner (personal communication, 2011).

T. Valenzuela, M. Cristiani, H. Gothe, and J. Eschner, “Cold Ytterbium atoms in high-finesse optical cavities: cavity cooling and collective interactions,” in “Lasers and Electro-Optics 2009 and the European Quantum Electronics Conference. CLEO Europe-EQEC 2009. European Conference on,” (IEEE, 2009), p. 1.
[Crossref]

D. C. Montgomery, Theory of the unmagnetized plasma (Gordon & Breach, 1971).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1
Fig. 1 Schematic drawing of the physical system. The particles are indicated as the shaded cloud in the center. The resonator mode decays at a rate κ and is driven by a laser characterised by the pump parameter η.
Fig. 2
Fig. 2 Normalized steady state photon number versus effective detuning δ/κ for N = 105, U 0 = 0.04κ, η =18κ, ωT = κ ωR = 0.5 · 10−3 κ and a Gaussian velocity distribution.
Fig. 3
Fig. 3 Response curve (red) and regions of instability (shaded) for η = 13κ (a) and η = 18κ (b). Those parts of the response curve that lie inside the unstable region (red dashed), correspond to linearly unstable steady states. The black dashed line separates the parameter space into positive and negative Δ. The intervalls designated A correspond to bistability, the intervall designated B supports no stable steady state at all. All other parameters as in Fig. (2).
Fig. 4
Fig. 4 a: Response curve (stable parts red, unstable parts black dotted) and time averaged normalized photon number (blue circles) versus effective detuning. A deviation from the response curve occurs only inside the unstable intervall [−3.3,0], where the system exhibits limit cycle oscillations. b: Limit cycle frequency versus detuning. Parameters as usual.
Fig. 5
Fig. 5 (a) Photon number as a function of time for δ = −κ (upper row) and δ = −3κ (lower row) and (b) normalized power spectrum. (c) Shape of the limit cycle. All additional parameters as in Fig. (2).
Fig. 6
Fig. 6 Comparison between the reduced (blue) and the Vlasov model Eq. (64) (red). Apart from a phase shift, the results are almost identical. All parameters as in Fig. (2) and Fig. (5).
Fig. 7
Fig. 7 Response curve and numerically determined limit cycle frequency (blue triangles) vs. theoretical frequency from Eq. (53) (green dashed). The blue circles show the time averaged photon number. Parameters: N = 105, kvT = κ, ωR = 10−3 κ and U 0 = 0.05κ.

Equations (66)

Equations on this page are rendered with MathJax. Learn more.

f ( x , v , t ) = m 2 π ħ e i z m v / ħ ρ p , 1 ( x + z 2 , x z 2 , t ) d z ,
f t + v f x + U 0 | α | 2 2 sin ( 2 k x ) ( f ( x , v + v R ) f ( x , v v R ) ) = 0 ,
α ˙ = [ κ + i ( Δ c NU 0 / 2 ) ] α i NU 0 2 α d v 0 L f ( x , v , t ) cos ( 2 k x ) d x + η .
f ( x , v , t ) = f 0 ( v ) + δ f ( x , v , t ) ,
( t + v x ) δ f = U 0 | α | 2 2 sin ( 2 k x ) ( f 0 ( v + v R ) f 0 ( v v R ) ) .
ϕ t + 2 i k v ϕ = i U 0 | α | 2 4 ( F 0 ( v + v R ) F 0 ( v v R ) ) ,
α ˙ = [ κ + i ( Δ c NU 0 / 2 ) ] α i NU 0 2 α Re { ϕ ( v , t ) } d v + η .
ϕ 0 ( v ) = v R U 0 | α 0 | 2 4 1 k v Δ F 0 Δ v ,
Δ F 0 Δ v : = F 0 ( v + v R ) F 0 ( v v R ) 2 v R .
0 = ( κ + i δ ) α 0 i v R NU 0 2 8 k | α 0 | 2 α 0 1 v Δ F 0 Δ v d v + η ,
F ( x ) : = v T F 0 ( x v T ) ,
J ( v R / v T ) : = 1 2 1 x F ( x + v R / v T ) F ( x v R / v T ) v R / v T d x ,
I 0 [ κ 2 + ( δ I 0 ω R NU 0 2 J 4 ω T 2 ) 2 ] = η 2 .
ϕ ( v , t ) = ϕ 0 ( v ) + δ ϕ ( v , t ) ,
α ( t ) = α 0 + δ α ( t ) .
δ ϕ t + 2 i k v δ ϕ = i v R U 0 2 ( α 0 δ α * + α 0 * δ α ) Δ F 0 Δ v ,
δ α ˙ = ( κ + i Δ ) δ α i NU 0 2 α 0 Re { δ ϕ ( v , t ) } d v ,
χ ( s ) : = NU 0 2 v R 8 ( Δ F 0 / Δ v s + 2 i k v Δ F 0 / Δ v s 2 i k v ) d v ,
( s + κ i Δ I 0 χ ( s ) α 0 2 χ ( s ) α 0 * 2 χ ( s ) s + κ + i Δ + I 0 χ ( s ) ) · ( δ α δ α * ) = ( I 1 I 2 ) ,
D ( s ) = ( κ + s ) 2 + Δ 2 2 i Δ I 0 χ ( s ) .
F 0 ( v ) = δ ( v ) ,
ϕ ( v , t ) = ϕ + ( t ) δ ( v v R ) + ϕ ( t ) δ ( v + v R ) + ϕ 0 ( v , t ) ,
ϕ 0 ( v , 0 ) e 2 i k v t d v 0 , as t .
x ¨ + ( 4 ω R ) 2 x = ω R U 0 | α | 2 ,
α ˙ = [ k + i ( Δ c NU 0 / 2 ) ] α i NU 0 2 α x + η ,
χ ( s ) = i ω R NU 0 2 2 1 s 2 + ( 4 ω R ) 2 .
lim ɛ 0 + 1 x y ± i ɛ = P x y i π δ ( x y ) .
χ ( ω ) = NU 0 2 v R 8 k [ 1 i P v Δ F 0 / Δ v ( ω / 2 k 2 ) v 2 d v + π 2 ( Δ F 0 Δ v | ω 2 k Δ F 0 Δ v | ω 2 k ) ] .
F 0 ( v ) = 1 v T π e ( v / v T ) 2 ; J = 2.
Δ F 0 Δ v d F 0 d v
ħ U 0 | α | 2 m v T 2 2 , or 4 U 0 η 2 ω R κ 2 ω T 2 1 ,
D r ( ω ) = κ 2 + Δ 2 ω 2 + NU 0 2 Δ I 0 ω R ω T 2 ( 1 ω ω T Da ( ω 2 ω T ) ) ,
D i ( ω ) = 2 κ ω π NU 0 2 Δ I 0 ω R 2 ω T 2 ω ω T e ( ω 2 ω T ) 2 .
Da ( x ) = e x 2 0 x e t 2 d t .
D r ( ω 0 ) > 0
D i ( ω 0 ) = 0.
I 0 = I 0 crit ( Δ )
δ = Δ NU 0 2 ω R 2 ω T 2 I 0 ,
Δ ( Δ NU 0 2 ω R 2 ω T 2 I 0 crit ( Δ ) , I 0 cirt ( Δ ) )
χ ( t ) : = Re { ϕ ( v , t ) } d v χ 0 + χ 1 cos ( ω t )
| α ( t ) | 2 = : I ( t ) = n I n e i n ω t
I n ( χ 0 , χ 1 ) = η 2 m = J m ( NU 0 χ 1 2 ω ) J m n ( NU 0 χ 1 2 ω ) [ κ + i Δ i m ω ] [ κ i Δ i ( n m ) ω ] ,
( ω + 2 k v ) ϕ n = v R U 0 I n 2 F 0 v .
ϕ n ( v ) = v R U 0 I n 2 ( F 0 ( v ) v P n ω + 2 k v + h n ( v ) δ ( n ω + 2 k v ) ) ,
e i n ω t lim γ 0 + e ( i n ω + γ ) t ,
h n ( v ) = i π F 0 ( v ) v .
χ 0 = v R U 0 I 0 2 k 1 v F 0 ( v ) v d v = U 0 ω R J 2 ω T 2 I 0
χ 1 = U 0 v R 4 k ( I a F 0 ( v ) v d v ω / 2 k + v I b π F 0 v ( ω / 2 k ) ) ,
I b F 0 ( v ) v d v ω / 2 k + v = π I a F 0 v ( ω / 2 k ) .
| NU 0 χ 1 2 ω | 1.
χ 0 = χ 0 ( 0 ) + χ 0 ( 2 ) ( ω ) ( NU 0 χ 1 4 ω ) 2 + ,
X 0 ( 0 ) = U 0 ω R J 2 ω T 2 η 2 κ 2 + Δ 0 2 ,
I 1 = I 1 ( 1 ) ( ω , χ 0 ( 0 ) ) NU 0 χ 1 4 ω + I 1 ( 3 ) ( ω , χ 0 ( 0 ) ) ( NU 0 χ 1 4 ω ) 3 + ,
I 1 ( 1 ) ( ω , χ 0 ( 0 ) ) = 2 η 2 ω Δ 0 ( κ 2 + Δ 0 2 ) ( κ 2 + Δ 0 2 ω 2 + 2 i κ ω ) .
2 κ ω F 0 ( v ) v d v ω / 2 k + v = ( κ 2 + Δ 0 2 ω 2 ) π F 0 v ( ω / 2 k ) ,
ω = κ 2 + Δ 0 2 + 4 κ ω T π .
t ψ ^ = i ħ 2 m 2 x 2 ψ ^ + 1 i ħ V ^ ψ ^ ,
f ^ ( x , p ) = 1 2 π ħ e i z p / ħ ψ ^ ( x z 2 ) ψ ^ ( x + z 2 ) d z ,
f ^ t + p m f ^ x + U ^ ( x , q ) f ^ ( x , p ħ q ) d q = 0 ,
U ^ ( x , q ) : = i 2 π ħ d z [ V ^ ( x + z / 2 ) V ( x z / 2 ) ] e i q z .
f ^ t + p m f ^ x + U 0 2 a ^ a ^ sin ( 2 k x ) [ f ^ ( x , p + ħ k ) f ^ ( x , p ħ k ) ] = 0.
a ^ ˙ = ( κ + i δ ) a ^ i U 0 a ^ 0 L d x f ^ ( x , p ) cos ( 2 k x ) d p + η + a ^ i n .
f t + p m f x + U 0 2 | α | 2 sin ( 2 k x ) [ f ( x , p + ħ k ) f ( x , p ħ k ) ] = 0 ,
α ˙ = ( κ + i δ ) α i U 0 α 2 0 L d x f ( x , p ) cos ( 2 k x ) d p + η .
f ( x , p ± ħ k ) f ( x , p ) ± ħ k f p ( x , p ) + ħ 2 k 2 2 2 f p 2 ( x , p ) .
f t + p m f x + ħ k U 0 | α | 2 sin ( 2 k x ) f p = 0 ,

Metrics

Select as filters


Select Topics Cancel
© Copyright 2022 | Optica Publishing Group. All Rights Reserved