Theory of an optical dipole trap for cold atoms

Size: px
Start display at page:

Download "Theory of an optical dipole trap for cold atoms"

Transcription

1 Theory of an optical ipole trap for col atoms Article (Publishe Version) Garraway B M an Minogin V G () Theory of an optical ipole trap for col atoms. Physical Review A 6 (4) ISSN This version is available from Sussex Research Online: This ocument is mae available in accorance with publisher policies an may iffer from the publishe version or from the version of recor. If you wish to cite this item you are avise to consult the publisher s version. Please see the URL above for etails on accessing the publishe version. Copyright an reuse: Sussex Research Online is a igital repository of the research output of the University. Copyright an all moral rights to the version of the paper presente here belong to the iniviual author(s) an/or other copyright owners. To the extent reasonable an practicable the material mae available in SRO has been checke for eligibility before being mae available. Copies of full text items generally can be reprouce isplaye or performe an given to thir parties in any format or meium for personal research or stuy eucational or not-for-profit purposes without prior permission or charge provie that the authors title an full bibliographic etails are creite a hyperlink an/or URL is given for the original metaata page an the content is not change in any way.

2 Theory of an optical ipole trap for col atoms B. M. Garraway Sussex Centre for Optical an Atomic Physics School of Chemistry Physics an Environmental Sciences University of Sussex Falmer Brighton BN 9QJ Englan V. G. Minogin Institute of Spectroscopy Russian Acaemy of Sciences 49 Troitsk Moscow Region Russia Receive 6 April ; publishe 5 September The theory of an atom ipole trap compose of a focuse far re-etune trapping laser beam an a pair of re-etune counterpropagating cooling beams is evelope for the simplest realistic multilevel ipole interaction scheme base on a moel of a 35-level atom. The escription of atomic motion in the trap is base on the quantum kinetic equations for the atomic ensity matrix an the reuce quasiclassical kinetic equation for atomic istribution function. It is shown that when the etuning of the trapping fiel is much larger than the etuning of the cooling fiel an with low saturation the one-photon absorption emission processes responsible for the trapping potential can be well separate from the two-photon processes responsible for sub-doppler cooling atoms in the trap. Two conitions are erive that are necessary an sufficient for stable atomic trapping. The conitions show that stable atomic trapping in the optical ipole trap can be achieve when the trapping fiel has no effect on the two-photon cooling process an when the cooling fiel oes not change the structure of the trapping potential but changes only the numerical value of the trapping potential well. It is conclue that the separation of the trapping an cooling processes in a pure optical ipole trap allows one to cool trappe atoms own to a minimum temperature close to the recoil temperature keeping simultaneously a eep potential well. PACS numbers: 3.8.Pj 4.5.Vk PHYSICAL REVIEW A VOLUME I. INTRODUCTION In recent years there has been a growing interest in eveloping traps for col atoms with applications to funamental physical problems nanotechnology high-resolution microwave an optical spectroscopy frequency stanars an atom optics 5. Among ifferent types of atom traps stuie in recent years of practical importance is a far-offresonance optical ipole trap FORT base on a single focuse re-etune laser beam 6. The FORT prouces a nearly conservative potential well for atoms but incorporates the inevitable heating ue to the photon recoil associate with the scattere laser light 89. Although the heating rate may be very small at very large etuning from the resonance the photon recoil heating inevitably introuces an upper limit on the lifetime of atoms in the trap. It was previously propose that the heating mechanism in the FORT might be suppresse by aing the cooling laser fiels to a focuse trapping laser beam 78. Latest experiments with ifferent types of cooling laser fiels have shown that the aition of the cooling fiel can increase the lifetime of atoms an even the atomic ensity in the trap. The aition of the cooling fiel changes not only the Boltzmann factor that efines the lifetime of atoms in the FORT but may also have a profoun effect on all the basic parameters of the trap since the cooling fiel may strongly influence the atomic populations an coherences. Typically any cooling laser fiel operates at a etuning less than the etuning of the trapping laser beam. The cooling fiel may thus be responsible not only for perturbation of the steaystate internal atomic state but the perturbation of the trapping potential as well. Up to now the FORT with an aitional cooling laser fiel was theoretically iscusse for the simplest moel of a two-level ipole interaction scheme 73. A two-level moel has however a very limite connection with real experimental techniques that typically explore multilevel ipole interaction schemes. Physically there is a major ifference between the moels of a two-level atom an a multilevel atom in applications relate to the trapping an cooling atoms. In a two-level moel both trapping an cooling fiels excite atoms on the same atomic transition. As a result for a two-level atomic scheme in the FORT the epth of the potential well an the cooling limit have generally the same orer of magnitue 73 efine by a so-calle Doppler temperature 4. In multilevel atomic schemes the trapping an cooling laser fiels can prouce principally ifferent atomic transitions. The trapping fiel can prouce a potential well ue to the one-photon transitions while the cooling laser fiel can cool atoms own to the sub-doppler temperatures 5 8 ue to the two-photon transitions 9. It is thus important that in multilevel atomic schemes the optical processes use for trapping atoms in the FORT an those for sub-doppler cooling can have the ifferent physical origin. The use of the ifferent optical processes for trapping an cooling multilevel atoms thus raises new questions on basic parameters of the trap an the lifetime of atoms achievable in the FORT with a superimpose sub-doppler cooling process. In this paper we present a theoretical analysis of the FORT compose of a single re-etune trapping laser beam an a re-etune cooling fiel compose of the counterpropagating circular-polarize laser waves. To be specific we evelop an analytical theory for a simplest realistic moel of a 35-level atom interacting with one-imensional 5-947//64/4346/$ The American Physical Society

