RESUMEN

A dramatic improvement is reported in the stability of colloidal particles when stabilizing surface grafts are systematically shortened from small polymers to single monomers. The colloidal dispersions consist of fluorinated latex particles, exhibiting a weak van der Waals attraction, with grafted steric layers of poly(ethylene glycol) (PEG) of different chain lengths. Using an effective salting-out electrolyte, Na2CO3, particle aggregates are detected above a threshold salt concentration that is independent of the particle concentration. The results are interpreted in terms of a sudden onset of nondispersibility of single particles, triggered by the solvent not completely wetting particle surfaces. By decreasing the PEG chain length, the threshold salt concentration is found to increase sharply. For grafts with just a single ethylene glycol group, dispersions remain stable up to exceedingly high concentrations of Na2CO3. However, on removal of the surface coverage altogether, the classical stability behavior of charge-stabilized dispersions is recovered. The behavior can be captured by a simple model that incorporates effective polymer-solvent interactions in the presence of an electrolyte.

RESUMEN

This study explores the thermodynamic and vibrational properties of water in the three-dimensional environment of solvated ions and small molecules using molecular simulations. The spectrum of intermolecular vibrations in liquid solvents provides detailed information on the shape of the local potential energy surface, which in turn determines local thermodynamic properties such as the entropy. Here, we extract this information using a spatially resolved extension of the two-phase thermodynamics method to estimate hydration water entropies based on the local vibrational density of states (3D-2PT). Combined with an analysis of solute-water and water-water interaction energies, this allows us to resolve local contributions to the solvation enthalpy, entropy, and free energy. We use this approach to study effects of ions on their surrounding water hydrogen bond network, its spectrum of intermolecular vibrations, and resulting thermodynamic properties. In the three-dimensional environment of polar and nonpolar functional groups of molecular solutes, we identify distinct hydration water species and classify them by their characteristic vibrational density of states and molecular entropies. In each case, we are able to assign variations in local hydration water entropies to specific changes in the spectrum of intermolecular vibrations. This provides an important link for the thermodynamic interpretation of vibrational spectra that are accessible to far-infrared absorption and Raman spectroscopy experiments. Our analysis provides unique microscopic details regarding the hydration of hydrophobic and hydrophilic functional groups, which enable us to identify interactions and molecular degrees of freedom that determine relevant contributions to the solvation entropy and consequently the free energy.

RESUMEN

Based on the dressed-ion theory and a simple physical argument regarding the conductivity of the solution, we derive a relation between the ionic strength and dielectric constant of an electrolyte solution. At its simplest, this model gives the dielectric constant at low ionic strength I as εr(I) = εr(0)(1 + αI)-1, where α (the excess polarization) is directly related to the dressed-ion charge. One contribution to the origin of the dielectric decrement is thus seen to stem from the electrostatic screening of the ions in solution, with no solvent contributions necessary.

RESUMEN

Fluctuating hydrodynamics (FHD) is a general framework of mesoscopic modeling and simulation based on conservational laws and constitutive equations of linear and nonlinear responses. However, explicit representation of electrical forces in FHD has yet to appear. In this work, we devised an Ansatz for the dynamics of dipole moment densities that is linked with the Poisson equation of the electrical potential ϕ in coupling to the other equations of FHD. The resulting ϕ-FHD equations then serve as a platform for integrating the essential forces, including electrostatics in addition to hydrodynamics, pressure-volume equation of state, surface tension, and solvent-particle interactions that govern the emergent behaviors of molecular systems at an intermediate scale. This unique merit of ϕ-FHD is illustrated by showing that the water dielectric function and ion hydration free energies in homogeneous and heterogenous systems can be captured accurately via the mesoscopic simulation. Furthermore, we show that the field variables of ϕ-FHD can be mapped from the trajectory of an all-atom molecular dynamics simulation such that model development and parametrization can be based on the information obtained at a finer-grained scale. With the aforementioned multiscale capabilities and a spatial resolution as high as 5 Å, the ϕ-FHD equations represent a useful semi-explicit solvent model for the modeling and simulation of complex systems, such as biomolecular machines and nanofluidics.

