GTFermiSurfaceCut

GTFermiSurfaceCut[Hamiltonian,Fermi energy,list of bands,ndel,plane, area] calculates a cut throught the Fermi surface corresponding to a Hamiltonian and Fermi energy if the Fermi surface contains parts from list of bands. The electronic structure is calculated in a cube at ndel points per spatial dimension. plane defines orientation and position of the cutting plane. area defines a quadratic region in the cutting plane.

DetailsDetails

  • GTFermiSurfaceCut is strongly related to GTFermiSurface. Instead of the whole Fermi surface a cut is created, i.e. a plane is defined and the figure resulting from the cut of the plane with the Fermi surface is presented.
  • The following options can be given:
  • GOVerboseFalsecontrol of put of additional information.
    GOPlotTruedefines the out output of the cuts.
    ContourStyle{}defines the style of the contour lines.
  • See: W. Hergert, M. Geilhufe, Group Theory in Solid State Physics and Photonics. Problem Solving with Mathematica, chapter 9.7.1.

ExamplesExamplesopen allclose all

Basic Examples  (1)Basic Examples  (1)

First load the package:

In[1]:=
Click for copyable input
In[2]:=
Click for copyable input

Read the TB parameterset for Cu.

In[3]:=
Click for copyable input
Out[3]=

Read the predefined Hamiltonian for the fcc structure.

In[4]:=
Click for copyable input

The Hamiltonian for Cu is prepared by insertion of the parameter set.

In[5]:=
Click for copyable input

To see that all works correctly, the band structure is calculated. This is not necessary for the construction of the XCrySDen files.

Bandstructure of Cu

In[6]:=
Click for copyable input
Out[6]=

The exact position of the Fermi energy can be found from the DOS. We calculate the isoenergetic surface to E = 0.6 Ryd for band number 6

Fermi surface

In[7]:=
Click for copyable input
Out[7]=

Now we create a cut throught the Fermi surface. We want to see the cut in the kk-ky-plane. Thus the normal vector will be (0,0,1) amd the shift vector (0,0,0).

In[8]:=
Click for copyable input
In[9]:=
Click for copyable input
Out[9]=

The Brillouin zone boundaries in this plane can be added to the plot.

In[10]:=
Click for copyable input
In[11]:=
Click for copyable input
In[12]:=
Click for copyable input
Out[12]=