In [3]:

    

import meep as mp
from meep import mpb
import math


    

In [4]:

    

lattice = mp.Lattice(size=mp.Vector3(5, 5),
                          basis1=mp.Vector3(math.sqrt(3)/2, 0.5),
                          basis2=mp.Vector3(math.sqrt(3)/2, -0.5))


    

In [5]:

    

real = mp.Vector3(1, 0)


    

In [7]:

    

ms = mpb.ModeSolver(geometry_lattice=lattice,
                    )


    

In [8]:

    

real = mp.Vector3(1, 0)
rec = mp.cartesian_to_reciprocal(real, lattice)
bril = ms.first_brillouin_zone(rec)
bril


    

    
Out[8]:





Vector3<0.33012701892219454, 0.33012701892219454, 0.0>

In [9]:

    

real = mp.Vector3(0, 1)
rec = mp.cartesian_to_reciprocal(real, lattice)
bril = ms.first_brillouin_zone(rec)
bril


    

    
Out[9]:





Vector3<0.5000000000000004, 0.49999999999999956, 0.0>

In [11]:

    

rec=mp.Vector3(-1/3, 1/3)
bril = ms.first_brillouin_zone(rec)
bril


    

    
Out[11]:





Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>

In [ ]: