Kinetics of adsorption on activated carbon: application of heterogeneous vacancy solution theory

The kinetics of single component adsorption on activated carbon is investigated here using a heterogeneous vacancy solution theory (VST) of adsorption. The adsorption isotherm is developed to account for the adsorbate non-ideality due to the size difference between the adsorbate molecule and the vacant site, while incorporating adsorbent heterogeneity through a pore-width-related potential energy. The transport process in the bidisperse carbon considers coupled mass transfer in both macropore and micropore phases simultaneously. Adsorbate diffusion in the micropore network is modeled through effective medium theory, thus considering pore network connectivity in the adsorbent, with the activation energy for adsorbate diffusion related to the adsorption energy, represented by the Steele 10-4-3 potential for carbons. Experimental data of five hydrocarbons, CO2 and SO2 on Ajax carbon at multiple temperatures, as well as three hydrocarbons on Norit carbon at three temperatures are first fitted by the heterogeneous VST model to obtain the isotherm parameters, followed by application of the kinetic model to uptake data on carbon particles of different sizes and geometry at various temperatures. For the hydrocarbons studied, the model can successfully correlate the experimental data for both adsorption equilibrium and kinetics. However, there is some deviation in the fit of the desorption kinetics for polar compounds such as CO2 and SO2, due to the inadequacy of the L-J potential model in this case. The significance of viscous transport in the micropores is also considered here and found to be negligible, consistent with recent molecular simulation studies. (C) 2002 Elsevier Science Ltd. All rights reserved.

Vacancy solution theory for binary adsorption equilibria in heterogeneous carbon

A heterogeneous modified vacancy solution model of adsorption developed is evaluated. The new model considers the adsorption process through a mass-action law and is thermodynamically consistent, while maintaining the simplicity in calculation of multicomponent adsorption equilibria, as in the original vacancy solution theory. It incorporates the adsorbent heterogeneity through a pore-width-related potential energy, represented by Steele's 10-4-3 potential expression. The experimental data of various hydrocarbons, CO2 and SO2 on four different activated carbons - Ajax, Norit, Nuxit, and BPL - at multiple temperatures over a wide range of pressures were studied by the heterogeneous modified VST model to obtain the isotherm parameters and micropore-size distribution of carbons. The model successfully correlates the single-component adsorption equilibrium data for all compounds studied on various carbons. The fitting results for the vacancy occupancy parameter are consistent with the pressure change on different carbons, and the effect of pore heterogeneity is important in adsorption at elevated pressure. It predicts binary adsorption equilibria better than the IAST scheme, reflecting the significance of molecular size nonideality.

Density functional theory analysis of the influence of pore wall heterogeneity on adsorption in carbons

Density functional theory for adsorption in carbons is adapted here to incorporate a random distribution of pore wall thickness in the solid, and it is shown that the mean pore wall thickness is intimately related to the pore size distribution characteristics. For typical carbons the pore walls are estimated to comprise only about two graphene layers, and application of the modified density functional theory approach shows that the commonly used assumption of infinitely thick walls can severely affect the results for adsorption in small pores under both supercritical and subcritical conditions. Under supercritical conditions the Henry's law coefficient is overpredicted by as much as a factor of 2, while under subcritical conditions pore wall heterogeneity appears to modify transitions in small pores into a sequence of smaller ones corresponding to pores with different wall thicknesses. The results suggest the need to improve current pore size distrubution analysis methods to allow for pore wall heterogeneity. The density functional theory is further extended here to allow for interpore adsorbate interactions, and it appears that these interaction are negligible for small molecules such as nitrogen but significant for more strongly interacting heavier molecules such as butane, for which the traditional independent pore model may not be adequate.

Cold Storage Preservation and Warm Ischaemic Injury to Isolated Arterial Segments: Endothelial Cell Injury

Injury to endothelial calls is thought to be important to the development of the vascular lesion of chronic rejection. It was the aim of this study to develop a semiquantitative method to assess endothelial injury in arterial grafts and to document the injury produced by cold storage preservation and additional warm ischaemia. Twelve- and 24-h cold preservation of rat aortic segments, together with an additional 1 h of warm ischaemia, were assessed. Electron micrographs of representative endothelial cells were scored for cytoplasmic, nuclear and mitochondrial injury. The overall injury score was obtained by addition of the individual scores. Storage for up to 24 h in University of Wisconsin (UW) and Terasaki did not produce any injury. Twenty-four hours of storage in Euro-Collins resulted in endothelial cell death. Injury occurred after 12 h of storage in Ross, Collins and normal saline, and the injury increased following 24 h of storage. One hour of warm ischaemia did not increase the injury. Injury to endothelial cells varies with the preservation solution used and the time of cold storage, so that both the type of solution and the storage time should be taken into account in clinical studies looking at the influence of cold ischaemia time and graft outcome.

Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