3 B. M. GARRAWAY AND V. G. MINOGIN PHYSICAL REVIEW A processes while a less re-etune - laser fiel is responsible for laser cooling ue to the two-photon optical processes. In the chosen interaction scheme the fiel of a linearly polarize trapping laser beam E t e E t rcosky t t is efine by a unit polarization vector e e z spatially nonuniform amplitue E t (r) an the wave vector k t ke y (k t /c). The cooling laser fiel is assume to be compose of two counterpropagating left circularly polarize waves E c E ce e i(kz c t) e e i(kz c t) E ce e i(kz c t) e e i(kz c t) FIG.. One-imensional scheme of the FORT compose of a re-etune trapping laser beam tr an the re-etune cooling laser fiel c in the - configuration a an the ipole interaction scheme of a 35-level atom with the trapping an cooling laser fiels b. The magnetic sublevels for the groun an excite states of the atom are enote as: j g mg m j e m e m. counterpropagating laser waves in a - configuration. In this moel the ipole potential is cause by the one-photon processes while the friction force comes from the twophoton processes. Our analysis shows that uner specific conitions formulate below the one-photon processes proucing the optical ipole potential can have a negligible influence on the twophoton processes responsible for sub-doppler cooling of atoms. On the other han the perturbation of the atomic internal state prouce by the two-photon cooling processes may not change the structure of the potential well but only the value of the potential well epth. The separation of the trapping an cooling processes achievable uner specifie conitions allows one to maintain the epth of the potential well in the FORT that excees consierably the kinetic energy of trappe atoms. As a result the lifetime of atoms associate with the photon scattering in the FORT can be extremely long thus leaving other processes such as atom collisions an the reabsorption of laser light 4 to be responsible for the atom storage time. II. MULTILEVEL INTERACTION SCHEME The simplest scheme for a FORT may be escribe by the ipole interaction of a 35-level atom with a linearly polarize focuse trapping laser beam an one-imensional cooling - laser-fiel configuration as shown in Fig.. In the scheme a far re-etune laser beam is assume to be responsible for atomic trapping ue to the one-photon optical where e (/)(e x ie y ) are spherical unit vectors c is the frequency of the cooling laser waves an k c /c is the magnitue of the wave vector. With respect to the quantization axis Oz the first wave in Eq. is a polarize wave an the secon one is a polarize wave. Note that in a nonrelativistic approach the magnitue of the wave vector for the counterpropagating cooling laser fiels is consiere to be the same. For the above interaction scheme the atomic Hamiltonian has a stanar form HH M D E where the Hamiltonian H escribes the internal atomic states with energy levels E g E g an E e E e E e an the last term escribes the ipole interaction between the atom an the total laser fiel EE t E c. III. UNPERTURBED TRAPPING POTENTIAL Consier first the trapping potential for the above interaction scheme in the absence of the cooling fiel. The quasiclassical atomic-ensity-matrix equations escribing the ipole interaction of a 35-level atom with the fiel in a rotating wave approximation RWA are liste in Appenix A. The set inclues only equations for the ensity matrix elements that o not ecay to zero for an interaction time t. The ensity matrix elements entering the equations of Appenix A are efine with respect to the timeepenent atomic wave functions. The equations inclue on the left-han sie the total convective time erivatives t t v r. The Rabi frequency for the trapping laser beam an the etuning are: 3 4 r E tr 3 t

4 THEORY OF AN OPTICAL DIPOLE TRAP FOR COLD ATOMS PHYSICAL REVIEW A where is the reuce ipole matrix element an is the atomic transition frequency. The Rabi frequency is efine with respect to the strongest -type ipole transition j g m g j e m e. The spontaneous ecay rate W sp is efine in the usual way as W sp c 3. 6 Eliminating an explicit time an position epenence in the equations for the ensity matrix elements with the substitutions g e g e e i(kyt) g e g e e i(kyt) g e g e e i(kyt) an solving next the equations for a steay state one can fin the ipole raiation force on a 35-level atom with the use of a well-justifie quasiclassical formula F D"E where DD ab ba is an expectation value of the atomic ipole moment an the graient is assume to act on the electric fiel E only. The raiation force on a 35-level atom in a laser beam inclues as usual two partial forces i.e. the graient force an the raiation pressure force FF gr F rp F gr rre g e 3 g e g e F rp ke y Im e g 3 e g e g The case of large negative etuning is of particular interest in this paper i.e. we will let (r) an for a motionless atom the forces are then F gr 3 7 r F rp 3 7 ke r y. The graient force prouces the potential well for col atoms that iffers from that for a two-level atom 7 by a numerical factor only that is U F gr r 3 r 7. 3 The raiation pressure force prouces an aitional asymmetric potential well that can be neglecte at large etuning. IV. QUANTUM KINETIC EQUATIONS In general to escribe the atomic motion we can use the approach base on the quantum kinetic equations for the atomic ensity matrix in the Wigner representation. This is relatively straightforwar for a scheme such as the one consiere in Fig. where the atom interacts with both laser fiels an. Since we wish to analyze the operation of the trap at large etunings it is sufficient for our purpose to consier the quantum kinetic equations at small optical saturation. We will take into account one- an two-photon optical processes an neglect higher-orer optical processes. To simplify the consieration of the equations of motion we write out the initial microscopic equations but we neglect first the spatial variation in the trapping laser beam amplitue E t (r). The effect of the spatially inhomogeneous trapping laser beam will be taken into account in Sec. VI. A set of the atomic ensity matrix equations in the Wigner representation an a rotating wave approximation RWA escribing an interaction of a 35-level atom with laser waves of constant amplitue at weak optical saturation is liste in Appenix B. The set only inclues the equations for the ensity matrix elements escribing one- an two-photon processes. As before the ensity matrix elements entering the equations of Appenix B are efine with respect to the time-epenent atomic wave functions. The upper inices for the momentum-shifte ensity matrix elements are chosen to have ifferent meaning for the terms coming from the interaction with a trapping an cooling fiel i.e. () ab a rp ke y tb () ab a rp ke z tb ab arpnktb 4 where n is a unit vector that efines the irection of the spontaneous photon emission. The total time erivative on the left-han sie of the equations is efine as before by Eq. 4. The Rabi frequency for the cooling laser fiel efine with respect to the strongest -type ipole transition an the etuning of the cooling fiel are E c 5 c. 5 Functions with efine the angular anisotropy of the spontaneous photon emission n 3 6 n z n 3 8 n z 6 where n z is the projection of the unit vector n on the quantization axis Oz. Numerical factors entering the ensity matrix equations are ue to the Clebsch-Goran coefficients efining the ipole matrix elements for specific atomic transitions. Our choice of sign of the reuce ipole moment can be seen from the values of the matrix elements for the strongest -type an -type atomic transitions

