Aggregation and sedimentation of shattered graphene oxide nanoparticles in dynamic environments: a solid-body rotational approach

Nanoparticle (NP) aggregation is typically investigated in either quiescent or turbulent mixing conditions; neither is fully representative of dynamic natural environments. In groundwater, complex interacting influences of advective – diffusive transport, pore tortuosity, and the arrival of aggregates from up-gradient pores impacts the aggregation behaviour of NPs, whereas in surface waters, continuous mixing of fresh particle and aged aggregate populations amends aggregation rates. To mimic such conditions, a cylinder reactor containing shattered graphene oxide NP ( < 100 nm) suspension was set to rotate with a Reynolds number (Re) close to one and with zero shear. Two main aggregation phases were then observed. Up to 250 – 350 min, NP remained near the rotational axis longer than in static conditions, giving rise to higher aggregation rates interpreted as an enhanced perikinetic aggregation and differential sedimentation due to mixing with resuspending aggregates. In this phase, a population-balance model estimated an attachment efficiency > 5 times in the rotating system than in the static system. Later (5 – 13 h) aggregates collided with extensively each other, broke, and reformed on the rotating cylinder wall giving rise to larger, denser aggregates ( > 1 cm). These results thus shed new light on the differences in aggregation behaviour between porous media and other natural environmental systems compared to quiescent batch experiments. fate of nanomaterials, to help developing predictive models for managing the release of these nanomaterials. To investigate the impact of environmental dynamics on the aggregation processes of shattered graphene oxide nanoparticles, this study uses a solid-body rotational approach to mimic dynamics of interacting populations of nanoaggregate with different sizes and structures. This sheds light on the greater aggregation rates observed in the aquatic environments such as groundwater and surface water systems compared to those observed in simplifying laboratory batch experiments.


