Testing predictions of macroscopic binary diffusion coefficients using lattice models with site heterogeneity

Quantitatively predicting mass transport rates for chemical mixtures in porous materials is important in applications of materials such as adsorbents, membranes, and catalysts. Because directly assessing mixture transport experimentally is challenging, theoretical models that can predict mixture diffusion coefficients using Only single-component information would have many uses. One such model was proposed by Skoulidas, Sholl, and Krishna (Langmuir, 2003, 19, 7977), and applications of this model to a variety of chemical mixtures in nanoporous materials have yielded promising results. In this paper, the accuracy of this model for predicting mixture diffusion coefficients in materials that exhibit a heterogeneous distribution of local binding energies is examined. To examine this issue, single-component and binary mixture diffusion coefficients are computed using kinetic Monte Carlo for a two-dimensional lattice model over a wide range of lattice occupancies and compositions. The approach suggested by Skoulidas, Sholl, and Krishna is found to be accurate in situations where the spatial distribution of binding site energies is relatively homogeneous, but is considerably less accurate for strongly heterogeneous energy distributions.

A diffusion and T2 relaxation MRI study of the ovine lumbar intervertebral disc under compression in vitro

The ovine lumbar intervertebral disc is a useful model for the human lumbar disc. We present preliminary estimates of diffusion coefficients and T-2 relaxation times in a pilot MRI study of the ovine lumbar intervertebral disc during uniaxial compression in vitro, and identify factors that hamper the ability to accurately monitor the temporal evolution of the effective diffusion tensor at high spatial resolution.

NMR imaging of the diffusion of water at 310 K into semi-IPNs of PEM and Poly(HEMA-co-THFMA)

Magnetic resonance imaging has been used to monitor the diffusion of water at 310 K into a series of semi-IPNs of poly(ethyl methacrylate), PEM, and copolymers of 2-hydroxyethyl methacrylate, HEMA, and tetrahydrofurfuryl methacrylate, THFMA. The diffusion was found to be well described by a Fickian kinetic model in the early stages of the water sorption process, and the diffusion coefficients were found to be slightly smaller than those for the copolymers of HEMA and THFMA, P(HEMA-co-THFMA), containing the same mole fraction of HEMA in the matrix. A second stage sorption process was identified in the later stage of water sorption by the PEM/PTHFMA semi-IPN and for the systems containing a P(HEMA-co-THFMA) component with a mole fraction HEMA of 0.6 or less. This was characterized by the presence of Water near the surface of the cylinders with a longer NMR T-2 relaxation time, which would be characteristic of mobile water, such as water present in large pores or surface fissures. The presence of the drug chlorhexidine in the polymer matrixes at a concentration of 5.625 wt % was found not to modify the properties significantly, but the diffusion coefficients for the water sorption were systematically smaller when the drug was present.

NMR imaging of the diffusion of water at 37 degrees C into poly(2-hydroxyethyl methacrylate) containing aspirin or vitamin B-12

The ingress of water into poly(2-hydroxyethyl methacrylate), PHEMA, loaded with either one of two model drugs, vitamin B-12 or aspirin, was studied at 37 degreesC using three-dimensional NMR imaging. PHEMA was loaded with 5 and 10 wt % of the drugs. From the imaging profiles, it was observed that incorporation of vitamin B-12 into PHEMA resulted in enhanced crack formation on sorption of water and the crack healing behind the diffusion front was slower than for PHEMA without added drug. This was accounted for by the anti-plasticization of PHEMA by vitamin B-12. Crack formation was inhibited in the P-HEMA-aspirin systems because of the plasticizing effect of the aspirin on the PHEMA matrix. All of the polymers were found to absorb water according to an underlying Fickian diffusion mechanism. For PHEMA loaded with 5 wt % of aspirin or vitamin B-12, the best values of the water diffusion coefficients were both found to be 1.3 +/- 0.1 x 10(-11) m(2) s(-1) at 37 degreesC, while the values for the polymer loaded with 10 wt % of the drugs were slightly higher, 1.5 +/- 0.1 x 10(-11) m(2) s(-1).

Water diffusion in methacrylate based copolymer hydrogels of 2-hydroxyethyl methacrylate

The diffusion of water into cylinders of polyHEMA and copolymers of HEMA with THFMA, BMA and CHMA were studied over a range of copolymer compositions. The diffusion of water into the polymers was found to follow a Fickian, or case I mechanism. The diffusion coefficients of water were determined from mass measurements and NMR imaging studies. They were found to vary from 1.7 +/- 0.2 x 10(-11) m(2) s(-1) for polyHEMA at 37 degreesC to lower values for the copolymers. The mass of water absorbed at equilibrium relative to the mass of dry polymer varied from 52-58 wt% for polyHEMA to lower values for the copolymers.

Diffusion of linear paraffins in nanoporous silica

