| This is part of the opes module | |
| It is only available if you configure PLUMED with ./configure –enable-modules=opes . Furthermore, this feature is still being developed so take care when using it and report any problems on the mailing list. |
Linear expansion, according to a parameter lambda.
This can be used e.g. for thermodynamic integration, or for multibaric simulations, in which case lambda=pressure. It can also be used for multithermal simulations, but for simplicity it is more convenient to use ECV_MULTITHERMAL.
The difference in Hamiltonian \(\Delta U\) is expected as ARG.
\[ \Delta u_\lambda=\beta \lambda \Delta U\, . \]
Use the DIMENSIONLESS flag to avoid multiplying for the inverse temperature \(\beta\).
Typical multibaric simulation:
vol: VOLUME ecv: ECV_LINEAR ...ARG=volcompulsory keyword the label of the Hamiltonian difference \f$\Delta U\f$TEMP=300compulsory keyword ( default=-1 ) temperature.LAMBDA=0.06022140857*2000 #2 kbarcompulsory keyword ( default=0 ) the lambda at which the underlying simulation runsMIN_LAMBDA=0.06022140857 #1 bar( default=0 ) the minimum of the lambda rangeMAX_LAMBDA=0.06022140857*4000 #4 kbar ... opes: OPES_EXPANDED( default=1 ) the maximum of the lambda rangeARG=ecv.volcompulsory keyword the label of the ECVs that define the expansion.PACE=500compulsory keyword how often the bias is updated
Typical thermodynamic integration:
DeltaU: EXTRACVNAME=energy_difference ecv: ECV_LINEARcompulsory keyword name of the CV as computed by the MD engineARG=DeltaUcompulsory keyword the label of the Hamiltonian difference \f$\Delta U\f$TEMP=300 opes: OPES_EXPANDEDcompulsory keyword ( default=-1 ) temperature.ARG=ecv.*compulsory keyword the label of the ECVs that define the expansion.PACE=100compulsory keyword how often the bias is updated
Notice that by defauly LAMBDA=0, MIN_LAMBDA=0 and MAX_LAMBDA=1, which is the typical case for thermodynamic integration.
| TEMP | ( default=-1 ) temperature. If not specified tries to get it from MD engine |
| ARG | the label of the Hamiltonian difference \(\Delta U\) |
| LAMBDA | ( default=0 ) the lambda at which the underlying simulation runs |
| NUMERICAL_DERIVATIVES | ( default=off ) calculate the derivatives for these quantities numerically |
| DIMENSIONLESS | ( default=off ) ARG is considered dimensionless rather than an energy, thus is not multiplied by \(\beta\) |
| MIN_LAMBDA | ( default=0 ) the minimum of the lambda range |
| MAX_LAMBDA | ( default=1 ) the maximum of the lambda range |
| STEPS_LAMBDA | uniformly place the lambda values, for a total of STEPS_LAMBDA |
| SET_ALL_LAMBDAS | manually set all the lamdbas |