A combined 3D and 2D light scattering study on aqueous colloidal model systems with tunable interactions

In this article we report on the synthesis and characterization of a system of colloidal spheres suspended in an aqueous solvent which can be refractive index-matched, thus allowing for investigations of the particle near-wall dynamics by evanescent wave dynamic light scattering at concentrations up to the isotropic to ordered transition and beyond. The particles are synthesized by copolymerization of a fluorinated acrylic ester monomer with a polyethylene-glycol (PEG) oligomer by surfactant free emulsion polymerization. Static and dynamic light scattering experiments in combination with cryo transmission electron microscopy reveal that the particles have a core shell structure with a significant enrichment of the PEG chains on the particles surface. In index-matching DMSO/water suspensions the particles arrange in an ordered phase at volume fraction above 7%, if no additional electrolyte is present. The near-wall dynamics at low volume fraction are quantitatively described by the combination of electrostatic repulsion and hydrodynamic interaction between the particles and the wall. At volume fractions close to the isotropic to ordered transition, the near-wall dynamics are more complex and qualitatively reminiscent of the behaviour which was observed in hard sphere suspensions at high concentrations.


Introduction
The properties of colloidal particles at or close to interfaces between an aqueous and a non-aqueous phase have acquired growing attention from the scientific community and industry during recent years. 1 The spectrum of topics spans fundamental scientific subjects such as particle stabilized emulsions, 2 protein stability and adsorption 3 or self-assembly at water-oil interfaces 4 for the fabrication of functional nano-films, with industrial applications such as corrosion protection 5 or water-borne 'green' coatings, 6 to name a few examples.While there is a plethora of publications considering static properties such as interaction forces, structure formation and phase behaviour of particles at interfaces, [7][8][9] there are comparatively few contributions reporting on experimental investigations of particle dynamics in the ultimate vicinity of interfaces, [10][11][12][13][14][15][16][17][18] despite the fact that the dynamics of colloidal particles at interfaces is an important feature, especially when the dimension of confining geometries are approaching those of the particles, as in the case of microfluidics, 19 emulsion stabilization 20,21 as well as many biological processes. 22,23Interfaces between a solid and a liquid phase are of particular technological importance e.g. in coating processes or bio-film formation.In these cases, the motion of the particles close to the wall is slowed down due to hydrodynamic interaction and this hindrance differs for particles moving parallel and perpendicular to the wall. 24,25The interplay between static particle wall interaction and the slowing-down caused by the so-called wall drag effect has significant implications in potential applications.Geometrical constraints were used to quench Brownian rotational motion of spherical Janus particles providing a means to steer their trajectory. 26Furthermore, the trajectories of active swimmers, driven by flagellar motion, can be directed by tuning their interaction with a surface, and spatial confinements can be used to accumulate active swimmers. 27For a targeted exploitation of the combined effects of static interaction, the wall drag effect and particle concentration, a fundamental understanding of their mutual interplay is essential.In this paper we report on the synthesis and characterization of an aqueous system, which allows a systematic variation of the static interactions and at the same time enables the study of the near-wall dynamics of particles at arbitrary concentrations.
For experimental purposes, interfaces between transparent solids and a liquid phase are convenient, since in this case, the near-wall dynamics of colloidal particles can be investigated by evanescent wave dynamic light scattering (EWDLS).This technique This journal is © The Royal Society of Chemistry 2016 was first devised by Lan et al. in 1986, 28 and has been advanced constantly since then.The principle is as follows: an evanescent wave is created as an illumination source when a laser beam is totally reflected off a glass-solution interface; particles, which are located within this illumination profile will scatter light, which is then collected by a detecting unit and passed down to a correlator to generate the intensity correlation function.The penetration depth of the evanescent wave can be tuned by varying the incident angle of the laser beam.Being a light scattering technique in essence, EWDLS is most suited to the measurement of particles within a radius range around 100 nm, and at the same time provides good statistics.Modern instrumentation has made it possible to independently vary the scattering vector components parallel, Q J , and normal, Q > , to the interface and hence to distinguish between the anisotropic diffusivity of colloids parallel and perpendicular to the interface experimentally. 15EWDLS has also been employed to study the rotational diffusion of optically anisotropic particles, 29 to probe the particle dynamics at liquid-liquid interfaces, 30 and to measure the near-wall velocity profile of a streaming suspension with a resolution of tens of nanometers. 31The near-wall dynamics of colloidal hard spheres dispersed in an organic solvent have been systematically studied over a broad range of volume fractions. 16,17In the latest endeavor, experimental data were thoroughly compared with theoretical predictions based on a virial approximation and simulation results 32 where excellent agreement has been achieved. 18This agreement between experiments and theory has provided a framework for the understanding of near-wall dynamics in general, and paved the way for the study of more complex systems.
Despite its high potential, EWDLS has not been employed to study the near-wall dynamics of particles in aqueous solution, except for very dilute suspensions.The main challenge is to find a suitable aqueous model system, which is monodisperse and can be refractive index-matched to the aqueous solvent to avoid multiple scattering.To this end, the fluorinated monomer 2,2,3,3,4,4,4-heptafluorobutyl methacrylate (FBMA) has been successfully used to synthesize monodisperse particles with a low refractive index via emulsion polymerization and the system was shown to form well-ordered crystalline colloidal arrays. 33Further effort has led to the production of FBMA particles with radii between 50 nm to 700 nm. 34However, such particles are susceptible to irreversible aggregation upon salt addition or under compression by centrifugation, as they are solely stabilized through electrostatic repulsion (ER) in aqueous solution.This would bring great difficulty into volume variation and interaction tuning, which is essential for the study of nearwall dynamics in dense suspensions.This situation calls for a more robust aqueous model system.The first attempt in this direction was made by Wiemann et al., who designed a coreshell system with a fluoroacrylate core and a poly(ethylene glycol) (PEG) surface layer, for the study of glass transition and re-entry behaviour. 35These so-called PEGylated particles are sterically stabilized in aqueous solution by the PEG layer, and were shown to be stable against compression up to a particle volume fraction of 69% and high salt concentrations.Continued effort has made it possible to synthesize particles with a single and well-defined PEG layer. 36Moreover, this system bridges fundamental research to life science application, as PEG is a biocompatible material and e.g.PEGylation has already been widely used for drug delivery applications. 37,38n this paper we show that PEGylated FBMA particles are an interesting model system for the study of near-wall dynamics in aqueous solution for several reasons.First of all, they fulfill the requirements of being monodisperse, iso-refractive with aqueous solvents and stable, which ensures reliable EWDLS measurements; secondly, they are highly versatile, as their interactions can be tuned by varying the PEG chain length, the solvent composition, and salt concentrations where the PEG-layer should ensure colloidal stability even at high salt concentration.We describe the synthesis of the particles by emulsion polymerization and their bulk characterization, by cryo transmission electron microscopy (cryo-TEM) and light scattering.Finally we show that the near-wall dynamics of this model system is in very good agreement with current theoretical predictions for low particle volume fractions.However, due to the long range particle-particle interactions and those between the particles and the wall, the effect of increased concentration on the particles' near-wall dynamics is much more pronounced than for suspensions of colloidal hard spheres.

