Some properly rotated lemniscates should *fake* p orbitals quite well. BTW
you can always try with some parametric equations (not that easy):
2px is given by the following parametric function (I'm using polar
coordinates):
1/4*pow(1/(2*pi),1/2)*pow(Z/a,3/2)*(Z*r/a)*exp(-Z*r/(2*a))*sin(theta)*sen(ph
i)
where Z is the atomic number and a = 0.529 e-10 m.
2py:
1/4*pow(1/(2*pi),1/2)*pow(Z/a,3/2)*(Z*r/a)*exp(-Z*r/(2*a))*cos(theta)
2pz:
1/4*pow(1/(2*pi),1/2)*pow(Z/a,3/2)*(Z*r/a)*exp(-Z*r/(2*a))*sin(theta)*cos(ph
i)

The parametric equations are (obviously):
x = r*sin(theta)*sin(phi);
y = r*sin(theta)*cos(phi);
z = r*cos(theta);

Nothing that I would really try. Good luck.

--
Jonathan.


"Paul Jones" <pdj### [at] psuedu> ha scritto nel messaggio
news:3c176f32@news.povray.org...
> Hello,
>
> Does anybody know how to create an isosurface for the atomic orbitals?
> Specifically the s, p and d (I can't find the equations for the f
> orbitals....). I have forgotten most of my quatum physics and even more
math
> (calc and diffeq).
>
> Any guidance or scene files would be much appreciated  :-)
>
> -thanks a lot
>
>
> paul
>
>

Post a reply to this message