Computer modeling and analysis of heterogeneous structures of microporous carbonaceous materials

The aim of this work was to study the problems connected with computer modeling and analysis of heterogeneous structures of microporous carbonaceous materials. The research was focused on the numerical properties of original mathematical models for heterogeneous multilayer adsorption on microporous carbonaceous materials presented in our earlier papers and their applicability to examination of real microporous materials. These models are aimed at drawing information on pore structure and capacity on the basis of adsorption isotherms of small molecule adsorbates. They easily fit typical adsorption data in wide relative pressure ranges. In the theory presented, adsorption of small nearly spherical molecules in irregular pores of molecular size has been considered and side adsorbate–adsorbate interactions are neglected. The molecules mentioned are located in pores by forming aggregates, the size of which is limited by the geometry of the pores. The set of adsorbate molecules, which were adsorbed mainly due to adhesive interactions with the adsorbent matter, is treated as the first layer adsorption. Joining further molecules is viewed as the second, third,... layer adsorption. The main idea of the approach to modeling microporous structure presented, consists of introducing of realistic relationships between geometrical properties of pores and adsorption energy. Special attention was focused on the analysis of the influence of the number of model parameters on identification reliability and evaluation errors of porous structure parameters. This paper gives more information on properties of the identification technique presented in our earlier papers. The five-parameter and six-parameter identification reliability is analyzed in more detail, for different values of the system parameters. In this context, the efficiency of simultaneous examination of two isotherms is also studied.     

Mechanism of methane adsorption on groove space in organic matter surface

Organic matter plays an important role in methane adsorption in shale. Pore surface of organic matter is usually rough and uneven, which results in a large amount of groove space on the pore surface. Thus, the influence of groove space on the adsorption capacity of methane in shale cannot be neglected. Nanoscale pore structures of the organic-rich shale in the Longmaxi Formation were investigated by low-pressure nitrogen gas adsorption as a basis for constructing models. We simplified the internal groove space into triangular prisms with different angles. The grand canonical Monte Carlo simulation and molecular dynamics simulation were used to analyse the methane molecule adsorption behaviour in pores. The results showed that the pore morphology of organic-rich shale in the Longmaxi Formation was mainly slit-shaped pores. The excess adsorption isotherms showed good agreement between experiment and simulation, indicating that the model is suitable and reliable. Methane molecules can enter into the groove space with an opening size of 0.738 nm, while they fail to enter into groove spaces with an opening size less than 0.492 nm. This understanding has important significance for the study of the adsorption characteristics of organic pores which have undergone multiple evolutions in geological history.     

Modelling Gas Adsorption in Slit-Pores Using Monte Carlo Simulation

Abstract

We discuss the use of Monte Carlo simulation to model the equilibrium adsorption of gases in slit pores. Databases of adsorption isotherms have been calculated for nitrogen, carbon-monoxide, methane and carbon-dioxide for a range of pressures, pore widths and temperatures. We discuss the implications of these results for materials characterisation procedures based on gas adsorption data.     

On the existence of a hysteresis loop in open and closed end pores

We have studied the adsorption of argon at 87 K in slit pores of finite length with a smooth graphitic potential, open at both ends or closed at one end. Simulations were carried out using conventional GCMC (grand canonical Monte Carlo) or kMC (kinetic Monte Carlo) in the canonical ensemble with extremely long Markov chain, of at least 2 × 10⁸ configurations; selected simulations with much longer Markov chains do not show any change in the results. When the pore width is in the micropore range (0.65 nm), type I isotherms are obtained for both pore models and for both simulation methods. However, wider pores (1, 2 and 3 nm in width) all exhibit hysteresis loops in the GCMC simulations, while in the canonical ensemble simulations, the isotherms pass through a sigmoid van der Waals type loop in the transition region. This loop locates the true equilibrium transition. For the pores with one closed end, this transition is close to, or coincides with, the adsorption branch of the GCMC hysteresis loop, but for the open-ended pores, it is more closely associated with the desorption branch. In a separate study of adsorption hysteresis in an infinitely long slit pore, using both simulation techniques, the van der Waals loop follows the adsorption branch of the GCMC isotherm to the transition, then reverts to a long vertical section that falls midway between the two hysteresis branches and finally moves to the desorption transition close to the evaporation pressure. An examination of molecular distributions inside the pores reveals two coexisting phases in the canonical simulations, whereas in the grand canonical simulations, the molecules are uniformly distributed along the length of the pores.     

