|
|
Line 1: |
Line 1: |
| ''' Welcome to the Quantum Mechanics A PHY5645 Fall2008''' | | [[Image:Quantum.png]] |
| [[Image:SchrodEq.png|thumb|Schrodinger equation. The most fundamental equation of quantum mechanics which describes the rule according to which a state <math>|\Psi\rangle</math> evolves in time. | | |
| ]] | | ''' Welcome to the Quantum Mechanics A PHY5645 Fall2008/2009''' |
| | [[Image:SchrodEq.png|thumb|550px|<b>[[Schrödinger Equation]]</b><br/>The most fundamental equation of quantum mechanics; given a Hamiltonian <math>\mathcal{H}</math>, it describes how a state <math>|\Psi\rangle</math> evolves in time.]] |
| | |
| This is the first semester of a two-semester graduate level sequence, the second being [[phy5646|PHY5646 Quantum B]]. Its goal is to explain the concepts and mathematical methods of Quantum Mechanics, and to prepare a student to solve quantum mechanics problems arising in different physical applications. The emphasis of the courses is equally on conceptual grasp of the subject as well as on problem solving. This sequence of courses builds the foundation for more advanced courses and graduate research in experimental or theoretical physics. | | This is the first semester of a two-semester graduate level sequence, the second being [[phy5646|PHY5646 Quantum B]]. Its goal is to explain the concepts and mathematical methods of Quantum Mechanics, and to prepare a student to solve quantum mechanics problems arising in different physical applications. The emphasis of the courses is equally on conceptual grasp of the subject as well as on problem solving. This sequence of courses builds the foundation for more advanced courses and graduate research in experimental or theoretical physics. |
| | | |
| The key component of the course is the collaborative student contribution to the course Wiki-textbook. Each team of students is responsible for BOTH writing the assigned chapter AND editing chapters of others. | | The key component of the course is the collaborative student contribution to the course Wiki-textbook. Each team of students is responsible for BOTH writing the assigned chapter AND editing chapters of others. |
|
| |
|
| ----
| | '''Team assignments:''' [[Phy5645_Fall09_teams|Fall 2009 student teams]] |
| '''Outline of the course:''' | |
| | |
| | |
| == Physical Basis of Quantum Mechanics ==
| |
| | |
| === Basic concepts and theory of motion in QM ===
| |
| In Quantum Mechanics, all of the information of the system of interest is contained in a wavefunction , <math>\Psi\,\!</math>. Physical properties of the system such as position, linear and angular momentum, energy, etc. can be represented via linear operators, called observables. These observables are a complete set of commuting Hermitian operators, which means that the common eigenstates (in the case of quantum mechanics, the wavefunctions) of these Hermitian operators form an orthonormal basis. Through these mathematical observables, a set of corresponding physical values can be calculated.
| |
| | |
| In order to clarify the paragraph above, consider an analogous example: Suppose that the system is a book, and we characterize this book by taking measurements of the dimensions of this book and its mass (The volume and mass are enough to characterize this system). A ruler is used to measure the dimensions of the book, and this ruler is the observable operator. The length, width, and height (values) from the measurements are the physical values corresponding to that operator (ruler). For measuring the weight of the book, a balance is used as the operator. The measured mass of the book is the physical value for the corresponding observable. The two observable operators (the ruler and the mass scale) have to commute with each other, otherwise the system can not be characterized at the same time, and the two observables can not be measured with infinite precision.
| |
| | |
| In quantum mechanics, there are some measurements that cannot be done at the same time. For example, suppose we want to measure the position of an electron. What we would do is send a signal (a gamma ray, for example), which would strike the electron and return to our detectors. We have, then, the position of the electron. But as the photon struck the electron, the electron gained additional momentum, so then our simultaneous momentum measurement could not be precise. Therefore both momentum and position cannot be measured at the same time. These measurement are often called "incompatible observables." This is explained in the Heisenberg uncertainty principle and implies, mathematically, that the two operators do not commute.
| |
| | |
| This concept contrasts with classical mechanics, where the two observables that do not commute with each other can still be measured with infinite precision. This is because of the difference in dimension of the object: macroscopic (classical mechanics) and microscopic scale (quantum mechanics). However, the prediction of quantum mechanics must be equivalent to that of the classical mechanics when the energy is very large (classical region). This is known as the Correspondence Principle, formally expressed by Bohr in 1923.
| |
| | |
| We can explain this principle by the following:
| |
| In quantum mechanics, particles cannot have arbitrary values of energy, only certain discrete values of energy. There are quantum numbers corresponding to specific values of energy and states of the particle. As the energy gets larger, the spacing between these discrete values becomes relatively small and we can regard the energy levels as a continuum. The region where the energy can be treated as a continuum is what is called the classical region.
| |
| | |
| === UV Catastrophy (Blackbody Radiation) ===
| |
| Imagine a perfect absorber cavity (i.e. it absorbs all radiation at all wavelengths, so that its spectral radiance is only going to depend on the temperature). This emission is called the blackbody radiation. This blackbody radiation experiment shows an important failure of classical mechanics. Lord Rayleigh (John William Strutt) and Sir James Jeans applied classical physics and assumed that the radiation in this perfect absorber could be represented by standing waves with nodes at the ends. The result predicts that the spectral intensity will increase quadratically with the increasing frequency, and will diverge to infinite energy or intensity squared at a UV frequency, or so called "Ultraviolet Catastophy."
| |
|
| |
|
| In 1900, Max Planck offered a successful explanation for blackbody radiation. He also assumed the the radiation was due to oscillations of the electron, but the difference between his assumption and Rayleigh's was that he assumed that the possible energies of an oscillator were not continuous. He proposed that the energy of this oscillator would be proportional to the frequency of a constant, the Planck constant.
| | '''Fall 2009 Midterm is October 15''' |
|
| |
|
| <math>E=h\nu=\hbar\omega</math>
| | == Outline of the Course == |
|
| |
|
| Here E is energy, h is the Planck constant (<math> h=6.626*10^{-34} Joule-seconds \!</math>
| | <b>Chapter 1: [[Physical Basis of Quantum Mechanics]]</b> |
| ) and <math>\nu\!</math> is the frequency of the oscillator. With the concept of energy being discrete in mind, the result is that Planck's calculation avoided the UV catastrophy, and instead the energy approached zero as the frequency is increased. In summary, the energy of the electromagnetic radiation is proportional to frequency instead of the amplitude of the electromagnetic waves, defying the classical physics where the energy is proportional to the intensity.
| |
|
| |
|
| === Photoelectric Effect ===
| | * [[Basic Concepts and Theory of Motion]] |
| Consider a system composed of light hitting a metal plate. From experimental observations, a current can be measured as light is incident on the metal plate. This phenomenon was first studied by Albert Einstein (1905). During this period, the classical point of view is that an electron is bound inside of an atom, and an excitation energy is needed in order to release it from the atom. This energy can be brought forth in the form of light. The classical point of view also includes the idea that the energy of this light is proportional to its intensity. Therefore, if enough energy (light) is absorbed by the electron, the electron eventually will be released. However, this is not the case.
| | * [[UV Catastrophe (Black-Body Radiation)]] |
| | * [[Photoelectric Effect]] |
| | * [[Stability of Matter]] |
| | * [[Double Slit Experiment]] |
| | * [[Stern-Gerlach Experiment]] |
| | * [[The Principle of Complementarity]] |
| | * [[The Correspondence Principle]] |
| | * [[The Philosophy of Quantum Theory]] |
|
| |
|
| The determining factor here is not the intensity of the light, but the frequency used on the electron. If the frequency of the light is the "specific" frequency, the electron will be released. This specific frequency of light is in resonance with the energy "frequency" of the electron. Einstein made the conclusion that the energy of a single photon is proportional to its frequency, not the intensity.
| |
|
| |
|
| Einstein realized that the classical view that light is a wave was not true, but instead light must be a particle. If light were a wave, then the mechanism that the wave transmits its energy is due to its perturbation, which is the amplitude of the wave. The photoelectric effect clearly shows that <math>(amplitude)^2 = intensity </math> has no affect on the electron energy, instead, only a specific frequency will have an effect on the energy of the electron. Agreeing with the UV catastrophy conclusions, he stated that light is made of particles called "photons," with an energy equal to hv.
| | <b>Chapter 2: [[Schrödinger Equation]]</b> |
| | |
| === Stability of Matter ===
| |
| One of the most important problems to inspire the creation of Quantum Mechanics was the model of the Hydrogen Atom. After Thompson discovered the electron, and Rutherford, the nucleus (or Kern, as he called it), the model of the Hydrogen atom was refined to one of the lighter electron of unit negative elementary charge orbiting the larger proton, of unit positive elementary charge. However, it was well known that classical electrodynamics required that charges accelerated by an EM field must radiate, and therefore lose energy. For an electron that moves in circular orbit about the more massive nucleus under the influence of the Coulomb attractive force, here is a simple non-relativistic model of this classical system:
| |
| | | |
| Where <math>r\,\!</math> is the orbital radius, and we neglect the motion of the proton by assuming it is much much more massive than the electron.
| | * [[Brief Derivation of Schrödinger Equation]] |
| | * [[Relation Between the Wave Function and Probability Density]] |
| | * [[Stationary States]] |
| | * [[Heisenberg Uncertainty Principle]] |
| | * [[Some Consequences of the Uncertainty Principle]] |
|
| |
|
| '''So the question is: What determines the rate <math>\rho</math> of this radiation? and how fast is this rate?'''
| |
|
| |
|
| The electron in the Bohr's model involves factors of: radius <math>r_0\,\!</math>, angular velocity <math>\omega\,\!</math>, charge of the particle <math>e\,\!</math>, and the speed of light, <math>c\,\!</math>: <math>\rho=\rho(r_0,\omega,e,c)\,\!</math>
| | <b>Chapter 3: [[Operators, Eigenfunctions, and Symmetry]]</b> |
|
| |
| The radius and charge will not enter separately, this is because if the electron is far from the proton, then the result can only depend on the dipole moment, which is .
| |
|
| |
|
| Therefore the above parameters is now:<math> \rho(er_0, \omega, c) \!</math>
| | * [[Linear Vector Spaces and Operators]] |
| | * [[Commutation Relations and Simultaneous Eigenvalues]] |
| | * [[The Schrödinger Equation in Dirac Notation]] |
| | * [[Transformations of Operators and Symmetry]] |
| | * [[Time Evolution of Expectation Values and Ehrenfest's Theorem]] |
|
| |
|
| '''What are the dimensions of <math>\rho\,\!</math>?'''
| |
|
| |
|
| Essentially, since light is energy, we are looking for how much energy is passed in a given time: <math>[\rho]=\frac{energy}{time} \!</math>
| | <b>Chapter 4: [[Motion in One Dimension]]</b> |
| | |
| Knowing this much already imposes certain constraints on the possible dimensions. By using dimensional analysis, let's construct something with units of energy.
| |
| | |
| From potential energy for coulombic electrostatic attractions: <math>energy=\frac{e^2 }{length} \!</math>
| |
| | |
| <math>e</math> has to be with <math>r_0</math> , multiply by <math>r^2</math> , and divide <math>length^2</math>.
| |
| | |
| The angular velocity is in frequency, so to get the above equations in energy/time, just multiply it with the angular velocity, <math>energy=\frac{e^2 }{length}\frac{r^2}{length^2}*\omega </math>
| |
| | |
| (Here, it is seen that the acceleration of the electron will increase with decreasing orbital radius. The radiation due to the acceleration a is given by the Larmor Formula: <math>energy \sim \frac{e^2r_0^2 }{(c/w)^3} w = \frac{e^2r_0^2 }{c^3}w^4\sim\frac{1}{r_0^4 } \!</math>
| |
| | | |
| It was known that the hydrogen atom had a certain radius on the order of .5 angstroms. Given this fact it can easily be seen that the electron will rapidly spiral into the nucleus, in the nanosecond scale. Clearly, the model depicts an unstable atom which would result in an unstable universe. A better representation of of an electron in an atom is needed.
| | * [[One-Dimensional Bound States]] |
| | | * [[Oscillation Theorem]] |
| ===Double Slit Experiment===
| | * [[The Dirac Delta Function Potential]] |
| | * [[Scattering States, Transmission and Reflection]] |
| | * [[Motion in a Periodic Potential]] |
| | * [[Summary of One-Dimensional Systems]] |
|
| |
|
| '''Bullet'''[[Image:Gun.png|thumb|125px|Double slit thought experiment with classical bullets]]
| |
|
| |
|
| Imagine a rattling gun which is shooting bullets in all directions. A histogram of the bullet's location after it passes through the two slits is plotted. If slit 2 is closed, but the slit 1 is open, then the green peak is observed which is given by the distribution function <math>p_1</math>. Similarly, if the slit 1 is closed, but he slit 2 is open, the pink peak is observed which is given by the distribution function <math>p_2</math>. When both slits are open, <math>peak_{12}</math> (purple) is observed. This agrees with the classical view, where the bullet is the particle and <math>p_{12}</math> is simply a sum of <math>p_1</math> and <math>p_2</math>.
| | <b>Chapter 5: [[Discrete Eigenvalues and Bound States - The Harmonic Oscillator and the WKB Approximation]]</b> |
|
| |
|
| The equation describing the probability of the bullet arrival if both of the slit are open is therefore
| | * [[Harmonic Oscillator Spectrum and Eigenstates]] |
| | * [[Analytical Method for Solving the Simple Harmonic Oscillator]] |
| | * [[Coherent States]] |
| | * [[Charged Particles in an Electromagnetic Field]] |
| | * [[WKB Approximation]] |
|
| |
|
| <math>p_{12}=p_1+p_2.</math>
| |
|
| |
|
| | <b>Chapter 6: [[Time Evolution and the Pictures of Quantum Mechanics]]</b> |
|
| |
|
| '''Classical Waves'''[[Image:waves.png|thumb|125px|right|Double slit thought experiment with water waves]]
| | * [[The Heisenberg Picture: Equations of Motion for Operators]] |
| | * [[The Interaction Picture]] |
| | * [[The Virial Theorem]] |
|
| |
|
|
| |
|
| As waves are passed through the double slit, the intensity of the waves which are proportional to the squares of the height of the wave motion <math>H1^2</math> and <math>H2^2</math> are observed when slit 1 and 2 are closed respectively. These intensities are similar to the histograms for the bullets in the previous demonstration. However, an interference pattern of the intensity (<math>H12</math>) is observed when both slits are opened. This is due to constructive and destructive interferences of the two waves. The resultant interference is the square of the sum of the two individual wave heights
| | <b>Chapter 7: [[Angular Momentum]]</b> |
| <math>H12 = (H1 +H2 )^2</math>
| |
| | |
| === Stationary states and Heisenberg Uncertainty relations ===
| |
| | |
| == Schrodinger equation and motion in one dimension ==
| |
| === The time dependent Schrodinger equation is===
| |
| | |
| <math>i\hbar \frac{\partial}{\partial t}|\psi\rangle=\mathcal{H}|\psi\rangle</math>
| |
| | |
| where the state <math>|\psi\rangle</math> evolves in time according to the Hamiltonian operator <math>\mathcal{H}</math>.
| |
| | |
| ===Motion in 1D===
| |
| '''Overview'''
| |
| | |
| Let's consider the motion in 1 direction of a particle in the potential V(x). Supposing that V(x) has finite values when x goes to infinity
| |
| :<math>\lim_{x \to -\infty}V(x)=V_-, \lim_{x \to +\infty}V(x)=V_+</math>
| |
| and assuming that: <math>V_-<V_+/!</math>
| |
| Schrodinger equation:
| |
| :<math>[-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+V(x)]\psi(x)=E\psi(x)</math>
| |
| <math>\rightarrow \frac{d^2}{dx^2}\psi(x)+\frac{2m}{\hbar^2}(E-V(x))\psi(x)=0</math>
| |
| | |
| From this equation we can discuss some general properties of 1-D motion as follows:
| |
| | |
| If <math>E>V_+\!</math>:
| |
| | |
| <math>E-V(x)>0\!</math> at both <math>-\infty</math> and <math>+\infty</math>. Therefore, the solution of Schrodinger equation are trigonometric function (sine or cosine). The wave function is oscillating at both <math>-\infty</math> and <math>+\infty</math>. The particle is in unbound state. The energy spectrum is continous. Both oscillating solutions are allowed, the energy level are two-fold degenerate.
| |
| | |
| If <math>V_-\le E \le V_+</math>:
| |
| | |
| <math>E-V(x)>0\!</math> at <math>-\infty</math> but <math>E-V(x)<0</math> at <math>+\infty</math>. Therefore, the wave function is oscillating at <math>-\infty</math> but decaying exponentially at <math>+\infty</math>. The energy spectrum is still continous but no longer degenerate.
| |
| | |
| If <math>E<V_-\!</math>:
| |
| | |
| <math>E-V(x)\!<0</math> at both <math>-\infty</math> and <math>+\infty</math>. Therefore, the wave function decays exponentially at both <math>-\infty</math> and <math>+\infty</math>. The particle is in bound state. The energy spectrum is discrete and non-degenerate.
| |
| | |
| == Operators, eigenfunctions, symmetry, and time evolution ==
| |
| | |
| === Commutation relations and simulatneous eigenvalues ===
| |
| | |
| COMMUTATORS
| |
| | |
| The commutator of two operators A and B is defined as follows:
| |
| | |
| <math>[A,B]=AB-BA\,\!.</math>
| |
| | |
| If <math>[A,B]=0</math>, we say that the operators <math>A</math> and <math>B</math> commute. Conversely, if , we say that the operators and do not commute. We can think of the commutator between two operators as a measure of how badly the two fail to commute.
| |
| | |
| Identities:<br/>
| |
| <math>
| |
| [A,B]+[B,A]=0 \!</math><br/>
| |
| <math>[A,A]= 0 \!</math><br/>
| |
| <math>[A,B+C]=[A,B]+[A,C]\!</math><br/>
| |
| <math>[A+B,C]=[A,C]+[B,C]\!</math><br/>
| |
| <math>[AB,C]=A[B,C]+[A,C]B\!</math><br/>
| |
| <math>[A,BC]=[A,B]C+B[A,C]\!</math><br/>
| |
| <math>[A,[B,C]]+[C,[A,B]]+[B,[C,A]]=0\!</math><br/>
| |
| | |
| | |
| | |
| | |
| | |
| Some more complicated commutator identities can be found here
| |
| http://sites.google.com/site/phy5645fall2008/some-useful-commutator-identites
| |
| | |
| COMPATIBLE OBSERVABLES
| |
| | |
| An operator which corresponds to some physically measurable property of a system is called an observable. All observables are Hermitian. If two Obersvables have simultaneous eigenkets, meaning for two obervables A and B,
| |
| | |
|
| |
| | |
| Then we have that, also, .
| |
| | |
| So we can see that,
| |
| | |
| The same logic works in reverse. So if two observables, A & B commute, i.e. if , then they have simultaneous eigenkets, and they are said to be compatible observables. Conversely, if , we say that the operators and do not commute and are incompatible observables.
| |
| | |
| | |
| POSITION AND MOMENTUM OPERATORS
| |
| An extremely useful example is the commutation relation of and . In 1D position space, = x and .
| |
| | |
| Applying , to an arbitrary state ket we can see that:
| |
| | |
| | |
| | |
| | |
| | |
| The position and momentum operators are incompatible. This provides a fundamental contrast to classical mechanics in which x and p obviously commute.
| |
| | |
| We can further generalize this relation between a position coordinate with its own and other conjugate momenta with what are known as the canonical commutation relations:
| |
| | |
| | |
| | |
| | |
| | |
| | |
| where the indices stand for x, y, or z.
| |
| | |
| | |
| | |
| CONNECTION BETWEEN CLASSICAL MECHANICS AND QUANTUM MEACHNICS
| |
| | |
| There is a wonderful connection between Classical mechanics and Quantum Mechanics. The Hamiltonian is a concept in the frame of classical mechanics. In this frame, the Hamiltonian is defined as:
| |
| | |
| | |
| | |
| | |
| There are two possibilities.
| |
| | |
| 1. If the Lagrangian does not depend explicitly on time the quantity H is conserved.
| |
| 2. If the Potential and the constraints of the system are time independent, then H conserves and also H is the Energy of the system.
| |
| | |
| It's clear from the above equation that:
| |
| | |
| | |
| | |
| | |
| | |
| | |
| This pair of the equations is called Hamilton's equations of motions. now let's define the following classical commutator:
| |
| | |
| | |
| | |
| | |
| This commutator is called Poisson Bracket, and it has interesting properties. To see, let's calculate commutation relationships between coordinates and momenta.
| |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| This relations clearly shows how close are the quantum commutators with classical world. If we perform the following identification:
| |
| | |
| | |
| | |
| | |
| Then we get quantum commutators. This identifications is called canonical quantization. As a final an important remark, the fact that we have classical commutators doesn't mean that we will have Heisenberg uncertainty relation for conjugate classical variables. This is because in classical mechanics the object of study are points (or body as a collection of points). In quantum mechanics we object of study is the state of a particle or system of particles - which describes the probability of finding a particle, and not it's exact, pointlike, location of momentum.
| |
| | |
| | |
| HAMILTONIAN
| |
| | |
| In Quantum Mechanics, an important property that one needs to check on a given operator (let's say ) is if it commutes with the Hamiltonian . If commutes with , then the eigenfunctions of can always be chosen to be eigenfunctions of . If commutes with the Hamiltonian and does not explicitly depend on time, then is a constant of motion.
| |
| | |
| | |
| COMMUTATORS & SYMMETRY
| |
| | |
| We can define an operator called the parity operator, which does the following:
| |
| | |
| | |
| | |
| | |
| | |
| The parity operator commutes with the Hamiltonian if the potential is symmetric, . Since the two commute, the eigenfunctions of the Hamiltonian can be written as eigenfunctions of the parity operator. This means that if the potential is symmetric, the solutions will be linear combinations of functions which have definite parity (even and odd functions).
| |
| | |
| GENERALIZED HEISENBERG UNCERTAINTY RELATION
| |
| | |
| === Heisenberg and interaction picture: Equations of motion for operators ===
| |
| === Feynman path integrals ===
| |
| | |
| == Discrete eigenvalues and bound states. Harmonic oscillator and WKB approximation ==
| |
|
| |
| Harmonic oscillator spectrum and eigenstates
| |
| Coherent states
| |
| Feynman path integral evaluation of the propagator
| |
| Motion in magnetic field
| |
| WKB
| |
| | |
| == Angular momentum ==
| |
|
| |
| Commutation relations, angular momentum as generator of rotations in 3D, eigenvalue quantization
| |
| Orbital angular momentum eigenfunctions
| |
| | |
| == Central forces ==
| |
|
| |
| ===Free particle in spherical coordinates===
| |
| A free particle is a specific case when <math>V_0=0\!</math> of the motion in a uniform potential . So it's more useful to consider a particle moving in a uniform potential. The Schrodinger equation for the radial part of the wave function is:
| |
| | |
| there is a sign error here (I will correct soon)
| |
| | |
| let , this should be k^2 not k. Rearranging the equation gives
| |
| | |
| Letting gives the equation:
| |
| :<math>(-\frac{\partial^2}{\partial\rho^2}+\frac{l(l+1)}{\rho^2})u_l(\rho)=u_l(\rho)=d_ld_l^+u_l(\rho)</math>
| |
| where and become the raising and lowering operators:
| |
|
| |
| Being <math>d_l^{\dagger}d_l=d_{l+1}d_{l+1}^{\dagger}, it can be shown that
| |
| | |
| | |
| For l=0, <math>\frac{partial^2}{partial \rho^2 u_o(\rho)u_o(\rho)
| |
| :<math>u_0=A\sin(\rho)+B\cos(\rho)</math>
| |
| Now applying the raising operator to the ground state
| |
| :<math>d_o^\dagger u_o(\rho)=(-\frac{\partial}{\partial\rho}+\frac{l+1}{\rho})u_o(\rho)=c_o u_1(\rho)</math>
| |
| | |
| ===Hydrogen atom===
| |
| | |
| == Continuous eigenvalues and collision theory ==
| |
| ===Differential cross-section and the Green's function formulation of scattering===
| |
| Much of what we know about forces and interactions in atoms and nuclei has been learned from scattering experiments, in which say atoms in the target are bombarded with beams of particles. These particles are scattered by the target atoms and then detected as a function of a scattering angle and energy. From theoretical point of view, we are now concerned with the continuous part of the energy spectrum. We are free to choose the value of the incident particle energy and by a proper choice of the zero of energy, this corresponds to E > 0 an to eigenfunctions of the unbound states. Before, when we were studying the bound states, the focus was on the discrete energy eigenvalues which allow a direct comparison of theory and experiments. In the continuous part of the spectrum, as it comes into play in scattering, the energy is given by the incident beam, and intensities are the object of measurement and prediction. These being the measures of the likelihood of finding a particle at certain places, are of course related to the eigenfunctions, rather than eigenvalues. Relating observed intensities to calculated wave functions is the first problem in scattering theory.
| |
| | |
| | |
| | |
| If <math>I_0</math> is the number of particles incident from the left per unit area per time and <math>Id\Omega</math> the number of those scattered into the cone per time, then the differential cross section is defined as
| |
| :
| |
| | |
| ===Central potential scattering and phase shifts===
| |
| Recall that for scattering we have Green function method
| |
| :<math>\psi_k =\psi_k^0 +\int d^3 r' G(r,r')V(r')\psi_k (r')</math>
| |
| where the Green's function satisfies that
| |
| :<math>(\nabla^2+k^2)G_k(\mathbf{r},\mathbf{r'})=\delta(\mathbf{r}-\mathbf{r'})</math>
| |
| and the solution is chosen such that the second term in Eq.(1) corresponds to an outgoing
| |
| wave.
| |
| :<math>G_k(\mathbf{r},\mathbf{r'})=-\frac{1}{4\pi}\frac{e^ { ik| \mathbf{r}-\mathbf{r'}|}}{|\mathbf{r}-\mathbf{r'}|}</math>
| |
| and in the asymptotic limit of r come to infinity
| |
| :<math>\lim_{r \to \infty}| \mathbf{r}-\mathbf{r'}|=r-\mathbf{r}\cdot\mathbf{r'}</math>
| |
| Thus
| |
| :<math>
| |
| \begin{array}{lcl}
| |
| \lim_{r \to \infty} \psi_k(\mathbf{r})&=&\psi_k^0(\mathbf{r}) -\frac{m}{2\pi\hbar^2 }\int d^3\mathbf{r'}e^ { -ik\mathbf{\hat{r}}\cdot\mathbf{r'}}V(\mathbf{r'})\psi_k(\mathbf{r'})\frac{e^ { ik\mathbf{r}}}{r}\\
| |
| &=&\psi_k^0(\mathbf{r})+f_k(\theta,\phi)\frac{e^ { ik\mathbf{r}}}{r}
| |
| \end{array}
| |
| </math>
| |
| where
| |
| :<math>f_k(\theta,\phi)=-\frac{m}{2\pi\hbar^2 }\int d^3\mathbf{r'}e^ { -ik\mathbf{\hat{r}}\cdot\mathbf{r'}}V(\mathbf{r'})\psi_k(\mathbf{r'})</math>
| |
| and the angles <math>\theta</math> and <math>\phi</math> are the angles between (the vector defining the detector) <math>\hat{r}</math> and <math>\mathbf{k}</math>
| |
| (the vector defining the in the incoming wave).
| |
| | |
| For central potentials, i.e. if <math>V(r)=V(|\mathbf{r}|)</math>, then <math>f_k(\theta)\!</math>, i.e. the scattering amplitude does
| |
| not depend on the azimuthal angle <math>\phi</math>.
| |
| To determine <math>f_k(\theta)\!</math> we need to find the solution of the Schrodinger equation:
| |
| | |
| :<math>(-\frac{\hbar^2 }{2m}\nabla^2+V(|\mathbf{r}|) )\psi=\frac{\hbar^2 k^2 }{2m}\psi</math>
| |
| | | |
| we use spherical coordinates
| | * [[Commutation Relations]] |
|
| | * [[Angular Momentum as a Generator of Rotations in 3D]] |
| :<math>(-\frac{\hbar^2 }{2m}\frac{\partial^2 }{\partial r^2 }+\frac{h^2 l(l+1)}{2mr^2}+V(|\mathbf{r}|) )u_l(r) =\frac{\hbar^2 k^2 }{2m}u_l(r)</math>
| | * [[Spherical Coordinates]] |
| | | * [[Eigenvalue Quantization]] |
| For V(r) with a finite range d, we have shown that for r >> d we have
| | * [[Orbital Angular Momentum Eigenfunctions]] |
| :<math>(-\frac{\hbar^2 }{2m}\frac{\partial^2 }{\partial r^2 }+\frac{\hbar^2 l(l+1)}{2mr^2} )u_l(r) =\frac{\hbar^2 k^2 }{2m}u_l(r)</math>
| |
| | |
| and the solution is a combination of spherical Bessel and Neumann functions
| |
|
| |
|
| :<math>\frac{u(r)}{r}=A_l j^l(kr) +B_l n_l(kr)</math>
| |
| when r is large enough we use approximation of Bessel function and Neumann function.
| |
|
| |
|
| :<math>\frac{u(r)}{r}\rightarrow A_l \frac{sin(kr-l\frac{\pi}{2})}{kr} -B_l \frac{cos(kr-l\frac{\pi}{2})}{kr}</math>
| | <b>Chapter 8: [[Central Forces]]</b> |
|
| |
| Letting
| |
| :<math>\frac{B_l }{A_l }=-tan\delta_l</math> | |
| | | |
| here <math>\delta_l</math> is called phase shift.
| | * [[General Formalism]] |
| we can rewrite the above expression (up to a normalization constant) as
| | * [[Free Particle in Spherical Coordinates]] |
| | * [[Spherical Well]] |
| | * [[Isotropic Harmonic Oscillator]] |
| | * [[Hydrogen Atom]] |
| | * [[WKB in Spherical Coordinates]] |
|
| |
|
| <math>\frac{u_l(r) }{r}\rightarrow\frac{sin(kr-l\frac{\pi}{2})}{kr}</math>
| |
|
| |
|
| Now since we are seeking the scattering amplitude with azimuthal symmetry, we can write
| | <b>Chapter 9: [[The Path Integral Formulation of Quantum Mechanics]]</b> |
| the solution of the Schrodinger equation as a superposition of <math>m=0\!</math> spherical harmonics
| |
| only:
| |
| :<math>\psi=\sum_{l=0}^{\mathop{ \infty}}a_l(k)p_l(cos\theta) \frac{u_l(r) }{r}</math>
| |
| :<math>\psi=\sum_{l=0}^{\mathop{ \infty}}a_l(k)p_l(cos\theta) \frac{sin(kr-l\frac{\pi}{2}+\delta_l )) }{kr}</math>
| |
| where the Legendre polynomials are
| |
| :<math>p_l(x)= \frac{1}{2^ll!}\frac{d^l }{dx^l }(x^2-1 )^l</math>
| |
| :<math>p_0(x)=1;p_1(x)=x;p_2(x)=\frac{1}{2}(3x^2 -1)</math>
| |
| let us fix the coeffcients <math>a_l(k)</math> by
| |
| :<math>e^{ikcos\theta}+f_k(\theta,\phi)\frac{e^ { ik\mathbf{r}}}{r}=\sum_{l=0}^{\mathop{ \infty}}a_l(k)p_l(cos\theta) \frac{sin(kr-l\frac{\pi}{2}+\delta_l )) }{kr}</math>
| |
| which must hold at large r and where we chose the coordinates by letting the incident wave
| |
| propagate along z-direction. Note that (due to an entirely separate argument):
| |
| :<math>e^{ikrcos\theta}=\sum_{l=0}^{\mathop{ \infty}}(2l+1)i^l p_l(cos\theta) \frac{sin(kr-l\frac{\pi}{2} ) }{kr}</math>
| |
| so
| |
| :<math>\sum_{l=0}^{\mathop{ \infty}}(2l+1)i^l p_l(cos\theta) \frac{sin(kr-l\frac{\pi}{2} ) }{kr}+f_k(\theta,\phi)\frac{e^ { -ik\mathbf{r}}}{r}=\sum_{l=0}^{\mathop{ \infty}}a_l(k)p_l(cos\theta) \frac{sin(kr-l\frac{\pi}{2}+\delta_l )) }{kr}</math>
| |
| We fix the coefficients <math>a_l(k)</math> by matching the incoming spherical waves on both sides of the
| |
| above equation. Note that this does not involve <math>f_k(\theta)</math> since the scattering amplitude
| |
| controls the outgoing spherical wave. Since <math>\cos x=(e^{ix} -e^{-ix})/2\!</math> we have
| |
| :<math>a_l(k)=(2l+1)i^le^{i\delta_l}</math>
| |
| Therefore:
| |
| :<math>f_k(\theta)=\frac{1}{k}\sum_{l=0}^{\infty}(2l+1)e^{i\delta_l }sin\delta_lp_l(cos\theta)</math>
| |
| Note that <math>\delta_l\!</math> is a function of k and therefore a function of the incident energy. If <math>\delta_l\!</math> is
| |
| known we can reconstruct the entire scattering amplitude and consequently the differential
| |
| cross section. The phase shifts must be determined from the solution of the Schrodinger
| |
| equation.
| |
|
| |
|
| Physically, we expect <math>\delta_l < 0\!</math> for repulsive potentials and <math>\delta_l > 0\!</math> for attractive potentials.
| | * [[Feynman Path Integrals]] |
| Also, if l/k >> d, then the classical impact parameter is much larger than the range of the
| | * [[The Free-Particle Propagator]] |
| potential and in this case we expect <math>\delta_l\!</math> to be small.
| | * [[Propagator for the Harmonic Oscillator]] |
|
| |
|
| The differential scattering cross section is
| |
| :<math>\frac{d\sigma}{d\Omega}=|f_k(\theta) |^2=\frac{1}{k^2 }|\sum_{l=0}^{\infty}(2l+1)e^{i\delta_l }sin\delta_lp_l(cos\theta) |^2</math>
| |
| By integrating <math>\frac{d\sigma}{d\Omega}</math> over the solid angle we obtain the total scattering cross section
| |
| :<math>\sigma_{tot}=\frac{4\pi}{k^2 }\sum_{l=0}^{\mathop{ \infty}}(2l+1)sin^2\delta_l</math>
| |
| which follows from the orthogonality of the Legendre polynomials
| |
|
| |
|
| :<math>\int_{-1}^{1}dxp_l(x) p_l'(x)=\frac{2}{(2l+2)}\delta_{ll }</math>
| | <b>Chapter 10: [[Continuous Eigenvalues and Collision Theory]]</b> |
| | | |
| Finally note that since <math>p_l(1) = 1\!</math> for all <math>l\!</math>, we have
| | * [[Differential Cross Section and the Green's Function Formulation of Scattering]] |
| :<math>\sigma_{tot}=\frac{4\pi}{k}\Im mf(0)</math>
| | * [[Central Potential Scattering and Phase Shifts]] |
| here we take the imaginary part.
| | * [[Coulomb Potential Scattering]] |
| This relationship is known as the optical theorem.
| |
| | |
| ===Born approximation and examples of cross-section calculations===
| |
| Recall the scattering of a particle in a potential V(r) has a differential cross section of:
| |
| | |
| | |
| The scattering amplitude, fk(r), is the coefficient of the outgoing wave.
| |
| | |
| | |
| | |
| The Born approximation, often called the first Born approximation, is a technique for find solutions when V(r) is small.
| |
| | |
| The scattering amplitude can be approximated by:
| |
| where -ikr is the scattered portion and ik is the incident portion.
| |
| | |
| The scattering amplitude is defined as the coefficient of the outgoing wave in the asymptotic solution (for large r)
| |
| | |
| of the Schrodinger equation:
| |
| | |
| For a central-force potential, V(r), the Born scattering amplitude reduces to
| |
| | |
| | |
| | |
| Which leads to the Born cross section:
| |
| | |
| | |
| | |
| | |
| Born Approximation for Spherically Symmetric Potentials:
| |
| | |
| | |
| | |
| Given spherically symmetry we may define
| |
| | |
| and align the polar axis for the r integral lie along this quantity. We then have:
| |
| | |
| Our first Born integration then takes the form:
| |
| | |
| | |
| The phi integral introduces a trivial 2pi. For the theta integral we can use the following identity:
| |
| | |
| | |
| and get:
| |
| | |
| | |
| where the angular dependence of f is carried by kappa:
| |
| | |
| | |
| | |
| Example 1:
| |
| | |
| | |
| Consider the scattering amplitude from a Gaussian potential of the form:
| |
| | |
| | |
| Our scattering amplitude then becomes:
| |
| | |
| | |
| | |
| | |
| Integration by parts gives:
| |
| | |
| ===Coulomb potential scattering===
| |
| Example 1
| |
| | |
| | |
| Lets look at an example of a Screened Coulomb (Yukawa) Potential, where we have a potential:
| |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| thus we have the differential cross section:
| |
| | |
| | |
| | |
| | |
| Example 2
| |
| | |
| When we are considering scattering due to the Coulomb potential, we can not neglect the effect of this potential at large distances because it is only a 1/r potential.
| |
| | |
| Use a change of coordinates from Cartesian to parabolic coordinates:
| |
| | |
| ,
| |
| ,
| |
| | |
| | |
| The following is a picture of parabolic coordinates:
| |
| | |
| Phi represents rotation about the z-axis, xi represents the parabolas with their vertex at a minimum, and eta represents parabolas with their vertex at a maximum.
| |
| | |
| So now we can write the Schrodinger equation in parabolic coordinates:
| |
| | |
| | |
| | |
| So we will seek solutions which are independent of phi. Recall that the scattering amplitude is a function of only.
| |
| | |
| Look for solutions of the form:
| |
| | |
| | |
| We can tidy up the notation a little bit by using the following substitution:
| |
| | |
| Now let:
| |
| From this we can write:
| |
| | |
| | |
| | |
| Recall the confluent hypergeometric function:
| |
| | |
| | |
| We can then write the recursion formula as the following:
| |
| | |
| | |
| | |
| This implies that:
| |
| | |
| | |
| where the confluent geometric function is written in terms of three new variables, and A is a c-number.
| |
| | |
| Now we can write the wavefunction due to Coulomb scattering:
| |
| | |
| | |
| | |
| Now we should look at the limit where z is taken to go to infiniti, and our confluent hypergeometric function is rewritten as:
| |
| | |
| | |
| | |
| Now we can use this to rewrite our equation for Phi of xi:
| |
| | |
| | |
| Rewriting our wavefunction :
| |
| | |
| | |
| | |
| where f(theta) is the following:
| |
| | |
| | |
| | |
| We can then get our differential cross section from that by squaring it:
| |
| | |
| | |
| | |
| If we normalize the wavefunction to give unit flux at large distances, we must take the following for the constant C:
| |
| | |
| | |
| | |
| So the wavefunction at large distances is given by the following:
| |
| | |
| | |
| | |
| where
| |
| | |
| | |
| | |
| Plugging this in for our wavefunction squared:
| |
| | |
| | |
| | |
| Now let's use the following quantity to represent the velocity:
| |
| | |
| | |
| For small incident velocities, we can write:
| |
| | |
|
| |
| | |
| | |
| | |
| where the first equation is for an attractive Coulomb potential, and the second equation is for a repulsive Coulomb potential.
| |
| | |
| ===Two particle scattering===
| |
| Classically, if we wish to consider a collection of identical particles, say billiard balls, it is always possible to label all the balls such that we can follow a single ball throughout interacting with others. We could, for example, label each ball with a different color. Then, after an arbitrary number of interactions, we can distinguish, say, a red ball from any other. It is not, however, possible to attach such labels to quantum mechanical systems of, say, electrons. Quantum mechanical particles are far to small to attach such physical labels and there are not enough degrees of freedom to label each particle differently. Again considering the classical case, one could imagine simply recording the position of a given particle throughout its trajectory to distinguish it from any other particle. Quantum mechanically, however, we again fail in following a single particles trajectory since each time we make a measurement of position we disturb the system of particles in some uncontrollable fashion. If the wave functions of the particles overlap at all, then the hope of following a single particles trajectory is lost. We now attempt to study the consequences of such indistinguishably between identical quantum particles.
| |
| | |
| | |
| Scattering of Identical Particles | |
| | |
| Lets look at the case of two identical bosons (spin 0) from their center of mass frame. To describe the system, we must use a symmetrized wave function. Under the exchange r1 <----> r2, rcm = (r1 + r2)/2 is invariant while r = r1-r2 changes sign. So the center of mass wave function is already symmetric. Furthermore, the wave function has even parity. This implies that the only possible eigenstates of angular momentum of the two particles are those with even angular momentum quantum numbers. This is evident from the property of the associated Legendre polynomials:
| |
| | |
| | |
| | |
| But we have to symmetrize Psi(r) by hand:
| |
| | |
| | |
| Under the transformation:
| |
| | |
| The first two terms of the symmetrized wave function represent the incident waves corresponding to the center of mass frame. Note that because we are considering identical particles we cannot distinguish the target particle from the incident one. Thus, each particle has equal amplitude of being either one.
| |
| | |
| The scattering amplitude is:
| |
| | |
| | |
| and can then be associated with the angle through which each particle is scattered. The total amplitude for particles to emerge at each angle is then exactly the sum of amplitudes for emerging at each angle, which is given above. The scattering amplitude remains consistent with the fact that we have two identical particles, and this gives us the differential cross section:
| |
| | |
| | |
| | |
| Note:
| |
| | |
| The first two terms in the differential cross section is what we would get if we had two distinguishable particles, while the third term give the quantum mechanical interference that goes along with identical particles.
| |
| | |
| As an example, consider scattering through a 90 degree angle. We then have:
| |
| | |
| | |
| Now if the particles are distinguishable, the cross section for observing a scattered particle at 90 degrees is then:
| |
| | |
| | |
| | |
| Where if the particles are indistinguishable, we see above that we will have:
| |
| | |
| | |
| | |
| Thus the differential cross-section is exactly twice the distinguishable case when the particles are indistinguishable.
| |