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