Simulation of the load-unload response ratio and critical sensitivity in the lattice solid model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Solution Structures of the cis- and trans-Pro30 Isomers of a Novel 38-Residue Toxin from the Venom of Hadronyche Infensa sp. that Contains a Cystine-Knot Motif within Its Four Disulfide Bonds

The primary sequence and three-dimensional structure of a novel peptide toxin isolated from the Australian funnel-web spider Hadronyche infensa sp. is reported. ACTX-HI:OB4219 contains 38 amino acids, including eight-cysteine residues that form four disulfide bonds. The connectivities of these disulfide bonds were previously unknown but have been unambiguously determined in this study. Three of these disulfide bonds are arranged in an inhibitor cystine-knot (ICK) motif, which is observed in a range of other disulfide-rich peptide toxins. The motif incorporates an embedded ring in the structure formed by two of the disulfides and their connecting backbone segments penetrated by a third disulfide bond. Using NMR spectroscopy, we determined that despite the isolation of a single native homologous product by RP-HPLC, ACTX-HI:OB4219 possesses two equally populated conformers in solution. These two conformers were determined to arise from cis/trans isomerization of the bond preceding Pro30. Full assignment of the NMR spectra for both conformers allowed for the calculation of their structures, revealing, the presence of a triple-stranded antiparallel sheet consistent with the inhibitor cystine-knot (ICK) motif.

The Three-dimensional Solution Structure of NaD1, a New Floral Defensin from Nicotiana alata and its Application to a Homology Model of the Crop Defense Protein alfAFP

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of NaD1, a novel antifungal and insecticidal protein isolated from the flowers of Nicotiana alata. NaD1 is a basic, cysteine-rich protein of 47 residues and is the first example of a plant defensin from flowers to be characterized structurally. Its three-dimensional structure consists of an a-helix and a triple-stranded anti-parallel beta-sheet that are stabilized by four intramolecular disulfide bonds. NaD1 features all the characteristics of the cysteine-stabilized up motif that has been described for a variety of proteins of differing functions ranging from antibacterial insect defensins and ion channel-perturbing scorpion toxins to an elicitor of the sweet taste response. The protein is biologically active against insect pests, which makes it a potential candidate for use in crop protection. NaD1 shares 31% sequence identity with alfAFP, an antifungal protein from alfalfa that confers resistance to a fungal pathogen in transgenic potatoes. The structure of NaD1 was used to obtain a homology model of alfAFP, since NaD1 has the highest level of sequence identity with alfAFP of any structurally characterized antifungal defensin. The structures of NaD1 and alfAFP were used in conjunction with structure - activity data for the radish defensin Rs-AFP2 to provide an insight into structure-function relationships. In particular, a putative effector site was identified in the structure of NaD1 and in the corresponding homology model of alfAFP. (C) 2002 Elsevier Science Ltd. All rights reserved.

Synthesis and structural analysis of the N-terminal domain of the thyroid hormone-binding protein transthyretin

Transthyretin (TTR) is a 55 kDa protein responsible for the transport of thyroid hormones and retinol in human serum. Misfolded forms of the protein are implicated in the amyloid diseases familial amyloidotic polyneuropathy and senile systemic amyloidosis. Its folding properties and stabilization by ligands are of current interest due to their importance in understanding and combating these diseases. To assist in such studies we developed a method for the solid phase synthesis of the monomeric unit of a TTR analogue and its folding to form a functional 55 kDa tetramer. The monomeric unit of the protein was chemically synthesized in three parts, comprising amino acid residues 151, 5499 and 102127, and ligated using chemoselective thioether ligation chemistry. The synthetic protein was folded and assembled to a tetrameric structure in the presence of the TTRs native ligand, thyroxine, as shown by gel filtration chromatography, native gel electrophoresis, TTR antibody recognition and thyroid hormone binding. In the current study the solution structure of the first of these fragment peptides, TTR(151) is examined to determine its intrinsic propensity to form beta-sheet structure, potentially involved in amyloid fibril formation by TTR. Despite the presence of extensive beta-structure in the native form of the protein, the Nterminal fragment adopts an essentially random coil conformation in solution.

Effect of soil composition and dissolved organic matter on pesticide sorption

