Differences

This shows you the differences between two versions of the page.

--- exercises:2017_ethz_mmm:stm [2017/05/25 12:09] – dpasserone
+++ exercises:2017_ethz_mmm:stm [2017/05/25 12:25] – dpasserone
@@ Line 37: / Line 37: @@
 An empirical potential in teh form of C6/R^6 plus a pauli repulsion
 is added to couple the adsorbate/substrate systems.
@@ Line 71: / Line 70: @@
 </code>
 </note>
+Now go to the STM directory andsubmit the run script
+<code>
+qsub run
+</code>
+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:
+<code>
+./plotorbitals.sc
+</code>
+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:
+<code>
+qsub run_sumbias
+</code>
+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
+<code>
+qsub run_stm
+</code>
+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 p.xyz file that you find in the STM directory
+and execute:
+<code>
+./run_PP
+</code>
+It will take ~ 5 minutes, then you will find a dir containing the AFM simulated image.
+===TASK_2===
+Repeat all the instructions of TASK_1 for the scripts present in the dir TASK_2
+<note warning>
+Be carefulhere we do a spin polarized simulation,
+we have to distinguish the three C atoms of one terminus of the ribbon from the
+three of the opposite terminus calling them C1 and C2.
+When the file p.xyz is created in the STM dir (after running ./pos.sc)
+copy it immediateli to the AFM dir.
+Now, before executing the instructions for the STM dir
+edit the file p.xyz and modify it in such a way that
+the first three C atoms will be labelled as C1
+and the C atoms from 4 to 6 will be labelled as C2
+<code>
+  C1        6.0848407282        7.8280098155       21.6125989354
+  C1        6.0865671686       12.7633436664       21.6071222309
+  C1        6.1020007836       10.2957686990       21.6036624306
+  C2       56.3447906713       10.2958157091       21.6033852713
+  C2       56.3619529363        7.8280149623       21.6128774460
+  C2       56.3601930737       12.7634261117       21.6063533886
+  H         4.9837063610        7.8327959357       21.5912164696
+  H         4.9855872642       12.7623732365       21.5844580428
+</code>
+</note>