Motion Captured – Springing Forth

August 1, 2014

Share this pageShare on FacebookTweet about this on TwitterShare on LinkedInGoogle+

This is the beginning of a regular series of short articles about basic (and not so basic) physics featuring some computer experiments. Below is the first: many point masses connected by springs. However, these are not regular springs. Try playing with them, the parameters $k$, $\alpha$ and $\beta$ are explained below, $N$ is the number of points on a side, so be aware of your processor power if you want to set it very high. Click “restart” to begin, “random” to reset the points and you can click on the red dots to move them around.

Go get a spring. The harder you squash or stretch it the harder it tries to get back to its “natural length”. So you know that the force must be related to the difference between the natural length and the length it’s been forced to i.e.
\[F = f(dx)\]
where $dx$ is the difference between the natural length and the stretched or squashed length. The equation above summarizes this in mathematical form, the force is given by some function of $dx$. 17th century weirdo Robert Hooke gave the first guess for what this function was in the form of a Latin anagram ceiiinosssttuv which is solved to give Ut tensio, sic vis, “As the extension, so the force” or “Force is proportional to the extension.” This gives us
\[F = k dx \]
where $k$ is some constant that tells us how stiff the spring is. Of course, as advanced 21st century scientists we can try a more complicated function,
\[F = k dx + \alpha dx^2 + \beta dx^3 \]
$\alpha$ and $\beta$ are just some numbers. Any force law other than Hooke’s one is called “non-linear”.

Now you could take a bunch of springs and join them together, making something like a mattress, pull on it and let it oscillate. For springs obeying Hooke’s law it will oscillate in a pretty boring way. For non-linear springs the motion will be crazy. Joining a bunch of non-linear springs together is called the Fermi-Pasta-Ulam problem. Fermi and Ulam are famous guys. Pasta was a American computer scientist with an unfortunate name. Also involved in the first paper was a lady called Mary Tsingou, but for some reason her name is not associated with the problem. Their idea was that if you had a non-linear force the system would oscillate in such a crazy way that it would go through every possible arrangement. When they actually simulated it however they found it didn’t, it oscillated in a surprisingly regular way. The conclusion in fancy physics jargon: non-linearity does not imply ergodicity.

From the force we can get a potential energy, using the non-linear force above gives us the following four basic shapes of the potential, the force is the slope of the curve up to some changes in the constants $k$, $\alpha$ and $\beta$.



You can think of the system as a little ball rolling on the potential hill. The top left is Hooke’s law, small oscillations about the mean value, imagine a little ball rolling back and forth inside the “well”. If you go to the simulation you can see what that looks like by setting $\alpha = \beta = 0$ and $k = 1$, for example. Much more interesting is the top right potential. Imagine carefully placing a little ball on top of the hump in the middle. Unless you were extremely careful it would roll down to the left or right. This is in fact exactly the sort of thing that happens in the Higgs mechanism which explains the mass of the electron. In our model, if we start almost at the top of the hump (press the “random” button and try $k=-1$, $\beta = 80$, $\alpha = 0$) we find that the springs want to compress together making clumps. These clumps are like our Higgs bosons. The symmetry, where all the points were more or less equal, is broken and some roll down the hill one way and some the other way.

The third one is also interesting, again imagine a little ball rolling on the hill on the bottom left. If it starts near the middle it will roll happily inside the little well, but if it gets above the hill it will roll down forever, getting faster and faster. You can see this happen, and break the simulation, by choosing say $k = 1$ and $\beta = -250$. After a little while some of the springs should, by chance, accumulate enough energy to get over the hump and fall to their doom. This is a metastability, interestingly our own universe could be metastable, but the hill to climb is certainly very large: see here. The bottom right potential is what it looks like if you turn $\beta = 0$ and turn on $\alpha$. You have one stable direction and one unstable one.

Try clicking on the dots, you can move them around if you want a different starting position. You can set up some interesting looking oscillations by choosing a good starting position for all the points. If you really want to, you can break some springs. Easiest is to make a potential with positive parameters. The forces can become huge and the algorithm that evolves the positions forward in time will break as the numbers overflow the maximum allowed size. These masses will break off. By pulling really hard on a spring you might be able to get it to snap off also, depending on the parameters. If you do break some springs clicking “random” will regenerate them.

Rudy Arthur