User Tools

Site Tools


exercises:2016_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
Last revisionBoth sides next revision
exercises:2016_ethz_mmm:simple_stm [2016/05/20 08:20] – [1. Task: Running the job and looking at the orbitals] pshindeexercises:2016_ethz_mmm:simple_stm [2016/05/30 11:41] dpasserone
Line 48: Line 48:
 </code> </code>
  
-There will a lot of "**EIG**" files, but only the last ("**EIG-1_0.MOLog**") file is of use for the energy level diagram plotting. Remove the other "**EIG**" files by using the command "** rm *EIG-1_0_*.MOLog**". To plot the energy level diagram, copy and paste following lines into the python script **eldplot.py**. The file **energy_ref.dat** contains energy eigenvalues (in a.u.) in one column from the last "**EIG**" file. The Fermi energy (**Ef** [a.u], is the energy of the highest occupied level) must be entered in the **eldplot.py** script. Use the command **python eldplot.py** to get the energy level diagram as a png image. (use display to visualize it). Identify the occupied and unoccupied energy levels and name them. Feel free to change the png image names.+There will a lot of "**EIG**" files, but only the last ("**EIG-1_0.MOLog**") file is of use for the energy level diagram plotting. Remove the other "**EIG**" files by using the command "** rm *EIG-1_0_*.MOLog**". To plot the energy level diagram, extract the energies from the eigenvalues file, and copy them (single column) into a file with the same shape as the provided example **energy_ref.dat**. Copy and paste following lines into the python script **eldplot.py**. The file **energy_ref.dat** contains the energy eigenvalues (in a.u.) in one column from the last "**EIG**" file (you can also use two names for the two molecules). The Fermi energy (**Ef** [a.u], is the energy of the highest occupied level) must be entered in the **eldplot.py** script. Use the command **python eldplot.py** to get the energy level diagram as a png image. (use display command to visualize it). Identify the occupied and unoccupied energy levels and name them. Feel free to change the png image name.
 <code> <code>
 import matplotlib.pyplot as plt import matplotlib.pyplot as plt
Line 59: Line 59:
 # Open file # Open file
  
-f = open('energy.dat', 'r')+f = open('energy_ref.dat', 'r')
  
 lines = f.readlines() lines = f.readlines()
exercises/2016_ethz_mmm/simple_stm.txt · Last modified: 2020/08/21 10:15 by 127.0.0.1