Adsorption and dynamics of linear and mono-branched hexane isomers in MIL-140 metal-organic frameworks.

Recent breakthrough experiments revealed the iso-reticular Zr-MOFs, MIL-140B and MIL-140C, as promising sorbents for the separation of C6 isomers. Interestingly while the ultra-small pore MIL-140B exhibited hexane isomer sorption hierarchy according to the normal boiling point order (n-C6 > 3MP (3-methyl pentane)), an uncommon shift in the elution order was observed in the larger pore MIL-140C. It was only speculated that the flexibility of the MOFs might be the origin of this intriguing behavior. Herein, flexible force field hybrid osmotic Monte Carlo combined with molecular dynamics simulations were carried out to unravel the microscopic mechanism of the adsorption and dynamics of both C6 isomers in MIL140B and MIL140C. Thermodynamically preferred adsorption of n-C6 over 3MP was predicted for MIL-140B and to a slightly less extent for MIL-140C. Interestingly while the mobility of n-C6 was found to remain higher than that of 3MP in the whole range of loading for MIL-140B, 3MP becomes more mobile than n-C6 at saturation in MIL-140C. This suggests that this kinetics order is most probably the origin of the inversion of the elution order observed experimentally for MIL-140C. The translational and rotational dynamics of the two guests in MIL-140B and MIL-140C was further understood in-depth.

Associated molecular liquids at the graphene monolayer interface.

We report molecular simulations of the interaction between a graphene sheet and different liquids such as water, ethanol, and ethylene glycol. We describe the structural arrangements at the graphene interface in terms of density profiles, number of hydrogen bonds (HBs), and local structuration in neighboring layers close to the surface. We establish the formation of a two-dimensional HB network in the layer closest to the graphene. We also calculate the interfacial tension of liquids with a graphene monolayer and its profile along the direction normal to the graphene to rationalize and quantify the strengthening of the intermolecular interactions in the liquid due to the presence of the surface.

Tuning the hexane isomer separation performances of Zeolitic Imidazole Framework-8 using mechanical pressure.

Hybrid osmotic Monte Carlo simulations were performed to anticipate the tunability of the separation performance of the flexible Zeolitic Imidazole Framework-8 (ZIF-8) via the application of an external mechanical pressure. This synergistic combination of mechanical control of the pore aperture/cage dimension and guest adsorption was applied to the challenging hexane isomers separation processes of vital importance in the field of petrochemical industry. The application of a mechanical pressure above 1 GPa was predicted to boost the linear hexane/2-methylpentane and 2-methylpentane/2,3-dimethylbutane selectivity by 40% and 17%, respectively, as compared to the pristine ZIF-8. We further unraveled the microscopic origin of this optimized performance with an in-depth analysis of the critical interplay between the structural changes of the ZIF-8 framework and the conformational rearrangements of C6 isomers under mechanical pressure.

Dynamics of water confined in mesopores with variable surface interaction.

We have investigated the dynamics of liquid water confined in mesostructured porous silica (MCM-41) and periodic mesoporous organosilicas (PMOs) by incoherent quasielastic neutron scattering experiments. The effect of tuning the water/surface interaction from hydrophilic to more hydrophobic on the water mobility, while keeping the pore size in the range 3.5 nm-4.1 nm, was assessed from the comparative study of three PMOs comprising different organic bridging units and the purely siliceous MCM-41 case. An extended dynamical range was achieved by combining time-of-flight (IN5B) and backscattering (IN16B) quasielastic neutron spectrometers providing complementary energy resolutions. Liquid water was studied at regularly spaced temperatures ranging from 300 K to 243 K. In all systems, the molecular dynamics could be described consistently by the combination of two independent motions resulting from fast local motion around the average molecule position and the confined translational jump diffusion of its center of mass. All the molecules performed local relaxations, whereas the translational motion of a fraction of molecules was frozen on the experimental timescale. This study provides a comprehensive microscopic view on the dynamics of liquid water confined in mesopores, with distinct surface chemistries, in terms of non-mobile/mobile fraction, self-diffusion coefficient, residence time, confining radius, local relaxation time, and their temperature dependence. Importantly, it demonstrates that the strength of the water/surface interaction determines the long-time tail of the dynamics, which we attributed to the translational diffusion of interfacial molecules, while the water dynamics in the pore center is barely affected by the interface hydrophilicity.

Calculation of the interfacial tension of the graphene-water interaction by molecular simulations.

