Compute harmonic phonons, DOS, and band structure for a crystal structure using phonopy with NequIP-OAM-XL forces, returned as a .phonons file Ouro renders as an interactive dispersion beside the element-projected DOS. Useful for spotting imaginary modes before kappa or ZT estimation.