5 B. M. GARRAWAY AND V. G. MINOGIN PHYSICAL REVIEW A e g g e 5 e g g e e g g e 5. Note that when the cooling laser fiel an photon recoil are neglecte the equations of Appenix B reuce to the quasiclassical equations of Appenix A. The quantum kinetic equations of Appenix B are use below to erive the quasiclassical kinetic equation of atomic motion that efines the ipole raiation force on the atom in the total laser fiel the momentum iffusion tensor an the potential well for the motionless atom. V. QUASICLASSICAL KINETIC EQUATION The initial quantum kinetic equations for the atomic ensity matrix elements can be first simplifie by the stanar proceure of eliminating the explicit time an position epenence with proper substitutions. Since the initial equations inclue only terms escribing one- an two-photon processes we introuce substitutions that also take into account only one- an two-photon processes: g e g e e ikzit g e g e e ikyit g e g e e ikzit g e g e e ikzit g e g e e ikyit g e g e e ikzit g e g e e ikzit g e g e e ikyit g g g g g g g g e e e e e e e e g e g e e ikzit e ik(yz)i()t g g e ik(yz)i()t g g g g e ikz e ik(yz)i()t g g e ik(yz)i()t e e e e e ik(yz)i()t e e e e e ikz e ik(yz)i()t e e e ik(yz)i()t e e e e e ikz e ik(yz)i()t e e e ik(yz)i()t e e e e e ikz e e e e e ik(yz)i()t. 7 Having mae the above substitutions we neglect in the microscopic equations all terms that escribe three-photon an higher-orer multiphoton processes. This proceure finally results in a set of equations that oes not inclue explicit time an position epenence. An analysis of the ensity matrix equations that have no explicit time an coorinate epenence can be one in the usual way 3. Uner the conition that atom-laser fiel interaction time excees the characteristic relaxation times of the atomic-ensity-matrix the momentum with of the atomic ensity-matrix elements can be assume to excee the photon momentum k. This principal assumption which will be checke below allows one to expan the ensitymatrix elements in powers of the photon momentum k. Consiering next the time evolution of the equations expane in successively increasing orers of the photon momentum k one can conclue that the iagonal elements aa an the off-iagonal elements ab are functionals of the Wigner quasiprobability istribution function ww(rpt) w g g e e. ; 8 The structure of the functional epenence for the iagonal an off-iagonal ensity matrix elements can be irectly etermine from the structure of the expane equations: aa R aa kr aa w kq aa w p z ab S ab ks ab w kt ab w p z 9 where R aa S ab R aa S ab...q aa T ab... are functions of atomic momentum p or atomic velocity vp/m that have to be etermine by the solution proceure. In accorance with the efinition of the quasi-probability istribution function 8 the unknown iagonal functions satisfy the normalization conitions R aa R aa Q aa.... Substituting next the general solution 9 into the expane equations for the atomic ensity-matrix elements one can see that the Wigner function ww(rpt) satisfies an infinite Fokker-Planck equation. To the secon orer in the photon momentum k the equation for the quasiprobability istribution function ww(rpt) is then escribe by a secon-orer Fokker-Planck kinetic equation w t w v r p Fw p D iiw i where ixyz F is the ipole raiation force on the atom an D ii is the momentum iffusion tensor

