pmtm (Multitaper) calculates PSD from Fourier coefficients as abs(Xx)^2; shouldn't that be Xx.conj(Xx)?
Understand PMTM multitaper PSD calculation in MATLAB! Is abs(xx)^2 correct for power spectral density from Fourier coefficients? Explore the xx * conj(xx) a