ReaxFF acceleration & parametrization Tokyo 2018
ReaxFF acceleration & parametrization London 2016
The materials below are based on a two-day hands-on workshop in London In September 2016. Input files, templates and example scripts are available in the whole NSCCS ReaxFF workshop package. Please see also the notes on scripting on Windows and MacOS at the bottom of this page.
The ReaxFF parametrization tutorial shows how to build a very simple force field for C/H/Br with a (too) simplistic training set. It begins with how to get a good starting guess, in this case finding a C/H/Cl force field with reasonable C-Cl bonding behavior. Then a training set is built from running ADF optimizations and linear transits, using both new and available ADF command line (python) scripting tools (the entire workshop package is available). After running a few MCFF runs to optimize the new force field, reactive MD can be run on radical polymerization. Finally some hints and examples on how to check overfitting are described. This is a very minimalistic starting point; see also the MCFF introduction and overview talk by Adri van Duin for more tips and tricks.
The fbMC tutorial shows how to build a large graphene supercell and compares how defected graphene heals more quickly with the force bias Monte Carlo acceleration technique than with regular ReaxFF molecular dynamics.
A more advanced, command-line based GCMC tutorial shows how a thermodynamic equilibrium under different pressure and temperature can be established with ReaxFF. The Monte Carlo moves are inserting, deleting and moving oxygen atoms from an O2 reservoir, followed by quick ReaxFF minimization. A python scripts analyzes the Pd/O ratio for a nanoparticle as a function of MC steps, establishing equilibrium within a few thousand moves. No GUI support is available yet, but with some extra effort the structures under different p, T conditions can be visualized in the GUI.
The recommend environment for running these tutorial examples is Linux, but you can also script on MacOS, and -to a lesser extent- on Windows.
On any operating system make sure you set the proper ADF environment variables (source adfrc.sh) and use ADF’s startpython to avoid conflicts (see also the scripting manual):
With PLAMS you can automate workflows and manage jobs. In the scripting tutorial you’ll learn how you can use it to access properties on the binary files (‘KF files’) through the KFReader functionality.
Scripting on MacOS – running ADF in a terminal
If you want to run on MacOS from the command line, open a terminal and type a dot. Then open a GUI and go to Help: Show ADFHOME in Finder in the menu, drag that to a terminal and press enter (see video)
Scripting on Windows – ADF in command-line mode
If you want to run these tutorials on Windows (as was done during the London workshop), make sure you always run in the command line mode:
- double click adf_command_line.bat in your ADF installation directory (e.g. C:/ADF2017.106)
- type sh to go to a shell
Note that this shell is even more limited than the MacOS terminal.
Again: whenever you run python scripts, use startpython to use the ADF python stack.
Scripting in AMS2018 and beyond
In the Amsterdam Modeling Suite 2018, you can just use Help -> Terminal for MacOS or Help -> Command-line on Windows to open up a terminal with the correct environment variables. On Windows you should still type
to go to a basic shell.
Please let us know at email@example.com if you have any questions or problems about running these tutorials, scripting with ADF, or anything else.