We report the calculation of the solid-liquid interface tension of the graphene-water interaction by using molecular simulations. Local profiles of the interfacial tension are given through the mechanical and thermodynamic definitions. The dependence of the interfacial tension on the graphene area is investigated by applying both reaction field and Ewald summation techniques. The structure of the interfacial region close to the graphene sheet is analyzed through the profiles of the density and hydrogen bond number and the orientation of the water molecules. We complete this study by plotting the profiles of the components of the pressure tensor calculated by the Ewald summation and reaction field methods. We also investigate the case of a reaction field version consisting in applying a damped shifted force in the case of the calculation of the pressure components.

The extent of the glass transition from molecular simulation revealing an overcrank effect.

A deep understanding of the transition between rubber and amorphous state characterized by a glass transition temperature, Tg , is still a source of discussions. In this work, we highlight the role of molecular simulation in revealing explicitly this temperature dependent behavior. By reporting the specific volume, the thermal expansion coefficient and the heat capacity versus the temperature, we actually show that the glass transition domain extends to a greater range of temperature, compared with experiments. This significant enlargement width is due to the fast cooling rate, and actually explains the difficulty to locate Tg . This result is the manifestation of an overcranking effect used by high-speed cameras to reveal slow-motion. Accordingly, atomistic simulation offers the significant opportunity to show that the transition from the rubber state to the glass phase should be detailed in terms of the degrees of freedom freeze. © 2017 Wiley Periodicals, Inc.

Interactions of Organics within Hydrated Selective Layer of Reverse Osmosis Desalination Membrane: A Combined Experimental and Computational Study.

In this work we have examined a computational approach in predicting the interactions between uncharged organic solutes and polyamide membranes. We used three model organic molecules with identical molecular weights (100.1 g/mol), 4-aminopiperidine, 3,3-dimethyl-2-butanone (pinacolone) and methylisobutyl ketone for which we obtained experimental data on partitioning, diffusion and separation on a typical seawater reverse osmosis (RO) membrane. The interaction energy between the solutes and the membrane phase (fully aromatic polyamide) was computed from molecular dynamics (MD) simulations and the resulting sequence was found to correlate well with the experimental rejections and sorption data. Sorption of the different organic solutes within the membrane skin layer determined from attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR) nicely agreed with interaction energies computed from molecular simulations. Qualitative information about solute diffusivity inside the membrane was also extracted from MD simulations while ATR-FTIR experiments indicated strongly hindered diffusion with diffusion coefficients in the membrane about 10-15 m2/s. The computational approach presented here could be a first step toward predicting rejections trends of, for example, hormones and pharmaceuticals by RO dense membranes.

Importance of the tail corrections on surface tension of curved liquid-vapor interfaces.

We report molecular simulations of the liquid-vapor cylindrical interface of methane. We apply the truncated Lennard-Jones potential and specific long-range corrections for the surface tension developed especially for cylindrical interfaces. We investigate the impact of the cutoff on the radial density profile, the intrinsic and long-range correction parts to the surface tension, and Tolman length. We also study the curvature dependence of the surface tension as a function of the cutoff used. In this work we shed light that both density and Tolman length are cutoff-dependent whereas the total surface tension is slightly curvature and cutoff dependent.

Computer modelling of the surface tension of the gas-liquid and liquid-liquid interface.

This review presents the state of the art in molecular simulations of interfacial systems and of the calculation of the surface tension from the underlying intermolecular potential. We provide a short account of different methodological factors (size-effects, truncation procedures, long-range corrections and potential models) that can affect the results of the simulations. Accurate calculations are presented for the calculation of the surface tension as a function of the temperature, pressure and composition by considering the planar gas-liquid interface of a range of molecular fluids. In particular, we consider the challenging problems of reproducing the interfacial tension of salt solutions as a function of the salt molality; the simulations of spherical interfaces including the calculation of the sign and size of the Tolman length for a spherical droplet; the use of coarse-grained models in the calculation of the interfacial tension of liquid-liquid surfaces and the mesoscopic simulations of oil-water-surfactant interfacial systems.

Can we approach the gas-liquid critical point using slab simulations of two coexisting phases?

In this paper, we demonstrate that it is possible to approach the gas-liquid critical point of the Lennard-Jones fluid by performing simulations in a slab geometry using a cut-off potential. In the slab simulation geometry, it is essential to apply an accurate tail correction to the potential energy, applied during the course of the simulation, to study the properties of states close to the critical point. Using the Janecek slab-based method developed for two-phase Monte Carlo simulations [J. Janec̆ek, J. Chem. Phys. 131, 6264 (2006)], the coexisting densities and surface tension in the critical region are reported as a function of the cutoff distance in the intermolecular potential. The results obtained using slab simulations are compared with those obtained using grand canonical Monte Carlo simulations of isotropic systems and the finite-size scaling techniques. There is a good agreement between these two approaches. The two-phase simulations can be used in approaching the critical point for temperatures up to 0.97 TC∗ (T∗ = 1.26). The critical-point exponents describing the dependence of the density, surface tension, and interfacial thickness on the temperature are calculated near the critical point.

