User Tools

Site Tools


exercises:2018_ethz_mmm:pmf

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
exercises:2018_ethz_mmm:pmf [2018/05/25 11:09] dpasseroneexercises:2018_ethz_mmm:pmf [2020/08/21 10:15] (current) – external edit 127.0.0.1
Line 77: Line 77:
 timestep 1.0 timestep 1.0
 # (4) NVT Dynamics # (4) NVT Dynamics
-fix temp1 molecule nvt temp 300 300 100 +fix temp1 molecule nvt temp _TEMP_ _TEMP_ 100 
-fix temp2 substrate nvt temp 300 300 100 +fix temp2 substrate nvt temp _TEMP_ _TEMP_ 100 
-velocity                        all create 300 293288+velocity                        all create _TEMP_ 293288
 fix PMF molecule recenter NULL _Y_ NULL fix PMF molecule recenter NULL _Y_ NULL
 compute temp_molecule molecule temp compute temp_molecule molecule temp
Line 86: Line 86:
 thermo_modify flush yes thermo_modify flush yes
 thermo 50 thermo 50
-dump xyz all xyz 100000 mol_sub.xyz+dump xyz all xyz 100000000 mol_sub.xyz
 dump_modify xyz element C C H C N K Cl dump_modify xyz element C C H C N K Cl
 dump coord all dcd 5000 trajectory.dcd dump coord all dcd 5000 trajectory.dcd
-restart 500000 TCB_PMF.restart+restart _NSTEPS_ TCB_PMF.restart
 run _NSTEPS_              run _NSTEPS_             
 </code> </code>
Line 95: Line 95:
  
 <note important> <note important>
-  - Edit the input file to run 50 picoseconds. These are not a lot, but the simulation will last a few minutes per distance.+  - Edit the input file to run 50 picoseconds, to write the restart at the end, to thermalize at the wished temperature 
 +  - **REPLACE _NSTEPS_ WITH THE NUMBER OF STEPS AND _TEMP_ WITH THE DESIRED TEMPERATURE IN THE md_temp.inp FILE**  
 +  - The length of the simulation is small to get converged averages, but the run will thus last a few minutes per distance making the exercise feasible on 16 cores.
   - Now run the chain of simulations   - Now run the chain of simulations
 <code> qsub run_distance_loop </code> <code> qsub run_distance_loop </code>
   - Directories R_<distance> will be created (check with ** ls -ltr ** )   - Directories R_<distance> will be created (check with ** ls -ltr ** )
   - You can enter a certain directory (e.g., R_12.5), check the **log.lammps** for the evolution of a trajectory, and visualize the trajectory with: **vmd mol_sub.xyz trajectory.dcd**   - You can enter a certain directory (e.g., R_12.5), check the **log.lammps** for the evolution of a trajectory, and visualize the trajectory with: **vmd mol_sub.xyz trajectory.dcd**
-  - At the end of everything, you can perform averages and integrals by understanding (and using) the script ** get_pot_mean_force **+  - At the end of everything, you can perform averages and integrals by understanding (and using) the script ** get_pot_mean_force **  FROM WITHIN THE DIRECTORY T_300: 
 +<code> 
 +./get_pot_mean_force 
 +</code> 
 +  - THIS PRODUCES THE POTENTIAL OF MEAN FORCE pot_mean_force for that temperature
   - Repeat the same for other temperatures (e.g.) 10 K and 800 K.   - Repeat the same for other temperatures (e.g.) 10 K and 800 K.
   - Potentials can be visualized by adapting the **gnuplot** script in the exercise_12 directory: pot_mean_force.gnu. Edit and adapt, then:   - Potentials can be visualized by adapting the **gnuplot** script in the exercise_12 directory: pot_mean_force.gnu. Edit and adapt, then:
Line 107: Line 113:
 </code> </code>
 </note> </note>
 +<note important> 
 +Check the already finished runs in REF_300:  
 +<code> 
 +cd REF_300 
 +cp ../get_pot_mean_force . 
 +./get_pot_mean_force 
 +gnuplot 
 +gnuplot > plot "pot_mean_force" w lp  
 +</code> 
 +</note>
  
 <note warning> <note warning>
exercises/2018_ethz_mmm/pmf.1527246592.txt.gz · Last modified: 2020/08/21 10:15 (external edit)