Atom-Atom Interactions in Ultracold Quantum Gases

Size: px
Start display at page:

Download "Atom-Atom Interactions in Ultracold Quantum Gases"

Transcription

1 Atom-Atom Interactions in Utracod Quantum Gases Caude Cohen-Tannoudji Lectures on Quantum Gases Institut Henri Poincaré, Paris, 5 Apri 7 Coège de France 1

2 Lecture 1 (5 Apri 7) Quantum description of eastic coisions between utracod atoms The basic ingredients for a mean-fied description of gaseous Bose Einstein condensates Lecture (7 Apri 7) Quantum theory of Feshbach resonances How to manipuate atom-atom interactions in a utracod quantum gas

3 A few genera references 1 L.Landau and E.Lifshitz, Quantum Mechanics, Pergamon, Oxford (1977) A.Messiah, Quantum Mechanics, North Hoand, Amsterdam (1961) 3 C.Cohen-Tannoudji, B.Diu and F.Laoë, Quantum Mechanics, Wiey, New Yor (1977) 4 C.Joachin, Quantum coision theory, North Hoand, Amsterdam (1983) 5 J.Daibard, in Bose Einstein Condensation in Atomic Gases, edited by M.Inguscio, S.Stringari and C.Wieman, Internationa Schoo of Physics Enrico Fermi, IOS Press, Amsterdam, (1999) 6 Y. Castin, in Coherent atomic matter waves, Lecture Notes of Les Houches Summer Schoo, edited by R. Kaiser, C. Westbroo, and F. David, EDP Sciences and Springer-Verag (1) 7 C.Cohen-Tannoudji, Cours au Coège de France, Année C.Cohen-Tannoudji, Compéments de mécanique quantique, Cours de 3 ème cyce, Notes de cours rédigées par S.Haroche 9 T.Köher, K.Gora, P.Juienne, Rev.Mod.Phys. 78, (6) 3

4 1 - Introduction Outine of ecture 1 - Scattering by a potentia. A brief reminder Integra equation for the wave function Asymptotic behavior. Scattering ampitude Born approximation 3 - Centra potentia. Partia wave expansion Case of a free partice Effect of the potentia. Phase shifts S-Matrix in the anguar momentum representation 4 - Low energy imit Scattering ength a Long range effective interactions and sign of a 5 - Mode used for the potentia. Pseudo-potentia Motivation Determination of the pseudo-potentia Scattering and bound states of the pseudo-potentia Pseudo-potentia and Born approximation 4

5 Interactions between utracod atoms At ow densities, -body interactions are predominant and can be described in terms of coisions. We wi focus here on eastic coisions (athough ineastic coisions and 3-body coisions are aso important because they imit the achievabe spatia densities of atoms). Coisions are essentia for reaching therma equiibrium At very ow temperatures, mean-fied descriptions of degenerate quantum gases depend ony on a very sma number of coisiona parameters. For exampe, the shape and the dynamics of Bose Einstein condensates depend ony on the scattering ength Possibiity to contro atom-atom interactions with Feshbach resonances. This expains the increasing importance of utracod atomic gases as simpe modes for a better understanding of quantum many body systems Purpose of these ectures: Present a brief review of the concepts of atomic and moecuar physics which are needed for a quantitative description of interactions in utracod atomic gases. 5

6 Notation Two atoms, with r r 1 (mass )V m, interacting with a -body interaction potentia In ecture 1, we ignore the spins degrees of freedom. They wi be taen into account in ecture. Hamitonian p (1 )p H = + + Vr r (1.1) 1 m m Change of variabes R = ( r + r ) P = p + p Center of mass variabes 1 /G G 1 M = m + m = m Tota mass 1 /= r r p = ( p p ) Reative variabes 1 /1 /r μ = mm ( m + m ) = m Reduced mass 1 1 ( = + μ ) H P M p + V ( r) G (1.) /Hamitonian of a free partice with mass M H /CM H re Hamitonian of a fictitious partice with mass μ, moving in V(r) 6

7 Finite range potentia Simpe case where V ( r) = for r > b b is caed the range of the potentia One can extend the resuts obtained in this simpe case to potentias decreasing fast enough with r at arge distances. For exampe, for the Van der Waas interactions between atoms decreasing as C 6 / r 6 for arge r, one can define an effective range See for exampe Ref. 5 b VdW μ /C 6 = 1 4 (1.3) 7

8 1 - Introduction Outine of ecture 1 - Scattering by a potentia. A brief reminder Integra equation for the wave function Asymptotic behavior. Scattering ampitude Born approximation 3 - Centra potentia. Partia wave expansion Case of a free partice Effect of the potentia. Phase shifts S-Matrix in the anguar momentum representation 4 - Low energy imit Scattering ength a Long range effective interactions and sign of a 5 - Mode used for the potentia. Pseudo-potentia Motivation Determination of the pseudo-potentia Scattering and bound states of the pseudo-potentia 8

9 Scattering by a potentia. A brief reminder Shrödinger equation for the reative partice (with E>) Δ+ V ( r) ψ( r) = Eψ( r) E = μ μ μ Δ+ ψ ( ) = ( ) ψ ( ) r V r r (1.4) Green function of Δ + Δ+ G ( r) = δ ( r) The boundary conditions for G wi be chosen ater on (1.5) (1.6) Integra equation for the soution of Shrödinger equation ϕ ( ) r Soution of the equation without the right member : Δ+ ϕ ( ) r = (1.7) ( ) ( ) 3 ψ r = ϕ r + d r G ( r r ) V ( r ) ψ ( r ) (1.8) 9

