Carmelo Herdes

Computer simulation of volatile organic compound adsorption in atomistic models of molecularly imprinted polymers.

Molecularly imprinted polymers (MIPs) offer a unique opportunity to significantly advance volatile organic compound (VOC) sensing technologies and a number of other applications. However, the development of these applications using MIPs has been hindered by poor understanding of the microstructure of MIPs, geometry of binding sites, and the details of molecular recognition processes in these materials. This is further complicated by the vast number of optimization parameters such as building components and processing conditions. Computer simulations and molecular modeling can help us understand adsorption and binding phenomena in MIPs on the molecular level and thus provide a route to more efficient MIP design strategies. So far, molecular models have been either oversimplified or severely limited in length scale, essentially focusing on a single binding site. Here, we propose a more general, atomistically detailed model that describes the microstructure of MIPs. We apply this model to investigate adsorption of pyridine, benzene, and toluene in MIPs and demonstrate that it is able to capture a number of essential experimental features. Therefore, this model can serve as a starting point in computational design and optimization of MIPs.

On the Calculation of Solid-Fluid Contact Angles from Molecular Dynamics

A methodology for the determination of the solid-fluid contact angle, to be employed within molecular dynamics (MD) simulations, is developed and systematically applied. The calculation of the contact angle of a fluid drop on a given surface, averaged over an equilibrated MD trajectory, is divided in three main steps: (i) the determination of the fluid molecules that constitute the interface, (ii) the treatment of the interfacial molecules as a point cloud data set to define a geometric surface, using surface meshing techniques to compute the surface normals from the mesh, (iii) the collection and averaging of the interface normals collected from the post-processing of the MD trajectory. The average vector thus found is used to calculate the Cassie contact angle (i.e., the arccosine of the averaged normal z-component). As an example we explore the effect of the size of a drop of water on the observed solid-fluid contact angle. A single coarse- grained bead representing two water molecules and parameterized using the SAFT-γ Mie equation of state (EoS) is employed, meanwhile the solid surfaces are mimicked using integrated potentials. The contact angle is seen to be a strong function of the system size for small nano-droplets. The thermodynamic limit, corresponding to the infinite size (macroscopic) drop is only truly recovered when using an excess of half a million water coarse-grained beads and/or a drop radius of over 26 nm.

Modelling the interfacial behaviour of dilute light-switching surfactant solutions

The direct molecular modelling of an aqueous surfactant system at concentrations below the critical micelle concentration (pre-cmc) conditions is unviable in terms of the presently available computational power. Here, we present an alternative that combines experimental information with tractable simulations to interrogate the surface tension changes with composition and the structural behaviour of surfactants at the water-air interface. The methodology is based on the expression of the surface tension as a function of the surfactant surface excess, both in the experiments and in the simulations, allowing direct comparisons to be made. As a proof-of-concept a coarse-grained model of a light switching non-ionic surfactant bearing a photosensitive azobenzene group is considered at the air-water interface at 298 K. Coarse-grained molecular dynamic simulations are detailed based on the use of the SAFT force field with parameters tuned specifically for this purpose. An excellent agreement is obtained between the simulation predictions and experimental observations; furthermore, the molecular model allows the rationalization of the macroscopic behaviour in terms of the different conformations of the cis and trans surfactants at the surface.

Molecular Recognition Effects in Atomistic Models of Imprinted Polymers

In this article we present a model for molecularly imprinted polymers, which considers both complexation processes in the pre-polymerization mixture and adsorption in the imprinted structures within a single consistent framework. As a case study we investigate MAA/EGDMA polymers imprinted with pyrazine and pyrimidine. A polymer imprinted with pyrazine shows substantial selectivity towards pyrazine over pyrimidine, thus exhibiting molecular recognition, whereas the pyrimidine imprinted structure shows no preferential adsorption of the template. Binding sites responsible for the molecular recognition of pyrazine involve one MAA molecule and one EGDMA molecule, forming associations with the two functional groups of the pyrazine molecule. Presence of these specific sites in the pyrazine imprinted system and lack of the analogous sites in the pyrimidine imprinted system is directly linked to the complexation processes in the pre-polymerization solution. These processes are quite different for pyrazine and pyrimidine as a result of both enthalpic and entropic effects.

Predicting neopentane isosteric enthalpy of adsorption at zero coverage in MCM-41

The isosteric enthalpy of adsorption for neopentane at relative pressures down to 3 × 10(-8) in MCM-41 was predicted for the temperature range from -15 to 0 °C. At such low pressures and temperatures, experimental measurements become problematic for this system. We used an atomistic model for MCM-41 obtained by means of a kinetic Monte Carlo method mimicking the synthesis of the material. The model was parametrized to represent experimental nitrogen adsorption isotherms at 77 K using grand canonical Monte Carlo simulations. The simulated isosteric enthalpy of adsorption shows very good agreement with available experimental data, demonstrating that GCMC simulations can predict heats of adsorption for conditions that are challenging for experimental measurements. Additional insights into the adsorption mechanisms, derived from energetic analysis at the molecular level, are also presented.

Thermodynamic properties and aggregate formation of surfactant-like molecules from theory and simulation

