Dear all,

When i was comparing the thermal conductivity of an anisotropic system calculated by emd(gree-kubo) method and nemd method, a question occured to me that does the nemd method calculated thermal conductivity has any meaning for the anisotropic system?

As we all know that the fourier heat transform formular has the form

[qx] [k11 k12 k13] [ ∂T / ∂X]

[qy]= -[k21 k22 k23] * [∂T / ∂y] (1)

[qz] [k31 k32 k33] [∂T / ∂z],

Imagining that under a coordinate system ,the thermal conductivity tensor component k12,k13,k31,k32 not equal zero, and under this coorinate system ,along the x axis, we perform the nemd method.

For performing the nemd method in x axis,so the heat flux along y,z ,that is qy and qz euqal 0, also the temperature gradient in y,z direction ∂T / ∂y and ∂T / ∂z equal 0 as well.Bring these into the equation(1)

,then the equation (1)take the form :

[qx] [k11 k12 k13] [ ∂T / ∂X] [k11* ∂T / ∂X]

[0]= - [k21 k22 k23] * [0] =[k21*∂T / ∂X]
[0] [k31 k32 k33] [0] [k31* ∂T / ∂X],

for ∂T / ∂X not equal zero,so k21=k31=0,which is not consistent with the premise we make before.So does the thermal conductivity λ calculated by using nemd methods make any means if the system is anisotropic?

Sincerely,

XiaonengRan