units metal
boundary p p p
atom_style atomic
lattice fcc 3.61
region wor block 0 50 0 10 0 50
region box block 0 50 -10 50 0 50
create_box 1 box
create_atoms 1 region wor
pair_style eam
pair_coeff* * Cu_u3.eam
run 0
neighbor 0.5 bin
neigh_modify every 5 delay 0 check yes
regionhot block 24 26 9 10 24 26 units box
grouphot region hot
thermo100
timestep 0.005
compute 110 all ke/atom
dump 200 all custom 10 dump.po xs ys zs c_110
fix 40 all nve
run 1000
fix 30 hot temp/rescale 10 300 1000 0.05 0.5
run 10000
目的是想在这个Cu片中间加热，观察热的传递，是不是看原子的动能/震动就好了？
还有fix temp/rescle 后两个参数应该怎么理解呀 谢啦

最后两个值，manuel里面说的很明白了吧
window = only rescale if temperature is outside this window (temperature units)
fraction = rescale to target temperature by this fraction
Rescaling is only performed if the difference between the current and desired temperatures is greater than the window value. The amount of rescaling that is applied is a fraction (from 0.0 to 1.0) of the difference between the actual and desired temperature. E.g. if fraction = 1.0, the temperature is reset to exactly the desired value.

units real
dimension 3
boundary p p p
atom_style charge
#neighbor 2.0 bin
# read topology
read_data sioh.data
# Force-field parameters
# O=1 Si=2 H=3
pair_style reax/c NULL checkqeq yes
pair_coeff * * ffield.reax.SIOH O Si H
fix 5 all nve
fix 10 all qeq/reax 1 0.0 10.0 1.0e-6 reax/c
run 1000