On the entropy of mixing of a binary mixture of hard spheres

Thermodynamically Stable Colloidal Solids: Interfacial Thermodynamics from the Particle Size Distribution

True thermodynamic stability of a solid colloidal dispersion is generally unexpected, so much that thorough experimental validation of proposed stable systems remains incomplete. Such dispersions are under investigated and would be of interest due to their long-term stability and insensitivity to preparation pathway. We apply classical nucleation theory (CNT) to such colloidal systems, providing a relationship which links the size-dependent interfacial free energy density of the particles to their size distribution, and use this expression in the fitting of previously reported size distributions for putatively thermodynamically stable nanoparticles. Experimental data from a gold-thiol system exhibiting inverse coarsening or "digestive ripening" can be well-described in terms of a power-law dependence of the interfacial free energy γ on radius based on capacitive charging of the nanoparticles, going asr - 3 , as suggested by prior authors. Data from magnetite nanoparticles in highly basic solutions also can be well-fit using the CNT relation, but with γ going asr - 2 . Slightly better fits are possible if the power of the radius is non-integral, but we stress that more complex models of γ will require richer data sets to avoid the problem of overfitting. Some parameters of the fits are still robustly at odds with earlier models that implicitly assumed absolute thermodynamic stability: first, the extrapolated free energy density of the flat surface in these systems is small and positive, rather than strongly negative; second, the shape of the distributions indicates the solution phase to be supersaturated in monomer relative to the bulk, and thus that these two systems may only be metastable. For future work, we derive expressions for the important statistical thermodynamic and chemical parameters of the interface energy in terms of 1) the surfactant concentration, 2) the temperature dependence, and 3) the concentrations of particles in the tail of the distribution.

Finite-Size Charged Species Diffusion and pH Change in Nanochannels

Molecular transport through nanofluidic structures exhibits properties that are unique at the nanoscale. The high surface-to-volume ratio of nanometer-sized confined spaces renders particle interactions with the surface of central importance. The electrical double layer (EDL) at the solid-liquid interface of charged surfaces generates an enrichment of counterions and the exclusion of co-ions that lead to a change in their diffusivity. In addition, the diffusive transport is altered by steric and hydrodynamic interactions between fluid molecules and the boundaries. An extensive body of literature investigates molecular transport at the nanoscale. However, most studies account for ionic species as point charges, severely limiting the applicability of the results to "large" nanofluidic systems. Moreover, and even more importantly, the change of pH in the nanoconfined region inside nanochannels has been completely overlooked. Corroborated by experimental data, here we present an all-encompassing analysis of molecular diffusion from the micro- to the ultra-nanoscale. While accounting for finite-size ions, we compute self-consistently the pH inside the channels. Surprisingly, we found that the concentration of ions H⁺ can change by more than 2 orders of magnitude compared to the bulk, hugely affecting molecular transport. Further, we found that counterions exhibit both enrichment and exclusion, depending on the size of nanochannels. Achieving a greater understanding of the effective transport properties of fluids at the nanoscale will fill the gap in knowledge that still limits development of innovative systems for medicine and industrial applications alike.

Accurate schemes for calculation of thermodynamic properties of liquid mixtures from molecular dynamics simulations

We explore different schemes for improved accuracy of entropy calculations in aqueous liquid mixtures from molecular dynamics (MD) simulations. We build upon the two-phase thermodynamic (2PT) model of Lin et al. [J. Chem. Phys. 119, 11792 (2003)] and explore new ways to obtain the partition between the gas-like and solid-like parts of the density of states, as well as the effect of the chosen ideal "combinatorial" entropy of mixing, both of which have a large impact on the results. We also propose a first-order correction to the issue of kinetic energy transfer between degrees of freedom (DoF). This problem arises when the effective temperatures of translational, rotational, and vibrational DoF are not equal, either due to poor equilibration or reduced system size/time sampling, which are typical problems for ab initio MD. The new scheme enables improved convergence of the results with respect to configurational sampling, by up to one order of magnitude, for short MD runs. To ensure a meaningful assessment, we perform MD simulations of liquid mixtures of water with several other molecules of varying sizes: methanol, acetonitrile, N, N-dimethylformamide, and n-butanol. Our analysis shows that results in excellent agreement with experiment can be obtained with little computational effort for some systems. However, the ability of the 2PT method to succeed in these calculations is strongly influenced by the choice of force field, the fluidicity (hard-sphere) formalism employed to obtain the solid/gas partition, and the assumed combinatorial entropy of mixing. We tested two popular force fields, GAFF and OPLS with SPC/E water. For the mixtures studied, the GAFF force field seems to perform as a slightly better "all-around" force field when compared to OPLS+SPC/E.

Entropy of mixing for Ag(x)In(1-x) and Ag(x)Sn(1-x) liquid binary alloys

The entropy of mixing for Ag(x)In(1-x) and Ag(x)Sn(1-x) liquid binary alloys has been systematically investigated by using the perturbation theory. The interionic interactions as one of the basic ingredients are described by a local model pseudopotential. Since the metals forming the concerned alloys are less simple in nature the effect of the sp-d hybridization is appropriately taken into account through the interionic interactions. Results of our calculations across the full range of Ag concentrations are found to be good in agreement with the available experimental data.

Variational approach to hard sphere segregation under gravity

It is demonstrated that the minimization of the free energy functional for hard spheres and hard disks yields the result that excited granular materials under gravity segregate not only in the widely known "Brazil nut" fashion, i.e., with the larger particles rising to the top, but also in reverse "Brazil nut" fashion. Specifically, the local density approximation is used to investigate the crossover between the two types of segregation occurring in the liquid state, and the results are found to agree qualitatively with previously published results of simulation and of a simple model based on condensation.