Introduction
With the increasing global proliferation of commercial applications for nanomaterials, they are becoming increasingly released into the natural environment either accidentally in waste streams or deliberately within environmental clean-up, agronomy, and petroleum reservoir recovery applications. 1 Among various nanomaterials, the use of graphene nanomaterials is rapidly growing due to their broad constructional, industrial, environmental, and medical applications and arising from their use with other nanomaterials as nanohybrids. [2][3][4][5][6][7] Graphene nanosheets may end up in the environment in various forms of particulate matter such as crumpled graphene, [8][9][10] multilayer graphene, 4,5,11 shattered nanosheets, 12 and fractal aggregates/hetero-aggregates. 13,14 Once they have entered aquatic environments, aggregation of such particles can significantly affect their functionality and transport behaviour, particularly in aqueous and porous media. 4,5,15 Despite the existence of many studies on homo-and hetero-aggregation of various nanoparticles (NP), 14,[16][17][18][19][20][21] and abundant reports on the impacts of various factors on the aggregation behaviour of these NP, [22][23][24][25][26] it still remains a problem of how system dynamics modify aggregation. It is particularly of paramount importance in groundwater (GW) transport to understand how complex multi-cascade processes of advective and diffusive transport, [27][28][29][30] tortuosity in porous media, 31 and the arrival of aggregates from upgradient pores 32,33 impact the aggregation behaviour of NP. Additionally, it is critical to investigate how resuspension phenomenon occurring within surface water (SW) bodies, such as lakes, river, and sea at various scales, influences the aggregation of NP in these environments. [34][35][36][37] Typically, current experimental techniques investigating NP aggregation dynamics are limited to simple quiescent or turbulent mixing batch experiments which may not adequately capture aggregation conditions in real systems. 38,39 Whilst research into incorporating aggregation models within porous media transport models continues, 38,40,41 steps must be taken to validate the individual aspects of aggregation models in dynamic media. One of these impacts can be the heterogeneous or anisotropic mixing in the population of nanoaggregates induced by the environmental dynamics such as cascading/cycling conditions within GW and SW systems [27][28][29][30]32,33,36,37,[42][43][44] that can modify the particle size distribution (PSD) evolution trends under aggregation/sedimentation conditions. To investigate such systems, we used a slow rotating cylinder in which these impacts are simulated by slow continuous revolutions of the cylinder. This simple apparatus allows monitoring the particle size dynamics during aggregation which is cumbersome or impossible in other experimental approaches used for studying the NP fate, such as packed column experiments 45 or mesocosm tests. 46 Rolling cylinders have been used extensively in marine science to mimic the natural condition of marine snow aggregate formation. 47 This system consists of a cylinder which is saturated with water and is placed horizontally to rotate around its axis slowly so that no significant shear is induced in the fluid inside. Since the initial condition is static, i.e., the cylinder is initially at rest at zero time, the initial rotations can cause shear in the fluid. Once this has passed, the system rotates as a solid body. 48,49 This technique, 50,51 has been also utilized to study the interactions of bacteria, 52 organic carbon accumulation in marine sediments, 53 aggregation of diatom, 54 and their interactions with resuspended sediments, 55 larvae interactions with anisotropic fluid motions, 56 and interactions of microplastic with phytoplankton aggregates. 57 However, to the best of our knowledge there has not been any systematic application of rolling cylinders to determine the aggregation of engineered NP in aqueous environmental systems to date. To mimic the condition of porous media, the system was designed with a Reynolds number (Re) below one. The main objective of using such a system was to investigate how NP aggregation behaviour changes over time when the local mass concentration is kept constant, a phenomenon that can occur in pores of subsurface porous media and in resuspending surface waters. Using such a rotating cylinder, the angle between NP displacement directions and the gravity or removal direction is continuously changing without significant mechanical shear, thus preventing particle removal via settling within the axial area of the cylinder. This allows aggregates to develop to very large dimensions until gravity forces impacting on large grown aggregates which are probably more compact than initial aggregates overcome rotational displacements, and the system enters an unstable situation where aggregates tend to settle. The primary aim therefore is to shed light on unexpected higher aggregation rates of NP in GW systems when compared with many standard experimental and modelling results, and on how NP aggregation can be influenced by resuspension and population mixing in SW systems. 34,38,58,59 This may underpin development of new mathematical descriptions of aggregation mechanisms for a more accurate consideration of these mechanisms in models of NP fate and transport in the environment which take other complex concurrent transport/transformation phenomena into account. 38,41,60 Since the largest use of graphene nanomaterials may result from their high strength and flexibility (∼100 stronger than steel), 4 it is possible that they release into the environment as shattered graphene because they may undergo extreme physical or chemical stresses during the usage or upon entering waste streams. Therefore, in this study we use shattered graphene oxide (SGO) particles with a primary hydrodynamic size (D H ) of 90 nm. In the scope of this study, which lumps the impact of particle surface characteristics and interactions in one free parameter, i.e., attachment efficiency, within the fitting process, SGO may also act as a representative for other NP types due to their small hydrodynamic size (D H < 100 nm) and since in terms of general aggregation behaviour, e.g., hydrodynamic size and mass concentration trends over time, they exhibit similar patterns to other NP, such as hydroxyapatite (HAp) NP. 40 It should be mentioned that the pristine GO sheets may not represent NP due to their large aspect ratio and being considered as a 2-D material.

