Hi there,
I am using command fix thermal/conductivity to compute the thermal conducitivity of carbon nanotube. However, I was unable to get the correct temperature profile after arriving at the steady state. Can anybody tell me what part I used was wrong? My guesssing is the problem occurs at fix ave/spatial part to output the temp profile. or can anybody give me your already checked code (as simple as it can) to let me learn? Thanks for your help.
units metal
dimension 3
boundary p p p
atom_style atomic
pair_style tersoff
read_data data.dat
pair_coeff * * SiC.tersoff C
mass 1 12
velocity all create 300 873586443 dist gaussian
fix 1 all nvt 300.0 300.0 0.5 drag 0.4
#output of varaibles
thermo_style custom step temp pe ke etotal pxx pyy pzz
thermo 1000
timestep 0.001
run 10000
unfix 1
fix 1 all nve
fix 2 all thermal/conductivity 100 z 40
compute ke all ke/atom
variable temp atom c_ke[]/1.5
fix 3 all ave/spatial 10000 1 100000 z lower 0.25 v_temp &
file tmp.profile units reduced
dump 1 all custom 100000 dump.* tag x y z
run 1000000
Best,
AC