Transport properties such as interfacial thermal conductance and thermal conductivity can be easily studied using reverse non-equilibrium molecular dynamics (RNEMD) methods. A minimally perturbative velocity shearing and scaling RNEMD method (VSS-RNEMD) has been developed which allows for periodic scaling of particle velocities, and it has been extended for use in non-periodic simulations. VSS-RNEMD was most recently used to study the solvent thermal conductivity and interfacial thermal conductance in a simulation of gold nanoarrays capped with phenylethanethiol in toluene and dichloromethane solvents. However, polarization of the gold sur-face was not considered in this work as the systems did not contain charged species. Interfacial thermal conductance has also been shown to increase with increasing surface curvature, but citrate's effect on this trend has not previously been studied. This study aims to determine the effect of gold surface polarization and curvature on heat transport in systems containing sodium citrate.