Quantum Mechanics Beginners ******************************* Fundamental Concepts ================================== What're the most important tricks in QM calculations? --------------------------------------------------------- * **Remember what basis we are working in** * **Identity** First Three Postulates ------------------------- * Physical state is described by kets in a Hilbert space. We need to specify a complete basis {:math:`\ket{i}`} to do calculations. .. math:: \ket{\psi} = \sum_i \ket{i}\bra{i}\ket{\psi} = \sum_i C_i \ket{i} * Operators are given by Hermitian operators; A measurement of the variable :math:`\hat \Omega` will yield one of the eigenvalues :math:`\omega` with the probability .. math:: \left|\braket{\omega}{\psi}\right|^2 . And the state of the system will change to :math:`\ket{\omega}`. * The state vector obeys the Schrödinger equation, .. math:: i\hbar \frac{\d}{\d t}\ket{\psi(t)} = \hat H \ket{\psi(t)} , where :math:`\hat H` is the Hamiltonian operator. .. admonition:: Comments :class: note The logic here is that we first find the way to describe a system, then think about how to find out the information we need from the state vector and also find the evolution of the state vector. Then we need the operator and Schrodinger equation. Finally, we would like to relate the theory to experiments, and it comes the measurement postulate. Later we will need the relation between position and momentum, which becomes the fourth postulate. * How to solve the evolution of a system? We just define a magical operator, propagator .. math:: \hat U \ket{\psi(t_0)} = \ket{\psi(t)} . This operator just gives us the evolution of state vector! Wait, can we write down the explicit expression of it? Let's find out. The only thing we know about the evolution of a state vector is the third postulate up there. .. math:: i\hbar \frac{\d }{\d t}\ket{\psi(t)} & = \hat H \ket{\psi(t)} \\ i\hbar \frac{\d }{\d t}\hat U \ket{\psi(t_0)} & = \hat H \hat U \ket{\psi(t_0)} \\ i\hbar \frac{\d }{\d t}\hat U & = \hat H \hat U Looks familiar? This just gives us a exponential result, **if the Hamiltonian is time independent**. .. math:: \hat U = e^{- i \hat H (t-t_0)/\hbar} **We can prove that this operator is Unitary because** :math:`\hat H` **is Hermitian.** This is just the abstract representation, we work in some basis, and the most convenient basis is the eigenstates of Hamiltonian, { :math:`\ket{\epsilon_i}` }, .. math:: \hat U \ket{\phi} & = e^{- i \hat H (t-t_0)/\hbar} \ket{\psi} \\ \hat U \ket{\phi} & = \sum_i e^{- i \hat H (t-t_0)/\hbar} \ket{\epsilon_i}\bra{\epsilon_i} \ket{\psi} \\ \hat U \ket{\phi} & = \sum_i e^{- i \epsilon_i (t-t_0)/\hbar} \ket{\epsilon_i}\bra{\epsilon_i} \ket{\psi} And we are going to use .. math:: \hat U = \sum_i e^{- i \epsilon_i (t-t_0)/\hbar} \ket{\epsilon_i}\bra{\epsilon_i} from now on. (Well, only on discrete eigenvalues ones.) **(See that? Identity does the work again.)** Position and Momentum Space ----------------------------- .. admonition:: Summary :class: note * Position 1. Define {:math:`\ket{x}`} basis. 2. Define :math:`\hat x` operator. 3. Find wave function in this basis. 4. Find measurement. * Evolution 1. Need propagator :math:`\hat U`. 2. Propagator needs the solution of Hamiltonian eigensystem. 3. (Free particles) Hamiltonian needs the solution of momentum eigensystem. * Momentum 1. Before we define some arbitrary momentum space, we should check the relation between momentum and position. And it turns out to be related by a commutator.(Postulate IV) 2. Use the postulate to momentum operator. 3. Find eigenstates. 4. (Calculate the propagator.) Position Space """""""""""""""" 1. Define :math:`\ket{x}` basis. **Orthonormal condition is** .. math:: \braket{x}{x'}=\delta(x-x'). **Completeness condition is** .. math:: \int \braket{x'}{x'} \d x' = \mathbb{I} 2. Define position operator. The position operator is defined as .. math:: \hat x \ket{x} = x \ket{x} And in :math:`\{\ket{x}\}` basis, this operator becomes a function, which is .. math:: &\bra{x}\hat x \ket{x'} \\ =& \left(\bra{x}\hat x \right)\ket{x'} \\ =& x \braket{x}{x'} \\ =& x \delta(x-x') 3. Find state vector in :math:`\{\ket{x}\}` basis. .. math:: \psi(t,x) = \braket{x}{\psi(t)} **Normalized** .. math:: \int \left| \psi(t,x) \right|^2 \d x = 1. And we are interpreting :math:`\left| \psi(t, x)\right|^2` as probability density. 4. Calculate probability of a measurement. Taking :math:`\hat x` as an example. .. math:: &\bra{\psi} \hat x \ket{\psi} \\ &= \iint \braket{\psi}{x}\bra{x} \hat x \ket{x'} \braket{x'}{\psi} \d x \d x' \\ &= \iint \psi^ * (t,x) x\delta(x-x') \psi(t,x') \d x \d x' \\ &= \int \left| \psi(t,x) \right|^2 x \d x Momentum Space """"""""""""""" To find the momentum operator, we need to check the relation between momentum and position before we just randomly define one. Truth is, we have a fourth postulate states the relation between them. Postulate IV ^^^^^^^^^^^^^^ The commutator of :math:`\hat x`, :math:`\hat p` is .. math:: \left[ \hat x, \hat p \right] = i \hbar Two comments: * Why i ? Eigenvalue of Anti-Hermitian operator. * Why :math:`\hbar`? Because people define the dimensions of position and momentum differently before they know this commutator. We would like to assign them the same dimension if we already know this relation. Momentum Space ^^^^^^^^^^^^^^^^ 1. Find momentum operator in position basis {:math:`\ket{x}`}. .. math:: \bra{x} \left[ \hat x, \hat p\right] \ket{x'} = i\hbar \delta(x-x') And write out the commutator and use the relation of delta function :math:`x\delta'(x) = -\delta(x)`, we find out the momentum operator in {:math:`\ket{x}`} basis, .. math:: \bra{x}\hat p \ket{x'} = -i\hbar \frac{\d }{\d t} \delta(x-x') **Let's talk physics.** What does that operator mean? We need to see what the result is when momentum operator is applied to a state. And remember we would work in {:math:`\ket{x}`} basis. .. math:: &\bra{x} \hat p \ket{\psi} \\ & = \iint \braket{x}{x'} \bra{x'} \hat p \ket{x''}\braket{x''}{\psi} \d x' \d x'' \\ & = \int \bra{x}\hat p \ket{x''}\psi(t,x'') \d x'' \\ & = \int \left( -i\hbar \frac{\d}{\d x} \delta(x-x') \psi(t,x') \right) \d x' \\ & = \int \left( -i\hbar \frac{\d}{\d x'} \delta(x'-x) \psi(t,x') \right) \d x' **Integrate by parts, we will find the expression.** (I am having a problem finding the right answer.) .. math:: \bra{x} \hat p \ket{\psi} = - i\hbar \frac{\d }{\d x}\psi(x) . 2. Eigenfunction for momentum. .. math:: \hat p \ket{p} = p \ket{p} . Again, we are going to project it on the {:math:`\ket{x}`} basis, .. math:: \bra{x}\hat p\ket{p} = \bra{x} p \ket{p} , where :math:`\braket{x}{p}` is the eigenstates in {:math:`\ket{x}`} basis, we call it :math:`\phi_p(x)`. .. math:: \bra{x}\hat p\ket{p} & = p \phi_p(x) \\ \int \bra{x}\hat p \ket{x'}\braket{x'}{p}\d x' & = p \phi_p(x) \\ -i\hbar \frac{\d }{\d x} \phi_p(x) & = p \phi_p(x) The solution is .. math:: \phi_p(x) = \mathrm{C} e^{i p x/\hbar} This constant C is found by the normalization condition, .. math:: \braket{p}{p'}=\int \phi_p^*(x)\phi_{p'}(x)\d x = \delta(p-p') **The final results should be** .. math:: \phi_p(x)=\frac{1}{\sqrt{2\pi \hbar}} \exp{(i p x/\hbar)} 3. Find the dynamics of free particles in quantum mechanics. **Find the propagator and everything solves.** The hamiltonian for a free particle is .. math:: \hat H = \frac{\hat p^2}{2m} . We argue here that the eigenvectors of momentum are also the eigenvectors of this hamiltonian. And we can easily guess the eigenvalues are :math:`p^2/2m`. So the propagator is .. math:: \hat U = \int e^{-i p^2 t/2m\hbar} \ket{p}\bra{p} \d p But that is too abstract to use, we can find the expression in {:math:`\ket{x}`} basis. .. math:: \bra{x}\hat U\ket{x} & = \int e^{-i p^2 t/2 m \hbar} \braket{x}{p}\braket{p}{x} \d p \\ & = \int e^{-i p^2 t/2 m \hbar} \left| \phi_p \right|^2 \d p Quantum in 1D ============== General ---------- Always start with the propagator for time independent Hamiltonian. .. math:: \ket{\psi(t)} = \hat U \ket{\psi(0)} For cases that Hamiltonian with discrete eigenvalues, .. math:: \ket{\psi(t)} = \sum _ n e^{-i \epsilon _ n t/ \hbar } \ket{n}\braket{n}{\psi(0)} If the initial state is just one of the eigenstates of Hamiltonian, say the mth one (normalized), .. math:: \ket{\psi(t)} = e^{- i \epsilon _ m t/\hbar} \ket{ m } Well, that phase factor doesn't have any effect for the topic we discuss. So our time evolution will stay on the same state forever. The same thing happens for continuous cases. So our task is simplified to solve the eigensystem of Hamiltonian, which is .. math:: \hat H \ket{\epsilon} = \epsilon \ket{\epsilon} Infinite Barriers """""""""""""""""" Math ^^^^^ Setup ~~~~~~ * Potential in a box .. math:: V(x) = & 0, 0< x & = \left( \frac{\d}{\d t}\bra{\psi(t)}\right) \hat \Omega \ket{\psi(t)} + \bra{\psi(t)} \dot{\hat \Omega} \ket{\psi(t)} + \bra{\psi(t)} \hat \Omega \left( \frac{\d}{\d t}\ket{\psi(t)}\right) \\ & = \frac{1}{i\hbar} \left ( - \bra{\psi(t)} \hat H \hat\Omega \ket{\psi(t)} +\bra{\psi(t)} \hat\Omega \hat H \ket{\psi(t)} \right) + \bra{\psi(t)} \dot{\hat \Omega} \ket{\psi(t)} \\ & = \frac{1}{i\hbar} \bra{\psi(t)}\left[\hat\Omega,\hat H\right] \ket{\psi(t)}+\bra{\psi(t)} \dot{\hat \Omega} \ket{\psi(t)} This is called Ehrenfest's Theorem. Simple Example of Ehrenfest's Theorem """"""""""""""""""""""""""""""""""""""" Suppose we have a system with Hamiltonian .. math:: \hat H = \frac{\hat p^2}{2m} + V(\hat x) We need to figure some commutators first. .. math:: 2m \left[ \hat x, \hat H \right] =\left[\hat x, \hat p^2\right] = \hat x \hat p\hat p - \hat p \hat p \hat x = \hat x \hat p \hat p -\hat p \hat x \hat p + \hat p \hat x \hat p - \hat p \hat p \hat x = \left[\hat x , \hat p\right]\hat p + \hat p \left[ \hat x,\hat p\right] = 2 i \hbar \hat p .. math:: \left[\hat p, \hat H\right] = \left[\hat p, V(\hat x) \right] = \left[\hat p, \sum_0^\infty \frac{V^{(n)}}{n!}\hat x^n\right] =\cdots =-i\hbar V'(\hat x) 1. Position average .. math:: \frac{\d }{\d t} \left< \hat x \right> & = \frac{1}{i\hbar} \bra{\psi(t)} \left[ \hat x, \hat H \right]\ket{\psi(t)} \\ & = \frac{\left< \hat p \right>}{m} We are familiar with this in classical mechanics. 2. Momentum average .. math:: \frac{\d}{\d t} \left<\hat p\right> & = \frac{1}{i\hbar} \bra{\psi(t)} \left[\hat p, \hat H\right] \ket{\psi(t)} \\ & = \frac{1}{i\hbar} \bra{\psi(t)} (-i\hbar V'(\hat x)) \ket{\psi(t)} \\ & = -\left< V'(\hat x) \right> In classical mechanics, the derivative of potential is force. And the result is just like Newton's 2n Law except the right hand side is not exactly like a force which should be :math:`-\frac{\d}{\d x} \left< V(\hat x) \right>`. What does :math:`-\left< V'(\hat x)\right>` mean """"""""""""""""""""""""""""""""""""""""""""""""""" Suppose the potential area is fairly small and distributed around some coordinate :math:`x_0=\left< \hat x \right>`, we can do Taylor expansion around :math:`x_0`. .. math:: < V(\hat x)> & = V(x_0) + V'(x_0) < (x - x_0)> + V''(x_0)<(x-x_0)^2> /2 + \cdots \\ & = V(x_0) + 0 + V''(x_0) (\Delta x)^2 + \cdots If the uncertainty is small enough, every term except the first one becomes small. So to the lowest order, average of potential is approximately the potential at :math:`x_0`. Similarly, the average of first derivative of potential :math:`` is approximately :math:`V'(x_0)`. These gives us a hint for the previous result we got for the time evolution of average momentum. The result reduces to classical mechanics one as long as we keep the lowest order of Taylor expansion. Those higher order terms show the quantum effect. Picture ----------- We can see deeper into Ehrenfest's Theorem through Heisenberg Picture of quantum mechanics. Schrödinger & Heisenberg Pictures """""""""""""""""""""""""""""""""""" Pictures are the ways we look at the evolution of systems. Schrödinger Picture ^^^^^^^^^^^^^^^^^^^^ In Schrödinger picture the states are evolving with time. .. math:: i\hbar \frac{\d}{\d t} \ket{\psi} _ S = \hat H \ket{\psi} _ S And for time independent Hamiltonian, .. math:: \ket{\psi}_S = U^\dagger \ket{\psi _ 0} _ S Heisenberg Picture ^^^^^^^^^^^^^^^^^^^^ In Heisenberg Picture, the states do not change with time. .. math:: \ket{\psi} _ H = \ket{\psi_0} _ H , and of course the initial is the same with Schrödinger Picture, .. math:: \ket{\psi_0} _ H = \ket{\psi _ 0} _ S . How do we relate to Heisenberg Picture to Schrödinger Picture? Through investigation of observables. We should have the same observation results in both Pictures. .. math:: {} _ H \bra{\psi} \hat \Omega _ H \ket{\psi} _ H & = {} _ S \bra{\psi} \hat \Omega _ S \ket{\psi} _ S \\ {} _ H \bra{\psi} \hat \Omega _ H \ket{\psi} _ H & = {} _ S \bra{\psi _ 0} \hat U^\dagger \hat \Omega _ S \hat U \ket{\psi _ 0} _ S \\ \hat \Omega _ H & = \hat U^\dagger \hat \Omega _ S \hat U So the operators change with time in Heisenberg Picture. Ehrenfest's Theorem in Heisenberg Picture """""""""""""""""""""""""""""""""""""""""""" .. math:: \frac{\d }{\d t} \hat \Omega _ H = \frac{1}{i\hbar } \left[ \hat \Omega _ H, \hat H \right] + \hat U ^ \dagger \frac{\partial }{\partial t} \Omega _ H \hat U This can be easily proved by throwing every definition need in to it. We also need the following equations. .. math:: \frac{\d }{\d t} \hat U = \frac{\d }{\d t} e^{-i\hat H t /\hbar} = \frac{\hat H}{i\hbar} \hat U And REMEMBER that propagator commute with time independent Hamiltonian, so .. math:: \hat H = \hat U^\dagger \hat U \hat H = \hat U^ \dagger \hat U \hat U \equiv \hat H _ H So this Ehrenfest's Theorem can also be written as .. math:: \frac{\d }{\d t} \hat \Omega _ H = \frac{1}{i\hbar } \left[ \hat \Omega _ H, \hat H _ H \right] + \hat U ^ \dagger \frac{\partial }{\partial t} \Omega _ H \hat U We can **define** .. math:: \frac{\partial}{\partial t}\hat \Omega _ H \equiv \hat U^\dagger \frac{\partial }{\partial t}\hat \Omega _ S \hat U , which is the time derivative of operator in Heisenberg Picture. **Reminder: The time derivative of an observable (average) depends not only the time derivative of itself, but also the commutator of the observable and Hamiltonian.** Example of Ehrenfest's Theorem in Heisenberg Picture ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ We will show why it is better to work in Heisenberg Picture to show the meanings of Ehrenfest's Theorem. Suppose we have a Hamiltonian in Heisenberg Picture, .. math:: \hat H_H = \frac{\hat p _ H^2 }{2m} + V(\hat x _ H) . Time derivative of position operator .. math:: \frac{\d}{\d t} \hat x _ H = \frac{1}{i\hbar} \left[\hat x _ H, \hat H _ H \right ] = \frac{\hat p _ H}{m} Time derivative of momentum operator .. math:: \frac{\d}{\d t} \hat p_H = \frac{1}{i\hbar } \left[ \hat p _ H, \hat H \right] = - V'(\hat x_H) So the operator in Heisenberg Picture just have a sense of the physical quantities in classical mechanics. That's why we like it. Comparison of Picutres """""""""""""""""""""""""""""" .. figure:: assets/images/picutres-comparison.png :align: center Comparison of different pictures. Notice that in Dirac picture, :math:`\hat W_I = \hat U_0 ^{-1}\hat W \hat U_0`, :math:`\ket{\psi(t)}_I = \hat U_I \ket{\psi(0)}_I`. A markdown file that is used to make this table can be downloaded `here `_ . Conservation --------------- We say a observable is conserved if the corresponding operator commutes with Hamiltonian, .. math:: \left[ \hat \Omega, \hat H \right]=0 1. Energy Hamiltonian always commutes with itself. .. math:: \frac{\d}{\d t} \left<\epsilon \right> = \bra{\psi} \left( \frac{\partial }{\partial t} \hat H \right) \ket{\psi} If Hamiltonian is time independent, then energy is conserved. (If Hamiltonian is tide dependent, energy is not conserved. This is kind of obvious in classical mechanics.) What is the nature of time dependence """"""""""""""""""""""""""""""""""""""" We can see this by looking at a simple example. Assume we have a system with energy eigenstates :math:`\ket{\epsilon _ n}`, and initially, .. math:: \ket{\psi _ 0} = \sum_n C _ n \ket{\epsilon _ n} . So .. math:: \ket{\psi(t)} = \sum _ n C _ n e^{-i\epsilon _ n t/\hbar} \ket{\epsilon _ n} . We can calculate the expectation value of some operator :math:`\hat \Omega`, .. math:: \left< \omega (t) \right> = \sum _ {n,m} \left( C _ n^ * e^{i\epsilon _ n t/\hbar } \bra{\epsilon _ n} \right) \hat \Omega \left( C _ m e^{-i \epsilon _ m t/\hbar} \ket{\epsilon _ m} \right) = \sum _ {n,m} C _ n ^* C _ m e^{-i(\epsilon _ m - \epsilon _ n) t/\hbar} \bra{\epsilon _ n} \hat \Omega \ket{\epsilon _ m} If :math:`\ket{\epsilon _ n}` are also the eigenvectors of :math:`\hat \Omega`, then .. math:: \bra{\epsilon _ n} \hat \Omega \ket{\epsilon _ m} = \omega _ m \delta _ {n,m} And the expectation value .. math:: \left< \omega (t) \right> = \sum _ {n} C _ n^* C _ n \omega _ n **The important thing is that the time dependence of this expectation value actually arise from this term** .. math:: e^{-i(\epsilon _ m - \epsilon _ n)t/\hbar} . As it is so important, we call .. math:: (\epsilon _ m - \epsilon _ n)/\hbar **Bohr frequency**. Harmonic Oscillators ===================== Why Harmonic Oscillators -------------------------- Many systems can reduce to it. Use Taylor expansion for the potential and redefine parameters we will find harmonic oscillators in the potential. Hamiltonian for 1D is .. math:: \hat H = \frac{\hat p^2}{2m} + \frac{1}{2} k \hat x^2 Standard Solution -------------------- We can use polynomial expansion for part of the solution. Dimension Schrodinger Equation """""""""""""""""""""""""""""""" First step is always finding out the characteristic length scale and characteristic energy scale. Assume we have an characteristic length :math:`\eta` and characteristic energy scale :math:`\epsilon_0`. Through uncertainty principle we know only for dimensional analysis .. math:: \left[\hat p\right]=\frac{\hbar}{\eta} Kinetic energy and potential energy have the same dimension .. math:: \frac{\hbar^2}{\eta^2 m}=k \eta^2 , so we have .. math:: \eta = \sqrt{\frac{\hbar}{m\omega}} with :math:`\omega^2=k/m`. A dimensional analysis shows that :math:`\epsilon_0=\hbar\omega`. Now we can define dimensionless variables, .. math:: z=x/\eta, e=\epsilon/\epsilon_0 The time independent Schrodinger equation in position basis is .. math:: -\hbar^2 \frac{\mathrm d^2}{\mathrm dx^2}\psi'' /m + k x^2 = 2\epsilon \psi . Using those characteristic scales, we can rewrite this equation into a dimensionless one, which is .. math:: \psi''+(2e-z^2)\psi = 0 in which :math:`\psi'=\frac{\mathrm d}{\mathrm dz}\psi`. Take Limits """""""""""""""""""" We need to look at the behavior of the solutions before we can guess a proper general solution. :math:`z\rightarrow \infty`, we have :math:`\psi''-z^2\psi=0`. Solution to this equation is :math:`\psi(z)~ e^{-z^2/2}`. The solution of the the equation should be in the form .. math:: \psi(z) = u(z) e^{-z^2/2} . Insert this to time independent Schrodinger equation, we can get the equation of :math:`u(z)`. .. math:: u'' - 2 z u' +(2e-1) u = 0 Polynomial Method """"""""""""""""""""" The simplest form of :math:`u(z)` is polynomial, .. math:: u(z) = \sum _ {n=0}^{\infty} u _ n z^n . Put this back to equation of u, we can get the recursion relation, .. math:: (n + 2)(n+1) u _ {n+2} = \left[ 2n - (2e - 1) \right] u _ n . If :math:`u_0` and :math:`u_1` are given, we can get the whole polynomial. Notice that we have definite parity here. So :math:`u _ 1` branch vanish because they are even. :math:`u_0` is set by the normalization condition. Terminate The Series """""""""""""""""""""" The series blow up if it doesn't terminate. So we need to terminate the series using the following relation, .. math:: 2e - 1 = 2n . Then we have the energy levels, which is :math:`e=n+1/2`. Complete Series """""""""""""""" By picking proper normalization factor, we can write down the energy levels and corresponding wave functions. In fact, this polynomial can be found in mathematical phyisics books. .. math:: H _ {n+1} = 2 z H _n -n H _ {n-1} Tricky Solution ----------------- Find out the characteristic length and energy .. math:: \eta = \sqrt{\frac{\hbar }{m\omega }} \\ \epsilon = \hbar \omega \\ \omega = \sqrt{\frac{k}{m}} One way to get the intrinsic length without writing down the dimensions of each quantity is to use the following relation .. math:: \left[ E \right] = \left[ m \omega^2 \hat x^2 \right] \\ \hbar \omega = m \omega^2 \eta^2 \\\ \eta = \sqrt{ \frac{\hbar}{m\omega} } Or if we are given the Hamiltonian in terms of :math:`k`, .. math:: \left[ \frac{\hat p^2}{2m} \right] = \left[ k \hat x^2 \right] \\ \frac{\hbar^2 / \eta^2 }{m} = k\eta^2 \\ \eta = \sqrt{\hbar}{ \sqrt{m k} } = \sqrt{ \hbar }{ m \omega } Rewrite the Hamiltonian .. math:: \hat H & = \frac{1}{2m} \left[ \left(\frac{\hat p}{\hbar/\eta}\right)^2 \left(\frac{\hbar}{\eta}\right)^2 + \frac{1}{2} m \omega^2 \left( \frac{\hat x}{\eta} \right)^2 \right] \\ & = \frac 1 2 \hbar \omega \left[ \left(\frac{\hat p}{\hbar/\eta}\right)^2 + \left(\frac{\hat x}{\eta}\right)^2 \right] \\ & = \frac 1 2 \hbar \omega \left( \frac{\hat x}{\eta} - i \frac{\hat p}{\hbar/\eta} \right) \left( \frac{\hat x}{\eta} + i\frac{\hat p}{\hbar/\eta} \right) - \frac{i}{\hbar} \left[\hat x, \hat p\right] \\ & = \frac 1 2 \hbar \omega (\sqrt 2 \hat a^\dagger \sqrt 2 \hat a + 1) \\ & = \hbar \omega \left( \hat a^\dagger \hat a + \frac 1 2\right) Now we can define :math:`\hat a^\dagger \hat a = \hat N`, which is just like an operator for (energy) quanta numbers. An impoertan relation is .. math:: \left[\hat a, \hat a^\dagger\right] = 1 \\ \left[\hat a, \hat N\right] = \hat a The eigen equation for this weird energy quanta number operator is .. math:: \hat N \ket{n} = n \ket{n} To find out the eigen state of :math:`\hat a` and :math:`\hat a^\dagger`, we try this, .. math:: \hat N (\hat a \ket{n}) = (n-1) (\hat a \ket{n}) \\ \hat N (\hat a^\dagger \ket{n}) = (n+1) (\hat a^\dagger \ket{n}) This means :math:`\hat a \ket{n}` and :math:`\hat a^\dagger \ket{n}` are also eigen states of :math:`\hat N`. The next step is very crucial. Since :math:`\hat a \ket{n}` and :math:`\hat a^\dagger \ket{n}` are eigen states of :math:`\hat N`, we know that .. math:: \hat a \ket{n} = C1 \ket{n} \\ \hat a^\dagger \ket{n} = C2 \ket{n} Then our next step is to find out what are :math:`C1` and :math:`C2` exactly. They way of finding them is to use invariant quantities, such as the inner product. Here we use average of :math:`\hat N` operator. .. math:: \hat a \ket{n} = \sqrt n \ket{n-1} \\ \hat a^\dagger \ket{n} = \sqrt{n+1} \ket{n+1} Final step is to constrain on :math:`n`, which should be integrals. This is true because we need a cut off for the eigen equation of :math:`\hat N`, whose avarage is n and it should be non negative. .. math:: \bra{n}\hat N \ket{n} \ge 0 leads to :math:`n\ge 0`. To get this proper cut off, :math:`n` should be integer because if it's not, according to .. math:: \hat a \ket{n} = \sqrt n \ket{n-1} n can go to negative numbers. If n is positive integer, .. math:: \hat a \ket{1} = \ket{0} \\ \hat a \ket{0} = 0 \ket{0} show an cut off at 0. We can even find out the wave functions of these :math:`\ket{n}` by finding the ground state first and apply :math:`\hat a^\dagger` to the ground state. Ground state in :math:`{\ket{x}}` basis can be found by solving the differential equation, .. math:: \bra{x} \hat a \ket{0} = 0 .. Very important: * The Hermitian conjugate of :math:`\hat a \ket{n}` is :math:`\bra{n} \hat a^\dagger`. * Hermitian conjugate of :math:`\hat a \hat a^\dagger` is :math:`\hat a \hat a^dagger`. This can be a trap. Hermitian conjugate is the complex conjugate AND TRANSPOSE! Semiclassical ------------------- Classical """""""""""""""""""""" In phase space, the trajectory of phase space points ( {:math:`x/\eta` and :math:`p/(\hbar/\eta)`} ) is on a circle of radius :math:`x_{max}/\eta`. Quantum semiclassical """"""""""""""""""""""" Key points: 1. What is the trajectory of :math:`\left<\hat x/\eta\right>` and :math:`\left<\hat p/(\hbar/\eta)\right>` 2. Can we make the trajectory just like the classical case by choosing some special conditions? 3. What do these special cases mean? * Expectation value of creation and annihilation operators Apply Ehrenfest theorem to annihilation operator, .. math:: i\hbar \frac{\mathrm d}{\mathrm d t} \avg{\hat a(t)} = \bra{\psi} \left[ \hat a(t), \hat H \right] \ket{\psi} = \hbar \omega \avg{\hat a(t)} Excellent. Now we can solve out :math:`\avg{\hat a(t)}`, which is .. math:: \avg{\hat a(t)} = \alpha_0 \exp(-i\omega t) Take the hermitian conjugate, .. math:: \avg{\hat a^\dagger (t)} = \alpha_0^* \exp(i\omega t) * Expectation value of position and momentum With these two operators, we can find out the average of :math:`\hat x` and :math:`\hat p` because .. math:: \hat x = \eta \frac{1}{\sqrt 2} \left( \hat a^\dagger + \hat a\right)\\ \hat p = \frac{\hbar}{\eta} i \frac{1}{\sqrt 2} \left(\hat a^\dagger - \hat a \right) , we have .. math:: \avg{\hat x(t)} = \eta \frac{1}{\sqrt 2} \left( \avg{\hat a^\dagger (t)} + \avg{\hat a(t)} \right) \\ \avg{\hat p(t)} = \frac{\hbar}{\eta} i \frac{1}{\sqrt {2} } \left( \avg{\hat a^\dagger (t) - \avg{\hat a(t)}} \right) We can have a look at these two averages, .. math:: \frac{\avg{\hat x(t)} }{\eta} = \frac{1}{\sqrt{2} } \left[ (\alpha_0 + \alpha_0^*)\cos(\omega t) + i (\alpha_0^* - \alpha_0 ) \sin(\omega t) \right] \\ \frac{\avg{\hat p(t)}}{\hbar/\eta} = \frac{1}{\sqrt{2}} \left[ (\alpha_0 + \alpha_0^*) \sin(\omega t) + i( \alpha_0 - \alpha_0^*)\cos(\omega t) \right] It is obvious that the average reduces to classical case if :math:`\alpha_0 = \alpha_0^*`. **But this is too strong for a semiclassical limit.** * Coherent state **Coherent state is the eigenstate of creation operator. Its wave package has the smallest spread allowed by quantum mechanics.** **The most special part about coherent state is that the system stays on coherent state if it start with coherent state.** .. math:: \hat a \ket{\alpha(t)} = \alpha(t) \ket{\alpha(t)} Take the hermitian conjugate, .. math:: \bra{\alpha(t)} \hat a^\dagger = \bra{\alpha(t)}\alpha(t)^* At :math:`t=0`, we have .. math:: \bra{\psi(0)} N \ket{\psi(0)} = \vert \alpha_0 \vert ^2 That is to say, energy should be .. math:: \bra{\psi(0)} \hat H \ket{\psi(0)} = \hbar \omega \left( \vert \alpha_0 \vert^2 + \frac{1}{2} \right) Initially, we also have .. math:: \bra{\psi(0)} (\hat a - \alpha_0)^\dagger (\hat a-\alpha_0) \ket{\psi(0)} = 0 This means .. math:: \hat a \ket{\psi(0)} = \alpha_0 \ket{\psi(0)} * Coherent state expanded using energy eigenstates (This result) (To Be Finished...)