GNUPLOT

- not so Frequently Asked Questions -

update 2005/11/19

Deformed Nucleus (Legendre Expansion)

[ver.4] ONLY !

The spherical harmonics Y[lm](theta,phi) is reduced to a simple Legendre function p[l](cos(theta)) scaled by a constant when m=0, which is independent of phi. The following equations and figure are Y[l](theta) for several l-values.

eq2 yfunc

Using this Y[l](theta) with l=lambda=even terms, shape of a deformed nucleus can be expanded as follows:

eq2

where the beta a parameter of deformation. If beta=0, the nucleus is spherical. The 3-dim. shape given by this equation is shown with gnuplot. As it is already shown in the previous section, we express the (x,y,z) coordinate with the angles u,v and radius r.

   x = R(theta)*cos(u)*cos(v)
   y = R(theta)*sin(u)*cos(v)
   z = R(theta)*sin(v)

where the theta is the angle measured from the Z-axis, so that the relation between theta and v is theta = pi/2-v. To draw the surface, the parameters u,v are varied from 0 to 360 deg. In the case of beta_2 = 0.3, beta_4 = 0.1, and R_0 = 1:

gnuplot> set parametric

        dummy variable is t for curves, u/v for surfaces
gnuplot> set angle degree
gnuplot> set urange [0:360]
gnuplot> set vrange [0:360]
gnuplot> set isosample 16,16
gnuplot> set ticslevel 0
gnuplot> set view 75,25
gnuplot> set size 0.7,1.0
gnuplot> set xrange [-2:2]
gnuplot> set yrange [-2:2]
gnuplot> set zrange [-2:2]
gnuplot> set urange [0:360]
gnuplot> set vrange [0:360]
gnuplot> y0(t)=1.0
gnuplot> y2(t)=sqrt(5.0/(4*pi))*( 3.0*cos(t)**2 -  1.0             )/2.0
gnuplot> y4(t)=sqrt(9.0/(4*pi))*(35.0*cos(t)**4 - 30*cos(t)**2 +3.0)/8.0
gnuplot> b2=0.3
gnuplot> b4=0.1
gnuplot> r(t) = 1 + b2*y2(0.5*pi-t) + b4*y4(0.5*pi-t)
gnuplot> fx(u,v)=cos(u)*cos(v)
gnuplot> fy(u,v)=sin(u)*cos(v)
gnuplot> fz(v)=sin(v)
gnuplot> set pm3d
gnuplot> splot r(v)*fx(u,v),r(v)*fy(u,v),r(v)*fz(v) with lines
legendre

The deformation parameters, beta_2 and beta_4 can be positive or negative. Here are some examples for some combinations of beta_2 and beta_4. The beta_2 parameters are taken to be -0.4 or 0.4, and for each beta_2, we changed the beta_4 value from -0.2 to 0.2. When beta_2 is positive the shape of nucleus is prolate, while it becomes oblate if beta_2 is negative.

gnuplot> set border 0
gnuplot> unset xtics
gnuplot> unset ytics
gnuplot> unset ztics
gnuplot> unset colorbox
gnuplot> b2 = -0.4 ; b4 = -0.2 ; replot ; pause -1
gnuplot> b2 = -0.4 ; b4 =  0.0 ; replot ; pause -1
gnuplot> b2 = -0.4 ; b4 =  0.2 ; replot ; pause -1
gnuplot> b2 =  0.4 ; b4 = -0.2 ; replot ; pause -1
gnuplot> b2 =  0.4 ; b4 =  0.0 ; replot ; pause -1
gnuplot> b2 =  0.4 ; b4 =  0.2 ; replot ; pause -1

beta_4 = -0.2 beta_4 = 0 beta_4 = 0.2
beta_2 = -0.4 legendre1 legendre2 legendre3
beta_2 = 0.4 legendre4 legendre5 legendre6
up