Germán Pérez-Sánchez

Modeling self-assembly of silica/surfactant mesostructures in the templated synthesis of nanoporous solids.

A novel coarse-grained (CG) model to study the self-assembly of silica/surfactant mesostructures during the synthesis of periodic mesoporous silica is reported. Molecular dynamics simulations of hexadecyltrimethylammonium bromide (also called cetyltrimethylammonium bromide, or CTAB) surfactants in water and in aqueous silicate solutions have been performed to understand micelle formation, micelle growth, and their size evolution during the synthesis of surfactant-templated mesoporous materials. Direct comparison of density profiles obtained for preassembled micelles employing an all-atom description, AA, with those calculated with the CG model has been carried out for checking the validity of the latter model. Good agreement between AA and CG approaches was found, demonstrating the potential of the CG approximation for modeling these highly complex systems. The micelle formation and micelle fusion/fission processes were analyzed after performing long CG simulations for surfactant and ionized silica-surfactant aqueous solutions. We observed the formation of rodlike micelles in the case of silica-surfactant solutions, while spherical micelles were stable under the same conditions for the CTAB+H(2)O system. This demonstrates that the interaction of anionic silicates with cationic surfactants promotes a sphere-to-rod transition in surfactant solutions, a key step in the synthesis of nanoporous silica materials.

Fluid-solid equilibrium of carbon dioxide as obtained from computer simulations of several popular potential models: the role of the quadrupole.

In this work the solid-fluid equilibrium for carbon dioxide (CO2) has been evaluated using Monte Carlo simulations. In particular the melting curve of the solid phase denoted as I, or dry ice, was computed for pressures up to 1000 MPa. Four different models, widely used in computer simulations of CO2 were considered in the calculations. All of them are rigid non-polarizable models consisting of three Lennard-Jones interaction sites located on the positions of the atoms of the molecule, plus three partial charges. It will be shown that although these models predict similar vapor-liquid equilibria their predictions for the fluid-solid equilibria are quite different. Thus the prediction of the entire phase diagram is a severe test for any potential model. It has been found that the Transferable Potentials for Phase Equilibria (TraPPE) model yields the best description of the triple point properties and melting curve of carbon dioxide. It is shown that the ability of a certain model to predict the melting curve of carbon dioxide is related to the value of the quadrupole moment of the model. Models with low quadrupole moment tend to yield melting temperatures too low, whereas the model with the highest quadrupole moment yields the best predictions. That reinforces the idea that not only is the quadrupole needed to provide a reasonable description of the properties in the fluid phase, but also it is absolutely necessary to describe the properties of the solid phase.

Thermal properties of ionic systems near the liquid-liquid critical point

Isobaric heat capacity per unit volume, C(p), and excess molar enthalpy, h(E), were determined in the vicinity of the critical point for a set of binary systems formed by an ionic liquid and a molecular solvent. Moreover, and, since critical composition had to be accurately determined, liquid-liquid equilibrium curves were also obtained using a calorimetric method. The systems were selected with a view on representing, near room temperature, examples from clearly solvophobic to clearly coulombic behavior, which traditionally was related with the electric permittivity of the solvent. The chosen molecular compounds are: ethanol, 1-butanol, 1-hexanol, 1,3-dichloropropane, and diethylcarbonate, whereas ionic liquids are formed by imidazolium-based cations and tetrafluoroborate or bis-(trifluromethylsulfonyl)amide anions. The results reveal that solvophobic critical behavior-systems with molecular solvents of high dielectric permittivity-is very similar to that found for molecular binary systems. However, coulombic systems-those with low permittivity molecular solvents-show strong deviations from the results usually found for these magnitudes near the liquid-liquid phase transition. They present an extremely small critical anomaly in C(p)-several orders of magnitude lower than those typically obtained for binary mixtures-and extremely low h(E)-for one system even negative, fact not observed, up to date, for any liquid-liquid transition in the nearness of an upper critical solution temperature.

Critical behavior of static properties for nitrobenzene-alkane mixtures.

We present experimental data of the isobaric heat capacity per unit volume C(p,x)V(-1) for mixtures containing nitrobenzene and an alkane (C(N)H(2N+2), with N ranging from 6 to 15) upon approaching their liquid-liquid critical points along a path of constant composition. Values for the critical amplitude A(+) have been determined. They have been combined with the previously reported ones for the leading term of the coexistence-curve width to obtain, with the aid of well-known universal relations, the critical amplitudes of the correlation length and of the osmotic susceptibility. The trends of all these critical parameters, which exhibit anomalous behavior in the low N region, are discussed in terms of particular microscopic phenomena characterizing NB-C(N)H(2N+2) mixtures. The work is completed with an analysis of the analog of the Yang-Yang anomaly in liquid-liquid criticality: the behavior of the partial molar heat capacities of the two liquid components is found to illustrate previously uncovered asymmetry effects.

Thermodynamic consistency near the liquid-liquid critical point.

The thermodynamic consistency of the isobaric heat capacity per unit volume at constant composition C(p,x) and the density rho near the liquid-liquid critical point is studied in detail. To this end, C(p,x)(T), rho(T), and the slope of the critical line (dT/dp)(c) for five binary mixtures composed by 1-nitropropane and an alkane were analyzed. Both C(p,x)(T) and rho(T) data were measured along various quasicritical isopleths with a view to evaluate the effect of the uncertainty in the critical composition value on the corresponding critical amplitudes. By adopting the traditionally employed strategies for data treatment, consistency within 0.01 K MPa(-1) (or 8%) is attained, thereby largely improving the majority of previous results. From temperature range shrinking fits and fits in which higher-order terms in the theoretical expressions for C(p,x)(T) and rho(T) are included, we conclude that discrepancies come mainly from inherent difficulties in determining the critical anomaly of rho accurately: specifically, to get full consistency, higher-order terms in rho(T) are needed; however, the various contributions at play cannot be separated unambiguously. As a consequence, the use of C(p,x)(T) and (dT/dp)(c) for predicting the behavior of rho(T) at near criticality appears to be the best choice at the actual experimental resolution levels. Furthermore, the reasonably good thermodynamic consistency being encountered confirms that previous arguments appealing to the inadequacy of the theoretical expression relating C(p,x) and rho for describing data in the experimentally accessible region must be fairly rejected.

Heat capacity anomalies along the critical isotherm in fluid-fluid phase transitions

The behavior of the isochoric heat capacity of pure fluids and the isobaric heat capacity at constant composition of binary mixtures along isothermal paths of approach to liquid-gas and liquid-liquid critical points is studied. From the complete scaling formulation of fluid-fluid criticality, explicit expressions for the critical amplitudes of the leading /Y-Y(c)/(-alpha/beta) (where Y can be the density or the mole fraction) contributions are found to reveal previously discovered features of the scaling function, whereas the nature of the most important asymmetry-related terms is characterized. Data for pure toluene and for the binary mixture nitromethane-isobutanol are described within experimental uncertainty using the /Y-Y(c)/(-alpha/beta) singularity plus a linear term. Extensive data for mixtures allow proper visualization of the topological features of the heat capacity-density-temperature surface in the critical region.

CHAPTER 11:Critical Behaviour: Pure Fluids and Mixtures

Experimental findings on critical behaviour for fluids and fluid mixtures over the past few years have been summarized. This includes the analysis of the complete scaling formulation of asymmetric criticality for the coexistence curves of liquid-gas phase transitions and also for liquid-liquid phase separation in binary mixtures. For these phase transitions, we also reviewed results for ionic solutions as well as the historically elusive check of the thermodynamic consistency between the density and the heat capacity.