Rotating cylinder setup
Acrylic cylinders were manufactured with inner diameters of 5 cm and length of 19.1 cm to rotate at 4.71 rotations per hour (rph) corresponding to a Reynolds number of 0.91. Methylene blue was used to investigate both the dispersion of matter inside the cylinder and the flow regime as a solidbody rotation. In this case the dye was injected through a needle as an instant point source into the cylinder to see how its diffusion/dispersion is affected by the rotations. A beltpulley was set to reduce, by >200 times, the rotational velocity of a laboratory stirrer (Eyela NZ, Japan) with an adjustable revolution rate as illustrated in Fig. 1 and S1. † The connection of cylinder to the rotational apparatus was filled with soft plastic, preventing any transfer of shaking from the apparatus to the cylinder. The prepared NP dispersions were poured into the cylinder in the beginning of each experiment to set up each case with an initially homogeneous dispersion. Sampling for particle size measurement was challenging due to the possibility of aggregate breakage resulting from any shear force. To alleviate this, we selected all the channels through which the samples had to pass including syringe tips and a glass pipe inside the cylinder, with an inner diameter larger than ∼2 mm. Sampling was performed by syringe removal (without needle) with a minor flow rate to minimize the impact of shear force on aggregates. With great care they were then directly transferred to the dynamic light scattering (DLS) cuvette for size measurement. Samples were taken using a straight glass pipe centred on one side of the cylinder up to approximately one-third of the length (Fig. 1). In order to keep the reactor saturated at all times during the experiment and to prevent bubble formation, a syringe with a long needle was used to inject the stock dispersion into the cylinder simultaneous with sampling, with the needle leaned towards the wall of the cylinder and entered up to approximately two-third of the length of the cylinder, to be far enough from the sampling point in the middle (see Fig. 1). The volume of each sample was <1.5 mL and an average of 16 samples were taken in each experimental setup. Therefore, the total volume of reactor exchanged with the initially provided stock dispersion was around 24 mL which is less than 10% when compared with the total volume of the cylinder (∼375 mL). After obtaining, the samples were immediately analysed for size and PSD using DLS (Malvern Zetasizer, Nano ZS model, UK) selecting the number of runs to 5 each with a duration of 10 sec and adjusting the beam attenuator to 11 and the position of measurement at 6.5 mm. The refractive index of graphene oxide (GO) used in the calculation of volume-based PSD was assumed as 1.333 and the absorption as 0.01. It should be noted that the sensitivity of volumebased PSD to these factors was minor despite their being expected to be influential. Based on preliminary batch experiment, ionic strength (IS) was selected at 0.5 and 0.75 mM CaCl 2 with a pH of 6. These IS values are less than or close to the critical coagulation concentration (CCC) (results not shown).
The derived count rate (DCR) of DLS output has been used as a measure of mass concentration. [61][62][63] However, since DCR significance under the reaction limited aggregation (RLA) remains unvalidated, we also determined the mass concentration in each sample using UV-vis spectroscopy (Hitachi, U-4100 model, Japan) at a wavelength of 230 nm, with a calibrated absorption curve. With such a curve (intercept = 0) we used measured absorbance normalized by the initial absorbance in each experiment as representative of normalized mass concentrations.

SGO synthesis
Graphene oxide was synthesized by the thermal exfoliation of natural graphite flakes using a modified Hummers method. 64 This was conducted through constant stirring of graphite (2.0 g) in a solution of 225 mL sulfuric acid and 50 mL phosphoric acid, followed by gradual addition of 5 g potassium permanganate at a constant temperature of 35°C and continuous stirring, which was maintained for 10 h. The mixture was then diluted with 225 mL of deionized (DI) water and kept in an ice bath, subsequently, 3 mL of H 2 O 2 was added for permanganate residual removal, with the mixture centrifuged and washed with 1 M hydrochloric acid. This was followed by washing with phosphate buffer (pH 11). Finally, the synthesized GO was washed with DI water repeatedly to reach the neutral pH. Shattered GO was prepared via intensive ultra-sonication of a GO dispersion (2 g L −1 ) to achieve initial uniform hydrodynamic size below 100 nm. This was conducted using a probe sonication at a power of 40 W for 2 h with 30 sec stops following 30 sec sonication intervals. This dispersion was then centrifuged for 30 min at 19 500 ± 500 rpm to remove the larger fraction of particles. Finally, the dispersion was passed through a 0.45 μm syringe filter, and the filtrate kept in the dark at 4°C as the stock dispersion. Sodium azide was added (10 mM) to each experimental dispersion to prevent the growth of bacteria during analysis and due to its buffering ability in maintaining a constant pH during the course of experiments. 40

