Normal mode analysis is a typical topic in junior/senior mechanics courses. Ours suffers from a lack of linear algebra as a prerequisite so I’ve worked to find ways to engage students with this material without that background. My typical approach is to model the system with a ton of friction so that it settles down to an equilibrium setup, then turn friction off and move the particles a little bit away from their equilibrium, then analyze the motions, looking for peaks in the Fourier transforms.
Last week I was a part of a chemistry honors defense that was looking at Raman spectra of bath salts. It was mostly a conversation about normal modes! So I wanted to see if I could model a benzene ring using this approach. Here’s my first try at it (I know, I know, I haven’t added the H’s around the edge. Babysteps!)
You can see that it works pretty well, but that the top “breathing” mode doesn’t look perfectly symmetrical.
Ok, here’s where I need your help. To find the resonant frequencies of the system, I simply Fourier transformed the time sequence of the distance one of the particles was away from its equilibrium position. What I’m wondering is whether there’s a better time sequence to use to do that Fourier analysis. I could certainly:
- do that trick for all the particles and, I don’t know, add the results
- fourier transform the time series of all the coordinates (that would be 12 for the case above – x and y for all particles). Then add them? That seems problematic because symmetric modes will cancel. Maybe add their amplitudes?
- Look for a measure of the whole system (like the standard deviation of all the displacements or something) that could then be Fourier transformed
I guess I don’t know what’s best. Hence this blog post. Any help you could give would be great. Here are some starters for you:
- You’re an idiot to teach this without linear algebra. But, if you have to, here’s what I would suggest . . .
- That’s not a benzene ring it’s a polysyllabatesupercomplicatedchemistryword and those modes don’t look right.
- How do you deal with degeneracies?
- How long does it take to run that?
- Can you share your code?
- Those animations are wrong because they all take the same amount of time to cycle through. Instead you should . . .
- I’ve tried modeling something like this before and all the particles stack up. The lengths are all at their equilibrium length, but the angles are all basically zero. How did you avoid that?
- quick answer: angular springs
- It’s obvious! You should fourier transform the time series of . . .