Monte Carlo simulations for the estimation of molecule pair interatcion
In this exercise you will perform a MC simulation for different coverages of “sumanene” mlecules adsorbed on a Ag(111) substrate hypothesizing different posisble values for the nearest neighbor energy in the molecule-molecule interaction
you@eulerX ~$ wget http://www.cp2k.org/_media/exercises:2017_ethz_mmm:exercise_3.zip you@eulerX ~$ unzip exercises:2017_ethz_mmm:exercise_3.zip you@eulerX ~$ cd exercise_3
In an experiment performed at Empa, sumanene molecules were adsorbed on a Ag(111) surface. It was found that at very low coverage (0.02) 30% of the molecules were weakly bound into dimers. [ http://dx.doi.org/10.1021/ja504126z J. Am. Chem. Soc. 2014, 136, 13666−13671]
While in execution, the code will show you snapshots of the positions of the molecules on the lattice on the left panel.
In the central panel average values for the number of isolated molecules, the number of dimers and the number of clusters is plotted. On the right panel the average energy of the system is plotted. Intermediate average values are also written in a output file data_de_T.out where de and T are the input values of the dimer energy and of the Temperature also a final snapshot of the graphic panels is saved as png image
coverage 0.02 number of cycles 200 binding energy in eV 0.0 Temperature in K 200
Estimate the dimer binding energy as DE=kT * ln(n0/nexp) where k is Boltzmann's constant T is the simulation (and experiment) temperature n0 is the concentration of dimers in the case of zero interaction nexp is the concentration of dimers found in the experiment To compute the concentration consider that at coverage 0.02, in the simulation, the total number of molecules is 50
TASK3 Repeat the simulation with coverage 0.1 and DE=-0.02 and DE=-0.1 Describe what you obtain. Now try coverage=0.1 T=400 DE=-0.1 Comment the result
TASK4 have a look at the pyhton code, identify the main MC steps in the MAIN part of the code
clusters_plot= for i in range(nouter): and also the section #### DECIDE whether to accept or not the move
What do you think it would happen if you replace the condition
with the condition
which finds out all teh molecules that are connected to a given one
which finds out all the molecules that are 1st neighbors to a given one.
The function “allconnected” is quite intuitive and inefficient. Can you imagine roughly a more efficient function to perform the same task?