10 Choice of boundary conditions We choose for ϕ a pane wave.with wave vector ϕ ( ) r i.= e = e κ r /i r κ = (1.9) and we choose, for the Green function G, boundary conditions corresponding to an outgoing spherica wave (see Ref., Chap.XIX) 1 e G + ( r r ) = 4π r ir r r (1.1) We thus get the foowing soution for Schrödinger equation r r (+ )1 3 e ( )(+ )i ir ψ r = e d ψ 4π r.v r r r r If V has a finite range b, the integra over r is restricted to a finite range and we can write: If r b, r r r r. n with n= r / r ir + i. r e ψ ( r) e f(, κ, n) r μ 3 - in. r + f(, κ, n) = ψ 4π d r e V( r ) ( r ) (1.11) (1.1) 1

11 Scattering state with an outgoing spherica wave Asymptotic behavior for arge r + The state ψ )r is a soution of the Schrodinger equation behaving exp.for arge r as the sum of an incoming pane wave ( ir) and of (,, ()exp()/ an outgoing spherica wave f κ n i r r Scattering ampitude f (,, ).μ (3 )+ ) r (i f κ n = d e ψ 4π r V r r -/We have pu t = n = r,, )r( κ n is the ampitude of /f the outgoing spherica wave in the direction of = n = r r. It depends ony on and on the poar anges θ and ϕ of with respect to (1.13) Differentia cross section Comparing the fuxes aong and, one gets : (,, / dσ dω = κ n )f (1.14) 11

12 Born approximation In the scattering ampitude, the potentia V appears expicity (,, )-.μ (3 )( r + )i f κ n = d e ψ 4π r V r )r(to owest order in V, one can thus repace ψ + r by the zeroth exp(. ) order soution of the Schrodinger equation ir (,, )().μ (3 r )i f κ n = d e 4π r V r This is the Born approximation In this approximation, the scattering ampitude is proportiona to the spatia Fourier transform of the potentia (1.15) (1.16) 1

13 Low energy imit The presence of V(r ) in the scattering ampitude (,,)-. μ ()(3 r + f κ n = r V r ψ r 4π d e restricts the integra over r to a finite range r < b i. r If b 1, one can repace e by 1. The scattering ampitude μ 3 + f(, κ, n) = d ψ 4π r V( r ) ( r ) then no onger depends on the direction of the scattering vector It is sphericay symmetric even if V( r) is not.,(,, )When f κ n a (.e ψ + ) ir ir a r e a 1 r r a is a constant, caed scattering ength, which wi be discussed in more detais ater on )i (1.17) (1.18). (1.19) 13

14 Another interpretation of the outgoing scattering state Another expression for this state (see refs. 4 and 8) ψ = ϕ + Lim V ψ T = p μ ε / + E T + iε + For ε non zero but very sma, ψ appears as the state obtained at t = by starting from the free state ϕ at t = and by switching on sowy V on a time interva on the order of ε Ingoing scattering state ( )ir 1 ( ) ( ψ r = e d ψ 4π r V r r r r 1 ψ = ϕ + Lim V ψ ε + E T iε If one starts from such a state at t = and if one switches off V sowy on a time interva on the order of ε, one gets the free state ϕ at t = + /3 e. r r )i / (1.) (1.1) 14

15 Definition S - Matrix (,)S = ϕ S ϕ = Lim ϕ U t t ϕ ji 1,j i t1 j i t + : U evoution operator in interaction representation + One can show that = ψ ψ Quaitative interpretation S ji j i (1.) (1.3) V is switched on sowy (time scae ħ/ε) between - and, and then switched off sowy (time scae ħ/ε) between and + One starts from ϕ i at t = - and one oos for the probabiity ampitude to be in ϕ j at t = + From t.= to t = the initia free state ϕ is transformed i + into ψ Since the evoution operator is unitary, and since ψ i j, + transforms into ϕ from t = to t = +,ψ ψ is the j j i ampitude to find the system in the free state ϕ at = + if one starts from ϕ at =.t j t i 15

16 1 - Introduction Outine of ecture 1 - Scattering by a potentia. A brief reminder Integra equation for the wave function Asymptotic behavior. Scattering ampitude Born approximation 3 - Centra potentia. Partia wave expansion Case of a free partice Effect of the potentia. Phase shifts S-Matrix in the anguar momentum representation 4 - Low energy imit Scattering ength a Long range effective interactions and sign of a 5 - Mode used for the potentia. Pseudo-potentia Motivation Determination of the pseudo-potentia Scattering and bound states of the pseudo-potentia 16

17 Centra potentia V depends ony on r 1D radia Schrödinger equation One oos for soutions of the form ϕ ( r) = R ( r) Y ( n) n= r / r If we put m m u( r) R( r) = r with the boundary condition u ( ) = (1.4) (1.5) (1.6) one gets for u the foowing 1D radia equation ) d + 1 μ () ( + V r = u r ) dr r ( 1D Schrödinger equation for a partice moving in a potentia which is the sum of V and of the centrifuga barrier (+ 1 ) μ r (1.7) (1.8) 17

18 Case of a free partice (V=) The soutions of the Schrödinger equation are: )( ) ()( (ϕ r = j r Y n m m ) π where the j are the spherica )Besse functions of order ()r ()1 sin(π )j r ()!j r (r r + 1 r r For arge r we thus have sin(/))( ) ( )r π ϕ r Y n m r m π (/r )(/)i r π i r π ( () = e e Y n m π ir = outgoing spherica wave + ingoing spherica wave These functions form )an ()orthonorma set (see Appendix) ϕ ϕ = δ ( δ δ ( m m ) mm (1.9) (1.3) (1.31) (1.3) 18

19 Expansion of a pane wave in free spherica waves Pane wave ()/. ( r = π e = δ 3 r )i (1.33) The factor (π) -3/ is introduced for the orthonormaization One can show that: ().()/* =+ ()( 3 3 )( )()ir π e = π 4π i Y n Y κ j r m m = m= *() m=+ 1 ()()( ) = i Y κ ϕ r m m = m= / /m with κ = (1.34) ({ () ϕm } The transformation from the orthonorma basis { } he orthonorma basis is given by the matrix ) ()* 1 ( ϕ = δ Y κ m m ) to t (1.35) 19

20 Effect of a potentia. Phase shifts We come bac to the Schrödinger equation with V. Consider, for r arge, an incoming wave exp[-i(r-π/)]. Since the refection coefficient of V is 1 (conservation of the norm), the refected outgoing wave has the same moduus and has just accumuated a phase shift with respect to the V= case. The superposition of the waves is thus a shifted sinusoid. We concude that there is a set of soutions of the Schrödinger equation with V which behave for arge r as: sin/()( ) ( ) π + δ ϕ r r Y n (1.36) m r m π r One can show that these functions are orthonormaized (see Appendix) ϕ ϕ = δ ( δ δ m m ) mm (1.37) They don t form a basis if there are aso bound states in the potentia V.

21 Partia wave expansion of the outgoing scattering state Consider the inear superposition of the states ϕ with the same m coefficients as those appearing in the expansion (1.34) of the pane wave exp(on ).the, each state being mutipied by the phase factor iδ We wi show that such a state is nothing but the outgoing state ψ (mutipied by for having an orthonormaized state) ) ϕ(m (/+ 3 π m=+ + 1 ( = = m ) = )( ) iδ ψ i Y κ e ϕ = m=+ = m= ϕ () m Before demonstrating this identity, et us discuss its physica meaning.the outgoing scattering state is obtained by switching ( on sowy V on the free state. Each spherica wave ϕ of the )m expansion of is transformed into ϕ, but in addition it acquires m iδ a phase factor e which depends on and which thus varies from one spherica wave of the expansion to another one.m * e iδ ϕ m m (1.38) 1

22 Demonstration For arge r, the inear superposition introduced in (1.38) behaves as: ( π + δ ) ( π ) m= + *i r 1 ()( ) ( ) i r e e i Y κ Y n m m = m= π ir exp/exp/ Using i r i r we get: ( δ 1) / =+ * sin(/)i i m 1 ()( ) ( ) π κ r i Y Y n + m m = m= π r i r ( ) ( ) ( δ π δ π i + = 1 + e 1) // π e e ir e (1.39) The contribution of the first term of the bracet is nothing but the asymptotic expansion of the pane wave in spherica waves. The second term gives an outgoing spherica wave 3.ir ()ir (,, /)e π e + f κ n (1.4) r This demonstrates that the state given in (1.38) is an outgoing scattering state and gives in addition the expression of the ampitude f (1.4) (1.41)

23 Scattering ampitude in terms of the phase shifts (, =+, )m 4π sin(δ )( *)i f κ n = e δ Y n Y κ m m = m= 1 ()sin(cos(1.43) ) i δ = + 1 e δ P θ (cos) = where P θ is a Legendre poynomia and where θ. is the ange between n and κ Integrating f over the poar anges of κ gives the scattering cross section ) ()()4π ()sin( σ = (σ σ = + 1 δ (1.44) ) = Scattering of identica partices Quantum interference between different paths )()()θ π - θ f θ (f θ + ε f π θ ) ε =+ 1 1 for bosons (fermions)(π sin()( σ = π θ dθ θ + ε π θ tota /f f (1.45) ) 3

24 Partia wave expansion of the ingoing scattering state ψ is given by a inear superposition of the states ϕ anaogous m + to the one introduced for ψ each state being now mutipied by exp() exp().the phase factor iδ instead of iδ () =+ =+ 1 * δ δ ψ = m ()() i κ e ϕ = m i i Y ϕ e ϕ m m m m = m= = m= The demonstration of this identity is simiar to the one given above for the outgoing scattering state (1.46) If we start from this state and if we switch off V sowy, it transforms ( into the free state. Each wave ϕ is transformed into ϕ, but m m iδ ) in addition its phase factor changes from e to 1 w hich corresponds +iδ to acquiring a phase factor e Finay, when we go from t = to t = + switching on and then )() switching off V sowy, we start from ϕ and we end in ϕ acquiring m m + iδ + iδ + iδ, a goba phase factor e e = e(4

25 S Matrix in the anguar momentum representation S = S = ψ ψ + We use the expansion of ψ + and ψ in spherica waves m=+ ( m =+ + ψ ϕ ) i e δ ϕ ψ ϕ ϕ = = m= m m = m= = m = = = m = ) iδ e( m m This gives a first expression of S m=+ m =+ + S = ψ ψ = ϕ ϕ ϕ ϕ (On the other hand, a change of basis gives for S = m= = m = )() + iδ iδ e e m m m m ( ) δ δ δ mm =+ m =+ (S S ϕ ϕ S ϕ ϕ = = Comparing the expressions obtained for we get, S )() + i δ () ϕ S ϕ = e(δ δ δ m m mm (1.48) (1.49) ))m () () ( m m m m which shows that the S-matrix is diagona exin p(the anguar )momentum representation, with diagona eements iδ ceary showing the unitarity of S (1.47) (1.5) (1.51) 5

26 1 - Introduction Outine of ecture 1 - Scattering by a potentia. A brief reminder Integra equation for the wave function Asymptotic behavior. Scattering ampitude Born approximation 3 - Centra potentia. Partia wave expansion Case of a free partice Effect of the potentia. Phase shifts S-Matrix in the anguar momentum representation 4 - Low energy imit Scattering ength a Long range effective interactions and sign of a 5 - Mode used for the potentia. Pseudo-potentia Motivation Determination of the pseudo-potentia Scattering and bound states of the pseudo-potentia 6

27 Centra potentia. Low energy imit Suppose first V=. The centrifuga barrier in the 1D Schrödinger equation prevents the partice from approaching near the region r= ) ) (= μr μ μ r ()/ r = + 1 μ r r () r = + ( 1 db If the range b of the potentia is sma enough, i.e. if.b a partice with cannot fee the potentia Ony = wave wi fee V " s-wave scattering" d μ () ( + V r = u r ) dr db (1.5) (1.53) 7

28 Scattering ength sin()for r arge enough, u = u varies as r + δ sin (). Let v be the function r + δ extending for a u r Let P be the intersection point of v with the r-axis which is the cosest from the origin.by definition, the scattering ength a is the imit of the abscissa of P when (see figure) Expansion of v in powers of r near r = v r r r Abscissa of : a ) sin()sin()cos(= + δ δ δ + r 1 tan()p δ tan()im δ π () π = δ (+ ) (1.54) (1.55) 8

29 Limit = Scattering ength (continued) d dr μ () ( r u r )V = Far from r=o, the soution of the S.E. is a straight ine and )v r r (a (1.57) The abscissa of Q is equa to a (1.56) Scattering cross section ()4π ()sin()()σ = δ σ π + 1 = 4 = ()() δ a σ = 4πa = () For identica bosons σ = 8π a = sin(δ ) (1.58) (1.59) (1.6) 9

30 Scattering ength for square potentias Square potentia barriers()v r Square barrier of height /V and width b = μ V a )u (r b )v (r r a r For r)> b (and = r = v r )u r (a For ()< b and ()r = u r = u r The curvature of (u is )positive and u r = = We concude that the scattering ength is aways positive and smaer than the range b of the potentia a b When V (hard sphere potentia) a b 3

31 For and (Scattering ength for square potentias (continued) Square potentia wes V r a V )v (r r a u r b r = /V μ r)r )u > b ( = = v r r a For )< b and ()= u r = ( u r The curvature of (u is )r negative and u r = = If V is sma enough so that there is no bound state in the potentia we, the curvature of u for r < b is sma and a is negative When increases, the curvature of u for r < bincreases in absoute vaue and a Then a switches suddeny to +.V and decreases. This divergence of a corresponds to the appearance of the first bound state in the potentia we 31

32 Square potentia wes (continued) Variations of a with figure taen from Ref.5 When the depth of the potentia we increases, divergences )of occur for a vaues of V such that b = n +(1π /a corresponding to the appearances of successive bound states in the potentia we. These divergences of a which goes from to + are caed "zero-energy" resonances 3

33 Long range effective interactions and sign of a The scattering ength determines how the ong range behavior of the wave functions is modified by the interactions. To understand how the sign of a is reated to the sign of the effective ong range interactions, it wi be usefu to consider the partice encosed in a spherica box with radius R, so that we have the boundary condition u R )=( (1.61) eading to a discrete energy spectrum In the absence of interactions (V=), the normaized eigenstates and the eigenvaues of the 1D Schrödinger equation are: sin(/))() 1 π π ψ = N r R 1 π = (N μ =,,r E N N R r R (/N r R )(1.6)..N inπsfigure corresponding to N=3 R r 33

34 V = R V a > a The dotted ine is the sinusoid outside the range of the potentia It has a shorter waveength than for V=, and thus a arger wave number. The inetic energy in this region, which is aso the tota energy, is arger R V a < a R The dotted ine is the sinusoid outside the range of the potentia It has a onger waveength than for V=, and thus a smaer wave number. The inetic energy in this region, which is aso the tota energy, is smaer 34

35 Correction to the energy to first order in a a R )/ For the state ψ, we have R = N λ N Because of the interactions, these N haf waveengths occupy now a ength /R a so that ()/λ = R /N λ = R/ a N ) = π λ = π λ = R R a ///) a E = μ E = E = E R R a (E 1 + N N N (N N R a π Finay, we have δ = = = N E E E E a (1.63) N N N N 3 R μ R Long range effective interactions are - repusive if a > /(- attractive if a < 35

36 1 - Introduction Outine of ecture 1 - Scattering by a potentia. A brief reminder Integra equation for the wave function Asymptotic behavior. Scattering ampitude Born approximation 3 - Centra potentia. Partia wave expansion Case of a free partice Effect of the potentia. Phase shifts S-Matrix in the anguar momentum representation 4 - Low energy imit Scattering ength a Long range effective interactions and sign of a 5 - Mode used for the potentia. Pseudo-potentia Motivation Determination of the pseudo-potentia Scattering and bound states of the pseudo-potentia 36

37 Mode used for the potentia V(r) Why not using the exact potentia? The interaction potentia is very difficut to cacuate exacty. A sma error in V can introduce a very arge error on the scattering ength deduced from this potentia. Mean fied description of utracod quantum gases require in genera a first order treatment of the effect of V (Born approximation). But Born approximation cannot be in genera appied to the exact potentia Approach foowed here The motivation here is not to cacuate the scattering ength. This parameter is supposed nown experimentay. We are interested in the derivation of the macroscopic properties of the gas from a mean fied description using a singe parameter which is a. The ey idea is to repace the exact potentia by a pseudo-potentia simper to use than the exact one and obeying conditions: - It has the same scattering ength as the exact potentia - It can be treated with Born approximation so that mean fied descriptions of its effects are possibe 37

38 Determination of the pseudo-potentia Derivation à a H.Bethe and R.Peiers (Y.Castin, private communication) We add to the 3D Schrödinger equation of a free partice (V=) a term proportiona to a deta function Δ μ ψ )+ δ ( ) (= μ ψ ( ) r C r r (1.64) To determine the coefficient, we impose to the soution of this equation to coincide with the extension to a r of the asymptotic in()/()/c behavior r + δ of the true wave function sr u r r In particuar, for sma enough, one shoud have: )()ψ r (B r a r (1.65) /r )(/)()Inserting (1.65) into 1.64 and using Δ 1 = 4π δ r,r we get an equation containing a deta function mutipied by a coefficient 4π (C a B μ which must vanish. This gives the coefficient C appearing in (1. 64) 4π 4π C = g B where g = a = a μ m (1.66) 38

39 Determination of the pseudo-potentia (continued) It wi be more convenient to express C=gB, not in terms of the coefficient B appearing in the wave function ψ = B(r-a)/r of equation (1.65), but in terms of the wave function ψ itsef. We use for that d B = rψ d r r = Equation (1) can be rewritten as: Δ ψ( r) + Vpseudoψ( r) = ψ( r) μ μ d where Vpseudoψ( r) = gδ( r) [ rψ( r) ] dr (/)is caed the pseudo-potentia. The term d d reguarizes pseudo r r ( ) /the action.v of δ r when it acts on functions behaving as 1 r near r = For functions which are reguar in r =, V has the pseudo ):same (effect as δ )()/r ()()()( )ψ r = u r r with u V ψ r = g u δ r (pseudo )()()( ()g ψ r reguar in r = V ψ r = gψ δ r pseudo (1.67) (1.68) (1.69) (1.7) 39

40 Scattering states of the pseudo-potentia We are ooing for soutions of equation (1.68) with E > For the centrifuga,barrier prevents the (partice ),from approaching r = and one (can show that ψ = so that, ). according to (1.7), V ψ r = V gives ony s-scattering pseudo pseudo,and one can write ψ ( r ) = ()/where ( )r r u can be.u )()()() r u u r u ()()u 1 d u Δ = Δ + = 4π u δ r + r r r (r dr The Schrödinger equation for u becomes: ()() ()( )u r ( )() 4π δ δ u r u r + + g r u = μ r μ r ( )Canceing (the term proportiona to δ and the term independant ),r of δ r we get equations: ) μ ()() ()u = g = a u r + u r = u π ( (1.71) (1.7) (1.73) 4

41 Scattering states of the pseudo-potentia (continued) The soution of the second equation can be written: )sir = nu r +(δ ( ) Inserting this soution into the first equation gives: anδ = ta On the other hand, the s-wave scattering ampitude is equa to: )1 siδ ni f =( e δ Using equation (1.75) giving tanδ finay gives after simpe agebra: )a f = 1 +(ia is proportiona to a A first order treatment of V pseudo pseudo thus gives the correct resut for ).V the scattering ampitude in the zero energy imit ( a =. This shows that Born approximation can be used with V for utracod atoms pseudo The conditions imposed above on V pseudo are thus fufied (1.74) (1.75) (1.76) (1.77) 41

42 The unitary imit From the expression of the scattering ampitude obtained above, we deduce the scattering ampitude for identica bosons ) 8π a σ = 1 +(a which is vaid for a. )The ow energy imit a 1 gives the we nown resut: ( σ( ) a 1 8πa There is another interesting imit, corresponding to high energy, )or strong interaction a 1 eading to resut independent of a: ( 8π σ ( ) a 1 This is the so caed unitary imit (1.78) (1.79) (1.8) 4

43 Bound state of the pseudo-potentia The cacuation is the same as for the /scattering states, except that we repace the positive energy μ by a negative one - κ / μ The equations derived from the Schrodinger equation are now: ) ()() ()u = a u r κ u r = (1.81) u ( The soution of the second equation (finite for r ) is: u r )=(e κ r which inserted into the first equation gives: κ =/1 a The pseudo-potentia thus has a bound state with an energy E = μa and a wave function: xp e r a (1.8) (1.83) (1.84) (1.85) 43

44 Energy shifts produced by the pseudo-potentia We come bac to the probem of a partice in a box of radius R. We have cacuated above the energy shifts of the discrete energy eves of this partice produced by a potentia characterized by a scattering ength a. To first order in a, we found: π δ = N E a N 3 μ R This resut was deduced directy from the modification induced by the interaction on the asymptotic behavior of the wave functions and not from a perturbative treatment of V. We show now that: ()( δe = ψ V ψ to first order in V N pseudo )N N pseudo (1.86) (1.87) which is another evidence for the fact that the effect of the pseudo potentia can be cacuated perturbativey, which is not the case for the rea potentia. For exampe, a hard core potentia (V= for r< a) cannot obviousy be treated perturbativey, but its scattering ength is a, and using a pseudo-potentia with scattering ength a aows perturbative cacuations. 44

45 Demonstration The unperturbed normaized eigenfunctions of the partice in the spherica box are: )sin(/)( 1 Nπ r R ψ r =()N )π R r () ψ r is reguar in r (= ) and ()N 1 1 ψ =(Nπ N 3 π R (1.47) (1.88) so that (()( )(V ψ r = gδ r ψ )() N N )pseudo (1.89) We deduce that ()( N π π δ = ψ = = N E g g a N )N 3 3 πr μr (1.9) which coincides with the resut obtained above in (1.63). 45

46 Concusion Eastic coisions between utracod atoms are entirey characterized by a singe number, the scattering ength Effective ong distance interactions are attractive if a< and repusive if a> Giving the same scattering ength as the rea potentia, the pseudo-potentia gives the good asymptotic behavior for the wave function describing the reative motion of atoms, and thus correcty describes their ong distance interactions In a diute gas, atoms are far apart. The pseudo-potentia is proportiona to a and can be treated perturbativey. A first order treatment of the pseudo-potentia is the basis of mean fied description of Bose Einstein condensates where each atom moves in the mean fied produced by a other atoms. Next step: can one change the scattering ength? 46

Atom-atom interactions in ultracold gases

Atom-atom interactions in ultracold gases Atom-atom interactions in ultracold gases Claude Cohen-Tannoudji To cite this version: Claude Cohen-Tannoudji. Atom-atom interactions in ultracold gases. DEA. Institut Henri Poincaré, 5 et 7 Avril 7, 7.

More information

1.2 Partial Wave Analysis

1.2 Partial Wave Analysis February, 205 Lecture X.2 Partia Wave Anaysis We have described scattering in terms of an incoming pane wave, a momentum eigenet, and and outgoing spherica wave, aso with definite momentum. We now consider

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation 22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements

More information

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p

More information

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18 Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms

More information

221B Lecture Notes Notes on Spherical Bessel Functions

221B Lecture Notes Notes on Spherical Bessel Functions Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,

More information

Elements of Kinetic Theory

Elements of Kinetic Theory Eements of Kinetic Theory Statistica mechanics Genera description computation of macroscopic quantities Equiibrium: Detaied Baance/ Equipartition Fuctuations Diffusion Mean free path Brownian motion Diffusion

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Physics 506 Winter 2006 Homework Assignment #6 Solutions

Physics 506 Winter 2006 Homework Assignment #6 Solutions Physics 506 Winter 006 Homework Assignment #6 Soutions Textbook probems: Ch. 10: 10., 10.3, 10.7, 10.10 10. Eectromagnetic radiation with eiptic poarization, described (in the notation of Section 7. by

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

Formulas for Angular-Momentum Barrier Factors Version II

Formulas for Angular-Momentum Barrier Factors Version II BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A

More information

Midterm 2 Review. Drew Rollins

Midterm 2 Review. Drew Rollins Midterm 2 Review Drew Roins 1 Centra Potentias and Spherica Coordinates 1.1 separation of variabes Soving centra force probems in physics (physica systems described by two objects with a force between

More information

Scattering of Particles by Potentials

Scattering of Particles by Potentials Scattering of Partices by Potentias 2 Introduction The three prototypes of spectra of sef-adjoint operators, the discrete spectrum, with or without degeneracy, the continuous spectrum, andthe mixed spectrum,

More information

hole h vs. e configurations: l l for N > 2 l + 1 J = H as example of localization, delocalization, tunneling ikx k

hole h vs. e configurations: l l for N > 2 l + 1 J = H as example of localization, delocalization, tunneling ikx k Infinite 1-D Lattice CTDL, pages 1156-1168 37-1 LAST TIME: ( ) ( ) + N + 1 N hoe h vs. e configurations: for N > + 1 e rij unchanged ζ( NLS) ζ( NLS) [ ζn unchanged ] Hund s 3rd Rue (Lowest L - S term of

More information

Elements of Kinetic Theory

Elements of Kinetic Theory Eements of Kinetic Theory Diffusion Mean free path rownian motion Diffusion against a density gradient Drift in a fied Einstein equation aance between diffusion and drift Einstein reation Constancy of

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

Physics 566: Quantum Optics Quantization of the Electromagnetic Field

Physics 566: Quantum Optics Quantization of the Electromagnetic Field Physics 566: Quantum Optics Quantization of the Eectromagnetic Fied Maxwe's Equations and Gauge invariance In ecture we earned how to quantize a one dimensiona scaar fied corresponding to vibrations on

More information

Joel Broida UCSD Fall 2009 Phys 130B QM II. Homework Set 2 DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS.

Joel Broida UCSD Fall 2009 Phys 130B QM II. Homework Set 2 DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS. Joe Broida UCSD Fa 009 Phys 30B QM II Homework Set DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS. You may need to use one or more of these: Y 0 0 = 4π Y 0 = 3 4π cos Y

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated

More information

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n. Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

Lecture 17 - The Secrets we have Swept Under the Rug

Lecture 17 - The Secrets we have Swept Under the Rug Lecture 17 - The Secrets we have Swept Under the Rug Today s ectures examines some of the uirky features of eectrostatics that we have negected up unti this point A Puzze... Let s go back to the basics

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The

More information

Lecture 8 February 18, 2010

Lecture 8 February 18, 2010 Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some

More information

Parallel-Axis Theorem

Parallel-Axis Theorem Parae-Axis Theorem In the previous exampes, the axis of rotation coincided with the axis of symmetry of the object For an arbitrary axis, the paraeaxis theorem often simpifies cacuations The theorem states

More information

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential Lecture 6 Povh Krane Enge Wiiams Properties of -nuceon potentia 16.1 4.4 3.6 9.9 Meson Theory of Nucear potentia 4.5 3.11 9.10 I recommend Eisberg and Resnik notes as distributed Probems, Lecture 6 1 Consider

More information

Elements of Kinetic Theory

Elements of Kinetic Theory Eements of Kinetic Theory Statistica mechanics Genera description computation of macroscopic quantities Equiibrium: Detaied Baance/ Equipartition Fuctuations Diffusion Mean free path Brownian motion Diffusion

More information

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Physics 505 Fall Homework Assignment #4 Solutions

Physics 505 Fall Homework Assignment #4 Solutions Physics 505 Fa 2005 Homework Assignment #4 Soutions Textbook probems: Ch. 3: 3.4, 3.6, 3.9, 3.0 3.4 The surface of a hoow conducting sphere of inner radius a is divided into an even number of equa segments

More information

arxiv: v1 [hep-lat] 23 Nov 2017

arxiv: v1 [hep-lat] 23 Nov 2017 arxiv:1711.08830v1 [hep-at] 23 Nov 2017 Tetraquark resonances computed with static attice QCD potentias and scattering theory Pedro Bicudo 1,, Marco Cardoso 1, Antje Peters 2, Martin Pfaumer 2, and Marc

More information

Agenda Administrative Matters Atomic Physics Molecules

Agenda Administrative Matters Atomic Physics Molecules Fromm Institute for Lifeong Learning University of San Francisco Modern Physics for Frommies IV The Universe - Sma to Large Lecture 3 Agenda Administrative Matters Atomic Physics Moecues Administrative

More information

Sound-Particles and Phonons with Spin 1

Sound-Particles and Phonons with Spin 1 January, PROGRESS IN PHYSICS oume Sound-Partices Phonons with Spin ahan Minasyan aentin Samoiov Scientific Center of Appied Research, JINR, Dubna, 498, Russia E-mais: mvahan@scar.jinr.ru; scar@off-serv.jinr.ru

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations .615, MHD Theory of Fusion ystems Prof. Freidberg Lecture : The Moment Equations Botzmann-Maxwe Equations 1. Reca that the genera couped Botzmann-Maxwe equations can be written as f q + v + E + v B f =

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),

More information

Math 124B January 17, 2012

Math 124B January 17, 2012 Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia

More information

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes... : Soution Set Seven Northwestern University, Cassica Mechanics Cassica Mechanics, Third Ed.- Godstein November 8, 25 Contents Godstein 5.8. 2. Components of Torque Aong Principa Axes.......................

More information

Nuclear Size and Density

Nuclear Size and Density Nucear Size and Density How does the imited range of the nucear force affect the size and density of the nucei? Assume a Vecro ba mode, each having radius r, voume V = 4/3π r 3. Then the voume of the entire

More information

17 Lecture 17: Recombination and Dark Matter Production

17 Lecture 17: Recombination and Dark Matter Production PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was

More information

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL

More information

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

More information

arxiv:gr-qc/ v2 10 Apr 1997

arxiv:gr-qc/ v2 10 Apr 1997 AN ALGEBRAIC INTERPRETATION OF THE WHEELER-DEWITT EQUATION arxiv:gr-qc/960900v 0 Apr 997 John W. Barrett Louis Crane 9 Apri 997 Abstract. We make a direct connection between the construction of three dimensiona

More information

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS

More information

arxiv:hep-ph/ v1 26 Jun 1996

arxiv:hep-ph/ v1 26 Jun 1996 Quantum Subcritica Bubbes UTAP-34 OCHA-PP-80 RESCEU-1/96 June 1996 Tomoko Uesugi and Masahiro Morikawa Department of Physics, Ochanomizu University, Tokyo 11, Japan arxiv:hep-ph/9606439v1 6 Jun 1996 Tetsuya

More information

arxiv:gr-qc/ v1 12 Sep 1996

arxiv:gr-qc/ v1 12 Sep 1996 AN ALGEBRAIC INTERPRETATION OF THE WHEELER-DEWITT EQUATION arxiv:gr-qc/960900v Sep 996 John W. Barrett Louis Crane 6 March 008 Abstract. We make a direct connection between the construction of three dimensiona

More information

David Eigen. MA112 Final Paper. May 10, 2002

David Eigen. MA112 Final Paper. May 10, 2002 David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.

More information

Cluster modelling. Collisions. Stellar Dynamics & Structure of Galaxies handout #2. Just as a self-gravitating collection of objects.

Cluster modelling. Collisions. Stellar Dynamics & Structure of Galaxies handout #2. Just as a self-gravitating collection of objects. Stear Dynamics & Structure of Gaaxies handout # Custer modeing Just as a sef-gravitating coection of objects. Coisions Do we have to worry about coisions? Gobuar custers ook densest, so obtain a rough

More information

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions Physics 116C Hemhotz s an Lapace s Equations in Spherica Poar Coorinates: Spherica Harmonics an Spherica Besse Functions Peter Young Date: October 28, 2013) I. HELMHOLTZ S EQUATION As iscusse in cass,

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary

More information

Charge Density from X-ray Diffraction. Methodology

Charge Density from X-ray Diffraction. Methodology Charge Density from X-ray Diffraction. Methodoogy Ignasi Mata imata@icmab.es Master on Crystaography and Crystaization, 2012 Outine I. Charge density in crystas II. The mutipoar refinement III. Methodoogy

More information

Electron-impact ionization of diatomic molecules using a configuration-average distorted-wave method

Electron-impact ionization of diatomic molecules using a configuration-average distorted-wave method PHYSICAL REVIEW A 76, 12714 27 Eectron-impact ionization of diatomic moecues using a configuration-average distorted-wave method M. S. Pindzoa and F. Robicheaux Department of Physics, Auburn University,

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

Lecture VIII : The pseudopotential

Lecture VIII : The pseudopotential Lecture VIII : The pseudopotentia I. KOHN-SHAM PROBLEM FOR AN ISOLATED ATOM For a one-eectron atom, the Couombic potentia, V ( r) = V (r) = Z/r is sphericay symmetric. The soutions may then be spit into

More information

Chemical Kinetics Part 2

Chemical Kinetics Part 2 Integrated Rate Laws Chemica Kinetics Part 2 The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates the rate

More information

AAPT UNITED STATES PHYSICS TEAM AIP 2012 DO NOT DISTRIBUTE THIS PAGE

AAPT UNITED STATES PHYSICS TEAM AIP 2012 DO NOT DISTRIBUTE THIS PAGE 2012 Semifina Exam 1 AAPT UNITED STATES PHYSICS TEAM AIP 2012 Semifina Exam DO NOT DISTRIBUTE THIS PAGE Important Instructions for the Exam Supervisor This examination consists of two parts. Part A has

More information

Methods for Ordinary Differential Equations. Jacob White

Methods for Ordinary Differential Equations. Jacob White Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa

More information

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,

More information

Candidate Number. General Certificate of Education Advanced Level Examination January 2012

Candidate Number. General Certificate of Education Advanced Level Examination January 2012 entre Number andidate Number Surname Other Names andidate Signature Genera ertificate of Education dvanced Leve Examination January 212 Physics PHY4/1 Unit 4 Fieds and Further Mechanics Section Tuesday

More information

Introduction to LMTO method

Introduction to LMTO method 1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations

More information

International Journal of Mass Spectrometry

International Journal of Mass Spectrometry Internationa Journa of Mass Spectrometry 280 (2009) 179 183 Contents ists avaiabe at ScienceDirect Internationa Journa of Mass Spectrometry journa homepage: www.esevier.com/ocate/ijms Stark mixing by ion-rydberg

More information

arxiv:quant-ph/ v3 6 Jan 1995

arxiv:quant-ph/ v3 6 Jan 1995 arxiv:quant-ph/9501001v3 6 Jan 1995 Critique of proposed imit to space time measurement, based on Wigner s cocks and mirrors L. Diósi and B. Lukács KFKI Research Institute for Partice and Nucear Physics

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

NEW PROBLEMS. Bose Einstein condensation. Charles H. Holbrow, Editor

NEW PROBLEMS. Bose Einstein condensation. Charles H. Holbrow, Editor NEW PROBLEMS Chares H. Hobrow, Editor Cogate University, Hamiton, New York 3346 The New Probems department presents interesting, nove probems for use in undergraduate physics courses beyond the introductory

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

APPENDIX C FLEXING OF LENGTH BARS

APPENDIX C FLEXING OF LENGTH BARS Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or

More information

Quantum Electrodynamical Basis for Wave. Propagation through Photonic Crystal

Quantum Electrodynamical Basis for Wave. Propagation through Photonic Crystal Adv. Studies Theor. Phys., Vo. 6, 01, no. 3, 19-133 Quantum Eectrodynamica Basis for Wave Propagation through Photonic Crysta 1 N. Chandrasekar and Har Narayan Upadhyay Schoo of Eectrica and Eectronics

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t. Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes

More information

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing

More information

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006 expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische

More information

Collapse of a Bose gas: Kinetic approach

Collapse of a Bose gas: Kinetic approach PRAMANA c Indian Academy of Sciences Vo. 79, No. 2 journa of August 2012 physics pp. 319 325 Coapse of a Bose gas: Kinetic approach SHYAMAL BISWAS Department of Physics, University of Cacutta, 92 A.P.C.

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

Thermophoretic interaction of heat releasing particles

Thermophoretic interaction of heat releasing particles JOURNAL OF APPLIED PHYSICS VOLUME 9, NUMBER 7 1 APRIL 200 Thermophoretic interaction of heat reeasing partices Yu Doinsky a) and T Eperin b) Department of Mechanica Engineering, The Pearstone Center for

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary

More information

Chemical Kinetics Part 2. Chapter 16

Chemical Kinetics Part 2. Chapter 16 Chemica Kinetics Part 2 Chapter 16 Integrated Rate Laws The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates

More information

14-6 The Equation of Continuity

14-6 The Equation of Continuity 14-6 The Equation of Continuity 14-6 The Equation of Continuity Motion of rea fuids is compicated and poory understood (e.g., turbuence) We discuss motion of an idea fuid 1. Steady fow: Laminar fow, the

More information

Thermal Leptogenesis. Michael Plümacher. Max Planck Institute for Physics Munich

Thermal Leptogenesis. Michael Plümacher. Max Planck Institute for Physics Munich Max Panck Institute for Physics Munich Introduction Introduction Probem #1: the universe is made of matter. Baryon asymmetry (from nuceosynthesis and CMB): η B n b n b n γ 6 10 10 must have been generated

More information

Data Search Algorithms based on Quantum Walk

Data Search Algorithms based on Quantum Walk Data Search Agorithms based on Quantum Wak Masataka Fujisaki, Hiromi Miyajima, oritaka Shigei Abstract For searching any item in an unsorted database with items, a cassica computer takes O() steps but

More information

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier

More information

1) For a block of mass m to slide without friction up a rise of height h, the minimum initial speed of the block must be

1) For a block of mass m to slide without friction up a rise of height h, the minimum initial speed of the block must be v m 1) For a bock of mass m to side without friction up a rise of height h, the minimum initia speed of the bock must be a ) gh b ) gh d ) gh e ) gh c ) gh P h b 3 15 ft 3) A man pus a pound crate up a

More information

2M2. Fourier Series Prof Bill Lionheart

2M2. Fourier Series Prof Bill Lionheart M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier

More information

Supplemental Information

Supplemental Information Suppementa Information A Singe-eve Tunne Mode to Account for Eectrica Transport through Singe Moecue- Sef-Assembed Monoayer-based Junctions by A. R. Garrigues,. Yuan,. Wang, E. R. Muccioo, D. Thompson,

More information

1. Measurements and error calculus

1. Measurements and error calculus EV 1 Measurements and error cacuus 11 Introduction The goa of this aboratory course is to introduce the notions of carrying out an experiment, acquiring and writing up the data, and finay anayzing the

More information

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations) Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =

More information

This shoud be true for any, so the equations are equivaent to the Schrodingerequation (H E) : () Another important property of the (2) functiona is, t

This shoud be true for any, so the equations are equivaent to the Schrodingerequation (H E) : () Another important property of the (2) functiona is, t I. Variationa method Introduction Usefu for the determination of energies and wavefunctions of different atomic states. H the time-independent Hamitonian of the system E n eigenenergies Ψ n - eigenfunctions

More information

Chemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer. Lecture 12, October 1, 2003

Chemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer. Lecture 12, October 1, 2003 Cemistry 350 Pysica Cemistry II (Quantum Mecanics) 3 Credits Fa Semester 003 Cristoper J. Cramer Lecture 1, October 1, 003 Soved Homework We are asked to demonstrate te ortogonaity of te functions Φ(φ)

More information