# Add the root path of the pygra library import os import sys sys.path.append(os.environ['PYGRAROOT']) from pygra import geometry from pygra import heterostructures import numpy as np import matplotlib.pyplot as plt g = geometry.chain() h = g.get_hamiltonian() h1 = h.copy() # copy h2 = h.copy() # copy h1.add_swave(.0) # add electron hole symmetry h2.add_swave(.01) # pairing gap of 0.01 ht = heterostructures.create_leads_and_central(h1, h2, h1) # create the junction ht.delta = 1e-6 # imaginary part of the energy es = np.linspace(-.05, .05, 100) # grid of energies ts = [ht.didv(e) for e in es] # calculate transmission plt.plot(es, ts, marker="o") plt.ylim([0, 4.1]) plt.show()
# Add the root path of the pygra library import os ; import sys ; sys.path.append(os.environ['PYGRAROOT']) from pygra import geometry from pygra import heterostructures import numpy as np import matplotlib.pyplot as plt g = geometry.honeycomb_zigzag_ribbon(5) h = g.get_hamiltonian() h.add_kane_mele(0.1) #h.get_bands() #exit() h1 = h.copy() # copy h2 = h.copy() # copy ht = heterostructures.create_leads_and_central(h1,h2,h1) # create the junction ht.delta = 1e-3 # imaginary part of the energy es = np.linspace(-.5,.5,20) # grid of energies ts = [ht.didv(e) for e in es] # calculate transmission plt.plot(es,ts,marker="o") plt.ylim([0,4.1]) plt.show()
# Add the root path of the pygra library import os ; import sys sys.path.append(os.path.dirname(os.path.realpath(__file__))+"/../../../src") from pygra import geometry from pygra import heterostructures import numpy as np g = geometry.chain() h = g.get_hamiltonian() ht = heterostructures.create_leads_and_central(h,h,h) h.get_bands() # get bandstructure es = np.linspace(-1.,1.,50) ts = [ht.landauer(e) for e in es] np.savetxt("TRANSPORT.OUT",np.matrix([es,ts]).T)
# Add the root path of the pygra library import os ; import sys ; sys.path.append(os.environ['PYGRAROOT']) from pygra import geometry from pygra import heterostructures import numpy as np g = geometry.chain() h = g.get_hamiltonian() ht = heterostructures.create_leads_and_central(h,h,h) h.get_bands() # get bandstructure es = np.linspace(-1.,1.,50) ts = [ht.landauer(e) for e in es] np.savetxt("TRANSPORT.OUT",np.matrix([es,ts]).T)