Dear Lammps Users,
I want to compute the temperature distribution in a nanochannel filled with diatomic gas molecules. My system is equilibrated at 300 K. When I use “compute temp” or “compute temp/com”, it indeed shows that the overall temperature in the gas domain is 300 K.
As you can see in the attached lammps input, I tried to derive the temperature distribution across the nanochannel using “compute chunk/atom”+ “fix ave/chunk” and also using “compute temp/chunk”, but using both approaches show that the temperature across the channel varies around 250 K.
The thing is previously I used the same approach for monoatomic gas and I haven’t faced such a problem in that case.
I highly appreciate it if someone could tell me how I can get the correct temperature distribution in a system with diatomic gas molecules.
#------------------------------------------input file--------------------------------------------
read_restart restart.Eq_after_thermos
group NiB region Ni_B
group NiT region Ni_T
#----------------------------------------------FORCE FIELD--------------------------------------------------------------