Thursday, June 26, 2014

Phi-4 model for Nose-Hoover thermostatting scheme

For the last couple of weeks I have been working on the $\phi^4$ conduction model. The model is very unique. Many works have followed Aoki & Kusnezov's work on finding thermal conductivity of this system. Prof. William Hoover has used this model to show that the Hamiltonian thermostats fail to promote heat flow.

I have written a code for those in need for $\phi^4$ model in thermal equilibrium. Thermal equilibrium is implemented using Nose-Hoover thermostat. Only a part of the system is thermostatted. The c-code is attached for those in need.

Feel free to email/comment if you find any error or have any doubts.

Link is: Nose_Hoover_Phi4_Model_C_File