Experimental procedure
The cylinder tests were performed as follows: (1) prepare particle dispersions in DI water with a final SGO concentration 50 mg L −1 , and sodium azide concentration 10 mM; (2) adjust pH at 6 ± 0.05 with NaOH/HCl (100 mM) while stirring with a magnet stirrer; (3) ultrasonicate for 5 min in sonication bath, add electrolyte, (CaCl 2 ) to reach final concentrations 0.5 and 0.75 mM, stop stirring 5 sec after adding the electrolyte, and take the first sample; (4) fill the cylinder with the prepared dispersion, emplace the syringe and the sampling pipe, seal, and leave at rest for static cases or on the rotating mode for dynamic cases (in total taking ∼10-15 min from the time of stopping the stirrer until the complete setup of the cylinder); and (5) take samples (<1.5 mL) at certain times to be analysed immediately by DLS. After obtaining each sample the glass pipe was evacuated using a syringe to prevent interference of the solution inside the pipe with the next sampling. After analysis by DLS, the same samples were carefully preserved at 4°C for analysis using UV-vis, which were carried out following dilution by 7.5 times and sonication for 5 min to obtain a uniform dispersion. Selected experiments which could have been prone to uncertainties (e.g., due to higher IS) were carried out in duplicate. Control experiments were undertaken in a similar condition to cylinder experiment but inside 3 mL DLS cuvettes with time resolved online measurement.
In these experiments after ultrasonication of the suspension for 5 min, the electrolyte was added, and the sample was immediately vortexed for 5 sec. Then it was transferred to a disposable cuvette, after which time-resolved DLS measurement was started immediately. The whole process duration, from ultrasonication until the start of the first measurement, was 70 ± 20 sec. The sample volume of 3 mL corresponds to a measurement depth of ∼23.3 mm which is close to the radius of the cylinder. In control experiments, time-resolved DCR data normalized to the initially observed DCR was used as an indicator of mass concentration. 40,[61][62][63] Theory and modelling Flow regime. A simplified form of Navier-Stokes equation in cylindrical coordinate can be expressed as: 48 where u θ is the angular component of the linear velocity, r is the radial distance from the centre of the cylinder, and ν is the kinematic viscosity [L 2 T −1 ]. For the initial condition, u θ = 0, and boundary condition, u θ = ωa cyl at r = a cyl , where a cyl is the radius of the cylinder and ω is the angular velocity of the cylinder. The analytical solution of eqn (1) is: 48,65 where A i = 2ω(λ i J 2 Ĳλ i a cyl )) −1 , J 1 and J 2 are Bessel functions of the first kind and orders 1 and 2, respectively, and λ i ≈ iπa cyl −1 . The horizontal and vertical rectangular components of velocity, u x and u y are: where U k is settling velocity of aggregates in size class k [L T −1 ]. Among various approaches to settling velocity, 66-69 a power-law model can best describe the settling velocity of NP aggregates as derived in the ESI: where a 0 is the radius of the primary particles, a k is the radius of aggregates, D f is the fractal dimension of aggregates, g is the gravitational acceleration, μ is the viscosity of fluid, ρ 0 is the density of primary particles [M L −3 ], and ρ w is the density of water [M L −3 ].
Once the system reaches the equilibrium conditions, i.e. solid-body rotation, particle trajectories can be calculated as: 48 where r 0 and θ 0 represent the initial position of a particle in cylindrical coordinates and x and y represent the position of a particle in rectangular coordinates at time t. The shear rate, γ, resulted from the gradient of velocity across the cylinder radius is given as: 48

