Gold nanoparticles
(GNPs) are versatile materials suitable for
various biological applications due to their tunable surface properties,
but structure–function relationships between specific GNP components
and GNP behavior are largely lacking. In this work, atomistic molecular
dynamics simulations were used to study the influence of gold core
morphology, size, and ligand length on the structure of uniformly
nonpolar alkanethiol monolayer-protected GNPs in water. By use of
a generalized system preparation workflow, three gold core models
were selected for this study: (1) a uniformly spherical hollow gold
core, (2) a spherical gold core cut from a bulk gold lattice, and
(3) a faceted gold core obtained from variance-constrained semigrand-canonical
simulations. Independent of the gold core morphology, we found that
long alkanethiol ligands exhibit increased ligand order and form quasi-crystalline
domains, or bundles, in which ligands orient in the same direction,
leading to asymmetric monolayer structures. Faceted gold cores promote
the formation of ligand bundles for short ligand lengths, but the
influence of the gold core is diminished for long ligand lengths.
We used a clustering algorithm to identify the subpopulation of bundled
ligands and found that bundling leads to heterogeneous surface properties,
whereby bundled ligands have a higher ligand order and lower surface
area accessible for solvent interactions. These findings demonstrate
the importance of GNP compositional features on monolayer structure,
which could impact GNPs interactions with other molecules.