Skip to content

Compare two descriptions of crystal structures to see if the have the same underlying structure

License

Notifications You must be signed in to change notification settings

alliepiper/XtalComp

Repository files navigation

======================================================================
What is XtalComp for?

Computationally identifying duplicate crystal structures taken from
the output of modern solid state calculations is a non-trivial
exercise for many reasons. The translation vectors in the description
are not unique -- they may be transformed into linear combinations of
themselves and continue to describe the same extended structure. The
coordinates and cell parameters contain numerical noise. The periodic
boundary conditions at the unit cell faces, edges, and corners can
cause very small displacements of atomic coordinates to result in very
different representations. The positions of all atoms may be uniformly
translated by an arbitrary vector without modifying the underlying
structure. Additionally, certain applications may consider
enantiomorphic structures to be identical.

The XtalComp algorithm overcomes these issues to detect duplicate
structures regardless of differences in representation. It begins by
performing a Niggli reduction on the inputs, standardizing the
translation vectors and orientations. A transform search is performed
to identify candidate sets of rotations, reflections, and translations
that potentially map the description of one crystal onto the other,
solving the problems of enantiomorphs and rotationally degenerate
lattices. The atomic positions resulting from each candidate transform
are then compared, using a cell-expansion technique to remove periodic
boundary issues. Computational noise is treated by comparing
non-integer quantities using a specified tolerance.


======================================================================
How to use XtalComp

mkdir build
cd build
ccmake ..
make

This will build the XtalComp static library and a sample testing
script. To use XtalComp in your application, link to the statically
built XtalComp library in the "build" directory and include
xtalcomp.h in your application. See test.cpp for example usage.

To run the sample test script, execute

./test

The output will indicate if any tests have failed.


======================================================================
Files included:

CMakeLists.txt     : Build system config
xtalcomp.*         : Implementation of the XtalComp algorithm
stablecomparison.* : Convenience functions for comparing floating
                       point number to a tolerance.
xcvector.*         : Simple column 3-vector class
xcmatrix.*         : Simple 3x3 matrix class
xctransform.*      : Simple class transform class
test.*             : Sample test script

About

Compare two descriptions of crystal structures to see if the have the same underlying structure

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Contributors 3

  •  
  •  
  •  

Languages