Soft Spheres

This is a Brownian Dynamics simulation of soft spheres that interact with a Lennard-Jones potential in a periodic, cubic box. Just hit the "Start" button to begin the simulation, which runs in your browser.

You can adjust the simulation parameters by changing the JavaScript code below. Chiefly, you can experiment with the number of particles (and thereby the system density) and the interaction strength (parameter `epsilon`). You can also think of the interaction strength as being the inverse of temperature; i.e., a stronger interaction is consistent with a lower temperature.

Within the simulation we compute the radial pair density correlation function, g(r). You can see the plot of g(r) throughout time at the bottom of the page.

Steps: 0

// Note the box is a cube of length 1.0 var r = 0.05; SoftSpheres.run({ Nparticles: 1000, particleRadius: r, // LJ interaction strength epsilon: 0.00002, // LJ length scale sigma: 2 * r, // BD update length scale randomScale: 0.5 * r, // How frequently to update the neighbor lists in simulation cycleDuration: 50, // How frequently to sample the pair density distribution analysisRate: 10 });