This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
probe_particle_model [2017/01/02 16:00] krejcio |
probe_particle_model [2017/01/02 16:16] krejcio |
||
---|---|---|---|
Line 4: | Line 4: | ||
New code is written in C/Python and can operate in framework of Lennard-Jones forces as well as electrostatic forces, if necessary. | New code is written in C/Python and can operate in framework of Lennard-Jones forces as well as electrostatic forces, if necessary. | ||
+ | |||
+ | {{:ptcda_df.png|}} | ||
===== Older Fortran Version ===== | ===== Older Fortran Version ===== | ||
Line 61: | Line 63: | ||
Symbol/or/Z-of-element x y z charge-optional | Symbol/or/Z-of-element x y z charge-optional | ||
- | The Lennard-Jones potential needed for calculations can be created by running: | + | The Lennard-Jones (L-J) potential needed for calculations can be created by running: |
python PATH_TO_YOUR_PROBE_PARTICLE_MODEL/generateLJFF.py -i YOUR_INPUT_FILE.xyz | python PATH_TO_YOUR_PROBE_PARTICLE_MODEL/generateLJFF.py -i YOUR_INPUT_FILE.xyz | ||
If charges are also in the input file add "-q" flag to create electrostatic field. | If charges are also in the input file add "-q" flag to create electrostatic field. | ||
Line 71: | Line 73: | ||
- | x, y and z components of Lennard-Jones forces are stored in __LJFF_x.xsf__, __LJFF_y.xsf__ and __LJFF_z.xsf__ files, respectively. These files that can be viewed e.g. via XCrySDen (http://www.xcrysden.org/) or VESTA (http://jp-minerals.org/vesta/en/). | + | x, y and z components of L-J forces are stored in __LJFF_x.xsf__, __LJFF_y.xsf__ and __LJFF_z.xsf__ files, respectively. These files that can be viewed e.g. via XCrySDen (http://www.xcrysden.org/) or VESTA (http://jp-minerals.org/vesta/en/). |
Line 97: | Line 99: | ||
scanStep 0.1 0.1 0.10 # steps of scan (dx, dy, dz) | scanStep 0.1 0.1 0.10 # steps of scan (dx, dy, dz) | ||
Amplitude 1.0 # [Å] oscillation amplitude for conversion Fz->df | Amplitude 1.0 # [Å] oscillation amplitude for conversion Fz->df | ||
+ | |||
+ | If you want to make a scan for different probe, you have to change the probeType in __params.ini__ and to recompute L-J forces. | ||
===== Simulating AFM ===== | ===== Simulating AFM ===== | ||
Line 102: | Line 106: | ||
With having L-J and electrostatic forces made by generating scripts, you can try to run an AFM scan via: | With having L-J and electrostatic forces made by generating scripts, you can try to run an AFM scan via: | ||
python PATH_TO_YOUR_PROBE_PARTICLE_MODEL/relaxed_scan.py | python PATH_TO_YOUR_PROBE_PARTICLE_MODEL/relaxed_scan.py | ||
- | which run the scan with charge Q and lateral stiffness K written in __params.ini__. The result - Fz force acting on the tip -- is saved in __Q?.??K?.??__ directory as an __OutFz.xsf__ file. | + | which run the scan with charge Q and lateral stiffness K written in __params.ini__. Please note, that electrostatic forces are necessary only if Q ≠ 0.0. The result - Fz force acting on the tip -- is saved in __Q?.??K?.??__ directory as an __OutFz.xsf__ file. |
The df results for constant height scans can be plotted by: | The df results for constant height scans can be plotted by: | ||
python PATH_TO_YOUR_PROBE_PARTICLE_MODEL/plot_results.py --df | python PATH_TO_YOUR_PROBE_PARTICLE_MODEL/plot_results.py --df | ||
Line 112: | Line 116: | ||
If a flag - -pos is applied for both commands (scanning & plotting) than xy positions of the relaxing Probe Particle (PP) are shown in __xy_???.png__ as a red dots, while the gray scale on the background maps represent z position of the PP (brighter - higher). | If a flag - -pos is applied for both commands (scanning & plotting) than xy positions of the relaxing Probe Particle (PP) are shown in __xy_???.png__ as a red dots, while the gray scale on the background maps represent z position of the PP (brighter - higher). | ||
- | ===== Test ===== | + | ===== Tests ===== |
Examples of df simulation is already in the downloaded/cloned repository in folder __examples/__. You should try to run these examples before going to your own stuff, in order to see that the code is working on your machine. | Examples of df simulation is already in the downloaded/cloned repository in folder __examples/__. You should try to run these examples before going to your own stuff, in order to see that the code is working on your machine. | ||
Line 120: | Line 124: | ||
===== Scans with different charge (Q), lateral stiffness (K) or oscillation amplitude (A) ===== | ===== Scans with different charge (Q), lateral stiffness (K) or oscillation amplitude (A) ===== | ||
- | A scan with different charge (Q) and/or lateral stiffness (K) than those written in __params.ini__ can be calculated via running: | + | A scan with different charge (Q) and/or lateral stiffness (K) than those written in __params.ini__ can be calculated via running (be aware, that for Q ≠ 0.0, you need to have precalculated electrostatic forces): |
python PATH_TO_YOUR_PROBE_PARTICLE_MODEL/relaxed_scan.py -q (Q) -k (K) | python PATH_TO_YOUR_PROBE_PARTICLE_MODEL/relaxed_scan.py -q (Q) -k (K) |