Try
Laplacian[Ez[\[Rho], \[Phi], z], {\[Rho], \[Phi], z}, "Cylindrical"] + n^2 \[Omega] k0^2 Ez[\[Rho], \[Phi], z] // Expand
Indeed! The solution is simple and correct. Thanks.