Experimental
2.1 Materials and synthesis 2.1.1Materials.The fluorinated monomer 2,2,3,3,4,4,4heptafluorobutyl methacrylate (HFBMA, 97%) was purchased from Alfa Aesar.Prior to use, the inhibitor hydroquinone was removed by washing the monomer three times with a 5% NaOH (w/v) solution.The macromonomer methoxy poly(ethylene glycol) acrylate with a molecular weight of M n = 950 g mol À1 was purchased from Sigma-Aldrich.To remove the inhibitor from the PEG macromonomer, it was dissolved in water and passed through an inhibitor remover packed column (CAS, Sigma-Aldrich).Dimethylsulfoxide (DMSO) was purchased from Merck.The initiator, potassium peroxodisulfate (K 2 S 2 O 8 ), was obtained from Sigma-Aldrich.
2.1.2Synthesis.The PEGylated fluorinated particles were prepared by surfactant free emulsion polymerization, following the procedure given in ref. 36.In a typical synthesis, 1000 mL H 2 O are deoxygenated three times in a 3 L three-necked round bottom flask by firstly applying a vacuum and afterwards flooding with nitrogen.After deoxygenation, 34.9091 g (130 mmol) HFBMA and 69.7 mL of an acryl-PEG solution (10 g dissolved in 100 mL water; 6.97 g; 7.34 mmol) were added under inert gas countercurrent.The emulsion was heated to 70 1C under reflux cooling, and vigorous stirring (we used a magnetic stirrer bar which was driven by a motor with a nominal speed of 310 rpm).After the injection is finished the reaction mixture turns turbid with a bluish colour.The input adapter is removed and the reaction is allowed to proceed for twenty hours.After cooling, the mixture is filtered through a Bu ¨chner funnel.To remove any remaining starting material, the dispersion is dialyzed for 7 days against deionized water from a Milli-Q water purification system with a resistivity of 18.2 MO cm À1 , changing the water twice a day.The solid content of the dialyzed dispersion was determined gravimetrically by drying three 1 mL aliquots at 100 1C under reduced pressure, weighing the solid remainder and calculating the average solid content.
2.2 Bulk characterization methods 2.2.1 Static and dynamic light scattering.The light-scattering measurements were performed on a commercial instrument by ALV-Laservertriebsgesellschaft (Langen, Germany) which is equipped with a HeNe Laser as a light source, an automated attenuator device, a Perkin Elmer avalanche diode and an ALV-6000 multiple tau digital correlator.Time-averaged static scattering intensities were obtained by measuring count rates of the diode and correcting them for attenuator transmission, scattering volume and solvent background.Furthermore, the scattering curves were normalized, using toluene as a reference and dividing out the particle concentration.Static intensities were recorded in an angular range between 301 and 1201 with an increment of two degrees.Time auto correlation functions of the scattered intensity g 2 (Q,t) were recorded in a range of scattering angles 301 r y r 1201 in steps of ten degree.Here Q = 4pn S sin(y/2)/l 0 is the scattering vector with n S the refractive index of the solution and l 0 the vacuum wave length of the laser.To determine the diffusion constant of the particles, D 0 , and their polydispersity index, U, a cumulant analysis of the correlation functions was performed.The particles' hydrodynamic radius, R H , was calculated via the Stokes-Einstein relation, using the appropriate solvent viscosity.
2.2.2 Cryo transmission electron microscopy: cryo-TEM.Transmission cryo electron microscopy (cryo-TEM) was used to evaluate the morphology and size of the colloidal particles.3 mL of a solution was applied on a hydrophilic Quantifoil TEM grid.In order to obtain a thin vitreous ice layer, the excess of solution was blotted for 2.5 seconds with a filter paper and the grid was plunged into liquid ethane in a FEI Vitrobot.The specimen was maintained at a temperature of approximately À196 1C using a cryo holder.Bright field TEM images were acquired on a FEI Spirit TEM operated at 120 kV.Since conventional TEM images correspond to twodimensional (2D) projections of three-dimensional (3D) objects, electron tomography was performed to investigate the morphology of the colloidal particles.A series of 2D projections was acquired while tilting the specimen around an axis perpendicular to the electron beam over the range from À651 to 601 with an increment of 51.After alignment, the 2D images were reconstructed by the Simultaneous Iterative Reconstruction Technique (SIRT) algorithm in the Astra Toolbox. 393 Evanescent wave dynamic light scattering 2.3.1 EWDLS set-up.EWDLS experiments were performed with an instrument built in house, based on a triple axis diffractometer by Huber Diffraktionstechnik, Rimsting, Germany, which has been described in detail elsewhere.15 The setup is equipped with a frequency doubled Nd/Yag Laser (Excelsior; Spectra Physics) with a vacuum wavelength of l 0 = 532 nm and a nominal power output of 300 mW as a light source.Scattered light is collected with an optical enhancer system by ALV Lasververtriebsgesellschaft, Langen, Germany, and fed into two avalanche photo diodes by Perkin Elmer via an ALV fiber splitter.The TTL signals of the diodes are cross-correlated in time using an ALV-6000 multiple tau correlator.The scattering geometry and the definition of the scattering vector and its components parallel and normal to the interface are sketched in Fig. 1.The sample cell (custom-made by Hellma GmbH, Mu ¨llheim, Germany) consists of a hemispherical lens as the bottom part, made of SF10 glass, with an index of refraction n 1 = 1.736 at l 0 = 532 nm.The suspension of FBMA particles is contained in a hemispherical dome which covers the lens from the top.The incident beam is reflected off the interface between the semispherical lens and the suspension, which is considered to be a flat wall.The resulting evanescent wave, which extends into the particle suspension is used as the illumination for the scattering experiment.It has a wave vector k e pointing in the direction of the reflected beam's projection onto the reflecting interface.The scattering vector is given by Q = k s À k e , where the scattered light wave vector, k s , is set by the two angles y and a r describing the position of the detecting unit.
The inverse penetration depth of the evanescent wave is determined by the angle of incidence a i as k=2 ¼ l 0 and the refractive index of the glass n 1 and the suspension n 2 , respectively.The magnitude of the scattering vector component parallel to the interface is given by 2 sin a r /l 0 .The two components can be varied independently of each other by changing y and a r .To distinguish between particle dynamics parallel and normal to the wall, series of correlation functions are recorded at fixed 2 EWDLS data analysis.The analysis of the scattered intensity time auto-correlation function g 2 (t) from EWDLS is more complicated than in conventional bulk dynamic light scattering (DLS), mainly for two reasons.The first major complication occurs from the fact that a simple quadratic Siegert relation between the measured g 2 (t) and the normalized correlation function of the scattered field g ˆ1(t), which is usually assumed in DLS, does not apply in most cases in EWDLS.The scattering signal from the sample will be accompanied almost inevitably by a static scattering signal originating from surface corrugations, which can lead to large variations of the amplitude of the experimental correlation functions, as is shown in the top part of Fig. 2. Therefore, a mixed homodyne/heterodyne detection scheme has to be taken into account and a generalized Siegert relation 15 has to be applied for the conversion from g 2 (t) to g ˆ1(t).Here, , with A the experimental intercept of g 2 (t).Furthermore, in many cases EWDLS intensity-auto correlation functions exhibit a very slow decay at large times, as shown in the top part of Fig. 2. The physical origin of this slow relaxation is not clear yet.While Garnier et al. conjecture that it is due to a slow reversible adsorption of the particles to the wall due to van der Waals attraction, 11 Lisicki et al. argue that it is also caused by the unavoidable stray-light from surface defects, which is scattered by colloids in the bulk of the suspension into the detector. 40Since these slow modes are in general well separated from the relaxation rates of interest, we approximate their contribution by an additional baseline B 1 to g ˆ1(t).Thus, to determine the initial slope G of g ˆ1(t), which is related to the dynamic properties of interest, we chose to fit the experimental correlation functions, using a non-linear least squares algorithm, to eqn (1), where g ˆ1(t) is modeled as a decaying single exponential function in time. ( According to eqn (1) and ( 2), B 1 is related to B 2 , the baseline of Consequently there are three fit parameters A, G and B 2 .Since an erroneous baseline value will lead to a deviation of G from the true relaxation rate, due to a normalization error, we fitted g 2 (t) repeatedly, starting with a number of data points, N p .After the convergence of a single fit, two data points at the long time end of g 2 (t) were removed, and the remaining points were fitted again with reduced N p .This procedure was repeated until N p o 20.With this technique it was possible to identify a limited range of N p 's where the best fitting parameters are essentially independent of the number of points.The G values determined in this range are considered to be the initial slope or the first cumulant of g ˆ1(t).Repeated measurements on the same sample show that G can be determined by this procedure with an uncertainty which is generally smaller than 10% as shown by the error bars displayed in Fig. 9.
2.3.3Near-wall dynamics of particles with long range interactions.The mobility of colloids close to a wall is smaller than in bulk and depends on their distance, z from the wall due to the particles' hydrodynamic interaction with the wall.The resulting slowing down of Brownian motion is more pronounced in the direction normal to the wall than parallel to the wall.This has been described by two position dependent friction coefficients l J (z) and l > (z). 24,25In the simplest case of hard spherical particles the initial slope of the field correlation function G is related to these friction coefficients by (3) where hD J,> i k refer to position-dependent diffusion coefficients averaged over the illumination profile of the evanescent wave 14 as The superscript k indicates that the averaged quantities depend slightly on the penetration depth of the evanescent wave.For particles, which are subject to a long-ranging interaction potential with the wall V(z), the Boltzmann-factor of the potential will enter into the average diffusion constants 15 as where b = 1/k B T is the inverse thermal energy.
3 Results and discussion 3.1 Particle characterization 3.1.1Synthesis and general properties.We prepared FBMA latex particles which are charge stabilized in aqueous media and which are additionally stabilized sterically by chemically attached PEG chains, according to a semi-batch emulsion polymerization procedure. 36Electrophoresis measurements using a Malvern Zetasizer 2000, show that the particles behave as charged spheres at low ionic strengths with a zeta potential of À30 mV.The density of particles was determined to be 1.5266 g mL À1 with a density meter (Anton Paar DMA 4500).Given a particle radius of R E 130 nm (see section on cryo-TEM and bulk light scattering) this density results a sedimentation length L s E 300R, indicating that sedimentation should not significantly influence the results of the EWDLS experiments which were performed to investigate particle near-wall dynamics.This is experimentally verified by the typical standard deviations of repeated EWDLS measurements, as shown by the error bars in Fig. 9.
The particles yield turbid suspensions, if dispersed in water at volume fractions above f \ 0.01, since the refractive index of the fluorinated polymer is higher than that of water.In order to conduct light scattering measurement at higher concentrations, the refractive index of the particles needs to be matched to the surrounding medium by adding a co-solvent.DMSO was chosen as co-solvent because of its miscibility with water in all proportions, and its very low toxicity.A series of DMSO/water solutions was prepared, with refractive indices ranging from 1.33 to 1.40.To keep the particle concentration constant in all cases, 0.2 mL of a dialyzed stock suspension were added to 0.8 mL DMSO/water mixtures of various compositions.The scattering intensities of the resulting samples were measured at a scattering angle of y = 901, as the scattering intensity at constant volume fraction is approximately proportional to the square of the refractive index difference between the solvent and the particles.At a volume ratio of water to DMSO, f(H 2 O) : f(DMSO) of 69 : 31, the scattering intensity shows a minimum, which corresponds to the smallest obtainable refractive index difference between particles and solvent.The refractive index at this solvent composition was determined to be n sol = 1.3847, with an ABBEMAT RXA 156 refractometer from Anton Paar.
To obtain index-matched suspensions for further scattering experiments, the aqueous suspension was centrifuged at low speed and the sediments re-dispersed in the index-matching solvent mixture, for three times.The resulting dispersion is completely transparent and used as stock solution.
3.1.2Cryo-TEM.The shape and size of the particles were observed by cryo-TEM.In Fig. 3 we show representative images of single particles and some clusters which formed due to the blotting.The interfaces between the voids of the clusters and the spheres appear blurred, which we speculate might be due to the PEG chains.A total of 132 spheres were investigated to determine the mean particle radius and the relative standard deviation, which were found to be R = 118 nm and s r = 0.05, respectively.
To confirm the spherical shape of the particles they were investigated by cryo-electron tomography.A series of 2D projections was acquired while tilting the specimen around an axis perpendicular to the electron beam (x-axis) over a range from À651 to +601 with a tilt increment of 51.Example images from this series are shown in Fig. 4, where the perfectly circular shape of the 2D projection can be observed for all tilt angles.This observation is supported by the calculated 3D reconstruction which is displayed in Fig. 5.
The elongation of the 3D-reconstruction, which is observed in the zx and xy planes, is an artifact, which is caused by the limited number of acquisitions and the so-called missing wedge problem.The missing wedge is caused by the inability to tilt the sample over AE901 range.This is related to the limited gap between the pole pieces of the objective lens and shadowing effects of the TEM grid at high angles. 41As the 3D reconstruction is   mathematically speaking an ill-defined inverse problem, multiple solutions correspond to the incomplete data set that was collected.This eventually results in a blurring/elongation of the reconstruction in the direction of the missing wedge.
3.1.3Static light scattering and phase behaviour.For a characterization of the particles in solution, combined static and dynamic light scattering experiments were performed.All the scattering data from bulk samples which will be discussed later were obtained from suspensions without any added salt.Therefore the salinity is about 10 À5 mol L À1 due to carbon dioxide adsorption, corresponding to a Debye screening length of l D E 95 nm.The time-averaged static scattering intensity, expressed in absolute intensities as the Rayleigh ratios R(Q), from three samples with low volume fractions is shown in Fig. 6.The data obtained from the samples with different volume fractions, f, in DMSO/water could be superimposed perfectly by dividing out the particle concentration, which indicates that at least up to f = 0.014 there is no discernible effect of the solution structure factor on the scattering curve.Intriguingly the shape of the scattering curve from the suspension in pure water differs significantly from the curves obtained from samples which were dispersed in the index-matching solvent.The latter two curves show a distinct minimum which could be misinterpreted to indicate a particle radius of about 200 nm.For homogeneous spheres the position of the first minimum in the particle form factor is related to the particle radius by RQ min E 4.49.Thus for particles with a radius of about 120 nm, as indicated by cryo-TEM, the first minimum should occur at Q min E 0.037 nm À1 , while here we find Q min E 0.023 nm À1 implying a particle radius significantly larger than observed in cryo-TEM.This is a strong indication that the particles are not optically homogeneous spheres but rather exhibit a variation of the scattering contrast with distance from the particle centre. 42,43We therefore fitted the static scattering data using the form factor of a core-shell particle with a core radius R c and an outer radius R which is given by where and V e = 4p(De cs R c 3 + De ss R 3 )/3 and De cs = e c À e s is the difference between the dielectric constants of the core and the shell while De ss = e s À e sol is the difference between the shell and the solvent with the corresponding refractive indexes n i ¼ ffiffiffi ffi e i p of the core, the shell and the solvent respectively.To account for the particles' polydispersity we applied a logarithmic normal distribution p(R c ) of the core radii assuming a constant shell thickness d shell = R À R c .With this the model function used for the fitting is PðQÞdR c (7)   with the wave vector of the incident beam, k I = 2pn S /l 0 , and the particles' average dielectric constant e = (R c 3 De cs À R 3 e s )/R 3 .The lines in Fig. 6 were obtained by simultaneously fitting the three sets of data with eqn (7), where the global parameters R c , R, s r were forced to be the same for all three sets of data, while the local parameters n c and n s were allowed to vary between the data sets.The resulting best fitting parameter values are collected in Table 1.
The geometrical parameters are in very good agreement with the information deduced from cryo-TEM.Further we find n c E 1.38 which is qualitatively in agreement the reported refractive index values for pure FBMA particles n FBMA E 1.37. 34The parameters associated with the particle shell are well in line with the picture of a thin PEG layer, which is swollen by the solvent.The molar mass of the PEG chain corresponds to an average degree of polymerization of N E 22, yielding an all-trans contour length of the chain l c E 8 nm and a radius of gyration R g E 0.8 nm if a length of about 0.13 nm per backbone atom 44,45 and a monomer segment length of approximately 0.38 nm is assumed.This implies that a shell consisting of PEG should have a thickness between approximately 1.5 nm and 8 nm in accordance with our findings from SLS hRi À hR c i = 5.9 nm.Assuming an ideal mixture, we can obtain an estimate for the PEG volume fraction in the shell as f PEG = (n s À n sol )/ (n PEG À n sol ).Applying n PEG E 1.358 46 this yields a polymer volume fraction of about 20% in both solvents.These findings strongly indicate that our particles exhibit a pronounced core shell structure with a significant enrichment of the PEG chains at the surface of an FBMA core, which was not reported before.
Generally, with increasing volume fraction of the particle suspension, a transition from an isotropic solution to a colloidal crystal is expected.Indeed, FBMA-particle suspensions with volume fractions f 4 10 are opalescent, indicating the formation of an   ordered structure.With lower volume fractions this colourful effect vanishes gradually and the determination of the transition volume fraction by optical inspection is rather ambiguous.Therefore we applied the Hansen-Verlet criterion to estimate the lowest volume fraction at which coexistence between an isotropic solution and an ordered structure occurs.According to this empirical rule, it is generally accepted that the maximum of the structure factor S(Q max ) = 2.85 at the freezing point of colloidal hard sphere suspensions, 47 while it slightly increases for particles interacting by electrostatic repulsion (ER), depending weakly on their surface charge density and the solvent salinity.9][50][51] In Fig. 7 we plotted the structure factors of suspensions with different volume fractions, obtained from the static scatting data according to SðQÞ % IðQ; fÞ where I(Q,f 0 ) was measured at f 0 = 1.4 Â 10 À2 at which the structure factor is unity at all Q to a very good approximation.At the given salinity S(Q max ) t 3 for the sample with f = 0.07, while the structure factor maximum is significantly larger than the threshold value for the sample with 8% volume fraction.Thus, at f = 0.08 there is definitely a coexistence between an isotropic solution and an ordered structure.
3.1.4Dynamic light scattering.Before conducting experiments on near-wall dynamics, it is important to check if and how the index-matching solvent influences the diffusion behaviour of the system in bulk other than by viscous friction.For this purpose, we performed dynamic light scattering experiments on dilute particle dispersions both in water and in the index-matching DMSO/water mixture.Fig. 8 shows a plot of the first cumulant G 1 vs. Q 2 for particles dispersed in pure water and in the mixed solvent, respectively.The diffusion coefficient of the particles is the slope of the observed linear relations, D 0 = 1630 nm 2 ms À1 in water and 746 nm 2 ms À1 in the mixed solvent.The hydrodynamic radius, calculated through the Stokes-Einstein relation is R H = 131 nm in both cases, when the correct values for the solvent viscosity are used, i.e.Z = 1 mPa s for water and 2.187 mPa s for the mixed solvent.This shows that the solvent exchange has no effect on the particles' bulk dynamic except for the increased viscous friction.Moreover, the linearity of G 1 versus Q 2 plots in both cases is a first hint to a narrow particle size distribution indicating the absence of aggregates in the dispersions.A more thorough estimate of the particle size distribution can be obtained from a cumulant analysis of the correlation functions.The polydispersity index is given by U = G 2 /G 1 2 where G 2 is the second cumulant.We find U = 0.03 for the particles in pure water and 0.09 in the solvent mixture.Assuming a logarithmic normal distribution of the particle radii these figures correspond to a relative standard deviation s r E 0.05 in the first case, which is in perfect agreement with the observed standard deviation from cryo-TEM.The relative standard deviation measured in the indexmatching solvent mixture is 0.1 again showing there is no significant aggregation of particles in this case.The hydrodynamic radius obtained from DLS is larger than the geometric radius obtained from cryo-TEM, which is R = 120 nm.This is explained by the fact that the two methods measure different moments of the size distribution.For a logarithmic normal distribution the ratio of both quantities is related to the standard deviation by R H /R = (1 + s 2 ) 5 which for the present system implies a relative standard deviation of the order of ten percent in good agreement with the values from the cumulant analysis of DLS data.

