annotate nvt.xml @ 3:6944c73fe6f9 draft default tip

Uploaded
author fmercuri
date Mon, 28 Jun 2021 10:24:56 +0000
parents c549135d400b
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
2
c549135d400b Uploaded
fmercuri
parents: 0
diff changeset
1 <tool id="gmx_nvt" name="GROMACS NVT equilibration daimoners' version" version="@VERSION@">
0
666dc804011d Uploaded
fmercuri
parents:
diff changeset
2 <description>- constant-temperature and -volume equilibration of a system</description>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
3 <macros>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
4 <import>macros.xml</import>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
5 </macros>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
6
666dc804011d Uploaded
fmercuri
parents:
diff changeset
7 <expand macro="requirements" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
8
666dc804011d Uploaded
fmercuri
parents:
diff changeset
9 <command detect_errors="exit_code"><![CDATA[
666dc804011d Uploaded
fmercuri
parents:
diff changeset
10 #if $mdp.mdpfile == "custom":
666dc804011d Uploaded
fmercuri
parents:
diff changeset
11 ln -s '$mdp.mdp_input' ./nvt.mdp &&
666dc804011d Uploaded
fmercuri
parents:
diff changeset
12 #end if
666dc804011d Uploaded
fmercuri
parents:
diff changeset
13 #if $mdp.mdpfile == "default":
666dc804011d Uploaded
fmercuri
parents:
diff changeset
14 ln -s '$nvt' ./nvt.mdp &&
666dc804011d Uploaded
fmercuri
parents:
diff changeset
15 #end if
666dc804011d Uploaded
fmercuri
parents:
diff changeset
16
666dc804011d Uploaded
fmercuri
parents:
diff changeset
17 ln -s '$gro_input' ./em.gro &&
666dc804011d Uploaded
fmercuri
parents:
diff changeset
18 ln -s '$top_input' ./top_input.top &&
666dc804011d Uploaded
fmercuri
parents:
diff changeset
19 ln -s '$itp_input' ./posres.itp &&
666dc804011d Uploaded
fmercuri
parents:
diff changeset
20
3
6944c73fe6f9 Uploaded
fmercuri
parents: 2
diff changeset
21 gmx grompp -f ./nvt.mdp -c ./em.gro -r ./em.gro -p ./top_input.top -o nvt.tpr -maxwarn 4 &>> verbose.txt &&
0
666dc804011d Uploaded
fmercuri
parents:
diff changeset
22 gmx mdrun -deffnm nvt &>> verbose.txt
666dc804011d Uploaded
fmercuri
parents:
diff changeset
23
666dc804011d Uploaded
fmercuri
parents:
diff changeset
24 #if $str == 'pdb' or $str == 'both'
666dc804011d Uploaded
fmercuri
parents:
diff changeset
25 && gmx editconf -f nvt.gro -o nvt.pdb &>> verbose.txt
666dc804011d Uploaded
fmercuri
parents:
diff changeset
26 #end if
666dc804011d Uploaded
fmercuri
parents:
diff changeset
27 ]]></command>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
28 <configfiles>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
29 <!-- .mdp file for the gromacs simulation -->
666dc804011d Uploaded
fmercuri
parents:
diff changeset
30 <configfile name="nvt">
666dc804011d Uploaded
fmercuri
parents:
diff changeset
31 #if $mdp.mdpfile == 'default':
666dc804011d Uploaded
fmercuri
parents:
diff changeset
32 title = NVT equilibration
666dc804011d Uploaded
fmercuri
parents:
diff changeset
33 define = -DPOSRES ; position restrain the protein
666dc804011d Uploaded
fmercuri
parents:
diff changeset
34 ; Run parameters
666dc804011d Uploaded
fmercuri
parents:
diff changeset
35 integrator = $mdp.integrator ; leap-frog integrator
666dc804011d Uploaded
fmercuri
parents:
diff changeset
36 nsteps = $mdp.md_steps ; 2 * 50000 = 100 ps
666dc804011d Uploaded
fmercuri
parents:
diff changeset
37 dt = $mdp.step_length ; 2 fs
666dc804011d Uploaded
fmercuri
parents:
diff changeset
38 ; Output control
666dc804011d Uploaded
fmercuri
parents:
diff changeset
39 nstxout = $mdp.write_freq ; save coordinates every n ps
666dc804011d Uploaded
fmercuri
parents:
diff changeset
40 nstvout = $mdp.write_freq ; save velocities every n ps
666dc804011d Uploaded
fmercuri
parents:
diff changeset
41 nstenergy = $mdp.write_freq ; save energies every n ps
666dc804011d Uploaded
fmercuri
parents:
diff changeset
42 nstlog = $mdp.write_freq ; update log file every n ps
666dc804011d Uploaded
fmercuri
parents:
diff changeset
43 nstxout-compressed = $mdp.write_freq ; save compressed coordinates every 10.0 ps
666dc804011d Uploaded
fmercuri
parents:
diff changeset
44 ; Bond parameters
666dc804011d Uploaded
fmercuri
parents:
diff changeset
45 continuation = no ; first dynamics run
666dc804011d Uploaded
fmercuri
parents:
diff changeset
46 constraint_algorithm = lincs ; holonomic constraints
666dc804011d Uploaded
fmercuri
parents:
diff changeset
47 constraints = $mdp.constraints ; all bonds (even heavy atom-H bonds) constrained
666dc804011d Uploaded
fmercuri
parents:
diff changeset
48 lincs_iter = 1 ; accuracy of LINCS
666dc804011d Uploaded
fmercuri
parents:
diff changeset
49 lincs_order = 4 ; also related to accuracy
666dc804011d Uploaded
fmercuri
parents:
diff changeset
50 ; Neighborsearching
666dc804011d Uploaded
fmercuri
parents:
diff changeset
51 cutoff-scheme = $mdp.cutoffscheme
666dc804011d Uploaded
fmercuri
parents:
diff changeset
52 ns_type = grid ; search neighboring grid cells
666dc804011d Uploaded
fmercuri
parents:
diff changeset
53 nstlist = 10 ; 20 fs, largely irrelevant with Verlet
666dc804011d Uploaded
fmercuri
parents:
diff changeset
54 rcoulomb = $mdp.rcoulomb ; Short-range electrostatic cut-off
666dc804011d Uploaded
fmercuri
parents:
diff changeset
55 rlist = $mdp.rlist ; Cut-off distance for the short-range neighbor list.
666dc804011d Uploaded
fmercuri
parents:
diff changeset
56 rvdw = $mdp.rvdw ; Short-range Van der Waals cut-off
666dc804011d Uploaded
fmercuri
parents:
diff changeset
57 ; Electrostatics
666dc804011d Uploaded
fmercuri
parents:
diff changeset
58 coulombtype = $mdp.coulombtype ; method for electrostatics calculations e.g. PME
666dc804011d Uploaded
fmercuri
parents:
diff changeset
59 pme_order = 4 ; cubic interpolation
666dc804011d Uploaded
fmercuri
parents:
diff changeset
60 fourierspacing = 0.16 ; grid spacing for FFT
666dc804011d Uploaded
fmercuri
parents:
diff changeset
61 ; Temperature coupling is on
666dc804011d Uploaded
fmercuri
parents:
diff changeset
62 tcoupl = V-rescale ; modified Berendsen thermostat
666dc804011d Uploaded
fmercuri
parents:
diff changeset
63 tc-grps = $mdp.sysname ; two coupling groups - more accurate
666dc804011d Uploaded
fmercuri
parents:
diff changeset
64 tau_t = 0.1 ; time constant, in ps
666dc804011d Uploaded
fmercuri
parents:
diff changeset
65 ref_t = $mdp.temperature; reference temperature, one for each group, in K
666dc804011d Uploaded
fmercuri
parents:
diff changeset
66 ; Pressure coupling is off
666dc804011d Uploaded
fmercuri
parents:
diff changeset
67 pcoupl = no ; no pressure coupling in NVT
666dc804011d Uploaded
fmercuri
parents:
diff changeset
68 ; Periodic boundary conditions
666dc804011d Uploaded
fmercuri
parents:
diff changeset
69 pbc = xyz ; 3-D PBC
666dc804011d Uploaded
fmercuri
parents:
diff changeset
70 ; Dispersion correction
666dc804011d Uploaded
fmercuri
parents:
diff changeset
71 DispCorr = EnerPres ; account for cut-off vdW scheme
666dc804011d Uploaded
fmercuri
parents:
diff changeset
72 ; Velocity generation
666dc804011d Uploaded
fmercuri
parents:
diff changeset
73 gen_vel = yes ; assign velocities from Maxwell distribution
666dc804011d Uploaded
fmercuri
parents:
diff changeset
74 gen_temp = 300 ; temperature for Maxwell distribution
666dc804011d Uploaded
fmercuri
parents:
diff changeset
75 gen_seed = -1 ; generate a random seed
666dc804011d Uploaded
fmercuri
parents:
diff changeset
76 #end if
666dc804011d Uploaded
fmercuri
parents:
diff changeset
77 </configfile>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
78
666dc804011d Uploaded
fmercuri
parents:
diff changeset
79 </configfiles>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
80 <inputs>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
81 <param argument="gro_input" type="data" format='gro' label="GRO structure file."/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
82 <param argument="top_input" type="data" format='top' label="Topology (TOP) file."/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
83 <param argument="itp_input" type="data" format='itp' label="Position restraint file."/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
84
666dc804011d Uploaded
fmercuri
parents:
diff changeset
85 <expand macro="md_inputs"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
86 </inputs>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
87 <outputs>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
88 <data name="output1" format="gro" from_work_dir="nvt.gro">
666dc804011d Uploaded
fmercuri
parents:
diff changeset
89 <filter>str == 'gro' or str == 'both'</filter>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
90 </data>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
91 <data name="output5" format="pdb" from_work_dir="nvt.pdb">
666dc804011d Uploaded
fmercuri
parents:
diff changeset
92 <filter>str == 'pdb' or str == 'both'</filter>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
93 </data>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
94 <data name="output2" format="cpt" from_work_dir="nvt.cpt"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
95 <data name="output3" format="trr" from_work_dir="nvt.trr">
666dc804011d Uploaded
fmercuri
parents:
diff changeset
96 <filter>traj == 'trr' or traj == 'both'</filter>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
97 </data>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
98 <data name="output4" format="xtc" from_work_dir="nvt.xtc">
666dc804011d Uploaded
fmercuri
parents:
diff changeset
99 <filter>traj == 'xtc' or traj == 'both'</filter>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
100 </data>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
101 <data name="report" format="txt" from_work_dir="verbose.txt">
666dc804011d Uploaded
fmercuri
parents:
diff changeset
102 <filter>capture_log</filter>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
103 </data>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
104 </outputs>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
105 <tests>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
106 <test>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
107 <param name="gro_input" value="em.gro" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
108 <param name="top_input" value="topol_solv.top" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
109 <param name="itp_input" value="posres.itp" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
110 <param name="traj" value="xtc"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
111 <param name="str" value="both"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
112
666dc804011d Uploaded
fmercuri
parents:
diff changeset
113 <expand macro="test_params"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
114
666dc804011d Uploaded
fmercuri
parents:
diff changeset
115 <output name="output1" file="nvt.gro" ftype="gro" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
116 <output name="output2" file="nvt.cpt" ftype="cpt" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
117 <output name="output4" file="nvt.xtc" ftype="xtc" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
118 <output name="output5" file="nvt.pdb" ftype="pdb" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
119
666dc804011d Uploaded
fmercuri
parents:
diff changeset
120 </test>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
121 <test>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
122 <param name="gro_input" value="em.gro" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
123 <param name="top_input" value="topol_solv.top" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
124 <param name="itp_input" value="posres.itp" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
125 <param name="traj" value="xtc"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
126 <param name="str" value="pdb"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
127
666dc804011d Uploaded
fmercuri
parents:
diff changeset
128 <expand macro="test_params"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
129
666dc804011d Uploaded
fmercuri
parents:
diff changeset
130 <output name="output2" file="nvt.cpt" ftype="cpt" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
131 <output name="output4" file="nvt.xtc" ftype="xtc" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
132 <output name="output5" file="nvt.pdb" ftype="pdb" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
133
666dc804011d Uploaded
fmercuri
parents:
diff changeset
134 </test>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
135 <test>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
136 <param name="gro_input" value="em.gro" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
137 <param name="top_input" value="topol_solv.top" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
138 <param name="itp_input" value="posres.itp" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
139 <param name="str" value="gro"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
140 <param name="traj" value="trr"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
141
666dc804011d Uploaded
fmercuri
parents:
diff changeset
142 <expand macro="test_params"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
143
666dc804011d Uploaded
fmercuri
parents:
diff changeset
144 <output name="output1" file="nvt.gro" ftype="gro" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
145 <output name="output2" file="nvt.cpt" ftype="cpt" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
146 <output name="output3" file="nvt.trr" ftype="trr" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
147 </test>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
148 <test>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
149 <param name="gro_input" value="em.gro" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
150 <param name="top_input" value="topol_solv.top" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
151 <param name="itp_input" value="posres.itp" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
152 <param name="mdpfile" value="custom" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
153 <param name="mdp_input" value="nvt.mdp" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
154 <param name="traj" value="none"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
155 <param name="str" value="gro"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
156 <output name="output1" file="nvt.gro" ftype="gro" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
157 <output name="output2" file="nvt.cpt" ftype="cpt" compare="sim_size"/>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
158 </test>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
159 </tests>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
160 <help><![CDATA[
666dc804011d Uploaded
fmercuri
parents:
diff changeset
161
666dc804011d Uploaded
fmercuri
parents:
diff changeset
162 Upload GRO, TOP and ITP (position restraint) files for equilibration under an NVT ensemble. To take advantage of all GROMACS features, upload an MDP file with simulation parameters. Otherwise, choose parameters through the Galaxy interface. See http://manual.gromacs.org/documentation/2018/user-guide/mdp-options.html for more information on the options.
666dc804011d Uploaded
fmercuri
parents:
diff changeset
163 ]]></help>
666dc804011d Uploaded
fmercuri
parents:
diff changeset
164
666dc804011d Uploaded
fmercuri
parents:
diff changeset
165 <expand macro="citations" />
666dc804011d Uploaded
fmercuri
parents:
diff changeset
166 </tool>