Cell vectors are passed correctly
PLUMED must receive the cell vectors from the MD code in order to calculate CVs correctly.
To test that cell vectors are passed correctly to PLUMED we run a short trajectory and output the cell vectors
that are passed to PLUMED using the following command:
c:CELLCalculate the components of the simulation cell More details
Trajectory
Input and output files for the test calculation are available in this zip archive
Results
The table below includes some of the results from the calculation. The columns contain:
- The values for the cell vectors that were obtained from the MD code, $x_{md}$.
- The values for the cell vectors that were obtained from PLUMED, $x_{pl}$.
- The tolerances that were used when comparing these quantities, $\delta$.
- The values of $100\frac{\vert x_{md} - x_{pl}\vert }{ \delta }$.
If the PLUMED interface is working correctly the first two sets of numbers should be identical and the final column should be filled with zeros.
MD code output | PLUMED output | Tolerance | % Difference |
---|---|---|---|
-0.2694 0.0000 0.2694 0.0000 0.2694 0.2694 -0.2694 0.2694 0.0000 | -0.2694 0.0000 0.2694 0.0000 0.2694 0.2694 -0.2694 0.2694 0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0220 0.0000 0.0220 0.0000 0.0220 0.0220 0.0220 0.0220 0.0000 |
-0.2694 -0.0000 0.2694 0.0000 0.2694 0.2694 -0.2694 0.2694 -0.0000 | -0.2694 -0.0000 0.2694 0.0000 0.2694 0.2694 -0.2694 0.2694 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0323 0.0000 0.0323 0.0426 0.0103 0.0103 0.0323 0.0323 0.0000 |
-0.2694 0.0000 0.2694 0.0000 0.2694 0.2694 -0.2694 0.2694 -0.0000 | -0.2694 0.0000 0.2694 0.0000 0.2694 0.2694 -0.2694 0.2694 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0185 0.0000 0.0185 0.0480 0.0293 0.0295 0.0180 0.0181 0.0005 |
-0.2694 0.0000 0.2694 0.0000 0.2694 0.2694 -0.2694 0.2694 -0.0000 | -0.2694 0.0000 0.2694 0.0000 0.2694 0.2694 -0.2694 0.2694 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0378 0.0001 0.0377 0.0140 0.0242 0.0237 0.0357 0.0361 0.0020 |
-0.2695 0.0000 0.2695 0.0000 0.2695 0.2695 -0.2695 0.2695 -0.0000 | -0.2695 0.0000 0.2695 0.0000 0.2695 0.2695 -0.2695 0.2695 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0082 0.0002 0.0083 0.0374 0.0447 0.0458 0.0133 0.0122 0.0049 |
-0.2695 0.0000 0.2695 0.0000 0.2695 0.2695 -0.2695 0.2695 -0.0000 | -0.2695 0.0000 0.2695 0.0000 0.2695 0.2695 -0.2695 0.2695 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0170 0.0003 0.0171 0.0139 0.0292 0.0313 0.0270 0.0249 0.0097 |
-0.2696 0.0000 0.2696 0.0000 0.2696 0.2696 -0.2696 0.2696 -0.0000 | -0.2696 0.0000 0.2696 0.0000 0.2696 0.2696 -0.2696 0.2696 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0147 0.0006 0.0145 0.0385 0.0210 0.0246 0.0025 0.0009 0.0166 |
-0.2697 0.0000 0.2697 0.0000 0.2697 0.2697 -0.2697 0.2697 -0.0000 | -0.2697 0.0000 0.2697 0.0000 0.2697 0.2697 -0.2697 0.2697 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0095 0.0009 0.0096 0.0052 0.0103 0.0157 0.0365 0.0313 0.0261 |
-0.2698 0.0000 0.2698 0.0000 0.2698 0.2698 -0.2698 0.2698 -0.0000 | -0.2698 0.0000 0.2698 0.0000 0.2698 0.2698 -0.2698 0.2698 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0153 0.0014 0.0150 0.0072 0.0140 0.0065 0.0245 0.0172 0.0385 |
-0.2699 0.0000 0.2699 0.0000 0.2699 0.2699 -0.2699 0.2699 -0.0000 | -0.2699 0.0000 0.2699 0.0000 0.2699 0.2699 -0.2699 0.2699 -0.0000 | 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 0.0010 | 0.0056 0.0019 0.0060 0.0375 0.0354 0.0454 0.0386 0.0483 0.0461 |