EWDLS measurements
To investigate the PEGylated FBMA particles near-wall dynamics at low volume fractions, we performed EWDLS measurements on suspensions with f = 0.01.As the colloids and the glass surface are negatively charged, the addition of electrolyte will make a difference to the particles' near-wall dynamics, by changing the range of the electrostatic repulsion (ER).To illustrate this influence, Fig. 9 shows values of G plotted versus the square of the parallel scattering vector component obtained during parallel scans, at two different salt concentrations.Here and in the following we show normalized quantities, where d = 2R H is the  This journal is © The Royal Society of Chemistry 2016 particle diameter, to facilitate comparison between the results obtained from different systems, in particular with a hard sphere system as shown in Fig. 12.The red symbols represent a sample with 10 mM NaCl, which corresponds to a Debye length of l D = 3 nm.In this situation the electrostatic repulsion (ER) is almost completely screened and the measured G agree within experimental error with the prediction for a hard sphere with a hydrodynamic radius of R H = 130 nm, represented by the black line in Fig. 9.The error bars assigned to these data, which represent the standard deviation of three repeated measurements, are representative for all EWDLS data discussed in the following.For the sake of clarity we will not display error bars in the following graphs.The data represented by the blue symbols were obtained from the sample without extra salt addition, which corresponds to an ionic strength B10 À5 M and a Debye length of l D = 95 nm.In this suspension, the particles move significantly faster than those dispersed at a high salt concentration.This can be explained intuitively, since the long range electrostatic repulsion (ER) will effectively push the particles away from the interface thereby reducing the hydrodynamic wall drag effect.To quantitatively analyze these data, we calculated the G vs. Q J 2 dependence, by numerically integrating eqn ( 5) into which we introduced the electrostatic interaction potential which may be approximated within the Debye-Hu ¨ckel-Derjaguin limit by The amplitude B ER scales with particle radius and depends on the surface charge density of particles and the wall.In the model calculations, this amplitude is the only adjustable parameter and we could achieve good agreement between all experimental data discussed in this paper and the predictions with a fixed value of B ER = 2k B T. In Fig. 9 and 10, the lines labeled ER represent predictions for l D = 95 nm (blue), and for l D = 3 nm (red).The black line shows the hard sphere prediction as a reference, where no electrostatic interaction is taken into account.The systematic deviation of the experimental data at high Q J is due to the fact that the primary laser beam is reflected upon exiting of the semi-spherical lens, by that giving rise to a second scattering event with an in-plane angle y r = 180 À y.This effect is frequently observed, even in bulk experiments, with index-matched particles.For the case of EWDLS on hard sphere samples it has been discussed in detail in a previous paper. 18As a consequence of this back reflection, data obtained from EWDLS experiments at y 4 901 are often not reliable, which in the present case corresponds to (Q J d) 2 \ 30, given a r = 601 and the refractive index of the solution n S E 1.38.
The agreement between experiments and predictions is less satisfying in the case of the Q > -scans, as shown in Fig. 10, since the scatter of the experimental data is significantly larger, especially at long Debye lengths.At high salt concentrations the experimental data are reasonably well represented by the prediction for a Debye length of 3 nm, while the experimental data obtained from the sample without additional salt follow the prediction only qualitatively.
To investigate the influence of the evanescent wave penetration depth on the measured values of G we performed experiments at different angles of incidence.Fig. 11 shows data obtained from Q J -scans at two different penetration depths together with the corresponding theoretical prediction.Although the agreement is very good, the scatter of the experimental data is larger than the difference between the predictions for different penetration depths.Thus, the influence of the penetration depth is too small to be reliably detected experimentally.
So far, we have shown that in dilute suspensions, PEGylated FBMA particles behave like charged spheres, and their nearwall dynamics can be suitably described by a combination of hydrodynamic wall effect and the particles' static interaction with the wall.However, a more interesting and challenging question concerns the near-wall dynamics at increased volume fractions, at which the system changes from an isotropic solution to an ordered structure, stabilized by long-ranged particle  interactions.As a first step to tackle this question we show data obtained from a suspension with a particle volume fraction of 7% in Fig. 12.The graph displays the normalized G versus Q J 2 data, obtained at kR = 1.81 from a suspension without added salt.As a general trend we observe that the relaxation rates are lager at f = 0.07 than at infinite dilution, which is represented by the dashed dotted line.Further, there is a significant deviation from the linear Q J 2 behaviour, which was observed at low volume fraction.The relaxation rates rather undulate with Q J 2 which is very reminiscent of the behaviour we reported earlier for hard spheres at high volume fractions.For comparison we also show experimental data obtained from a suspension of hard spheres at a volume fraction of 30% together with a set of simulated data for the same system.Note that these data were obtained at a r = 27.81 and n S E 1.49, thus the effect of back reflections will only become visible at (Q J d) 2 \ 50.Based on the observed qualitative similarity we conjecture that also the near-wall dynamics of the PEGylated FBMA particles at f = 0.07 are significantly affected by the suspension's structure factor, which is in line with our observation that suspensions with a volume fraction of 7% are close to the disorder/order coexistence of the system.A systematic experimental investigation and a quantitative theoretical analysis of this effect will be the subject of a future contribution.

