8.1.8. Bonded particle modelsΒΆ

The BPM package implements bonded particle models which can be used to simulate mesoscale solids. Solids are constructed as a collection of particles which each represent a coarse-grained region of space much larger than the atomistic scale. Particles within a solid region are then connected by a network of bonds to provide solid elasticity.

Unlike traditional bonds in molecular dynamics, the equilibrium bond length can vary between bonds. Bonds store the reference state. This includes setting the equilibrium length equal to the initial distance between the two particles but can also include data on the bond orientation for rotational models. This produces a stress free initial state. Furthermore, bonds are allowed to break under large strains producing fracture. The examples/bpm directory has sample input scripts for simulations of the fragmentation of an impacted plate and the pouring of extended, elastic bodies.


Bonds can be created using a read data or create bonds command. Alternatively, a molecule template with bonds can be used with fix deposit or fix pour to create solid grains.

In this implementation, bonds store their reference state when they are first computed in the setup of the first simulation run. Data is then preserved across run commands and is written to binary restart files such that restarting the system will not reset the reference state of a bond. Bonds that are created midway into a run, such as those created by pouring grains using fix pour, are initialized on that timestep.

As bonds can be broken between neighbor list builds, the special_bonds command works differently for BPM bond styles. There are two possible settings which determine how pair interactions work between bonded particles. First, one can turn off all pair interactions between bonded particles. Unlike bond quartic, this is not done by subtracting pair forces during the bond computation but rather by dynamically updating the special bond list. This is the default behavior of BPM bond styles and is done by updating the 1-2 special bond list as bonds break. To do this, LAMMPS requires newton bond off such that all processors containing an atom know when a bond breaks. Additionally, one must do either (A) or (B).

  1. Use the following special bond settings

    special_bonds lj 0 1 1 coul 1 1 1
    

    These settings accomplish two goals. First, they turn off 1-3 and 1-4 special bond lists, which are not currently supported for BPMs. As BPMs often have dense bond networks, generating 1-3 and 1-4 special bond lists is expensive. By setting the lj weight for 1-2 bonds to zero, this turns off pairwise interactions. Even though there are no charges in BPM models, setting a nonzero coul weight for 1-2 bonds ensures all bonded neighbors are still included in the neighbor list in case bonds break between neighbor list builds.

  2. Alternatively, one can simply overlay pair interactions such that all bonded particles also feel pair interactions. This can be accomplished by using the overlay/pair keyword present in all bpm bond styles and by using the following special bond settings

    special_bonds lj/coul 1 1 1
    

See the Howto page on broken bonds for more information.


Currently there are two types of bonds included in the BPM package. The first bond style, bond bpm/spring, only applies pairwise, central body forces. Point particles must have bond atom style and may be thought of as nodes in a spring network. Alternatively, the second bond style, bond bpm/rotational, resolves tangential forces and torques arising with the shearing, bending, and twisting of the bond due to rotation or displacement of particles. Particles are similar to those used in the granular package, atom style sphere. However, they must also track the current orientation of particles and store bonds and therefore use a bpm/sphere atom style. This also requires a unique integrator fix nve/bpm/sphere which numerically integrates orientation similar to fix nve/asphere.

To monitor the fracture of bonds in the system, all BPM bond styles have the ability to record instances of bond breakage to output using the dump local command. Additionally, one can use compute nbond/atom to tally the current number of bonds per atom.

In addition to bond styles, a new pair style pair bpm/spring was added to accompany the bpm/spring bond style. This pair style is simply a hookean repulsion with similar velocity damping as its sister bond style.


While LAMMPS has many utilities to create and delete bonds, only the following are currently compatible with BPM bond styles:

Note create_bonds requires certain special_bonds settings. To subtract pair interactions, one will need to switch between different special_bonds settings in the input script. An example is found in examples/bpm/impact.