CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace: • gapw • ls_scf • simple_stm • bs • basic_electronic_structure • surface_au • basic_electronic_structure • tio2_gap • hybrid • surface_cu
exercises:2015_ethz_mmm:simple_stm

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:2015_ethz_mmm:simple_stm [2015/05/29 10:47] – pshindeexercises:2015_ethz_mmm:simple_stm [2020/08/21 10:15] (current) – external edit 127.0.0.1
Line 12: Line 12:
 Download, as usual, the files from the media manager: {{exercise-10.2.tar.gz|}}. The 1h.1.5.inp file is **commented**). Download, as usual, the files from the media manager: {{exercise-10.2.tar.gz|}}. The 1h.1.5.inp file is **commented**).
  
-Please use command **module load cp2k/trunk.2.5.13191** to do this particular exercise.+Please use command **module load cp2k/2.5.14075** , as well as **module load python/2.7.9** to do this particular exercise.
 </note> </note>
 ===== 1. Task: Running the job and looking at the orbitals ===== ===== 1. Task: Running the job and looking at the orbitals =====
Line 45: Line 45:
     &END     &END
  
-      &MO +
-! molecular orbitals eigenvalues and molecular ocuupations after each scf +
-       FILENAME EIG +
-! add the last iteration with description +
-       ADD_LAST NUMERIC +
-! printing of eigenvalues of molecular orbitals +
-       EIGENVALUES +
-! also print occupation of the molecular orbitals +
-       OCCUPATION_NUMBERS +
-      &END +
-    &END+
 </code> </code>
  
-There will be an output file with the energy levels and their occupation. The last one can be easily found... +There will be an output file with the unoccupied energy levels and the last "**EIG**" file with occupied energy levels. 
- +To plot the energy level diagram, copy and paste following lines into the python script **eldplot.py**. The file **energy.dat** contains energy eigenvalues (in a.u.) in one column from the output file and from the EIG file. The Fermi energy (**Ef** in a.u., is the energy of the highest occupied value) must be entered in the **eldplot.py** script. Use the command **python eldplot.py** to get the energy level diagram as a postscript image. (use gs to visualize it). Identify the occupied and unoccupied energy levels and name them. Feel free to change the png image names.
-To plot the energy level diagram, copy and paste following lines into the python script **eldplot.py**. The file **energy.dat** contains energy eigenvalues in a column from the output file at the previous step. The Fermi energy (Ef) must be entered in the **eldplot.py** script. Use command **python eldplot.py** to get the energy level diagram as a png image. Identify the occupied and unoccupied energy levels and name them.+
 <code> <code>
 import matplotlib.pyplot as plt import matplotlib.pyplot as plt
Line 81: Line 70:
  
 # Fermi energy # Fermi energy
- 
 Ef=0.0 Ef=0.0
  
Line 110: Line 98:
 #plt.show() #plt.show()
  
-# save plot in a png file +# save plot in a eps file 
-plt.savefig('ELD.png')+plt.savefig('ELD.eps')
  
 </code> </code>
 An example of the energy level diagram is as follows, An example of the energy level diagram is as follows,
  
-{{ eld.png?direct&700 |}}+{{ eld1.png?direct&700 |}}
  
 <note tip> <note tip>
exercises/2015_ethz_mmm/simple_stm.1432896471.txt.gz · Last modified: (external edit)