Probing the Limits of Conventional Extended X-ray Absorption Fine Structure Analysis Using Thiolated Gold Nanoparticles

We present a method for quantifying the accuracy of extended X-ray absorption fine structure (EXAFS) fitting models. As a test system, we consider the structure of bare Au-147 nanoparticles as well as particles bound with thiol ligands, which are used to systematically vary disorder in the atomic structure of the nanoparticles. The accuracy of the fitting model is determined by comparing two distributions of bond lengths: (1) a direct average over a molecular dynamics (MD) trajectory using forces and energies from density functional theory (DFT) and (2) a fit to the theoretical EXAFS spectra generated from that same trajectory. Both harmonic and quasi-harmonic EXAFS fitting models are used to characterize the first-shell Au-Au bond length distribution. The harmonic model is found to significantly underestimate the coordination number, disorder, and bond length. The quasi-harmonic model, which includes the third cumulant of the first-shell bond length distribution, yields accurate bond lengths, but incorrectly predicts a decrease in particle size and little change in the disorder with increasing thiol ligands. A direct analysis of the MD data shows that the particle surfaces become much more disordered with ligand binding, and the high disorder is incorrectly interpreted by the EXAFS fitting models. Our DFT calculations compare well with experimental EXAFS measurements of Au nanoparticles, synthesized using a dendrimer encapsulation technique, showing that systematic errors in EXAFS fitting models apply to nanoparticles 1-2 nm in size. Finally we show that a combination of experimental EXAFS analysis with candidate models from DFT is a promising strategy for a more accurate determination of nanoparticle structures.