From: Karson Bader (Qphysicist_at_hotmail.com)
Date: Fri May 13 2016 - 12:07:52 CDT
I have a 6 layer graphite system. The bottom three are completely fixed, the 4th layer Carbon atoms are attached to a heat bath with the outer most atoms of the 4th later fixed. The 5th and 6th layers have their outer most atoms fixed to prevent them from drifting away and the inner atoms are free to move. A basic diagram is below, where Xs are fixed, Os are free but allowed to move and 1s are the atoms coupled to the heat bath.
When I look at the temperature data from the output file, the system temperature agrees with what I have in the configuration file, but when I look at the temperature of each individual layer (by using a custom script) the coupled layer is much higher (140-150K) with layers 5 and 6 being between 50-70K. The overal average of of the layers agrees with the configuration file temp.
I'm guessing that when the program checks the temperature it looks at ALL moving atoms, not just the layer that is coupled to the heat bath.
Is there a way assign a temperature to only a subset of atoms and have ONLY those atoms looked at when the temperature is being checked?
I am using two different files to assign fixed atoms and tcoupled atoms.
My configuration file is below:
#namd configuration file
set temperature 80
#Periodic Boundary Conditions
cellBasisVector1 42.6 0.00 0.00
cellBasisVector2 0.00 36.9 0.00
#Fixed Atoms Constraint
#Berendsen Temperature Control
This archive was generated by hypermail 2.1.6 : Tue Dec 27 2016 - 23:22:07 CST