6 THEORY OF AN OPTICAL DIPOLE TRAP FOR COLD ATOMS PHYSICAL REVIEW A In this section we consier only homogeneous laser waves where the ipole raiation force has the clear meaning of a raiation pressure force FF rp e z F rp. This force an the momentum iffusion tensor entering Eq. are F rp k ImS e g S e g 6 S e g S e g D ii k ii R e e iz k R e e ii ImT e g S e g S e g R e e 3 R e e R e e R e e ImT g e T g e 3 R e e R e e T e g 3 ImT g e T g e. 3 In the above equations the coefficients ii ii etermine the angular anisotropy of the spontaneous photon emission ii nn i n xx yy 3 zz 5 ; ii nn i n xx yy 5 4 zz 5. 5 The values of the force an momentum iffusion coefficients F rp an D ii which govern the time evolution of the quasiclassical istribution function can be explicitly etermine by solving the steay-state equations that follow from the expane equations for the atomic-ensity-matrix elements. These are consiere separately in zero an first orer in photon momentum k. Explicit equations for the quantities F rp an D ii are given in Sec. VIII. VI. DIPOLE GRADIENT FORCE We next take into consieration the ipole graient force F gr associate with the graient of the trapping laser beam amplitue E t E t (r). The erivation of the graient force can be one in a way that generalizes the proceure consiere in the preceing section. Representing the trapping laser beam amplitue in the form of a Fourier expansion ir ab p 3/ iqe iq"r ab p q 3 q. 7 For a laser beam amplitue that varies in space on a scale that is large compare to the size of the atomic wave packet we use an expansion of the ensity-matrix elements to first orer in the small momentum q ab p q ab p q p abp. This transforms the terms in equations of Appenix B into ir ab p ir ab p r p abp. 8 9 When the substitutions 9 are introuce into the equations of Appenix B the Fokker-Planck equation inclues in the above approximation the total raiation force as a sum of the ipole graient force an the raiation pressure force FF gr F rp. 3 The ipole graient force is then etermine by the steaystate optical coherences as F gr rre S g e 3 S g e S g e. 3 Note that formal mathematical expression for this part of the total raiation force coincies with that given for the graient force by Eq. 9. The explicit expressions given by Eqs. 3 an 9 are generally ifferent since the functions S ab escribe the interaction of the atom with the total laser fiel efine by the equations of Appenix B while the functions ab escribe the interaction of the atom with only the trapping laser beam efine by the equations of Appenix A. We will see in Sec. VIII that Eq. 3 gives the graient force 4 that iffers by a numerical factor from the graient force given by Eq. 9. VII. LOW-SATURATION SOLUTION We will etermine the steay-state atomic ensity matrix elements which efine the partial forces an the momentum iffusion tensor in a regime characterize by two important limits. That is we consier large etunings where an we consier low optical saturation 3 E t r 3/ E t qe iq"r 3 q 6 s t s c 33 one shoul introuce into the equations of Appenix B the following substitutions: when the one-photon an two-photon processes play a ominant role in the time evolution of the atomic-ensity-matrix

7 B. M. GARRAWAY AND V. G. MINOGIN PHYSICAL REVIEW A elements. In aition to the above approximations we restrict our treatment to the case of slowly moving atoms where where the steay-state groun-state populations N R g g N R g g an N R g g connecte to the groun-state two-photon coherence g g are y kv y z kv z. 34 Uner all these conitions the sets of equations for the functions R aa S ab an T ab which follow from the equations of Appenix B expane to the first orer in the photon momentum have a simple analytical solution. When solving the last equations one can note that the efficiencies of the twophoton optical processes between the groun-state magnetic sublevels epen crucially on two parameters an N D 3 N D N D 6 k v z an the common enominator is kv z9k v z kv z9k v z D 7 4 k v z which efine the frequency withs of the two-photon resonance structures relate to the groun-state coherences g g an g g g g. Physically the origin of the twophoton frequency withs can be unerstoo in the following way. The laser fiel connects the groun-state probability amplitues with the upper-state probability amplitues through one-photon excitation processes. These one-photon absorption processes perturb the groun-state probability amplitues an then lea to ecay rates for the groun-state coherences of the orer of the rates of the one-photon transitions i.e. of the orer of an. These quantities play accoringly the role of the frequency withs for the twophoton resonance processes relate to the groun-state coherences. Since the two-photon processes inuce by the cooling laser fiel are only important for eep sub-doppler cooling of atoms in the trap 9 the physical meaning of parameters an allows us to introuce two conitions necessary for separating the cooling an trapping processes: i.e.. 36 The left inequality guarantees that the trapping fiel oes not influence the two-photon cooling process while the right inequality is neee for the two-photon friction coefficient to be greater than the one-photon friction coefficient. Uner the conition 36 an in lowest orer of the small parameters s t s c y an z the low-saturation lowvelocity solutions give the one-photon coherences entering Eq. 3 as S g e i3/ i N S g e i i N S g e i3/ i N 37 The other steay-state coherences an populations entering Eqs. an 3 are efine by similar equations S g e i i N S g e i/6 i N S g e i/ i N... 4 Qualitatively the behavior of the groun-state populations at low velocities is efine by the groun-state coherence g g that is inuce by the cooling fiel. The velocity position of the narrow two-photon structures in the grounstate populations can be seen from the energy conservation law. In an atom rest frame the absorption of a photon from one traveling wave an emission of a photon into the other traveling wave results in a two-photon transition between groun-state sublevels g g that oes not change the atom energy (kv)(kv). Energy conservation thus irectly shows that the two-photon resonance structures which are ue to the cooling fiel are locate at zero velocity kv. VIII. MODIFIED OPTICAL POTENTIAL DEPTH AND KINETIC ENERGY The atomic coherences that we have erive Eqs. 37 together with the values of the groun-state populations give a new value of the graient force at zero velocity an negative etuning F gr r. 4 In a corresponing way they give a new value of the potential Ur F gr r 55 r r

