CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace: • magnus-lih-hfx • band • rtp_field_xas • xas_tdp • c2h2_bond_energy_2018 • faq
exercises:2017_uzh_cmest:band

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
exercises:2017_uzh_cmest:band [2017/10/17 10:51] – jglanexercises:2017_uzh_cmest:band [2020/08/21 10:15] (current) – external edit 127.0.0.1
Line 1: Line 1:
-======= Getting the band structure of graphene =======+======= Getting the band structure of WO$_3$ Lattice ======= 
 + 
 +In this exercise, you will carry out band structure calculation using K-point sampling for Cubic lattice WO$_3$. The reference band structure you can find in [[http://pubs.acs.org/doi/abs/10.1021/cm3032225|this paper]] 
 + 
 +{{:exercises:2017_uzh_cmest:wo3.jpeg?1200|}}
  
-In this exercise, you will carry out band structure calculation using K-point sampling for Cubic lattice WO3. The reference band structure you can find in [[http://pubs.acs.org/doi/abs/10.1021/cm3032225|this paper]] 
-{{:exercises:2017_uzh_cmest:wo3.jpeg?400|}} 
 To get the band structure for WO3, only a few changes are required compared to the previous example for [[PDOS|calculating the PDOS]]: To get the band structure for WO3, only a few changes are required compared to the previous example for [[PDOS|calculating the PDOS]]:
  
Line 70: Line 72:
                SPECIAL_POINT ???   #R                SPECIAL_POINT ???   #R
                SPECIAL_POINT ???   #M                SPECIAL_POINT ???   #M
-               NPOINTS 5+               NPOINTS ???
             &END             &END
          &END BAND_STRUCTURE          &END BAND_STRUCTURE
Line 117: Line 119:
   * The ''SPECIAL_POINT'' keyword is used to specify the start-, mid- and endpoints of the line. Those points usually denote special points in the reciprocal lattice/unit cell, like the $\Gamma$, $M$ or $K$ point. You can find the definition for these in the appendix section of [[http://www.sciencedirect.com/science/article/pii/S0927025610002697|this paper]]. This keyword can also be specified multiple times, making it possible to directly get the band structure for a complete //path//.   * The ''SPECIAL_POINT'' keyword is used to specify the start-, mid- and endpoints of the line. Those points usually denote special points in the reciprocal lattice/unit cell, like the $\Gamma$, $M$ or $K$ point. You can find the definition for these in the appendix section of [[http://www.sciencedirect.com/science/article/pii/S0927025610002697|this paper]]. This keyword can also be specified multiple times, making it possible to directly get the band structure for a complete //path//.
  
-Now, when you run this input file you will get in addition the the output file, a file named ''graphene.bs'' which will look similar to the following:+Now, when you run this input file you will get in addition the the output file, a file named ''WO3.bs'' which will look similar to the following:
  
 <code> <code>
Line 149: Line 151:
 Your tasks: Your tasks:
  
-  * Lookup the special points for the $\Gamma$, $X$,$M$,$\Gamma$,$R$,$M$ points in the [[http://pubs.acs.org/doi/abs/10.1021/cm3032225|mentioned paper]] (make sure you choose the right lattice). Calculate and plot the band structure for graphene from $\Gamma$ over $M$, $K$ back to $\Gamma$ (you are free to decide whether to use multiple K-Point sets are multiple special points in a single set). Mark the special points. Choose an appropriate number of interpolation points to get a smooth plot.+  * Lookup the special points for the $\Gamma$, $X$,$M$,$R$ points in the [[http://pubs.acs.org/doi/abs/10.1021/cm3032225|mentioned paper]] (make sure you choose the right lattice). Calculate and plot the band structure for WO3 lattice along $\Gamma$, $X$,$M$,$\Gamma$,$R$,$M$ (you are free to decide whether to use multiple K-Point sets are multiple special points in a single set). Mark the special points. Choose an appropriate number of interpolation points to get a smooth plot.
   * Compare your plot with plots from literature. What is different?   * Compare your plot with plots from literature. What is different?
-  * Why do you get 8 orbital energies? Try to change the input to get more unoccupied orbitals.+  * How many orbital energies do you get and why? Try to change the input to get more unoccupied orbitals.
  
-To convert the band structure file to a file which can be loaded directly into MATLAB for example, you can use the script ''cp2k_bs2csv.py'' from below, which when passed a band structure file ''graphene.bs'' as an argument will write files ''graphene.bs-setN.csv'' for each set containing the K-Point coordinates and the energies in one line.+To convert the band structure file to a file which can be plotted directly, you can use the script ''cp2k_bs2csv.py'' from below, which when passed a band structure file ''WO3.bs'' as an argument will write files ''WO3.bs-set-1.csv'' for each set containing the K-Point coordinates and the energies in one line.
  
 +To plot the ''WO3.bs-set-1.csv'' file, you can either load it into $MATLAB$ or use $GNUPLOT$ command line.
 +<code>gnuplot>set yrange [-8:14]
 +gnuplot>plot for [i=4:23] "WO3.bs.set-1.csv" u 0:i w l t ""
 + </code>
 <file python cp2k_bs2csv.py> <file python cp2k_bs2csv.py>
 #!/usr/bin/env python #!/usr/bin/env python
exercises/2017_uzh_cmest/band.1508237481.txt.gz · Last modified: (external edit)