Normal mode calculations of trigonal selenium