Theoretical investigation of the ionic selectivity of polyelectrolyte multilayer membranes in nanofiltration.

Polyelectrolyte multilayer membranes have proven to be promising materials for ion nanofiltration. In this work, we implement a continuum mesoscopic transport model developed in previous works (Szymczyk, A.; Zhu, H.; Balannec, B. Langmuir 2010, 26, 1214; Szymczyk, A.; Zhu, H.; Balannec, B. J. Phys. Chem. B 2010, 114, 10143) to investigate the pressure-driven transport of electrolyte mixtures through this kind of membrane. The model accounts for an inhomogeneous distribution of the fixed charge through an arbitrary number of polyelectrolyte bilayers. We show that accounting for the multiple bipolar charge distribution resulting from the layer-by-layer assembly of polyelectrolytes with opposite charge is responsible for the increase in the Na(+)/Mg(2+) selectivity reported experimentally with respect to conventional nanofiltration membranes. The model also allows the rationalizing of the seemingly contradictory experimental results reported in the literature (i.e., the increase or decrease in the selectivity with the number of bilayers or the existence of an optimum number of bilayers). It is shown, however, that the nonmonotonous variation of the ionic selectivity does not originate from the multibipolar distribution of the fixed charge through polyelectrolyte multilayer membranes but from the existence of an optimum skin layer thickness.

Superpermittivity of nanoconfined water.

Nowadays, it is well established that the physical properties of confined liquids strongly differ from those in bulk phase. While dynamical and structural properties were strongly explored, dielectric properties are poorly studied despite their importance in the understanding and the modelling of molecular mechanism in a number of nano-applications such as nanofluidics, nanofiltration, and nanomedicine. Among them, the dielectric permittivity is probably one of the most important. The lack of knowledge about it strongly limits our ability to model fluid-material interactions and more generally our understanding of the behaviour of confined fluids. Recently, the dielectric permittivity of confined water in silica, Metal Organic Frameworks, and graphene materials was found to be slightly higher than the permittivity of water in bulk phase. In this work, the permittivity of water and dichloromethane confined in carbon nanotubes was predicted by means of molecular dynamics simulations. The static dielectric constant was found to be 700, i.e., 10-fold higher than the bulk value. This superpermittivity has, for origin, the excluded volume and the presence of an unconfined direction leading to a pre-orientation of water molecules close to the pore wall and an increase in dipolar fluctuations.

Thermotropic luminescent clustomesogen showing a nematic phase: a combination of experimental and molecular simulation studies.

Octahedral Mo6 nanoclusters are functionalized with two organic ligands containing cyanobiphenyl (CB) units, giving luminescent hybrid liquid crystals (LC). Although the mesogenic density around the bulky inorganic core is constant, the two hybrids show different LC properties. Interestingly, one of them shows a nematic phase, which is particularly rare for this kind of supermolecular system. This surprising result is explained by using large-scale molecular dynamic simulations.

Water confinement in nanoporous silica materials.

The influence of the surface polarity of cylindrical silica nanopores and the presence of Na(+) ions as compensating charges on the structure and dynamics of confined water has been investigated by molecular dynamics simulations. A comparison between three different matrixes has been included: a protonated nanopore (PP, with SiOH groups), a deprotonated material (DP, with negatively charged surface groups), and a compensated-charge framework (CC, with sodium cations compensating the negative surface charge). The structure of water inside the different pores shows significant differences in terms of layer organization and hydrogen bonding network. Inside the CC pore the innermost layer is lost to be replaced by a quasi bulk phase. The electrostatic field generated by the DP pore is felt from the surface to the centre of pore leading to a strong orientation of water molecules even in the central part of the pore. Water dynamics inside both the PP and DP pores shows significant differences with respect to the CC pore in which the sub-diffusive regime of water is lost for a superdiffusive regime.

Discotic columnar liquid crystal studied in the bulk and nanoconfined states by molecular dynamics simulation.

A prototypical Gay Berne discotic liquid crystal was studied by means of molecular dynamics simulations both in the bulk state and under confinement in a nanoporous channel. The phase behavior of the confined system strongly differs from its bulk counterpart: the bulk isotropic-to-columnar transition is replaced by a continuous ordering from a paranematic to a columnar phase. Moreover, a new transition is observed at a lower temperature in the confined state, which corresponds to a reorganization of the intercolumnar order. It reflects the competing effects of pore surface interaction and genuine hexagonal packing of the columns. The translational molecular dynamics in the different phases has been thoroughly studied and discussed in terms of collective relaxation modes, non-Gaussian behavior, and hopping processes.

Crossover in structure and dynamics of a primary alcohol induced by hydrogen-bonds dilution.

