Hello everyone,

I am trying to use `rNEMD`

method (Muller-Plathe) to get thermal conductivity for my system. However, I am not sure after getting the results, how should I convert it to W/(m.k) units.

In my simulations, I am using `real`

units, so the involved factors are in the following unit:

```
Q: KCal/mol
t: fs (10^-15 s)
Z: Angstrom (10^-10 m)
```

So, I want to convert thermal conductivity from

```
(KCal/mol)/(fs.Ang.k) to (W/m.k)
```

I have tried the same procedure as mentioned in this link. However, I failed to get reasonable results. Can anybody help with this?

These are the results that I got from my simulation:

```
Kinetic energy: 135997 KCal/mol
Simulation time: 10^6 fs
A: 59 (Ang)*59 (Ang)
dT/dz: 8.366 (k/Ang)
No. of atoms: 15660
```

Also, here is the temperature profile that I got from running the system.

TC_2.xlsx (16.4 KB)

I really appreciate any help in this regard,

Cheers,

Mahsheed