Why do the two methods give different answers for the Fourier transform of $Y = \cos(\omega_0 t + \phi)$?
They are the same.

Using the Dirac impulse sifting property $f(x) \delta(x-a) = f(a)\delta(x-a)$ you can also verify that the second method produces the same output as the first method.