When using DFPT or the finite displacement (also known as the frozen phonon) method to compute surface phonon spectra, it is generally necessary to perform high-precision optimization on the slab structure.
My question is, when optimizing this slab structure, should the atoms in the bottom layers of the slab be fixed as usual?