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.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

N-aminopyrroledione-hydrazonoketene-pyrazolium oxide-pyrazolone rearrangements and pyrazolone tautomerism

Flash vacuum thermolysis (FVT) of 1-(dimethylamino)pyrrole-2,3-diones 5 causes extrusion of CO with formation of transient hydrazonoketenes 7. The transient ketenes 7 are observable in the form of weak bands at 2130 (7a) or 2115 cm(-1) (7b) in the Ar matrix IR spectra resulting from either FVT or photolysis of either 5 or 1,1- dimethylpyrazolium-5- oxides 8, and these absorptions are in excellent agreement with B3LYP/6-31G* frequency calculations. Under FVT conditions the ketenes 7 cyclize to pyrazolium oxides 8, which undergo 1,4-migration of a methyl group to yield 1,4-dimethyl-3-phenylpyrazole-5(4H)-one 9a and 1,4,4-trimethyl-3-phenylpyrazole-5(4H)-one 9b. All three tautomers of 9a have been characterized, viz. the CH form 9a (most stable form in the gas phase, the solid state and solvents of low polarity), the OH form 9a' (metastable solid at room temperature) and the NH form 9a (stable in aprotic dipolar solvents). The isomeric 1,4-dimethyl-5-phenylpyrazole-3(2H)-one 12 tautomerizes to the 3-hydroxypyrazole 12'. The crystal structure of the hydrochloride 14 of 9a'/9a is reported, representing the first structurally characterised example of a protonated 5-hydroxypyrazole.

Analysis of multicomponent adsorption kinetics on activated carbon

An integrated mathematical model for the kinetics of multicomponent adsorption on microporous carbon was developed. Transport in this bidisperse solid is represented by balance equations in the macropore and micropore phases, in which gas-phase diffusion dominates the mass transfer in the macropores, with the phenomenological diffusivities represented by the generalized Maxwell-Stefan (GMS) formulation. Viscous flow also contributes to the macropore fluxes and is included in the MS expressions. Diffusion of the adsorbed phase controls the mass transfer in the micro ore phase, p which is also described in a similar way by the MS method. The adsorption isotherms are represented by a new heterogeneous modified vacancy solution theory formulation of adsorption, which has proved to be a robust method for adsorption on activated carbons. The model is applied to the coadsorption and codesorption of C2H6 and C3H8 on Ajax and Norit carbon, as well as the displacement on Ajax carbon. The effect of the viscous flow in the macropore phase is not significant for the cases studied. The model accurately predicts the overshoot behavior and rollup of C2H6 during coadsorption. The prediction for the heavier compound C3H8 is always satisfactory, though at higher C3H8 mole fraction, the overshoot extent of C2H6 is overpredicted, possibly due to neglect of heat effects.

Recent advances in catalysts for methanol synthesis via hydrogenation of CO and CO2

Since the start of last century, methanol synthesis has attracted great interests because of its importance in chemical industries and its potential as an environmentally friendly energy carrier. The catalyst for the methanol synthesis has been a key area of research in order to optimize the reaction process. In the literature, the nature of the active site and the effects of the promoter and support have been extensively investigated. In this updated review, the recent progresses in the catalyst innovation, optimization of the reaction conditions, reaction mechanism, and catalyst performance in methanol synthesis are comprehensively discussed. Key issues of catalyst improvement are highlighted, and areas of priority in R&D are identified in the conclusions.

Structure of thermolysin cleaved microcin J25: Extreme stability of a two-chain antimicrobial peptide devoid of covalent links

sThe structure of a two-chain peptide formed by the treatment of the potent antimicrobial peptide microcin J25 (MccJ25) with thermolysin has been characterized by NMR spectroscopy and mass spectrometry. The native peptide is 21 amino acids in size and has the remarkable structural feature of a ring formed by linkage of the side chain of Glu8 to the N-terminus that is threaded by the C-terminal tail of the peptide. Thermolysin cleaves the peptide at the Phe10-Val11 amide bond, but the threading of the C-terminus through the N-terminal ring is so tight that the resultant two chains remain associated both in the solution and in the gas phases. The three-dimensional structure of the thermolysin-cleaved peptide derived using NMR spectroscopy and simulated annealing calculations has a well-defined core that comprises the N-terminal ring and the threading C-terminal tail. In contrast to the well-defined core, the newly formed termini at residues Phe10 and Val11 are disordered in solution. The C-terminal tail is associated to the ring both by hydrogen bonds stabilizing a short beta-sheet and by hydrophobic interactions. Moreover, unthreading of the tail through the ring is prevented by the bulky side chains of Phe19 and Tyr20, which flank the octapeptide ring. This noncovalent two-peptide complex that has a remarkable stability in solution and in highly denaturing conditions and that survives in the gas phase is the first example of such a two-chain peptide lacking disulfide or interchain covalent bonds.

