Simulation and reference interaction site model theory of methanol and carbon tetrachloride mixtures

Competition between the hydrogen bond and the halogen bond in a [CH3OH-CCl4] complex: a matrix isolation IR spectroscopy and computational study

Methanol (CH3OH) is the simplest alcohol and carbon tetrachloride (CCl4) is widely used as a solvent in the chemical industry. CH3OH and CCl4 are both important volatile substances in the atmosphere and CCl4 is an important precursor for atmospheric ozone depletion. Moreover, mixtures of CH3OH and CCl4 are an important class of non-aqueous mixtures as they exhibit a large deviation from Raoult's law. The specific interaction between CH3OH and CCl4 is not yet investigated experimentally. The interaction between CH3OH and CCl4 at the molecular level can be twofold: hydrogen bond (O-HCl) and halogen bond (C-ClO) interaction. One halogen bonded minimum and two hydrogen bonded minima are identified in the dimer potential energy surface. Herein, the 1 : 1 complex of [CH3OH-CCl4] has been characterised using matrix-isolation infrared spectroscopy and electronic structure calculations to investigate the competition between hydrogen bonded and halogen bonded complexes. Vibrational spectra have been monitored in the C-Cl, C-O, and O-H stretching regions. The exclusive formation of halogen bonded 1 : 1 complexes in argon and nitrogen matrices is confirmed by a combination of experimental and simulated vibrational frequency, stabilisation energy, energy decomposition analysis, and natural bond orbital and atoms-in-molecules analyses. This investigation helps to understand the specific interactions in the [CH3OH-CCl4] mixture and also the possibilities of formation of halogen bonded atmospheric complexes that may influence the atmospheric chemical activities, and enhance aerosol formation and deposition of CCl4.

The effect of surface roughness on the phase behavior of colloidal particles

Shape anisotropy of colloidal particles can give rise to complex intermolecular interactions that determine particle packing and phase behavior. The vapor-liquid coexistence curves of attractive rough particles display a shift when compared to attractive smooth spherical particles. We use Integral Equation Theory (IET) to determine the vapor-liquid spinodal phase diagram of smooth and rough colloidal particles interacting through square-well attraction. Additionally, we use Gibbs Ensemble Monte Carlo (GEMC) simulations to locate their vapor-liquid coexistence curves. We model a rough colloidal particle as a spherical core with small beads embedded on its surface. The critical point of smooth spherical particle systems predicted by theory and simulations is in quantitative agreement. An increase in surface roughness due to an increase in either the number of beads or the diameter of the beads has a modest effect on the local structure of the system in the supercritical region. In contrast, increasing surface roughness consistently shifts the vapor-liquid coexistence curves to higher temperatures. The critical temperature is found to be a quadratic function of the number of beads. At a fixed bead size and number of beads, the critical temperature does not vary with the arrangement of beads on the core. Both IET and GEMC simulations predict that unlike critical temperatures, critical packing fractions vary non-monotonically with surface roughness. We find that the feasibility and accuracy of the integral equation theory depend sensitively on the chosen closure combination.

Evidence of Structural Inhomogeneities in Hard-Soft Dimeric Particles without Attractive Interactions

We perform Monte Carlo simulations of a simple hard-soft dimeric model constituted by two tangent spheres experiencing different interactions. Specifically, two hard spheres belonging to different dimers interact via a bare hard-core repulsion, whereas two soft spheres experience a softly repulsive Hertzian interaction. The cross correlations are soft as well. By exploring a wide range of temperatures and densities we investigate the capability of this model to document the existence of structural inhomogeneities indicating the possible onset of aggregates, even if no attraction is set. The fluid phase behavior is studied by analyzing structural and thermodynamical properties of the observed structures, in particular by computing radial distribution functions, structure factors and cluster size distributions. The numerical results are supported by integral equation theories of molecular liquids which allow for a finer and faster spanning of the temperature-density diagram. Our results may serve as a framework for a more systematic investigation of self-assembled structures of functionalized hard-soft dimers able to aggregate in a variety of structures widely oberved in colloidal dispersion.

