Spin texture calculation (Oxygen version)

 

by Benedikt Maurer for exciting oxygen

Purpose: In this tutorial you will learn how to calculate the spin texture with exciting for GaAs.

Fold Table of Contents 0. Define relevant shell variables 1. Theoretical background 2. Spin texture of GaAs Literature

---

0. Define relevant shell variables

Read the following tutorials before starting with the next sections!

• How to set environment variables for tutorials scripts.
• Spin-orbit coupling.

From now on the symbol $ will indicate the shell prompt.

---

1. Theoretical background

An interesting property of the electronic structure with spin-orbit coupling (SOC) is the so-called spin texture. It can be used, e.g., to differentiate between the Rashba and the Dresselhaus effect (see Kepenekian2015).

The spin texture s(n,k) is defined as the expectation value of a vector of Pauli matrices on a two dimensional k-grid. At k for the state n it is defined as:

(1)

s(n,k)=⟨n,k|σ|n,k⟩

---

2. Spin texture of GaAs

Before you start, create a working directory for this tutorial and move into it:

$ mkdir spin-texture-tutorial && cd spin-texture-tutorial

To calculate the spin texture for GaAs, create the file input.xml with your editor of choice and copy and paste the following lines into it:

<input>
   <title>Ga As spin texture</title>
 
   <structure speciespath="$EXCITINGROOT/species">
      <crystal scale="10.683" >
         <basevect>  0.0     0.5     0.5 </basevect>
         <basevect>  0.5     0.0     0.5 </basevect>
         <basevect>  0.5     0.5     0.0 </basevect>
      </crystal>
 
      <species speciesfile='Ga.xml' rmt='2.0'>
         <atom coord="0.00 0.00 0.00"/>
      </species>
 
      <species speciesfile='As.xml' rmt='2.0'>
         <atom coord="0.25 0.25 0.25"/>
      </species>
   </structure>
 
   <groundstate
      do="fromscratch"
      xctype="GGA_PBE"
      rgkmax="6"
      gmaxvr="15"
      ngridk="4 4 4"
      swidth="1d-6"
      epsengy="1d-6"
      maxscl="40"
      nempty="15">
 
      <spin spinorb="true"/>
   </groundstate>
 
   <properties>
      <spintext bands="23 30">
         <plot2d>
            <parallelogram grid = "12 12">
               <origin coord = "-0.042506  -0.042506  -0.085013"/>
               <point  coord = "-0.042506   0.050235   0.007728"/>
               <point  coord = " 0.050235  -0.042506   0.007728"/>
            </parallelogram>
         </plot2d>
      </spintext>
   </properties>
</input>

Set up the correct species path and run exciting to perform the groundstate calculation:

$ SETUP-excitingroot.sh
$ time exciting_smp &

Note that spin-orbit coupling (SOC) is enabled spinorb="true". Spin texture calculations are only allowed for ground states with enabled SOC.
Note also the spintext element in the properties element:

• The bands attribute defines the band range that is considered. In this case we are interested in the bands around the band gap. You can extract the band numbers from the EIGVAL.OUT.
• The parallelogram element defines the k-grid for the spin texture calculation. The k-points are homogeneously distributed on a parallelogram in the reciprocal space. In the example of this tutorial, we consider a square with sidelength 0.1 [a.u.] (a.u. = Bohr^-1) in the (k_x, k_y) plane. The centre of this square is the Γ point and the borders are parallel to the Cartesian k_x and k_y axes.
• grid defines the number of k-points per direction.
• The parallelogram is spanned between origin and the points defined by point.

Be aware of that coord is given in reciprocal coordinates.

When the calculation is done, you can plot the spin texture for the first conduction band using the script PLOT-spintext.py (for a detailed description of the script arguments see The python script "PLOT-spintext.py") as follows:

$ PLOT-spintext.py -b 29 -c energy

• -b defines the band that is considered for the plot which is in this case the first conduction band. You can extract the numbers of the states from EIGVAL.OUT.
• -c specifies if the energy or the z component of the spin orientation is considered for the contour plot. The choices are "energy" or "spin_z".

Note that the script is designed for rectangular spin textures!

The resulting PNG file PLOT-spintext.png can be should look like the following:

In the plot, k₁ and k₂ correspond to the k-vectors (k₁,0,0)^T and (0,k₂,0)^T, respectively. The color map represents the energy of the chosen state. Note that the unit of k₁ and k₂ are given in atomic units (Bohr^-1).

---

Literature

• Kepenekian2015: The form of the spin texture and more information on Rashba and Dresselhaus effects is given in: Kepenekian, M., Robles, R., Katan, C., Sapori, D., Pedesseau, L., & Even, J. (2015). Rashba and Dresselhaus Effects in Hybrid Organic-Inorganic Perovskites: From Basics to Devices. ACS Nano, 9(12), 11557–11567. https://doi.org/10.1021/acsnano.5b04409

 

 

http://exciting.wikidot.com/oxygen-spintexture