Sohaib
Mohammed
a,
Ivan
Kuzmenko
b and
Greeshma
Gadikota
*a
aSchool of Civil and Environmental Engineering, Cornell University, Ithaca, New York 14853, USA. E-mail: gg464@cornell.edu; Tel: +1 607-255-4796
bAdvanced Photon Source, Argonne National Laboratory, Lemont, Illinois 60439, USA
First published on 29th November 2021
Achieving reversible and tunable assembly of silica nanoparticles at liquid–liquid interfaces is vital for a wide range of scientific and technological applications including sustainable subsurface energy applications, catalysis, drug delivery and material synthesis. In this study, we report the mechanisms controlling the assembly of silica nanoparticles (dia. 50 nm and 100 nm) at water–heptane and water–toluene interfaces using sodium dodecyl sulfate (SDS) surfactant with concentrations ranging from 0.001–0.1 wt% using operando ultrasmall/small-angle X-ray scattering, cryogenic scanning electron microscopy imaging and classical molecular dynamics simulations. The results show that the assembly of silica nanoparticles at water–hydrocarbon interfaces can be tuned by controlling the concentrations of SDS. Silica nanoparticles are found to: (a) dominate the interfaces in the absence of interfacial SDS molecules, (b) coexist with SDS at the interfaces at low surfactant concentration of 0.001 wt% and (c) migrate toward the aqueous phase at a high SDS concentration of 0.1 wt%. Energetic analyses suggest that the van der Waals and electrostatic interactions between silica nanoparticles and SDS surfactants increase with SDS concentration. However, the favorable van der Waals and electrostatic interactions between the silica nanoparticles and toluene or heptane decrease with increasing SDS concentration. As a result, the silica nanoparticles migrate away from the water–hydrocarbon interface and towards bulk water at higher SDS concentrations. These calibrated investigations reveal the mechanistic basis for tuning silica nanoparticle assembly at complex interfaces.
Furthermore, advances in the use of engineered silica nanoparticles to architect new matter for applications related to catalysis,4,5 drug delivery,6,7 materials synthesis,8,9 separation,10,11 sensors,12 adsorption,13 enhanced oil recovery14 and food processing15 motivate us to develop mechanistic insights into the assembly of silica nanoparticles at complex fluidic interfaces. For example, Pickering emulsions have been used as polymerization vessels for preparing hollow microcapsules and composites capsules.16 These materials are widely used in food, pharmaceutical and cosmetic products. Advancing fundamental insights into the interfacial behavior and the energetic interactions underlying the organization of silica nanoparticles and surfactants at complex interfaces is crucial for synthesizing novel materials. With this context, the following questions are addressed in this study: (i) What is the influence of surfactant concentrations on the organization or displacement of silica nanoparticles at water–hydrocarbon interfaces? (ii) How does the interfacial behavior of surfactants and silica nanoparticles evolve with the size of nanoparticles? (iii) What are the energetic interactions underlying the organization or displacement of silica nanoparticles away from water–hydrocarbon interfaces?
Prior experimental and computational studies utilized various techniques, including Pendant drop tensiometer and molecular dynamics simulations17–22 to probe the interfacial behavior of surfactants and nanoparticles at various water–hydrocarbon interfaces. While some studies concluded that the synergetic interactions of nanoparticles and surfactants on the interface enhance emulsion stability,22–27 others reported that the surfactant concentrations higher than the critical micelle concentration (CMC) leads to the displacement of the nanoparticles from the interface.28,29 Although it is poorly understood, the mechanism governing the displacement of nanoparticles from the interface is related to the interfacial tension of the associated emulsion and the contact angle between the nanoparticles and the liquid phases.17
The energy required to displace a nanoparticle from the interface can be quantified as E = πr2γOW(1 ± cosφ)2 where r is the radius of the nanoparticle, γOW is the interfacial tension of the oil-water emulsion, and φ is the contact angle of the adsorbed nanoparticles at the interface.30 The extent to which the nanoparticles adsorb at the interface depends on the wettability of these nanoparticles. Nanoparticles with highly hydrophilic or hydrophobic surfaces tend to disperse in water or hydrocarbon phases, respectively. According to the equation above, the maximum energy required to displace a nanoparticle from the interface is at a contact angle of 90°. The synergetic interactions of surfactants and nanoparticles on the interface can alter the wetting properties of the assembled nanoparticles and eventually enhance or inhibit the assembly of nanoparticles on the interface. Vashisth and co-workers29 showed that adding high surfactant concentrations can completely disrupt the assembly of nanoparticles from the interface because the presence of the surfactant at the water–oil interface is energetically favored compared to nanoparticles. This displacement was also attributed to the influence of the added surfactants on the nanoparticle's contact angle at the interface,31,32 the extent of nanoparticles flocculation,33,34 the dynamics of nanoparticles at the interface35 and the competitive assembly of surfactant molecules at the interface.36 Resolving the synergetic effects of surfactant and nanoparticles enables us to tune the self-assembly behavior of nanoparticles at water–hydrocarbon interfaces in addition to the ability to control the interfacial properties of the associated emulsion including the coalescence of small drops into larger drops that facilitate the separation of water and oil phases.
In this study, we introduce an experimental and computational approach to elucidate the interfacial behavior underlying the co-existence of silica nanoparticles with diameters of 50 nm and 100 nm and sodium dodecyl sulfate (SDS) concentrations ranging from 0.001–0.1% at water–heptane and water–toluene interfaces. Ultrasmall- and Small-Angle X-Ray Scattering (USAXS/SAXS) measurements, Cryogenic Scanning Electron Microscopy (Cryo-SEM) imaging, and Molecular Dynamics (MD) simulations are harnessed to probe the coexistence of hydrophilic silica nanoparticles and SDS molecules on the interfaces and their influence on the interfacial properties of the studied emulsions. The interfacial behavior of the silica nanoparticles is investigated at SDS concentrations of 0.001, 0.01 and 0.1 wt%. The experiments and simulations are performed at 1 bar and 298 K. This experimental and simulation approach enables us to develop fundamental insights into the structure, dynamics and morphology of nanoparticles and surfactant assembled at water–hydrocarbon interfaces across spatial and temporal scales.
The nonbridging surface oxygens on the nanoparticle's surface are functionalized with OH groups with a surface density of about 7.3 OH nm−2. The rationale behind choosing this level of hydroxylation is driven by the chemistry of the nanoparticles used in the experiments which are dominated by silanol groups (Si–OH) on the surface and the typical density of OH groups on the silica surfaces used in previous studies.43–45 Further, our reactive molecular dynamics simulations (ReaxFF) shows that nonbridging oxygens on the surface of 3 nm sized silica nanoparticles undergo continuous hydroxylation when dissolved in 1 g cm−3 of water at 298 K and 1 bar (see Fig. S1 and S2†).
Five silica nanoparticles are randomly distributed in a 7 nm × 7 nm × 10 nm cell with a concentration of 1 wt% to match the experimental concentration of the silica nanoparticles suspensions. The simulation cell is extended in the z-direction by 5 nm and SDS–heptane and SDS–toluene mixtures are added with surfactant concentrations of 0 wt%, 0.001 wt%, 0.01 wt%, and 0.1 wt% (see Fig. 1(b)). The dimensions of the initial simulation cells are 7 nm × 7 nm × 15 nm in x, y and z directions, respectively. Periodic boundary conditions are used in the x, y and z directions.
Silica nanoparticles and water molecules are modeled using parameters from CLAYFF46 and TIP4P,47 respectively, while heptane, toluene and SDS molecules are modeled using the OPLS/AA forcefield48,49 (see the forcefield parameters in Table S1†). The combination of OPLS/AA and CLAYFF forcefields have been used extensively in investigating the interactions of hydrocarbons and surfactants with silica surfaces.50–52 50000 steps of energy minimization using “steepest descent” method are performed on the initial configurations to decrease their high energy and remove the inappropriate geometries. The minimized configurations underwent semi-isotropic compression in z-direction, the direction normal to the interface for 50 ns using the constant number of molecules, constant pressure, constant interfacial area, and constant temperature (NPAT) ensemble. The simulations are performed at 1 bar and 298 K to be consistent with the experimental conditions. Nose–Hoover thermostat53,54 with a relaxation time of 1 ps and Parrinello–Rahman barostat55 are used to maintain the required temperature and pressure. The equation of motion is integrated using the leapfrog algorithm with a time step of 1 fs. The short-range interactions are calculated within a cutoff of 1.4 nm, while the long-range electrostatic interactions are treated using Particle Mesh Ewald (PME).56 The nonbonded van der Waals and electrostatic interactions are modeled using 12-6 Lennard-Jones and coulombic models, respectively. The bonded interactions account for bonds stretching, angles bending and dihedrals, except for silica nanoparticles where only OH bond stretching are accounted for. All the simulations are conducted using GROningen Machine for Chemical Simulations (GROMACS 2018) simulation package.57
As 0.1 wt% SDS is added to the system, the water–heptane and water–toluene interfaces are completely occupied by SDS molecules and no scattering from nanoparticles was observed, suggesting that the interfacial nanoparticles are completely displaced from the interfaces and emerged in the water phase. The competitive assembly of silica nanoparticles and SDS molecules on water–hydrocarbon interfaces is consistent with prior experimental work showing that the small quantities of surfactant stabilized nanoparticle emulsions.29,59
To determine the size distribution and the number of interfacial silica nanoparticles at water–heptane and water–toluene interfaces, we modeled the USAXS/SAXS curves that represent the scattering behavior of silica nanoparticles and surfactant at these interfaces. Quantitative analyses on the nanoparticles’ dominated interfaces are performed by modeling the structure factor and the form factor of the intensity curves on the Q scale. The intensity scattered from identical spherical nanoparticles can be expressed as the product of and such that:60
(1) |
The size distribution and the number of interfacial silica nanoparticles at water–heptane and water–toluene interfaces (Fig. 3) are obtained from fitting at the high Q regions. At 0 wt% SDS, the number of 50 nm and 100 nm nanoparticles assembled at water–toluene surfaces is higher compared to water–heptane interfaces. The number of 50 nm and 100 nm silica nanoparticles at water–heptane and water–toluene interfaces decreased substantially as the SDS concentration increased to 0.001 wt%, indicating that the added SDS molecules preferentially occupy the interfacial area and thus reduce the available area for the nanoparticles to adsorb on the interface. It is interesting to note that the number of interfacial silica nanoparticles at water–heptane interface remains similar when the concentrations of surfactants are 0.001 wt% and 0.01 wt% irrespective of the size of the nanoparticles (Fig. 3(a) and (c)). In contrast, the number of silica nanoparticles at the water–toluene interface dropped to zero on increasing the concentration of surfactants to 0.01 wt% and 0.1 wt% (Fig. 3(b) and (d)). These observations agree with the findings of Vashisth and co-workers that interfacial surfactants displace the adsorbed nanoparticles after accounting for the difference in the sizes of nanoparticles.29
The observations noted from Fig. 3 are consistent with the scattering data shown in Fig. 2. The scattering pattern from water–heptane and water–toluene interfaces at high SDS concentrations (i.e., 0.1 wt% at water–heptane and ≥ 0.01 wt% at water–toluene) suggests the absence of silica nanoparticles at these interfaces and the dominance of surfactant assemblies (see Fig. 2). At 0.1 wt% SDS, the scattering curves from the water–heptane interfacial regions showed peaks at about 0.05 Å−1 indicating scattering from self-assembled surfactant molecules. The scattering from self-assembled surfactant at water–toluene are noticed at 0.01 wt% and at 0.1 wt% SDS concentrations, whereas peaks are detected at about 0.031 Å−1 and 0.026 Å−1, respectively. The locations of the peaks in water–heptane and water–toluene interfaces show the dependence of the formed assembles sizes and morphologies on the interfacial properties of the corresponding emulsions. Porod analysis of the interfacial surfactants on water–heptane and water–toluene interfaces are performed at high Q region to investigate the morphology of the self-assembled SDS molecules (see Fig. 2). The power law slope values of the scattering intensity curves on water–heptane interfaces in the presence of 0.1 wt% are 2.4 and 2.3 for emulsions with 50 nm and 100 nm silica nanoparticles, respectively. These values indicate the formation of 2D disk-like morphologies by the interfacial surfactants. However, SDS molecules at the interfaces of water–toluene forms both 2D disk-like (Q−2) and 1D rod-like (Q−1) morphologies. The formation of surfactant ensembles on water–hydrocarbon interfaces are evident from prior experimental and computational studies.61–64 Czajka and Armes61 concluded that the presence of SDS surfactants results in the formation of spherical latex emulsions based on SAXS measurements. Shi and Guo64 performed MD simulations on the effect of SDS molecular structure on its aggregation morphology at water–trichloroethylene interfaces. The results indicated that the surfactants form a continuous monolayer at the interface as the interfacial coverage increases.
The morphology of the self-assembled silica nanoparticles and SDS ensembles are further analyzed by calculating the power law slope in the USAXS (low Q) region (see Fig. 4). The morphology of self-assembled nanoparticles and SDS molecules change significantly with the SDS concentration. The low power law slope (< 1) in the absence (0 wt%) and presence of low SDS concentrations (0.001 wt%) indicate that the interfacial nanoparticles and SDS molecules are well dispersed at the interfaces. As the SDS concentration increases, the power law slope values increase substantially, implying the emergence of SDS aggregates with varying morphologies. A power law slope of < 3 indicates the formation of surface fractals.65,66 Further, the displaced silica nanoparticles to the bulk water are well dispersed with no indication of self-assembly as the power low slope is < 1 in the USAXS region (see Fig. S3†).
The displacement of interfacial nanoparticles from water–heptane and water–toluene interfaces are confirmed by cryo-SEM imaging (see Fig. 5). At low SDS concentrations of 0.001 wt%, the nanoparticles are abundant in the interfacial region of water–heptane and water–toluene emulsions. The nanoparticles are mostly found as individual nanoparticles with minor traces as aggregates on water–heptane and water–toluene interfaces (see Fig. 5(a) and (c)). However, the number of interfacial nanoparticles on water–heptane and water–toluene interfaces were significantly reduced as the SDS concentration increased to 0.1 wt%. Only traces of nanoparticles are observed at the interfaces due to the displacement into the aqueous phase by the interfacial SDS assembles. The cryo-SEM data agrees with the observation of USAXS/SAXS scattering measurements (Fig. 2).
Additional insights into the assembly of silica nanoparticles and SDS molecules on water–heptane and water–toluene interfaces are obtained from the density profiles governed by MD simulations (Fig. 6). The density profiles are averaged over the last 10 ns of the simulation time and calculated along the axis normal to the interface (z-axis). Water and hydrocarbons formed two separated phases, while the silica nanoparticles and SDS molecules are dispersed in the water and hydrocarbon phases, respectively, with favorable assembly at the interfaces (Fig. 6(a) and (b)). The initially dispersed nanoparticles in the water phase migrated and assembled on the interfaces to extents that differ with the concentration of SDS molecules. SDS density profiles showed peaks at water–heptane and water–toluene interfaces, indicating preferential assembly at water–hydrocarbon interfaces compared to the dispersion in bulk hydrocarbons. Higher peaks in SDS profiles are observed on water–toluene interfaces compared to water–heptane interfaces.
Interestingly, the density profiles of silica nanoparticles at water–heptane interfaces showed a systematic shift toward the water phase as the SDS concentrations increase in the system (Fig. 6(c)). The nanoparticles’ density peaks at water–heptane interfaces are shifted from 8.7 nm to 7.6 nm as the SDS concentration increased from 0 wt% to 0.1 wt%, respectively. Similarly, the density peaks of nanoparticles at water–toluene are shifted from 8.5 nm to 8.3 nm as the SDS concentrations increased from 0 wt% to 0.1 wt%, respectively. The shift toward the aqueous phase suggests that the assembled SDS molecules on the interfaces displace the assembled nanoparticles away from the interface. This observation is consistent with the USAXS/SAXS intensities (Fig. 2 and 3) and agrees with previous studies.29
The displacement of nanoparticles from the interface is also evident from the density of water molecules around the surface of nanoparticles. These densities are quantified by calculating the radial distribution function (g(r)) of oxygen atoms in water molecules from the hydroxyl groups (–OH) on the surface of silica nanoparticles as a reference atom (see Fig. 7). The peaks represent the first coordination shells and second coordination shells of the oxygen in the water molecule around silica –OH groups. The density of the first and second coordination shells increased systematically with the increase in SDS concentrations. The denser water shells stem from the exposure of nanoparticle surface to water molecules that cover the nanoparticle as it is displaced from the interface and emerges into the aqueous phase to solvate the surfaces of nanoparticles and reduce its surface free energy.
(2) |
In this expression, γWH is the interfacial tension of water–heptane and water–toluene mixtures, n is the number of interfaces in the simulation cell, Pii is the pressure (where i = x, y and z directions) and Lz is the equilibrium cell length in the direction normal to the interface (z direction).
The calculated interfacial tensions of water–heptane and water–toluene at 1 bar and 298 K are 51.8 mN m−1 and 36.4 mN m−1, respectively, which are consistent with prior experimental and simulation results.67,71 Introducing silica nanoparticles to the interface resulted in a 4.7% and 9.4% reduction in the interfacial tension of water–heptane and water–toluene mixtures, respectively. Increases in the SDS concentrations to 0.01 wt% resulted in a slight reduction in the interfacial tension of both emulsions due to the addition of limited quantities of surfactant molecules that accumulate at the water–hydrocarbon interfaces. A substantial reduction in the interfacial tensions of both emulsions is observed on adding 0.1 wt% of SDS to the systems. The interfacial tension of water–heptane and water–toluene decreased to 24.3 and 9.1 mN m−1, respectively. These results suggest that high concentrations of SDS are more effective in stabilizing the emulsion of interest compared to hydrophilic silica nanoparticles. The simulations and experiments reveal that the stabilization of water–hydrocarbon emulsions in the presence of high concentrations of SDS surfactants is accompanied by limiting the assembly of silica nanoparticles at water–hydrocarbon interfaces.
The mechanism by which silica nanoparticles and SDS molecules reduce the interfacial tension of water–heptane and water–toluene emulsions involves the assembly and dispersion of these nanoparticles and surfactant at the interfaces (see Fig. 9). Silica nanoparticles form a cluster of five nanoparticles at water–heptane and water–toluene in the absence of SDS molecules. Interestingly, the aggregation number of the nanoparticles decreases systematically as the SDS concentration increases (Fig. 9A and B) due to the preferential assembly of the SDS molecules at the interfacial region (Fig. 9C and D). Thus, the reduction in the interfacial tensions is driven by three scenarios (Fig. 9E): (a) the presence of nanoparticles at the interface results in slight reduction in the interfacial tension, (b) the coexistence of nanoparticles and surfactant at the interface results in a greater reduction in the interfacial tension compared to the first scenario, and (c) the excess SDS molecules at the interfaces of water–heptane and water–toluene reduces the interfacial tension significantly.
Interestingly, the interactions of nanoparticles with water in water–heptane and water–toluene systems increase systematically with the SDS concentrations. The increase in nanoparticles–water interactions is associated with a systematic decrease in the magnitude of nanoparticle–heptane and nanoparticle–toluene interaction energies. The trends of nanoparticles–water and nanoparticles–hydrocarbons energetic interactions as SDS concentrations increase confirm the displacement of nanoparticles from the interfaces toward the aqueous phase by SDS molecules. The displacement of the nanoparticles toward the aqueous phase is also confirmed by the increase of the average number of water molecules in the first coordination shell of the hydroxyl group on the surface of the nanoparticles (see Table 1). The number of Owater in the first coordination shell of OHsilica increases systematically with the increase of SDS concentrations in the hydrocarbon phases, indicating the increase of water density around the displaced nanoparticles. The enhanced solvation of silica nanoparticles due to the displacement into the aqueous phase is associated with an increase in the number of the hydrogen bonds between the hydroxyl groups on the surface of the nanoparticles and the surrounding water molecules (see Fig. S4†).
SDS concentration (wt%) | Water–heptane | Water–toluene |
---|---|---|
0 | 0.43 ± 0.01 | 0.43 ± 0.02 |
0.001 | 0.45 ± 0.04 | 0.44 ± 0.01 |
0.01 | 0.47 ± 0.02 | 0.47 ± 0.03 |
0.1 | 0.48 ± 0.02 | 0.54 ± 0.01 |
van der Waals interactions contribute significantly to the coexistence of silica nanoparticles and SDS surfactant at the interfaces of water–heptane and water–toluene emulsions, with electrostatic interactions having minor contribution (see Table 2). van der Waals and electrostatic interactions increase systematically with the concentrations of SDS due to the assembly of SDS at the interfaces. Further, van der Waals and electrostatic interactions between the interfacial nanoparticles and surfactants are higher at water–toluene compared to at water–heptane at all SDS concentrations. The higher intermolecular interactions of nanoparticles–SDS at the interface of water–toluene explains the significant reduction in the interfacial tension driven by the synergetic interactions between the interfacial nanoparticles and surfactants. However, the increase in the hydrogen bonding of nanoparticles–water (Fig. S4†) and SDS–water (Fig. S5†) with the increase of SDS concentrations confirms the favorable adsorption of SDS surfactant at the interface over the silica nanoparticles.
SDS concentration (wt%) | Water–heptane interface | Water–toluene interface | ||
---|---|---|---|---|
van der Waals | Electrostatic | van der Waals | Electrostatic | |
0.001 | −2.52 ± 0.04 | −0.43 ± 0.02 | −2.60 ± 0.01 | −0.45 ± 0.04 |
0.01 | −16.10 ± 0.12 | −2.41 ± 0.08 | −23.05 ± 0.27 | −3.35 ± 0.09 |
0.1 | −43.25 ± 1.39 | −5.19 ± 0.42 | −69.207 ± 2.03 | −11.02 ± 0.21 |
Interestingly, attractive electrostatic interactions of nanoparticles with water are noted while the van der Waals interactions between the OH-terminated nanoparticles and water molecules are repulsive (see Table 3). Hydrogen bonding between the hydroxyl groups on the nanoparticles surfaces and water molecules contributed positively to the strong attraction between the nanoparticles and the aqueous phase. Nanoparticle interactions with heptane are dominated by van der Waals interactions with minor contributions of Coulomb attraction, while the interaction with toluene includes a substantial contribution from van der Waals and electrostatic interactions. The interaction energies elucidate the systematic displacement of silica nanoparticles from water–heptane and water–toluene interfaces that observed from USAXS/SAXS measurements, cryo-SEM imaging and MD simulations.
SDS concentration (wt%) | Silica nanoparticles–water interactions | Silica nanoparticles–hydrocarbon interactions | ||
---|---|---|---|---|
van der Waals | Electrostatic | van der Waals | Electrostatic | |
Water–heptane | ||||
0 | 1435.4 ± 11 | −18805.4 ± 49 | −1150.67 ± 11 | −45.46 ± 5 |
0.001 | 1461.8 ± 12 | −20155.9 ± 48 | −874.658 ± 16 | −33.65 ± 4 |
0.01 | 1489.8 ± 19 | −20495.0 ± 31 | −748.29 ± 8 | −26.18 ± 1 |
0.1 | 1499.3 ± 20 | −20545.6 ± 19 | −634.7 ± 12 | −22.14 ± 2 |
Water–toluene | ||||
0 | 1378.1 ± 10 | −19614.7 ± 53 | −884.7 ± 14 | −970.2 ± 12 |
0.001 | 1454.0 ± 23 | −20135.9 ± 49 | −774.7 ± 18 | −953.4 ± 11 |
0.01 | 1476.7 ± 23 | −20224.0 ± 51 | −648.3 ± 13 | −951.1 ± 13 |
0.1 | 1527.3 ± 16 | −21426.6 ± 21 | −534.7 ± 9 | −353.4 ± 6 |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1nr06807e |
This journal is © The Royal Society of Chemistry 2022 |