On the orientation of N2and CO2molecules adsorbed in slit pore models with oxidised graphitic surface

Numerical Analysis of the Porous Structure of Activated Carbons Derived from Synthetic Polymers

This paper presents original results from the unique analysis of the porous structure of activated carbons (ACs) produced through the chemical activation of polyethylene terephthalate (PET) and polyacrylonitrile (PAN), as well as from a physical mixture of both polymers. An advanced method of adsorbent surface analysis-more specifically, the new method of numerical clustering-based adsorption analysis regarding the surface heterogeneity, pore geometry and adsorption energy distribution parameters-allowed us to obtain information about the porous structure of the ACs from the synthetic polymers mentioned above. As the results showed, ACs obtained with PAN were characterised by a first adsorbed layer with the highest volume. When the surface heterogeneity, highly desirable in most advanced adsorption processes, is taken into account, the materials with the best surface properties in both potassium carbonate (K2CO3) and potassium hydroxide (KOH) activation processes were the ACs obtained with a mass proportion of PET to PAN of 1:3, which were characterised by a low degree of surface heterogeneity and a first adsorbed layer presenting a relatively large volume.

Controlling the Thermal Conductivity of Monolayer Graphene with Kirigami Structure

In this work, the thermal conductivity performance of graphene kirigami (GK) was systematically investigated via molecular dynamics (MD) simulations. The results indicate that the degree of defects (DD) on GK has a significant influence on thermal conductivity. Reducing the DD is the most effective way to decrease the thermal conductivity of GK. For zigzag-incised GK sheets, the change rate of thermal conductivity can reach up to 1.86 W/mK per 1% change in DD by tuning the incision length. The rate of changing thermal conductivity with DD can be slowed down by changing the width among incisions. Compared with the zigzag-incised GK sheets, heat transfer across the armchair-incised GK comes out more evenly, without significant steep and gentle stages along the heat transfer routes. More importantly, the GK structure can adjust the thermal conductivity by stretching, which the previously reported nanoporous graphene does not have. The change rate of thermal conductivity achieves about 0.17 W/mK with 1% stretching strain for simulated GK and can be further reduced at high tensile strain rates, benefiting the precise and variable control of the thermal conductivity of the monolayer graphene.

Shape Transformations and Self-Assembly of Hairy Particles under Confinement

Molecular dynamics simulations are used to investigate the behavior of polymer-tethered nanoparticles between two inert or attractive walls. The confinement in pores creates new possibilities for controlling the shape transformation of individual hairy particles and their self-organization. We introduce a minimalistic model of the system; only chain-wall interactions are assumed to be attractive, while the others are softly repulsive. We show how the shape of isolated particles can be controlled by changing the wall separation and the strength of the interaction with the surfaces. For attractive walls, we found two types of structures, “bridges” and “mounds”. The first structures are similar to flanged spools in which the chains are connected with both walls and form bridges between them. We observed various bridges, symmetrical and asymmetrical spools, hourglasses, and pillars. The bridge-like structures can be “nano-oscillators” in which the cores jump from one wall to the other. We also study the self-assembly of a dense fluid of hairy particles in slit-like pores and analyze how the system morphology depends on interactions with the surfaces and the wall separation. The hairy particles form layers parallel to the walls. Different ordered structures, resembling two-dimensional crystalline lattices, are reported. We demonstrate that hairy particles are a versatile soft component forming a variety of structures in the slits.

Umbrella Sampling Simulations of Carbon Nanoparticles Crossing Immiscible Solvents

We use molecular dynamics to compute the free energy of carbon nanoparticles crossing a hydrophobic–hydrophilic interface. The simulations are performed on a biphasic system consisting of immiscible solvents (i.e., cyclohexane and water). We solvate a carbon nanoparticle into the cyclohexane layer and use a pull force to drive the nanoparticle into water, passing over the interface. Next, we accumulate a series of umbrella sampling simulations along the path of the nanoparticle and compute the solvation free energy with respect to the two solvents. We apply the method on three carbon nanoparticles (i.e., a carbon nanocone, a nanotube, and a graphene nanosheet). In addition, we record the water-accessible surface area of the nanoparticles during the umbrella simulations. Although we detect complete wetting of the external surface of the nanoparticles, the internal surface of the nanotube becomes partially wet, whereas that of the nanocone remains dry. This is due to the nanoconfinement of the particular nanoparticles, which shields the hydrophobic interactions encountered inside the pores. We show that cyclohexane molecules remain attached on the concave surface of the nanotube or the nanocone without being disturbed by the water molecules entering the cavity.

Binding Free Energy Calculations of Bilayer Graphenes Using Molecular Dynamics

Bilayer graphenes are dimeric assemblies of single graphene layers bound together by π-complexation interactions. Controlling these assemblies can be complicated, as the layered compounds disperse in solvents or aggregate into higher columnar configurations and clusters. One way to assess the interactions that contribute to the stability of the layered compounds is to use molecular simulation. We perform pulling molecular dynamics on bilayer graphenes with different sizes and obtain the normal and shear force profiles of dissociation. We generate pathways of dissociation along the two directions and calculate the binding free energies of the structures with umbrella sampling simulations. We show that the dissociation process is direction-dependent. Along the shear direction, we compute the same free energy for the different samples, which validates the consistency of our simulations. We notice that the dissociation is less adiabatic on the normal than the shear direction, having an entropic contribution to the Gibbs energy. This contribution is more enhanced for the larger bilayer graphenes.

Calculating adsorption isotherms using Lennard Jones particle density distributions

The computation of particle density distributions in pore channels is a fundamental post process practice in molecular adsorption simulations. The distributions, although not appropriate for direct experimental interrogation, when expressed in variable temperature, may be used to evaluate thermodynamic properties. As with molecular simulations, we can spotlight any frame or region of interest inside the computational cell, the distributions and subsequently the thermodynamic evaluations can be pore or site specific. This allows us to establish correlations for the adsorption capacity, on different pore partitions and surface textures and compare strengths of explicit interactions. We confirm this assumption using adsorption simulations of Lennard Jones particles such as argon in the dual cylindrical cavity of ZIF-69. We initiate the computations from the density distribution functions and reproduce full curves of adsorption isotherms which correlate well with the direct simulations. The computations are fast and implemented using linear van' t Hoff plots for the surface coverage at a sequence of pressures.

The effect of gme topology on multicomponent adsorption in zeolitic imidazolate frameworks

We employ a simulation approach to study the adsorption of single, binary and ternary mixtures on eight gme Zeolitic Imidazolate Frameworks (ZIFs) at 298 K. Four adsorbate fluids were considered; carbon dioxide, methane, nitrogen and water. We compute the high pressure adsorption density profiles inside the micropore channels of each crystal. The profiles are compared directly for the different structures and adsorbate components and used to highlight the influence of the imidazolate ligands on pure and competitive adsorption. ZIFs with long ligands reveal an additional, accessible to the fluid space detected for the first time. This is a wedged volume on one direction of the pore walls, shaped thus because the long ligands tilt in order to be connected. We estimate the pressure required for water to become the dominating competing adsorbate within different crystal cavities. The simulated data for CO₂ adsorption in ZIF69 strongly correlate with a set of Raman spectroscopy intensity values that correspond to the same adsorbate-adsorbent system.