Poomiwat Phadungbut

A GCMC simulation and experimental study of krypton adsorption/desorption hysteresis on a graphite surface

Adsorption isotherms and isosteric heats of krypton on a highly graphitized carbon black, Carbopack F, have been studied with a combination of Monte Carlo simulation and high-resolution experiments at 77K and 87K. Our investigation sheds light on the microscopic origin of the experimentally observed, horizontal hysteresis loop in the first layer, and the vertical hysteresis-loop in the second layer, and is found to be in agreement with our recent Monte Carlo simulation study (Diao et al., 2015). From detailed analysis of the adsorption isotherm, the latter is attributed to the compression of an imperfect solid-like state in the first layer, to form a hexagonally packed, solid-like state, immediately following the first order condensation of the second layer. To ensure that capillary condensation in the confined spaces between microcrystallites of Carbopack F does not interfere with these hysteresis loops, we carried out simulations of krypton adsorption in the confined space of a wedge-shaped pore that mimics the interstices between particles. These simulations show that, up to the third layer, any such interference is negligible.

On the phase transition in a monolayer adsorbed on graphite at temperatures below the 2D-critical temperature

Monte Carlo simulations in the grand ensemble and meso-canonical ensemble in which the adsorbent is connected to a finite reservoir have been used to study adsorption isotherms for monolayer argon adsorption on graphite at temperatures below the 2D-critical temperature in order to elucidate the microscopic details of the 2D-transitions: vapour–solid, vapour–liquid and liquid–solid. An S-shaped van der Waals (vdW) loop was found when a small square surface was used; however, for large square surfaces and rectangular surfaces the isotherms exhibit a vdW-type loop with a vertical segment which indicates the coexistence of two phases separated by a boundary that changes its shape with the loading. This coexistence occurs at the same chemical potential as determined by the mid-density scheme, developed by Do and co-workers (Z. Liu, L. Herrera, V.T. Nguyen, D.D. Do, and D. Nicholson, A Monte Carlo scheme based on mid-density in a hysteresis loop to determine equilibrium phase transition. Mol Simul. 37(11):932–939, 2011; Z. Liu, D.D. Do, and D. Nicholson, A thermodynamic study of the mid-density scheme to determine the equilibrium phase transition in cylindrical pores. Mol Simul. 38(3):189–199, 2011).

Adsorption on ordered and disordered duplex layers of porous anodic alumina

We have carried out systematic experiments and numerical simulations of the adsorption on porous anodic aluminum oxide (AAO) duplex layers presenting either an ordered or a disordered interconnecting interface between the large (cavity) and small (constriction) sections of the structured pores. Selective blocking of the pore openings resulted in three different pore topologies: open structured pores, funnel pores, and ink-bottle pores. In the case of the structured pores having an ordered interface, the adsorption isotherms present a rich phenomenology characterized by the presence of two steps in the condensation branch and the opening of one (two) hysteresis loops during evaporation for the ink-bottle (open and funnel) pores. The isotherms can be obtained by summing the isotherms measured on uniform pores having the dimensions of the constrictions or of the cavities. The numerical analysis of the three different pore topologies indicates that the shape of the junction between the two pore sections is only important for the adsorption branch. In particular, a conic junction which resembles that of the AAO pores represents the experimental isotherms for the open and funnel pores better, but the shape of the junction in the ink bottle pores does not matter. The isotherms for the duplex layers with a disordered interface display the same general features found for the ordered duplex layers. In both cases, the adsorption branches coincide and have two steps which are shifted to lower relative pressures compared to those for the ordered duplex. Furthermore, the desorption branches comprise hysteresis loops much wider than those of the ordered duplex layers. Overall, this study highlights the important role played by morphologies where there are interconnections between large and small pores.

Computational methodology for determining textural properties of simulated porous carbons

We have refined and improved the computational efficiency of the TriPOD technique, used to determine the accessible characteristics of porous solids with a known configuration of solid atoms. Instead of placing a probe molecule randomly, as described in the original version of the TriPOD method (Herrera et al., 2011), we implemented a scheme for dividing the porous solid into 3D-grids and computing the solid-fluid potential energies at these grid points. We illustrate the potential of this technique in determining the total pore volume, the surface area and the pore size distribution of various molecular models of porous carbons, ranging from simple pore models to a more complex simulated porous carbon model; the latter is constructed from a canonical Monte Carlo simulation of carbon microcrystallites of various sizes.