The ACEMD3 input file may contain any of the commands listed below.
Each command must start on a new line and be followed by its arguments. Commands are case-insensitive, but their arguments may be case-sensitive.
Blank lines are permitted. Lines beginning with a “#” are treated as comments and ignored.
The filename of an Amber PRMTOP file
The filename of a CHARMM PAR file. May be specified more than once.
The filename of a CHARMM PSF file
Mandatory initial system geometry in PDB format.
Optional initial system geometry in NAMD BINCOOR format. If specified, overrides “coordinates”
Optional initial system velocity field in NAMD BINCOOR format. If specified, overrides field generated by “temperature”
Optional initial system velocity field in NAMD BINCOOR format. If specified, overrides field generated by “temperature” and “velocities”
Filename of a NAMD XSC format file giving the periodic cell dimensions. Overrides “celldimension” and any dimensions in the “coordinates” PDB
The temperature in Kelvin for which to generate an initial velocity field.
The dimensions of the unit cell in Angstrom. Note that the unit cell must be cuboid. Overrides any dimension given in the “coordinates” PDB.
ACEMD3 will always output
output.xsc at the end of a simulation. Additionally, it will
periodically write out the coordinate state to an XTC format trajectory
file. Note that DCD format is no longer supported.
The name of the trajectory file to write to.
The number of steps between outputting coordinates to the trajectory file
Use particle-mesh Ewald summation for long-range electrostatics.
Real-space cutoff in Angstroms for electrostatics and van der Waals.
Apply switching function to the van der Waals potential.
Distance in Angstrom at which to begin applying the switching function.
ACEMD3 continues to use a Langevin thermostat. Old commands with the “langevin” prefix will be converted to their new equivalents and a warning issued.
Enable thermostatic control
Target temperature (K) for thermostatic control
Damping constant for the Langevin thermostat in ps^-1
The barostat algorithm in ACEMD3 has changed from Berendsen to a non-virial Monte Carlo scheme. Old commands with the “berendsen” prefix will be converted to their new equivalents and a warning issued.
Enable pressure control.
The target pressure in bar.
Allow X, Y and Z unit cell dimensions to vary independently
Constrain the X,Y dimensions of the unit cell. Allow Z to vary independently.
Constrain the X:Y radio of the unit cell dimensions. Allow Z to vary independently.
The filename of a PLUMED2 input file. If set, enables PLUMED metadynamics
For more details on the use of these new commands please read Applying Restraints
Specify a restraining potential centred on specific atom centres.
Specify a restraining potential centred on a group of atoms.
The integration timestep in femtoseconds. Note that if timestep >= 2.0 then constraints will be applied to all bonded terms including an H. Additionally, if timestep >= 4.0 hydrogen mass repartitioning is performed, increasing the mass of H to 4au.
If set, attempt to restart a simulation from checkpoint files. Restart
files are always written every
trajectoryFreq steps irrespective of this
The number of energy minimization steps to perform before commencing dynamics.
The length of simulation to run. May be specified as a number of steps or as a time if one of the suffices “us”, “ns”, “ps”, “fs” is used.