Thiago Colla

Thermodynamic and dynamic anomalies for a three-dimensional isotropic core-softened potential.

Using molecular-dynamics simulations and integral equations (Rogers-Young, Percus-Yevick, and hypernetted chain closures) we investigate the thermodynamics of particles interacting with continuous core-softened intermolecular potential. Dynamic properties are also analyzed by the simulations. We show that, for a chosen shape of the potential, the density, at constant pressure, has a maximum for a certain temperature. The line of temperatures of maximum density (TMD) was determined in the pressure-temperature phase diagram. Similarly the diffusion constant at a constant temperature, D, has a maximum at a density rho(max) and a minimum at a density rho(min) < rho(max). In the pressure-temperature phase diagram the line of extrema in diffusivity is outside of the TMD line. Although this interparticle potential lacks directionality, this is the same behavior observed in simple point charge/extended water.

Structural anomalies for a three dimensional isotropic core-softened potential

Using molecular dynamics simulations we investigate the structure of a system of particles interacting through a continuous core-softened interparticle potential. We found for the translational order parameter t a local maximum at a density rho(t-max) and a local minimum at rho(t-min)>rho(t-max). Between rho(t-max) and rho(t-min), the t parameter anomalously decreases upon increasing pressure. For the orientational order parameter Q(6) a maximum was observed at a density rho(t-max)<rho(Qmax)<rho(t-min). For densities between rho(Qmax) and rho(t-min), both the translational (t) and orientational (Q(6)) order parameters have anomalous behavior. We know that this system also exhibits density and diffusion anomalies. We found that the region in the pressure-temperature phase diagram of the structural anomaly englobes the region of the diffusion anomaly that is larger than the region limited by the temperature of maximum density. This cascade of anomalies (structural, dynamic, and thermodynamic) for our model has the same hierarchy as that observed for the simple point charge/extended water.

Equilibrium properties of charged microgels: A Poisson-Boltzmann-Flory approach

The equilibrium properties of ionic microgels are investigated using a combination of the Poisson-Boltzmann and Flory theories. Swelling behavior, density profiles, and effective charges are all calculated in a self-consistent way. Special attention is given to the effects of salinity on these quantities. In accordance with the traditional ideal Donnan equilibrium theory, it is found that the equilibrium microgel size is strongly influenced by the amount of added salt. Increasing the salt concentration leads to a considerable reduction of the microgel volume, which therefore releases its internal material - solvent molecules and dissociated ions - into the solution. Finally, the question of charge renormalization of ionic microgels in the context of the cell model is briefly addressed.

Equation of state of charged colloidal suspensions and its dependence on the thermodynamic route

The thermodynamic properties of highly charged colloidal suspensions in contact with a salt reservoir are investigated in the framework of the renormalized Jellium model (RJM). It is found that the equation of state is very sensitive to the particular thermodynamic route used to obtain it. Specifically, the osmotic pressure calculated within the RJM using the contact value theorem can be very different from the pressure calculated using the Kirkwood-Buff fluctuation relations. On the other hand, Monte Carlo simulations show that both the effective pair potentials and the correlation functions are accurately predicted by the RJM. It is suggested that the lack of self-consistency in the thermodynamics of the RJM is a result of neglected electrostatic correlations between the counterions and coions.

The renormalized Jellium model of colloidal suspensions with multivalent counterions

An extension of the renormalized Jellium model which allows to study colloidal suspensions containing trivalent counterions is proposed. The theory is based on a modified Poisson-Boltzmann equation which incorporates the effects of counterion correlations near the colloidal surfaces using a new boundary condition. The renormalized charges, the counterion density profiles, and osmotic pressures can be easily calculated using the modified renormalized Jellium model. The results are compared with the ones obtained using the traditional Wigner-Seitz (WS) cell approximation also with a new boundary condition. We find that while the thermodynamic functions obtained within the renormalized Jellium model are in a good agreement with their WS counterpart, the effective charges predicted by the two theories can be significantly different.

