Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Wyckoff Positions and Standard Cell #8

Open
Br0kenSmi1e opened this issue Dec 12, 2024 · 0 comments
Open

Wyckoff Positions and Standard Cell #8

Br0kenSmi1e opened this issue Dec 12, 2024 · 0 comments

Comments

@Br0kenSmi1e
Copy link

Br0kenSmi1e commented Dec 12, 2024

We process the C2DB database with spglib.get_layergroup. Positions of atoms are recorded as std_positions in spglib.
And things went wrong for some materials, such as https://c2db.fysik.dtu.dk/material/2AgNaTe-1 with layer group No. 64.
For Wyckoff position 2a, the coordinates given in https://www.cryst.ehu.es/ are (3/4, 1/4, 0) and (1/4, 3/4, 0).
However, the 2a Wyckoff positions in spglib are (1/2, 0, 0) and (0, 1/2, 0).

We can check the layer group operations of spglib.
The operations on input cell are stored as rotations and translations.
And the transformation between input cell and standard cell is stored as transformation_matrix and origin_shift.
So we can calculate the operations on standard cell, they are listed as following:
(x, y, z) (-x, y, z) (x, -y, z) (-x, -y, z) (y, x, z) (-y, x, z) (y, -x, z) (-y, -x, z) (x+1/2, -y+1/2, -z) (-y+1/2, -x+1/2, -z) (-x+1/2, y+1/2, -z) (y+1/2, x+1/2, -z) (-x+1/2, -y+1/2, -z) (y+1/2, -x+1/2, -z) (x+1/2, y+1/2, -z) (-y+1/2, x+1/2, -z).
This is different from the operations recorded in https://www.cryst.ehu.es/, which is
(x,y,z) (-x+1/2,-y+1/2,z) (-y+1/2,x,z) (y,-x+1/2,z) (-x,y+1/2,-z) (x+1/2,-y,-z) (y+1/2,x+1/2,-z) (-y,-x,-z) (-x,-y,-z) (x+1/2,y+1/2,-z) (y+1/2,-x,-z) (-y,x+1/2,-z) (x,-y+1/2,z) (-x+1/2,y,z) (-y+1/2,-x+1/2,z) (y,x,z).

These two groups of operations may be isomorphic to each other, i.e., they are just expressed in two different coordinate systems.
But the problem is how to check they are isomorphic, and find the isomorphism (the transformation between the two set of coordinate systems)?
Because CrystalFormer works under the coordinate system given by https://www.cryst.ehu.es/, so we need to perform the transformation so that CrystalFormer recognizes the output of spglib.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant