N-Heterocyclic Carbenes as Key Intermediates in the Synthesis of Fused, Mesoionic, Tricyclic Heterocycles.

Coupling between 5-bromoimidazo[1,5-a]pyridinium salts and malonate or arylacetate esters leads to a facile and straightforward access to the new mesoionic, fused, tricyclic system of imidazo[2,1,5-cd]indolizinium-3-olate. Mechanistic studies show that the reaction pathway consists of nucleophilic aromatic substitution on the cationic, bicyclic heterocycle by an enolate-type moiety and in the nucleophilic attack of a transient free N-heterocyclic carbene (NHC) species on the ester group; the relative order of these two steps depends on the nature of the starting ester. This work highlights the valuable implementation of free NHC species as key intermediates in synthetic chemistry, beyond their classical use as stabilizing ligands or organocatalysts.

Self-assembly of polymer-grafted nanoparticles in thin films.

We use large-scale molecular dynamics simulations with a coarse-grained model to investigate the self-assembly of solvent-free grafted nanoparticles into thin free-standing films. Two important findings are highlighted. First, for appropriately chosen values of system parameters the nanoparticles spontaneously assemble into monolayer thick films. Further, the nanoparticles self-assemble into a variety of morphologies ranging from dispersed particles, finite stripes, long strings, to percolating networks. The main driving force for these morphologies is the competition between strong short-range attractions of the particle cores and long-range entropic repulsions of the grafted chains. The grafted nanoparticle systems provide practical means to realize two-length-scale systems that have been previously seen, using a simple two-dimensional model [G. Malescio and G. Pellicane, Nat. Mater., 2003, 2, 97], to generate a variety of morphologies. However, there are only relatively narrow ranges of interaction strengths and chain lengths for which anisotropic self-assembly is possible.

Group contribution methodology based on the statistical associating fluid theory for heteronuclear molecules formed from Mie segments.

A generalization of the recent version of the statistical associating fluid theory for variable range Mie potentials [Lafitte et al., J. Chem. Phys. 139, 154504 (2013)] is formulated within the framework of a group contribution approach (SAFT-γ Mie). Molecules are represented as comprising distinct functional (chemical) groups based on a fused heteronuclear molecular model, where the interactions between segments are described with the Mie (generalized Lennard-Jonesium) potential of variable attractive and repulsive range. A key feature of the new theory is the accurate description of the monomeric group-group interactions by application of a high-temperature perturbation expansion up to third order. The capabilities of the SAFT-γ Mie approach are exemplified by studying the thermodynamic properties of two chemical families, the n-alkanes and the n-alkyl esters, by developing parameters for the methyl, methylene, and carboxylate functional groups (CH3, CH2, and COO). The approach is shown to describe accurately the fluid-phase behavior of the compounds considered with absolute average deviations of 1.20% and 0.42% for the vapor pressure and saturated liquid density, respectively, which represents a clear improvement over other existing SAFT-based group contribution approaches. The use of Mie potentials to describe the group-group interaction is shown to allow accurate simultaneous descriptions of the fluid-phase behavior and second-order thermodynamic derivative properties of the pure fluids based on a single set of group parameters. Furthermore, the application of the perturbation expansion to third order for the description of the reference monomeric fluid improves the predictions of the theory for the fluid-phase behavior of pure components in the near-critical region. The predictive capabilities of the approach stem from its formulation within a group-contribution formalism: predictions of the fluid-phase behavior and thermodynamic derivative properties of compounds not included in the development of group parameters are demonstrated. The performance of the theory is also critically assessed with predictions of the fluid-phase behavior (vapor-liquid and liquid-liquid equilibria) and excess thermodynamic properties of a variety of binary mixtures, including polymer solutions, where very good agreement with the experimental data is seen, without the need for adjustable mixture parameters.

Accurate statistical associating fluid theory for chain molecules formed from Mie segments.

A highly accurate equation of state (EOS) for chain molecules formed from spherical segments interacting through Mie potentials (i.e., a generalized Lennard-Jones form with variable repulsive and attractive exponents) is presented. The quality of the theoretical description of the vapour-liquid equilibria (coexistence densities and vapour pressures) and the second-derivative thermophysical properties (heat capacities, isobaric thermal expansivities, and speed of sound) are critically assessed by comparison with molecular simulation and with experimental data of representative real substances. Our new EOS represents a notable improvement with respect to previous versions of the statistical associating fluid theory for variable range interactions (SAFT-VR) of the generic Mie form. The approach makes rigorous use of the Barker and Henderson high-temperature perturbation expansion up to third order in the free energy of the monomer Mie system. The radial distribution function of the reference monomer fluid, which is a prerequisite for the representation of the properties of the fluid of Mie chains within a Wertheim first-order thermodynamic perturbation theory (TPT1), is calculated from a second-order expansion. The resulting SAFT-VR Mie EOS can now be applied to molecular fluids characterized by a broad range of interactions spanning from soft to very repulsive and short-ranged Mie potentials. A good representation of the corresponding molecular-simulation data is achieved for model monomer and chain fluids. When applied to the particular case of the ubiquitous Lennard-Jones potential, our rigorous description of the thermodynamic properties is of equivalent quality to that obtained with the empirical EOSs for LJ monomer (EOS of Johnson et al.) and LJ chain (soft-SAFT) fluids. A key feature of our reformulated SAFT-VR approach is the greatly enhanced accuracy in the near-critical region for chain molecules. This attribute, combined with the accurate modeling of second-derivative properties, allows for a much improved global representation of the thermodynamic properties and fluid-phase equilibria of pure fluids and their mixtures.

SAFT-γ force field for the simulation of molecular fluids: 2. Coarse-grained models of greenhouse gases, refrigerants, and long alkanes.

In the first paper of this series [C. Avendaño, T. Lafitte, A. Galindo, C. S. Adjiman, G. Jackson, and E. A. Müller, J. Phys. Chem. B2011, 115, 11154] we introduced the SAFT-γ force field for molecular simulation of fluids. In our approach, a molecular-based equation of state (EoS) is used to obtain coarse-grained (CG) intermolecular potentials that can then be employed in molecular simulation over a wide range of thermodynamic conditions of the fluid. The macroscopic experimental data for the vapor-liquid equilibria (saturated liquid density and vapor pressure) of a given system are represented with the SAFT-VR Mie EoS and used to estimate effective intermolecular parameters that provide a good description of the thermodynamic properties by exploring a wide parameter space for models based on the Mie (generalized Lennard-Jones) potential. This methodology was first used to develop a simple single-segment CG Mie model of carbon dioxide (CO2) which allows for a reliable representation of the fluid-phase equilibria (for which the model was parametrized), as well as an accurate prediction of other properties such as the enthalpy of vaporization, interfacial tension, supercritical density, and second-derivative thermodynamic properties (thermal expansivity, isothermal compressibility, heat capacity, Joule-Thomson coefficient, and speed of sound). In our current paper, the methodology is further applied and extended to develop effective SAFT-γ CG Mie force fields for some important greenhouse gases including carbon tetrafluoride (CF4) and sulfur hexafluoride (SF6), modeled as simple spherical molecules, and for long linear alkanes including n-decane (n-C10H22) and n-eicosane (n-C20H42), modeled as homonuclear chains of spherical Mie segments. We also apply the SAFT-γ methodology to obtain a CG homonuclear two-segment Mie intermolecular potential for the more challenging polar and asymmetric compound 2,3,3,3-tetrafluoro-1-propene (HFO-1234yf), a novel replacement refrigerant with promising properties. The description of the fluid-phase behavior and the prediction of the other thermophysical properties obtained by molecular simulation using our SAFT-γ CG Mie force fields are found to be of comparable quality (and sometimes superior) to that obtained using the more sophisticated all-atom (AA) and united-atom (UA) models commonly employed in the field. We should emphasize that though the focus of our current work is on simple homonuclear models, the SAFT-γ methodology is based on a group contribution methodology which is naturally suited to the development of more sophisticated heteronuclear models.

Pressure and density scaling for colloid-polymer systems in the protein limit.

Grand canonical Monte Carlo and histogram reweighting techniques are used to study the fluid-phase behavior of an athermal system of colloids and nonadsorbing polymers on a fine lattice in the "protein limit," where polymer dimensions exceed those of the colloids. The main parameters are the chains' radius of gyration, R_{g}, the diameter of the colloids, σ_{c}, and the monomer diameter, σ_{s}. The phase behavior is controlled by the macroscopic size ratio, q_{r}=2R_{g}/σ_{c}, and the microscopic size ratio, d=σ_{s}/σ_{c}. The latter ratio is found to play a significant role in determining the critical monomer concentration for q_{r}≲4 and the critical colloid density for all chain lengths. However, the critical (osmotic) pressure is independent of the microscopic size ratio at all macroscopic size ratios studied. Quantitative agreement is observed between our simulation results and experimental data. We scale our results based on the polymer correlation length, which has previously been suggested to universally collapse these binodals [Bolhuis et al., Phys. Rev. Lett. 90, 068304 (2003); Fleer and Tuinier, Phys. Rev. E 76, 041802 (2007)]. While the density binodals exhibit universal characteristics along the low-colloid-density branch, such features are not present in the corresponding high-density phase. However, pressure binodals do collapse nicely under such a scaling, even far from the critical point, which allows us to produce a binodal curve whose shape is independent of either size ratio.

SAFT-γ force field for the simulation of molecular fluids. 1. A single-site coarse grained model of carbon dioxide.

An application of the "top-down" concept for the development of accurate coarse-grained intermolecular potentials of complex fluids is presented. With the more common "bottom-up" procedure, coarse-grained models are constructed from a suitable simplification of a detailed atomistic representation, and small adjustments to the intermolecular parameters are made by comparison with limited experimental data where necessary. In contrast, in the top-down approach, a molecular-based equation of state is used to obtain an effective coarse-grained intermolecular potential that reproduces the macroscopic experimental thermophysical properties over a wide range of conditions. These coarse-grained intermolecular potentials can then be used in conventional molecular simulation to obtain properties (such as structure or dynamics) that are not directly accessible from the equation of state or at extreme conditions where the theory is expected to fail. To demonstrate our procedure, a coarse-grained model for carbon dioxide (CO(2)) is obtained from a recent implementation of the statistical associating fluid theory of variable range (SAFT-VR) employing a Mie (generalized Lennard-Jones) potential; the parameters of this single-site Mie model of CO(2) are estimated by optimizing the equation of state's description of the experimental vapor-pressure and saturated liquid density data. This approach is only viable because of the excellent agreement of the SAFT-VR Mie EoS with simulation data. Our single-site SAFT-γ coarse-grained model for CO(2) is used in Monte Carlo molecular simulation to assess the adequacy of the description of the fluid-phase behavior and properties that were not used to develop the potential model such as the enthalpy of vaporization, interfacial tension, density profiles, supercritical densities, and second-derivative thermodynamic properties (thermal expansivity, isothermal compressibility, heat capacity, Joule-Thompson coefficient, and speed of sound). The accuracy of the description with the single-site SAFT-γ model of CO(2) is found to be of similar quality to that of more sophisticated intermolecular potentials such as a six-site (three LJ centers and three charged sites) all-atom model. The SAFT-γ top-down approach to coarse-graining resolves a key challenge with coarse-graining techniques: the provision of a direct robust link between the microscopic and macroscopic scales.

Simultaneous application of the gradient theory and Monte Carlo molecular simulation for the investigation of methane/water interfacial properties.

This work is dedicated to the simultaneous application of the gradient theory of fluid interfaces and Monte Carlo molecular simulations for the description of the interfacial behavior of the methane/water mixture. Macroscopic (interfacial tension, adsorption) and microscopic (density profiles, interfacial thickness) properties are investigated. The gradient theory is coupled in this work with the SAFT-VR Mie equation of state. The results obtained are compared with Monte Carlo simulations, where the fluid interface is explicitly considered in biphasic simulation boxes at both constant pressure and volume (NPT and NVT ensembles), using reliable united atom molecular models. On one hand, both methods provide very good estimations of the interfacial tension of this mixture over a broad range of thermodynamic conditions. On the other hand, microscopic properties computed with both gradient theory and MC simulations are in very good agreement with each other, which confirms the consistency of both approaches. Interfacial tension minima at high pressure and prewetting transitions in the vicinity of saturation conditions are also investigated.

On the isobaric thermal expansivity of liquids.

The temperature and pressure dependence of isobaric thermal expansivity, α(p), in liquids is discussed in this paper. Reported literature data allow general trends in this property that are consistent with experimental evidence to be established. Thus, a negative pressure dependence is to be expected except around the critical point. On the other hand, α(p) exhibits broad regions of negative and positive temperature dependence in the (T, p) plane depending on the nature of the particular liquid. These trends are rationalized here in terms of various molecular-based equations of state. The analysis of the Lennard-Jones, hard sphere square well and restricted primitive model equations allows understanding the differences in the α(p) behavior between liquids of diverse chemical nature (polar, nonpolar, and ionic): broader regions of negative temperature and positive pressure dependencies are obtained for liquids characterized by larger ranges of the interparticle potential. Also, using the statistical associating fluid theory (SAFT) allowed the behavior of more complex systems (basically, those potentially involving chain and association effects) to be described. The effect of chain length is rather simple: increasing it is apparently equivalent to raise the interaction range. By contrast, association presents a quite complex effect on α(p), which comes from a balance between the dispersive and associative parts of the interaction potential. Thus, if SAFT parameters are adjusted to obtain low association ability, α(p) is affected by each mechanism at clearly separate regions, one at low temperature, due to association, and the other to dispersive forces, which has its origin in fluctuations related with vapor-liquid transition.

Interfacial properties of water/CO2: a comprehensive description through a Gradient Theory-SAFT-VR Mie approach.

The Gradient Theory of fluid interfaces is for the first time combined with the SAFT-VR Mie EOS to model the interfacial properties of the water/CO(2) mixture. As a preliminary test of the performance of the coupling between both theories, liquid-vapor interfacial properties of pure water have been determined. The complex temperature dependence of the surface tension of water can be accurately reproduced, and the interfacial thickness is in good agreement with experimental data and simulation results. The water/CO(2) mixture presents several types of interfaces as the liquid water may be in contact with gaseous, liquid, or supercritical CO(2). Here, the interfacial tension of the water/CO(2) mixture is modeled accurately by the gradient theory with a unique value of the crossed influence parameter over a broad range of thermodynamic conditions. The interfacial density profiles show a systematic adsorption of CO(2) in the interface. Moreover, when approaching the saturation pressure of CO(2), a prewetting transition is highlighted. The adsorption isotherm of CO(2) is computed as well in the case of a gas/liquid interface and compared with experimental data. The good agreement obtained is an indirect proof of the consistency of interfacial density profiles computed with the gradient theory for this mixture and confirms that the gradient theory is suitable and reliable to describe the microstructure of complex fluid interfaces.

Interfacial properties of the Mie n-6 fluid: Molecular simulations and gradient theory results.

In a first part, interfacial properties of a pure monoatomic fluid interacting through the Mie n-6 potential (n=8, 10, 12, and 20) have been studied using extensive molecular simulations. Monte Carlo and molecular dynamics simulations have been employed, using, respectively, the test area approach and the mechanic route. In order to yield reference values, simulations have been performed with a cutoff radius equal to 10sigma, which is shown to be sufficient to avoid long range corrections. It is shown that both approaches provide results consistent with each other. Using the molecular simulations results, it is demonstrated that a unique scaling law is able to provide an accurate estimation of the surface tension whatever the repulsive exponent n, even far from the critical point. Furthermore, it is shown that the surface tension of the Mie n-6 fluid is as well accurately described by a unique Parachor's law. Density profiles are shown to be well represented by the tanh mean field profile, with slight deviations for the lowest temperatures and the smallest n. In addition, the interfacial width is shown to increase when n decreases (for a given reduced temperature) and to follow the usual scaling behavior for not too low temperature. In a second part, interfacial properties of the Mie n-6 fluid computed by the gradient theory, coupled with an equation of state based on the Barker-Henderson perturbation theory, have been compared with those obtained by molecular simulations. It is demonstrated that, even far from the critical point, the gradient theory is efficient to compute surface tensions and density profiles of this model fluid, provided the equation of state accurately model the phase behavior of the fluid involved (which is not the case for n=8 in this study).

Thermodynamic properties of the Mie n-6 fluid: a comparison between statistical associating fluid theory of variable range approach and molecular dynamics results.

Molecular dynamics (MD) simulations of direct and derivative thermodynamic properties of the Mie n-6 fluid (n=8, 10, and 12) have been performed for liquid to supercritical states. Using the results, an in depth test of the monomer-monomer interaction estimation of a recently derived statistical associating fluid theory of variable range (SAFT-VR) equation of state [Lafitte et al., J. Chem. Phys., 124, 024509 (2006)] has been carried out based on the Mie n-6 potential. For pure fluids, using an appropriate scaling, MD simulations show that density and isometric heat capacity are nearly independent of n, whereas sound velocity and thermal pressure coefficient tend to increase with n. In addition, the results show that predictions provided by the equation of state are consistent with those coming from MD and catch correctly the trends of each property with n except for the heat capacity. The comparison is next extended to binary mixtures with components differing only in the value of the n parameter and which demonstrate the reliability of the scheme (MX1b) used by Lafitte et al. to deal with this parameter in the SAFT-VR equation of state. In addition, a new empirical one-fluid approximation of the n parameter is proposed thanks to MD simulations, which very favorably compare with the one-fluid model on n previously proposed in the literature. The consistency of this approximation is addressed by making use of it in combination with the SAFT-VR Mie equation of state. It is shown that using such an approach, which is easier to handle than the MX1b one, yields slightly improved results compared to those of the MX1b.

A comprehensive description of chemical association effects on second derivative properties of alcohols through a SAFT-VR approach.

A recently derived version of the statistical associating fluid theory (SAFT), denoted as SAFT-VR Mie, which incorporates the Mie potentials within the SAFT-VR framework to model the monomer segment interactions (Lafitte et al. J. Chem. Phys. 2006, 124, 024509), is used for the study of second-order derivative properties and phase equilibria of alcohols and 1-alcohol + n-alkane binary mixtures. For this purpose, a variable repulsive potential is used to induce nonconformal interactions in the reference nonbonded fluid. These features have a significant influence on the chain and association contributions through the contact value of the radial distribution function, and they enhance the SAFT theory performance in the application to associating substances. When dealing with pure alcohols and 1-alcohol + n-alkane binary mixtures, an accurate description of both phase equilibria and second-order derivatives is obtained with a single set of molecular parameters. To explore the predictive ability limit of the model we have particularly focused our attention on secondary derivative properties, which display singularities due to the formation of aggregates. With this approach, we have found that the model is able to reproduce accurately the complex behavior of the isobaric heat capacity of alcohols as, for instance, the maximum versus temperature in the compressed liquid region. Furthermore, in the case of 1-hexanol + n-hexane binary mixtures, the proposed equation is found to capture the association effects on the pressure and temperature dependence of the isobaric thermal expansivity. These two special features, which to our knowledge have never been described by a theoretical model, emphasize both the validity of the changes in the model proposed and the physical meaning of the molecular parameters obtained in this study.

Simultaneous estimation of phase behavior and second-derivative properties using the statistical associating fluid theory with variable range approach.

A modified statistical associating fluid theory (SAFT) with variable range version is presented using the family of m-n Mie potentials. The use of this intermolecular potential for modeling repulsion-dispersion interactions between the monomer segments, together with a new method for optimizing the molecular parameters of the equation of state, is found to give a very accurate description of both vapor-liquid equilibria and compressed liquid bulk properties (volumetric and derivative properties) for long-chain n-alkanes. This new equation improves other SAFT-like equations of state which fail to describe derivative properties such as the isothermal compressibility and the speed of sound in the condensed liquid phase. Emphasis is placed on pointing out that the key for modeling the latter properties is the use of a variable repulsive term in the intermolecular potential. In the case of the n-alkanes series, a clear dependence of the characteristic molecular parameters on increasing chain length is obtained, demonstrating their sound physical meaning and the consistency of the new fitting procedure proposed. This systematic method for optimizing the model parameters includes data on the saturation line as well as densities and speed of sound data in the condensed liquid phase, and the results show undoubtedly that the model performance is enhanced and its range of applicability is now widened, keeping in any case a good balance between the accuracy of the different estimated properties.