Difference between revisions of "AutoMeKin"

From AutoMeKin
Jump to navigation Jump to search
(License)
(Program execution)
Line 30: Line 30:
 
==[[Installation instructions]]==
 
==[[Installation instructions]]==
  
==[[Program execution]]==
+
==Program execution==
 +
To start using any of the scripts of the program, you have to load tsscds/2018 module:
 +
 
 +
<code>module load tsscds/2018</code>
 +
 
 +
To run the low-level calculations use:
 +
 
 +
<code>nohup llcalcs.sh molecule.dat ntasks niter runningtasks >llcalcs.log 2>&1 &</code>
 +
 
 +
where:<br />
 +
<code>molecule</code> is  the name of your molecule<br />
 +
<code>ntasks</code> is the number of tasks<br />
 +
<code>niter</code> is the number of iterations<br />
 +
<code>runningtasks</code> is the number of simultaneous tasks
 +
 
 +
To run the high-level calculations use:
 +
 
 +
<code>nohup hlcalcs.sh molecule.dat runningtasks >hlcalcs.log 2>&1 &</code>
 +
 
 +
For more details, follow the instructions given in the [http://forge.cesga.es/frs/?group_id=77&release_id=21 tutorial]
  
 
==[[References and citations]]==
 
==[[References and citations]]==

Revision as of 16:37, 14 April 2019

Automated Reaction Mechanisms and Kinetics (AutoMeKin)

AutoMeKin has been designed to discover reaction mechanisms and solve the kinetics in an automated fashion, using chemical dynamics simulations. AutoMeKin obtains transition state (TS) guess structures from trajectory simulations performed at very high energies or temperatures. From the obtained TS structures, minima and product fragments are determined following the intrinsic reaction coordinate (IRC). Then, with all the stationary points, the reaction network is constructed. Finally, the kinetics is solved using the Kinetic Monte Carlo (KMC) method. The program is interfaced with MOPAC2016 and Gaussian 09 (G09), but work is in progress to incorporate more electronic structure programs.

Details of the method are described in these two publications: 1 2

Users are encouraged to read the tutorial, which is thought to guide you through the various steps necessary to predict reaction mechanisms and kinetics of unimolecular decompositions. To facilitate the presentation, we consider, as an example, the decomposition of formic acid (FA).

The program has been tested on the following Linux distros: CentOS 7, Red Hat Enterprise Linux and Ubuntu 16.04.3 LTS

Authors

George L. Barnes, Sabine Kopec, Daniel Peláez, Aurelio Rodríguez, Roberto Rodríguez-Fernández, James J. P. Stewart, Pablo Tahoces, Saulo A. Vázquez and
Emilio Martínez-Núñez

Address
Departamento de Química Física, Facultade de Química, Avda. das Ciencias s/n, 15782 Santiago de Compostela, SPAIN
email me

License

A computer program for finding reaction mechanisms and solving the kinetics. Copyright (C) 2018 Emilio Martinez-Nunez, Aurelio Rodriguez, Roberto Rodriguez-Fernandez, Saulo Vazquez-Rodriguez, George L. Barnes and James J. P. Stewart

This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Download source code & tutorial

Installation instructions

Program execution

To start using any of the scripts of the program, you have to load tsscds/2018 module:

module load tsscds/2018

To run the low-level calculations use:

nohup llcalcs.sh molecule.dat ntasks niter runningtasks >llcalcs.log 2>&1 &

where:
molecule is the name of your molecule
ntasks is the number of tasks
niter is the number of iterations
runningtasks is the number of simultaneous tasks

To run the high-level calculations use:

nohup hlcalcs.sh molecule.dat runningtasks >hlcalcs.log 2>&1 &

For more details, follow the instructions given in the tutorial

References and citations

Changelog

News

Use threads=1 in your input file (see the tutorial) to avoid multhreading in MOPAC calculations for much better performance