I want to calculate cluster sizes for a box of particles using the following method:

# Compute cluster IDs for atoms

compute myCluster all cluster/atom 0.9

# Define chunks based on cluster IDs

compute myChunks all chunk/atom c_myCluster compress yes

# Compute the number of atoms in each chunk

compute sizeChunks all property/chunk myChunks count

fix chunkSizes all ave/time 10000 1 10000 c_sizeChunks file chunk_sizes.dat mode vector

However as LAMMPS documentation said,

A cluster is defined as a set of atoms, each of which is within the cutoff distance from one or more other atoms in the cluster. If an atom has no neighbors within the cutoff distance, then it is a 1-atom cluster.

I want to discard one atom clusters. More precisely I want to define cluster only if number of particles in cluster becomes more than a cutoff (say, N=10).

More Shuvadip Dutta's questions See All
Similar questions and discussions