Juan D.
Olarte-Plata
,
Gøran
Brekke-Svaland
and
Fernando
Bresme
Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, White City Campus, 80 Wood Lane, London W12 0BZ, UK. E-mail: j.olarte@imperial.ac.uk; g.svaland15@imperial.ac.uk; f.bresme@imperial.ac.uk
First published on 8th May 2020
We investigate the impact of nanoparticle roughness on the phase behaviour of suspensions in models of calcium carbonate nanoparticles. We use a Derjaguin approach that incorporates roughness effects and interactions between the nanoparticles modelled with a combination of DLVO forces and hydration forces, derived using experimental data and atomistic molecular dynamics simulations, respectively. Roughness effects, such as atomic steps or terraces appearing in mineral surfaces result in very different effective inter-nanoparticle potentials. Using stochastic Langevin Dynamics computer simulations and the effective interparticle interactions we demonstrate that relatively small changes in the roughness of the particles modify significantly the stability of the suspensions. We propose that the sensitivity of the phase behavior to the roughness is connected to the short length scale of the adhesive attraction arising from the ordering of water layers confined between calcite surfaces. Particles with smooth surfaces feature strong adhesive forces, and form gel fractal structures, while small surface roughness, of the order of atomic steps in mineral faces, stabilize the suspension. We believe that our work helps to rationalize the contrasting experimental results that have been obtained recently using nanoparticles or extended surfaces, which provide support for the existence of adhesive or repulsive interactions, respectively. We further use our model to analyze the synergistic effects of roughness, pH and ion concentration on the phase behavior of suspensions, connecting with recent experiments using calcium carbonate nanoparticles.
Calcite surfaces are prone to dissolution and crystallization, as a consequence the surfaces can develop different levels of roughness. Experimental studies of extended calcium carbonate surfaces using the surface forces apparatus,7 indicate that roughness effects can play a significant role in determining the interactions between calcite surfaces. Roughness effects are manifested in the measurement of repulsive interactions when the surfaces are immersed in water. Experimental studies of CN aqueous suspensions, with nanoparticle diameters of 60–70 nm, demonstrated the formation of gel phases,4 and therefore the existence of strong adhesive interparticle interactions.8 The mechanical strengthening of calcium carbonate pastes obtained from recrystallization of amorphous calcium carbonate and vaterite mixtures, has been rationalized considering the smoothing of the nanoparticle surfaces at grain contacts.9 Recent molecular dynamics simulations of flat calcium carbonate surfaces immersed in water do also provide support for attractive interactions. Strong adhesive minima were observed at inter-surface separations of about 1 nm. This strong adhesion is mostly driven by the layers of water molecules adsorbed at the calcite surfaces. Shifts in the relative position of the surface planes was shown to influence the interaction strength too.10 While the simulation results and the behaviour of nanoparticle suspensions seem consistent with each other regarding the observation of adhesive interactions, the experiments on extended surfaces demonstrated on the other hand the existence of repulsive interactions. We believe that the consideration of surface roughness might help to understand these results.
Understanding the role of nanoparticle roughness on the phase behavior of a suspension requires an extension of the existing theoretical models, such as the DLVO theory,11 to incorporate roughness effects. The latter effects have indeed been considered in previous works. Surface and chemical heterogeneities were identified as potential contributors to the discrepancies between theory and experiments,12–14 prompting the implementation of surface roughness in theoretical models.15–21 These works demonstrated that the roughness can lead to interactions that differ from the DLVO potential. The double layer repulsion arising from the calcite surface charge is expected to be much smaller than the solvation forces at nanometer inter-surface separations, given the low surface potential of calcite (−0.02 C m−2 (ref. 22)), and therefore should contribute little at this range of distances. Previous studies have highlighted the impact of surface roughness on the local surface charge, such as steps or terraces, leading to a ten-fold increase with respect to flat calcite surfaces.23 Particle size is also known to modify the surface charge in titanium dioxide nanoparticles.24 However, experiments on calcium carbonate nanoparticles have reported the formation of gel structures,4,8 which support the existence of adhesive interactions. This suggest that the double layer repulsion arising from surface charges may not be dominant at very short intersurface distances.
We adopt in this work a convolution approach to investigate the influence of nanoparticle roughness on the particle interactions. The approach builds on the model developed by Parsons et al.,21 whereby the force between flat surfaces is convoluted with a probability distribution that defines the roughness of the surface of interest. A key input for this approach is the solvent mediated interactions between flat calcite surfaces, which were computed recently using state of the art forcefields and molecular simulations. We use the resulting effective potentials to investigate the phase behavior of colloidal suspensions by means of Langevin Dynamics computer simulations. The surface roughness is shown to have a strong impact on the adhesive minimum found in flat surfaces, which disappears in rough colloids, even for small RMSD roughness of the surface ∼2.0 Å, rendering repulsive interactions dominant. We investigate the interplay of adhesion, DLVO forces and roughness effects in defining the interactions between CN, as a function of the ionic strength of the aqueous solution.
We reproduce in Fig. 1, the solvent mediated interaction free energy of two flat calcite surfaces corresponding to the (104) plane. The free energy features a strong adhesive minimum at r < 1 nm.10 We note that surface hydration can result in long range repulsion. Such effect has been observed in soft interfaces too.25 However, the calcite (104) surface induces the ordering of the interfacial water molecules, which adopt a checkerboard structure. This ordering has been observed both in experiments and computer simulations.10,26–28 When the two surfaces are put in close contact (nanometer distance) the ordering of the water layers spans the whole confined region, leading to strong adhesive force. Changes in the inter-surface distance disrupt the structure of the water layers, resulting in either narrow adhesive minima or repulsion at short distances (see r < 7.5 Å in Fig. 1). Strong adhesive minima have also been reported in mesoscopic models of Calcium Silicate Hydrate nanoparticles studies of cement matrices.29
Fig. 1 Free energy interaction between two atomically flat calcite (104) surfaces. The inset shows the structuring of water confined between two calcium carbonate (104) surfaces in registry at D = 7.5 Å intersurface separation. Results taken from ref. 10. |
In our work we propose that the adhesive forces operate between calcite nanoparticles (several 10's of nm in diameter), where smooth surfaces might be present. A significant amount of hydrated ions at the surfaces could lead to rough surfaces and disrupt the adhesive force. However, previous studies concluded that ion-ion correlations do not play a significant role in the formation of calcite pastes.8 The estimated charge density is fairly low, ∼0.1 charges per nm2. Furthermore, it has been reported that calcite nanoparticles feature smooth surfaces (see TEM images in ref. 8). The existence of such smooth surfaces in these small nanoparticles might well be compatible with the estimated, much longer ∼1 μm, length scales of growth-induced roughening.30
We demonstrate below that small RMS surface roughnesses ρ ∼ 2.0 Å, defined as the standard deviation of the surface height distribution with respect to the average surface height, modify the attractive interaction, leading to strong hydration repulsion shor interparticle distances. Combining the roughened hydration interaction with the DLVO theory then yields effective potentials that can be used to investigate the phase behaviour of the colloidal suspensions. We consider below experimental conditions corresponding to different ionic strengths and pH, which result in different double layer repulsions, which are modelled with the DLVO theory and the Derjaguin approximation.
F(r) = 2πReffW(r) | (1) |
(2) |
Fig. 2 AFM images of a calcite (104) surface in (A) deionized water and (B) supersaturated solution with respect to calcite. The arrows indicate the crystallographic directions of terrace growth. Reprinted with permission from Vavouraki et al., Cryst. Growth Des. 2010, 10(1), 60–69.31 Copyright 2010 American Chemical Society. (C) Schematic representation of our model for a rough nanoparticle. The different surface heights result in an average inter-particle interaction given by the “roughened” Derjaguin approximation in eqn (4), which averages the interaction of different surface heights (see detail of surface in zoomed image) by their probability distribution. (D) Inter-particle interactions for varying concentration of Ca(OH)2, for surface roughness of ρ = 1.5 Å (longer range included in the ESI†). (E) Inter-particle interactions as a function of the particle roughness (increasing ρ from bottom to top), for the zero ion conditions (c = 0 mM). |
To calculate the roughness on the inter-particle interactions we introduce the height function, fh, which is define by the Gaussian distribution:
(3) |
(4) |
This approach has been pioneered by Parsons et al. to incorporate surface roughness in the theoretical calculation of surface forces.21 These authors showed that roughness amplifies the long range behaviour of DLVO forces, and shifts the repulsive branch detected in surface force measurements to a longer distance. The shift scales with the Root Mean Square roughness of the surfaces.
(5) |
(6) |
(7) |
To define the DLVO interactions we use information from the speciation analysis reported in ref. 8 which yields the Debye length and zeta potential of the calcium carbonate surfaces as a function of the initial concentration of Ca(OH)2. The values used in the present study are summarized in Table 1.
c (mM) | I (mM) | λ (nm) | ζ (mV) |
---|---|---|---|
0 | 0.73 | 11.1 | 9.3 |
3 | 0.47 | 13.8 | 10.9 |
15 | 1.7 | 7.3 | 12.5 |
30 | 10.2 | 3.0 | 18.6 |
50 | 43.2 | 1.4 | 18.6 |
We used the viscosity of the fluid to select the damping parameter τ = m/γ required for the thermostat employed in the LD simulations. For water, η = 8.90 × 10−4 Pas at 298 K. Using the definition of the friction coefficient, γ = 3πησ, gives in reduced units which corresponds to 0.8 ns in SI units, using the density of calcite 2710 kg m−3, and nanoparticles of diameter 70 nm.
We used the thermostat for the Langevin equations of motion proposed by Bussi and Parrinello,32 as implemented in LAMMPS.33 Due to the very steep and short range interactions, a very small timestep must be used to ensure accurate integration. To select the simulation time step, we monitored the conservation of effective energy,32 defined as = Etot − ΔEtstat, where Etot is the total energy of the system and ΔEtstat is the increment in the energy due to the thermostat. was monitored for different values of timestep and damping parameters, for the inter-particle interaction corresponding to surfaces with roughness ρ = 1.5 Å of particles with σ = 70 nm. We find good energy conservation for Δt = 1 × 10−5 (see Fig. 1 in the ESI†), in a range of damping parameters 1–100 τwater. Larger damping parameters correspond to lower viscosities, thus effectively increasing the efficiency of the simulation. For this reason, we set τ = 100 τwater. While the change in the effective viscosity of the suspension modifies the dynamics (not investigated in this work), it does not influence the final structure of the colloidal suspension.
The addition Ca(OH)2 has been shown to increase the zeta potential as well as inducing the electrostatic screening of the solution, thereby tuning the DLVO interaction of the colloidal suspension. The results are represented in reduced units, namely U* = U/kBT with T = 298 K, and r* = r/σ, where σ is the nanoparticle diameter.
Relatively low surface roughnesses, ρ = 1.5 Å result in deep energy minima, ∼20–30kBT. This interactions are similar to those considered before in models of adhesive hard-spheres and patchy colloidal potentials,34–38 which included short-ranged attractive wells of the order of 10's kBT. Based on these works, we expect that the effective interaction employed here should lead to irreversible and diffusion-limited cluster aggregation (DLCA),35 characterized by the formation of a gel phase at low particle packing fractions, likely in the interval ϕ = 0.01–0.10.39 Surface roughnesses below 1.5 Å should also result in irreversible aggregation.
Fig. 2D shows that the strongest attraction is obtained in the case of pure calcite. As the concentration of Ca(OH)2 increases to 30 mM, the interaction strength decreases, and for higher concentrations, c = 50 mM, the interaction strength increases again. The effect of surface roughness is shown in Fig. 2E. Increasing the surface roughness results in a shift to longer distances of the repulsive branch associated to the hydration forces, with a concomitant reduction of the interaction strength (given by the potential well depth). The interaction depends strongly on the roughness and for small roughnesses (ρ > 2.0 Å), it becomes purely repulsive, while for smooth surfaces (ρ = 0), adhesive forces due to the solvent structure lead to stable suspensions. Our analysis therefore shows that the conditions required to find a stable CN suspension will depend both on the concentration of Ca(OH)2 (the concomitant change of the pH) and the surface roughness. Experimental studies indicated that at concentrations of the order 30 mM the suspensions are stable, as indicated by the minimization of the storage modulus.8
Depending on the surface roughness, the CN formed a stable suspension (larger roughness ρ > 2.0 Å, see Fig. 3A) or aggregated into fractal clusters (smaller roughness ρ < 2.0 Å, see Fig. 3B and C). The interplay of the roughness, the adhesive and DLVO forces led to the formation of compact spherical aggregates (see Fig. 3D). The latter appear at high roughness in a wide range of volume fractions ϕ. Intermediate roughnesses favour the formation of fractal and percolating clusters (see Fig. 3C).
To quantify the nanostructure of the nanoparticle suspensions we computed the mean cluster size, employing a distance criterion to construct the clusters. Two nanoparticles i and j were assigned to the same cluster if their distance dij < 1.02 ma. The attractive minima are contained within this distance, as seen in Fig. 2D and E, and thus corresponds to the characteristic distance between aggregated nanoparticles. The mean cluster size distribution, s, was computed using the equation:40,41
(8) |
(9) |
The fractal dimension of the suspension was calculated using clusters containing more than 20 colloids, with a similar expression to eqn (8):
(10) |
Following the work by Griffiths et al., we calculate the fractal dimension at two scales. The first scale corresponds to box sizes with a size similar to the characteristic size of the particle (l = [1σ,5σ]), hence probing the local structure of the cluster (local fractal dimension), and characterizing the degree of compactness of individual clusters. A second scale corresponds to box counting cell sizes in the range l = [5σ,L], where L is the size of the simulation box. This length scale quantifies the amount of volume fraction occupied by each cluster, and therefore it is related with the percolation of the clusters (global fractal dimension). See ESI† for more information with respect to local and global fractal dimensions.
We show in Fig. 4 a colour map that illustrates the range of mean cluster size and fractal dimensions in the roughness/packing fraction plane. The mean cluster size as well as global the fractal dimension show that the cluster fill the space more efficiently at high volume fractions (large fractal dimension), and form percolating structures with a percolation threshold between ϕ = 0.04–0.05, as evidenced by an increase in the mean cluster size 〈s〉 converging to ∼N, the number of particles in the simulation cell. For volume fractions ϕ > 0.05 and roughness values ∼2.0 Å, the structures transition from a percolating gel to a fluid phase. This region of the phase diagram is characterized by an increase in the local fractal dimension, df ∼ 1.8. The increase in the fractal dimension reflects an increase in the degree of compaction of the clusters, which arises from weaker attractive interaction. The reduction in the interactions allows the clusters to rearrange themselves into more compact structures, specifically crystalline structures, as shown in Fig. 3D. Aggregation into local crystalline structures as opposed to the formation of percolating gels, has also been observed by Griffiths and coworkers39 using the Morse pair potential with different interaction strengths. Those authors reported local fractal dimensions similar to the ones we find here. Similar fractal dimensions have also been reported in experimental studies of gold colloidal aggregates formed via irreversible kinetic aggregation.45 This supports our view that the percolating phases form following diffusion-limited cluster aggregation.
Our simulations show that at fixed roughness, the addition of the Ca(OH)2 and therefore the increase in the solution pH, results in non-monotonic changes of the fractal dimension with volume fraction. For instance we observe a minimum in the fractal dimension at 30 mM and ρ = 1.7, which reflects the formation of a stable suspension, while at this roughness most of the suspensions at other concentrations feature adhesive behaviour.
Our results show that nanoparticle surface roughness has a large impact on the inter-colloidal interactions and on the stability of the suspensions. Our model predicts aggregation and gel formation (as reported in experiments of calcite nanoparticles8) when the surface roughness is small (ρ < 2.0 Å). The fractal dimension of the clusters obtained in this way is comparable to the diffusion-limited cluster aggregation fractal dimension of ∼1.75,42 where particles aggregate when undergoing a random walk, without further reorganization of the structure upon joining the nanoparticle cluster. Increasing the surface roughness leads to a shift of the effective repulsion to longer inter-particle distances and to the stabilization of the suspension. We have observed a synergistic effect between nanoparticle roughness and DLVO forces (mediated by changes in the concentration of Ca(OH)2). At conditions corresponding to high base concentrations, 50 mM of Ca(OH)2, attractive interactions may be important, even for relatively large roughness, leading to the aggregation of the suspensions. At intermediate concentrations, 30 mM, the suspensions are stable, even for relatively low surface roughness, and the suspensions features a minimum in the roughness required for gel formation. Our model predicts that at lower concentrations, <30 mM, gel phases are formed for small roughness. High pH conditions and high ionic concentration (high Ca(OH)2) may disrupt the particle morphology inducing an increase in the roughness. According to our model such effects would eliminate the adhesive minimum, making aggregation less likely. The lack of gel phases at 30 mM concentration is consistent with the experimental measurements of calcium carbonate nanoparticle suspensions, which reported a minimum in the storage modulus at this concentration,8 signalling the formation of a fluid phase. Our model indicates that this phase can be formed because the nanoparticles are sufficiently rough (according to our model ρ > 1.5 Å), since roughness below this value result in the formation a gel phase. This result should motivate additional experimental analyses to resolve the surface structure of calcite nanoparticles.
Overall, our result highlight the importance of the surface structure of nanoparticles on the phase behavior of suspensions. We provide theoretical predictions of the relationship between the stability of the suspension and the inter-colloidal energy profiles, with particular emphasis on the modification of the adhesive minimum. We have illustrated the interplay between roughness and adhesive effects in a models of calcite nanoparticle suspensions, which are relevant in industrial applications, building materials and biomineralization processes. Our results indicate that small changes in the surface roughness of the nanoparticles (of the order of calcite atomic steps) might modify significantly the interparticle forces. Indeed, we demonstrate that attraction or repulsion between the surfaces can be observed with the same underlying interactions, when the roughness of the surface is taken into account. We expect that our work will serve to rationalize existing contrasting observations, as well as to connect the phase behavior of suspensions to the surface topography and volume fraction of the suspension. The phase diagram reported in this work (Fig. 3E) provides a route to establish this connection. While, we have focused on calcite, by considering a combination of interactions obtained from atomistic simulations and experimental studies, the theoretical approach discussed here to predict the relationship between the stability of the suspension and the inter-colloidal interactions, could be extended to other colloidal suspensions.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0nr00834f |
This journal is © The Royal Society of Chemistry 2020 |