Particle aggregation
A flexible population-balance model known as fixed pivot (FP), 70 capable of considering initial PSD and consequent evolution of PSD in early and late stages of aggregation is used to model the dynamic behaviour of hydrodynamic diameter, overall mass concentration, and PSD for SGO aggregation experiments. Conserving two properties of mass and number, the model equation for aggregation combined with sedimentation term is as follows: 70 where n k is the particle number concentration of aggregates in size class k [L −3 ], β is the collision frequency, α is the attachment efficiency factor, v is the representative volume of each size class in the grid, δ is Kronecker's delta, Z s is sedimentation depth or measurement depth [L], and η k is: An explicit forward Euler scheme was used for the time discretization of eqn (9) with an adjustable time-step, as this can be more efficient than higher-order schemes for solving "stiff" problems such as population balance models. 71,72 Since the Brinkman permeability-based model [73][74][75] performs best for modelling collision frequencies, this approach was used in the present study. All collision kernels were expressed based on aggregate volume so that the impact of aggregate shapes which are unknown will interfere less in the model results. The measured initial PSD was directly used as the initial condition. A new optimization algorithm code 40 was also used here to estimate the parameter values (α and D f ) by matching the modelled with experimental data for mean hydrodynamic diameter (D H ) and normalized concentration data (C/C 0 ). In these fittings, the contribution weight of the D H was assumed to be twice C/C 0 . Aggregation modelling was performed only at a single spatial point, i.e., the centre of the cylinder cross-section, corresponding to the point of measurement and considering the sedimentation depth in eqn (9), Z s , as the radius of the cylinder. Full model equations together with details of extensive model testing and the MATLAB codes and optimization algorithm are available in the ESI † and in Babakhani et al. 40 A summary of the model parameters and experimental characteristics are provided in Table 1.

Flow regime and particle trajectory analysis
The experimental results of the cylinder with a drop of dye injected in the beginning of the experiment (Fig. S2 †) shows that generally within the first ∼10-20 min a thin layer of dye forms at the rotating wall. After ∼60 min, however, this thin layer is still noticeable, but begins to slowly diffuse into the bulk solution and after ∼90-100 min it is barely visible. These results suggest that the rotational regime can be considered as solid-body and the diffusion of dye is the only mechanism of mixing in the cylinder after ∼10-20 min. Consistent with these outcomes, the calculations of rotation velocity and shear rate for the cylinder, shown in Fig. S3, † demonstrate that within 10 min rotational velocity reaches a steady-state (solid-body rotation) and the shear stress produced in cylinder fluid totally disappears (blue lines). This duration is negligible compared to the total course of experi-ment >5 h and is not captured in the sampling intervals when NP dispersions are used in the cylinder. Both experimental and modelling results confirm the validity of assuming steady-state and zero-shear conditions in the system. Fig. 2 shows results of particle trajectory analysis modelled using eqn (1)-(8) with varying fractal dimension and particle radius for two scenarios: (1) particle is initially placed in the middle of the cylinder radius and (2) particle is initially placed near the wall of the cylinder. Based on this analysis, the movement of NP and their aggregates ranging in radius from 20 nm to 100 μm depends significantly on the fractal dimension. For aggregates up to a radius of ∼1 μm and for fractal dimension up to 2.5, their trajectories follow the "fixed" solid-body rotation regime as illustrated by Tagawa et al. 76 However, for aggregate radii in range of 10-100 μm with D f ≥ 2.5, trajectories start to deviate from this regime and show a tendency towards the so-called "cascading" or "suspending" regimes. 76 Sedimentation velocity is so rapid for aggregates of radius >100 μm and D f = 2.9 that they interact with the cylinder wall instead of rotating with the solid body. In agreement with these results, Fig. S4 † shows that velocity vectors exhibit a nonuniform distribution for aggregates with 10 μm radius and D f = 2.9, highlighting the importance of aggregate fractal dimensionality and size in controlling particle trajectories.

