Chemistry 365: Normal Mode Analysis David Ronis McGill University
|
|
- Virginia Preston
- 5 years ago
- Views:
Transcription
1 Chemistry 365: Normal Mode Analysis David Ronis McGill University 1. Quantum Mechanical Treatment Our starting point is the Schrodinger wave equation: Σ h 2 2 2m i N i=1 r 2 i + U( r 1,..., r N ) Ψ( r 1,..., r N ) = E Ψ( r 1,..., r N ), (1.1) where N is the number of atoms in the molecule, m i is the mass of the i th atom, and U( r 1,..., r N )isthe effective potential for the nuclear motion, e.g., as is obtained in the Born- Oppenheimer approximation. If the amplitude of the vibrational motion is small, then the vibrational part of the Hamiltonian associated with Eq. (1.1) can be written as: H vib Σ N h 2 2 2m i i=1 2 i + U N Σ K i, j : i j, (1.2) i, j=1 R i is the equilibrium posi- where U is the minimum value of the potential energy, i r i R i, tion of the i th atom, and K i, j 2 U (1.3) r i r j r k = R k is the matrix of (harmonic) force constants. Henceforth, we will shift the zero of energy so as to make U =. Note that in obtaining Eq. (1.2), we have neglected anharmonic (i.e., cubic and higher order) corrections to the vibrational motion. The next and most confusing step is to change to matrix notation. We introduce a column vector containing the displacements as: [ x 1, y 1, z 1,..., x N, y N, z N ]T, (1.4) where "T "denotes a matrix transpose. Similarly, we can encode the force constants or masses into 3N 3N matrices, and thereby rewrite Eq. (1.2) as H vib = h 2 M K, (1.5) where all matrix quantities are emboldened, the superscript " " denotes the Hermitian conjugate (matrix transpose for real matrices), and where the mass matrix, M, is a diagonal matrix with the masses of the given atoms each appearing three times on the diagonal. The Hamiltonian given by Eq. (1.5) is the generalization of the usual harmonic oscillator Hamiltonian to include more particles and to allow for "springs" between arbitrary particles. Unless K is diagonal (and it usually isn t) Eq. (1.5) would seem to suggest that the vibrational problem for a polyatomic is non-separable (why?); nonetheless, as we now show, it can be separated. This is first shown using using the quantum mechanical framework we ve just set up.
2 Chemistry Normal Mode Analysis Later, we will discuss the separation using classical mechanics. This is valid since, for harmonic oscillators, you get the same result for the vibration frequencies. We now make the transformation i m 1/2 i i,which allows Eq. (1.5) to be reexpressed in the new coordinates as where H vib = h K, (1.6) K M 1/2 KM 1/2. (1.7) Since the matrix K is Hermitian (or symmetric for real matrices), it is possible to find a unitary (also referred to as an orthogonal matrix for real matrices) matrix which diagonalizes it; i.e., you can find a matrix P which satisfies where λ is diagonal (with real eigenvalues) and where P KP = λ or K = PλP, (1.8) PP = 1 (1.9) where 1 is the identity matrix. Note that P is the unitary matrix who s columns are the normalized eigenvectors. (See a good quantum mechanics book or any linear algebra text for proofs of these results). We now make the transformation in Eq. (1.5). This gives P (1.1) H vib = h 2 P P λ, (1.11) Finally, Eq. (1.8) allows us to cancel the factors of P in this last equation and write H vib = 3N Σ 1 2 h 2 i=1 2 2 i λ i 2 i. (1.12) Thus the two transformations described above have separated the Hamiltonian into 3N uncoupled harmonic oscillator Hamiltonians, and are usually referred to as a normal mode transformation. Remember that in general the normal mode coordinates are not the original displacements from the equilibrium positions, but correspond to collective vibrations of the molecule. 2. Normal Modes in Classical Mechanics The starting point for the normal mode analysis in classical mechanics are Newton s equations for a system of coupled harmonic oscillators. Still, using the notation that led to Eq. (1.12), we can write the classical equations of motion as: M (t) = K (t), (2.1) where components of the left-hand-side of the equation are the rates of change of momentum of the nuclei, while the right-hand-side contains the harmonic forces.
3 Chemistry Normal Mode Analysis Since we expect harmonic motion, we ll look for a solution of the form: (t) = cos(ω t + φ)y, (2.2) where φ is an arbitrary phase shift. If Eq. (2.2) is used in Eq. (2.1) it follows that: (Mω 2 K)Y =, (2.3) where remember that Y is a column vector and M, and K are matrices. We want to solve this homogeneous system of linear equations for Y. In general, the only way to get a nonzero solution is to make the matrix [Mω 2 K] singular; i.e., we must set det(mω 2 K) = det(ω 2 M 1/2 KM 1/2 )det(m) =. (2.4) Since, det(m) = (M 1 M 2...M N ) 3 wesee that the second equality is simply the characteristic equation associated with the eigenvalue problem: Ku = λu, (2.5) with λ ω 2 and with K given byeq. (2.7) above. The remaining steps are equivalent to what was done quantum mechanically. 3. Force Constant Calculations Here is an example of a force constant matrix calculation. We will consider a diatomic molecule, where the two atoms interact with a potential of the form: U(r 1, r 2 ) K r 1 r 2 R ; (3.1) i.e., a simple Hookian spring. It is easy to take the various derivatives indicated in the preceding sections; here, however, we will explicitly expand the potential in terms of the atomic displacements, i. Bywriting, r i = R i + i (where R i is the equilibrium position of the i th nucleus), Eq., (3.1) can be rewritten as: U(r 1, r 2 ) = 1 1/2 2 2 K R R 12 ( 1 2 ) R, (3.2) where R 12 R 1 R 2. Clearly, the equilibrium will have R 12 = R. Moreover, we expect that the vibrational amplitude will be small, and thus, the terms in the s in the square root in Eq. (3.2) will be small compared with the first term. The Taylor expansion of the square root implies that we can write which can be rewritten as A + B A + U(r 1, r 2 ) = 1 2 K R R 12 ( 1 2 ) B +..., (3.3) 2 A 2R 12 2 R, (3.4) U(r 1, r 2 ) = 1 2 K[ ˆR 12 ( 1 2 )] 2, (3.5)
4 Chemistry Normal Mode Analysis where the ˆ denotes a unit vector and where all terms smaller than quadratic in the nuclear displacements have been dropped. If the square is expanded, notice the appearance of cross terms in the displacements of 1 and 2. It is actually quite simple to finish the normal mode calculation in this case. To do so, define the equilibrium bond to point in along the x axis. Equation (3.5) shows that only x-components of the displacements cost energy, and hence, there will no force in thy y or z directions (thereby resulting in 4 zero eigen-frequencies). For the x components, Newton s equations become: m 1 m 2 1 x = K K x 2 K K 1 x 2 x, (3.6) where m i is the mass of the i th nucleus. This in turn leads to the following characteristic equation for the remaining frequencies: and = det m 1 m 2 ω 2 K K = (m 1 ω 2 K)(m 2 ω 2 K) K 2 = ω 2 (µω 2 K), where µ m 1 m 2 /(m 1 + m 2 )isthe reduced mass. Thus we pick up another zero frequency and a nonzero one with ω = K/µ, which is the usual result. (Note that we don t count ± roots twice--why?). 4. Normal Modes in Crystals The main result of the preceding sections is that the characteristic vibrational frequencies are obtained by solving for eigenvalues, cf., Eq. (2.5). For small to mid-size molecules this can be done numerically, onthe other hand, this is not practical for crystalline solids where the matrices are huge (3N 3N with N O(N A ))! Crystalline materials differ from small molecules in one important aspect; they are periodic structures made up of identical unit cells; specifically, each unit cell is placed at K K R = n 1 a 1 + n 2 a 2 + n 3 a 3, (4.1) where n i is an integer and a i is a primitive lattice vector. The atoms are positioned within each unit cell at positions labeled by an index, α. Thus, the position of any atom in entire crystal is determined by specifying R and α.hence, Eq. (2.5) can be rewritten as λu R,α = Σ R α Σ K R,α ;R,α u R,α. (4.2) The periodicity of the lattice implies that only the distance between unit cells can matter; i.e., K R,α ;R,α = K R R ;α,α,and this turns Eq. (4.2) into a discrete convolution that can be simplified by introducing a discrete Fourier transform, i.e., we let ũ k,α Σ e ik R u R,α. (4.3) R
5 Chemistry Normal Mode Analysis When applied to both sides of Eq. (2.5), this gives * where λũ k,α K k;α,α Σ K k;α,α ũ k,α, (4.4) α Σ e ik R K R;α,α. (4.5) R The resulting eigenvalue problem has rank 3N cell,where N cell is the number of atoms in a unit cell, and is easily solved numerically. We haven t specified the values for the k s. Itturns out that a very convenient choice is to use k s expressed in terms of the so-called reciprocal lattice vectors, G,i.e., G m 1 b 1 + m 2 b 2 + m 3 b 3, (4.6) where the m i s are integers, and where the reciprocal lattice basis vectors are defined as a b 1 2 a 3 2π a 1 ( a 2 a 3 ), a b 2 3 a 1 2π a 2 ( a 3 a 1 ), and b a 3 1 a 2 2π a 3 ( a 1 a 2 ). (4.7) Note the following: a) The denominators appearing in the definitions of the b i are equal; i.e., a 1 ( a 2 a 3 ) = a 2 ( a 3 a 1 ) = a 3 ( a 1 a 2 ) = ν cell, (4.8) where ν cell is the volume of the unit cell. The volume of the reciprocal lattice unit cell is (2π ) 3 /ν cell. b) By construction, where δ i, j is a Kronecker-δ. a i b j = 2πδ i, j, (4.9) c) The real and reciprocal lattices need not be the same, even ignoring how the basis vectors are normalized; e.g., they are for the SCC, but the reciprocal lattice for the BCC lattice is a FCC lattice. d) For any lattice vector, R,and reciprocal lattice vector, G, e i R G = e 2π i(m 1n 1 +m 2 n 2 +n 3 m 3 ) = 1, (4.1) cf. Eqs. (4.1), (4.6), and (4.9). Hence, adding any reciprocal a lattice vector to the k in Eq. (4.4) changes nothing, and so, we restrict k to what is known as the First Brillouin Zone; specifically, k = k 1 b 1 + k 2 b 2 + k 3 b 3, (4.11) where k i m i /N i with m i =, 1, 2,..., N i 1. Note that N 1 N 2 N 3 = N, the total number of * Strictly speaking, in obtaining Eq. (4.5) we have assumed some special properties of the lattice. We don t use a finite lattice, rather one that obeys periodic boundary conditions (see below). Here we re using the definition of C. Kittel, Introduction to Solid State Physics (Wiley, 1966), p. 53. There are others, see, e.g., M. Born and K. Huang, Dyamical Theory of Crystal Lattices (Oxford, 1968), p. 69.
6 Chemistry Normal Mode Analysis cells in the crystal. The N i s are the number of cells in the a i direction. With this choice *, consider 1 N k Σ e i k Rũ k,α = 1 N k Σ Σ e i k ( R R) 3 u = R,α Σ u Π 1 R,α The sums in parenthesis are geometric series and give 1 N i N i 1 m i = R Σ e 2π im i(n i n i )/N i which when used in Eq. (4.12) shows that 1 N k R i=1 N i N i 1 Σ e 2π im i(n i n i )/N i. (4.12) m i = = 1 e 2π i(n i n ) i 1 N i e 2π i(n i n i )/N i = δ 1 ni,n i, (4.13) Σ e i k Rũ k,α = u R,α. (4.14) The sums over k are really sums over the m i s, and thus, k i hardly changes as we go from m i m i + 1for large N i,cf. Eq. (4.11). This allows the sums over m i to be replaced by integrals; i.e., Σ... k 3 Π N i 1 i=1 dm i...= where V = Nν cell is the volume of the system. V (2π ) d k..., (4.15) 3 1st Brillouin The main goal of this discussion is to obtain the vibrational density of states. By using Eq. (4.15), it is easy to show that g(ω ) = zone V (2π ) d k 3 δ (ω ω (k)), (4.16) 1st Brillouin zone where δ (x) is the Dirac δ -function. Once the frequencies are known, it is relatively easy to numerically bin them by frequency asafunction of k. 5. Normal Modes in Crystals: An Example Consider a crystal with one atom of mass m per unit cell and nearest neighbor interactions of the type considered in Sec. 3. Newton s equations of motion for this atom become m,, = K[ê x ( x 1,, + x 1,, 2 x,,) + ê y ( y,1, + y, 1, 2 y,, ) +ê z ( x,,1 + x,,1 2 x,,)], (5.1) where ê i is a unit vector in the i direction. Note that this model implies that vibrations in the x, y, z, directions are separable. By using this in Eq. (4.5) we see that * The more common choice for the range of the m i s is m i = N i /2,...,,..., N i /2, which, cf. Eq. (4.11), means that 1/2 k i 1/2.
7 Chemistry Normal Mode Analysis sin 2 (k 1 π ) K k = 4ω 2 sin 2 (k 2 π ) sin 2 (k 3 π ), (5.2) where ω (K/m) 1/2 and where the identity, 1 cos(x) = 2sin 2 (x/2) was used. Since K k is already diagonal, we see that ω i ( k) = 2ω sin(k i π ), i = 1, 2, 3, which is also the result for a one dimensional chain (as was expected, given the separability of the vibrations in the x, y, and z directions). The normal mode is just ê i,which shows that our model potential is too simple. Consider the first eigen-vector. Itcorresponds to an arbitrary displacement in the x direction, with an eigenvalue that is independent of k y and k z. When k x = wedon t get a single zero frequency, asexpected, but one for any ofthe N y N z values of (k y, k z ). Thus, with 3N 2/3 zero frequencies, not 6, the crystal is unstable! by The problem is easily fixed by modifying the interaction potential, i.e., we replace Eq. (3.5) U K 2 ( 1 2 ) 2, (5.3) which is invariant under rigid translations and rotations, and allows us to rewrite Eq. (5.1) as m,, = K( 1,, + 1,, 2,, +,1, +, 1, 2,, +,,1 +,,1 2,, ). (5.4) With this, it follows that K k = 4ω 2 [sin 2 (k 1 π ) + sin 2 (k 2 π ) + sin 2 (k 3 π )] 1, (5.5) where 1 is a 3 3identity matrix and where we have expressed k in terms of the reciprocal lattice basis vectors, cf. Eq (5.11). The triply degenerate vibrational frequencies are simply ω i ( k) = 2ω [sin 2 (k 1 π ) + sin 2 (k 2 π ) + sin 2 (k 3 π )] 1/2. (5.6) If all the k i, Eq. (5.6) becomes ω ( k) 2ω kπ,where k (k k k 2 3) 1/2 which is the expected linear dispersion law for long wevelength, acoustic vibrations. Some constant frequency surfaces in the First Brillouin Zone are shown in Figs
8 Chemistry Normal Mode Analysis Fig A constant normal mode frequency surface for ω /ω = 1. The reciprocal lattice basis was used which need not be orthogonal (although it is for the SCC lattice). Also note that we ve switched to the other definition of the first Brillouin zone, with 1/2 k i 1/2. The surface is roughly spherical, as expected for long wavelength acoustic phonons. Fig As in Fig. 5.1 but with ω /ω = 2. Notice the C 4 axises through the centers of any ofthe unit cell faces. Fig ω /ω = 3. As in Fig. 5.1 but with Normal mode frequencies for this model were computed numerically for a uniform sample of k i s in the First Brillouin Zone and binned in order to get the vibrational density of of states as shown in Fig. 5.4
9 Chemistry Normal Mode Analysis Fig Vibrational density of states for the isotropic model defined by Eq. (5.3). Any single atom per unit cell lattice will give the same result. The data was obtained by numerically binning the frequencies in the First Brillouin Zone. Note that the maximum frequency for this model is 2 3ω. By noting that where f (x) has zeros at x = x i, i = 1,..., and that we can rewrite Eq. (4.16) as g(ω ) 3N δ ( f (x)) = Σ δ (x x i), (5.7) i df (x i ) dx i δ (x) = ds 2π eixs, (5.8) = 2ω 1 dk 1 1 dk 2 1 dk 3 δ (ω 2 ω 2 (k)) = ω π = ω π ds 1 dk 1 1 dk 2 1 dk 3 e is(ω 2 ω 2 (k)) ds eisω 2 Φ 3 (s), where the extra factor of 3 is due to the triple degeneracy ofeach mode and where Φ(s) 1 dk e 4iω 2 s sin2 (π k) When we use Eq. (5.7) for the frequencies, it follows that the wavevector integrations separate and
10 Chemistry Normal Mode Analysis where g(ω ) = 2Vω Φ(s) dk π e is4ω sin 2 (π k ) 1 = e 2ω is 2π 2 ds 2π eisω 2 Φ 3 (s), (5.9) π dz e2ω 2 is cos(z) 1 2π dz e 2ω 2 isz (1 z 2 ) = e 2ω 1/2 2π = e 2ω 2 is is J (2ω 2 s), (5.1) where J (x) isabessel function of the first kind. * When this is used in Eq. (5.1), we see that g(ω ) = 2Vω (2π ) 4 ds eis(ω 2 6ω 2 ) J 3 (2ω 2 s) {5.11) * The integral leading to the last equality can be found in I.S. Gradsheyn and I.M. Ryzhik, Table of Integrals, Series, and Products, A.Jeffrey editor, (Academic Press, 198), Eq. ( ) on p. 321.
Phonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationPhonons I - Crystal Vibrations (Kittel Ch. 4)
Phonons I - Crystal Vibrations (Kittel Ch. 4) Displacements of Atoms Positions of atoms in their perfect lattice positions are given by: R 0 (n 1, n 2, n 3 ) = n 10 x + n 20 y + n 30 z For simplicity here
More informationClassical Theory of Harmonic Crystals
Classical Theory of Harmonic Crystals HARMONIC APPROXIMATION The Hamiltonian of the crystal is expressed in terms of the kinetic energies of atoms and the potential energy. In calculating the potential
More informationSolid State Physics. Lecturer: Dr. Lafy Faraj
Solid State Physics Lecturer: Dr. Lafy Faraj CHAPTER 1 Phonons and Lattice vibration Crystal Dynamics Atoms vibrate about their equilibrium position at absolute zero. The amplitude of the motion increases
More informationTopic 5-1: Introduction to Phonons Kittel pages: 91, 92
Topic 5-1: Introduction to Phonons Kittel pages: 91, 92 Summary: In this video we introduce the concept that atoms are not rigid, fixed points within the lattice. Instead we treat them as quantum harmonic
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationShort Sample Solutions to the Sample Exam for (3rd Year Course 6) Hilary Term 2011
Short Sample Solutions to the Sample Exam for (3rd Year Course 6) Hilary Term 0. [4] A Lattice is an infinite set of points in space where the environment of any given point is identical to the enivironment
More informationPhonons (Classical theory)
Phonons (Classical theory) (Read Kittel ch. 4) Classical theory. Consider propagation of elastic waves in cubic crystal, along [00], [0], or [] directions. Entire plane vibrates in phase in these directions
More informationChemistry 593: The Semi-Classical Limit David Ronis McGill University
Chemistry 593: The Semi-Classical Limit David Ronis McGill University. The Semi-Classical Limit: Quantum Corrections Here we will work out what the leading order contribution to the canonical partition
More informationPH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5
PH575 Spring 2009 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 Spring 2009 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 Spring 2009 POP QUIZ
More informationRotations and vibrations of polyatomic molecules
Rotations and vibrations of polyatomic molecules When the potential energy surface V( R 1, R 2,..., R N ) is known we can compute the energy levels of the molecule. These levels can be an effect of: Rotation
More informationIntroduction to Condensed Matter Physics
Introduction to Condensed Matter Physics The Reciprocal Lattice M.P. Vaughan Overview Overview of the reciprocal lattice Periodic functions Reciprocal lattice vectors Bloch functions k-space Dispersion
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More informationKronig-Penney model. 2m dx. Solutions can be found in region I and region II Match boundary conditions
Kronig-Penney model I II d V( x) m dx E Solutions can be found in region I and region II Match boundary conditions Linear differential equations with periodic coefficients Have exponentially decaying solutions,
More informationQuantum Condensed Matter Physics Lecture 5
Quantum Condensed Matter Physics Lecture 5 detector sample X-ray source monochromator David Ritchie http://www.sp.phy.cam.ac.uk/drp2/home QCMP Lent/Easter 2019 5.1 Quantum Condensed Matter Physics 1. Classical
More information! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»
Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable
More informationCoupled Oscillators Monday, 29 October 2012 normal mode Figure 1:
Coupled Oscillators Monday, 29 October 2012 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l m k θ 2
More informationGeometry of Crystal Lattice
0 Geometry of Crystal Lattice 0.1 Translational Symmetry The crystalline state of substances is different from other states (gaseous, liquid, amorphous) in that the atoms are in an ordered and symmetrical
More informationMP464: Solid State Physics Problem Sheet
MP464: Solid State Physics Problem Sheet 1) Write down primitive lattice vectors for the -dimensional rectangular lattice, with sides a and b in the x and y-directions respectively, and a face-centred
More informationeigenvalues eigenfunctions
Born-Oppenheimer Approximation Atoms and molecules consist of heavy nuclei and light electrons. Consider (for simplicity) a diatomic molecule (e.g. HCl). Clamp/freeze the nuclei in space, a distance r
More information5.1 Classical Harmonic Oscillator
Chapter 5 Harmonic Oscillator 5.1 Classical Harmonic Oscillator m l o l Hooke s Law give the force exerting on the mass as: f = k(l l o ) where l o is the equilibrium length of the spring and k is the
More informationIn which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid.
Coupled Oscillators Wednesday, 26 October 2011 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l k θ 2
More informationQuantum Theory of Matter
Imperial College London Department of Physics Professor Ortwin Hess o.hess@imperial.ac.uk Quantum Theory of Matter Spring 014 1 Periodic Structures 1.1 Direct and Reciprocal Lattice (a) Show that the reciprocal
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 9, February 8, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 9, February 8, 2006 The Harmonic Oscillator Consider a diatomic molecule. Such a molecule
More informationElectrons in a weak periodic potential
Electrons in a weak periodic potential Assumptions: 1. Static defect-free lattice perfectly periodic potential. 2. Weak potential perturbative effect on the free electron states. Perfect periodicity of
More informationVibrations and Rotations of Diatomic Molecules
Chapter 6 Vibrations and Rotations of Diatomic Molecules With the electronic part of the problem treated in the previous chapter, the nuclear motion shall occupy our attention in this one. In many ways
More informationMP464: Solid State Physics Problem Sheet
MP464: Solid State Physics Problem Sheet 1 Write down primitive lattice vectors for the -dimensional rectangular lattice, with sides a and b in the x and y-directions respectively, and a face-centred rectangular
More information5.80 Small-Molecule Spectroscopy and Dynamics
MIT OpenCourseWare http://ocw.mit.edu 5.80 Small-Molecule Spectroscopy and Dynamics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.80 Lecture
More informationStructure and Dynamics : An Atomic View of Materials
Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic
More informationOpinions on quantum mechanics. CHAPTER 6 Quantum Mechanics II. 6.1: The Schrödinger Wave Equation. Normalization and Probability
CHAPTER 6 Quantum Mechanics II 6.1 The Schrödinger Wave Equation 6. Expectation Values 6.3 Infinite Square-Well Potential 6.4 Finite Square-Well Potential 6.5 Three-Dimensional Infinite- 6.6 Simple Harmonic
More informationReciprocal Lattice. Let A(r) be some function featuring discrete translation symmetry:
Reciprocal Lattice From Quantum Mechanics we know that symmetries have kinematic implications. Each symmetry predetermines its specific quantum numbers and leads to certain constraints (conservation laws/selection
More informationCh 125a Problem Set 1
Ch 5a Problem Set Due Monday, Oct 5, 05, am Problem : Bra-ket notation (Dirac notation) Bra-ket notation is a standard and convenient way to describe quantum state vectors For example, φ is an abstract
More informationFloquet Topological Insulator:
Floquet Topological Insulator: Understanding Floquet topological insulator in semiconductor quantum wells by Lindner et al. Condensed Matter Journal Club Caltech February 12 2014 Motivation Motivation
More informationIntroduction to density functional perturbation theory for lattice dynamics
Introduction to density functional perturbation theory for lattice dynamics SISSA and DEMOCRITOS Trieste (Italy) Outline 1 Lattice dynamic of a solid: phonons Description of a solid Equations of motion
More informationPrinciples of Molecular Spectroscopy
Principles of Molecular Spectroscopy What variables do we need to characterize a molecule? Nuclear and electronic configurations: What is the structure of the molecule? What are the bond lengths? How strong
More informationHandout 10. Applications to Solids
ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem
More informationAddition of Angular Momenta
Addition of Angular Momenta What we have so far considered to be an exact solution for the many electron problem, should really be called exact non-relativistic solution. A relativistic treatment is needed
More informationTopic 9-1: Bloch Theorem and the Central Equation Kittel Pages:
Topic 9-1: Bloch Theorem and the Central Equation Kittel Pages: 167-174 Summary: We begin here by postulating Bloch s theorems which develop the form of the wavefunction in a periodic solid. We then show
More informationRotational motion of a rigid body spinning around a rotational axis ˆn;
Physics 106a, Caltech 15 November, 2018 Lecture 14: Rotations The motion of solid bodies So far, we have been studying the motion of point particles, which are essentially just translational. Bodies with
More informationLecture contents. A few concepts from Quantum Mechanics. Tight-binding model Solid state physics review
Lecture contents A few concepts from Quantum Mechanics Particle in a well Two wells: QM perturbation theory Many wells (atoms) BAND formation Tight-binding model Solid state physics review Approximations
More informationVibrational Motion. Chapter 5. P. J. Grandinetti. Sep. 13, Chem P. J. Grandinetti (Chem. 4300) Vibrational Motion Sep.
Vibrational Motion Chapter 5 P. J. Grandinetti Chem. 4300 Sep. 13, 2017 P. J. Grandinetti (Chem. 4300) Vibrational Motion Sep. 13, 2017 1 / 20 Simple Harmonic Oscillator Simplest model for harmonic oscillator
More informationUnderstanding Phonon Dynamics via 1D Atomic Chains
Understanding Phonon Dynamics via 1D Atomic Chains Timothy S. Fisher Purdue University School of Mechanical Engineering, and Birck Nanotechnology Center tsfisher@purdue.edu Nanotechnology 501 Lecture Series
More informationP3317 HW from Lecture and Recitation 10
P3317 HW from Lecture 18+19 and Recitation 10 Due Nov 6, 2018 Problem 1. Equipartition Note: This is a problem from classical statistical mechanics. We will need the answer for the next few problems, and
More information1 Infinite-Dimensional Vector Spaces
Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationBorn-Oppenheimer Approximation
Born-Oppenheimer Approximation Adiabatic Assumption: Nuclei move so much more slowly than electron that the electrons that the electrons are assumed to be obtained if the nuclear kinetic energy is ignored,
More informationEigenvectors and Hermitian Operators
7 71 Eigenvalues and Eigenvectors Basic Definitions Let L be a linear operator on some given vector space V A scalar λ and a nonzero vector v are referred to, respectively, as an eigenvalue and corresponding
More informationIntroduction to Vibrational Spectroscopy
Introduction to Vibrational Spectroscopy Harmonic oscillators The classical harmonic oscillator The uantum mechanical harmonic oscillator Harmonic approximations in molecular vibrations Vibrational spectroscopy
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationChapter 5 Phonons II Thermal Properties
Chapter 5 Phonons II Thermal Properties Phonon Heat Capacity < n k,p > is the thermal equilibrium occupancy of phonon wavevector K and polarization p, Total energy at k B T, U = Σ Σ < n k,p > ħ k, p Plank
More informationTime dependent perturbation theory 1 D. E. Soper 2 University of Oregon 11 May 2012
Time dependent perturbation theory D. E. Soper University of Oregon May 0 offer here some background for Chapter 5 of J. J. Sakurai, Modern Quantum Mechanics. The problem Let the hamiltonian for a system
More informationPhonons II - Thermal Properties (Kittel Ch. 5)
Phonons II - Thermal Properties (Kittel Ch. 5) Heat Capacity C T 3 Approaches classical limit 3 N k B T Physics 460 F 2006 Lect 10 1 Outline What are thermal properties? Fundamental law for probabilities
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationAn Introduction to Lattice Vibrations
An Introduction to Lattice Vibrations Andreas Wacker 1 Mathematical Physics, Lund University November 3, 2015 1 Introduction Ideally, the atoms in a crystal are positioned in a regular manner following
More informationPhysics with Neutrons I, WS 2015/2016. Lecture 11, MLZ is a cooperation between:
Physics with Neutrons I, WS 2015/2016 Lecture 11, 11.1.2016 MLZ is a cooperation between: Organization Exam (after winter term) Registration: via TUM-Online between 16.11.2015 15.1.2015 Email: sebastian.muehlbauer@frm2.tum.de
More informationLecture 10. Central potential
Lecture 10 Central potential 89 90 LECTURE 10. CENTRAL POTENTIAL 10.1 Introduction We are now ready to study a generic class of three-dimensional physical systems. They are the systems that have a central
More informationLecture 4: Basic elements of band theory
Phys 769 Selected Topics in Condensed Matter Physics Summer 010 Lecture 4: Basic elements of band theory Lecturer: Anthony J. Leggett TA: Bill Coish 1 Introduction Most matter, in particular most insulating
More informationLecture 6. Tight-binding model
Lecture 6 Tight-binding model In Lecture 3 we discussed the Krönig-Penny model and we have seen that, depending on the strength of the potential barrier inside the unit cell, the electrons can behave like
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 8: Lattice Waves in 1D Monatomic Crystals Outline Overview of Lattice Vibrations so far Models for Vibrations in Discrete 1-D Lattice Example of Nearest
More informationNormal modes. where. and. On the other hand, all such systems, if started in just the right way, will move in a simple way.
Chapter 9. Dynamics in 1D 9.4. Coupled motions in 1D 491 only the forces from the outside; the interaction forces cancel because they come in equal and opposite (action and reaction) pairs. So we get:
More informationLecture 14 The Free Electron Gas: Density of States
Lecture 4 The Free Electron Gas: Density of States Today:. Spin.. Fermionic nature of electrons. 3. Understanding the properties of metals: the free electron model and the role of Pauli s exclusion principle.
More informationSolutions for Homework 4
Solutions for Homework 4 October 6, 2006 1 Kittel 3.8 - Young s modulus and Poison ratio As shown in the figure stretching a cubic crystal in the x direction with a stress Xx causes a strain e xx = δl/l
More informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
More informationPH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5
PH575 Spring 2014 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 POP QUIZ Phonon dispersion relation:
More informationPhysics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension
Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension In these notes we examine Bloch s theorem and band structure in problems with periodic potentials, as a part of our survey
More informationPhysics 541: Condensed Matter Physics
Physics 541: Condensed Matter Physics In-class Midterm Exam Wednesday, October 26, 2011 / 14:00 15:20 / CCIS 4-285 Student s Name: Instructions There are 23 questions. You should attempt all of them. Mark
More informationTime-Independent Perturbation Theory
4 Phys46.nb Time-Independent Perturbation Theory.. Overview... General question Assuming that we have a Hamiltonian, H = H + λ H (.) where λ is a very small real number. The eigenstates of the Hamiltonian
More informationChem120a : Exam 3 (Chem Bio) Solutions
Chem10a : Exam 3 (Chem Bio) Solutions November 7, 006 Problem 1 This problem will basically involve us doing two Hückel calculations: one for the linear geometry, and one for the triangular geometry. We
More informationPHYS 3313 Section 001 Lecture # 22
PHYS 3313 Section 001 Lecture # 22 Dr. Barry Spurlock Simple Harmonic Oscillator Barriers and Tunneling Alpha Particle Decay Schrodinger Equation on Hydrogen Atom Solutions for Schrodinger Equation for
More informationIsotropic harmonic oscillator
Isotropic harmonic oscillator 1 Isotropic harmonic oscillator The hamiltonian of the isotropic harmonic oscillator is H = h m + 1 mω r (1) = [ h d m dρ + 1 ] m ω ρ, () ρ=x,y,z a sum of three one-dimensional
More informationLattice Vibrations. Chris J. Pickard. ω (cm -1 ) 200 W L Γ X W K K W
Lattice Vibrations Chris J. Pickard 500 400 300 ω (cm -1 ) 200 100 L K W X 0 W L Γ X W K The Breakdown of the Static Lattice Model The free electron model was refined by introducing a crystalline external
More informationThe Quantum Heisenberg Ferromagnet
The Quantum Heisenberg Ferromagnet Soon after Schrödinger discovered the wave equation of quantum mechanics, Heisenberg and Dirac developed the first successful quantum theory of ferromagnetism W. Heisenberg,
More informationKeble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering
Tomi Johnson Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering Please leave your work in the Clarendon laboratory s J pigeon hole by 5pm on Monday of
More informationPHYS3113, 3d year Statistical Mechanics Tutorial problems. Tutorial 1, Microcanonical, Canonical and Grand Canonical Distributions
1 PHYS3113, 3d year Statistical Mechanics Tutorial problems Tutorial 1, Microcanonical, Canonical and Grand Canonical Distributions Problem 1 The macrostate probability in an ensemble of N spins 1/2 is
More information5.5. Representations. Phys520.nb Definition N is called the dimensions of the representations The trivial presentation
Phys50.nb 37 The rhombohedral and hexagonal lattice systems are not fully compatible with point group symmetries. Knowing the point group doesn t uniquely determine the lattice systems. Sometimes we can
More informationNearly-Free Electrons Model
Nearly-Free Electrons Model Jacob Shapiro March 30, 2013 Abstract The following is my attempt to explain to myself what is going on in the discussion of the Nearly-Free Electron Model which can be found
More informationCHAPTER 6 Quantum Mechanics II
CHAPTER 6 Quantum Mechanics II 6.1 6.2 6.3 6.4 6.5 6.6 6.7 The Schrödinger Wave Equation Expectation Values Infinite Square-Well Potential Finite Square-Well Potential Three-Dimensional Infinite-Potential
More informationPhysics 607 Final Exam
Physics 607 Final Exam Please be well-organized, and show all significant steps clearly in all problems. You are graded on your work, so please do not just write down answers with no explanation! Do all
More information1 Reduced Mass Coordinates
Coulomb Potential Radial Wavefunctions R. M. Suter April 4, 205 Reduced Mass Coordinates In classical mechanics (and quantum) problems involving several particles, it is convenient to separate the motion
More informationModel for vibrational motion of a diatomic molecule. To solve the Schrödinger Eq. for molecules, make the Born- Oppenheimer Approximation:
THE HARMONIC OSCILLATOR Features Example of a problem in which V depends on coordinates Power series solution Energy is quantized because of the boundary conditions Model for vibrational motion of a diatomic
More informationLecture 12: Phonon heat capacity
Lecture 12: Phonon heat capacity Review o Phonon dispersion relations o Quantum nature of waves in solids Phonon heat capacity o Normal mode enumeration o Density of states o Debye model Review By considering
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More information5.62 Physical Chemistry II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 5.62 Physical Chemistry II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.62 Spring 2008 Lecture
More informationEigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator.
PHYS208 spring 2008 Eigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator. 07.02.2008 Adapted from the text Light - Atom Interaction PHYS261 autumn 2007 Go to list of topics
More informationDispersion relation for transverse waves in a linear chain of particles
Dispersion relation for transverse waves in a linear chain of particles V. I. Repchenkov* It is difficult to overestimate the importance that have for the development of science the simplest physical and
More informationChemistry 483 Lecture Topics Fall 2009
Chemistry 483 Lecture Topics Fall 2009 Text PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon A. Background (M&S,Chapter 1) Blackbody Radiation Photoelectric effect DeBroglie Wavelength Atomic
More informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationdynamics computer code, whiccan simulate the motion of atoms at a surface and We coembininewton's equation for the nuclei,
DUI FILE COpy T DTIC ANNUAL LETTER REPORT BY 0. F. ELEC33 20KE AD-A224 069 P.1. FOR N00014-85-K-0442 U Theory of electronic states and formations energies of.defect comp exes, interstitial defects, and
More informationLattice dynamics. Javier Junquera. Philippe Ghosez. Andrei Postnikov
Lattice dynamics Javier Junquera Philippe Ghosez Andrei Postnikov Warm up: a little bit of notation Greek characters ( ) refer to atoms within the unit cell Latin characters ( ) refer to the different
More information(2) A two-dimensional solid has an electron energy band of the form, . [1]
(1) The figure shows a two-dimensional periodic lattice, containing A atoms (white) and B atoms (black). The section of lattice shown makes a 3a 4a rectangle, as shown (measured from A atom to A atom).
More informationChemistry 795T. NC State University. Lecture 4. Vibrational and Rotational Spectroscopy
Chemistry 795T Lecture 4 Vibrational and Rotational Spectroscopy NC State University The Dipole Moment Expansion The permanent dipole moment of a molecule oscillates about an equilibrium value as the molecule
More informationPhysics 828 Problem Set 7 Due Wednesday 02/24/2010
Physics 88 Problem Set 7 Due Wednesday /4/ 7)a)Consider the proton to be a uniformly charged sphere of radius f m Determine the correction to the s ground state energy 4 points) This is a standard problem
More informationPotentially useful reading Sakurai and Napolitano, sections 1.5 (Rotation), Schumacher & Westmoreland chapter 2
Problem Set 2: Interferometers & Random Matrices Graduate Quantum I Physics 6572 James Sethna Due Friday Sept. 5 Last correction at August 28, 2014, 8:32 pm Potentially useful reading Sakurai and Napolitano,
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationAssignment 6. Using the result for the Lagrangian for a double pendulum in Problem 1.22, we get
Assignment 6 Goldstein 6.4 Obtain the normal modes of vibration for the double pendulum shown in Figure.4, assuming equal lengths, but not equal masses. Show that when the lower mass is small compared
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More informationChapter 2 The Group U(1) and its Representations
Chapter 2 The Group U(1) and its Representations The simplest example of a Lie group is the group of rotations of the plane, with elements parametrized by a single number, the angle of rotation θ. It is
More informationChemistry 881 Lecture Topics Fall 2001
Chemistry 881 Lecture Topics Fall 2001 Texts PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon MATHEMATICS for PHYSICAL CHEMISTRY, Mortimer i. Mathematics Review (M, Chapters 1,2,3 & 4; M&S,
More information