diff options
Diffstat (limited to 'figs/bands_third.fig')
-rw-r--r-- | figs/bands_third.fig/Makefile | 23 | ||||
-rw-r--r-- | figs/bands_third.fig/bands_third.gnuplot | 41 | ||||
-rw-r--r-- | figs/bands_third.fig/bands_third.py | 86 |
3 files changed, 150 insertions, 0 deletions
diff --git a/figs/bands_third.fig/Makefile b/figs/bands_third.fig/Makefile new file mode 100644 index 0000000..fc025fa --- /dev/null +++ b/figs/bands_third.fig/Makefile @@ -0,0 +1,23 @@ +PROJECTNAME=bands_third + +all: $(PROJECTNAME).pdf + +$(PROJECTNAME).pdf: + python $(PROJECTNAME).py > $(PROJECTNAME).dat + gnuplot $(PROJECTNAME).gnuplot + pdflatex -jobname $(PROJECTNAME) $(PROJECTNAME)_plot.tex + +install: $(PROJECTNAME).pdf + cp $(PROJECTNAME).pdf $(INSTALLDIR)/ + + +clean-aux: + rm -f $(PROJECTNAME).aux + rm -f $(PROJECTNAME).log + rm -f $(PROJECTNAME).dat + rm -f $(PROJECTNAME)_plot.tex + +clean-pdf: + rm -f $(PROJECTNAME).pdf + +clean: clean-aux clean-pdf diff --git a/figs/bands_third.fig/bands_third.gnuplot b/figs/bands_third.fig/bands_third.gnuplot new file mode 100644 index 0000000..2845518 --- /dev/null +++ b/figs/bands_third.fig/bands_third.gnuplot @@ -0,0 +1,41 @@ +#!/usr/bin/env gnuplot + +datafile="bands_third.dat" +# output +outfile="bands_third_plot.tex" +set output outfile + +# terminal +set term lua tikz size 10,6 standalone tightboundingbox + +# format axes +set key off +unset colorbox +unset border +unset xtics +unset ytics +unset ztics + +# colors +## FF143C (bright red) +## 32CD32 (bright green) + +# angle camera +set view 80,50 + +# style for mesh lines +set style line 1 linetype rgbcolor "#000000" + +# splot mode +set pm3d hidden3d depthorder + +# set colors +set palette defined (0 "#FF143C", 1 "#32CD32") + +# set data range +set zrange [-0.0001:0.0001] + +# plot +splot datafile using 1:2:3:4 with pm3d linestyle 1 , \ + datafile using 5:6:7:8 with pm3d linestyle 1 + diff --git a/figs/bands_third.fig/bands_third.py b/figs/bands_third.fig/bands_third.py new file mode 100644 index 0000000..f4d2290 --- /dev/null +++ b/figs/bands_third.fig/bands_third.py @@ -0,0 +1,86 @@ +#!/usr/bin/env python + +## compute the bands for bilayer graphene (no g4 or delta) +## third regime + +from math import * +import cmath +import numpy +from scipy import optimize + +g1=0.1 +g3=0.33*g1 + +def Omega(x,y): + return(1+2*cmath.exp(-3j/2*x)*cos(sqrt(3)/2*y)) + +# Hamiltonian +def H(x,y): + return(numpy.array(\ + [[0,g1,0,Omega(x,y).conjugate()],\ + [g1,0,Omega(x,y),0],\ + [0,Omega(x,y).conjugate(),0,g3*Omega(x,y)*cmath.exp(3j*x)],\ + [Omega(x,y),0,g3*(Omega(x,y).conjugate())*cmath.exp(-3j*x),0]]\ + )) + +# eigenvalues +def eigsH(x,y): + return(numpy.linalg.eigvals(H(x,y))) + +# resolution +nrpoints=101 +# xrange +xmin,xmax=2*pi/3-0.004,2*pi/3+0.004 +# yrange +ymin,ymax=2*pi/3/sqrt(3)-0.004,2*pi/3/sqrt(3)+0.004 + +# sample points +x=numpy.linspace(xmin,xmax,nrpoints) +y=numpy.linspace(ymin,ymax,nrpoints) +x,y=numpy.meshgrid(x,y) + + +## We want a data point at each Fermi point, which must first be computed + +# A function that vanishes only at the Fermi points +def fermpt(v): + eigs=numpy.sort(numpy.real(eigsH(v[0],v[1]))) + return([eigs[1]-eigs[2],0]) +# compute Fermi points +## approximate starting value +pf_app=numpy.zeros((4,2)) +pf_app[0]=[2*pi/3,2*pi/3/sqrt(3)] +pf_app[1]=[2*pi/3,2/sqrt(3)*acos(1/2-g1*g3/2)] +pf_app[2]=[acos(g1*g3/2-1/2),acos(g1*g3/2-1/2)/sqrt(3)] +pf_app[3]=[4*pi/3-acos(g1*g3/2-1/2),acos(g1*g3/2-1/2)/sqrt(3)] +pf=pf_app +for i in range(0,4): + pf[i]=optimize.root(fermpt,pf_app[i],tol=1.e-10).x + +# Reset the data point closest to each Fermi point to the fermi point +for k in range(0,4): + mini,minj=0,0 + minval=100 + for i in range(0,nrpoints): + for j in range(0,nrpoints): + if(sqrt((x[i,j]-pf[k,0])**2+(y[i,j]-pf[k,1])**2)<minval): + mini,minj=i,j + minval=sqrt((x[i,j]-pf[k,0])**2+(y[i,j]-pf[k,1])**2) + [x[mini,minj],y[mini,minj]]=pf[k] + + +# data points +z=numpy.zeros((4,nrpoints,nrpoints)) +for i in range(0,nrpoints): + for j in range(0,nrpoints): + eigs=numpy.sort(numpy.real(eigsH(x[i,j],y[i,j]))) + for k in range(0,4): + z[k,i,j]=(eigs[k]) + +# output +for i in range(0,nrpoints): + for j in range(0,nrpoints): + for k in [1,2]: + print("%10f %10f %10f %d " % (x[i,j],y[i,j],z[k,i,j],k),end='') + print() + print() |