Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aers_v1_0.tar.gz(75333 Kbytes)
Manuscript Title: The VENUS/NWChem Software Package. Tight Coupling between Chemical Dynamics Simulations and Electronic Structure Theory
Authors: Upakarasamy Lourderaj, Rui Sun, Swapnil C. Kohale, George L. Barnes, Wibe A. de Jong, Theresa L. Windus, William L. Hase
Program title: VENUS/NWChem
Catalogue identifier: AERS_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 185(2014)1074
Programming language: Fortran 77 with some C in NWChem, MPI.
Computer: All Linux based workstations and parallel supercomputers.
Operating system: Linux.
Has the code been vectorised or parallelized?: Venus is a sequential code; NWChem can run in parallel.
Keywords: Direct Dynamics, Classical Trajectories, Molecular Simulation.
Classification: 16.8.

Subprograms used:
Cat Id Title Reference
AEGI_v1_0 NWChem CPC 181(2010)1477

Nature of problem:
Direct dynamics simulations play an important role in investigating and understanding atomic-level chemical dynamics information such as atomistic reaction mechanisms, unimolecular and bimolecular rate constants, intramolecular vibrational energy redistribution rates, etc. The ability to couple direct dynamics with electronic structure methods brings a level of fidelity to the simulations that is important for complex systems. However, a tight coupling between two codes that have their own development teams and schedules can be challenging.

Solution method:
The VENUS/NWChem interface is designed to link the general electronic structure program (NWChem) and classical chemical dynamics simulation program (VENUS) to perform direct dynamics simulation in which the trajectories "on the fly" with the potential and its derivatives obtained directly from electronic structure theory. One of the design goals is to build interfaces that require as little interference in NWChem and VENUS as possible so that each of the code developments can continue independently. This is especially important since VENUS is currently a sequential code and NWChem is a parallel code and being able to compute the energies, gradients, and Hessian in parallel is an important aspect of making the software useful to users. In this manuscript, the tight coupling interface between the two codes is described and examples of its use are given. In the classical chemical dynamics simulation an ensemble of trajectories is calculated, and the initial sampling represents the conditions of the reactants for the chemical reaction under investigation. Each trajectory is evaluated by numerically integrating either Hamilton's or Newton's equations of motion. The Schrödinger equation is solved and the energy and energy gradient is calculated in the electronic structure program (NWChem), this information is passed to the classical trajectory program (VENUS) to solve the equations of motion.

Additional comments:
Full documentation is provided in the distribution file. This includes a README file giving the names and brief description of all the files that make up the package and instructions on the installation and execution of the program. Sample input and output data for test run will also be provided.

The software is free to download and use once a signed license agreement has been received. The agreement will be displayed when the program is requested.

Running time:
The running time depends on the size of the chemical system, simulation time, complexity of the ab initio method and number of CPU's. The ab initio method is the most time consuming part of each step in the calculations and scaling for different types of systems and levels of theory are available in M. Valiev, E. J. Bylaska, D. Wang, K. Kowalski, N. Govind, T. P. Straatsma, J. Nieplocha, E. Apra, T. L. Windus, and W. A. deJong, Comp. Phys. Commun. 181 (2010) 1477. Again, there are many factors that affect the running time for the full simulation and it can range from several hours for simulations of a few atoms with DFT and a small basis set running on a single compute node to several days for simulation of tens of heavy atoms with larger basis set running parallel.