R软件绘制3D量子轨道（球谐）

Question

我在R中使用surf3D或persp3D试图重现Maple的以下绘图功能：

>plot3d(Re(sqrt(15/(8*Pi))*(sin(theta)*cos(theta)*exp(I*phi)))^2,phi=0..2*Pi,theta=0..Pi,coords=spherical,scaling=constrained,style=patchcontour,numpoints=5000,axes=frame)

要获得与以下可见图相似（或更优）的东西：

在我开始自己对这个问题进行个人调查之前，有人能快速回答吗？

提前向您寻求帮助。

Answer 1

在这里，我认为解决方案是一块（或完整的？）（也许我们可以做得更好和更简单，例如，等距曲线和应该从中心到远而不是红色的颜色，而不仅仅是在垂直轴z上）：

library("plot3D")
theta<-seq(from = 0, to = pi, by = pi/320)
phi<-seq(from = 0, to = 2*pi, by = pi/160)

m<-length(phi); n=length(theta)
Phi<-matrix(rep(phi,each=n),nrow=n)
Theta<-matrix(rep(theta,m),nrow=n)


r = Re(1/2*sqrt(15/(2*pi))*exp(-1i*Phi)*sin(theta)*cos(Theta))
x = r*cos(Phi)
y = r*sin(Phi)
z = r*cos(Theta)

surf3D(x, y, z,colkey=TRUE,box=TRUE,bty="b2",main="A Spherical Harmonic")

结果给出：