Open kbarros opened 1 year ago
Another caution is that the Spglib documentation states this:
[Site symmetry symbols] are valid for the standard settings. For different settings and choices belonging to the same space group type, the same set of the symbols is returned.
Source. It's not clear to me whether this means that the returned symbols could be invalid for non-standard settings?
The same problem arose in a surprising way for Na4Ir3O8 (the iridium are arranged in a hyperkagome here)
a = 8.981
latvecs = lattice_vectors(a, a, a, 90, 90, 90)
x = 0.141
p = [0.125, x, x+1/4]
cryst = Crystal(latvecs, [p], 213) # Wyckoff 12d, point group ..2
@assert !isnothing(cryst.sitesyms)
x = 0
p = [0.125, x, x+1/4]
cryst = Crystal(latvecs, [p], 213) # Site symmetry is the same, but can't be inferred
@assert isnothing(cryst.sitesyms)
In the special case of x = 0
, it appears that spglib needs more crystallographic information to infer the spacegroup symmetry 213, and our "hack" breaks.
Just noting I mentioned the symmetry/Wyckoff analysis functions in doped
which could possibly be used for this goal in https://github.com/spglib/spglib/issues/135
This looks great, thanks for sharing! I wonder if the logic could be ported into Sunny/Julia without too much trouble.
If the chemical cell is known to be "conventional" for the given space group (in the sense of the Bilbao server), then it is possible to look up the data from tables, e.g., as available from the wyckoffs
function of Crystalline.jl: https://thchr.github.io/Crystalline.jl/dev/api/#Crystalline.wyckoffs-Union{Tuple{Integer},%20Tuple{D},%20Tuple{Integer,%20Val{D}}}%20where%20D
In this example,
Sunny builds the correct crystal, but cannot infer the site symmetry. The output is:
However, if we also include the iodine atoms, then
correctly outputs:
which indicates the point group of
-3m.
for Fe.The latter works by asking Spglib to infer all symmetry information about the crystal cell. From what I can tell, Spglib has no functionality to report point group symmetry information for a position and a given Hall number. See, e.g., https://github.com/spglib/spglib/issues/135.