One of the most common means of gold nanoparticle (AuNP) biofunctionalization involves the manipulation of precursor citrate-capped AuNPs via ligand displacement. However, the molecular-level structural characteristics of the citrate overlayer adsorbed at the aqueous Au interface at neutral pH remain largely unknown. Access to atomistic-scale details of these interfaces will contribute much needed insight into how AuNPs can be manipulated and exploited in aqueous solution. Here, the structures of such citrate overlayers adsorbed at the aqueous Au(111) interface at pH 7 are predicted and characterized using atomistic molecular dynamics simulations, for a range of citrate surface densities. We find that the overlayers are disordered in the surface density range considered, and that many of their key characteristics are invariant with surface density. In particular, we predict the overlayers to have 3-D, rather than 2-D, morphologies, with the anions closest to the gold surface being oriented with their carboxylate groups pointing away from the surface. We predict both striped and island morphologies for our overlayers, depending on the citrate surface density, and in all cases we find bare patches of the gold surface are present. Our simulations suggest that both citrategold adsorption and citratecounterion pairing contribute to the stability of these citrate overlayer morphologies. We also calculate the free energy of adsorption at the aqueous Au(111) interface of a single citrate molecule, and compare this with the corresponding value for a single arginine molecule. These findings enable us to predict the conditions under which ligand displacement of surface-adsorbed citrate by arginine may take place. Our findings represent the first steps toward elucidating a more elaborate, detailed atomistic-scale model relating to the biofunctionalization of citrate-capped AuNPs.