Imidoylketene-alpha-oxoketenimine and alpha-oxoketene-alpha-oxoketene rearrangements. 1,3-Shifts of substituted phenyl groups

1,3-Phenyl shifts interconvert imidoylketenes 1 and alpha-oxoketenimines 2 and, likewise, alpha-oxoketenes 3 automerize by this 1,3-shift. These rearrangements usually take place in the gas phase under conditions of. ash vacuum thermolysis. Energy profiles calculated at the B3LYP/6-31G(d, p) and B3LYP/6311 + G(3df,2p)//B3LYP/6-31G(d,p) levels demonstrate that electron donating substituents ( D) in the migrating phenyl group and electron withdrawing ones ( W) in the non-migrating phenyl group, can stabilise the transition states TS1 and TS2 to the extent that activation barriers of ca. 100 kJ mol(-1) or less are obtained; i.e. enough to make these reactions potentially observable in solution at ordinary temperatures. The calculated transition state energies Delta G(TS1) show an excellent correlation with the Hammett constants sigma(p)(W) and sigma(p) +(D).

Permeability of subcritical hydrocarbons in activated carbon

A new diffusion and flow model is presented to describe the behavior of hydrocarbon vapors in activated carbon. The micro/mesopore size distribution (PSD) is obtained according to Do's method which consists of two sequential processes of pore layering and pore filling. This model uses the micro/meso PSD obtained from each adsorbate equilibrium isotherm, which reflects the dynamics behavior of adsorbing molecules through the solid. The initial rise in total permeability is mainly attributed to adsorbed-phase diffusion (that is, surface diffusion), whereas the decrease over reduced pressure of about 0.9 is attributed to the reduction of pore space available for gas phase diffusion and flow. A functional form of surface diffusivity is proposed and validated with experimental data. This model predicts well the permeability of condensable hydrocarbon vapors in activated carbon. (C) 2005 American Institute of Chemical Engineers.

Anaerobic and aerobic metabolism of glycogen-accumulating organisms selected with propionate as the sole carbon source

In the microbial competition observed in enhanced biological phosphorus removal (EBPR) systems, an undesirable group of micro-organisms known as glycogen-accumulating organisms (GAOs) compete for carbon in the anaerobic period with the desired polyphosphate-accumulating organisms (PAOs). Some studies have suggested that a propionate carbon source provides PAOs with a competitive advantage over GAOs in EBPR systems; however, the metabolism of GAOs with this carbon source has not been previously investigated. In this study, GAOs were enriched in a laboratory-scale bioreactor with propionate as the sole carbon source, in an effort to better understand their biochemical processes. Based on comprehensive solid-, liquid- and gas-phase chemical analytical data from the bioreactor, a metabolic model was proposed for the metabolism of propionate by GAOs. The model adequately described the anaerobic stoichiometry observed through chemical analysis, and can be a valuable tool for further investigation of the competition between PAOs and GAOs, and for the optimization of the EBPR process. A group of Alphaproteobacteria dominated the biomass (96% of Bacteria) from this bioreactor, while post-fluorescence in situ hybridization (FISH) chemical staining confirmed that these Alphaproteobacteria produced poly-beta-hydroxyalkanoates (PHAs) anaerobically and utilized them aerobically, demonstrating that they were putative GAOs. Some of the Alphaproteobacteria were related to Defluvicoccus vanus (16% of Bacteria), but the specific identity of many could not be determined by FISH. Further investigation into the identity of other GAOs is necessary.

On the equilibrium and dynamic behavior of alcohol vapors in activated carbon

As alcohol molecules such as methanol and ethanol have both polar and non-polar groups, their adsorption behavior is governed by the contributions of dispersion interaction (alkyl group) and hydrogen bonding (OH group). In this paper, the adsorption behavior of alcohol molecules and its effect on transport processes are elucidated. From the total permeability (B-T) of alcohol molecules in activated carbon, an adsorption mechanism is proposed, describing well the experimental data, by taking combination effects of clustering, entering micropores, layering and pore filling processes. Unlike the case of non-polar compounds, it was found that at low pressures there are two rises in the BT of alcohol molecules in activated carbon. The first rise is due to the major contribution of surface diffusion to the transport (which is the case of non-polar molecules) and the second one may be associated with cluster formation at the edge of micropores and entering micropores when the clusters are sufficiently large enough to induce a dispersive energy. In addition the clusters formed may enhance surface diffusion at low pressures and hinder gas phase diffusion and flow in meso/macropores. (c) 2006 Elsevier Ltd. All fights reserved.

Oxoketene-oxoketene, imidoylketene-imidoylketene and oxoketenimine-imidoylketene rearrangements. 1,3-shifts of phenyl groups

