Dear Johannes,
I found that hexagonal unit vectors do not work as expected with
ucf. First you need to generate a new cubic unit cell with
orthogonal vectors. You can do this with cif2cell tool (you can get
the cif file from any crystallographic database). The transformation
you need from hexagonal to cubic is (2x+y,2y,z). the command for
this is
cif2cell orig.cif --no-reduce --supercell=[[2,1,0],[0,2,0],[0,0,1]]
-p cif -o cubic.cif
Then you can use the new lattice constants and atom positions from
cubic.cif into ucf file with (1,0,0),(0,1,0),(0,0,1) unit vectors.
The difficulty comes with filling the interactions part of ucf file
with the increased number of atoms. For a single species, this can
be on a spreadsheet program. For more, you need more automation.
Please refer to gen_nbrs_example.py script I sent before.
Dear Richard,
Sorry I did not report this before, because I did not follow it
throughly. System generation from non-orthogonal unit cell
(hexagonal or rhombohedral) does not work as expected. So we need to
generate a cubic unit cell to make it work properly.
Best wishes,
Muftah
--
=============================================================
Postdoctoral fellow,
Crossing-Correlation Spin Electronics Lab (Sahashi group).
Dept. of Electronic Engineering, Tohoku University.
Research Center for Rare Metal and Green Innovation (RaMGI),
468-1, Aza-Aoba, Aramaki, Aoba-ku, Sendai 980-8579, JAPAN.
Tel/Fax: +81-22-752-2172
=============================================================