Skip to main content

Decoupling harmonic oscillators

Wow!  Here's the short version of what's going on.  My quantum II final was a bit rough because I didn't know how to decouple a pair of coupled harmonic oscillators.  However, I found an excellent article on decoupling in my second favorite journal the American Journal of Physics.  I wanted to capture a few notes here on the process and on one of the math operations left out of the original article, but hashing out the level of background to provide has taken awhile.  I decided to go with the terser representation of the subject.  It may be extremely sketchy since I'm more than a bit scattered studying for finals.  Basically, these are my notes, and I hope they're helpful :)  Any ideas, comments, corrections, or suggestions, as always, are more than welcome! For more detail and an excellent article on how to decouple coupled oscillators see [1].

OK, here goes.

What does 'decoupling harmonic oscillators' mean and what is it used for?
Pretty frequently in physics, we come across problems where you have to decouple a system of coupled harmonics oscillators.  For folks who haven't had analytical mechanics yet, this can be a bit dicy.  I'm one of those folks.  Decoupling means taking the generic Lagrangian that describes the system of two coupled harmonic oscillators that looks like this (picture 1)



and modifying/rewriting it so that it looks like this (picture 2)


In doing this, we get rid of the terms with the b and q coefficients that mix the characteristics of the two oscillators and arrive at an equation where the energy of the system can be calculated easily and directly from the two new terms lambda_sub_1 and lambda_sub_2.

Why does everyone say 'Just diagonalize it'?
What they're referring to is that the above two equations can be written in matrix notation as (picture 3)



If you do the matrix multiplies through, you'll wind up back at the first equation above.  If we actually manage to modify the equation so that the b and q coefficients go to zero, then the two square matrices will be diagonal, they'll only have non-zero entries on a diagonal going from upper left to lower right.  Hence, diagonalizing is just another way of saying that we want to rearrange things so that b and q go to zero.

Why is the next thing everyone says 'Rotate the coordinates'
Looking back at the first equation, you might notice that it can be broken into two parts, each of which describe an ellipse when they're set to zero (picture 4).


If the b and q terms aren't zero, then you wind up with diagonal ellipses like those shown in (a) and (b) in picture 6 [1].  Paradoxically, diagonalizing the matrices and setting b and q to zero will force the ellipses to not be diagonal.  Aravind points out a very cool way to visualize the entire coupling process using the diagram to the left (picture 5).

First, we perform a rotation in the eta_dot plane to force the kinetic energy ellipse to no longer be diagonal.  Cool, the b term of our matrix that describes the kinetic energy of the system is now zero!  

As you can see however, the ellipse in the eta plane is still tilted off diagonal.  We need to rotate it as well.  There's a bit of a problem though.  When we rotate one of the ellipses we rotate them both.  To get around this, we rotate the potential energy ellipse and scale both ellipses at the same time in the same way so that the kinetic energy ellipse becomes a circle.  You can rotate a circle to your heart's content in any way you please and it's still a circle, and, (Cool!), it's still parallel to the coordinate axes in the way we want so that the b entries in the kinetic energy matrix are still zero.  When we're done, we've rotated the potential energy ellipse as well so that the q terms are now zero!

The step that Aravind kind of glossed over was the details on how to calculate the angle of rotation.  My version of the calculation follows.

The matrix equation that describes rotating the eta axes is shown in picture 6.  The squiggles on the right that are hard to distinguish are primes that indicate the new set of axes.


The matrix in the middle is called the rotation matrix and the angle theta is what we need to calculate.  We know that we want to rotate the axes so that the b components will become zero in the square matrix shown in picture 3.  We also know that the rotation matrix has to be applied to the square matrix from picture 3 in the following way, (picture 7), to bring it into the new coordinate system.  The square matrix with the a, b, and c components is denoted by T and the rotation matrix is denoted by F.  


We know that in the new version of T, we want the matrix entries where the b coefficients go to be zero.  So, to calulate the angle, we perform the matrix multiply, then set the resulting equation that winds up in that position to zero and solve.  Here goes (picture 8)


The first multiply of T times F is shown above and it gives the result below, (picture 9), that is multiplied by the transpose of F, (the term on the left).


The resulting equation for the coefficient in the lower left-hand corner of the T' matrix, (the position where the coefficient resides that we want to set to zero), and its solution for zero are shown below, (picture 10).  Gotta get back to studying now.



References:
1.  Geometrical interpretation of the simultaneous diagonalization of two
quadratic forms
http://dx.doi.org/10.1119%2F1.16069
Aravind P.K. (1989). Geometrical interpretation of the simultaneous diagonalization of two quadratic forms, American Journal of Physics, 57 (4) 309. DOI:

Comments

Popular posts from this blog

Cool Math Tricks: Deriving the Divergence, (Del or Nabla) into New (Cylindrical) Coordinate Systems

The following is a pretty lengthy procedure, but converting the divergence, (nabla, del) operator between coordinate systems comes up pretty often. While there are tables for converting between common coordinate systems, there seem to be fewer explanations of the procedure for deriving the conversion, so here goes!

What do we actually want?

To convert the Cartesian nabla



to the nabla for another coordinate system, say… cylindrical coordinates.



What we’ll need:

1. The Cartesian Nabla:



2. A set of equations relating the Cartesian coordinates to cylindrical coordinates:



3. A set of equations relating the Cartesian basis vectors to the basis vectors of the new coordinate system:



How to do it:

Use the chain rule for differentiation to convert the derivatives with respect to the Cartesian variables to derivatives with respect to the cylindrical variables.

The chain rule can be used to convert a differential operator in terms of one variable into a series of differential operators in terms of othe…

Lab Book 2014_07_10 More NaI Characterization

Summary: Much more plunking around with the NaI detector and sources today.  A Pb shield was built to eliminate cosmic ray muons as well as potassium 40 radiation from the concreted building.  The spectra are much cleaner, but still don't have the count rates or distinctive peaks that are expected.
New to the experiment?  Scroll to the bottom to see background and get caught up.
Lab Book Threshold for the QVT is currently set at -1.49 volts.  Remember to divide this by 100 to get the actual threshold voltage. A new spectrum recording the lines of all three sources, Cs 137, Co 60, and Sr 90, was started at approximately 10:55. Took data for about an hour.
Started the Cs 137 only spectrum at about 11:55 AM

Here’s the no-source background from yesterday
In comparison, here’s the 3 source spectrum from this morning.

The three source spectrum shows peak structure not exhibited by the background alone. I forgot to take scope pictures of the Cs137 run. I do however, have the printout, and…

Unschooling Math Jams: Squaring Numbers in their own Base

Some of the most fun I have working on math with seven year-old No. 1 is discovering new things about math myself.  Last week, we discovered that square of any number in its own base is 100!  Pretty cool!  As usual we figured it out by talking rather than by writing things down, and as usual it was sheer happenstance that we figured it out at all.  Here’s how it went.

I've really been looking forward to working through multiplication ala binary numbers with seven year-old No. 1.  She kind of beat me to the punch though: in the last few weeks she's been learning her multiplication tables in base 10 on her own.  This became apparent when five year-old No. 2 decided he wanted to do some 'schoolwork' a few days back.

"I can sing that song... about the letters? all by myself now!"  2 meant the alphabet song.  His attitude towards academics is the ultimate in not retaining unnecessary facts, not even the name of the song :)

After 2 had worked his way through the so…