The numerical investigation of the statics and dynamics of systems in non-equilibrium in general, and under shear flow in particular, has become more and more common. However, not all the numerical methods developed to simulate equilibrium systems can be successfully adapted to out-of-equilibrium cases. This is especially true for thermostats. Indeed, even though thermostats developed to work under equilibrium conditions sometimes display good agreement with rheology experiments, their performance rapidly degrades beyond weak dissipation and small shear rates. Here we focus on gauging the relative performances of three thermostats, Langevin, dissipative particle dynamics, and Bussi-Donadio-Parrinello under varying parameters and external conditions. We compare their effectiveness by looking at different observables and clearly demonstrate that choosing the right thermostat (and its parameters) requires a careful evaluation of, at least, temperature, density and velocity profiles. We also show that small modifications of the Langevin and DPD thermostats greatly enhance their performance in a wide range of parameters.
On the effect of the thermostat in non-equilibrium molecular dynamics simulations / Ruiz-Franco, José; Rovigatti, Lorenzo; Zaccarelli, Emanuela. - In: THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER. - ISSN 1292-8941. - 41:7(2018), p. 80. [10.1140/epje/i2018-11689-4]
On the effect of the thermostat in non-equilibrium molecular dynamics simulations
Ruiz-Franco, José
;Rovigatti, Lorenzo;Zaccarelli, Emanuela
2018
Abstract
The numerical investigation of the statics and dynamics of systems in non-equilibrium in general, and under shear flow in particular, has become more and more common. However, not all the numerical methods developed to simulate equilibrium systems can be successfully adapted to out-of-equilibrium cases. This is especially true for thermostats. Indeed, even though thermostats developed to work under equilibrium conditions sometimes display good agreement with rheology experiments, their performance rapidly degrades beyond weak dissipation and small shear rates. Here we focus on gauging the relative performances of three thermostats, Langevin, dissipative particle dynamics, and Bussi-Donadio-Parrinello under varying parameters and external conditions. We compare their effectiveness by looking at different observables and clearly demonstrate that choosing the right thermostat (and its parameters) requires a careful evaluation of, at least, temperature, density and velocity profiles. We also show that small modifications of the Langevin and DPD thermostats greatly enhance their performance in a wide range of parameters.File | Dimensione | Formato | |
---|---|---|---|
Ruiz-Franco_On_the_effect_2018.pdf
solo gestori archivio
Tipologia:
Versione editoriale (versione pubblicata con il layout dell'editore)
Licenza:
Tutti i diritti riservati (All rights reserved)
Dimensione
505.49 kB
Formato
Adobe PDF
|
505.49 kB | Adobe PDF | Contatta l'autore |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.