def get_zigzag_bilayer(n=10):
"""Get hamiltonian of a zigzag bilayer"""
import geometry
g = geometry.honeycomb_zigzag_ribbon(40)
drs = [np.array([0., 0., -1.]), np.array([0., 1., 1.])]
g = geometry.apilate(g, drs=drs)
h = g.get_hamiltonian()
return h
def get_zigzag_bilayer(n=10): """Get hamiltonian of a zigzag bilayer""" import geometry g = geometry.honeycomb_zigzag_ribbon(40) drs = [np.array([0.,0.,-1.]),np.array([0.,1.,1.])] g = geometry.apilate(g,drs=drs) h = g.get_hamiltonian() return h
def get_armchair_bilayer(n=10):
"""Get hamiltonian of an armchair bilayer"""
import geometry
g = geometry.honeycomb_armchair_ribbon(30)
drs = [np.array([0., 0., -1.]), np.array([np.sqrt(3.) / 2., 0., 1.])]
g = geometry.apilate(g, drs=drs)
h = g.get_hamiltonian()
return h
def get_armchair_bilayer(n=10): """Get hamiltonian of an armchair bilayer""" import geometry g = geometry.honeycomb_armchair_ribbon(30) drs = [np.array([0.,0.,-1.]),np.array([np.sqrt(3.)/2.,0.,1.])] g = geometry.apilate(g,drs=drs) h = g.get_hamiltonian() return h