GTPlotStateWeights

GTPlotStateWeights[cluster,distance,basis,wave function,scale]
demonstrates at which atoms of a cluster, the wave function (a solution of a real space tight-binding Hamiltonian) has the heighest weights.

DetailsDetails

  • The real space tight-binding Hamiltionian of a cluster of ( - number of basis atoms, - number of atoms of -th set in the cluster) atoms has the size ( - number of orbitals of -th sort). The Hamiltonian has therefore eigenvalues and the corresponding number of eigenvectors. The weight of the orbitals of an atom and sort is defined as:
  • ( - corresponding components in the eigenvector).
  • The wave function weights at the atoms are depicted by spheres of different size.
  • clustercluster of atoms, used to construct the Hamiltonian
    distanceminimal distance of atoms, used to construct bonds for the plot
    basisinformation about the basis used to construct the Hamiltonian
    wave functionan eigenvector of the Hamiltonian
    scalescaling factor, defines the radius of a sphere of maximal weight
  • The following options can be given:
  • GOColorScheme"ElementData"Defines the colors used to plot the atoms in the structure graph
    GOPlotTruePlots the graph
  • See: W. Hergert, M. Geilhufe, Group Theory in Solid State Physics and Photonics. Problem Solving with Mathematica

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

For comparison the DOS is calculated in k-space

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

To demonstrate GTPlotStateWeights a small cluster is used.

The Hamiltonian in real space has the dimension 370x370. The corresponding cluster consists of 37 Ga and 37 As atoms. This real space Hamiltonian is also already prepared:

In[4]:=
Click for copyable input

The cluster is also already calculated:

In[5]:=
Click for copyable input

The density of states can be calculated.

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

Due to the small cluster the real space DOS deviates from the k-space results. The comparison of both DOS calculations shows, that due to the surface of the cluster surface states appear in the band gap in the real space calculation.

GTFindStateNumbers can be used to find the numbers of eigenstates in an energy interval. To investigate the weights of the states on the different atoms, the eigenvectors have to be calculated too.

In[7]:=
Click for copyable input
In[8]:=
Click for copyable input
Out[8]//TableForm=
In[9]:=
Click for copyable input
In[10]:=
Click for copyable input
Out[10]=

The plot demonstrates the state number 335, lying in the gap, is located on the surface of the cluster.