In this simulation the bonding is done with a Lennard-Jones (LJ) potential. LJ potentials have extremely large repelling forces when the particles get too close and longer range attractive forces when the particles are farther apart. The LJ potential is of the form
`V(r)=4epsilon[(sigma/r)^2b-(sigma/r)^b]`
where `epsilon` is the potential minimum,`sigma` is a scaling parameter, `r` is the particle separation, and `b` is an integer that is adjustable by the learner. The larger the value of `b`, the more abrupt is the force change with distance, r. The solid is initially set up so that all the particles are separated by an `r` that causes the nearest neighbor forces to be zero. An obvious result of that set up is that alternate rows of particles have to be staggered by 1/2 spacing in the horizontal direction and the vertical spacing of rows is `sqrt(3)/2` times the horizontal spacing of particles in any row. The longitudinal acoustic wave is initiated by oscillating the vertical position of the top row by a given amplitude at a selected frequency. The transverse acoustic wave is initiated by oscillating the horizontal position of the top row by a given amplitude at a selected frequency.The details of the resulting waves in the rest of the solid are shown by plots of both the displacement sums of either the rows (blue) or the columns (green). For a stable solid, it is seen that about 40 rows are required to see a full wavelength of displacement. Since 50 rows is about one wavelength, `lambda`, and wave speed is frequency*`lambda` that means that the wave speed is about
`c=40omega_0r_(avg)/(2*pi) ~ 20`
pixels/time increment, where `omega_0` is the driving frequency expressed in radians per second, an `r_(avg)` is the initial separation between particles.The learner may find that certain combinations of the selectable parameters do not result in a stable solid. If that is the case, a restart can be accomplished by pressing the "Refresh Page" button.