Project ID: plumID:22.028
Source: plumed_FM5.dat
Originally used with PLUMED version: 2.7
Stable: zipped raw stdout - zipped raw stderr - stderr
Master: zipped raw stdout - zipped raw stderr - stderr
######### PLUMED file for performing RECT simulations on N-glycans ######### #MOLINFOThis command is used to provide information on the molecules that are present in your system. More detailsSTRUCTURE=reference_FM5.pdb # # Define torsion angles for FM5 # phi1_2:a file in pdb format containing a reference structureTORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O5-2,the O5 atom in residue 2. Click here for more information.@C1-2,the C1 atom in residue 2. Click here for more information.@O4-1,the O4 atom in residue 1. Click here for more information.@C4-1psi1_2:the C4 atom in residue 1. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@C1-2,the C1 atom in residue 2. Click here for more information.@O4-1,the O4 atom in residue 1. Click here for more information.@C4-1,the C4 atom in residue 1. Click here for more information.@C3-1# phi2_3:the C3 atom in residue 1. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O5-3,the O5 atom in residue 3. Click here for more information.@C1-3,the C1 atom in residue 3. Click here for more information.@O4-2,the O4 atom in residue 2. Click here for more information.@C4-2psi2_3:the C4 atom in residue 2. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@C1-3,the C1 atom in residue 3. Click here for more information.@O4-2,the O4 atom in residue 2. Click here for more information.@C4-2,the C4 atom in residue 2. Click here for more information.@C3-2# phi3_4:the C3 atom in residue 2. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O5-4,the O5 atom in residue 4. Click here for more information.@C1-4,the C1 atom in residue 4. Click here for more information.@O3-3,the O3 atom in residue 3. Click here for more information.@C3-3psi3_4:the C3 atom in residue 3. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@C1-4,the C1 atom in residue 4. Click here for more information.@O3-3,the O3 atom in residue 3. Click here for more information.@C3-3,the C3 atom in residue 3. Click here for more information.@C2-3# phi3_5:the C2 atom in residue 3. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O5-5,the O5 atom in residue 5. Click here for more information.@C1-5,the C1 atom in residue 5. Click here for more information.@O6-3,the O6 atom in residue 3. Click here for more information.@C6-3psi3_5:the C6 atom in residue 3. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@C1-5,the C1 atom in residue 5. Click here for more information.@O6-3,the O6 atom in residue 3. Click here for more information.@C6-3,the C6 atom in residue 3. Click here for more information.@C5-3omega3_5:the C5 atom in residue 3. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O6-3,the O6 atom in residue 3. Click here for more information.@C6-3,the C6 atom in residue 3. Click here for more information.@C5-3,the C5 atom in residue 3. Click here for more information.@O5-3# phi5_7:the O5 atom in residue 3. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O5-7,the O5 atom in residue 7. Click here for more information.@C1-7,the C1 atom in residue 7. Click here for more information.@O6-5,the O6 atom in residue 5. Click here for more information.@C6-5psi5_7:the C6 atom in residue 5. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@C1-7,the C1 atom in residue 7. Click here for more information.@O6-5,the O6 atom in residue 5. Click here for more information.@C6-5,the C6 atom in residue 5. Click here for more information.@C5-5omega5_7:the C5 atom in residue 5. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O6-5,the O6 atom in residue 5. Click here for more information.@C6-5,the C6 atom in residue 5. Click here for more information.@C5-5,the C5 atom in residue 5. Click here for more information.@O5-5# phi5_6:the O5 atom in residue 5. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O5-6,the O5 atom in residue 6. Click here for more information.@C1-6,the C1 atom in residue 6. Click here for more information.@O3-5,the O3 atom in residue 5. Click here for more information.@C3-5psi5_6:the C3 atom in residue 5. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@C1-6,the C1 atom in residue 6. Click here for more information.@O3-5,the O3 atom in residue 5. Click here for more information.@C3-5,the C3 atom in residue 5. Click here for more information.@C2-5# phi1_8:the C2 atom in residue 5. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O5-8,the O5 atom in residue 8. Click here for more information.@C1-8,the C1 atom in residue 8. Click here for more information.@O6-1,the O6 atom in residue 1. Click here for more information.@C6-1psi1_8:the C6 atom in residue 1. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@C1-8,the C1 atom in residue 8. Click here for more information.@O6-1,the O6 atom in residue 1. Click here for more information.@C6-1,the C6 atom in residue 1. Click here for more information.@C5-1omega1_8:the C5 atom in residue 1. Click here for more information.TORSIONCalculate a torsional angle. More detailsATOMS=the four atoms involved in the torsional angle@O6-1,the O6 atom in residue 1. Click here for more information.@C6-1,the C6 atom in residue 1. Click here for more information.@C5-1,the C5 atom in residue 1. Click here for more information.@O5-1# # Define puckering for each monosaccharide only to print it to file # puck1:the O5 atom in residue 1. Click here for more information.PUCKERINGCalculate sugar pseudorotation coordinates. More detailsATOMS=the five or six atoms of the sugar ring in the proper order@O5-1,the O5 atom in residue 1. Click here for more information.@C1-1,the C1 atom in residue 1. Click here for more information.@C2-1,the C2 atom in residue 1. Click here for more information.@C3-1,the C3 atom in residue 1. Click here for more information.@C4-1,the C4 atom in residue 1. Click here for more information.@C5-1puck2:the C5 atom in residue 1. Click here for more information.PUCKERINGCalculate sugar pseudorotation coordinates. More detailsATOMS=the five or six atoms of the sugar ring in the proper order@O5-2,the O5 atom in residue 2. Click here for more information.@C1-2,the C1 atom in residue 2. Click here for more information.@C2-2,the C2 atom in residue 2. Click here for more information.@C3-2,the C3 atom in residue 2. Click here for more information.@C4-2,the C4 atom in residue 2. Click here for more information.@C5-2puck3:the C5 atom in residue 2. Click here for more information.PUCKERINGCalculate sugar pseudorotation coordinates. More detailsATOMS=the five or six atoms of the sugar ring in the proper order@O5-3,the O5 atom in residue 3. Click here for more information.@C1-3,the C1 atom in residue 3. Click here for more information.@C2-3,the C2 atom in residue 3. Click here for more information.@C3-3,the C3 atom in residue 3. Click here for more information.@C4-3,the C4 atom in residue 3. Click here for more information.@C5-3puck4:the C5 atom in residue 3. Click here for more information.PUCKERINGCalculate sugar pseudorotation coordinates. More detailsATOMS=the five or six atoms of the sugar ring in the proper order@O5-4,the O5 atom in residue 4. Click here for more information.@C1-4,the C1 atom in residue 4. Click here for more information.@C2-4,the C2 atom in residue 4. Click here for more information.@C3-4,the C3 atom in residue 4. Click here for more information.@C4-4,the C4 atom in residue 4. Click here for more information.@C5-4puck5:the C5 atom in residue 4. Click here for more information.PUCKERINGCalculate sugar pseudorotation coordinates. More detailsATOMS=the five or six atoms of the sugar ring in the proper order@O5-5,the O5 atom in residue 5. Click here for more information.@C1-5,the C1 atom in residue 5. Click here for more information.@C2-5,the C2 atom in residue 5. Click here for more information.@C3-5,the C3 atom in residue 5. Click here for more information.@C4-5,the C4 atom in residue 5. Click here for more information.@C5-5puck6:the C5 atom in residue 5. Click here for more information.PUCKERINGCalculate sugar pseudorotation coordinates. More detailsATOMS=the five or six atoms of the sugar ring in the proper order@O5-6,the O5 atom in residue 6. Click here for more information.@C1-6,the C1 atom in residue 6. Click here for more information.@C2-6,the C2 atom in residue 6. Click here for more information.@C3-6,the C3 atom in residue 6. Click here for more information.@C4-6,the C4 atom in residue 6. Click here for more information.@C5-6puck7:the C5 atom in residue 6. Click here for more information.PUCKERINGCalculate sugar pseudorotation coordinates. More detailsATOMS=the five or six atoms of the sugar ring in the proper order@O5-7,the O5 atom in residue 7. Click here for more information.@C1-7,the C1 atom in residue 7. Click here for more information.@C2-7,the C2 atom in residue 7. Click here for more information.@C3-7,the C3 atom in residue 7. Click here for more information.@C4-7,the C4 atom in residue 7. Click here for more information.@C5-7puck8:the C5 atom in residue 7. Click here for more information.PUCKERINGCalculate sugar pseudorotation coordinates. More detailsATOMS=the five or six atoms of the sugar ring in the proper order@O5-8,the O5 atom in residue 8. Click here for more information.@C1-8,the C1 atom in residue 8. Click here for more information.@C2-8,the C2 atom in residue 8. Click here for more information.@C3-8,the C3 atom in residue 8. Click here for more information.@C4-8,the C4 atom in residue 8. Click here for more information.@C5-8# # Apply 1D metadynamics on each torsion angle. Different bias potentials are applied in different replicas, according to the RECT keyword. # metaD_phi1_2:the C5 atom in residue 8. Click here for more information.METADUsed to performed metadynamics on one or more collective variables. More detailsARG=phi1_2the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_phi1_2 metaD_psi1_2:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=psi1_2the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_psi1_2 # metaD_phi2_3:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=phi2_3the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_phi2_3 metaD_psi2_3:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=psi2_3the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_psi2_3 # metaD_phi3_4:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=phi3_4the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_phi3_4 metaD_psi3_4:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=psi3_4the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_psi3_4 # metaD_phi3_5:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=phi3_5the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_phi3_5 metaD_psi3_5:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=psi3_5the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_psi3_5 metaD_omega3_5:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=omega3_5the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_omega3_5 # metaD_phi5_6:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=phi5_6the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_phi5_6 metaD_psi5_6:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=psi5_6the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_psi5_6 # metaD_phi5_7:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=phi5_7the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_phi5_7 metaD_psi5_7:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=psi5_7the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_psi5_7 metaD_omega5_7:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=omega5_7the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_omega5_7 # metaD_phi1_8:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=phi1_8the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_phi1_8 metaD_psi1_8:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=psi1_8the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_psi1_8 metaD_omega1_8:a file in which the list of added hills is storedMETADUsed to performed metadynamics on one or more collective variables. More detailsARG=omega1_8the input for this action is the scalar output from one or more other actionsTAU=4.0in well tempered metadynamics, sets height to (k_B Delta T*pace*timestep)/tauPACE=500the frequency for hill additionSIGMA=0.35the widths of the Gaussian hillsGRID_MIN=-pithe lower bounds for the gridGRID_MAX=pithe upper bounds for the gridGRID_BIN=200the number of bins for the gridTEMP=310.15the system temperature - this is only needed if you are doing well-tempered metadynamicsRECT=1,1.2,1.46,1.82,2.3,2.94,3.78,4.89,6.34,8.23,10.7,14list of bias factors for all the replicasFILE=HILLS_omega1_8 # # PRINT torsion angle values and puckering coordinates to separate files #a file in which the list of added hills is storedPrint quantities to a file. More detailsARG=phi1_2,psi1_2,phi2_3,psi2_3,phi3_4,psi3_4,phi3_5,psi3_5,omega3_5,phi5_6,psi5_6,phi5_7,psi5_7,omega5_7,phi1_8,psi1_8,omega1_8the input for this action is the scalar output from one or more other actionsSTRIDE=200the frequency with which the quantities of interest should be outputFILE=COLVAR #the name of the file on which to output these quantitiesPrint quantities to a file. More detailsARG=puck1.theta,puck2.theta,puck3.theta,puck4.theta,puck5.theta,puck6.theta,puck7.theta,puck8.theta,puck1.phi,puck2.phi,puck3.phi,puck4.phi,puck5.phi,puck6.phi,puck7.phi,puck8.phithe input for this action is the scalar output from one or more other actionsSTRIDE=200the frequency with which the quantities of interest should be outputFILE=COLVAR_thetathe name of the file on which to output these quantities