We report the results of molecular dynamics simulations of the properties of a pseudo-atom model of dodecane thiol ligated 5-nm diameter gold nanoparticles (AuNP) in vacuum as a function of ligand coverage and particle separation in three state of aggregation, namely the isolated AuNP, isolated pair of AuNPs and a square assembly of AuNPs. Our calculations show that the ligand density along a radius emanating from the core of an isolated AuNP has the same gross features for all values of the coverage; it oscillates around a constant value up to a distance along the chain corresponding to the position of the fourth pseudo-atom, then smoothly decays to zero, reflecting both the restricted conformations of the chain near the core surface and the larger numbers of conformations available further from the core. Interaction between two AuNPs generates changes in the ligand distributions of each. We examine the structure and general shape of the ligand envelope as a function of the coverage and demonstrate that the equilibrium structure of the envelope and the deformation of that envelope generated by interaction between the NPs is coverage-dependent, so that the shape, depth and position of the minimum of the potential of mean force displays a systematic dependence on the ligand coverage.We propose an accurate analytical description of the calculated potential of mean force as a function of a set of parameters that scale linearly with the ligand coverage. Noting that the conformational freedom of the ligands implies that multiparticle induced deviations from additivity of the pair potential of mean force are likely important, we define and calculate an effective pair potential of mean force for a square configuration of particles; our definition contains, implicitly, both the three-and four-particle contributions to deviation from additivity. We find that the effective pair potential of mean force in this configuration has a different minimum and a different well depth than the isolated pair potential of mean force. Previous work has found that the three-particle contribution to deviation from additivity is