Modeling micro-heterogeneity in mixtures: The role of many body correlations

A two-component interaction model is introduced herein, which allows us to describe macroscopic miscibility with various modes of tunable micro-segregation, ranging from phase separation to micro-segregation, and is in excellent agreement with structural quantities obtained from simulations and the liquid state hypernetted-chain like integral equation theory. The model is based on the conjecture that the many-body correlation bridge function term in the closure relation can be divided into one part representing the segregation effects, which are modeled herein, and the usual part representing random many body fluctuations. Furthermore, the model allows us to fully neglect these second contributions, thus increasing the agreement between the simulations and the theory. The analysis of the retained part of the many body correlations gives important clues about how to model the many body bridge functions for more realistic systems exhibiting micro-segregation, such as aqueous mixtures.

Development of molecular closures for the reference interaction site model theory with application to square-well and Lennard-Jones homonuclear diatomics

Inspired by significant improvements obtained for the performances of the polymer reference interaction site model (PRISM) theory of the fluid phase when coupled with 'molecular closures' (Schweizer and Yethiraj 1993 J. Chem. Phys. 98 9053), we exploit a matrix generalization of this concept, suitable for the more general RISM framework. We report a preliminary test of the formalism, as applied to prototype square-well homonuclear diatomics. As for the structure, comparison with Monte Carlo shows that molecular closures are slightly more predictive than their 'atomic' counterparts, and thermodynamic properties are equally accurate. We also devise an application of molecular closures to models interacting via continuous, soft-core potentials, by using well established prescriptions in liquid state perturbation theories. In the case of Lennard-Jones dimers, our scheme definitely improves over the atomic one, providing semi-quantitative structural results, and quite good estimates of internal energy, pressure and phase coexistence. Our finding paves the way to a systematic employment of molecular closures within the RISM framework to be applied to more complex systems, such as molecules constituted by several non-equivalent interaction sites.

Determination of size of molecular clusters of ethanol by means of NMR diffusometry and hydrodynamic calculations

The microscopic structure of ethanol in the liquid state is characterized as a dynamic equilibrium of hydrogen-bonded clusters of different sizes and topologies. We have developed a novel method for determination of the average size of the clusters that combines the measurement of diffusion coefficient by means of NMR diffusometry technique and hydrodynamic simulations. The approach includes the use of HydroNMR [J. Garcı̀a de la Torre, M. L. Huertas, and B. Carrasco, J. Magn. Reson. 147, 2000, 138] for small molecules, which is attained here by the calibration procedure using a dilute solution of tetramethylsilane. It is thus possible to correlate the experimentally determined diffusion coefficient of ethanol with calculated diffusion coefficients of the modeled clusters of different sizes. We found that average size of the clusters in 0.16 M solution of ethanol in n-hexane corresponds to the monomer above 300 K and to the pentamer/hexamer below 240 K. The clusters in the case of 0.44 M solution are generally slightly larger, from the average size corresponding to the dimer at 320 K and the hexamer at 210 K.

Association effects in the {methanol + inert solvent} system via Monte Carlo simulations. I. Structure

In this work, the clusters residing in the {methanol + inert solvent} binary system have been characterized using a specific methodology in the framework of Monte Carlo molecular simulations. The cluster classification scheme considered distinguishes into five types: linear chains, cyclic clusters or isolated rings, branched linear chains, branched cyclic clusters, and composite rings. The procedure allows one to compute the next rich structural information: the fraction of molecules in the monomer or associated state, the fraction of each type of aggregate with a given size (and of molecules belonging to them), and the most probable and average cluster size for each type; likewise, the degree of branching in branched linear chains and the size distribution of the inner ring in branched cyclic clusters can be quantified. Specifically, all these properties were obtained for the {Optimized Potential for Liquid Simulation methanol + Lennard-Jones spheres} system at 298.15 K and 1 bar throughout the composition range. The results have provided a complete structural picture of this mixture describing comprehensively the effect of dilution into the hydrogen-bonded network of the pure associated fluid.