A new and effective method for thermostatting confined fluids
- DNRF Center “Glass and Time,” IMFUFA, Department of Science, Systems and Models, Roskilde University, DK-4000 Roskilde (Denmark)
- School of Applied Sciences, RMIT University, Melbourne, Victoria 3001 (Australia)
We present a simple thermostatting method suitable for nanoconfined fluid systems. Two conventional strategies involve thermostatting the fluid directly or employing a thermal wall that couples only the wall atoms with the thermostat. When only a thermal wall is implemented, the temperature control of the fluid is true to the actual experiment and the heat is transferred from the fluid to the walls. However, for large or complex systems it can often be computationally prohibitive to employ thermal walls. To overcome this limitation many researchers choose to freeze wall atoms and instead apply a synthetic thermostat to the fluid directly through the equations of motion. This, however, can have serious consequences for the mechanical, thermodynamic, and dynamical properties of the fluid by introducing unphysical behaviour into the system [Bernardi et al., J. Chem. Phys. 132, 244706 (2010)]. In this paper, we propose a simple scheme which enables working with both frozen walls and naturally thermostatted liquids. This is done by superimposing the walls with oscillating particles, which vibrate on the edge of the fluid control volume. These particles exchange energy with the fluid molecules, but do not interact with wall atoms or each other, thus behaving as virtual particles. Their displacements violate the Lindemann criterion for melting, in such a way that the net effect would not amount to an additional confining surface. One advantage over standard techniques is the reduced computational cost, particularly for large walls, since they can be kept rigid. Another advantage over accepted strategies is the opportunity to freeze complex charged walls such as β-cristobalite. The method furthermore overcomes the problem with polar fluids such as water, as thermalized charged surfaces require higher spring constants to preserve structural stability, due to the effects of strong Coulomb interactions, thus inevitably degrading the thermostatting efficiency.
- OSTI ID:
- 22255174
- Journal Information:
- Journal of Chemical Physics, Vol. 140, Issue 5; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
A stochastic thermostat algorithm for coarse-grained thermomechanical modeling of large-scale soft matters: Theory and application to microfilaments
Mechanism of Thermal Reversal of the (Fulvalene)tetracarbonyldiruthenium Photoisomerization: Toward Molecular Solar-Thermal Energy Storage