Compute harmonic phonons, DOS, and band structure for a crystal structure using phonopy with NequIP-OAM-XL forces. Useful for spotting imaginary modes before kappa or ZT estimation.