I want to heat a region of the system with a pulse duration energy of 1ns.
I tried to use “fix heat”,but it increased energy at once when N steps is run. Then I tried to set “N” to 10ps(timestep=1fs,N=10000) and set “eflux” 100 times smaller, Then ran 1000000(heat the group 100 times in 1ns with a 100 times smaller energy).
However,the target is group but not region. Because some atoms ran out of the region,when I heated the group first time. It’s not my purpose.