Project ID: plumID:19.018
Source: ES/plumed_excitedstates.dat
Originally used with PLUMED version: 2.3
Stable: zipped raw stdout - zipped raw stderr - stderr
Master: zipped raw stdout - zipped raw stderr - stderr
WHOLEMOLECULESThis action is used to rebuild molecules that can become split by the periodic boundary conditions. More details STRIDE the frequency with which molecules are reassembled=1 ENTITY0the atoms that make up a molecule that you wish to align=1-107 T1 : TORSIONCalculate one or multiple torsional angles. More details ATOMSthe four atoms involved in the torsional angle=4,3,36,37 T2 : TORSIONCalculate one or multiple torsional angles. More details ATOMSthe four atoms involved in the torsional angle=58,57,25,26 T3 : TORSIONCalculate one or multiple torsional angles. More details ATOMSthe four atoms involved in the torsional angle=79,78,46,47 # this the lambda value for current window lambdaCONSTANTCreate a constant value that can be passed to actions This action has hidden defaults. More details VALUEthe single number that you would like to store=0.00 :
# first we define V0 MATHEVALAn alias to the CUSTOM function that can also be used to calaculate combinations of variables using a custom expression. More details ... LABELa label for the action so that its output can be referenced in the input to other actions=V0 ARGthe values input to this function=T3 VARthe names to give each of the arguments in the function=x FUNCthe function you wish to evaluate=13.68-0.057*cos(x-pi)-48.06*(cos(x-pi))^2+0.17*(cos(x-pi))^3+42.2*(cos(x-pi))^4-0.0809*(cos(x-pi))^5 PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO ... MATHEVAL
# then V1 MATHEVALAn alias to the CUSTOM function that can also be used to calaculate combinations of variables using a custom expression. More details ... LABELa label for the action so that its output can be referenced in the input to other actions=V1 ARGthe values input to this function=T3 VARthe names to give each of the arguments in the function=x FUNCthe function you wish to evaluate=341.-0.463*cos(x-pi)-92.73*(cos(x-pi))^2+1.23*(cos(x-pi))^3+55.6*(cos(x-pi))^4-0.803*(cos(x-pi))^5 PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO ... MATHEVAL
# This is Vlambda for current window Vlambda : MATHEVALAn alias to the CUSTOM function that can also be used to calaculate combinations of variables using a custom expression. More details ARGthe values input to this function=V0,V1,lambda VARthe names to give each of the arguments in the function=x,y,z FUNCthe function you wish to evaluate=(1.0-z)*x+z*y PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO # These are Vlambda for nearest neighbors Vlambdap : MATHEVALAn alias to the CUSTOM function that can also be used to calaculate combinations of variables using a custom expression. More details ARGthe values input to this function=V0,V1,lambda VARthe names to give each of the arguments in the function=x,y,z FUNCthe function you wish to evaluate=(1.0-z-0.02)*x+(z+0.02)*y PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO Vlambdam : MATHEVALAn alias to the CUSTOM function that can also be used to calaculate combinations of variables using a custom expression. More details ARGthe values input to this function=V0,V1,lambda VARthe names to give each of the arguments in the function=x,y,z FUNCthe function you wish to evaluate=(1.0-z+0.02)*x+(z-0.02)*y PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO # activate Vlambda bias for current window Vlbias : BIASVALUETakes the value of one variable and use it as a bias More details ARGthe labels of the scalar/vector arguments whose values will be used as a bias on the system=Vlambda METADUsed to performed metadynamics on one or more collective variables. More details ARGthe labels of the scalars on which the bias will act=T3 SIGMAthe widths of the Gaussian hills=0.2 HEIGHTthe heights of the Gaussian hills=1.2 PACEthe frequency for hill addition=500 FILE a file in which the list of added hills is stored=HILLS_0.00 BIASFACTORuse well tempered metadynamics and use this bias factor=25.0 TEMPthe system temperature - this is only needed if you are doing well-tempered metadynamics=300.0 GRID_MINthe lower bounds for the grid=-pi GRID_MAXthe upper bounds for the grid=pi # print out useful stuff PRINTPrint quantities to a file. More details ARGthe labels of the values that you would like to print to the file=T3,Vlambdam,Vlambda,Vlambdap,T1,T2 STRIDE the frequency with which the quantities of interest should be output=500 FILEthe name of the file on which to output these quantities=COLVAR_0.00