The diffusion of hexane, heptane, octane, and decane in nanoporous MCM-41 silica at various temperatures is investigated by the zero-length-column method. The diffusion coefficients are derived by a complete-time-range analysis of desorption curves at different purge flow rates and temperatures. The results show that the calculated low-coverage diffusivity values decrease monotonically, and the derived Henry's law constants increase, as the carbon number of paraffins increases. The study reveals that transport is strongly influenced by intracrystalline diffusion and dominated by the sorbate-sorbent interaction. The diffusion activation energy and adsorption isosteric heat at zero loading increase monotonically with the carbon number of linear paraffins, but their ratio is essentially constant for each adsorbate compound.

Mass uptake study of the diffusion of water and SBF into poly(2-hydroxyethyl methacrylate-co-tetrahydrofurfuryI methacrylate) containing aspirin or vitamin B-12

The ingress of water and Kokubo simulated body fluid (SBF) into poly (2-hydroxyethyl methacrylate) (PHEMA), and its co-polymers with tetrahydrofurduryl methacrylate (THFMA), loaded with either one of two model drugs, vitamin 1312 or aspirin, was studied by mass uptake over the temperature range 298-318 K. The polymers were studied as cylinders and were loaded with either 5 wt% or 10 wt% of the drugs. From DSC studies it was observed that vitamin B-12 behaved as a physical cross-linker restricting chain segmental mobility, and so had a small anti-plasticisation effect on PHEMA and the co-polymers rich in HEMA, but almost no effect on the T-g of co-polymers rich in THFMA. On the other hand, aspirin exhibited a plasticising effect on PHEMA and the copolymers. All of the polymers were found to absorb water and SBF according to a Fickian diffusion mechanism. The polymers were all found to swell to a greater extent in SBF than in water, which was attributed to the presence of Tris buffer in the SBF. The sorptions of the two penetrants were found to follow Fickian kinetics in all cases and the diffusion coefficients at 310 K for SBF were found to be smaller than those for water, except for the polymers containing aspirin where the diffusion coefficients were higher than for the other systems. For example, for sorption into PHEMA the diffusion coefficient for water was 1.41 X 10(-11) m(2)/s and for SBF was 0.79 x 10-11 m(2)/s, but in the presence of 5 wt% aspirin the corresponding values were 1.27 x 10(-1)1 m(2)/s and 1.25 x 10(-11) m(2)/s, respectively. The corresponding values for PHEMA loaded with 5 wt% B-12 were 1.25 x 10(-11) m(2)/s and 0.74 x 10(-11) m(2)/s, respectively.

NMR imaging of the diffusion of water at 310 K into poly[(2-hydroxyethyl methacrylate)-co-(tetrahydrofurfuryl methacrylate)] containing vitamin B12 or aspirin

The ingress of water into copolymers of 2-hydroxyethyl methacrylate (HEMA) and tetrahydrofurfuryl methacrylate (THFMA) loaded with either one of two model drugs, ie vitamin B-12 or aspirin, was studied at 310 K using three-dimensional nuclear magnetic resonance (3D NMR) imaging. The poly(HEMA) was loaded with 5 wt% of the drugs. From the imaging profiles it was observed that incorporation of vitamin B-12 into the polymers rich in HEMA resulted in crack formation at the interface between the rubbery region and the glassy core on sorption of water, although these cracks were 'healed' behind the diffusion front. However, for the copolymers with low HEMA contents and for those containing aspirin, no evidence for similar crack formation was found. For the copolymers loaded with 5 wt% of aspirin or vitamin B-12 the values of the water diffusion coefficients, determined by curve-fitting the relative water concentration profiles from magnetic resonance imaging (MRI) measurements, were found to be smaller than those obtained from a mass uptake study. (C) 2004 Society of Chemical Industry.

Transport diffusion of gases is rapid in flexible carbon nanotu

Molecular dynamics simulations of rigid, defect-free single-walled carbon nanotubes have previously suggested that the transport diffusivity of gases adsorbed in these materials can be orders of magnitude higher than any other nanoporous material (A. I. Skoulidas et al., Phys. Rev. Lett. 2002, 89, 185901). These simulations must overestimate the molecular diffusion coefficients because they neglect energy exhange between the diffusing molecules and the nanotube. Recently, Jakobtorweihen et al. have reported careful simulations of molecular self-diffusion that allow nanotube flexibility (Phys. Rev. Lett. 2005, 95, 044501). We have used the efficient thermostat developed by Jakobtorweihen et al. to examine the influence of nanotube flexibility on the transport diffusion of CH4 in (20,0) and (15,0) nanotubes. The inclusion of nanotube flexibility reduces the transport diffusion relative to the rigid nanotube by roughly an order of magnitude close to zero pressure, but at pressures above about I bar the transport diffusivities for flexible and rigid nanotubes are very similar, differing by less than a factor or two on average. Hence, the transport diffusivities are still extremely large compared to other known materials when flexibility is taken into account.

NMR microscopy: A tool for measuring monomer diffusion in supercritical CO2

