exercises:2015_ethz_mmm:simple_stm
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionLast revisionBoth sides next revision | ||
exercises:2015_ethz_mmm:simple_stm [2015/05/28 21:01] – pshinde | exercises:2015_ethz_mmm:simple_stm [2015/06/01 14:43] – dpasserone | ||
---|---|---|---|
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/ |
</ | </ | ||
===== 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 | + | |
- | | + | |
- | ! add the last iteration with description | + | |
- | | + | |
- | ! printing of eigenvalues of molecular orbitals | + | |
- | | + | |
- | ! also print occupation of the molecular orbitals | + | |
- | | + | |
- | &END | + | |
- | &END | + | |
</ | </ | ||
- | 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 |
- | + | To plot the energy level diagram, copy and paste following lines into the python script | |
- | To plot energy level diagram copy and paste following lines into a file **eld_plot.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 **eld_plot.py** file. Use command python | + | |
< | < | ||
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 88: | Line 76: | ||
y1.append(float((p[0]-Ef)*27.212)) | y1.append(float((p[0]-Ef)*27.212)) | ||
- | yv = np.array(y1) | + | for yv in y1: |
+ | x = [0.2,0.8] | ||
+ | y = [yv, yv] | ||
- | x = [0.2,0.8] | + | if yv <= 0: |
- | y = [yv, yv] | + | plt.plot(x,y,color=' |
+ | else: | ||
+ | plt.plot(x,y, | ||
# add label to y-axis | # add label to y-axis | ||
Line 102: | Line 94: | ||
# remove x-axis label | # remove x-axis label | ||
plt.gca().xaxis.set_major_locator(plt.NullLocator()) | plt.gca().xaxis.set_major_locator(plt.NullLocator()) | ||
- | plt.plot(x, | ||
+ | # show plot | ||
#plt.show() | #plt.show() | ||
- | plt.savefig(' | + | |
+ | # save plot in a eps file | ||
+ | plt.savefig(' | ||
</ | </ | ||
+ | An example of the energy level diagram is as follows, | ||
+ | |||
+ | {{ eld1.png? | ||
<note tip> | <note tip> |
exercises/2015_ethz_mmm/simple_stm.txt · Last modified: 2020/08/21 10:15 by 127.0.0.1