Primary alcohols show a prominent Debye process in the dielectric relaxation located at a timescale longer than the main structural relaxation. Böhmer and co-workers studied dilution effects of the hydrogen bonding network of n-butanol (BuOH) with n-bromobutane (BuBr) to better understand the origin of this process. Interestingly, this work has evidenced a crossover in Debye relaxation time (τD) for a critical concentration in BuBr xc = 0.5. By using molecular dynamics simulations and NMR experiments we propose here to explore further dilution effects on the dipolar and translational dynamics. Moreover, we discuss the relation between structural and dynamical properties in the context of a detailed study of the microstructure and the H-bond network. The overall results are consistent with the existence of a topological change in the liquid structure occurring at about xc = 0.5 from a hypernetted percolating network to independent nanodomains of n-butanol molecules embedded in the n-bromobutane phase.

Multiscale modelling of transport in polymer-based reverse-osmosis/nanofiltration membranes: present and future.

Nanofiltration (NF) and reverse osmosis (RO) processes are physical separation technologies used to remove contaminants from liquid streams by employing dense polymer-based membranes with nanometric voids that confine fluids at the nanoscale. At this level, physical properties such as solvent and solute permeabilities are intricately linked to molecular interactions. Initially, numerous studies focused on developing macroscopic transport models to gain insights into separation properties at the nanometer scale. However, continuum-based models have limitations in nanoconfined situations that can be overcome by force field molecular simulations. Continuum-based models heavily rely on bulk properties, often neglecting critical factors like liquid structuring, pore geometry, and molecular/chemical specifics. Molecular/mesoscale simulations, while encompassing these details, often face limitations in time and spatial scales. Therefore, achieving a comprehensive understanding of transport requires a synergistic integration of both approaches through a multiscale approach that effectively combines and merges both scales. This review aims to provide a comprehensive overview of the state-of-the-art in multiscale modeling of transport through NF/RO membranes, spanning from the nanoscale to continuum media.

Importance of the Electrostatic Correlations in Surface Tension of Hydrated Reline Deep Eutectic Solvent from Combined Experiments and Molecular Dynamics Simulations.

In this study, the surface tension and the structure of hydrated reline are investigated by using diverse methods. Initially, the surface tension displays a nonlinear pattern as water content increases, decreasing until reaching 45 wt %, then gradually matching that of pure water. This fluctuation is associated with strong electrostatic correlations present in pure reline, which decrease as more water is added. Changes in surface tension reflect a shift from charge layering in pure reline to an increased interfacial hydrogen bonding as the water content rises. This shift causes the segregation of urea molecules into the bulk phase and a gradual anchoring of water molecules to the air-reline interface. An interesting observation is the antisurfactant effect, where heightened interfacial anchoring results in an unexpected increase in real contribution of surface tension. This, along with weakened electrostatic correlations beyond 45 wt % due to reinforced interfacial hydrogen bonding, contributes to the complex behavior of surface tension observed in this study.

Recent advances in Many Body Dissipative Particles Dynamics simulations of liquid-vapor interfaces.

Many Body Dissipative Particles Dynamics (MDPD) simulation is a novel promising mesoscopic method to model the liquid-vapor interfaces. Based upon works of Paganobarraga and Frenkel (J. Chem. Phys. 15, 5015 (2001)) and Trofimov (J. Chem. Phys. 117, 9383 (2002)) and of Warren (Phys. Rev. E 68, 066702 (2003)) this method has been critically reviewed during this last decade. We propose here to give an overview of the Many Body Dissipative Particles Dynamic simulation within the framework of the liquid-vapor interfaces. We recall the theoretical background of MDPD and we present some recent results of systems of interest such as water liquid-vapor interfaces and salt effect on water surface tension. Additionally we discuss the ability of MDPD to capture the mechanisms at the mesoscopic scale through the formation of micelles and the coalescence of a nanodroplet water on water surface.

Prediction of the concentration dependence of the surface tension and density of salt solutions: atomistic simulations using Drude oscillator polarizable and nonpolarizable models.

Molecular simulations using Drude oscillator polarizable and nonpolarizable models for water and ions are carried out to predict the dependence of the surface tension on salt concentration. The polarizable water and ion models are based only on the classical Drude oscillators. The temperature dependence of the surface tension of water is examined for different water models. The dependence of salt densities on salt concentration is investigated through the nonpolarizable and Drude oscillator polarizable models. Finally, the reproduction of the surface tension of salt solution over a large range of concentrations is analyzed through a number of combinations between ions and water force fields. The structure of the interface is then discussed as a function of polarization effects. We establish here the inability of the Drude oscillator polarizable force fields to reproduce the salt concentration dependence of surface tension of NaCl aqueous solutions.