Daniel P.
Szekrényes
a,
Cyrille
Hamon
b,
Doru
Constantin
*bc and
András
Deák
*a
aCentre for Energy Research, 1121, Budapest, Hungary. E-mail: andras.deak@ek-cer.hu
bUniversité Paris-Saclay, CNRS, Laboratoire de Physique des Solides, 91405 Orsay, France
cInstitut Charles Sadron, CNRS and Université de Strasbourg, 67034 Strasbourg, France. E-mail: constantin@unistra.fr
First published on 20th October 2022
Gold nanoparticles coated with polyethylene glycol (PEG) are able to form clusters due to the collapse of the surface-grafted polymer chains when the temperature and ion concentration of the aqueous medium are increased. The chain collapse reduces the steric repulsion, leading to particle aggregation. In this work, we combine small angle X-ray scattering (SAXS) and visible light spectroscopy to elucidate the structure of the developing clusters. The structure derived from the SAXS measurements reveals a decrease in interparticle distance and drastic narrowing of its distribution in the cluster, indicating restricted particle mobility and displacement within the cluster. Surprisingly, instead of forming a large crystalline phase, the evolving clusters are composed of about a dozen particles. The experimental optical extinction spectra measured during cluster formation can be very well reproduced by optical simulations based on the SAXS-derived structural data.
We have shown earlier that the bulk phase clustering of PEG covered gold nanoparticles can be triggered by raising the temperature at elevated background ion concentration, leading to the formation of large and compact nanoparticle clusters with ordered internal structure.5 The driving force of the process is the developing soft-sphere type colloidal pair interaction energy profile as a result of surface grafted polymer chain collapse,13 and has been shown theoretically to be critical in the preparation of ordered nanoparticle clusters.14 Kinetic control of the clustering process can be also realized by controlling the concentration and temperature levels, leading to particle aggregates with different sizes.15 Our concept has been confirmed by others, preparing 2D and 3D aggregates of PEGylated gold nanoparticles, indicating a hexagonal (FCC) arrangement of the particles in the nanoparticle solid.4,16–18 It is important to note that this clustering approach differs substantially from creating ab ovo unstable systems or aggregation schemes relying on (macro)molecular linkers,8,19 but belongs to the broader family of solvent quality tuning based aggregation approaches.6,7,20
In this work, we aimed to investigate the development of nanoparticle clusters, where the driving force for clustering is rather moderate, favoring the formation of smaller particle aggregates. We investigated the assembly of PEGylated nanoparticles by a combined SAXS and visible light spectroscopy approach in order to consistently link these independent measurement results to the structural changes. The measurements were interpreted by theoretical calculations both of the structure factor and of the extinction spectra. We showed that SAXS provides a detailed picture of the cluster structure and can capture subtle structural changes as a function of temperature. By full-curve fitting the measured SAXS data, we found compact clusters composed of only a moderate number of particles, instead of the formation of a large crystalline nanoparticle solid phase. Optical extinction measurements to follow the clustering process of gold nanoparticles are far easier to perform, but they are less sensitive to fine details of the cluster structure as aggregation-related plasmon coupling dominates the observed spectral changes. Realistic agreement between simulated and measured optical spectra can be obtained, however, if the SAXS derived structural data is used as input for the simulations and the surface damping of the localized nanoparticle surface plasmons is also taken into account. Besides demonstrating the advantages in combining SAXS and spectroscopy approaches for the study of particle clustering, the results also underline that special care is required when only simple ensemble optical measurements are employed to monitor the clustering of gold nanoparticles – an approach that is broadly applied in various detection schemes.
The spherical nanoparticles were prepared by the traditional Turkevich method,21 while their PEGylation was carried out as reported earlier.5 Visible light extinction spectroscopy was realized using a fiber coupled stabilized tungsten halogen light source (Thorlabs SL201L) and spectrometer (Thorlabs CCS200) at different temperatures in a cuvette holder (Thorlabs CVH100) with laboratory built external temperature control. Three different temperature values (35 °C, 45 °C and 55 °C) and two different potassium sulfate concentrations (0.125 M and 0.25 M) were considered. The gold concentration is always 7.2 nM, corresponding to a volume fraction φP = 10−5.
For recording the optical extinction spectra, a custom written Labview software was used. Before recording the spectra, the temperature of both the nanoparticle sol and the salt solution was set to the desired value and the dilution of the salt solution upon injection was taken into account. The temperature was kept constant during the measurement. Optical simulations of the nanoparticles’ and nanoparticle clusters’ extinction spectra were performed by the MNPBEM MATLAB module,22 using the dielectric function of gold given by Olmon et al.23 This dielectric function was modified to include an interfacial damping term, which accounts for surface scattering of the conduction electrons and for chemical interface damping.24 The simulation algorithm accounts for retardation and the effective refractive index of the embedding medium was set to 1.467 due to the presence of PEG. The particles and the particle assemblies were illuminated along the x, y and z principal axes each with two perpendicular polarizations and the extinction spectrum was obtained by averaging over these 6 excitations. Ten different, randomized nanoparticle cluster structures were simulated and averaged in each case, generated in the same way as for the SAXS evaluation (see ESI†).
SAXS experiments were carried out using the SWING beamline at SOLEIL synchrotron facility, at an energy of 16 keV. The sample-detector distance was set to 6.223 m. The scattered intensity as a function of q vector modulus was obtained by radially averaging the measured scattering pattern using the FOXTROT software developed at the beamline. This SAXS signal (I(q)) can be written in the form of I(q) = nV2Δρ2P(q)S(q), where V is the volume of the particle and Δρ is the difference in scattering length density between the particle and the surrounding medium. The form factor, denoted by P(q), normalized such that P(q → 0) = 1, is due to the size and the shape of the particles, while the structure factor S(q) provides information on their configuration. To determine the experimental structure factor, the SAXS signal of assembled nanospheres was divided by a polydisperse sphere model fit to the initial (well-dispersed) sample: I0(q) S(q) = I(q)/Imodel(q). As shown in the ESI,†Imodel(q) is very close to I0(q), confirming that interparticle interactions can be neglected at room temperature. Note also that, due to this polydispersity, S(q) is strictly speaking an effective structure factor and is only equal to the true structure factor in the decoupling approximation: 〈V2P(q)S(q)〉 = 〈V2P(q)〉〈S(q)〉: considering the low relative polydispersity (about 10%), we neglect this distinction in the following.
Combined with the gold concentration, this parameter yields a number concentration of particles nP = 6.5 × 10−9 nm−3. The prepared samples are kinetically stable, as in addition to the electric double layer repulsion the end-grafted polymer chains also counteracts the van der Waals attraction. When the PEGylated particles are placed in a salt solution, however, their aggregation can be induced at elevated temperature due to the collapse of the surface-grafted polymer chains and the decrease of the associated steric repulsion between particles (Scheme 1).5,18 Depending on the resulting colloidal interaction profile, destabilization can lead either to relatively slow nanoparticle cluster formation, or to rapid aggregation,25 similar to the results obtained on ‘hairy’ and plain silica nanoparticles.26 When a soft-sphere type interaction curve is present due to a strong short range repulsion (due to PEG in the present case) with only a small (few kT) attractive potential minimum, the particles will assemble into fairly ordered colloidal clusters.14 Larger attractive potentials on the other hand results in fast aggregation rates and the diffusion limited growth produces fractal-like aggregates.5,27
For the PEGylated particles used in the present study, the decreasing stability of the particles with increasing temperature can be clearly detected in their time-dependent visible extinction spectra (Fig. 1). At the lowest investigated salt concentration and temperature (0.125 M; 35 °C) the sample remains stable and no change in the extinction maximum nor in the width of the peak can be observed. As the temperature is increased to 45 or 55 °C, however, a small redshift is observed and the extinction measured in the visible wavelength range decreases continuously. As the extinction value around 400 nm is dominated by interband transition instead of the localized plasmon resonance, its value correlates with the amount of Au0 in the light path.28 Hence, this observed decrease already indicates that some particles leave the detection volume, that is, the samples sediment to some extent. The effect of sedimentation can be separated from the changes induced by plasmon coupling by performing a renormalization of the extinction spectra at 400 nm. The renormalized spectra (for details see ESI,† Fig. S2) at 0.125 M for the 45 and 55 °C case show a small redshift and (especially for the 55 °C case) peak broadening can be observed, with slightly increasing extinction in the near infrared region. This plasmon-coupling-related effect confirms that the nanoparticles approach each other. On the other hand, the extent of the near-field coupling is small, indicating a rather large coordination distance. More pronounced spectral changes can be observed at higher, 0.25 M salt concentration. There is already a clear indication of cluster formation at 35 °C and the NIR extinction becomes dominant at higher temperature. At 55 °C a qualitatively different behavior is observed, with the appearance of a broad extinction peak in the NIR wavelength region, typical for the case of aggregation due to stronger interparticle attraction.5,11 It appears that, even though clustering takes place at different concentration/temperature combinations, there is a fundamental difference in the extent of clustering and/or in the resulting structure for the 0.25 M; 55 °C case. This is also indicated by SEM images taken from drop-casted samples after 60 minutes, which show that smaller clusters are formed at lower temperatures and ion concentrations (Fig. 2). Note, however, that the structure of dried-in SEM samples might differ from that of the as-formed cluster structure in the liquid phase. The spectroscopy measurements also indicate that after 60 minutes the clustering process comes to a hold, and no further aggregation related spectral shifts can be anticipated, except for sedimentation of large aggregates.
![]() | ||
Fig. 2 Scanning electron microscopy images of drop-casted samples taken after 60 minutes clustering time at 45 °C and 0.125 M (a) and 55 °C 0.25 M (b). |
SAXS investigation of the same systems was carried out to obtain more direct information about the structure of the evolving nanoparticle clusters. A locally compact nanosphere arrangement was assumed when evaluating the SAXS data,29 locally similar to the FCC structure of the extended nanoparticle aggregates obtained with large (6 kDa) PEG grafted nanoparticles.16 Cluster configurations are generated starting from an initial arrangement with equal interparticle distances, by adding along the three space directions random position fluctuations drawn from a Gaussian distribution with a given standard deviation. The structure factor of these objects is computed using the Debye equation.30 Our approach is similar to Hansen's Monte Carlo method.31
The structure factor depends sensitively on the average separation distance between the particles, its distribution, the number of spheres in the cluster as well as on the fraction of free vs. clustered nanoparticles (see ESI,† Fig. S3 for details). It is noteworthy that the calculated structure factors show the average of 1000 individual configurations, the corresponding particle–particle separation distribution functions are shown as insets in Fig. 3. The main parameters obtained from the SAXS measurements are listed in Table 1; particle surface separation refers to gold particle surface-to-surface distance, α is the ratio of clustered to all particles. The time evolution of the SAXS measurements (ESI,† S4 and S5.) corresponds to the spectral changes observed using visible light spectroscopy, that is, cluster formation proceeds rather slowly on the hour scale, and the clusters do not completely disintegrate when the system is cooled back to 25 °C.
Particle surface separation (Å) | Fraction of clustered particles – α | |
---|---|---|
0.125 M; 45 °C | 27.9 ± 13.3 | 0.63 |
0.125 M; 55 °C | 27.5 ± 13.8 | 0.78 |
0.25 M; 45 °C | 26.0 ± 12.4 | 0.83 |
0.25 M; 55 °C | 18.0 ± 1.7 | 0.96 |
The structure factor calculations provide good agreement with the measurements in Fig. 3(a)–(c), whereas for the sample measured at 0.25 M and 55 °C there is a clear difference at the peak around q = 0.09. The experimental spectrum exhibits additional features around q = 0.04 and q = 0.08, not captured by the calculation. They might originate from the presence of more complex cluster structures. It is important to emphasize that the best agreement with the experimental data was obtained assuming N = 13 particles in the cluster, regardless of salt concentration and temperature (ESI,† Fig. S6 and S7 for the demonstration). This also suggests that, in spite of the obvious nanoparticle association, the system is kinetically trapped in the few-particle cluster state, without developing huge particle aggregates. This conclusion is coherent with the results of the time-dependent spectroscopy measurements shown above. Even for the highest investigated salt concentration and temperature (0.25 M and 55 °C), the loose coordination of the smaller clusters can be anticipated, without leading to very large, compact aggregates. This is in agreement with the spectroscopic data (Fig. 1) and the SEM images. The average gap between the particles in the clusters shows generally a decreasing trend as the ion concentration and the temperature increase (Table 1). This is in qualitative agreement with recent results involving gold nanoparticles grafted with larger polymer chains.17 More importantly, there is a fundamental difference in the particle separation distribution. While for (a)–(c) the gap between the particles covers a fairly broad range, in the case of (d) there is a clear reduction of the particle separation and a pronounced narrowing of the gap distribution (also see ESI,† Fig. S8). We attribute this difference to the collapse of the surface-grafted PEG chains, which results in a deeper potential well in terms of colloidal interactions due to reduction of the steric repulsion range,5,15 and the attractive forces effectively lock the particles in place as the cluster is formed. This also explains the trend observed during the visible light extinction measurements: as long as the PEG only shrinks, the particles experience only moderate attraction and are separated from each other due to the still effective steric repulsion, so both the particle coordination within the cluster and the plasmon coupling will be rather weak. When PEG collapses, on the other hand, it leads to aggregation with small interparticle separations and hence to pronounced plasmon coupling. This conclusion is also supported by the calculated nanoparticle pair interaction potentials (ESI,† Fig. S9). An initial PEG brush thickness of 1.85 nm can be inferred from dynamic light scattering experiments.5 As a rough estimate, the SAXS-based surface separation values listed in Table 1 divided by a factor of 2 indicate ca. 25% smaller values (ca. 50% in the 0.25 M; 55 °C case); the separation data in Table 1 were used as input to obtain the steric interaction contribution. As expected, at these high ion concentrations the electric double layer interaction is basically absent with Debye lengths well below 1 nm (ESI,† Fig. S9), the interaction is dominated by the balance of the attractive dispersion and repulsive steric interactions.5
While there is a steady decrease of the interaction energy minimum, for the 0.25 M; 55 °C sample its depth drops well below −5 kT. This also means that while the particle separation is reduced, the strong interparticle attraction impedes reorganization of the particles within the clusters, enabling the growth of non-centrosymmetric clusters.14 The observed trend of α is also in line with the derived structural data and the colloidal interaction calculations. It indicates that the fraction of nanoparticles incorporated in a cluster gradually increases with ion concentration and temperature.
To confirm this scenario, data derived from the SAXS measurements was used as input for optical simulations and the clusters’ optical extinction spectra were compared with the experimental ones. As the particles used in this study are well below the mean free electron path length (ca. 37.7 nm for gold),32 the conduction electrons will scatter at the particle surface. Additionally, PEG is attached to the particles’ surface via a thiol moiety, resulting in an additional chemical surface damping (CID) of the plasmons. For spherical particles, the combined effect of surface scattering and chemical interface damping can be given as the interfacial damping Γint = AvF/R, where vF is the Fermi velocity (1.4 nm fs−1)33 and R is the particle radius in nanometers. A is the combined surface and CID – or interfacial – parameter, which is of order unity and can be taken into account in simulations by modifying the bulk dielectric function of gold.24 As shown in Fig. 4(a), the increase of the interfacial damping contribution results in a blueshift and broadening of the individual particle's resonance. While in theory the size of the nanoparticles as well as the polydispersity can affect both the position and width of the dipole resonance, for the given particle size their influence is negligible (ESI,† Fig. S10). As a consequence, the measured ensemble optical extinction spectrum can be compared to a particle's calculated spectrum. Following the approach of Berciaud et al.,34 setting the interfacial (combined surface and CID) parameter to A = 0.35 yields very good agreement with the experimental spectrum (Fig. 4(b)). Our A parameter is somewhat larger than the previously reported value (A = 0.25) for similarly sized gold nanospheres,34 which is the result of chemical linkage between the gold surface and the thiolated PEG in the present case, leading to a larger CID contribution,35 close to the value obtained for CTAB-capped gold nanoparticles (A ≈ 0.35).36,37
The good agreement of the measured and calculated spectra allows the calculation of the nanoparticle clusters’ spectra, with clusters constructed based on the SAXS-derived data. The cluster extinction spectrum simulation not only includes the separation distance and distribution values, but the fraction of particles incorporated in the clusters as well. This latter is done by dividing first the calculated cluster spectrum with the corresponding number of particles, yielding the average spectral contribution of a single clustered particle. The final simulated spectrum is then calculated as the sum of the spectral contributions from clustered and free particles with weighting factors α and 1 − α for the clustered and free particles, respectively, where the value of α is taken from the corresponding SAXS evaluation (Table 1). The comparison of the starting and final (after 60 minutes) measurement with the optical calculations is shown in Fig. 5. It is noteworthy that we achieve excellent agreement for (a)–(c), implying that the optical calculation based on the SAXS-derived structure effectively captures the small changes occurring during clustering.
A similar spectral evolution for clustered gold nanoparticles was observed earlier for 3-mercatopropionate-induced clustering in the presence of PEG.19 Note, however, that for the system indicating PEG collapse (0.25 M; 55 °C – Fig. 5(d)), the calculation fails to reproduce the experimental data: the emerging peak in the NIR wavelength region and the extreme broadening are not reproduced in the optical simulation. The measured spectrum for (d) indicates that, as a result of strong interparticle attraction, elongated structures (e.g. chain-like aggregates) are also present, and extinction from these optically dominant modes determines the measured spectra.4,11 Note that when the nanoparticle design limits the growth of the clusters (e.g. when Janus-type gold nanospheres are used where the clustering is self-limited due to anisotropic particle structure or encapsulation of the small clusters in a polymer shell), the optical spectrum does not show this extreme broadening.9,38
Footnote |
† Electronic supplementary information (ESI) available: SAXS measurement evaluation details, uncorrected experimental optical extinction spectra, scanning electron microscopy images of drop-casted samples, calculated structure factors and derived particle separation distributions, time evolution of the experimental structure factors (PDF). See DOI: https://doi.org/10.1039/d2sm01257j |
This journal is © The Royal Society of Chemistry 2022 |