Jaap Frank wrote:

> the attached scene-files produce pictures of the
> quantum wave functions of hydrogen.
>
> [...]

Thanks, I see it's too complicated for a quick look, but I'll study it when 
I have a little more time. In case you're interested I attach my own 
version of hydrogen-like orbitals, I use a macro to define the orbital 
functions, and these functions can be used with media or isosurfaces.
It seemed to work faster if I defined the functions and the normalization 
factors separately.

A simple example:

DefineOrbital(2,0,0,2,0.05) // 2s orbital for Z=2, threshold=0.05

#declare c1 = <-10,-10,0>;
#declare c2 = <10,10,10>;

isosurface {
  function { Orb_pos(x,y,z) }
  contained_by { box { c1,c2 } }
  accuracy 0.001
  evaluate 0.7, 1.0, 1.0
  all_intersections
  pigment { color rgb <0,1,0> }
  finish { brilliance 2.0 phong 0.5 phong_size 100 }
}

isosurface {
  function { Orb_neg(x,y,z) }
  contained_by { box { c1,c2 } }
  accuracy 0.001
  evaluate 0.7, 1.0, 1.0
  all_intersections
  pigment { color rgb <1,0,0> }
  finish { brilliance 2.0 phong 0.5 phong_size 100 }
}

P.S. I neglected the 4*pi factor in the probability density.
P.P.S. Can an isosurface be declared and then the function be altered?

-- 

Linux User #289967 (counter.li.org)
PGP Pub Key ID: 0x01A95F99 (pgp.escomplinux.org)

Post a reply to this message