From Hv=μv, v†Hv=μ∥v∥2 is real, so μ∈R. For distinct eigenvalues, μivi†vj=vi†Hvj=μjvi†vj, hence orthogonality. Normalising the eigenvectors and taking them as the rows of U gives H=U†ΛU.
Pu=0, hence HPu=0. In the diagonal basis, f(λ)=1−∑juj2μj/(μj−λ)=−λ∑juj2/(μj−λ). Between consecutive μj its nonzero factor is strictly monotone from one infinite sign to the other, giving one root; zero weights give equality at endpoints. Hence μ1≤σ1≤μ2≤⋯≤σn−1≤μn.