Advances in Optics, Photonics, Spectroscopy & Applications <strong>VI</strong> <strong>ISSN</strong> <strong>1859</strong> - <strong>4271</strong>Here the quasistationary conditions have been assumed and the notation δ Xˆ≡ Xˆ− 〈 Xˆ〉 hasbeen introduced.Fig. 5: Diffusion coefficient D as a function of (a) the atom-to-surface distance r − a and (b) theaxial position z for the parameters of Fig. 2, where ωp= ωc= ω0.We again consider the case where the probe field is so weak that the excited state is hardlyoccupied and there is at most one <stro<strong>ng</strong>>ph</stro<strong>ng</strong>>oton in the cavity. We plot the diffusion coefficient D inFig. 5 for the parameters of Fig. 2. Figure 5(a) shows that D has a peak when the atom-tosurfacedistance r − a is small enough but not too small. The appearance of such a peak is aresult of the balance between the compe<s<stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>>>tin</s<stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>>>g effects of the van der Waals potential and the atomcavitycoupli<strong>ng</strong>.V. CENTER-OF-MASS MOTION OF THE ATOMThe classical motion of the center of mass of the atom is governed by the equationsM & z= F z+ p&z(D)(M r& = Mr & ϕ + F r+ p&r,2 D),(D)Mr & ϕ= − 2Mr&& ϕ + F + p&.(16)ϕϕHere M is the mass of the atom, Fr, F ϕ , and Fzare the cylindrical-coordinate components(D) (D)(D)of the force F , and p&r, p&ϕ , and p&zare the random variables that characterize themomentum diffusion.We use the method of Ref. [6] to solve Eqs. (??). In every time step ∆ t of the numerical( D)integration of p , a stochastic term ∆p = rˆR ( 6 ) ˆrDr∆t+ ϕRϕ( 6Dϕ∆t) + zˆRz( 6Dz∆t) isadded. Here rˆ , ϕˆ , and ẑ are the unit vectors for the cylindrical coordinates r , ϕ , and z ,respectively. The quantity Rj( p), where j = r,ϕ , z , produces a random number between − pand p . The coefficients Dr, D ϕ , and Dzare the cylindrical-coordinate components of thediffusion coefficient D = Dr+ D ϕ + Dz, which are responsible for the motion alo<strong>ng</strong> the rˆ , ϕˆ ,and ẑ directions, respectively.68
Nhữ<strong>ng</strong> tiến bộ <stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>> <stro<strong>ng</strong>>Qua<strong>ng</strong></stro<strong>ng</strong>> học, <stro<strong>ng</strong>>Qua<strong>ng</strong></stro<strong>ng</strong>> <stro<strong>ng</strong>>ph</stro<strong>ng</strong>>ổ và Ứ<strong>ng</strong> dụ<strong>ng</strong> <strong>VI</strong> <strong>ISSN</strong> <strong>1859</strong> - <strong>4271</strong>Fig. 6: Time dependences of (a) the radial distance r , the azimuthal a<strong>ng</strong>le ϕ , (c) the axial positionz , and (d) the mean intracavity <stro<strong>ng</strong>>ph</stro<strong>ng</strong>>oton number N for a movi<strong>ng</strong> atom. The initial position and velocity ofthe atom are ( r = 3a= 600 nm,ϕ = 0, z = 0 ) and ( vr = 0, vϕ= 0, vz= 20vrec ≅ 7 cm/s ), respectively.The probe field, the cavity, and the atom are at exact resonance, that is, ωp= ωc= ω0. Otherparameters are as in Fig. 2.We plot in Fig. 6 the time dependences of the atomic position ( r,ϕ , z)and the mean†intracavity <stro<strong>ng</strong>>ph</stro<strong>ng</strong>>oton number N = 〈 a a〉. In these calculations, we assume that the initial velocityof the atom is aligned alo<strong>ng</strong> the direction z of the fiber axis and that the atom stops movi<strong>ng</strong> afterit reaches the fiber surface. Figure 6(c) shows that the velocity of the atom alo<strong>ng</strong> the direction zof the fiber axis is almost constant. The reason is the followi<strong>ng</strong>: Since the probe field is weak andis at exact resonance with the cavity and the atom, the optical force is small. Consequently, thetranslational motion of the atom in the case of Fig. 6 is mainly determined by the van der Waalsforce, which acts only in the radial direction. Figure 6(d) shows that the mean intracavity <stro<strong>ng</strong>>ph</stro<strong>ng</strong>>otonnumber N oscillates in time, with peaks and dips at the nodes and an<s<stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>>>tin</s<stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>>>odes of the cavitystandi<strong>ng</strong>-wave field, respectively. Such oscillations closely follow the motion of the atomthrough the nodes and an<s<stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>>>tin</s<stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>>>odes of the cavity field.<strong>VI</strong>. CONCLUSIONSWe have calculated numerically and analytically the force, the friction coefficients, and themomentum diffusion. We have found that the spatial dependences of the force, the frictioncoefficients, and the momentum diffusion are very complicated due to the evanescent-wavenature of the atom-field coupli<strong>ng</strong> and the effect of the van der Waals potential. We have shownthat the time development of the mean number of <stro<strong>ng</strong>>ph</stro<strong>ng</strong>>otons in the cavity closely follows thetranslational motion of the atom through the nodes and an<s<stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>>>tin</s<stro<strong>ng</strong>>tro<strong>ng</strong></stro<strong>ng</strong>>>odes of the fiber-guided cavitystandi<strong>ng</strong>-wave field even though the cavity finesse is moderate, the cavity is lo<strong>ng</strong>, and thedrivi<strong>ng</strong>-field power is very weak.69