Aggregation experiment results
Experiments were undertaken in the cylinder under both rotating and static (non-rotating) conditions, and as a control within a cuvette under static conditions with continuous monitoring using dynamic light scattering (DLS).  shows that the averaged hydrodynamic diameter (D H ) follows a log-linear trend with time at 0.5 mM CaCl 2 in all experiments. However, at 0.75 mM CaCl 2 (Fig. 3c) the gradient de-creases over time after a rapid initial rise, which we interpret as a transition from early (fresh, monodisperse particle population) to late (aged, polydisperse particle population) stage with virtually varying D f and particle radii considering two scenarios: (1) for particle initially placed in the middle of the cylinder radius (left-hand side panels) and (2) particle initially placed near the wall of the cylinder (right-hand side panels). The rotational direction is counter clockwise. conditions within a coupled aggregation/sedimentation system. Under the solutions present and within the selected periods for size measurement (up to 300-350 min), the sedimentation is not generally significant ( Fig. 3b and d). The trend of D H measurement at 0.5 mM CaCl 2 within the static cylinder closely matches that in the control, demonstrating the validity of the design of sampling in the cylinder experiments. This was also the case at 0.75 mM CaCl 2 but only in the early stages <120 min, as shown in Fig. 3c after this time D H in the control continues to rise, whereas D H in the static cylinder tends towards an asymptote over time. The reason for this behaviour is not clear, it could be due to minor differences in the conditions between the two experiments, such as variation in temperature or dimensionality between the cylinder and the cuvette. One might attribute this behaviour to the breakage of aggregates larger than ∼500 nm during sampling from the cylinder experiments. However, this is less likely to be a reason for the observed lower growth of aggregates in the static cylinder compared to the control system, because the control experiment shows a level of sedimentation which is not observed for the static case (Fig. 3d) whereas sedimentation measurement, is not expected to be affected by the breakage. This suggests that in this case there might be a level of uncertainty for the control test rather than the occurrence of breakage in the cylinder samples as also indicated by the larger standard deviations in the control after ∼230 min.
Generally, the experimental observations using the rotating cylinder initially show higher aggregation rates than that of the static cylinder (Fig. 3a and c). In the case of 0.5 mM CaCl 2, normalized mass concentration curves show increases in concentration to above the initial uniform concentration at the point of measurement (Fig. 3b). Interestingly, these results are consistent for both measurement techniques, i.e., UV-vis used for the cylinder samples and DLS derived count rate used for the control data, and are also in agreement with previous measurements of DCR for HAp NP aggregation at IS below CCC. 40 Babakhani et al. 40 used the FP model already to validate the use of DCR data as representative of mass concentration during NP aggregation at IS above CCC leading to a decline in mass over time. The present study's results tend to agree with the use of DCR data as representative of mass concentration at IS below CCC as well.
The experimental photos of the rotating cylinder at 0.75 mM CaCl 2 at different times are shown in Fig. S5 † and the magnified images of aggregates formed at each time are shown in Fig. 4. It appears that visible aggregates form on the wall after ∼290 min, corresponding to the start of removal in the sedimentation curve after 250-300 min (Fig. 3d). Large networks of gel-like aggregates form at the bottom after ∼315 min which as a result of extensive sliding on the wall and consequent collisions, breakage, and reforming, become denser over time until finally they form larger (>1 cm) granular dense aggregates as indicated by their darker colour compared to those formed at earlier times ( Fig. 4 and S5 †). 77,78 When finally, the sizes of aggregates grow larger and their structure become more compact to induce their higher settling velocity due to their higher bulk density, the settling velocity overcomes the resuspension due to the cylinder rotation. This results in millimetre-sized aggregates appear immediately on the cylinder wall sliding and colliding consequently producing compact centimetre-sized SGO granules appear on the wall of the rotating cylinder in longer time (5-13 h), and the suspension becomes clear of the SGO NP.

