Simulations of dipolar fluids using effective many-body isotropic interactions

Arrested dynamics of the dipolar hard sphere model

We report the combined results of molecular dynamics simulations and theoretical calculations concerning various dynamical arrest transitions in a model system representing a dipolar fluid, namely, N (soft core) rigid spheres interacting through a truncated dipole-dipole potential. By exploring different regimes of concentration and temperature, we find three distinct scenarios for the slowing down of the dynamics of the translational and orientational degrees of freedom: at low (η = 0.2) and intermediate (η = 0.4) volume fractions, both dynamics are strongly coupled and become simultaneously arrested upon cooling. At high concentrations (η≥ 0.6), the translational dynamics shows the features of an ordinary glass transition, either by compressing or cooling down the system, but with the orientations remaining ergodic, thus indicating the existence of partially arrested states. In this density regime, but at lower temperatures, the relaxation of the orientational dynamics also freezes. The physical scenario provided by the simulations is discussed and compared against results obtained with the self-consistent generalized Langevin equation theory, and both provide a consistent description of the dynamical arrest transitions in the system. Our results are summarized in an arrested states diagram which qualitatively organizes the simulation data and provides a generic picture of the glass transitions of a dipolar fluid.

A Solvation/Desolvation Indicator Based on van der Waals Interactions between Solvents and Porphyrins

Solvation is a ubiquitous phenomenon associated with molecules in solutions. It often determines the equilibria of molecular systems and the rates of chemical reactions. Van der Waals interactions (a general term) includes weak interactions among noncharged compounds and it contributes significantly to solvation. The distinct observation of van der Waals interaction between solvent and porphyrin derivatives is reported herein. Bis(imidazolylporphyrinatozinc) structures connected through a 1,3-butadiyne moiety give two types of coordination polymers, E (extended) and S (stacked) polymers, exclusively. E polymers have larger solvent-accessible surface areas than the corresponding S polymers. Therefore, E polymers are better solvated than S polymers, providing an indicator of solvation and desolvation for the solvents used. A simple method (like a litmus test) was developed to evaluate the solvation ability of various solvents. Sixty-seven solvents and liquid compounds were tested, under the same conditions, using a conventional UV/Vis spectrometer. The results revealed a new liquid group with high solvation ability towards the porphyrins, and clarified van der Waals interaction assisted by secondary interaction on the substituents. The indicator system should contribute to the solution chemistry of molecules and materials, and to supramolecular chemistry interactions among hetero components.

Separating the effects of repulsive and attractive forces on the phase diagram, interfacial, and critical properties of simple fluids

Molecular simulations in the canonical and isothermal-isobaric ensembles were performed to study the effect of varying the shape of the intermolecular potential on the phase diagram, critical, and interfacial properties of model fluids. The molecular interactions were modeled by the Approximate Non-Conformal (ANC) theory potentials. Unlike the Lennard-Jones or Morse potentials, the ANC interactions incorporate parameters (called softnesses) that modulate the steepness of the potential in their repulsive and attractive parts independently. This feature allowed us to separate unambiguously the role of each region of the potential on setting the thermophysical properties. In particular, we found positive linear correlation between all critical coordinates and the attractive and repulsive softness, except for the critical density and the attractive softness which are negatively correlated. Moreover, we found that the physical properties related to phase coexistence (such as span of the liquid phase between the critical and triple points, variations in the P-T vaporization curve, interface width, and surface tension) are more sensitive to changes in the attractive softness than to the repulsive one. Understanding the different roles of attractive and repulsive forces on phase coexistence may contribute to developing more accurate models of liquids and their mixtures.