Conclusion
We have synthesized PEGylated fluorinated particles by an emulsion polymerization, for the study of near-wall dynamics of colloidal particles in aqueous solution.Bulk characterizations show that the particles have a core shell structure, a negligible polydispersity, are stable against compression and addition of salt, and can be easily index-matched to aqueous solutions, which makes them a good model system for scattering experiments in aqueous suspensions at increased concentrations up to the transition from an isotropic solution to an ordered phase and beyond.By static light scattering we have shown that a coexistence of isotropic and ordered suspensions already occurs at a volume fraction of 8% in DMSO/water, if no additional electrolyte is present.In future we will use this system for the systematic investigation of near-wall dynamics of particles with long-ranging interactions.In this contribution we showed that EWLDS measurements on this system produce reliable and reproducible data.At low volume fractions the near-wall dynamics can be well predicted taking into account the particles hydrodynamic and electrostatic interaction with the wall.However, even at moderate volume fractions the nearwall dynamics are reminiscent of the behaviour which was observed for hard spheres at much higher volume fractions, which is attributed to the long range particle-particle repulsion.A thorough experimental investigation of the interaction effects on near-wall dynamics and the design of a theoretical framework for their prediction is the objective of ongoing work.The dashed line represents simulation data for the hard sphere system, the dashed dotted line is the prediction for charged spheres at low f according to eqn ( 5) and ( 9) and the dotted line is a guide for the eyes.
After 1 hour of homogenization the stirring speed is reduced to nominally 150 rpm and an adapter for a PTFE tube is added to the flask.A solution of This journal is © The Royal Society of Chemistry 2016 Soft Matter, 2016, 12, 8485--8494 | 8487

