Real Spherical Harmonics Visualization

Taking the  function as the example, the full mathematical form of the function is,

To visualize the function, we first need to make it clear what we are trying to visualize. First, physically, this functions gives the distribution of, e.g., the electron density, as the function of the polar angle  and the azimuthal angle  -- here it should be noted that it is actually the square of the function that gives the actual density but not the function value itself. Second, the distribution is only about the angular part of the distribution -- in another word, the radial part of the distribution is not of a concern here with the spherical harmonics functions. Given a certain combo of  and  values, we can a value of the function -- again, this is associated with, e.g., the electron density. But before squaring, the value can be positive or negative, representing the phase of the wavefunction. In 3D space, when we try to visualize the function, we get a series of  and  angles and for each of the combo of  and , we obtain a value which will then be taken as the radius value so to give a coordinate in the 3D space, i.e., .
Down below is the 3D visualization following such a scheme, for the  function presented above. The red color is for the positive lobe and the cyan color is for the negative lobe.

Ψ[θ_,ϕ_]:=Sin[θ]^3*Cos[θ]*Cos[3*ϕ]
plot3D=SphericalPlot3D[
Abs[Ψ[θ,ϕ]],{θ,0,π},{ϕ,0,2 π},
ColorFunction->Function[{x,y,z,θ,ϕ,r},If[Ψ[θ,ϕ]>0,Red,Cyan]],
ColorFunctionScaling->False,
Mesh->None,PlotPoints->50,PlotRange->All,
Boxed->False,Axes->False,
Lighting->"Neutral"
]

Viewing it from the top, we can clearly see the 6 lobes along the  dimension, and this is a typical shape of the -wave Fermi surface.

==============================================================================
To showcase the gradual change of the function value, we can use other color functions, e.g.,
==============================================================================

Ψ[θ_,ϕ_]:=Sin[θ]^3*Cos[θ]*Cos[3*ϕ]
plot3D=SphericalPlot3D[
Abs[Ψ[θ,ϕ]],{θ,0,π},{ϕ,0,2 π},
ColorFunction->Function[{x,y,z,θ,ϕ,r},ColorData["TemperatureMap"][Ψ[θ,ϕ]]],
ColorFunctionScaling->False,
Mesh->None,PlotPoints->50,PlotRange->All,
Boxed->False,Axes->False,
Lighting->"Neutral"
]

=====================================================
N.B. To get a full list of the available color functions, we can do,
=====================================================

ColorData["Gradients"]