\(\renewcommand{\AA}{\text{Å}}\)

delete_atoms command

Syntax

delete_atoms style args keyword value ...
  • style = group or region or overlap or random or variable

    group args = group-ID
    region args = region-ID
    overlap args = cutoff group1-ID group2-ID
      cutoff = delete one atom from pairs of atoms within the cutoff (distance units)
      group1-ID = one atom in pair must be in this group
      group2-ID = other atom in pair must be in this group
    random args = ranstyle value eflag group-ID region-ID seed
      ranstyle = fraction or count
        for fraction:
          value = fraction (0.0 to 1.0) of eligible atoms to delete
          eflag =  no for fast approximate deletion, yes for exact deletion
        for count:
          value = number of atoms to delete
          eflag = no for warning if count > eligible atoms, yes for error
      group-ID = group within which to perform deletions
      region-ID = region within which to perform deletions
                  or NULL to only impose the group criterion
      seed = random number seed (positive integer)
    variable args = variable-name
  • zero or more keyword/value pairs may be appended

  • keyword = compress or bond or mol

    compress value = no or yes
    bond value = no or yes
    mol value = no or yes

Examples

delete_atoms group edge
delete_atoms region sphere compress no
delete_atoms overlap 0.3 all all
delete_atoms overlap 0.5 solvent colloid
delete_atoms random fraction 0.1 yes all cube 482793 bond yes
delete_atoms random fraction 0.3 no polymer NULL 482793 bond yes
delete_atoms random count 500 no ions NULL 482793
delete_atoms variable checkers

Description

Delete the specified atoms. This command can be used, for example, to carve out voids from a block of material or to delete created atoms that are too close to each other (e.g., at a grain boundary).

For style group, all atoms belonging to the group are deleted.

For style region, all atoms in the region volume are deleted. Additional atoms can be deleted if they are in a molecule for which one or more atoms were deleted within the region; see the mol keyword discussion below.

For style overlap pairs of atoms whose distance of separation is within the specified cutoff distance are searched for, and one of the 2 atoms is deleted. Only pairs where one of the two atoms is in the first group specified and the other atom is in the second group are considered. The atom that is in the first group is the one that is deleted.

Note that it is OK for the two group IDs to be the same (e.g., group all), or for some atoms to be members of both groups. In these cases, either atom in the pair may be deleted. Also note that if there are atoms which are members of both groups, the only guarantee is that at the end of the deletion operation, enough deletions will have occurred that no atom pairs within the cutoff will remain (subject to the group restriction). There is no guarantee that the minimum number of atoms will be deleted, or that the same atoms will be deleted when running on different numbers of processors.

For style random a subset of eligible atoms are deleted. Which atoms to delete are chosen randomly using the specified random number seed. Which atoms are deleted may vary when running on different numbers of processors.

For ranstyle = fraction, the specified fractional value (0.0 to 1.0) of eligible atoms are deleted. If eflag is set to no, then the number of deleted atoms will be approximate, but the operation will be fast. If eflag is set to yes, then the number deleted will match the requested fraction, but for large systems the selection of deleted atoms may take additional time to determine.

For ranstyle = count, the specified integer value is the number of eligible atoms are deleted. If eflag is set to no, then if the requested number is larger then the number of eligible atoms, a warning is issued and only the eligible atoms are deleted instead of the requested value. If eflag is set to yes, an error is triggered instead and LAMMPS will exit. For large systems the selection of atoms to delete may take additional time to determine, the same as for requesting an exact fraction with pstyle = fraction.

Which atoms are eligible for deletion for style random is determined by the specified group-ID and region-ID. To be eligible, an atom must be in both the specified group and region. If group-ID = all, there is effectively no group criterion. If region-ID is specified as NULL, no region criterion is imposed.

New in version 4May2022.

For style variable, all atoms for which the atom-style variable with the given name evaluates to non-zero will be deleted. Additional atoms can be deleted if they are in a molecule for which one or more atoms were deleted within the region; see the mol keyword discussion below. This option allows complex selections of atoms not covered by the other options listed above.


Here is the meaning of the optional keywords.

If the compress keyword is set to yes, then after atoms are deleted, then atom IDs are re-assigned so that they run from 1 to the number of atoms in the system. Note that this is not done for molecular systems (see the atom_style command), regardless of the compress setting, since it would foul up the bond connectivity that has already been assigned. However, the reset_atoms id command can be used after this command to accomplish the same thing.

Note that the re-assignment of IDs is not really a compression, where gaps in atom IDs are removed by decrementing atom IDs that are larger. Instead the IDs for all atoms are erased, and new IDs are assigned so that the atoms owned by individual processors have consecutive IDs, as the create_atoms command explains.

A molecular system with fixed bonds, angles, dihedrals, or improper interactions, is one where the topology of the interactions is typically defined in the data file read by the read_data command, and where the interactions themselves are defined with the bond_style, angle_style, etc. commands. If you delete atoms from such a system, you must be careful not to end up with bonded interactions that are stored by remaining atoms but which include deleted atoms. This will cause LAMMPS to generate a “missing atoms” error when the bonded interaction is computed. The bond and mol keywords offer two ways to do that.

It the bond keyword is set to yes then any bond or angle or dihedral or improper interaction that includes a deleted atom is also removed from the lists of such interactions stored by non-deleted atoms. Note that simply deleting interactions due to dangling bonds (e.g., at a surface) may result in a inaccurate or invalid model for the remaining atoms.

It the mol keyword is set to yes, then for every atom that is deleted, all other atoms in the same molecule (with the same molecule ID) will also be deleted. This is not done for atoms with molecule ID = 0, since such an ID is assumed to flag isolated atoms that are not part of molecules.

Note

The molecule deletion operation is invoked after all individual atoms have been deleted using the rules described above for each style. This means additional atoms may be deleted that are not in the group or region, that are not required by the overlap cutoff criterion, or that will create a higher fraction of porosity than was requested.

Restrictions

The overlap styles requires inter-processor communication to acquire ghost atoms and build a neighbor list. This means that your system must be ready to perform a simulation before using this command (force fields setup, atom masses set, etc.). Since a neighbor list is used to find overlapping atom pairs, it also means that you must define a pair style with the minimum force cutoff distance between any pair of atoms types (plus the neighbor skin) \(\ge\) the specified overlap cutoff.

If the special_bonds command is used with a setting of 0, then a pair of bonded atoms (1–2, 1–3, or 1–4) will not appear in the neighbor list, and thus will not be considered for deletion by the overlap styles. You probably do not want to delete one atom in a bonded pair anyway.

The bond yes option cannot be used with molecular systems defined using molecule template files via the molecule and atom_style template commands.

Default

The option defaults are compress = yes, bond = no, mol = no.