Hello developer, my model setup is shown below. There are vacuum layer, Newton layer, constant temperature layer and fixed layer. I use ttm/mod command to simulate femtosecond laser with two-temperature model。
compute ke all ke/atom
variable temp atom c_ke/1.5/8.625e-5
compute layers all chunk/atom bin/1d z lower 0.025 units reduced
fix 2 all ave/chunk 2 5 100 layers v_temp file temp.txt
This temperature calculation I selected all the areas, including the fixed layer, but the speed of the fixed layer is 0, so I think there are some problems in temperature calculation.
。
But if I choose to calculate temperature using only Newtonian layers, does his chunk/atom layer divide only Newtonian layers?。
And his output file shows that atoms and temperatures are also present in the vacuum layer region.。