Abstract:Chemical reactions are often associated with an energy barrier along the reaction pathway which hinders the spontaneity of the reaction. Changing the energy barrier along the reaction pathway allows one to modulate the performance of a reaction. We present a module, Python Algorithms for Searching Transition stAtes (PASTA), to calculate the energy barrier and locate the transition state of a reaction efficiently. The module is written in python and can perform nudged elastic band, climbing image nudged elastic band and automated nudged elastic band calculations. These methods require the knowledge of the potential energy surface (and its gradient along some direction). This module is written such that it works in conjunction with density functional theory (DFT) codes to obtain this information. Presently it is interfaced with three well known DFT packages: Vienna Ab initio Simulation Package (VASP), Quantum Espresso and Spanish Initiative for Electronic Simulations with Thousands of Atoms (SIESTA). This module is easily extendable and can be interfaced with other DFT, force-field or empirical potential based codes. The uniqueness of the module lies in its user-friendliness. For users with limited computing resources, this module will be an effective tool as it allows to perform the calculations image by image. On the other hand, users with plentiful computing resources (such as users in a high performance computing environment) can perform the calculations for large number of images simultaneously. This module gives users complete flexibility, thereby enabling them to perform calculations on large systems making the best use of the available resources. Program summary Program Title: PASTA Program files doi: http://dx.doi.org/10.17632/rv7fdm5gkf.1 Licensing provisions: BSD 3-clause Programming language: Python External routines/libraries: numpy, matplotlib Nature of problem: Most of the reactions have an energy barrier on their reaction pathway. This energy barrier affects the progress of the reaction. Solution method: We implement the NEB, CI-NEB and AutoNEB method to locate transition state and estimate the energy barrier. Our module works with density functional theory codes: VASP, SIESTA and Quantum Espresso presently.

The circumstellar envelope of the C-rich post-AGB star HD 56126

Towards a Topological Classification of Nonadiabaticity in Chemical Reactions

Phase space barriers and dividing surfaces in the absence of critical points of the potential energy: Application to roaming in ozone

Distance measures and evolution of polymer chains in their topological space

Analyzing Complex Reaction Mechanisms Using Path Sampling

Metastability, conformation dynamics, and transition pathways in complex systems

Constraint methods for determining pathways and free energy of activated processes

Topology of Classical Molecular Optimal Control Landscapes in Phase Space

PASTA: Python Algorithms for Searching Transition stAtes

Integrated Reaction Path Processing from Sampled Structure Sequences

Target search of active agents crossing high energy barriers

Transition Pathways in Complex Systems: Reaction Coordinates, Isocommittor Surfaces, and Transition Tubes

From A to B in free energy space

Reactive Vortexes in a Naturally Activated Process: Non-Diffusive Rotational Fluxes at Transition State Uncovered by Persistent Homology

Topological methods for exploring low-density states in biomolecular folding pathways

Toward an Accurate Black-Box Tool for the Kinetics of Gas-Phase Reactions Involving Barrier-less Elementary Steps

Combined Molecular Dynamics and Coordinate Driving Method for Automatically Searching Complicated Reaction Pathways

Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias

Path Topology in Molecular and Materials Sciences

Roaming dynamics in ion-molecule reactions: phase space reaction pathways and geometrical interpretation

A Roadmap for Simulating Chemical Dynamics on a Parametrically Driven Bosonic Quantum Device