Dibenzoylketene 5 undergoes degenerate 1,3-shifts of the phenyl group between acyl and ketene carbon atoms, thus interconverting it with 6 and 7. This 1,3-shift takes place in the gas phase under flash vacuum thermolysis (FVT) conditions, but not in solution at 110-145 degrees C. Imidoyl(benzoyl)ketene 13 undergoes degenerate 1,3-shift of the phenyl group on FVT, thus interconverting it with 14, but the ketenimine isomer 15 is not formed, and none of these shifts take place in the solid state at 250 degrees C. Imidoyl(p-toluoyl)ketene 21 undergoes a 1,3-p-tolyl shift, interconverting it with ketene 22 but not with ketenimine 23. The imidoyl(p-toluoyl)ketene rotamer 25 cyclizes to 4-toluoyloxyquinoline 28 and 4-quinolone 29. The cyclization of imidoyl(benzoyl)ketene 13 to 4-benzoyloxyquinoline 18, and of 25 to 28 involves 1,3-C-to-O shifts of benzoyl (toluoyl) groups. Calculations of the transition states for the transformations at the B3LYP/6-31G** level of theory are in agreement with the observed reaction preferences.

High-pressure adsorption of methane and carbon dioxide on coal

Adsorption isotherms of methane and carbon dioxide on two kinds of Australian coals have been measured at three temperatures up to pressures of 20 MPa. The adsorption behavior is described by three isotherm equations: extended three-parameter, Langmuir, and Toth. Among these, the Toth equation is found to be the most suitable, yielding the most realistic values of pore volume of the coals and the adsorbed phase density. Also, the surface area of coals obtained from CO2 adsorption at 273 K is found to be the meaningful parameter which captures the CO2 adsorption capacity. A maximum in the excess amount adsorbed of each gas appears at a lower pressure with a decrease in temperature. For carbon dioxide, after the appearance of the maximum, an inflection point in the excess amount adsorbed is observed close to the critical density at each temperature, indicating that the decrease in the gas-phase density change with pressure influences the behavior of the excess amount adsorbed. In the context of CO2 sequestration, it is found that CO2 injection pressures of lower than 10 MPa may be desirable for the CH4 recovery process and CO2-holding capacity.

Effect of the preparation technique on the catalytic properties of mesoporous V-HMS for the oxidation of toluene

Various mesoporous catalysts with vanadium loadings between 0.5 and 6 V wt.% and surface areas around 1300 m(2)/g were synthesized using the isomorphous substitution (IS) and molecular designed dispersion (MDD) techniques. Their catalytic properties were tested using toluene as a model VOC in a fixed bed reactor at temperatures between 300 and 550 degrees C. It was found that during the oxidation of toluene, over V-HMS synthesized via IS, conversion of toluene mainly results in carbon oxides, benzene, benzaldehyde and water. Total conversion is greatly improved when the vanadium content is increased from around 1.5 to 3.0 wt.%, but an increase in the textural porosity (V-TEX/V-MESO) from 0.3 to 0.6 had no discernable effect on the conversion. This can be explained by the fact that a V-TEX/V-MESO as low as 0.3 is sufficient to facilitate the access of toluene into the framework confined mesopores without any molecular transport limitations. However, when using V-HMS synthesized by MDD, conversion of toluene is greatly improved when the V-TEX/ V-MESO ratio is increased from 0.1 to 0.6. This is because the diffusion limitations are minimized by this increase. V-HMS synthesized via MDD does not exhibit selectivity to benzaldehyde, favoring total oxidation to CO and CO2. This different oxidation mechanism can be explained in terms of location, accessibility and number of active species on the surface of the HMS support. (c) 2005 Elsevier Inc. All rights reserved.

Aryliminopropadienone-C-amidoketenimine-amidinoketene-2-aminoquinolone cascades and the ynamine-isocyanate reaction

Imidoylketenes 11 and oxoketenimines 12 are generated by flash vacuum thermolysis of Meldrum's acid derivatives 9, pyrrolediones 17 and 18, and triazole 19 and are observed by IR spectroscopy. Ketenimine-3-carboxylic acid esters 12a are isolable at room temperature. Ketenes 11 and ketenimines 12 undergo rapid interconversion in the gas phase, and the ketenes cyclize to 4-quinolones 13. When using an amine leaving group in Meldrum's acid derivatives 9c, the major reaction products are aryliminopropadienones, ArN=C=C=C=O (15). The latter react with 1 equiv of nucleophile to produce ketenimines 12 and with 2 equiv to afford maIonic acid imide derivatives 16. N-Arylketenimine-C-carboxamides 12c cyclize to quinolones 13c via the transient amidinoketenes 11c at temperatures of 25-40 degrees C. This implies rapid interconversion of ketenes and ketenimines by a 1,3-shift of the dimethylamino group, even at room temperature. This interconversion explains previously poorly understood outcomes of the ynamine-isocyanate reaction. The solvent dependence of the tautomerism of 4-quinolones/4-quinolinols is discussed. Rotational barriers of NMe2 groups in amidoketenimines 12c and malonioc amides and amidines 16 (24) are reported.