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. Fa...