lammps中fix heat命令的疑问
大家好，小弟最近想用lammps中的fix heat命令模拟恒定热流的条件，但是看manual中的fix heat命令未能理解该命令的内涵。
Manual中有这样一句话“Note that since eflux has per-time units (i.e. it is a flux), this means that a larger value of N will add/subtract a larger amount of energy each time the fix is invoked.”其中eflux的单位是eV/ps（在metal单位下），我既然已经确定了热流的大小，那么N的取值大小有何意义呢？我实在是没有理解N在该条命令中的作用，还望大家指教。
另，如果是恒定热流，命令中的N应该如何取值呢？
欢迎lammps users前来讨论！ Manual中有这样一句话“Note that since eflux has per-time units (i.e. it is a flux), this means that a larger value of N will add/subtract a larger amount of energy each time the fix is invoked.”其中eflux的单位是eV/ps（在metal单位下），我既然已经确定了热流的大小，那么N的取值大小有何意义呢？我实在是没有理解N在该条命令中的作用，还望大家指教。
另，如果是恒定热流，命令中的N应该如何取值呢？
N的值确实不影响热流大小。它只影响加热流的频率和每次加的能量。假定热流大小是 1eV/ps，积分步长是 1 fs。
如果 N=1，则每一步加一次热流，每次向高温区输入 0.001 eV 的能量，从低温区带走同样的能量。
如果 N=10，则每 10 步加一次热流，每次向高温区输入 0.01 eV 的能量，从低温区带走同样的能量。
如果 N=100，则每 100 步加一次热流，每次向高温区输入 0.1 eV 的能量，从低温区带走同样的能量。
建议：取 N = 1，这样比较温和。 谢谢您的回复，这样一解释我就明白多了。
页:
[1]