Lecture 5 Spin Hamiltonians Angular momentum. Angular momentum operators in matrix representation. Angular momentum as a spin of higher size. Spin Hamiltonians. Addition of angular momenta. Clebsh‐Gordun coefficients. 1) Angular momentum Hydrogen atom and its es (eigenstates) R radial part of wave function. θ, ϕ – spherical harmonic, angular part of w.f. – principal quantum number, defines energy level. ℓ ‐ angular momentum quantum number, defines magnitude of the angular momentum. projection quantum number, defines projection of angular momentum on axis z (taken as a convention). – orbital degeneracy number. Classical angular momentum ⃗ ⃗⃗. Quantum angular momentum is a vector operator with three components All these operators can be represented in spherical coordinates θ, ϕ. See textbook. Spherical harmonics are the eigenstates of operators and . From definition of angular momentum (eq.1) it can be shown that they satisfy the following relation. For a given ℓ , takes the numbers ℓ,ℓ1,..0,..,ℓ resulting in degeneracy 2ℓ1. 1 We can also introduce rising and lowering operator, which are defined and act as following: (see textbook for details) From the above relation one can see that the set of angular wave functions, ,…, ,…., , is a closed set, meaning that result of action of any angular momentum operator on can be expressed as a linear superposition of spherical harmonics from this set. So it is customary to speak about such objects as higher spins. We can create the corresponding matrix description of angular momentum (spin) rd operators in the basis of . See problems in Griffiths (3 ed. 4.34, 4.61, 4.62) . Let’s do this for a particle with the spin 1, that is with ℓ1. First I switch, for consistency, from spherical harmonics to z‐components of spin 1. and then use them as basis vectors The matrix for operator can be easily constructed. You can check that it acts on the basis vectors in accordance with eq.5. In the next step we find how rising and lowering operator acts on the basis vectors using eq. (8,9) Next we express , via , and construct the matrices. 2 Continuing this procedure, we get At this point we can equate meaning by a higher spin Generalization: 1) From atomic orbitals (spherical harmonics) we can construct matrices for operators with rank 0,3,5,7.. corresponding to operators with 0,1,2, …. Purely .. mathematically, using expression for raising and lowering operators we can get matrices for 1/2, 3/2,5/2,.. Most importantly we get op erators with the lowest non‐zero rank 1/2, that is starting from , get matrices for ,. This is remarkable, because we do not know how to write the wave function corresponding to say / 2) The same algebra of matrix operators works when we deal with elementary particles in high energy physics and in atomic physics. 3 2. Spin Hamiltonians s 2.1. Spin in magnetic field In general, μ⃗ γ⃗ , where μ⃗ magnetic moment, γ gyromagnetic ratio, ⃗ spin vector operator. || For spin ½ of electron μ⃗ g ⃗, where 2 so‐called g‐factor. Another popular notation convention used both for spin ½ and higher spin J. μμ , where is the total angular momentum or spin measured in integer of half‐integer ||ℏ numbers : ‐1,0,+1 or ‐1/2 , +1/2. And μ is Bohr magneton. Notice that ℏ is moved from spin to μ. The minus sign in eq.2 and gyromagnetic ratio appears because magnetic moment is antiparallel to spin angular momentum for an electron. (electron has negative charge) The energy of magnetic moment placed in magnetic field in classical EM is μ⃗⃗ In QM, it is replaced with operator Hamiltonian (and this is our first spin Hamiltonian) Problem: ⃗ 0,0, , at 0, spin is in the state , that is it point into positive x‐ direction. Q.1 Describe time evolution of the spin. Here Ωγ is Larmor frequency. 4 The state of the spin at 0, can be represented as Then its time evolution is given by the time evolution of the eigenstates, Q.2 Find ⟨⟩ average value of as a function of time. ℏ You can further show that Ω. So the physical picture is similar to the classical behavior. Spin rotates (precesses) around z axis in x‐y plane. Energy level splitting picture #1: Arrows indicate direction of magnetic moment. This picture is used in literature all the time. Energy level splitting picture #2. Arrows indicate direction of the spin. It is different for electron (negative particle) and proton (positive particle) 5 2.2 Spin Hamiltonians with the product term. Spin‐orbit, Heisenberg ….. The are many examples of physical systems where Hamiltonians depend on interaction between spins. Let’s give some examples. 1) Spin‐orbit coupling. Fine structure of hydrogen and other atoms. Griffiths p.299. I do not give all details just basics facts that allow us to write the Hamiltonian. Static charge. Electrical field acts on the charge, magnetic field on magnetic moment. However, moving electron experience Lorentz force ⃗⃗. There is a mirror effect. A particle moving in electrical field “sees” magnetic field ⃗ ⃗⃗. Handwave explanation: In our usual frame, electron rotates around proton. However, in the reference frame of electron, proton rotates around it and produces magnetic field in the electron position. The orbital motion of electron is related to its angular momentum. The Hamiltonian has the following form (actually it is a small correction to Coulomb potential). ⃗ and ⃗ characterize the same electron. 2) Hyperfine interaction (splitting) – interaction between spin of electron and spin of proton (nuclei) 6 Here we have two magnetic moments interacting with each other via dipole‐dipole magnetic interaction. Additional “contact” interaction occurs for s electrons, because they have finite probability to be found at the origin, where the proton (nuclei) is located. When ℓ0 one can find: 3) Heisenberg model In some approximation, the Hamiltonian of two interacting electrons in, say, diatomic molecule can be represented as: (Comm: we will derive this result) ⃗⃗ This Hamiltonian is known as the Heisenberg model. is the exchange constant. It can be both positive and negative. The Heisenberg model can be extended to the case of many spins in 1,2,3 dimension. Spin can be ½, 1, …. The model plays important role in the theory of magnetism. To summarize, we have a rather generic model Hamiltonian, which has the product of two vector operators. Let’s use the matrix method to solve it for the simplest case of two interacting ½ spins. Importantly, here we start dealing with many‐body QM, so the approach has many generic features. First, we simplify Hamiltonian using Pauli matrices. Plank constant are absorbed in . The basis vectors are formed as direct product of states of the first spin and the second spin. Notation with arrows are often used in literature. I prefer more explicit notation which specify the state and particle. For completeness I also list the number of the state in the Hamiltonian basis and corresponding column. 7 To compute matrix elements of H such as we need to find an action of individual Pauli matrices on individual spin states Let’s demonstrate how we find matrix element for Heisenberg Hamiltonian. The rule is each operator acts on its own spin sate 1 on 1, 2 on 2. Some intermediate results needed for computation of matrix elements. ∗ Matrix is Hermitian so . Here is what we get 8 We see that our matrix of rank 4 splits into three square matrices A, B, C of rank 1,2, and 1, respectively. The procedure of finding eigenstates and eigenvalues for these matrices can be done independently. We see that the eigenstates of the Hamiltonian can be split into two groups. The group with form multiplet corresponding to the total spin equal 1 (in ℏ units) and projections m=+1 for , m=‐1 for , and m=0 for /√2. (You will be asked to prove the last statement in the homework problem. ) This multiplet has the special name triplet. The state with 3 is called singlet The picture of the energy splitting when (positive) exchange interaction is present can be summarized as following. 9 3. Addition of angular momentums. The result that we have obtained for the Heisenberg model can be obtain in a simpler, as some think, method of the addition of angular momentums. I still believe the matrix method is much superior – it uncovers all details, allows to compute all parameters and does not need analogy with classical physics case often used in the textbooks. Still for the sake of completeness I would like to consider the addition method as well. We consider the spin‐orbit coupling term, where ⃗ is operator of the electron angular momentum and ⃗ its spin operator Next we introduce total angular momentum operator and make some manipulations. Our starting choice of the basis vectors is all combinations of angular momentum and spin component such as for example: ⃗ ⃗ The is not the eigenstate of the operator . However, we can go around this problem if we notice the following relations So, we found that operators , , , commute with ⃗ ⃗ and therefore they have the same eigenstates and eigenvalues. The Hamiltonian in the basis | ,, , ⟩ takes diagonal form. Eigenvalues can be computed using relations Which can be generalized for total angular momentum 10 So the eigenvalues of the ⃗⃗ take the form Let’s see how it works for the case we just have solved using the matrix method, 1/2, 1/2. Here we have two options for total angular momentum 1, 0. We some dropped common factor we have got the same results from the matrix method.