Beat Frequency Slides Beat Frequency
Energy Exchange in Coupled Systems
Lesson 1: Coupled Oscillations
The Curious Case of Energy Transfer
Imagine two identical pendulums connected by a weak spring.
If you displace only one pendulum and release it:
It gradually loses amplitude until it stops.
Simultaneously, the second pendulum begins to swing.
The energy oscillates back and forth between them.
Why does this happen?
Periodic Energy Exchange
The Model System
Consider two masses \(m\) connected by three springs with constants \(k, \kappa, k\).
\(k\)
m
\(\kappa\)
m
\(k\)
Coordinates
\(x_1(t)\) and \(x_2(t)\) represent displacements from equilibrium.
Coupling
The central spring \(\kappa\) couples the motion of the two masses.
Equations of Motion
Newton's Second Law for each mass:
\[ m\ddot{x}_1 = -kx_1 - \kappa(x_1 - x_2) \]
\[ m\ddot{x}_2 = -kx_2 - \kappa(x_2 - x_1) \]
Note: These are coupled linear differential equations.
The Normal Modes
Symmetric (\(x_1 = x_2\))
Masses move together; the coupling spring length is constant.
\(\omega_s = \sqrt{\frac{k}{m}}\)
Anti-symmetric (\(x_1 = -x_2\))
Masses move in opposition; the coupling spring is maximally stressed.
\(\omega_a = \sqrt{\frac{k + 2\kappa}{m}}\)
Superposition and Beats
Any motion is a linear combination of the two normal modes:
\(x_1(t) = A\cos(\omega_s t) + B\cos(\omega_a t)\)
When the coupling is weak (\(\kappa \ll k\)), the frequencies are close:
Beat Frequency: \(\Delta\omega = \omega_a - \omega_s\)
This creates an envelope that modulates the amplitude of oscillation.
Amplitude Modulation (Beats)
Energy Flow
Initial
Mass 1 has all the energy. Mass 2 is at rest.
Exchange
Coupling spring does work, transferring power between masses.
The time for a complete energy transfer from \(M_1\) to \(M_2\) and back is the Beat Period.
Key Takeaways
01
Coupled systems exhibit collective behavior that can be decomposed into independent normal modes.
02
Normal modes occur when all parts of the system oscillate with the same frequency and phase.
03
Beats arise from the superposition of modes with slightly different frequencies.
Coupled Pendulums Teacher Guide Coupled Pendulums Facilitation Guide
Lesson 1: Beats and Energy Exchange
Teacher Resource
Lesson Objective
To provide a phenomenological and mathematical foundation for coupled oscillations. Students will move from observing energy exchange (beats) to identifying the underlying normal modes that remain invariant in their frequency of oscillation.
Essential Hook
"When two pendulums are connected by a spring, why does one stop completely while the other swings, only to reverse roles moments later? How can we predict the time it takes for this energy baton to be passed?"
Quick Specs
90 Minutes
Graduate / Upper-Undergrad
ODE, Linear Superposition
Facilitation Timeline
15 min
Observation & Hypothesis
Show a demonstration (physical or simulation) of two coupled pendulums. Ask students to describe the motion of one pendulum vs. the system as a whole. Identify the "Envelope" and the "Carrier" frequencies.
30 min
Formal Modeling
Derive the equations of motion on the board. Focus on the coupling term \(\kappa(x_1 - x_2)\). Emphasize that the force depends on the relative displacement.
30 min
Normal Mode Identification
Guide the class to "guess" the simple solutions: Symmetric and Anti-symmetric. Prove that these modes decouple the equations into two independent harmonic oscillators.
15 min
Synthesis & Beats
Show how the general solution is a sum of modes. Relate the difference in mode frequencies to the beat frequency observed in the initial hook.
Technical Guidance
Key Derivation: Decoupling
Standard coordinates \(x_1, x_2\) are frustrating. Introduce:
\(q_s = x_1 + x_2\) (Symmetric)
Show that adding and subtracting the equations of motion leads directly to:
\(\ddot{q}_s + \omega_s^2 q_s = 0\)
The Weak Coupling Limit
The beat frequency \(\Delta\omega = \omega_a - \omega_s\) is most apparent when \(\kappa \ll k\). Use the binomial expansion:
\(\omega_a \approx \omega_s (1 + \frac{\kappa}{k})\)
Common Misconceptions
Mode Identity: Students often think "modes" are things the system does sequentially. Clarify that modes are basis states; the system is usually in a superposition of all modes.
Beats vs. Modes: Students may confuse the beat frequency with a fundamental frequency of the system. Stress that \(\omega_{beat}\) is a derived frequency of the envelope, not an eigenfrequency.
Beats and Energy Worksheet Beats and Energy Dynamics
Problem Set 1.1: Coupled Linear Oscillators
Name:
Date:
"The problem of coupled oscillations is the first bridge from the mechanics of a single particle to the collective physics of many-body systems."
1
The Governing Equations
Consider two identical masses \(m\) connected to walls by springs of constant \(k\), and to each other by a central spring of constant \(\kappa\). Let \(x_1\) and \(x_2\) be the displacements from equilibrium.
A) Write the Lagrangian \(L = T - V\) for this system in terms of \(x_1, x_2, \dot{x}_1, \dot{x}_2\).
B) Derive the coupled equations of motion using the Euler-Lagrange equations.
2
Normal Modes and Energy Exchange
Assume the symmetric and anti-symmetric solutions: \(x_1 = x_2\) and \(x_1 = -x_2\).
Determine the eigenfrequencies \(\omega_s\) and \(\omega_a\) for these modes.
C) For the initial conditions \(x_1(0) = a, \dot{x}_1(0) = 0, x_2(0) = 0, \dot{x}_2(0) = 0\), show that the motion of the first mass can be written as:
\(x_1(t) = a \cos\left(\frac{\omega_a - \omega_s}{2}t\right) \cos\left(\frac{\omega_a + \omega_s}{2}t\right)\)
Identify the carrier frequency and the modulation frequency. How does the period of energy transfer relate to these frequencies?
D) Qualitative Sketch: On the axes below, sketch the displacement \(x_1(t)\) for the case where \(\kappa \ll k\). Label the "Beat Period" clearly.
Displacement (x)
Time (t) →
3
Physical Reflection
Explain why the energy is transferred completely between identical oscillators, but only partially transferred if the masses are different.
In the anti-symmetric mode, the central spring is compressed more than in the symmetric mode. Relate this physical observation to the relative magnitudes of \(\omega_s\) and \(\omega_a\).
Challenge Extension
If we add a damping force \(b\dot{x}\) to each mass, do the normal modes change? Does the energy transfer still occur indefinitely? Provide a qualitative argument for how the beat envelope would decay.
Matrix Mechanics Slides Matrix Mechanics
Formalizing N-Coupled Oscillators
Lesson 2: Lagrangian Eigenvalue Problems
Beyond Two Masses
In Lesson 1, we "guessed" the symmetric and anti-symmetric modes. This approach fails for:
Systems with unequal masses.
Asymmetric spring configurations.
Systems with \(N\) components.
We need a systematic tool: Linear Algebra.
The Quadratic Lagrangian
For small oscillations, the kinetic and potential energies are quadratic forms:
Kinetic Energy
\(T = \frac{1}{2} \sum_{ij} M_{ij} \dot{x}_i \dot{x}_j\)
Represented by the Mass Matrix (M).
Potential Energy
\(V = \frac{1}{2} \sum_{ij} K_{ij} x_i x_j\)
Represented by the Potential Matrix (K).
The matrices M and K are symmetric and positive definite.
The Secular Equation
Assuming a harmonic solution \(\vec{x}(t) = \vec{a} e^{i\omega t}\), the Euler-Lagrange equations become:
The Core Problem
\((K - \omega^2 M) \vec{a} = 0\)
Non-trivial solutions exist if and only if the determinant vanishes.
\(\det(K - \omega^2 M) = 0\)
Generalized Eigenvalue Recipe
1
Define generalized coordinates \(q_i\).
2
Write \(T\) and \(V\) to identify matrices \(M\) and \(K\).
3
Solve the Secular Equation for eigenfrequencies \(\omega_i^2\).
4
Substitute each \(\omega_i^2\) back to find eigenvectors \(\vec{a}_i\).
Example: Symmetric Case
Lesson 1 Review
Mass Matrix
\(M = \begin{pmatrix} m & 0 \\ 0 & m \end{pmatrix}\)
Potential Matrix
\(K = \begin{pmatrix} k+\kappa & -\kappa \\ -\kappa & k+\kappa \end{pmatrix}\)
Solving \(\det(K - \omega^2 M) = 0\) yields:
\(\omega_s^2 = k/m\)
\(\omega_a^2 = (k + 2\kappa)/m\)
Exactly the same result, but derived algorithmically!
Orthogonality & Normalization
The eigenvectors are orthogonal with respect to the mass matrix:
\(\vec{a}_i^T M \vec{a}_j = \delta_{ij}\)
This "Mass-Weighted Orthogonality" is the key to decoupling complex systems.
In Lesson 3, we will use this property to define Normal Coordinates.
Secular Equation Practice The Secular Equation
Problem Set 2.1: Matrix Formulation
Name:
1
Asymmetric Coupled Pendulums
Consider the same configuration as Lesson 1, but with different masses: \(m_1 = m\) and \(m_2 = 2m\). The spring constants remain \(k\) (outer) and \(\kappa\) (coupling).
A) Construct the Mass Matrix \(M\) and the Potential Matrix \(K\).
Matrix M
Matrix K
B) Set up the secular equation \(\det(K - \omega^2 M) = 0\). Expand the determinant to find the characteristic polynomial for \(\omega^2\).
2
Triatomic Linear Molecule
A simple model of a linear triatomic molecule (e.g., \(CO_2\)) consists of three masses in a line connected by two identical springs \(k\). Let the central mass be \(M\) and the outer masses be \(m\).
m
M
m
Determine the potential matrix \(K\) using the coordinates \(x_1, x_2, x_3\). Why is there no term for \(x_1\) and \(x_3\) coupling directly?
By inspection or derivation, identify the "Zero Frequency" mode (\(\omega = 0\)). What does this mode represent physically?
Solve for the two non-zero eigenfrequencies of the system.
Conceptual Challenge
"A system with \(N\) degrees of freedom will have \(N\) eigenfrequencies, some of which may be degenerate (identical)."
Suppose we add another spring to the triatomic molecule, connecting the two outer masses \(m\). How would this affect the \(M\) and \(K\) matrices? Would the zero-frequency mode still exist? Explain.
Matrix Eigenvalue Key Matrix Eigenvalue Key
Teacher Reference: Lesson 2
Answer Key
Problem 1: Asymmetric Masses (\(m, 2m\))
A) Matrices
\(M = \begin{pmatrix} m & 0 \\ 0 & 2m \end{pmatrix}\)
\(K = \begin{pmatrix} k+\kappa & -\kappa \\ -\kappa & k+\kappa \end{pmatrix}\)
B) Secular Equation
The determinant \(\det(K - \omega^2 M) = 0\) leads to:
\((k+\kappa - m\omega^2)(k+\kappa - 2m\omega^2) - \kappa^2 = 0\)
Expanding and solving for \(\omega^2\) yields a quadratic in \(\omega^2\):
\(2m^2\omega^4 - 3m(k+\kappa)\omega^2 + (k^2 + 2k\kappa) = 0\)
Encourage students to use the quadratic formula. The resulting frequencies will no longer correspond to purely symmetric/anti-symmetric motion. The masses will oscillate with different relative amplitudes.
Problem 2: Linear Triatomic Molecule
The Matrices (Coordinates \(x_1, x_2, x_3\))
\(M = \begin{pmatrix} m & 0 & 0 \\ 0 & M & 0 \\ 0 & 0 & m \end{pmatrix}\)
\(K = \begin{pmatrix} k & -k & 0 \\ -k & 2k & -k \\ 0 & -k & k \end{pmatrix}\)
Teacher Note: Coupling Inspection
Students should notice \(K_{13} = K_{31} = 0\). This is because the outer masses are not directly connected by a spring; they "feel" each other only through the central mass \(M\).
Eigenfrequencies
1
\(\omega_1 = 0\)
Pure Translation. The eigenvector is \((1, 1, 1)\). The molecule moves as a rigid body.
2
\(\omega_2^2 = k/m\)
Symmetric Stretch. The central mass remains stationary. Eigenvector \((1, 0, -1)\).
3
\(\omega_3^2 = \frac{k}{m}(1 + 2m/M)\)
Asymmetric Stretch. The outer masses move one way while the inner moves the other. Eigenvector \((1, -2m/M, 1)\).
Conceptual Challenge Response
Adding a spring between outer masses changes \(K_{11}\) and \(K_{33}\) (increases diagonal terms) and makes \(K_{13} = K_{31} = -\kappa_{new}\). However, the row/column sums of \(K\) will still be zero if the system is not anchored to a wall. This ensures the \(\omega=0\) translational mode is preserved.
Normal Modes Slides Normal Coordinates
The Art of Decoupling
Lesson 3: Eigenvector Basis & Mode Shapes
What is a Mode Shape?
An eigenvector \(\vec{a}_i\) associated with frequency \(\omega_i\) defines the relative amplitude and phase of all particles in that mode.
Property:
In a pure normal mode, all particles pass through equilibrium at the same time and reach their extrema simultaneously.
Example: Symmetric Mode
\(\vec{a}_s = \begin{pmatrix} 1 \\ 1 \end{pmatrix}\)
Example: Anti-symmetric Mode
\(\vec{a}_a = \begin{pmatrix} 1 \\ -1 \end{pmatrix}\)
The Coordinate Transformation
We define new coordinates \(\eta_i\) (normal coordinates) such that the physical displacement is a weighted sum of mode shapes:
\(\vec{x}(t) = \sum_{i=1}^N \eta_i(t) \vec{a}_i\)
\(\vec{x} = A \vec{\eta}\)
Where \(A\) is the modal matrix containing eigenvectors as columns.
Diagonalization
Substituting \(\vec{x} = A\vec{\eta}\) into the kinetic and potential energy forms:
Kinetic
\(T = \frac{1}{2} \dot{\vec{\eta}}^T (A^T M A) \dot{\vec{\eta}}\)
Potential
\(V = \frac{1}{2} \vec{\eta}^T (A^T K A) \vec{\eta}\)
The Magic Result
By orthogonality, both \(A^T M A\) and \(A^T K A\) become diagonal matrices!
The equations of motion for \(\eta_i\) are completely independent simple harmonic oscillators.
Visualization: CO2 Modes
Symmetric Stretch
Central mass stays still.
Asymmetric Stretch
Center moves opposite to ends.
Translation
System moves as one.
Key Insight
"Any arbitrary vibration of a complex system can be uniquely represented as a linear combination of its normal modes, each oscillating independently at its own natural frequency."
Molecule Mode Discovery The Modal Transformation
Activity 3.1: Discovering Independent Modes
Name:
1
Constructing the Modal Matrix
Recall the linear triatomic molecule from Lesson 2. Let the masses be \(m, M, m\) and the spring constants be \(k\). The eigenvectors (mode shapes) were:
\(\vec{a}_1 = \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}\)
\(\vec{a}_2 = \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix}\)
\(\vec{a}_3 = \begin{pmatrix} 1 \\ -2m/M \\ 1 \end{pmatrix}\)
A) Construct the modal matrix \(A = [\vec{a}_1 | \vec{a}_2 | \vec{a}_3]\).
A =
B) Verify the mass-weighting orthogonality: Calculate the product \(A^T M A\) and show that it is diagonal. (Note: Do not worry about normalizing to unity yet, just show it is diagonal).
2
Independent Harmonic Oscillators
After the transformation \(\vec{x} = A\vec{\eta}\), we get three independent equations for the normal coordinates \(\eta_1, \eta_2, \eta_3\).
Write the resulting general solution for \(\eta_i(t)\) in terms of their respective eigenfrequencies.
C) Physical Synthesis: If we excite the system with initial displacements \(\vec{x}(0) = (1, 0, 1)\) and zero initial velocity, which normal modes are present in the resulting motion? Solve for the coefficients \(\eta_i(0)\).
Discovery Question
Why is it practically impossible to excite only a single normal mode in a laboratory setting? What physical imperfections (non-linearities, air resistance, mass variations) would cause other modes to "leak" in over time?
Coordinate Transformation Guide Coordinate Transformation Guide
Teacher Resource: Lesson 3
Reference
The Logic of Transformation
The fundamental goal of this lesson is to show that a complex, coupled physical system is mathematically equivalent to a set of independent simple harmonic oscillators when viewed through the "lens" of normal coordinates.
Properties of the Modal Matrix \(A\)
The columns of \(A\) are the eigenvectors \(\vec{a}_i\).
\(A^T M A\) results in a diagonal matrix \(\mathcal{M}\) (the modal mass).
\(A^T K A\) results in a diagonal matrix \(\mathcal{K}\) (the modal stiffness).
The natural frequencies are \(\omega_i = \sqrt{\mathcal{K}_{ii}/\mathcal{M}_{ii}}\).
Pedagogical Tip
"Emphasize that while displacements \(x_i\) are the quantities we measure in the lab, the normal coordinates \(\eta_i\) are the quantities that describe the underlying physics of the system."
Activity 3.1 Solution Key
1A) The Modal Matrix
\(A = \begin{pmatrix} 1 & 1 & 1 \\ 1 & 0 & -2m/M \\ 1 & -1 & 1 \end{pmatrix}\)
1B) Diagonalization
Multiplying \(A^T M A\) should yield:
\(\mathcal{M} = \text{diag}(2m+M, 2m, M(1 + 2m/M))\)
Note: If student normalization is different, the values will vary, but the off-diagonals must be zero.
2C) Physical Synthesis
The vector \(\vec{x}(0) = (1, 0, 1)\) is a combination of translation and the asymmetric stretch. To find \(\eta_i(0)\), solve \(A\vec{\eta}(0) = \vec{x}(0)\). The symmetric mode (\(\vec{a}_2\)) will have zero amplitude because the end masses are displaced symmetrically.
Technical Deep Dive
Orthonormalization
Explain that if we normalize \(\vec{a}_i\) such that \(\vec{a}_i^T M \vec{a}_j = \delta_{ij}\), then:
\(A^T M A = I\) (Identity)
This is the "canonical" form of the transformation and is extremely useful for numerical simulations.
Research Connection
In quantum mechanics, this coordinate transformation is essentially the same as changing to a basis where the Hamiltonian is diagonal. In solid-state physics, normal modes of a crystal lattice are called phonons. The decoupling we perform here classically is the exact same math used to treat phonons as independent "particles" (quasiparticles).
Infinite Chain Slides The Infinite Chain
From Discrete to Continuum
Lesson 4: Dispersion Relations
The 1D Lattice
Consider an infinite chain of identical masses \(m\) connected by springs \(k\) with lattice spacing \(a\).
n-1
\(a\)
n
\(a\)
n+1
Equation of Motion for the \(n\)-th particle:
\(m\ddot{u}_n = k(u_{n+1} + u_{n-1} - 2u_n)\)
Traveling Wave Solution
Assume a solution in the form of a traveling wave:
\(u_n(t) = A e^{i(kna - \omega t)}\)
\(k\): Wavenumber
\(\omega\): Angular frequency
\(na\): Position of the \(n\)-th atom
The Substitution
Substituting this into the EOM and using the identity \((e^{ika} + e^{-ika} - 2) = -4\sin^2(ka/2)\):
\(-m\omega^2 = -4k \sin^2(ka/2)\)
The Dispersion Relation
Characteristic Frequencies
\(\omega(k) = 2\sqrt{\frac{k}{m}} \left| \sin\left(\frac{ka}{2}\right) \right|\)
The frequency \(\omega\) is a periodic function of the wavenumber \(k\).
1. Small \(k\) Limit (Continuum)
\(\sin(ka/2) \approx ka/2 \implies \omega \approx (a\sqrt{k/m})k\)
The relationship is linear, representing sound waves.
2. Short Wavelength Limit
At \(k = \pm \pi/a\), \(\omega\) reaches a maximum: \(\omega_{max} = 2\sqrt{k/m}\).
The system acts as a low-pass filter.
The First Brillouin Zone
-\(\pi\)/a \(\pi\)/a k \(\omega\)
Because atoms are at discrete points \(na\), any wave with \(k' = k + 2\pi m/a\) is physically indistinguishable from \(k\).
The First Brillouin Zone:
\(-\frac{\pi}{a} < k \leq \frac{\pi}{a}\)
This range contains all unique physical modes of the lattice.
The Wave Equation
As the lattice spacing \(a \to 0\) and number of atoms \(N \to \infty\):
The Classical Wave Equation
\(\frac{\partial^2 u}{\partial t^2} - v^2 \frac{\partial^2 u}{\partial x^2} = 0\)
The discrete coupled oscillators merge into a continuous medium.
Dispersion Derivation Notes Dispersion Derivation Notes
Teacher Resource: Lesson 4
Lecture Notes
1. The Difference Equation
The infinite chain is governed by a set of infinite, coupled linear differential equations. For mass \(n\):
\(m\ddot{u}_n = k(u_{n+1} - u_n) - k(u_n - u_{n-1}) = k(u_{n+1} + u_{n-1} - 2u_n)\)
2. The Traveling Wave Ansätz
Assume a periodic solution (Bloch wave) of the form \(u_n = A e^{i(kna - \omega t)}\). Substituting the second derivative:
\(-m\omega^2 e^{i(kna - \omega t)} = k [ e^{ik(n+1)a} + e^{ik(n-1)a} - 2e^{ikna} ] e^{-i\omega t}\)
Dividing by \(e^{i(kna - \omega t)}\) leaves:
\(-m\omega^2 = k [ e^{ika} + e^{-ika} - 2 ]\)
3. Trig Identity & Final Relation
Using \(e^{ika} + e^{-ika} = 2\cos(ka)\):
\(-m\omega^2 = 2k (\cos(ka) - 1)\)
Since \(1 - \cos(ka) = 2\sin^2(ka/2)\):
\(\omega(k) = 2\sqrt{\frac{k}{m}} \left| \sin\left(\frac{ka}{2}\right) \right|\)
Key Discussion Points
Phase vs. Group Velocity
Unlike waves in a continuous string, waves in a lattice are dispersive.
Phase Velocity: \(v_p = \omega/k\)
Group Velocity: \(v_g = d\omega/dk = a\sqrt{k/m}\cos(ka/2)\)
Point out that at the BZ boundary (\(k = \pi/a\)), \(v_g = 0\). This represents a standing wave.
Brillouin Zone Intuition
Why stop at \(\pi/a\)? Show students that if we increase \(k\) by \(2\pi/a\), the displacements \(u_n\) at the atomic sites do not change.
A lattice "samples" the wave at discrete intervals \(a\). Any wavelength shorter than \(2a\) (higher frequency than \(\pi/a\)) is aliased down to a longer wavelength in the First BZ.
The Cutoff Frequency
The lattice acts as a mechanical low-pass filter. It cannot transmit vibrations with frequencies higher than \(\omega_{max} = 2\sqrt{k/m}\). In the student worksheet, have them think about what happens to energy above this limit (evanescent waves).
Chain Limit Worksheet The Monatomic Chain
Problem Set 4.1: Dispersion in Lattice Solids
Name:
1
Deriving the Relation
Start from the equation of motion \(m\ddot{u}_n = k(u_{n+1} + u_{n-1} - 2u_n)\). Assume the Bloch solution \(u_n = A e^{i(kna - \omega t)}\).
A) Show the explicit steps that lead from the substitution to the trigonometric form: \(\omega^2 = \frac{4k}{m} \sin^2(ka/2)\).
2
Velocities and Dispersion
B) Derive the expressions for the phase velocity \(v_p = \omega/k\) and the group velocity \(v_g = d\omega/dk\).
C) Calculate the limit of \(v_p\) and \(v_g\) as \(k \to 0\). Are they equal? Explain the physical significance of this limit for acoustic sound waves.
D) What is the group velocity at the Brillouin Zone boundary \(k = \pi/a\)? Sketch the spatial configuration of the atoms for this specific mode.
3
Transition to Continuous Media
Consider the limit where the atomic spacing \(a\) becomes infinitesimally small (\(a \to 0\)). Let \(\rho = m/a\) be the linear mass density and \(Y = ka\) be the Young's modulus property of the spring chain.
Show that in this limit, the discrete difference equation reduces to the second-order partial differential equation for a continuous string.
Thinking Problem
Real crystals are three-dimensional. How would the dispersion relation change if the atoms could vibrate in three dimensions instead of just one? Do you expect the Brillouin Zone to be a simple interval or something more complex?
Phonon Physics Slides Phonon Physics
Acoustic & Optical Branches
Lesson 5: Diatomic Chains
Two Atoms per Unit Cell
Consider a chain with two masses \(M\) and \(m\) and distance \(a\) between identical masses.
M
m
M
m
M
A System of Two Equations:
\(M\ddot{u}_{2n} = k(u_{2n+1} + u_{2n-1} - 2u_{2n})\)
The Dispersion Relation
Solving the secular equation for the diatomic system yields two solutions for \(\omega^2\) at each \(k\):
\(\omega^2 = k\left(\frac{1}{m} + \frac{1}{M}\right) \pm k\sqrt{\left(\frac{1}{m} + \frac{1}{M}\right)^2 - \frac{4\sin^2(ka/2)}{mM}}\)
(-) Branch: Acoustic
Behaves like sound at small k.
(+) Branch: Optical
Interact with light (dipole moments).
The Forbidden Gap
The mass difference opens a Band Gap at the BZ boundary.
Optical: Atoms oscillate against each other.
Acoustic: Atoms move in the same direction.
No waves can propagate with frequencies in the gap region.
GAP
Why "Optical"?
In ionic crystals (like NaCl), the two atoms have opposite charges.
When they oscillate in opposition, they create a fluctuating dipole moment.
This dipole can absorb or emit electromagnetic radiation in the infrared spectrum.
-
The Phonon Picture
"By mastering the classical mechanics of coupled oscillators, we have built the foundation for Solid State Physics. In the quantum world, these independent modes are the particles that carry heat and sound."
Discrete
Atomic Strings
Continuum
Elastic Waves
Quantum
Phonons
Diatomic Lattice Guide Diatomic Lattice Facilitation
Teacher Resource: Lesson 5
Closing Lesson
The Complexity Factor
The transition from Lesson 4 to Lesson 5 is the introduction of a basis. Instead of one atom per unit cell, we now have two. This doubles the size of our secular equation and gives us two branches.
The Acoustic Branch
At \(k=0\), the frequency \(\omega=0\). In this limit, the unit cell moves as a whole. This is the sound limit. Both atoms in the cell move in the same direction.
\(\omega^2 \approx \frac{1}{2} \frac{k}{M+m} a^2 k^2\)
The Optical Branch
At \(k=0\), \(\omega\) is non-zero. The two atoms in the unit cell move in opposite directions, keeping the center of mass stationary. In an ionic crystal, this creates a dipole.
\(\omega^2 = 2k (\frac{1}{m} + \frac{1}{M})\)
The Band Gap
Highlight that the gap size depends on the mass ratio. If \(M = m\), the gap closes and we recover the monatomic chain (with a halved Brillouin zone). As \(M\) and \(m\) become very different, the gap grows.
Pedagogical Synthesis
Misconception: Gap Propagation
Students often think the gap means no waves exist. Clarify that traveling waves do not exist. If you drive the system at a gap frequency, the amplitude will decay exponentially with distance (evanescent waves).
Physical Analogy
Compare the optical mode to two children on a trampoline pushing against each other, vs the acoustic mode being them jumping together.
"The optical mode is the internal vibration of the basis, while the acoustic mode is the vibration of the lattice as a whole."
Final Question for Class
"If heat is carried by phonons, and optical phonons have a low group velocity (very flat \(\omega(k)\)), which branch is more responsible for thermal conductivity in a solid?"
Answer: Acoustic phonons carry most of the heat because their group velocity \(d\omega/dk\) is much higher.
Acoustic vs Optical Challenge The Diatomic Challenge
Problem Set 5.1: Acoustic and Optical Phonons
Name:
1
Dispersion Mapping
For a diatomic chain with masses \(M\) and \(m\) (where \(M > m\)), the dispersion relation has two solutions for every \(k\).
A) Sketch the dispersion relation \(\omega(k)\) over the first Brillouin Zone \((-\pi/a, \pi/a)\). Label the Acoustic branch, the Optical branch, and the Forbidden Gap.
\(-\pi/a\)
\(\pi/a\)
\(\omega\)
2
Boundary Dynamics
At the Brillouin Zone boundary \(k = \pm\pi/a\), the frequencies are given by:
\(\omega_{lower} = \sqrt{2k/M}\)
\(\omega_{upper} = \sqrt{2k/m}\)
B) Describe the motion of the atoms in each case. Which masses are moving and which are stationary at these boundary frequencies? Why?
C) As \(M \to m\), what happens to the size of the band gap? Relate your answer to the physical symmetry of the system.
3
The Big Picture
Synthesis Reflection
We have traced the journey from two coupled pendulums to the complex vibrations of a lattice.
How does the concept of a "Normal Mode" allow us to simplify a system of \(10^{23}\) coupled atoms into something manageable? Why is it mathematically significant that these modes are independent?
Finally, speculate on why the Optical branch is named "Optical". How would the presence of an electric field (light) interact differently with the Acoustic vs. Optical modes?
