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.

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
});