The effect of the solid and dissolved organic matter fractions, mineral composition and ionic strength of the soil solution on the sorption behaviour of pesticides were studied. A number of soils, chosen so as to have different clay mineral and organic carbon content, were used to study the sorption of the pesticides atrazine (6-chloro-N-2-ethyl-N-4-isopropyl-1,3,5-triazine-2,4-diamine), 2,4-D ((2,4-dichlorophenoxy) acetic acid), isoproturon (3-(4-isopropylphenyl)1,1-dimethylurea) and paraquat (1,1'-dimethyl-4,4'-bipyridinium) in the presence of low and high levels of dissolved organic carbon and different background electrolytes. The sorption behaviour of atrazine, isoproturon and paraquat was dominated by the solid state soil components and the presence of dissolved organic matter had little effect. The sorption of 2,4-D was slightly affected by the soluble organic matter in the soil. However, this effect may be due to competition for adsorption sites between the pesticide and the soluble organic matter rather than due to a positive interaction between the pesticide and the soluble fraction of soil organic matter. It is concluded that the major factor governing the sorption of these pesticides is the solid state organic fraction with the clay mineral content also making a significant contribution. The dissolved organic carbon fraction of the total organic carbon in the soil and the ionic strength of the soil solution appear to have little or no effect on the sorption/transport characteristics of these pesticides over the range of concentrations studied. (C) 2002 Elsevier Science B.V. All rights reserved.

Electron exchange coupling for single-donor solid-state spin qubits

Intervalley interference between degenerate conduction band minima has been shown to lead to oscillations in the exchange energy between neighboring phosphorus donor electron states in silicon [B. Koiller, X. Hu, and S. Das Sarma, Phys. Rev. Lett. 88, 027903 (2002); Phys. Rev. B 66, 115201 (2002)]. These same effects lead to an extreme sensitivity of the exchange energy on the relative orientation of the donor atoms, an issue of crucial importance in the construction of silicon-based spin quantum computers. In this article we calculate the donor electron exchange coupling as a function of donor position incorporating the full Bloch structure of the Kohn-Luttinger electron wave functions. It is found that due to the rapidly oscillating nature of the terms they produce, the periodic part of the Bloch functions can be safely ignored in the Heitler-London integrals as was done by Koiller, Hu, and Das Sarma, significantly reducing the complexity of calculations. We address issues of fabrication and calculate the expected exchange coupling between neighboring donors that have been implanted into the silicon substrate using an 15 keV ion beam in the so-called top down fabrication scheme for a Kane solid-state quantum computer. In addition, we calculate the exchange coupling as a function of the voltage bias on control gates used to manipulate the electron wave functions and implement quantum logic operations in the Kane proposal, and find that these gate biases can be used to both increase and decrease the magnitude of the exchange coupling between neighboring donor electrons. The zero-bias results reconfirm those previously obtained by Koiller, Hu, and Das Sarma.

Composition and quality of harvest residues and soil organic matter under windrow residue management in young hoop pine plantations as revealed by solid-state C-13 NMR spectroscopy

Solid-state C-13 nuclear magnetic resonance (NMR) with cross-polarisation (CP) and magic-angle-spinning (MAS) was used to: (a) examine the changes in carbon (C) composition of windrowed harvest residues during the first 3 years of hoop pine plantations in subtropical Australia; (b) assess the impacts of windrowed harvest residues on soil organic matter (SOM) composition and quality in the 0-10 cm soil layer. Harvest residues were collected from 0-, 1-, 2- and 3-year-old windrows of ca. 2.5 m width (15 m apart for 0-, 1- and 2-year-old sites and 10 m apart for 3-year-old site). Soils from the 0 to 10 cm soil layer were collected from the 1-, 2- and 3-year-old sites. The 13C NMR spectra of the harvest residues indicated the presence of lignin in the hoop pine wood, foliage and newly incorporated organic matter (NIOM). Condensed tannin structures were found in the decay-resistant bark, small wood and foliage, but were absent in other residue components and SOM. The NMR spectra of small wood samples contained condensed tannin structures because the outer layer of bark was not removed. NIOM showed a shift from foliage-like structures (celluloses) to lignin-type structures, indicating an incorporation of woody residues from the decomposing harvest residues. Suberins were also present in the small wood, foliage and bark. The 13C CP NMR spectra of SOM indicated that in areas where windrows were present, SOM did not show compositional changes. However, an increase in SOM quality under the windrows in the second year after their formation as characterised by the alkyl C/O-alkyl C (A/O-A) ratio was mainly due to inputs from the decomposition of the labile, readily available components of the windrowed harvest residues. (C) 2002 Published by Elsevier Science B.V.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.