Dynamics

Newtonian and stochastic dynamics backends for atooms.

Quick start

Run a molecular dynamics simulation of a Lennard-Jones system from an existing xyz file

from atooms.trajectory import Trajectory
from atooms.simulation import Simulation
from atooms.dynamics.netwonian import VelocityVerlet

# Start from the last frame of input.xyz
trajectory = Trajectory('input.xyz')
system = trajectory[-1]
system.interaction = Interaction('lennard_jones')
backend = VelocityVerlet(system, timestep=0.002)
sim = Simulation(backend, steps=200)
sim.run()

Do the same via the API, storing configurations in output.xyz,

from atooms.dynamics.api import md

md('input.xyz', 'output.xyz',
   method='velocity-verlet', model='lennard_jones',
   dt=0.002, nsteps=200, config_number=20)

or from the command line

md.py --method velocity-verlet -n 200 --dt 0.002 --config-number 20 input.xyz output.xyz

Features

Integration algorithms:

• Netwonian dynamics
  • velocity-Verlet
  • Nosé-Poincaré
  • Nosé-Hoover and Nosé-Hoover chains
  • Berendsen thermostat/barostat
  • event-driven
• Stochastic dynamics
  • Langevin
  • overdamped Langevin

Documentation

Check out the tutorial for more examples and the public API for full details.

Installation

From pypi

pip install atooms-dynamics

You can clone the code repository and install from source

git clone https://framagit.org/atooms/dynamics.git
cd dynamics
pip install .

Contributing

Contributions to the project are welcome. If you wish to contribute, check out these guidelines.

Authors

Daniele Coslovich: https://www.units.it/daniele.coslovich/