Flow dynamics in bioreactors containing tissue engineering scaffolds

Bioreactors are widely used in tissue engineering as a way to distribute nutrients within porous materials and provide physical stimulus required by many tissues. However, the fluid dynamics within the large porous structure are not well understood. In this study, we explored the effect of reactor geometry by using rectangular and circular reactors with three different inlet and outlet patterns. Geometries were simulated with and without the porous structure using the computational fluid dynamics software Comsol Multiphysics 3.4 and/or ANSYS CFX 11 respectively. Residence time distribution analysis using a step change of a tracer within the reactor revealed non-ideal fluid distribution characteristics within the reactors. The Brinkman equation was used to model the permeability characteristics with in the chitosan porous structure. Pore size was varied from 10 to 200 microm and the number of pores per unit area was varied from 15 to 1,500 pores/mm(2). Effect of cellular growth and tissue remodeling on flow distribution was also assessed by changing the pore size (85-10 microm) while keeping the number of pores per unit area constant. These results showed significant increase in pressure with reduction in pore size, which could limit the fluid flow and nutrient transport. However, measured pressure drop was marginally higher than the simulation results. Maximum shear stress was similar in both reactors and ranged approximately 0.2-0.3 dynes/cm(2). The simulations were validated experimentally using both a rectangular and circular bioreactor, constructed in-house. Porous structures for the experiments were formed using 0.5% chitosan solution freeze-dried at -80 degrees C, and the pressure drop across the reactor was monitored.     

Properties of Confined Square-well Fluids using the Gibbs Ensemble Simulation Technique

In this work we have used the extension of the Gibbs ensemble simulation technique to inhomogeneous fluids [Panagiotopoulos, A.Z. (1987) "Adsorption and capillary condensation of fluid in cylindrical pores by Monte Carlo simulation in the Gibbs ensemble", Mol. Phys. , 62 (3), 701-719], which has been applied to adsorption phenomena of confined fluids. Fluid molecules are described by spherical particles interacting via a square-well potential. The fluid is confined in two types of walls: symmetrical (two hard walls) and non-symmetrical (one square-well wall and one hard wall). In order to analyze the behavior of the confined fluid by varying the potential parameters, we evaluated the bulk and confined densities, the internal energies and the density profiles for different supercritical temperatures. A variety of adsorption profiles can be obtained by using this model. The simulation data reported here complements the available simulation data for this system and can be useful in the development of inhomogeneous fluid theories. Since the square-well parameters can be related to real molecules this system can also be used to understand real adsorption systems.     

Domains and anomalous adsorption isotherms of dipalmitoylphosphatidylcholine membranes and lipophilic ions: pentachlorophenolate, tetraphenylborate, and dipicrylamine.

Dipalmitoylphosphatidylcholine (DPPC) vesicles acquire negative surface charge on adsorption of negatively charged pentachlorophenolate (PCP-), and lipophilic ions tetraphenylborate (TPhB-), and dipicrylamine (DPA-). We have obtained (a) zeta-potential isotherms from the measurements of electrophoretic mobility of DPPC vesicles as a function of concentration of the adsorbing ions at different temperatures (25-42 degrees C), and (b) studied the effect of PCP- on gel-to-fluid phase transition by measuring the temperature dependence of zeta-potential at different PCP- concentrations. The zeta-potential isotherms of PCP- at 25, 32, and 34 degrees C correspond to adsorption to membrane in its gel phase. At 42 degrees C the zeta-potential isotherm corresponds to membrane in its fluid phase. These isotherms are well described by a Langmuir-Stern-Grahame adsorption model proposed by McLaughlin and Harary (1977. Biochemistry. 15:1941-1948). The zeta-potential isotherm at 37 degrees C does not follow the single-phase adsorption model. We have also observed anomalous adsorption isotherms for lipophilic ions TPhB- and DPA- at temperatures as low as 25 degrees C. These isotherms demonstrate a gel-to-fluid phase transition driven by ion adsorption to DPPC membrane during which the membrane changes from weakly to a strongly adsorbing state. The anomalous isotherm of PCP- and the temperature dependence of zeta-potential can be described by a two-phase model based on the combination of (a) Langmuir-Stern-Grahame model for each phase, (b) the coexistence of gel and fluid domains, and (c) depression of gel-to-fluid phase transition temperature by PCP-. Within the anomalous region the magnitude of zeta-potential rapidly increases concentration of adsorbing species, which was characterized in terms of a Esin-Markov coefficient. This effect can be exploited in membrane-based devices. Comments are also made on the possible effect of PCP, as an uncoupler, in energy transducing membranes.     

