User Tools

Site Tools


This is an old revision of the document!

Simulation of STM and AFM images for a graphene nanoribbons adsorbed on a metallic substrate

connect to hypatia:
ssh -X
ssh -X hypatia
module load python/2.7.12 

go to your scratch directory:

cd /mnt/scratch/your_username

and copy there the tar file of the exercise:

cp /home/cpi/exercise_12.tar ./
tar -xvf exercise_12.tar
cd exercise_12

We consider two possible chemical terminations for a finite size 7-AGNR. In TASK_1 the ribbon is terminated with a C-H2 bonding while in TASK_2 the termination is C-H The additional H atom present at the termini of the ribbon of TASK_1 will suppress the spin polarized edge states that are evident in the ribbon of TASK_2


Have a look to the cp2k input file cp2k.inp used to obtain quickly the optimized geometry of a ribbon adsorbed on a Au substrate. The ribbon is modelled within DFTB (similar to tight binding) while the substrate is modelled via Embedded Atom Model. An empirical potential in teh form of C6/R^6 plus a pauli repulsion is added to couple the adsorbate/substrate systems.

Two geometry fiels are present: and The input needs both of them.

Have a look at the geometry of the system using ASE:

In [1]: from import read

In [2]: from ase.visualize import view

In [3]: s=read("")

In [4]: view(s)

In [5]: exit()
submit the geometry optimization run
qsub run

After completion of the optimization you should extract the final coordinates of the molecule and copy them in the STM directory to compute the KS orbitals and to ocmpute the STM images you can extract the coordinates running the following script:


Now go to the STM directory andsubmit the run script

qsub run

The program will compute the 10 highest and 10 lowest KS orbitals. You can produce a contour plot of each orbital on a plane ~2A above the ribbon running a pyhton script:


I will also show you how to visualize the orbitals with VMD.

To obtain teh stm images you have to combine different KS orbitals (depending on the bias voltage applied) into a single cube file:

qsub run_sumbias

you will then obtain a cube file for each desired bias voltage (see the script run_sumbias)

Now you can compuyte a constant current STM image runnong the script

qsub run_stm

Please note that we are simulating a molecule, we do not include the electrons of the substrate thus we have a disceret spectrum of energies and it is quite likely that for values of the bias voltage that fall in the HOMO-LUMO gap we will obtain an empty image

Now we can simulate for teh same ribbon a AFM image: Go the the AFM directory of TASK_1 copy there the file that you find in the STM directory and execute:

exercises/2017_ethz_mmm/stm.1495714735.txt.gz · Last modified: 2017/05/25 12:18 by dpasserone