Domain decomposition of the two-temperature model in dl_poly_4

ABSTRACTIncluding the effects of excited electrons in classical simulations, at the level of the two-temperature model, involves the coupling of a grid-based finite-difference solver for a heat diffusion equation and classical molecular dynamics simulations with an inhomogeneous thermostat. Simulation of large systems requires domain decomposition of both particle-based and grid-based techniques. Starting with the CCP5 flagship code dl_poly_4 as the domain-decomposed molecular dynamics code, we devised a method to divide up temperature grids among processor cores in a similar fashion, including the appropriate communications between cores to deal with boundaries between grid divisions. This article gives the outline of how the domain decomposition of the temperature grids was achieved, as well as some example applications of the two-temperature model implementation in dl_poly_4.