Aggregation modelling results
Data from the first 250-350 min of experiments (before the evolution of larger, settling macro-aggregates) were used to fit models of the aggregation in the suspended phase. The estimated parameter values are shown in Fig. 5 and are also available in Table S1 † along with Nash-Sutcliffe 15 goodness-of-fit criterion, R 2 NS , values. The model set used in this study (FP population-balance technique with the power-law model for settling velocity and the Brinkman permeability model for collision frequencies) can fit both D H and normalized concentration (C/C 0 ) well with an average R 2 NS 0.758 ± 0.18 for D H data. Parameter values estimated for the model matched to the control data agree well with those of static cylinder measurements at electrolyte concentrations of 0.5 mM CaCl 2 and 0.75 mM CaCl 2 (<29% difference). The rolling cylinder experiment shows 5 and 7 times greater attachment efficiency (Fig. 5a) than that of the static system at 0.5 and 0.75 mM CaCl 2 , respectively. Likewise, D f estimated in the rolling cylinder experiment is on average 10% and 31% larger than that of the static cylinder at 0.5 and 0.75 mM CaCl 2 , respectively (Fig. 5b).
The results clearly indicate that conditions within the rotating cylinder enhance aggregation rate and leads to aggregates with more compact structures (via increasing D f ). For higher IS than those of present study's experimental conditions and longer times than those of aggregation simulations, the impact of system dynamics may well be higher (Fig. 4 and S5 †). Although the process of the aggregation at longer times (sliding on the wall) is different from that at early stages (aggregation of suspended particles/aggregates), occurrence of both phases is relevant to the natural environment. In GW system, large aggregates attached to porous medium solid surfaces can be translated, rolled over and collide with each other, thereby aggregating in the retained phase. 32,38 They can also be detached and transported to a new pore where they can interact with new NP or aggregates, 32,38 whether in the suspended phase or in the retained phase. Likewise, in a SW system the aggregates which have already settled onto deeper depths, can be transported with sediments or be resuspended where they could encounter new particles. It should be noted that in both modes, aggregates View Article Online become more compact, shown either by enhanced D f resulted from the model or by change in colour of aggregates to a darker colour observed in the cylinder. 77,78 The current cylinder condition ignores many aspects of convoluted groundwater pores, such the constricted geometry of pores, shear effect, nonuniform flow regime, and dynamic interactions between suspended particles and solid surfaces. [79][80][81][82] Nevertheless, ignoring such aspects may allow a better understanding of the dynamic impacts on the individual aggregation/ sedimentation mechanisms. For instance, including shear stress in the system 39,82,83 can overshadow the impact of system dynamics due to addition of other complexities such as break-up of aggregates. 84 This cylinder can be further modi-fied in future studies, e.g., by modifying the wall surface roughness or adding short edges/partitions on the inner wall in order to consider more complexities of natural conditions. The system developed different stages in respect to the interplay between aggregation and sedimentation mechanisms within the rotating cylinder. Up to 250-350 min, where the sliding of aggregates on the wall is still not significant, the maximum D f obtained in the present study is 2.34 (at 0.75 mM CaCl 2 within the rotational cylinder) and the maximum aggregate D H is below 1 μm (Fig. 3a and c). According to the particle tracking analysis results ( Fig. 2 and S4 †), for these values of D f and D H , the effect of gravity on aggregates may not be strong enough to deviate aggregate trajectories from a uniform solid-body regime. In agreement with this, the results from sedimentation over time (Fig. 3b) show minor or no removal of NP over the periods up to 250 or 350 min. Hence, within these periods the only difference between the static and rotating systems is seemingly that particle orientations with respect to the gravitational force vector are continuously changing for aggregates in the rotational system while they remain constant within the static system. This might maintain the local particle concentration within the axial area high at initial stages and therefore enhance the number of collisions resulting from Brownian motion, thereby increasing the overall rate of aggregation. At 0.75 mM CaCl 2 , at times <150 min normalized mass concentration data shows a slight decrease in the suspended mass in both the control and static cylinder experiments which is not noticeable in the case of rotational cylinder (Fig. 3d). However, after 250-300 min mass removal in the rotating cylinder becomes increasingly significant as the aggregates become large and dense enough for the gravity force on them to overcome the rotational force leading to their settlement on the cylinder wall. Enhanced aggregation and removal within the rotating cylinder compared to the static cylinder emanates also from differential sedimentation mechanism 85 of aggregation. In fact, attachment efficiency is generally expected to be a function of NP surface characteristics. 86 However, the present study indicates that system dynamics increase the effective attachment efficiency, even though chemical conditions controlling surface properties were consistent between rotating, static and control experiments. To investigate this behaviour, we calculated the collision frequencies between particles with radii of 50, 500, and 5000 nm and particles with all other size classes up to a radius of 15 μm and for various D f (1.5-2.9) as illustrated in Fig. S6. † Based on collision frequency formulations used in the present study, representing the state-of-art equations for environmentally-relevant NP, 40,72,87 the impact of differential sedimentation mechanism of aggregation, which is notably affected by the structure of aggregates (D f ), is around 3-5 orders-of-magnitude lower than that of perikinetic aggregation, which is impacted less by the structure of aggregates (D f ), for the range of aggregate sizes observed in this study in early stages of aggregation. Experimental photos (at 0.75 mM CaCl 2 , Fig. 4 and S5 †) suggest that at 290 min aggregates are sufficiently large (>10 μm, based on Fig. 2) and/or compact enough (D f > 2.5, based on Fig. 2) for the trajectories to deviate from solid body rotation and to settle on the wall, suggesting that in this condition sedimentation and therefore differential settling mechanism of aggregation is significant. However, the current collision frequency formulation may not be able to take this effect into account because for a similar range of aggregate size/fractal dimension they predict negligible contribution for differential sedimentation. This deviation may also arise from restructuring or change in the fractal dimension over the course of the experiment increases its significance in the later stages of aggregation. 40,68,88 Taking such a transient change in the fractal dimension, which causes a variation in the particle volume discretization, into account needs special computational consideration. Modelling the whole rotating cylinder domain for aggregation may shed light on the complex interactions of resuspending aggregates and the reasons for the observed differences in the attachment efficiency between the two systems of dynamic and static. Furthermore, additional experimental analysis on the structure of aggregates is needed to verify the estimated fractal dimension values, e.g., using static light scattering (SLS) 88 to monitor the dynamics of the aggregate structure online during the test. The impact of system dynamics can be even more crucial where the existence of natural colloids induces hetero-aggregation with NP. 16 In such systems, the initial diversities in density, size, and interactions of particles may affect both the initial availability of surfaces for Brownian-driven aggregation and differential settling aggregation. These are beyond the scope of current study and may be subject of future studies.

