Gold nanostructures have been studied extensively for use in photothermal therapies and other targeted biomedical applications. While it is known that particle morphology can influence thermal conductivity, there is no consensus yet on how heat transfer is affected by nanostructure curvature in more complex systems. Nanostructure viability for biomedical applications is also highly dependent on surface functionalization. Citrate is a protecting agent thought to work via creation of an electrical double-layer when adsorbed onto gold surfaces. Since it is also a common ion in physiological settings, it has shown promise as a biocompatible ligand. Molecular dynamics simulations of nanostructures are a useful method for recreating and explaining experimentally observed phenomena. The structures and binding of proteins and citrate overlayers on gold nanoparticles have both been studied via molecular simulation. In a re-cent set of papers, Chong et al. studied the different binding motifs of citrate in various protonation states on the surfaces of gold particles. Perfilieva et al. stressed the importance of including gold polarizability when simulating with charged ligands such as citrate. They determined that the nature of citrate adsorption is influenced by whether the gold surface is polarizable. The polarizable force fields previously used to study citrate capped nanoparticles included a rigid-rod dipole model in which a freely rotating dipole was attached to each gold atom. While the gold atoms were held fixed in space, the dipoles rotated when a charged species approached the gold surface. This is an elegant model, but the lack of atomic motion in the gold surfaces prevents simulation of vibrational contributions to thermal transport. Alternatively, the polarizable metal potential developed by Bhattarai et al. gives each metal atom an additional dynamic partial charge variable that represents fluctuations in local valence density. This partial charge is then propagated along with the nuclear degrees of freedom, allowing the charge distribution in the metal to respond to the local environment. This density readjusting embedded atom model (DR-EAM) is a superset of EAM-style potentials with the addition of Coulomb and self-potential terms. In a study of the Au(111)-water interface, polarization of the gold surface using DR-EAM increased thermal conductance at the interface relative to EAM. This increase was likely due to an alteration in water ordering at the gold surface, as well as the image charge and image dipole interactions present in the DR-EAM potential. The highly charged citrate anion is a natural application of DR-EAM, as the charges are likely to alter the surface behavior of gold when modeled using a polarizable potential.