Pressure effects in confined nanophases

In this article, we review how pressure effects in pores affect both the physics of the confined fluid and the properties of the host porous material. Molecular simulations in which high-pressure effects were observed are first discussed; we will see how the strong dependence on bulk phase pressure of the freezing temperature of a fluid confined in nanopores can be explained by important variations of the pressure within the pore. We then discuss recent works in which direct calculations of the pressure tensor of fluids confined in pores provide evidence for large pressure enhancements. Finally, practical applications of these pressure effects in which gas adsorption in microporous solids (pore size < 2 nm) was found to enhance their mechanical properties by increasing the elastic modulus by a factor 4 are discussed.     

Rate-limiting mass transfer in immunosorbents: characterisation of the adsorption of paraquat-protein conjugates to anti-paraquat Sepharose 4B

A series of paraquat-protein conjugates of different molecular size has been prepared by the coupling of paraquat hexanoate to the proteins lysozyme, ovalbumin, bovine serum albumin. The characteristics of the adsorption of these conjugates to an immunosorbent consisting of monoclonal anti-paraquat antibodies covalently immobilised to Sepharose 4B have been determined. Equilibrium adsorption isotherms were found to obey the Langmuir equation and indicated that 80% or more of the antibody binding sites were accessible to the conjugates. The rates of mass transfer of the conjugates to their adsorption sites on the immobilised antibodies was well described by a model in which mass transfer is controlled by transfer across the external film and diffusion within the porous adsorbent bead. The effective diffusivities of the conjugates within the immunosorbent were measured and has allowed the effect of the size of the adsorbing molecule on the rate of adsorption to be considered. The amount of paraquat that could be adsorbed and the rates of adsorption decreased as the size of the protein to which it is conjugated increased. The diffusivity of the conjugates within the pores of the adsorbent is reduced between two and five times compared to their diffusivities in free solution. The reduction is greater for the larger proteins and the variations of the effective diffusivities and the pore diffusivities with the molecular weight of the conjugate can be well described with simple correlations.     

Adsorption of phenol on formaldehyde-pretreated Pinus pinaster bark: equilibrium and kinetics

This work studies phenol adsorption on Pinus pinaster bark that has been previously treated with formaldehyde in acid medium. The influence of several variables such as solid/liquid ratio, pH and initial concentration of phenol in the solution on the adsorption capacity of the bark has been analysed. A kinetic model based on phenol diffusion within the pores of the adsorbent was in agreement with the results obtained for high initial concentrations of phenol, allowing the determination of diffusion coefficients. Adsorption equilibrium data were fitted by the Freundlich and BET isotherms. From their parameters phenol adsorption capacity and intensity, as well as the specific surface (BET) of the adsorbent, were determined.     

Adsorption isotherms for dilute solutions via the mean force method

The adsorption of solutes from a dilute liquid solution is of great technical importance but calculations of the local density of the solute and of the adsorption isotherm by standard molecular simulation yield large scattering with increasing dilution. As alternative the mean force (MF) method was suggested where the MF on a constrained solute molecule is integrated over a path from the bulk fluid to the wall. It has already been shown that the MF method gives reliable results for the relative local density, even at high dilution. Here, an extension of this method is introduced, where the absolute value of the bulk density is determined by particle balance. Thus, it is possible to calculate adsorption isotherms from the Henry regime to any finite concentration. Molecular dynamics simulations for the local density and the adsorption isotherm were performed for a model solution consisting of tetrahedral Lennard–Jones (LJ) solvent and linear LJ solute molecules in contact with a plane wall. It is found that the MF-results show less scattering than the results from standard simulations. Moreover, results for the orientation and the selectivity are given.     

