Getting started with PyRETIS

After installing PyRETIS, you can run simulations by using a text based input file or by using the PyRETIS library explicitly.

The former approach is perhaps the simplest, but you will first have to learn how to create input files and this is explained in the input description.

The latter option is more involved, but you are then given more freedom in defining, running and interacting with a simulation. In order to make use of the library, you will have to read about the structure of the PyRETIS library in the introduction to the library and in the detailed reference section.

A more extensive overview can be found in the full user guide. Since version 2.4, PyVisA is automatically installed with PyRETIS. Here, we report some examples showing the use of PyRETIS and PyVisA.


PyRETIS examples

In addition to reading the documentation, we have prepared several examples to highlight the usage and capabilities of PyRETIS. The full list of examples can be found in the example section.

../_images/md-nve.png

Molecular dynamics

In this example, we simply run a MD simulation. This is just intended as an example of how one can make use of PyRETIS as a library.

../_images/tis-1d-pot.png

TIS

This example shows how we can run a Transition Interface Sampling calculation and obtain a crossing probability. Here, we consider a 1D potential in which a single particle is moving.

../_images/retis-1d-pot.png

RETIS

This example show how we can run a Replica Exchange Transition Interface Sampling calculation and obtain a crossing probability and a rate constant. Here, we consider a 1D potential in which a single particle is moving.

../_images/sub_wf.png

Subtrajectory moves

This example shows how to use the subtrajectory monte carlo moves Stone Skipping, Web Throwing and Wire Fencing for sampling trajectories of a particle in a 1D well.

../_images/2dpot.png

RETIS 2D

In this example, we perform a simulation of a 2D potential which is constructed such that the selection of the order parameter is not so obvious.

../_images/pyretisrev.png

Extending with C/FORTRAN

This example shows how we can use FORTRAN or C to speed up PyRETIS calculations.

https://openmm.org/images/logo.svg

Using OpenMM

This example demonstrates how we interface between OpenMM and PyRETIS internal code.

../_images/retis-2d-wca.png

RETIS 2D WCA

Here we calculate the rate for the breaking of a bond between two particles in a fluid. We consider two cases: a low barrier and a high barrier case and we implement the force field in C.

../_images/hydrogen.png

Using CP2K

This example demonstrates how we can make use of CP2K for running the dynamics for PyRETIS. Here, we are just studying a toy example - breakage of the bond in hydrogen.

../_images/hydrate.png

Using GROMACS

This example demonstrates how we can make use of GROMACS for running the dynamics for PyRETIS. In this particular example, we study the diffusion of methane in a sI hydrate structure.

../_images/2d-lammps.png

Using LAMMPS

This example demonstrates how we can make use of LAMMPS for running the dynamics for PyRETIS. In this particular example, we revisit the 2D WCA example and use LAMMPS as our molecular dynamics engine.

../_images/contour-400x400.png

PyVisA

In this example, we show the usage of PyVisA. First, the optional requisite, PyQt5, is installed to enable PyVisA GUI. The compressor tool, and the visualization tool usage is then demostrated. A few sample pictures are reported.

../_images/corr_mat.png

Post processing using PyVisA

In this example, we show the usage of the features for post processing availeable in PyVisA. Here we will use the methane hydrate system from the “Using GROMACS” example, where we will add two new collective variables and perform PCA and clustering.