find a point-defect state

Dear all,

I have several questions to discuss you.
 
Q1. In the section of Finding a Point-defect State of  User Tutorial, I did almost what it said and get the reasonable result. It writes that we needn't to repeat the run to output the fields and we can just do (output-field-z 25). At this point, I can't understand. I put the function (output-field-z 25) in a control file, say, 1.ctl and write mpb 1.ctl in command line, but I failed. I didn't succeed when I just write mpb (output-field-z 25) command line, too. Could you help me?
 
Q2. Still, in the section of Finding a Point-defect State of  User Tutorial, you write that we can get the same result by adjust the control file, 
( set!  num-bands  1 )
( set!  target-freq  ( /  ( +  0.2812   0.4174 )  2 ) )
( set!  tolerance  1e-8 )
In this case, target-freq is 0.3493. The difference between 0.3493 and 0.378166( the result ) is not small. So why tolerance is so small, only 1e-8?
 
Q3. I want to have a trial of coordinate conversion functions, so I put  ( reciprocal->cartesian vector3 (0.1 0.3 0) ) in sq-rods.ctl which is attached. I can't find a result in sq-rods.out. Could you explain how to use these functions?

Thank you very much.

Best Regards,
Young Soo Lee

************************************************************************************************************
(define-param r 0.2) ; radius of the rods
(define-param eps 12) ; dielectric constant

(define GaAs (make dielectric (epsilon eps)))

(set! geometry-lattice (make lattice (size 5 5 no-size))) ; 2d cell

(set! geometry 
      (list
       (make cylinder 
     (material GaAs) 
     (center 0 0) (radius r) (height infinity))))

(set! geometry (geometric-objects-lattice-duplicates geometry))

(set! geometry (append geometry
                      (list (make cylinder (center 0 0)
                                  (radius 0.2) (height infinity)
                                  (material air)))))
      
(set-param! resolution 16)
(set! mesh-size 20)
(set-param! num-bands 26)

(set! k-points (list (vector3 0.5 0.5 0)))
(reciprocal->cartesian (vector3 0.1 0.3 0))

(run-tm)

__________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com

_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss