Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

Posts Tagged ‘anharmonic oscillator’

Final update of notes for PHY487 (condensed matter physics)

Posted by peeterjoot on January 20, 2014

Here is what will likely be the final update of my class notes from Winter 2013, University of Toronto Condensed Matter Physics course (PHY487H1F), taught by Prof. Stephen Julian.

Official course description: “Introduction to the concepts used in the modern treatment of solids. The student is assumed to be familiar with elementary quantum mechanics. Topics include: bonding in solids, crystal structures, lattice vibrations, free electron model of metals, band structure, thermal properties, magnetism and superconductivity (time permitting)”

This document contains:

• Plain old lecture notes. These mirror what was covered in class, possibly augmented with additional details.
• Personal notes exploring details that were not clear to me from the lectures, or from the texts associated with the lecture material.
• Assigned problems. Like anything else take these as is.
• Some worked problems attempted as course prep, for fun, or for test preparation, or post test reflection.
• Links to Mathematica workbooks associated with this course.
My thanks go to Professor Julian for teaching this course.

NOTE: This v.5 update of these notes is still really big (~18M).  Some of my mathematica generated 3D images result in very large pdfs.

Changelog for this update (relative to the first, and second, and third, and the last pre-exam Changelogs).

January 19, 2014 Quadratic Deybe

January 19, 2014 One atom basis phonons in 2D

January 07, 2014 Two body harmonic oscillator in 3D
Figure out a general solution for two interacting harmonic oscillators, then use the result to calculate the matrix required for a 2D two atom diamond lattice with horizontal, vertical and diagonal nearest neighbour coupling.

December 04, 2013 Lecture 24: Superconductivity (cont.)

December 04, 2013 Problem Set 10: Drude conductivity and doped semiconductors.

Advertisements

Posted in Math and Physics Learning. | Tagged: , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , | Leave a Comment »

An updated compilation of notes, for ‘PHY452H1S Basic Statistical Mechanics’, Taught by Prof. Arun Paramekanti

Posted by peeterjoot on March 27, 2013

Here’s my second update of my notes compilation for this course, including all of the following:

March 27, 2013 Fermi gas

March 26, 2013 Fermi gas thermodynamics

March 26, 2013 Fermi gas thermodynamics

March 23, 2013 Relativisitic generalization of statistical mechanics

March 21, 2013 Kittel Zipper problem

March 18, 2013 Pathria chapter 4 diatomic molecule problem

March 17, 2013 Gibbs sum for a two level system

March 16, 2013 open system variance of N

March 16, 2013 probability forms of entropy

March 14, 2013 Grand Canonical/Fermion-Bosons

March 13, 2013 Quantum anharmonic oscillator

March 12, 2013 Grand canonical ensemble

March 11, 2013 Heat capacity of perturbed harmonic oscillator

March 10, 2013 Langevin small approximation

March 10, 2013 Addition of two one half spins

March 10, 2013 Midterm II reflection

March 07, 2013 Thermodynamic identities

March 06, 2013 Temperature

March 05, 2013 Interacting spin

plus everything detailed in the description of my first update and before.

Posted in Math and Physics Learning. | Tagged: , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , | 1 Comment »

Heat capacity of perturbed harmonic oscillator

Posted by peeterjoot on March 11, 2013