Fig. 1
Fig. 1 Scattering geometry of an EWDLS experiment with spherical symmetry.

Fig. 2
Fig. 2 Time auto correlation functions recorded during a Q J scan in EWDLS at a r = 601 and values of y as indicated in the legend.Top: Intensity correlation functions g ˆ2(t) À 1 as measured.Bottom: Field correlation functions calculated by the generalized Siegert relation, eqn (1).

Fig. 4
Fig. 4 2D projections of PEGylated FBMA particles as observed during a tomographic tilt series.

Fig. 5
Fig.5Visualization of the 3D reconstruction.In the yz plane a spherical shape is observed, while in the zx and xy planes the reconstruction is elongated along the direction of the missing wedge.

Fig. 7
Fig. 7 Structure factors obtained from PEGylated FBMA particle suspensions in index-matching DMSO/water mixtures with different volume fractions as indicated in the legend and a screening length of l D E 95 nm.

Fig. 8
Fig. 8 Bulk DLS result of PEGylated FBMA particles dispersed in water (open red squares) and in DMSO/water mixture (open blue circles).Symbols represent experimental data and lines are linear least squares fits with slopes equal to the diffusion coefficients, D 0 , as indicated.

Fig. 9
Fig. 9 Normalized initial relaxation rates G as a function of the parallel scattering vector component squared.Symbols represent experimental data obtained from a f = 0.01 sample at different Debye screening lengths, as indicated in the legend.Lines represent theoretical predictions for infinite dilution according to eqn (5) and (9).The error bars assigned to the data obtained from the sample with added salt, represent standard deviations from three repeated measurements.

Fig. 10
Fig. 10 Normalized initial relaxation rates G as a function of the normal scattering vector component squared.Symbols represent experimental data obtained from a f = 0.01 sample at different Debye screening lengths, as indicated in the legend.Lines represent theoretical predictions for infinite dilution according to eqn (5) and (9).

Fig. 11
Fig. 11 Effect of evanescent wave penetration depth on measured G. Symbols represent experimental data obtained from a f = 0.01 sample at different penetration depths.Lines of equal color represent theoretical predictions for infinite dilution according to eqn (5) and (9).Penetration depths are indicated in the legend.

Fig. 12
Fig. 12 Comparison of Q J -scans on a PEGylated FBMA suspension with f = 0.07 (open circles) and a hard sphere system at f = 0.3 (open squares).The dashed line represents simulation data for the hard sphere system, the dashed dotted line is the prediction for charged spheres at low f according to eqn (5) and (9) and the dotted line is a guide for the eyes.

Table 1
1.Best fitting parameter values obtained from fitting the static light scattering data to eqn(7)