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.