The goal of this work is twofold: to predict the phase equilibria behavior of simplified surfactant models and to predict the population of aggregates as a function of pressure. We compare Monte Carlo simulation results of these systems with predictions from a modified version of the statistical associating fluid theory (soft-SAFT). Surfactant-like molecules are modeled as Lennard-Jones chains of tangent segments with one or two association sites. We study the influence of the number and location of the association sites on the thermodynamic properties and fraction of nonbonded molecules in all cases. The influence of the chain length is also investigated for a particular location of the sites. Results are compared with NPT Monte Carlo simulations to test the accuracy of the theory, and to study the molecular configurations of the system. Soft-SAFT is able to quantitatively predict the MC PVT results, independently of the location of the association sites. The theory is also able to capture the qualitative trend of the population of aggregates with pressure. Quantitative agreement is only obtained for specific locations of the sites.

Predicting the adsorption of n-perfluorohexane in BAM P109 standard activated carbon by molecular simulation using SAFT-gamma Mie coarse-grained force fields

This work is framed within the Eighth Industrial Fluid Properties Simulation Challenge, with the aim of assessing the capability of molecular simulation methods and force fields to accurately predict adsorption in porous media for systems of relevant practical interest. The current challenge focuses on predicting adsorption isotherms of n-perfluorohexane in the certified reference material BAM-P109 standard activated carbon. A temperature of T = 273  K and pressures of p / p 0 = 0 . 1 , 0.3, and 0.6 relative to the bulk saturation pressure p0 (as predicted by the model) are the conditions selected in this challenge. In our methodology we use coarse-grained intermolecular models and a top-down technique where an accurate equation of state is used to link the experimental macroscopic properties of a fluid to the force-field parameters. The state-of-the-art version of the statistical associating fluid theory (SAFT) for potentials of variable range as reformulated in the Mie group contribution incarnation (SAFT-γ Mie) is employed here. The parameters of the SAFT-γ Mie force field are estimated directly from the vapour pressure and saturated liquid density data of the pure fluids using the equation of state, and further validated by molecular dynamic simulations. The coarse-grained intermolecular potential models are then used to obtain the adsorption isotherm kernels for argon, carbon dioxide, and n-perfluorohexane in graphite slit pores of various widths using Grand Canonical Monte Carlo simulations. A unique and fluid-independent pore size distribution curve with total micropore volume of 0.5802 cm3/g is proposed for the BAM-P109. The pore size distribution is obtained by applying a non-linear regression procedure over the adsorption integral equation to minimise the quadratic error between the available experimental adsorption isotherms for argon and carbon dioxide and purpose-built Grand Canonical Monte Carlo kernels. The predicted adsorption levels of n-perfluorohexane at 273 K in BAM-P109 are 72.75 ± 0.01, 73.82 ± 0.01, and 75.44 ± 0.05 cm3/g at Standard Temperature and Pressure (STP) conditions for p / p 0 = 0 . 1 , 0.3, and 0.6, respectively.

Fundamental studies of methyl iodide adsorption in DABCO impregnated activated carbons

Methyl iodide capture from a water vapor stream using 1,4-diazabicyclo[2.2.2]octane (DABCO)-impregnated activated carbons is, for the first time, fundamentally described here on the atomic level by means of both molecular dynamics and grand canonical Monte Carlo simulations. A molecular dynamics annealing strategy was adopted to mimic the DABCO experimental impregnation procedure in a selected slitlike carbon pore. Predictions, restricted to the micropore region, are made about the adsorption isotherms of methyl iodide, water, and nitrogen on both impregnated and bare activated carbon models. Experimental and simulated nitrogen adsorption isotherms are compared for the validation of the impregnation strategy. Selectivity analyses of the preferential adsorption toward methyl iodide over water are also reported. These simulated adsorption isotherms sum up to previous experimental studies to provide an enhanced picture for this adsorption system of widespread use at nuclear plant HVAC facilities for the capture of radioactive iodine compounds.

A molecular-based equation of state for process engineering

Abstract We outline here some of the steps we are taking towards the development of reliable tools for quantitative predictions of thermodynamic properties of complex fluids with equations based on statistical mechanics. The long term objective is to provide a user-friendly computer code and a wide database of molecular parameters for different compounds, able to be implemented in a process simulator. We have observed that the keys of the success when using molecular modelling tools for predictions rely on the selection of the appropriate model, representative of the molecular structure, and the use of physically meaningful molecular parameters.

A multiscale method for simulating fluid interfaces covered with large molecules such as asphaltenes

Abstract The interface between two liquids is fully described by the interfacial tension only for very pure liquids. In most cases the system also contains surfactant molecules which modify the interfacial tension according to their concentration at the interface. This has been widely studied over the years, and interesting phenomena arise, e.g. the Marangoni effect. An even more complicated situation arises for complex fluids like crude oil, where large molecules such as asphaltenes migrate to the interface and give rise to further phenomena not seen in surfactant-contaminated systems. An example of this is the “crumpling drop” experiments, where the interface of a drop being deflated becomes non-smooth at some point. In this paper we report on the development of a multiscale method for simulating such complex liquid–liquid systems. We consider simulations where water drops covered with asphaltenes are deflated, and reproduce the crumpling observed in experiments. The method on the nanoscale is based on using coarse-grained molecular dynamics simulations of the interface, with an accurate model for the asphaltene molecules. This enables the calculation of interfacial properties. These properties are then used in the macroscale simulation, which is performed with a two-phase incompressible flow solver using a novel hybrid level-set/ghost-fluid/immersed-boundary method for taking the complex interface behaviour into account. We validate both the nano- and macroscale methods. Results are presented from nano- and macroscale simulations which showcase some of the interesting behaviour caused by asphaltenes affecting the interface. The molecular simulations presented here are the first in the literature to obtain the correct interfacial orientation of asphaltenes. Results from the macroscale simulations present a new physical explanation of the crumpled drop phenomenon, while highlighting shortcomings in previous hypotheses.