Charge neutrality breakdown in confined aqueous electrolytes: Theory and simulation

We study, using Density Functional theory (DFT) and Monte Carlo simulations, aqueous electrolyte solutions between charged infinite planar surfaces, in contact with a bulk salt reservoir. In agreement with recent experimental observations [Z. Luo et al., Nat. Commun. 6, 6358 (2015)], we find that the confined electrolyte lacks local charge neutrality. We show that a DFT based on a bulk-HNC expansion properly accounts for strong electrostatic correlations and allows us to accurately calculate the ionic density profiles between the charged surfaces, even for electrolytes containing trivalent counterions. The DFT allows us to explore the degree of local charge neutrality violation, as a function of plate separation and bulk electrolyte concentration, and to accurately calculate the interaction force between the charged surfaces.

Effective interactions in polydisperse systems of penetrable macroions

A coarse-graining approach based on the linear response theory is applied to deduce general expressions for the effective pair potentials in a multi-component system of soft macroions. Within the underlying approximations, the effective pair potentials can be written as simple convolutions between the intrinsic macroion charge distributions and a Yukawa-like potential which effectively contains the averaged contributions from the small ions. Two different charge distributions are assigned to the soft macroions: a Gaussian-like diffuse distribution and a uniform charge distribution inside the particle cores. The resulting effective pair potentials are then used in an effective model based on the hyppernetted chain approximation to investigate the structural properties of a two-component system of oppositely charged particles as a function of the various system parameters. It is found that the condensation of counterions is much stronger in the case of particles with a Gaussian charge distribution, leading to much weaker electrostatic interactions and less structured pair correlations in comparison with the system of uniformly charged macroions.

Lattice Model of an Ionic Liquid at an Electrified Interface

We study ionic liquids interacting with electrified interfaces. The ionic fluid is modeled as a Coulomb lattice gas. We compare the ionic density profiles calculated using a popular modified Poisson-Boltzmann equation with the explicit Monte Carlo simulations. The modified Poisson-Boltzmann theory fails to capture the structural features of the double layer and is also unable to correctly predict the ionic density at the electrified interface. The lattice Monte Carlo simulations qualitatively capture the coarse-grained structure of the double layer in the continuum. We propose a convolution relation that semiquantitatively relates the ionic density profiles of a continuum ionic liquid and its lattice counterpart near an electrified interface.

Self-Assembly of Ionic Microgels Driven by an Alternating Electric Field: Theory, Simulations, and Experiments

The structural properties of a system of ionic microgels under the influence of an alternating electric field are investigated both theoretically and experimentally. This combined investigation aims to shed light on the structural transitions that can be induced by changing either the driving frequency or the strength of the applied field, which range from string-like formation along the field to crystal-like structures across the orthogonal plane. In order to highlight the physical mechanisms responsible for the observed particle self-assembly, we develop a coarse-grained description, in which effective interactions among the charged microgels are induced by both equilibrium ionic distributions and their time-averaged hydrodynamic responses to the applied field. These contributions are modeled by the buildup of an effective dipole moment at the microgels backbones, which is partially screened by their ionic double layer. We show that this description is able to capture the structural properties of this system, allowing for very good agreement with the experimental results. The model coarse-graining parameters are indirectly obtained via the measured pair distribution functions and then further assigned with a clear physical interpretation, allowing us to highlight the main physical mechanisms accounting for the observed self-assembly behavior.

Yukawa particles in a confining potential

We study the density distribution of repulsive Yukawa particles confined by an external potential. In the weak coupling limit, we show that the mean-field theory is able to accurately account for the particle distribution. In the strong coupling limit, the correlations between the particles become important and the mean-field theory fails. For strongly correlated systems, we construct a density functional theory which provides an excellent description of the particle distribution, without any adjustable parameters.