Hi all users,
I am trying to apply Muller-Plathe algorithm for only one particle in a
simulation box, specifying the group of that particle atoms for fix
fix heat_swap hbn thermal/conductivity 100 y 6
'hbn' is the name of this group (I created it with type or region
options of 'group' but it didn't differ), the particle stands in the
middle of the simulation box and I am kind of fixing its position by 'mom yes
rot yes' options of 'velocity' command and also with 'fix angle hbn
momentum 10 angular'. Please also note that the other parts of the
simulation box, different than particle, is empty.
The problem is, I got 'f_heat_swap' values of 0. I think this command
divides whole simulation box to 6 slabs instead of dividing just group
of 'hbn' so it does not see any atoms in empty slabs and does not change
any velocity. But if it is the case, then why fix thermal/conductivity
command enables us to select any group that we want? I tried same
script, shrinking the box dimensions to particle dimensions and then it
started to show some energy change values, but this is not the case
that I want to see.
I am very curios that if there is someone who agrees with me about the reason of the
problem, is there a way that I can solve this or anyone who experienced it before? I will appreciate your comments.