Asunto(s)

RESUMEN

We introduce a simple improvement on the method to calculate equilibrium entropy differences between classical energy levels proposed by Davis [S. Davis, Phys. Rev. E 84, 050101 (2011)]. We demonstrate that the modification is superior to the original whenever the energy levels are sufficiently closely spaced or whenever the microcanonical averaging needed in the method is carried out by importance sampling Monte Carlo. We also point out the necessary adjustments if Davis's method (improved or not) is to be used with molecular dynamics simulations.

Asunto(s)

RESUMEN

Monte Carlo switching moves ("perturbations") are defined between two or more classical Hamiltonians sharing a common ground-state energy. The ratio of the density of states (DOS) of one system to that of another is related to the ensemble averages of the microcanonical acceptance probabilities of switching between these Hamiltonians, analogously to the case of Bennett's acceptance ratio method for the canonical ensemble [C. H. Bennett, J. Comput. Phys. 22, 245 (1976)]. Thus, if the DOS of one of the systems is known, one obtains those of the others and, hence, the partition functions. As a simple test case, the vapor pressure of an anharmonic Einstein crystal is computed, using the harmonic Einstein crystal as the reference system in one dimension; an auxiliary calculation is also performed in three dimensions. As a further example of the algorithm, the energy dependence of the ratio of the DOS of the square-well and hard-sphere tetradecamers is determined, from which the temperature dependence of the constant-volume heat capacity of the square-well system is calculated and compared with canonical Metropolis Monte Carlo estimates. For these cases and reference systems, the perturbation calculations exhibit a higher degree of convergence per Monte Carlo cycle than Wang-Landau (WL) sampling, although for the one-dimensional oscillator the WL sampling is ultimately more efficient for long runs. Last, we calculate the vapor pressure of liquid gold using an empirical Sutton-Chen many-body potential and the ideal gas as the reference state. Although this proves the general applicability of the method, by its inherent perturbation approach the algorithm is suitable for those particular cases where the properties of a related system are well known.

RESUMEN

With computational efficacy in mind, a one-center model for linear molecules is heuristically sketched. When parametrized for CO(2), all parameters save for two are supplied from literature quantum chemistry calculation or, in one case, heuristic argument. Using the remaining two adjustable parameters, the mean unsigned relative errors (predicted/observed) over the temperature range 220-290 K are 4.7% for the energy of vaporization, 0.6% for the liquid, and 8.0% for the vapor coexistence densities, respectively. The critical temperature is estimated at T(c) = 308 K, the critical density at ρ(c) = 0.460 g/cm(3), and the critical pressure at p(c) = 8.26 ± 0.11 MPa. This order of accuracy is comparable to that of many all-atom potential descriptions of CO(2) but is obtained at roughly nine times the speed. When supplied with the experimental bond length, somewhat worse agreement with experiment is exhibited for the neutron-weighted atomic pair distribution function of the liquid. This disparity is tentatively attributed to an overestimated electrostatic quadrupole-quadrupole interaction relative to the other forces present.

RESUMEN

A number of simple pair interaction potentials of the carbon dioxide molecule are investigated and found to underestimate the magnitude of the second virial coefficient in the temperature interval 220-448 K by up to 20%. Also the third virial coefficient is underestimated by these models. A rigid, polarizable, three-site interaction potential reproduces the experimental second and third virial coefficients to within a few percent. It is based on the modified Buckingham exp-6 potential, an anisotropic Axilrod-Teller correction, and Gaussian charge densities on the atomic sites with an inducible dipole at the center of mass. The electric quadrupole moment, polarizability, and bond distances are set to equal experiment. Density of the fluid at 200 and 800 bars pressure is reproduced to within some percent of observation over the temperature range 250-310 K. The dimer structure is in passable agreement with electronically resolved quantum-mechanical calculations in the literature, as are those of the monohydrated monomer and dimer complexes using the Gaussian charge polarizable model water potential. Qualitative agreement with experiment is also obtained, when quantum corrections are included, for the relative stability of the trimer conformations, which is not the case for the pair potentials.

Asunto(s)