[Click here for a PDF of this post with nicer formatting (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]

This problem was suggested as prep for the second midterm, but I spent too much time on my problem set. That’s pretty unfortunate since this showed exactly the approach that was expected for the second midterm problem. Not hard, just not obvious in the heat of the moment how to do that Taylor expansion.

Question: Anharmonic oscillator ([1] pr 3.29)

The potential energy of a one-dimensional, anharmonic oscillator may be written as

\begin{aligned}V(q) = c q^2 - g q^3 - f q^4,\end{aligned} \hspace{\stretch{1}}(1.1)

where c, g, and f are positive constant; quite generally, g and f may be assumed to be very small in value.

Show that the leading contribution of anharmonic terms to the heat cpacity of the oscillator, assumed classical, is given by

\begin{aligned}\frac{3}{2} k_{\mathrm{B}}^2 \left( {\frac{f}{c^2} + \frac{5}{4} \frac{g^2}{c^3} } \right) T,\end{aligned} \hspace{\stretch{1}}(1.2)

To the same order, show that the mean value of the position coordinate q is given by

\begin{aligned}\frac{3}{4} \frac{g k_{\mathrm{B}} T}{c^2}.\end{aligned} \hspace{\stretch{1}}(1.3)

Answer

Our partition function is

\begin{aligned}Z = \int dp dq e^{-\beta p^2/2m} e^{ -\beta \left( { c q^2 - g q^3 - f q^4} \right) }= \sqrt{\frac{2 \pi m}{\beta} }\int dq e^{ -\beta \left( { c q^2 - g q^3 - f q^4} \right) }\end{aligned} \hspace{\stretch{1}}(1.4)

How to expand this wasn’t immediately clear to me (as it wasn’t on the midterm either). We can’t Taylor expand in q, because there’s no single position q that is of interest to expand around (we are integrating over all q). What we can do though is Taylor expand about the values f and g, which are assumed to be small. Here’s the two variable Taylor expansion of this perturbed harmonic oscillator exponential. With

\begin{aligned}A(f, g) = e^{ -\beta \left( { c q^2 - g q^3 - f q^4} \right) }\end{aligned} \hspace{\stretch{1}}(1.5)

The expansion to second order is

\begin{aligned}A(f, g) &= A(0, 0) + f {\left. \frac{\partial A}{\partial f} \right\vert}_{f = 0} + g {\left. \frac{\partial A}{\partial g} \right\vert}_{g = 0} + \frac{1}{2} f^2 {\left. \frac{\partial^2 A}{\partial f^2} \right\vert}_{f = 0} + \frac{1}{2} g^2 {\left. \frac{\partial^2 A}{\partial g^2} \right\vert}_{g = 0} + f g {\left. \frac{\partial^2 A}{\partial g \partial f} \right\vert}_{f, g = 0} + \cdots \\ &= e^{ -\beta c q^2 } \left( 1 + g \beta q^3 + f \beta q^4+ \frac{1}{2} g^2 \left( \beta q^3 \right)^2 + \frac{1}{2} f^2 \left( \beta q^4 \right)^2 + f g \left( \beta q^3 \right) \left( \beta q^4 \right) + \cdots \right) \\ &= e^{ -\beta c q^2 } \left( 1 + g \beta q^3 + f \beta q^4+ \frac{1}{2} g^2 \beta^2 q^6+ f g \beta^2 q^7+ \frac{1}{2} f^2 \beta^2 q^8 + \cdots \right) \end{aligned} \hspace{\stretch{1}}(1.6)

This can now be integrated by parts, where any odd powers are killed. For even powers we have

\begin{aligned}\int q^{2 N} e^{-a q^2} dq &= \int q^{2 N - 1} d \frac{ e^{-a q^2} }{-2 a} \\ &= \frac{2 N - 1}{2a} \int q^{2 (N - 1)} e^{-a q^2} dq \\ &= \frac{(2 N - 1)!!}{(2a)^N} \sqrt{\frac{\pi}{a}}.\end{aligned} \hspace{\stretch{1}}(1.7)

This gives us

\begin{aligned}Z &= \sqrt{ \frac{\pi}{\beta c} }\sqrt{\frac{2 \pi m}{\beta} }\left( {1 + f \beta \frac{3!!}{(2 \beta c)^2}+ \frac{1}{{2}} g^2 \beta^2 \frac{5!!}{(2 \beta c)^3}+ \frac{1}{{2}} f^2 \beta^2 \frac{7!!}{(2 \beta c)^4}+ \cdots} \right) \\ &= \frac{\pi}{\beta}\sqrt{ \frac{2 m}{c} }\left( {1 + \frac{3 f}{4 c^2 \beta}+ \frac{15 g^2}{16 c^3 \beta}+ \frac{105 f^2}{32 c^4 \beta}+ \cdots} \right).\end{aligned} \hspace{\stretch{1}}(1.8)

Retaining only the first two terms of the expansion, we have

\begin{aligned}\boxed{Z \approx\frac{\pi}{\beta}\sqrt{ \frac{2 m}{c} }\left( {1 + \frac{3 f}{4 c^2 \beta}+ \frac{15 g^2}{16 c^3 \beta}} \right).}\end{aligned} \hspace{\stretch{1}}(1.9)

Our average energy, in this approximation

\begin{aligned}\left\langle{{H}}\right\rangle &= -\frac{\partial {}}{\partial {\beta}} \ln Z \\ &\approx -\frac{\partial {}}{\partial {\beta}}\left( {- \ln \beta + \ln \left( 1 + \frac{3 f}{4 c^2 \beta}+ \frac{15 g^2}{16 c^3 \beta} \right) } \right) \\ &= \frac{1}{{\beta}} + \frac{1}{{\beta^2}} \frac{\frac{3 f}{4 c^2} + \frac{15 g^2}{16 c^3} }{1 + \frac{1}{{\beta}}\left( {\frac{3 f}{4 c^2} + \frac{15 g^2}{16 c^3} } \right)} \\ &= k_{\mathrm{B}} T + k_{\mathrm{B}}^2 T^2\left( { \frac{3 f}{4 c^2} + \frac{15 g^2}{16 c^3} } \right)\left( {1 - k_{\mathrm{B}} T \left( \frac{3 f}{4 c^2} + \frac{15 g^2}{16 c^3}  \right) + \cdots} \right).\end{aligned} \hspace{\stretch{1}}(1.10)

So to first order in T our specific heat is

\begin{aligned}C_{\mathrm{V}} = \frac{1}{{k_{\mathrm{B}}}} \frac{\partial {\left\langle{{H}}\right\rangle}}{\partial {T}}\approx1+ 2 k_{\mathrm{B}} T \left( {\frac{3 f}{4 c^2} + \frac{15 g^2}{16 c^3} } \right)\end{aligned} \hspace{\stretch{1}}(1.11)

or

\begin{aligned}\boxed{C_{\mathrm{V}} = 1+ k_{\mathrm{B}} T \left( {\frac{3 f}{2 c^2} + \frac{15 g^2}{8 c^3} } \right)+ \cdots}\end{aligned} \hspace{\stretch{1}}(1.12)

For the coordinate

\begin{aligned}\left\langle{{q}}\right\rangle &= \sqrt{\frac{2 \pi m}{\beta} }\frac{1}{{Z}}\int q e^{ -\beta \left( { c q^2 - g q^3 - f q^4} \right) } dq \\ &= \sqrt{\frac{2 \pi m}{\beta} }\frac{   \int q    e^{ -\beta c q^2 }   \left( { 1 + g \beta q^3 + f \beta q^4 } \right) dq}{   \frac{\pi}{\beta}   \sqrt{ \frac{2 m}{c} }   \left( { 1 + \frac{3 f}{4 c^2 \beta} + \frac{15 g^2}{16 c^3 \beta} } \right)} \\ &\approx \sqrt{\frac{\not{{2 \pi m}}}{\not{{\beta}}} }g \not{{\beta}} \frac{3!!}{(2 \beta c)^2} \sqrt{\frac{\not{{\pi}}}{\not{{\beta}} \not{{c}}}}\frac{1}{{   \frac{\not{{\pi}}}{\beta}   \sqrt{ \frac{\not{{2 m}}}{\not{{c}}} }}},\end{aligned} \hspace{\stretch{1}}(1.13)

or

\begin{aligned}\boxed{\left\langle{{q}}\right\rangle \approx\frac{3 g k_{\mathrm{B}} T}{4 c^2}.}\end{aligned} \hspace{\stretch{1}}(1.14)

Compare this to the expectation of the coordinate for an unperturbed harmonic oscillator

\begin{aligned}\left\langle{{q}}\right\rangle = \frac{\int q e^{ -\beta c q^2}}{\int e^{-\beta c q^2}} = 0.\end{aligned} \hspace{\stretch{1}}(1.15)

We now have a temperature dependence to the expecation of the coordinate that we didn’t have for the harmonic oscillator.

References

[1] RK Pathria. Statistical mechanics. Butterworth Heinemann, Oxford, UK, 1996.

Posted in Math and Physics Learning. | Tagged: , , , , , , , , , | Leave a Comment »