“…Strictly speaking, the matrices −ω 2 M eα and −ω 2 M Iα should respectively be replaced by − i ωB eα − ω 2 M eα and − i ωB Iα −ω 2 M Iα in order to account for nonzero dissipation, though the dissipation matrices B eα and B Iα may be taken to be infinitesimal; also, once again, the diagonal blocks K Iαα entering Ẑ(0) α should actually include the effects of couplings to nuclear oscillators in other bodies as are present in the offdiagonal blocks K Iαβ for all β = α. With details explained in, 61,63 the RMB oscillator matrix parameters Q eα , M eα , M Iα , K eα , and K Iαα (the latter initially excluding couplings to nuclear oscillators in other bodies) along with the equilibrium atomic positions are all computed using density functional theory (DFT) for each body in isolation, while the matrices B eα and B Iα are assigned phenomenological values. These 2 × 2 block matrices can then be used in place of Ẑ(0) α and ∆ Ẑαβ in the formula for two components with general couplings (10) to find the combined heat transfer including PCHT and RHT: the couplings among valence electronic and nuclear DOFs through EM fields means that PCHT and RHT contributions are not separable, but in fact affect each other.…”