#
# mp_graphene.ipynb
#
from qtdraw import QtDraw
# initialization.
win = QtDraw() # create QtDraw application.
win.mp_set_group("D6h^1") # set space group.
win.set_model("graphene") # set model name.
# draw equivalent sites and bonds.
win.mp_add_site("[1/3,2/3,0]")
win.mp_add_bond("[2/3,1/3,0];[1/3,2/3,0]")
# draw vortex-like M-dipole structure.
lst = win.mp_vector_samb_list("[1/3,2/3,0]", "M")
win.mp_add_vector_samb_modulation("[[M03, 1/4, [1/2,0,0], cos], [M03, 1/2, [0,1/2,0], cos], [M04, -sqrt(3)/4, [1/2,0,0], cos], [M05, -sqrt(3)/4, [1/2,1/2,0], cos], [M06, 1/4, [1/2,1/2,0], cos]] : [3,3,1]", length=0.8)
# set cell and view.
win.set_cell("all")
win.set_view([0,0,1])