We present a first-principles investigation of the localization properties of vibrational modes in amorphous silicon nitride (a-Si 3 N 4 ). Our investigation further confirms that the vibrational modes underlaying the peak at ∼471 cm −1 in the infrared spectrum of silicon nitride mainly consist of nitrogen motion in the direction normal to the nearest Si neighbors plane. Inplane stretching of N-Si bonds becomes largely dominant above ∼700 cm −1 . In particular vibrational modes underlaying the infrared peak at 825 cm −1 arise from Si-N bonds stretching motions. If N-N homopolar bonds were present, we show that N-N bond stretching occur above ∼1100 cm −1 . Furthermore, we investigate the localization properties of vibrational modes by calculating their inverse participation ratio (IPR) and phase quotient. From this analysis we infer that modes above 600 cm −1 shows a progressive increase of the localization degree and optic-like behavior, especially above 1000 cm −1 . At about 650 cm −1 , given the considerable IPR value, and on the basis of projectional anlysis on siliconbreathing-like motions of the NSi 3 units, we suggest that vibrational modes may involve correlated motions of neighboring SiN 4 tetrahedra.