where ∆T is the temperature drop across the interface and J q is the heat flux, calculated using Equation . The temperature points near the heat source and heat sink regions and those adjacent to the interface are ignored in a linear fit to the solid and liquid temperature profiles used to calculate the ITR. For ClayFF, 45,76 ∆T = 17.9 K, while for Savio et al. ∆T = 5.9 K. for Savio et al. The value calculated for the strongest potential is similar to that reported from previous NEMD simulations of the lower limit of the ITC (R k ∼ 6 m 2 K GW -1 ) for ionic liquid-graphene interfaces. The range of ITR values is also similar to that observed experimentally (R k = 6-20 m 2 K GW -1 ) for the interface between water and chemically functionalized aluminum and gold surfaces. The ITC was found to increase linearly with ε F e-C and ε F e-H , as shown in Figure in the Supplementary Information. This linear relationship between G k and ε has been reported in previous NEMD studies for simple LJ systems and the Si/Water interface. We do not observe an exponential scaling of G k with ε F e-C and ε F e-H , as has been observed in some NEMD simulations at low interaction strengths. 98,99 We conducted further analysis to explore the temperature dependency of R k . Simulations