8 THEORY OF AN OPTICAL DIPOLE TRAP FOR COLD ATOMS PHYSICAL REVIEW A we can erive an explicit expression for the friction coefficient ue to the two-photon processes 6: 7 r 44 where r k /M is a recoil frequency. The last expression shows that at low effective saturation the friction ue to the two-photon processes oes not epen on intensity. The twophoton friction is however much higher than that ue to the one-photon processes. This higher friction is naturally explaine by the smaller value of the velocity with of the two-photon structure /k compare to the with of the onephoton structure /k see the right inequality 36. In a way that is similar to the raiation force the momentum iffusion tensor D ii also inclues a narrow two-photon velocity structure locate at zero velocity Fig.. The contribution of the two-photon structure can be shown to ecrease the value of the iffusion coefficient DD zz at zero velocity. The reason is that the analysis of the equations liste in Appenix B shows that the secon part of the iffusion coefficient in Eq. 3 is much bigger than the first one. However this secon part of the iffusion coefficient D manifests a narrow velocity ip at velocity vv z. In our basic approximations the iffusion coefficient D at zero velocity D can be estimate as FIG.. Velocity epenence of the raiation pressure force F rp soli line an iffusion coefficient DD zz ashe line in a small velocity region kv z at negative etuning an saturation parameter G5. which is about. times less than the unperturbe potential 3. The reuce value of the graient force an the potential in the total laser fiel is naturally explaine by the ominant role of the cooling fiel in proucing atomic populations in accorance with the leftmost conition in Eq. 36. While the trapping fiel alone prouces the groun-state atomic populations at zero velocity as N N 3/34 N 4/7 the cooling fiel reistributes the populations to the values N N 9/ N 4/. This reistribution multiplie by the relative strengths of the ipole transitions gives the ecrease in the potential. However in a threeimensional cooling fiel the ecrease of the potential shoul be less profoun. Substitution of the steay-state atomic coherences an groun-state populations into Eq. gives an explicit expression for the raiation pressure force in a low-saturation low-velocity approximation. Since the above analytical expressions 38 an 4 neglect a weak velocity epenence ue to the one-photon Doppler processes the raiation pressure force shown in Fig. is only escribe in our analysis near zero velocity i.e. where v /k. In the case of a re-etune cooling fiel where the low-velocity part of the raiation pressure force reuces to the cooling force. If we represent the force at a re etuning an at low velocities in the stanar form of a friction force F rp Mv 43 D 46 7 k. 45 The iffusion coefficient an the friction coefficient jointly efine the kinetic energy of the trappe atoms an thus an effective temperature foun accoring to the Einstein relation as Ek B T D 3 M 3. IX. CONDITIONS FOR STABLE TRAPPING 46 Assuming now that kinetic energy 46 of col trappe atoms is much less than the epth U() of the potential well 4 one can get a sufficient conition for stable atomic trapping. 47 Comparing the last conition with the conition for eep laser cooling 36 one can see that both conitions can be satisfie when the etuning of the trapping fiel is much larger than that of the cooling fiel. 48 The above two conitions efine by Eqs. 36 an 47 thus justify the iea of a stable ipole trap. The stable atomic trapping in the optical ipole trap can thus be achieve when the trapping fiel has no effect on the two-photon cooling process an when the cooling fiel oes not change the

9 B. M. GARRAWAY AND V. G. MINOGIN PHYSICAL REVIEW A structure of the trapping potential but changes only the numerical value of the trapping potential well. The lifetime of atoms in the trap associate with the iffusive heating can finally be estimate as e U()/E 49 where is the oscillation frequency of an atom in the trap an the Boltzmann factor U()/E consierably excees unity. The above evaluations of the conitions necessary an sufficient for stable atomic trapping can be illustrate by a specific example. Assume that the Rabi frequency an the etuning for a cooling fiel are accoringly an an the Rabi frequency an the etuning for the trapping fiel are an 4. For these parameters the one-photon withs 4 an.5 3 satisfy conitions 36 for eep laser cooling. In their turn the kinetic energy 46 estimate as E.4 an the potential-well epth estimate accoring to Eq. 4 as U().8 satisfy the sufficient conition 47 for stable atomic trapping giving for the Boltzmann factor a sufficiently large value U()/E. APPENDIX A Below is liste a set of quasiclassical equations for the atomic-ensity-matrix elements in the RWA that escribe the ipole interaction of a 35-level atom with the trapping laser fiel efine by Eq.. Note that the set inclues only the equations that involve ensity-matrix elements not vanishing at the interaction time t. t g g i3 ei(kyt) e g c.c. e e 3 e e t g g ie i(kyt) e g c.c. e e 4 3 e e e e t g g i3 ei(kyt) e g c.c. 3 e e e e X. CONCLUSION The above analysis shows thus that in an optical ipole trap one-photon processes responsible for trapping atoms can be well separate from the two-photon processes responsible for the eep sub-doppler cooling of atoms. This separation of the trapping an cooling processes ensures in turn a stable operation of the trap accoring to Eq. 49. This principal conclusion is also vali for the three-imensional case since if one as new - cooling fiels along the axes Ox an Oy all the above conclusions remain vali with only the numerical coefficients changing. This justifies the possibility of separating the trapping an cooling processes in a three-imensional case. Finally note that the above quasiclassical analysis cannot give a quantitative limit for the temperature of laser-coole atoms since it consiers atomic velocities higher than the recoil velocity. The cooling limit can however be seen from a qualitative estimation. Accoring to the above analysis the two-photon cooling process can operate when the velocity with of the two-photon resonance excees the recoil velocity v/kv r. That is when the atomic energy excees the recoil energy EE r k /M. These limits give the physical justification for the key assumption that the atomicensity-matrix elements can be expane in powers of the photon momentum k. ACKNOWLEDGMENT This work was supporte in part by the Royal Society an the Russian Fun for Basic Research Grant No t e e i3 ei(kyt) g e c.c. e e t e e ie i(kyt) g e c.c. e e t e e i3 ei(kyt) g e c.c. e e t g e i3 ei(kyt) g g e e g e t g e ie i(kyt) g g e e g e t g e ie i(kyt) g g e e g e. APPENDIX B Below are liste the atomic-ensity-matrix equations in the Wigner representation an the RWA. They escribe the ipole interaction of a 35-level atom with trapping an cooling laser fiels efine by Eqs. an. The equations inclue only terms escribing one-photon an twophoton optical processes

