Electric Potential of Citrate-Capped Gold Nanoparticles Is Affected by Poly(allylamine hydrochloride) and Salt Concentration

The structure near polyelectrolyte-coated gold nanoparticles (AuNPs) is of significant interest because of the increased use of AuNPs in technological applications and the possibility that the acquisition of polyelectrolytes can lead to novel chemistry in downstream environments. We use all-atom molecular dynamics (MD) simulations to reveal the electric potential around citrate-capped gold nanoparticles (cit-AuNPs) and poly(allylamine hydrochloride) (PAH)-wrapped cit-AuNP (PAH-AuNP). We focus on the effects of the overall ionic strength and the shape of the electric potential. The ionic number distributions for both cit-AuNP and PAH-AuNP are calculated using MD simulations at varying salt concentrations (0, 0.001, 0.005, 0.01, 0.05, 0.1, and 0.2 M NaCl). The net charge distribution (Z(r)) around the nanoparticle is determined from the ionic number distribution observed in the simulations and allows for the calculation of the electric potential (ϕ(r)). We find that the magnitude of ϕ(r) decreases with increasing salt concentration and upon wrapping by PAH. Using a hydrodynamic radius (RH) estimated from the literature and fits to the Debye-Hü̈ckel expression, we found and report the ζ potential for both cit-AuNP and PAH-AuNP at varying salt concentrations. For example, at 0.001 M NaCl, MD simulations suggest that ζ = -25.5 mV for cit-AuNP. Upon wrapping of cit-AuNP by one PAH chain, the resulting PAH-AuNP exhibits a reduced ζ potential (ζ = -8.6 mV). We also compare our MD simulation results for ϕ(r) to the classic Poisson-Boltzmann equation (PBE) approximation and the well-known Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. We find agreement with the limiting regimes─with respect to surface charge, salt concentration and particle size─in which the assumptions of the PBE and DLVO theory are known to be satisfied.