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).