Hi Onder,
I'm not sure if there is any difference between triangular lattice
and hexagonal lattice. They look exactly the same to me, geometrically.
Compare to the tri-rods.ctl, your basis is just rotated 30 degree
counter-clockwise.
For arbitrary lattice vectors,I remember that the basis of the k vector is
given in the header of the .out file. Eg,
Lattice vectors:
(0.939693, 0.34202, 0)
(0, 1, 0)
(0, 0, 1)
Cell volume = 0.939693
Reciprocal lattice vectors (/ 2 pi):
(1.06418, -0, 0)
(-0.36397, 1, -0)
(0, -0, 1)
Then the k-points are expressed in this set of reciprocal lattice basis. You
can just run some test files to see what's the result after you specify new
lattice vectors.
Best regards
Boshen
On Mon, Apr 25, 2011 at 12:00 PM, <[email protected]>wrote:
> 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. hex.PhC k-points (?zg?r ?nder KARAKILIN?)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 25 Apr 2011 16:20:32 +0300
> From: ?zg?r ?nder KARAKILIN? <[email protected]>
> Subject: [MPB-discuss] hex.PhC k-points
> To: [email protected]
> Message-ID:
> <[email protected]>
> Content-Type: text/plain;charset=iso-8859-9
>
> Dear mpb users
>
> I examine the tri-rods.ctl example in mpb wiki page.
> This real space structure is triangular and corresponded reciprocal
> structure is hexagonal, isn’t it?
> When I change the real space basis for hexagonal structure as the
> following code, what will the k-points become?
> Does the k-point remain still same?
> best...
> Onder
>
>
>
> (define-param npoly 1.46)
> (define-param epspoly (* npoly npoly))
> (define-param epsair 1)
>
> (define-param a 1)
> (define-param r (* 0.2 a))
>
> (set! num-bands 8)
>
> (set! default-material (make dielectric (epsilon epspoly)))
>
> (set! geometry-lattice (make lattice (size 1 1 no-size)
> (basis1 0.5 (/ (sqrt 3) 2))
> (basis2 1 0)))
> (set! geometry (list (make cylinder
> (center 0 0 0) (radius r) (height infinity)
> (material (make dielectric (epsilon epsair))))))
>
> (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)
>
>
>
>
>
> ------------------------------
>
> _______________________________________________
> mpb-discuss mailing list
> [email protected]
> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
>
> End of mpb-discuss Digest, Vol 62, Issue 1
> ******************************************
>
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