Pore network modelling of affinity chromatography: determination of the dynamic profiles of the pore diffusivity of beta-galactosidase and its effect on column performance as the loading of beta-galactosidase onto anti-beta-galactosidase varies with time.

A three-dimensional pore network model for diffusion in porous adsorbent particles was employed in a dynamic adsorption model that simulates the adsorption of a solute in porous particles packed in a chromatographic column. The solution of the combined model yielded the dynamic profiles of the pore diffusion coefficient of beta-galactosidase along the radius of porous adsorbent particles and along the length of the column as the loading of beta-galactosidase onto anti-beta-galactosidase immobilized on the surface of the pores of the particles occurred, and, the dynamic adsorptive capacity of the chromatographic column as a function of the design and operational parameters of the chromatographic system. It was found that for a given column length the dynamic profiles of the pore diffusion coefficient were influenced by (a) the superficial fluid velocity in the column, (b) the diameter of the adsorbent particles, and (c) the pore connectivity of the porous structure of the adsorbent particles. The effect of the magnitude of the pore connectivity on the dynamic profiles of the pore diffusion coefficient of beta-galactosidase increased as the diameter of the adsorbent particles and the superficial fluid velocity in the column increased. The dynamic adsorptive capacity of the column increased as (i) the particle diameter and the superficial fluid velocity in the column decreased, and (ii) the column length and the pore connectivity increased. In preparative affinity chromatography, it is desirable to obtain high throughputs within acceptable pressure gradients, and this may require the employment of larger diameter adsorbent particles. In such a case, longer column lengths satisfying acceptable pressure gradients with adsorbent particles having higher pore connectivity values could provide high dynamic adsorptive capacities. An alternative chromatographic system could be comprised of a long column packed with large particles which have fractal pores (fractal particles) that have high pore connectivities and which allow high intraparticle diffusional and convective flow mass transfer rates providing high throughputs and high dynamic adsorptive capacities. If large scale monoliths could be made to be reproducible and operationally stable, they could also offer an alternative mode of operation that could provide high throughputs and high dynamic adsorptive capacities.     

Differences in the adsorption and diffusion behaviour of water and non-polar gases in nanoporous carbon: role of cooperative effects of pore confinement and hydrogen bonding

We investigate the effect of pore confinement and molecular geometry on the adsorption and self-diffusion of H₂O, CO₂, Ar, CH₄, C₃H₆, SF₆ and C₅H₁₂, in a realistic model of nanoporous silicon carbide derived carbon (SiC-DC), constructed using hybrid reverse Monte Carlo simulation. Adsorption isotherms, adsorbate–adsorbate and adsorbate–adsorbent contributions to the isosteric heat of adsorption are determined to study the effect of pore confinement, microporosity and molecular geometry on adsorption of these gases. We describe the cooperative effect of pore confinement and hydrogen bonding on the formation of water clusters and anomalous adsorption behaviour of water compared with non-polar gases. We find that, in contrast to literature results based on the slit-pore model, pore-filling does not occur below the saturation pressure in hydrophobic amorphous carbon materials such as SiC-DC and activated carbon fibre. We also compare self-diffusivities and activation energy barriers of water and non-polar gases in the microporous structure of SiC-DC to identify underlying correlations with molecular properties. We demonstrate that the self-diffusivity of water deviates considerably from the correlation between diffusivity and molecular kinetic diameter observed for non-polar gases. This is attributed to the reduced diffusivity of water, and its relatively large energy barrier at high loadings despite its small kinetic diameter, which is due to the blocking effect of water clusters at pore entries.     

A grand canonical Monte Carlo simulation study of argon and krypton confined inside weakly attractive slit pores

