CAPSChain Assembly and Packing Suite

Run · module 16

Dynamics

Run molecular dynamics directly in CAPS with modern thermostats and barostats, bond constraints and multiple time steps — or write the same run as a deck for LAMMPS or GROMACS.

Open it Top bar›Dynamics

You start with
A structure with a force field, usually minimised or equilibrated.
You get
A trajectory, thermodynamic log and the final state; or input decks for another engine.
Dynamics · CAPS Dynamics with the ensemble, thermostat and barostat settings and live snapshots.

Step by step

  1. Choose the ensemble and time step

    NVE, NVT (default), NPT or NPH; Δt 1 fs. For 2 fs, constrain bonds to hydrogen (SHAKE, RATTLE (as LAMMPS) or LINCS (as GROMACS)) or use r-RESPA.

  2. Set the thermostat and barostat

    Thermostat: Bussi velocity rescaling (default), Langevin (BAOAB) or Nosé–Hoover chain; T 300 K, τ_T 100 fs.

    Barostat: Stochastic cell rescaling (default), Berendsen (early relaxation only) or MTK; P 1 atm, τ_P 1000 fs; coupling isotropic, per axis, only z, only x and y, or full shape.

  3. Check the interactions

    Electrostatics by damped shifted force or particle-mesh Ewald, the cut-off, and the mixing rule (the force field's own by default).

  4. Set the length and the output

    Steps 20 000, a frame every 1000 steps, a checkpoint interval, and Run on: this machine or a host from Settings.

  5. Press Run

    Snapshots and thermodynamics update live; Stop or pause at any time and continue later.

    Prefer another engine? Under Export to other engines pick LAMMPS, GROMACS or both and press Save input deck….

The same without the Studio

Every Studio page calls the same core as the caps command and the caps Python package, so a step you clicked can be repeated in a script.

caps md cell_eq.data -o cell_md.data --steps 20000 --dt 1 --temp 300 \
    --thermostat bussi --barostat crescale --pressure 1 \
    --dump traj.lammpstrj --every 1000 --log thermo.csv

Where to go next

Read more