Documentation

The EON program can run the following methods or tasks to explore the configuration space of molecular systems and to accelerate the simulation of their dynamics over long times. To run EON, a config.ini file must be generated using the options specified in the documentation. Each section header is denoted by square brackets, and is followed by the key/value pairs. For example, to set the job_type key of the [Main] section to the value process_search, your config.ini file would include the lines:

[Main]
job_type = process_search

There are specific options for each method, and a set of general options which are shared between methods. Examples of these general options include specificion of the interatomic potential and the parameters for doing structural optimizations and comparisons.

EON is designed to run in serial on one computer or in parallel using a communicator to send jobs from a server to clients and receive the results back. Several parallelization schemes have been implemented, including local communication through files, cluster communication via a queuing system or via mpi. Methods run in parallel are broken up by the eon server into tasks which are run by client program. The server then compiles the information sent back by the clients in a way that can be used by the sampling or dynamics methods.

Jobs

Adaptive kinetic Monte Carlo: uses saddle searches to find possible reactive events and KMC to determine the state-to-state kinetics.

  • Recycling: reactive events (processes) between states to speed up AKMC.

  • Kinetic Database: of reactive events so they can be reused between AKMC runs.

  • Coarse Graining: to identify and escape groups of states connected by low barriers.

Parallel Replica Dynamics: uses a set of trajectories to accerate the rate of escape from states.

Molecular Dynamics: standard classical dynamics.

  • Hyperdynamics: uses a bias potential to accelerate transitions between states.

Basin Hopping: lowers the barrier between states to accelerate Monte Carlo sampling.

Process Search: finds possible escape mechanisms from a state.

Saddle Search: finds a nearby saddle point.

  • Dimer: method to find the lowest curvature mode.

  • Lanczos: method to find the lowest curvature mode.

Minimization: optimize the geometry of a structure.

Nudged Elastic Band: Locate minimum energy pathways using nudged elastic band calculation.

General Options

Main: options which not specific to a single job type.

Communicator: methods by which the code can be run in parallel.

Potential: the interatomic potentials available, as well as interfaces to LAMMPS and DFT codes.

Optimizer: methods for optimizing atomic structures.

Hessian: options for controlling the calculation of Hessian matricies.

Prefactor: options for controlling the calculation of harmonic transition state theory prefactors.

Structure Comparison: options for determining if two structures are the same.

Paths: directories in which the calculations will be run.

Debug: options for printing more information to help debug calculations.