There is a newer version of the record available.

Published January 26, 2022 | Version v1
Dataset Open

On the effects of the degrees of freedom on calculating diffusion properties in nanoporous materials

  • 1. Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1951 Sion, Valais, Switzerland
  • 2. Department of Physics, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, Realfagbygget D5-149, 7491 Trondheim, Norway

* Contact person

Description

If one carries out a molecular simulation of N particles using periodic boundary conditions, linear momentum is conserved and hence the number of degrees of freedom is set to 3N-3. In most programs, this number of degrees of freedom is the default setting. However, if one carries out a molecular simulation in an external field, one needs to ensure that degrees of freedom are changed from this default setting to 3N, as in an external field the velocity of the center of mass can change. Using the correct degrees of freedom is important in calculating the temperature and in some algorithms to simulate at constant temperature. For sufficiently large systems, the difference between 3N and 3N-3 is negligible in the way. However, there are systems in which the comparison with experimental data requires molecular dynamics simulations of a small number of particles. In this work, we illustrate the effect of an incorrect setting of degrees of freedom in molecular dynamic simulations studying the diffusion properties of guest molecules in nanoporous materials. We show that previously published results have reported a surprising diffusion dependence on the loading, which could be traced back to an incorrect setting of the degrees of freedom. As the correct settings are convoluted and counter-intuitive in some of the most commonly used molecular dynamics programs, we carried out a systematic study on the consequences of the various commonly used (incorrect) settings. Our conclusion, is that for systems smaller than 50 particles the results are most likely unreliable as these are either performed at an incorrect temperature or the temperature is incorrectly used in some of the results. Furthermore, a novel and efficient method to calculate diffusion coefficients of guest molecules into nanoporous materials at zero loading conditions is introduced.

Files

File preview

files_description.md

All files

Files (31.0 MiB)

Name Size
md5:5cd1e9c15ce028c7c6e7f49729b8006e
894 Bytes Preview Download
md5:9fa11d0786662327f084c55b9df23d9d
5.4 MiB Download
md5:5b5a35b5b6fb9d336b387e34ffedf30c
2.0 MiB Download
md5:56b8c6a36bfb3dfe6594671d00c676d5
7.5 MiB Download
md5:0bee05ef4018dea8772af286580fd904
15.5 MiB Download
md5:fed082a2dba46d4210ed8390f9862598
100.0 KiB Download
md5:bd8edbb288b641fabf978b0204b90e6d
50.0 KiB Download
md5:067100af1c43313da00d61cd1ae91168
50.0 KiB Download
md5:674ec7c23cdad7f3e4a44846a5d49798
360.0 KiB Download
md5:50129b1059b5fc938682217eaadfb71d
1.8 KiB Preview Download

References

Preprint
H. Xu, R. Cabriolu, B. Smit, Journal of Chemical Theory and Computation, to be submitted (2022)