I need to generate a set of trajectory of the rotation of one dihedral in a polymer. AMBERTOOLS14 has a command “dihedralscan” to do that, but it only works for dihedrals of proteins. It’s not allowed to define your own dihedral. So I chose SMD to do this.
It’s not complicated to set up the system. The only confusing part is the value of coordinate of torsion is rad not degree.
The input file of NPT part is listed below in case you might need. The initial torsion was about -1rad. To rotate 360degree. I set the reaction coordinate change from the initial value to 5.3. The diehdral is on residue 10 ( atom 100, 196, 155, 175). The rest parts are restrained with 500kcal/mol/angstrom^2.