Re: mpb-discuss Digest, Vol 18, Issue 1

Why don't you transfer the vector from real space to k-space.
 
I've spent almost 4 months as your question.
 
Please check a solid state physics book and you see the reciprocal lattice.
 
when you chek it, set the vector a1,a2,a3 as 
a1 = ((sqrt 3)/2 , 0.5, 0)
a2 = ((sqrt 3)/2 , -0.5, 0)
a3 = ( 0, 0, c)
 
then 
b1 = (1/(sqrt3), 1, 0)
b2 = (1/(sqrt3), -1, 0)
b3 = this is not important. 
 
Good luck.
 
 
-----Original Message-----
From: [email protected] 
To: [email protected]
Cc: 
Sent: 07-08-02(목) 01:00:02
Subject: mpb-discuss Digest, Vol 18, Issue 1
Send mpb-discuss mailing list submissions to
[email protected]
To subscribe or unsubscribe via the World Wide Web, visit
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
or, via email, send a message with subject or body 'help' to
[email protected]
You can reach the person managing the list at
[email protected]
When replying, please edit your Subject line so it is more specific
than "Re: Contents of mpb-discuss digest..."
Today's Topics:
1. One simple question as a new MPB user (bing bai)
----------------------------------------------------------------------
Message: 1
Date: Wed, 1 Aug 2007 23:25:02 +0800 (CST)
From: bing bai <[email protected]>
Subject: [MPB-discuss] One simple question as a new MPB user
To: [email protected]
Message-ID: <[email protected]>
Content-Type: text/plain; charset="gb2312"
Hi,everyone
I am a new MPB user. When I try to simulate the ctl files in examples, I have a problem.
When I am trying the tri-rods ctl file, I could not understand why the K wavevector in irreducible Brillouin zone is (vector3 -1/3, 1/3, 0). I chceked the book,"Molding the flow of light",it seems that the K point should be (vector3 1/3, 1/3,0), and I can not find out where I was wrong. Please help me!
Thanks a lot!
Bing
; A triangular lattice of dielectric rods in air. (This structure has
; a band-gap for TM fields.) This file is used in the "Data Analysis
; Tutorial" section of the MPB manual.
(set! num-bands 8)
(set! geometry-lattice (make lattice (size 1 1 no-size)
(basis1 (/ (sqrt 3) 2) 0.5)
(basis2 (/ (sqrt 3) 2) -0.5)))
(set! geometry (list (make cylinder
(center 0 0 0) (radius 0.2) (height infinity)
(material (make dielectric (epsilon 12))))))
(set! k-points (list (vector3 0 0 0) ; Gamma
(vector3 0 0.5 0) ; M
(vector3 (/ -3) (/ 3) 0) ; K
(vector3 0 0 0))) ; Gamma
(set! k-points (interpolate 4 k-points))
(set! resolution 32)
(run-tm (output-at-kpoint (vector3 (/ -3) (/ 3) 0) 
fix-efield-phase output-efield-z))
(run-te)
---------------------------------
????????3.5G???20M??? 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://ab-initio.mit.edu/pipermail/mpb-discuss/attachments/20070801/b40aab54/attachment.html 
------------------------------
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
End of mpb-discuss Digest, Vol 18, Issue 1
******************************************

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