Conclusions
In the early stage of experiments reported here, whilst the majority of particles remain suspended, maintenance of mass concentration in the axial area of the rotating cylinder increased the number of collisions induced by Brownian motion and arising from differential sedimentation collisions with resuspending larger aggregates. This enhanced the aggregation rate which emerged relating to the attachment efficiency estimated based on model fit to experimental data. As the system evolves and aggregates grow sufficiently to show significant sedimentation, resuspension of settling aggregates in the rotating cylinder increased differential sedimentation, further enhancing the aggregation rate. This consequently leads to an immediate formation of aggregates on the cylinder wall that afterwards start rolling and sliding on the wall. These processes give rise to the formation of more compact aggregates, leading to an increase in fractal dimension in all stages as revealed by greater fractal dimension in the rotating system compared to that of static system and control and indicated by dark colour of aggregate formed on the cylinder wall at the end of experiments. Numerous studies have reported higher aggregation rates in groundwater, surface-water and other natural environmental systems when compared to quiescent batch experiments. 34,38,58,59 The results from this study demonstrate that this effect can be understood in terms of systems where, even in the absence of hydrodynamic shear, particle suspensions are to some extent prevented from sedimentation. Such effects are simulated here by continuous rotation in a slowmoving cylinder, which models slow flow through tortuous pore networks in natural systems which are characteristic of changing flow directions and mixing of old and new particle populations. It should be noted that considering the importance of hetero-aggregation (aggregation of NP with natural colloids) in the fate of NP in natural waters, 16,89-91 the present study serves as a first stage to investigate heteroaggregation in more realistic conditions. The system used in this study has also important applications for water treatment, e.g., as a settling unit that can be optimized through adjustment of the fractal dimension and aggregation rate where the adsorption of contaminants to NP is maximised whilst the NP removal due to settling is minimised. Further studies are required to consider a variable fractal dimension in the model for taking reorganization of aggregates into account and to investigate the spatial variations of aggregate sizes in both phases of suspended and settled particles. It also remains for future studies to experimentally verify the model estimated fractal dimension values through direct, online measurements.

Conflicts of interest
There are no conflicts to declare.