Curvature dependence of the effect of ionic functionalization on the attraction among nanoparticles in dispersion.

Solubilization of nanoparticles facilitates nanomaterial processing and enables new applications. An effective method to improve dispersibility in water is provided by ionic functionalization. We explore how the necessary extent of functionalization depends on the particle geometry. Using molecular dynamics/umbrella sampling simulations, we determine the effect of the solute curvature on solvent-averaged interactions among ionizing graphitic nanoparticles in aqueous dispersion. We tune the hydrophilicity of molecular-brush coated fullerenes, carbon nanotubes, and graphane platelets by gradually replacing a fraction of the methyl end groups of the alkyl coating by the ionizing -COOK or -NH3Cl groups. To assess the change in nanoparticles' dispersibility in water, we determine the potential-of-mean-force profiles at varied degrees of ionization. When the coating comprises only propyl groups, the attraction between the hydrophobic particles intensifies from spherical to cylindrical to planar geometry. This is explained by the increasing fraction of surface groups that can be brought into contact and the reduced access to water molecules, both following the above sequence. When ionic groups are added, however, the dispersibility increases in the opposite order, with the biggest effect in the planar geometry and the smallest in the spherical geometry. These results highlight the important role of geometry in nanoparticle solubilization by ionic functionalities, with about twice higher threshold surface charge necessary to stabilize a dispersion of spherical than planar particles. At 25%-50% ionization, the potential of mean force reaches a plateau because of the counterion condensation and saturated brush hydration. Moreover, the increase in the fraction of ionic groups can weaken the repulsion through counterion correlations between adjacent nanoparticles. High degrees of ionization and concomitant ionic screening gradually reduce the differences among surface interactions in distinct geometries until an essentially curvature-independent dispersion environment is created. Insights into tuning nanoparticle interactions can guide the synthesis of a broad class of nonpolar nanoparticles, where solubility is achieved by ionic functionalization.

Extent of Surface Force Additivity on Chemically Heterogeneous Substrates at Varied Orientations.

Surface interactions between chemically mixed surfaces, as well as those among dissolved biomolecules, comprise distinct contributions from polar and hydrophobic moieties. These contributions are often context dependent. Approximate compliance to the Cassie additivity equation for the wetting free energies on mixed surfaces in water is, however, indicative of similarly additive forces between individual surface elements, suggesting a quadratic interpolation model for total force from the forces between pure surfaces. We use molecular dynamics/umbrella sampling simulations of parallel and nonparallel mixed surfaces with demonstrable Cassie-like behavior to verify how well the total surface force between the heterogeneous, molecularly rough surfaces can be approximated as a combination of forces among the homogeneous ones. When accounting for dissimilar distances of approach between functional groups of different types, our results for graphene surfaces with mixed methyl and nitrile coating show such a superposition to provide a reasonable first order approximation of interactions between the platelets. Deviations from additivity are more prominent in parallel-plate configurations, at high content of hydrophobic groups, and small separations. The inclusion of water polarizability does not visibly alter the observed behavior regardless of platelet orientations. The outcome of this study determines the necessary molecular conditions for observing force additivity that emphasize the context dependence of hydrophobic interaction in the presence of polar groups. This notion provides guidelines for the syntheses of new, chemically heterogeneous materials with tailored function-oriented properties in aqueous media.

Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes.

In addition to the direct attraction, sizable hydrophobes in water experience an attractive force mediated by interfacial water. Using simple geometric arguments, we identify the conditions at which the water-induced interaction between curved hydrocarbon surfaces becomes repulsive. The repulsive contribution arises from the thermodynamic penalty due to the emergence of the liquid/vapor boundary created as water gets expelled between curved hydrophobes. By augmenting the mean field approach with atomistic simulations of pristine and alkyl-coated graphitic nanoparticles in three distinct geometries, spherical, cylindrical and planar, immersed in water, we show the macroscopic thermodynamics remarkably works down to the molecular scale. The new insights improve the prediction and control of wetting and dispersion properties for a broad class of nonpolar nanoparticles.

Fluctuation-driven anisotropy in effective pair interactions between nanoparticles: thiolated gold nanoparticles in ethane.

Fluctuations within the ligand shell of a nanoparticle give rise to a significant degree of anisotropy in effective pair interactions for low grafting densities [B. Bozorgui, D. Meng, S. K. Kumar, C. Chakravarty, and A. Cacciuto, Nano Lett. 13, 2732 (2013)]. Here, we examine the corresponding fluctuation-driven anisotropy for gold nanocrystals densely passivated with short ligands. In particular, we consider gold nanocrystals capped by alkylthiols, both in vacuum and in ethane solvent at high density. As in the preceding study, we show that the anisotropy in the nanoparticle pair potential can be quantified by an angle-dependent correction term to the isotropic potential of mean force (PMF). We find that the anisotropy of the ligand shells is distance dependent, and strongly influenced by ligand interdigitation effects as well as expulsion of ligand chains from the interparticle region at short distances. Such fluctuation-driven anisotropy can be significant for alkylthiol-coated gold nanoparticles, specially for longer chain lengths, under good solvent conditions. The consequences of such anisotropy for self-assembly, specially as a function of grafting density, solvent quality and at interfaces, should provide some interesting insights in future work. Our results clearly show that an isotropic two-body PMF cannot adequately describe the thermodynamics and assembly behavior of nanoparticles in this dense grafting regime and inclusion of anisotropic effects, as well as possibly many-body interactions, is necessary. Extensions of this approach to other passivated nanoparticle systems and implications for self-assembly are considered.