The diffusion of styrene into linear low density polyethylene in a solution of supercritical CO2 was investigated using NMR microimaging. For both pure styrene and styrene dissolved in supercritical CO2, the diffusion was found to follow Fickian kinetics. Supercritical CO2 was found to enhance the rate and extent of diffusion of styrene into the substrate by up to three times under the conditions of this investigation, compared to pure styrene. NMR imaging was used to measure the concentration profiles of the styrene penetrants in real time, and the results were fitted to a Fickian model for diffusion. At a CO2 pressure of 150 bar and temperature of 40 degrees C, the diffusion coefficient of a 30 wt-% solution of styrene into LLDPE was calculated to be 1 X 10(-11) m(2). s(-1). This is significantly faster than the diffusion coefficient measured for pure styrene diffusion at 40 degrees C (3 x 10(-12) m(2). s(-1)). The diffusion coefficients determined by gravimetric analysis were slightly higher than those determined by the imaging method. This was probably due to residual styrene and/or polystyrene adhering to the surface of the films in the gravimetric technique.

Estimation of concentration-dependent diffusion coefficient in drying process from the space-averaged concentration versus time with experimental data

The estimation of a concentration-dependent diffusion coefficient in a drying process is known as an inverse coefficient problem. The solution is sought wherein the space-average concentration is known as function of time (mass loss monitoring). The problem is stated as the minimization of a functional and gradient-based algorithms are used to solve it. Many numerical and experimental examples that demonstrate the effectiveness of the proposed approach are presented. Thin slab drying was carried out in an isothermal drying chamber built in our laboratory. The diffusion coefficients of fructose obtained with the present method are compared with existing literature results.

Differential diffusion: Often a finite-mising length effect

Many instances of differential diffusion, i e, different species having different turbulent diffusion coefficients in the same flow, can be explained as a finite mixing length effect. That is, in a simple mixing length scenario, the turbulent diffusion coefficient has the form 1 ( m )2 m m c l K w l OL   =  +    where, wm is the mixing velocity, lm the mixing length and Lc the overall distribution scale for a particular species. The first term represents the familiar gradient diffusion while the second term becomes important when lm/Lc is finite. This second term shows that different species will have different diffusion coefficients if they have different overall distribution scales. Such different Lcs may come about due to different boundary conditions and different intrinsic properties (molecular diffusivity, settling velocity etc) for different species. For momentum transfer in turbulent oscillatory boundary layers the second term is imaginary and explains observed phase leads of shear stresses ahead of velocity gradients.

Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation. (C) 2004 American Institute of Physics.

Heterogeneous model for gas transport in carbon molecular sieves

A dual resistance model with distribution of either barrier or pore diffusional activation energy is proposed in this work for gas transport in carbon molecular sieve (CMS) micropores. This is a novel approach in which the equilibrium is homogeneous, but the kinetics is heterogeneous. The model seems to provide a possible explanation for the concentration dependence of the thermodynamically corrected barrier and pore diffusion coefficients observed in previous studies from this laboratory on gas diffusion in CMS.(1.2) The energy distribution is assumed to follow the gamma distribution function. It is shown that the energy distribution model can fully capture the behavior described by the empirical model established in earlier studies to account for the concentration dependence of thermodynamically corrected barrier and pore diffusion coefficients. A methodology is proposed for extracting energy distribution parameters, and it is further shown that the extracted energy distribution parameters can effectively predict integral uptake and column breakthrough profiles over a wide range of operating pressures.

Weathering geochronology by (U-Th)/He dating of goethite

Nine samples of supergene goethite (FeOOH) from Brazil and Australia were selected to rest the suitability of this mineral for (U-Th)/He dating. Measured He ages ranged from 61 to 8 Ma and were reproducible to better than a few percent despite very large Variations in [U] and [Th]. In all Samples with internal stratigraphy or independent age constraints, the He ages corroborated the expected relationship's. These data demonstrate that internally consistent He ages can be obtained on goethite. but do not prove quantitative 4 He retention. To assess possible diffusive He loss, stepped-heating experiments were performed on two goethite samples that were subjected to proton irradiation to produce a homogeneous distribution of spallogenic He-3. The He-3 release pattern indicates the presence of at least two diffusion domains, one with high helium retentivity and the other with very low retentivity at Earth surface conditions. The low retentivity domain, which accounts for similar to 5% of He-3, contains no natural He-4 and may represent poorly crystalline or intergranular material which has lost all radiogenic He-4 by diffusion in nature. Diffusive loss of He-3 from the high retentivity domain is independent of the macroscopic dimensions of the analyzed polycrystalline aggregate, so probably represents diffusion from individual micrometer-size goethite crystals. The He-2/He-3 evolution during the incremental heating experiments shows that the high retentivity domain has retained 90%-95% of its radiogenic helium. This degree of retentivity is in excellent agreement with that independently predicted from the helium diffusion coefficients extrapolated to Earth surface temperature and held for the appropriate duration. Considering both the high and low retentivity domains, these data indicate that one of the samples retained 90% of its radiogenic He-4 over 47.5 Ma and the other retained 86% over 12.3 Ma. Thus while diffusive-loss corrections to supergene goethite He ages are required. these initial results indicate that the corrections are not extremely large and can be rigorously quantified using the proton-irradiation He-4/He-3 method. Copyright (C) 2005 Elsevier Ltd.