1) prepare your POSCAR file. 2) phonopy --symmetry --tolerance 0.01 -c POSCAR 3) POS2ABR > ABR.out (converting PPOSCAR to POSCAR_std and generating ABRs) (* paste PPOSCAR below or download the source code *) 4) cp POSCAR_std POSCAR; do VASP/QE calculations. 5) IRVSP -sg $sgn [ -nb $m $n ] > outir (generating tqc.txt/tqc.data when using only maximal HSKPs) 6) solve EBR and ABR decompositions (using tqc.data and PPOSCAR). (* paste tqc.data below *) Note: please fill in both boxes above! (only valid without spin-orbit coupling)