Grand canonical Monte Carlo simulations are performed to investigate the adsorption of argon and krypton inside weakly attractive slit pores. We examine the effects of confinement on these monoatomic fluids (modelled using the triangle-well potential) in a hard wall slit pore as also when the pore-fluid interactions are uniformly and weakly attractive. The effects of temperature and pressure on the adsorption isotherms of these confined fluids are found to be the same as those reported in literature. The equilibrium density profiles for argon and krypton exhibit both uniform distribution and layering under different conditions. In addition, for krypton, under specific conditions inside the narrow pores, we note the development of frustrated layering.     

Predicting helium and neon adsorption and separation on carbon nanotubes by Monte Carlo simulation

The adsorption of helium and neon mixtures on single-walled carbon nanotubes (SWCNTs) was investigated at various temperatures (subcritical and supercritical) and pressures using canonical Monte Carlo (CMC) simulation. Adsorption isotherms were obtained at different temperatures (4, 40, 77 and 130 K) and pressures ranging from 1 to 16 MPa. Separation factors and isosteric enthalpies of adsorption were also calculated. Moreover, the adsorption isotherms were obtained at constant specific temperatures (4 and 40 K) and pressures (0.2 and 1.0 MPa) as a function of the amount adsorbed. All of the adsorption isotherms for an equimolar mixture of helium and neon have a Langmuir shape, indicating that no capillary condensation occurs. Both the helium and the neon adsorption isotherms exhibit similar behavior, and slightly more of the helium and neon mixture is adsorbed on the inner surfaces of the SWCNTs than on their outer surfaces. More neon is adsorbed than helium within the specified pressure range. The data obtained show that the isosteric enthalpies for the adsorption of neon are higher than those for helium under the same conditions, which means that adsorption of neon preferentially occurs by (15, 15) SWCNTs. Furthermore, the isosteric enthalpies of adsorption of both gases decrease with increasing temperature.     

Local mobility in lipid domains of supported bilayers characterized by atomic force microscopy and fluorescence correlation spectroscopy

Fluorescence correlation spectroscopy (FCS) is used to examine mobility of labeled probes at specific sites in supported bilayers consisting of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid domains in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). Those sites are mapped beforehand with simultaneous atomic force microscopy and submicron confocal fluorescence imaging, allowing characterization of probe partitioning between gel DPPC and disordered liquid DOPC domains with corresponding topography of domain structure. We thus examine the relative partitioning and mobility in gel and disordered liquid phases for headgroup- and tailgroup-labeled GM1 ganglioside probes and for headgroup- and tailgroup-labeled phospholipid probes. For the GM1 probes, large differences in mobility between fluid and gel domains are observed; whereas unexpected mobility is observed in submicron gel domains for the phospholipid probes. We attribute the latter to domain heterogeneities that could be induced by the probe. Furthermore, fits to the FCS data for the phospholipid probes in the DOPC fluid phase require two components (fast and slow). Although proximity to the glass substrate may be a factor, local distortion of the probe by the fluorophore could also be important. Overall, we observe nonideal aspects of phospholipid probe mobility and partitioning that may not be restricted to supported bilayers.     

Pneumolysin-damaged cells benefit from non-homogeneous toxin binding to cholesterol-rich membrane domains

Nucleated cells eliminate lesions induced by bacterial pore-forming toxins, such as pneumolysin via shedding patches of damaged plasmalemma into the extracellular milieu. Recently, we have shown that the majority of shed pneumolysin is present in the form of inactive pre-pores. This finding is surprising considering that shedding is triggered by Ca²⁺-influx following membrane perforation and therefore is expected to positively discriminate for active pores versus inactive pre-pores.Here we provide evidence for the existence of plasmalemmal domains that are able to attract pneumolysin at high local concentrations. Within such a domain an immediate plasmalemmal perforation induced by a small number of pneumolysin pores would be capable of triggering the elimination of a large number of not yet active pre-pores/monomers and thus pre-empt more frequent and perilous perforation events. Our findings provide further insights into the functioning of the cellular repair machinery which benefits from an inhomogeneous plasmalemmal distribution of pneumolysin.