Relating structure, entropy, and energy of solvation of nanoscale solutes: application to gold nanoparticle dispersions.

Structural estimators for the entropy are combined with an analysis of the different contributions to the energy of solvation to understand the molecular basis of the thermodynamics of solvation of passivated nanoparticles. Molecular dynamics simulations of thiolated gold clusters in ethane are performed over a wide range of densities close to the critical isotherm. The entropic changes associated with solvent reorganization around the passivated nanoparticle are estimated from the nanoparticle-solvent pair correlation function, while the entropy of the ligand shell is estimated from the covariance in the positional fluctuations of the ligand atoms. The ligand-shell entropy (S(L)) is shown to be fairly insensitive to variations in solvent density ranging from vacuum to twice the critical density (ρ(c)). In contrast, the entropy change due to solvent reorganization (ΔS(ns)(ord)) shows a minimum around the critical point where the solvent excess shows a maximum. Combining the entropic estimates with the nanoparticle-solvent interaction energies, the free energy of solvation is shown to decrease with density once the critical point is crossed in a manner qualitatively consistent with available experimental data. The results suggest that such an approach to obtain structural insights into the thermodynamics of solvation of passivated nanoparticles could be useful in understanding the stability of nanoparticle dispersions of widely varying chemistries. This study also demonstrates that the theoretical analysis of solvation and self-assembly developed in the context of biomolecular hydration can be very usefully extended to understand the behavior of inorganic nanoparticle dispersions.

Water and other tetrahedral liquids: order, anomalies and solvation.

In order to understand the common features of tetrahedral liquids with water-like anomalies, the relationship between local order and anomalies has been studied using molecular dynamics simulations for three categories of such liquids: (a) atomistic rigid-body models for water (TIP4P, TIP4P/2005, mTIP3P, SPC/E), (b) ionic melts, BeF(2) (TRIM model) and SiO(2) (BKS potential) and (c) Stillinger-Weber liquids parametrized to model water (mW) and silicon. Rigid-body, atomistic models for water and the Stillinger-Weber liquids show a strong correlation between tetrahedral and pair correlation order and the temperature for the onset of the density anomaly is close to the melting temperature. In contrast, the ionic melts show weaker and more variable degrees of correlation between tetrahedral and pair correlation metrics, and the onset temperature for the density anomaly is more than twice the melting temperature. In the case of water, the relationship between water-like anomalies and solvation is studied by examining the hydration of spherical solutes (Na(+), Cl(-), Ar) in water models with different temperature regimes of anomalies (SPC/E, TIP4P and mTIP3P). For both ionic and nonpolar solutes, the local structure and energy of water molecules is essentially the same as in bulk water beyond the second-neighbour shell. The local order and binding energy of water molecules are not perturbed by the presence of a hydrophobic solute. In the case of ionic solutes, the perturbation is largely localized within the first hydration shell. The binding energies for the ions are strongly dependent on the water models and clearly indicate that the geometry of the partial charge distributions, and the associated multipole moments, play an important role. However the anomalous behaviour of the water network has been found to be unimportant for polar solvation.

Relating composition, structural order, entropy and transport in multi-component molten salts.

Molecular dynamics simulations of the LiF-BeF(2) molten salt mixture are used to establish relationships between composition, structural order, entropy, and transport properties of multi-component ionic liquids. A sharp rise in tetrahedral order associated with formation of the fluoroberyllate network occurs for compositions with BeF(2) concentrations greater than that of the Li(2)BeF(4)-BeF(2) eutectic. The excess entropy of the liquid in this regime, within the pair correlation approximation, is strongly correlated with the local tetrahedral order. The different degree of participation of beryllium, fluorine, and lithium ions in the cooperative dynamics of the fluoroberyllate network can be related to the degree of deviation from Rosenfeld-type excess entropy scaling, with the lithium ions remaining essentially unaffected by the liquid state network. We demonstrate that the deviations from Nernst-Einstein and Stokes-Einstein behaviour emerge only in temperature-composition regimes where tetrahedral order strongly correlates with the pair entropy. Implications for understanding structure-property relationships in other ionic liquids, such as molten salts, oxide melts, and RTILs are considered.

Tetrahedral order, pair correlation entropy, and waterlike liquid state anomalies: comparison of GeO2 with BeF2, SiO2, and H2O.

Molecular dynamics simulations of the Oeffner-Elliot model of germania (GeO(2)) are performed to identify nested regions of anomalous behavior in structural order, diffusivity, and pair entropy in the density-temperature plane, analogous to that seen in BeF(2), SiO(2), and H(2)O. The decreasing constraint of local tetrahedrality in GeO(2), compared to SiO(2) and BeF(2), substantially lowers the onset temperatures for anomalous behavior relative to the experimental melting temperatures (T(m)). Germania resembles water, more strongly than the ionic melts, in terms of temperatures for onset of anomalous behavior as well as in the order maps; for example, the structural anomaly sets in at 3.42T(m) in BeF(2), 3.09T(m) in SiO(2), 1.43T(m) in GeO(2), and 1.21T(m) in H(2)O. The detailed shapes of the anomalous regimes vary for different systems but the relative temperatures of onset for different anomalies are very similar in the different systems. The pair correlation entropy is shown to be a crucial and experimentally accessible quantity for relating structure, entropy, and diffusivity that could be potentially useful for a large class of inorganic ionic liquids.