An enhanced version of the heat exchange algorithm with excellent energy conservation properties
arXiv:1507.07081 · doi:10.1063/1.4931597
Abstract
We propose a new algorithm for non-equilibrium molecular dynamics simulations of thermal gradients. The algorithm is an extension of the heat exchange algorithm developed by Hafskjold and co-workers [Mol. Phys. 80, 1389 (1993); Mol. Phys. 81, 251 (1994)], in which a certain amount of heat is added to one region and removed from another by rescaling velocities appropriately. Since the amount of added and removed heat is the same and the dynamics between velocity rescaling steps is Hamiltonian, the heat exchange algorithm is expected to conserve the energy. However, it has been reported previously that the original version of the heat exchange algorithm exhibits a pronounced drift in the total energy, the exact cause of which remained hitherto unclear. Here, we show that the energy drift is due to the truncation error arising from the operator splitting and suggest an additional coordinate integration step as a remedy. The new algorithm retains all the advantages of the original one whilst exhibiting excellent energy conservation as illustrated for a Lennard-Jones liquid and SPC/E water.
8 pages, 6 figures
Cited by in corpus (9)
- Equivalence of the equilibrium and the nonequilibrium molecular dynamics methods for thermal conductivity calculations: From bulk to nanowire silicon
- Exploring high thermal conductivity polymers via interpretable machine learning with physical descriptors
- Thermo-osmosis in charged nanochannels: effects of surface charge and ionic strength
- Non-equilibrium simulations of thermally induced electric fields in water
- Thermal conductivity calculation in anisotropic crystals by molecular dynamics : application to alpha-Fe2O3
- Numerical Evidence for Thermally Induced Monopoles
- Theoretical prediction of thermal polarisation
- Deep Learning Inter-atomic Potential for Thermal and Phonon Behaviour of Silicon Carbide with Quantum Accuracy
- Quantum multifractality in thermal conduction across random interfaces