A tool for finding and analyzing atomic environments in crystal structures.
The purpose of this tool is to find atomic environments in crystal structures. These environments can then be used to define measures of similarity and collective variables for enhanced sampling simulations. The output from this tool can be used directly to create reference environments for the EnvironmentSimilarity collective variable in PLUMED. The tool can also be used to visualize chemical environments around an atom for general purposes.
- Choose from the examples or upload your own configuration.
- Determine unique environments in crystal structures. Generally, these environments correspond to the basis of the crystal structure.
- Two algorithms for finding unique environments are available.
- A sorting algorithm that is fast but in some pathological cases can be wrong.
- A permutation algorithm that is slow but always gives correct results. My advice is to use the fast version and change to the slow one if the algorithm finds more unique environments than expected.
- Visualize the environments.
- Output the environments in Protein Data Bank (pdb) format.
The app has four tabs:
- Choose configuration
- Define environments
- Analyze environments
- Output environments
that should be used sequentially.
The Choose configuration tab has a dropdown menu and an Upload button. By default the user can select a configuration from some examples that are useful to play with and get a grasp of the tool. The Upload button should be used to analyze a configuration stored in the user's hard drive. Once a configuration file has been uploaded the file should appear in the dropdown choices. Any file format supported by the ASE library can be used. A complete list can be found here. ASE determines the file format by its extension so be sure to name the file properly!
The Define environments tab is where the user defines the environments. There are three different ways to define the environments, Type, String, and Step that can be selected with the toggle buttons. Typically the most useful one is the Type option. In this case one chooses the atom type of the central atom, the atom type of the neighboring atoms, and a cutoff. These parameters are set using the Central atoms type:, Neighbor atoms type:, and Cutoff (Å): boxes. The app will find all environments with atoms of the second type around atoms of the first type inside the specified cutoff.
The tool then calculates the unique environments. For this purpose it compares all environments and keeps only the ones that are different from each other. Two environments are considered equal if the distance of every neighboring atom differs by less than some tolerance. The tolerance is specified in the Tolerance (Å): box. In order to compare environments, all permutations of the atoms have to be considered and that calculation can be slow. The Fast algorithm checkbox switches between a sorting algorithm that is fast but in some pathological cases can be wrong, and a permutation algorithm that is slow but always gives correct results. My advice is to use the fast version and change to the slow one if the algorithm finds more unique environments than expected.
When the calculation has finished the app prints the number of environments that have been found and the average number of neighbors in the environments.
- Cutoff should be given in angstrom
- Tolerance should be given in angstrom
The Analyze environments tab allows to visualize the calculated environments. One can toggle between the Unique environments and All of them. There is a slide bar to select the environment to visualize.
The Output environments tab prints the environments in Protein Data Bank (PDB) format. One can toggle between the Unique environments and All of them. There is a slide bar to select the environment to visualize.
- The app uses several python libraries, for instance ASE, NGLVIEW, and ipywidgets.
- I am grateful to Giovanni Pizzi and Dou Du for suggesting to deploy the tool using Binder+appmode.
- This tool was developed with support of the Swiss National Science Foundation (SNSF) through an Early Postdoc.Mobility fellowship.
- I also acknowledge funding from the NCCR MARVEL funded by the SNSF and from the CSI Computational Science Center funded by the Department of Energy of the USA.
How to cite
If you are using this tool to find environments for enhanced sampling simulations please read and cite: