I tried a nac calculated for In2O3 system, containing 40 atoms. When I used “amset phonon-frequency -o2 OUTCAR” to extracted POP frequency, “ValueError: cannot reshape array of size 120 into shape (3,3)” occurs. I found in phonon_frequency.py file, the eigenvectors calculated by VASP will be truned into a 33 matrix. However, in In2O3 system, the eigenvectors for a single q point will be a 1 120 vectors, so it cannot be turned to be 3*3 matrixs. It confused me, I hope I will get help from you. The OUTCAR, INCAR, POSCAR are attached.