Dear Li, This one can be done using matdyn.x, too.
>My system has hexagonal symmetry and its primitive cell is rhombohedral. If the difference in between number of atoms not too much, I would prefer hexagonal, as it has higher symmetry and this will make things easier.