10 THEORY OF AN OPTICAL DIPOLE TRAP FOR COLD ATOMS PHYSICAL REVIEW A t g g i3 ei(kyt) () e g ie ikzit () e g i 6 eikzit () e g c.c. n e e n e e 3 n e e n t g g ie i(kyt) () e g i eikzit () e g t g g i3 i eikzit () e g c.c. n e e 4 3 n e e ei(kyt) () e g n e e n ie ikzit () e g i 6 eikzit () e g c.c. 3 n e e n e e n e e n t e e ie ikzit () g e c.c. e e t e e i3 ei(kyt) () g e i eikzit () g e c.c. e e t e e ie i(kyt) () g e i 6 eikzit () g e e ikzit () g e c.c. e e t e e i3 ei(kyt) () g e i eikzit () g e c.c. e e t e e ie ikzit () g e c.c. e e t g e i3 ei(kyt) () e e ie ikzit () e e e ikzit () g g i 6 eikzit () e e g e t g e i3 ei(kyt) () g g () e e ie ikzit () e e i 6 eikzit () e e i eikzit () g g g e i(kyt) t () g e ie g g 3 e e () ie ikzit () e e i 6 eikzit () g g e ikzit () g g e ikzit () e e g e t g e i3 ei(kyt) g g (t) () e e i eikzit () g g e ikzit () e e e ikzit () e e g e t g e ie i(kyt) () g g e e () i eikzit (t) e e e ikzit () e e i 6 eikzit () g g e ikzit () g g g e t g e i3 ei(kyt) () g g () e e i eikzit () g g e ikzit () e e e ikzit () e e g e

11 B. M. GARRAWAY AND V. G. MINOGIN PHYSICAL REVIEW A () () 3 () i(kyt) t () g e ie g g 3 e e t g g i e i(kyt) e g ei(kyt) g e ie ikzit () e e i 6 eikzit () g g e ikzit () g g e ikzit () e e g e i 6 eikz e it () e g e it g e 3 n e e () 3 n e e t g e i3 t g e i3 ei(kyt) () g g () e e i 6 eikzit () e e ie ikzit () e e i eikzit () g g g e ei(kyt) () e e ie ikzit () g g e ikzit () e e i 6 eikzit () e e g e t g g i 3 t g g i ei(kyt) () e g e i(kyt) () g e i 6 eikz e it () e g e it () g e n e e 3 n e e 3 n e e n 6 eikzit () e g e ikzit () g e 3 n e e n e e 3 n e e n t e e i3 n e e n ei(kyt) () e g e e ie ikzit () g e t e e ie ikzit () g e i 6 eikzit () e g t e e i 3 t e e i e e ei(kyt) () g e e i(kyt) () e g i eikzit () g e i 6 eikzit () e g e e eikz e it () g e e it () e g e e t e e i e i(kyt) () g e 3 ei(kyt) () e g i eikzit () e g i 6 eikzit () g e e e t e e ie ikzit () e g i 6 eikzit () g e t e e i3 e e ei(kyt) () g e e e. ie ikzit () e g 4346-

12 THEORY OF AN OPTICAL DIPOLE TRAP FOR COLD ATOMS PHYSICAL REVIEW A Laser Manipulation of Atoms an Ions Proceeings of the International School of Physics Enrico Fermi Course CX- VIII Varenna Italy 99 eite by E. Arimono W. D. Phillips an F. Strumia North-Hollan Amsteram 99. H. Metcalf an P. van er Straten Phys. Rep V. S. Letokhov M. A. Olshanii an Yu. B. Ovchinnikov Quantum Semiclassic. Opt C. Aams an E. Riis Prog. Quantum Electron R. Grimm M. Weiemüller an Yu. B. Ovchinnikov Av. At. Mol. Opt. Phys A. Ashkin Phys. Rev. Lett J. P. Goron an A. Ashkin Phys. Rev. A S. Chu J. E. Bjorkholm A. Ashkin an A. Cable Phys. Rev. Lett J. D. Miller R. A. Cline an D. J. Heinzen Phys. Rev. A 47 R K. L. Corwin S. J. M. Kuppens D. Cho an C. E. Wieman Phys. Rev. Lett H. Lee C. Aams M. Kasevich an S. Chu Phys. Rev. Lett D. Boiron A. Michau J. M. Fournier L. Simar M. Sprenger G. Grynberg an C. Salomon Phys. Rev. A 57 R V. G. Minogin an V. S. Letokhov Laser Light Pressure on Atoms Goron an Breach New York V. S. Letokhov V. G. Minogin an B. D. Pavlik Zh. Éksp. Theor. Fiz Sov. Phys. JETP P. D. Lett R. N. Watts C. I. Westbrook W. D. Phillips P. L. Goul an H. J. Metcalf Phys. Rev. Lett J. Dalibar an C. Cohen-Tannouji J. Opt. Soc. Am. B D. Boiron C. Triché D. R. Meacher P. Verkerk an G. Grynberg Phys. Rev. A 5 R A. Hemmerich M. Weiemüller T. Esslinger C. Zimmermann an T. W. Hänsch Phys. Rev. Lett S. Chang T. Y. Kwon Ho S. Lee an V. Minogin Phys. Rev. A S. Chang T. Y. Kwon Ho S. Lee an V. G. Minogin Phys. Rev. A M. Prentiss A. Cable J.E. Björkholm S. Chu E.L. Raab an D.E. Pritchar Opt. Lett D. Sesko T. Walker C. Monroe A. Gallagher an C. Wieman Phys. Rev. Lett C. D. Wallace T. P. Dinneen K.-Y. N. Tan T. T. Grove an P. L. Goul Phys. Rev. Lett D. W. Sesko T. G. Walker an C. E. Wieman J. Opt. Soc. Am. B

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

