This is an old revision of the document!
Simple STM images
In this exercise we will consider different termination of two polyantryl molecules that are an intermediate step for the formation of a long armchair nanoribbon. 10.1021/ja311099k.
We will show how a simple change in the termination (1 vs. 2 Hydrogens) changes the state structure completely.
bsub -n 16
.
Copy, as usual, the files from the directory /cluster/home03/matl/danielep/LECTURE10/EXERCISE_10.2 (and later here on the media manager).
1. Task: Running the job and looking at the orbitals
This time we will not optimize the structure. With an ENERGY run, we run with cp2k the job 1h.1.5.inp and 2h.1.5.inp, meaning that there are here 1.5 units of the original molecule in the gas phase.
The interesting part of the code is in the &PRINT section of &DFT :
&PRINT &STM BIAS -2.0 -1.0 1.0 2.0 TH_TORB S &END STM &MO_CUBES NHOMO 10 NLUMO 10 STRIDE 2 2 2 WRITE_CUBE T &END &V_HARTREE_CUBE &END &MO FILENAME EIG ADD_LAST NUMERIC EIGENVALUES OCCUPATION_NUMBERS &END &END
There will be an output file with the energy levels and their occupation. The last one can be easily found…
- Look with vmd at the cube files corresponding to the most interesting levels (close to Fermi…). Comment on the distribution of the states.</note>
===== 2. Task: Producing a simple STM image =====
The section &STM shown above produces STM images at different bias (feel free to change), meaning, using the Tersoff-Hamann approximation, it integrates all the density of states with energies between Fermi energy and the Bias potential: this energy interval is involved in the tunnel current. The *STM*cube files are 3D maps of the integrated density of states. Imagine that we have a microscope with a feedback that can keep constant current between tip and sample, by changing the height of the tip on the surface. Since the current is proportional to the density of states, we move the tip on isosurfaces of our cubefile. The program stm ** (in the same working dir) allows to extract a 2D map of the height of a given isosurface.
Run the program in the following way: > stm -c 2h*STM*.cube --isovalues 1E-5 > stm.out
The resulting .igor files contain the z profile (in bohr) and may for example be plotted by gnuplot: