Dear Steve
Thank you!
I have changed the code, and get the temp profile
So “is it ok” means “is it right in RNEMD principle”
CL Wang
Dear Steve
Thank you!
I have changed the code, and get the temp profile
So “is it ok” means “is it right in RNEMD principle”
CL Wang
MP basically is a RNEMD method with controlled heat flux, so you need the periodic boundary condition, namely, doubling the domain size you are interested in. (for instance, http://iopscience.iop.org/0957-4484/12/1/305)
There are other methods utilizing adiabatic boundary condition, for example, http://prb.aps.org/abstract/PRB/v76/i8/e085424, but you need to modify the code more to meet the method’s requirement.
LC Liu