PHYS 414 Problem Set 2: Turtles all the way down

PHYS 414 Problem Set 2: Turtles all the way down PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field PHYSICAL REVIEW A VOLUME 60, NUMBER 2 AUGUST 1999 Quantum optics of a Bose-Einstein conensate couple to a quantize light fiel M. G. Moore, O. Zobay, an P. Meystre Optical Sciences Center an Department

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Average value of position for the anharmonic oscillator: Classical versus quantum results

Average value of position for the anharmonic oscillator: Classical versus quantum results verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

A Model of Electron-Positron Pair Formation

A Model of Electron-Positron Pair Formation Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Convective heat transfer

Convective heat transfer CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences

More information

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION Long-istance coherent coupling in a quantum ot array Floris R. Braakman 1, Pierre Barthelemy 1, Christian Reichl, Werner Wegscheier, Lieven M.K. Vanersypen 1 1 Kavli Institute of Nanoscience, TU Delft,

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

Extinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0

Extinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0 Extinction, σ/area 1.5 1.0 t = t 0 t = 0.7 t 0 t = t 0 t = 1.3 t 0 t = 1.5 t 0 0.7 0.9 1.1 Energy (ev) = 20 nm t 1.3 Supplementary Figure 1: Plasmon epenence on isk thickness. We show classical calculations

More information

Online Appendix for Trade Policy under Monopolistic Competition with Firm Selection

Online Appendix for Trade Policy under Monopolistic Competition with Firm Selection Online Appenix for Trae Policy uner Monopolistic Competition with Firm Selection Kyle Bagwell Stanfor University an NBER Seung Hoon Lee Georgia Institute of Technology September 6, 2018 In this Online

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

arxiv:quant-ph/ v1 29 Jun 2001

arxiv:quant-ph/ v1 29 Jun 2001 Atomic wave packet basis for quantum information Ashok Muthukrishnan an C. R. Strou, Jr. The Institute of Optics, University of Rochester, Rochester, New York 14627 (March 15, 2001) arxiv:quant-ph/0106165

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses.

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses. Physics 195a Course Notes The K 0 : An Interesting Example of a Two-State System 021029 F. Porter 1 Introuction An example of a two-state system is consiere. involve: 1. Treatment of a ecaying particle.

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

PLASMA ASSISTED CO 2 DISSOCIATION MODELS FOR ENVIRONMENT, ENERGY AND AEROSPACE APPLICATIONS

PLASMA ASSISTED CO 2 DISSOCIATION MODELS FOR ENVIRONMENT, ENERGY AND AEROSPACE APPLICATIONS PLASMA ASSISTED CO 2 DISSOCIATION MODELS FOR ENVIRONMENT, ENERGY AND AEROSPACE APPLICATIONS G. Colonna, L. D. Pietanza, G. D Ammano, A. Laricchiuta, an M. Capitelli CNR-IMIP, via Amenola 122/D, 70126 Bari

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

Advanced Partial Differential Equations with Applications

Advanced Partial Differential Equations with Applications MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.

More information

Three-Dimensional Theory of the Magneto-Optical Trap

Three-Dimensional Theory of the Magneto-Optical Trap ISSN 1063-7761, Journal of Experimental and Theoretical Physics, 015, Vol. 10, No. 4, pp. 587 594. Pleiades Publishing, Inc., 015. Original Russian Text O.N. Prudnikov, A.V. Taichenachev, V.I. Yudin, 015,

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In the usual geometric derivation of Bragg s Law one assumes that crystalline Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this

More information

Model for Dopant and Impurity Segregation During Vapor Phase Growth

Model for Dopant and Impurity Segregation During Vapor Phase Growth Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Electromagnetic surface and line sources under coordinate transformations

Electromagnetic surface and line sources under coordinate transformations PHYSICAL REVIEW A 8, 3382 29 Electromagnetic surface an line sources uner coorinate transformations Steven A. Cummer,* Nathan Kuntz, an Bogan-Ioan Popa Department of Electrical an Computer Engineering

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Patterns in bistable resonant-tunneling structures

Patterns in bistable resonant-tunneling structures PHYSICAL REVIEW B VOLUME 56, NUMBER 20 Patterns in bistable resonant-tunneling structures 15 NOVEMBER 1997-II B. A. Glavin an V. A. Kochelap Institute of Semiconuctor Physics, Ukrainian Acaemy of Sciences,

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Light scattering and dissipative dynamics of many fermionic atoms in an optical lattice

Light scattering and dissipative dynamics of many fermionic atoms in an optical lattice Light scattering an issipative ynamics of many fermionic atoms in an optical lattice S. Sarkar, S. Langer, J. Schachenmayer,, an A. J. Daley, 3 Department of Physics an Astronomy, University of Pittsburgh,

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

A variational approach to Bogoliubov excitations and dynamics of dipolar Bose Einstein condensates

A variational approach to Bogoliubov excitations and dynamics of dipolar Bose Einstein condensates IOP PUBLISHING JOURNAL OF PHYSICS B: ATOMIC, MOLECULAR AND OPTICAL PHYSICS J. Phys. B: At. Mol. Opt. Phys. () 5 (pp) oi:./5-5///5 A variational approach to Bogoliubov excitations an ynamics of ipolar Bose

More information

Open System Dynamics with Non-Markovian Quantum Trajectories

Open System Dynamics with Non-Markovian Quantum Trajectories Open System Dynamics with Non-Markovian Quantum Trajectories W. T. Strunz 1,,L.Diósi,anN.Gisin 3 1 Fachbereich Physik, Universität GH Essen, 45117 Essen, Germany Research Institute for Particle an Nuclear

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Delocalization of boundary states in disordered topological insulators

Delocalization of boundary states in disordered topological insulators Journal of Physics A: Mathematical an Theoretical J. Phys. A: Math. Theor. 48 (05) FT0 (pp) oi:0.088/75-83/48//ft0 Fast Track Communication Delocalization of bounary states in isorere topological insulators

More information

1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a

1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a Theory of the nerson impurity moel: The Schrieer{Wol transformation re{examine Stefan K. Kehrein 1 an nreas Mielke 2 Institut fur Theoretische Physik, uprecht{karls{universitat, D{69120 Heielberg, F..

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

Kinetics of radiationless energy transfer from upper excited states

Kinetics of radiationless energy transfer from upper excited states 1 August 1997 Ž Chemical Physics Letters 74 1997 171 176 Kinetics of raiationless energy transfer from upper excite states EN Bounov 1, MN Berberan-Santos, JMG Martinho Centro e Quımica-Fısica Molecular,

More information

Physics 2212 GJ Quiz #4 Solutions Fall 2015

Physics 2212 GJ Quiz #4 Solutions Fall 2015 Physics 2212 GJ Quiz #4 Solutions Fall 215 I. (17 points) The magnetic fiel at point P ue to a current through the wire is 5. µt into the page. The curve portion of the wire is a semicircle of raius 2.

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

Laser Cooling and Trapping of Atoms

Laser Cooling and Trapping of Atoms Chapter 2 Laser Cooling and Trapping of Atoms Since its conception in 1975 [71, 72] laser cooling has revolutionized the field of atomic physics research, an achievement that has been recognized by the

More information

Entanglement is not very useful for estimating multiple phases

Entanglement is not very useful for estimating multiple phases PHYSICAL REVIEW A 70, 032310 (2004) Entanglement is not very useful for estimating multiple phases Manuel A. Ballester* Department of Mathematics, University of Utrecht, Box 80010, 3508 TA Utrecht, The

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Alpha Particle scattering

Alpha Particle scattering Introuction Alpha Particle scattering Revise Jan. 11, 014 In this lab you will stuy the interaction of α-particles ( 4 He) with matter, in particular energy loss an elastic scattering from a gol target

More information

Self-focusing and soliton formation in media with anisotropic nonlocal material response

Self-focusing and soliton formation in media with anisotropic nonlocal material response EUROPHYSICS LETTERS 20 November 1996 Europhys. Lett., 36 (6), pp. 419-424 (1996) Self-focusing an soliton formation in meia with anisotropic nonlocal material response A. A. Zoulya 1, D. Z. Anerson 1,

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

arxiv: v1 [physics.class-ph] 20 Dec 2017

arxiv: v1 [physics.class-ph] 20 Dec 2017 arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.

More information

Spectral Flow, the Magnus Force, and the. Josephson-Anderson Relation

Spectral Flow, the Magnus Force, and the. Josephson-Anderson Relation Spectral Flow, the Magnus Force, an the arxiv:con-mat/9602094v1 16 Feb 1996 Josephson-Anerson Relation P. Ao Department of Theoretical Physics Umeå University, S-901 87, Umeå, SWEDEN October 18, 2018 Abstract

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

arxiv: v1 [math-ph] 5 May 2014

arxiv: v1 [math-ph] 5 May 2014 DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic

More information

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9 Oregon State University PH 213 Spring Term 2018 Prep 1 Suggeste finish ate: Monay, April 9 The formats (type, length, scope) of these Prep problems have been purposely create to closely parallel those

More information

The Three-dimensional Schödinger Equation

The Three-dimensional Schödinger Equation The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation

More information

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010 NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Physics 115C Homework 4

Physics 115C Homework 4 Physics 115C Homework 4 Problem 1 a In the Heisenberg picture, the ynamical equation is the Heisenberg equation of motion: for any operator Q H, we have Q H = 1 t i [Q H,H]+ Q H t where the partial erivative

More information

Unified kinetic model of dopant segregation during vapor-phase growth

Unified kinetic model of dopant segregation during vapor-phase growth PHYSICAL REVIEW B 72, 195419 2005 Unifie kinetic moel of opant segregation uring vapor-phase growth Craig B. Arnol 1, * an Michael J. Aziz 2 1 Department of Mechanical an Aerospace Engineering an Princeton

More information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.

More information

Impurities in inelastic Maxwell models

Impurities in inelastic Maxwell models Impurities in inelastic Maxwell moels Vicente Garzó Departamento e Física, Universia e Extremaura, E-671-Baajoz, Spain Abstract. Transport properties of impurities immerse in a granular gas unergoing homogenous

More information

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract

More information

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating droplets tracking by holographic high speed video in turbulent flow Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es

More information

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012 arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,

More information