TAP Reactor Development and Application in the Kinetic Characterization of Catalysts for Heterogeneously Catalyzed Gas Phase Reactions

This manuscript describes the application and further development of the TAP technique in kinetic characterization of heterogeneous catalysis. The major application of TAP systems is to study mechanisms, kinetics and transport phenomena in heterogeneous catalysis, all of which is made possible by the sub-millisecond time resolution. Furthermore, the kinetic information obtained can be used to gain an insight into the mechanism occurring over the catalyst system. This is advantageous as heterogeneous catalysts with an improved efficiency can be developed as a result. TAP kinetic studies are carried out at low pressure (~1x10-7 mbar) and TAP pulses are sufficiently small (1013-1015 molecules) so as to maintain this low pressure. The use of a small number of molecules in comparison to the total number of active sites means the state of the catalyst remains relatively unchanged. The use of the low intensity pulses also makes the pressure gradient negligible and so allows the TAP reactor system to operate in the Knudsen Diffusion regime, where gas-gas reactions are eliminated. Hence only gas-catalyst reactions are investigated and, by the use of moment analysis of observed exit flow, rate constants of elementary steps of the reaction can be obtained.

In this manuscript, two attempts to further the TAP technique are reported. Firstly, the work undertaken at QUB to attempt to control the number of molecules of condensable reagents that can be pulsed during a TAP pulse experiment is disclosed. Secondly, a collaborative project with SAI Ltd Manchester is discussed in a separate chapter, where technical details and validation of a customised time of flight mass spectrometer (ToF MS) for the QUB TAP-1 system are reported. A collaborative project with Cardiff Catalysis Institute focusing on the study of CO oxidation over hopcalite catalysts is also reported. The analysis of the experimental results has provided an insight into the possible mechanism of the oxidation of CO over these catalysts. A correction function has also been derived which accounts for the adsorption of reactant molecules over inert materials that are used for the reactor packing in TAP experiments. This function was then applied to the selective reduction of O2 in a H2 rich ethene feed, so that more accurate TAP moment based analysis could be conducted.

Possibility of designing catalysts beyond the traditional volcano curve: A theoretical framework for multi-phase surfaces

The current theory of catalyst activity in heterogeneous catalysis is mainly obtained from the study of catalysts with mono-phases, while most catalysts in real systems consist of multi-phases, the understanding of which is far short of chemists' expectation. Density functional theory (DFT) and micro-kinetics simulations are used to investigate the activities of six mono-phase and nine bi-phase catalysts, using CO hydrogenation that is arguably the most typical reaction in heterogeneous catalysis. Excellent activities that are beyond the activity peak of traditional mono-phase volcano curves are found on some bi-phase surfaces. By analyzing these results, a new framework to understand the unexpected activities of bi-phase surfaces is proposed. Based on the framework, several principles for the design of multi-phase catalysts are suggested. The theoretical framework extends the traditional catalysis theory to understand more complex systems.

Modeling reaction kinetics and mass transfer in ozonation in water solutions

This dissertation is based on four articles dealing with modeling of ozonation. The literature part of this considers some models for hydrodynamics in bubble column simulation. A literature review of methods for obtaining mass transfer coefficients is presented. The methods presented to obtain mass transfer are general models and can be applied to any gas-liquid system. Ozonation reaction models and methods for obtaining stoichiometric coefficients and reaction rate coefficients for ozonation reactions are discussed in the final section of the literature part. In the first article, ozone gas-liquid mass transfer into water in a bubble column was investigated for different pH values. A more general method for estimation of mass transfer and Henry’s coefficient was developed from the Beltrán method. The ozone volumetric mass transfer coefficient and the Henry’s coefficient were determined simultaneously by parameter estimation using a nonlinear optimization method. A minor dependence of the Henry’s law constant on pH was detected at the pH range 4 - 9. In the second article, a new method using the axial dispersion model for estimation of ozone self-decomposition kinetics in a semi-batch bubble column reactor was developed. The reaction rate coefficients for literature equations of ozone decomposition and the gas phase dispersion coefficient were estimated and compared with the literature data. The reaction order in the pH range 7-10 with respect to ozone 1.12 and 0.51 the hydroxyl ion were obtained, which is in good agreement with literature. The model parameters were determined by parameter estimation using a nonlinear optimization method. Sensitivity analysis was conducted using object function method to obtain information about the reliability and identifiability of the estimated parameters. In the third article, the reaction rate coefficients and the stoichiometric coefficients in the reaction of ozone with the model component p-nitrophenol were estimated at low pH of water using nonlinear optimization. A novel method for estimation of multireaction model parameters in ozonation was developed. In this method the concentration of unknown intermediate compounds is presented as a residual COD (chemical oxygen demand) calculated from the measured COD and the theoretical COD for the known species. The decomposition rate of p-nitrophenol on the pathway producing hydroquinone was found to be about two times faster than the p-nitrophenol decomposition rate on the pathway producing 4- nitrocatechol. In the fourth article, the reaction kinetics of p-nitrophenol ozonation was studied in a bubble column at pH 2. Using the new reaction kinetic model presented in the previous article, the reaction kinetic parameters, rate coefficients, and stoichiometric coefficients as well as the mass transfer coefficient were estimated with nonlinear estimation. The decomposition rate of pnitrophenol was found to be equal both on the pathway producing hydroquinone and on the path way producing 4-nitrocathecol. Comparison of the rate coefficients with the case at initial pH 5 indicates that the p-nitrophenol degradation producing 4- nitrocathecol is more selective towards molecular ozone than the reaction producing hydroquinone. The identifiability and reliability of the estimated parameters were analyzed with the Marcov chain Monte Carlo (MCMC) method. @All rights reserved. No part of the publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior permission of the author.

Mass transfer and reaction kinetics in single droplet systems

As reactive extraction grown more and more popular in a variety of technological applications, optimizing its performance becomes more and more important. The process of complex formation is affected by a great number of both physical and chemical properties of all the components involved, and sometimes their interference with one another makes improving the effectiveness of such processes very difficult. In this Master’s Theses, the processes of complex formation between the aqueous phase - represented by copper sulfate water solution, and organic phase – represented by Acorga M5640 solvent extractor, were studied in order to establish the effect these components have on reactive extraction performance and to determine which step is bottlenecking the process the most.

Oxidation kinetics of nickel nano crystallites obtained by controlled thermolysis of diaquabis(ethylenediamine) nickel(II) nitrate

The metal complex, [Ni(en)2(H2O)2](NO3)2 (en = ethylenediamine), was decomposed in a static furnace at 200 C by autogenous decomposition to obtain phase pure metallic nickel nanocrystallites. The nickel metal thus obtained was studied by XRD, IR spectra, SEM and CHN analysis. The nickel crystallites are in the nanometer range as indicated by XRD studies. The IR spectral studies and CHN analyses show that the surface is covered with a nitrogen containing species. Thermogravimetric mass gain shows that the product purity is high (93%). The formed nickel is stable and resistant to oxidation up to 350 C probably due to the coverage of nitrogen containing species. Activation energy for the oxidation of the prepared nickel nanocrystallites was determined by non-isothermal methods and was found to depend on the conversion ratio. The oxidation kinetics of the nickel crystallites obeyed a Johnson–Mehl–Avrami mechanism probably due to the special morphology and crystallite strain present on the metal.

Kinetics of Bacterial Colony Growth by Laser Induced Fluorescence

The growth kinetics of an aerial bacterial colony on solid agar media was studied using laser induced fluorescence technique. Fluorescence quenching of Rhodamin B by the bacterial colony was utilized for the study. The lag phase, log phase, and stationary phase of growth curve of bacterial colony was identified by measuring peak fluorescence intensity of dye doped bacterial colony.

Comment on "Kinetics of spinodal decomposition in the ising model with vacancy diffusion"

In a recent paper [Phys. Rev. B 50, 3477 (1994)], P. Fratzl and O. Penrose present the results of the Monte Carlo simulation of the spinodal decomposition problem (phase separation) using the vacancy dynamics mechanism. They observe that the t1/3 growth regime is reached faster than when using the standard Kawasaki dynamics. In this Comment we provide a simple explanation for the phenomenon based on the role of interface diffusion, which they claim is irrelevant for the observed behavior.

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Reply to "Comment on 'Nature and entropy content of the ordering transitions in RCo2'"

In the previous Comment, Forker and co-workers claim that perturbed angular correlation (PAC) data leave no alternative to the conclusion that the spontaneous magnetization of PrCo2 and NdCo2 undergoes a discontinuous, first-order phase transition at TC. We show here that their claim is in clear contradiction with a wealth of experimental evidence, including our own. Finally, we propose a possible origin for the disagreement between their interpretation of the PAC results and the literature on this subject.

Numerical model of solid phase transformations governed by nucleation and growth: microstructure development during isothermal crystallization

A simple numerical model which calculates the kinetics of crystallization involving randomly distributed nucleation and isotropic growth is presented. The model can be applied to different thermal histories and no restrictions are imposed on the time and the temperature dependences of the nucleation and growth rates. We also develop an algorithm which evaluates the corresponding emerging grain-size distribution. The algorithm is easy to implement and particularly flexible, making it possible to simulate several experimental conditions. Its simplicity and minimal computer requirements allow high accuracy for two- and three-dimensional growth simulations. The algorithm is applied to explore the grain morphology development during isothermal treatments for several nucleation regimes. In particular, thermal nucleation, preexisting nuclei, and the combination of both nucleation mechanisms are analyzed. For the first two cases, the universal grain-size distribution is obtained. The high accuracy of the model is stated from its comparison to analytical predictions. Finally, the validity of the Kolmogorov-Johnson-Mehl-Avrami model SSSR, is verified for all the cases studied

Pressure-jump-induced kinetics reveals a hydration dependent folding/unfolding mechanism of ribonuclease A

Pressure-jump (p-jump)-induced relaxation kinetics was used to explore the energy landscape of protein folding/unfolding of Y115W, a fluorescent variant of ribonuclease A. Pressure-jumps of 40MPa amplitude (5ms dead-time) were conducted both to higher (unfolding) and to lower (folding) pressure, in the range from 100 to 500MPa, between 30 and 50°C. Significant deviations from the expected symmetrical protein relaxation kinetics were observed. Whereas downward p-jumps resulted always in single exponential kinetics, the kinetics induced by upward p-jumps were biphasic in the low pressure range and monophasic at higher pressures. The relative amplitude of the slow phase decreased as a function of both pressure and temperature. At 50°C, only the fast phase remained. These results can be interpreted within the framework of a two-dimensional energy surface containing a pressure- and temperature-dependent barrier between two unfolded states differing in the isomeric state of the Asn-113–Pro-114 bond. Analysis of the activation volume of the fast kinetic phase revealed a temperature-dependent shift of the unfolding transition state to a larger volume. The observed compensation of this effect by glycerol offers an explanation for its protein stabilizing effect

Excision of anabaena PCC 7120 nifD element in escherichia coli: growth kinetics and RecA regulated xisA expression and DNA rearrangement

Anabaena PCC 7120 nifHDK operon is interrupted by an 11 kb DNA element which is excised during the development of heterocysts by Excisase A, encoded by the xisA gene residing on the element. The excision is a site-specific recombination event that occurs at the I I base pair direct repeats flanking the element. Earlier work showed the excision of the I I kb element in Escherichia coli at a frequency 0.3%. We report here the excision of this element at 1.1% and 1.98% in E. coli DH5 alpha, and 1.9% and 10.9% in E. coli JM 101 when grown on Luria broth and minimal media, respectively. Excision of nifD element in isogenic recA(-) (RK1) and recA(+) (RK2) E. coli JM101 P1 transductants, showed similar results to that of E. coli JM101 and DH5 alpha, respectively. A plasmid pMX32, carrying a xisA defective 11 kb element, showed no excision in E. coli RK2 strain. In contrast to Anabaena PCC 7120, excision of nifD element did not increase in E. call DH5 alpha grown in iron-deficient conditions. A PxisA::lacZ transcriptional fusion, used to detect the expression of elusive xisA gene, showed maximal beta-galactosidase activity in the stationary phase. The results suggest that the excision event in E. coli may involve additional factors, such as RecA and that the physiological status can influence the excision of nifD element. (C) 2007 Elsevier Ltd. All rights reserved.

Direct detection of chlorosilylene and time resolved study of some of its reactions in the gas-phase using a new photochemical precursor

Chlorosilylene, ClSiH, was prepared by 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene in the gas phase. ClSiH was monitored in real time at 457.9 nm using a CW argon ion laser. The kinetics of reactions of ClSiH with C2H4, CH2 = CHCMe3, C2H2, Me2O, SO2, HCl, MeSiH3, Me2SiH2, Me3SiH, MeGeH3, MeGeH3 and precursor have been studied at ambient temperatures for the first time. Addition reactions of ClSiH and reactions with lone pair donors are faster than insertion reactions. Surprisingly ClSiH inserts faster into Si-H than Ge-H bonds. ClSiH is intermediate in reactivity between SiH2 and SiCI2. Relative reactivities of CISiH and SiH2 vary considerably. (c) 2005 Elsevier B.V. All rights reserved.

Gas phase kinetic and quantum chemical studies of the reactions of silylene with the methylsilanes. Absolute rate constants, temperature dependences, RRKM modelling and potential energy surfaces

Time resolved studies of silylene, SiH2, generated by the 193 nm laser. ash photolysis of phenylsilane, have been carried out to obtain rate coefficients for its bimolecular reactions with methyl-, dimethyl- and trimethyl-silanes in the gas phase. The reactions were studied over the pressure range 3 - 100 Torr with SF6 as bath gas and at five temperatures in the range 300 - 625 K. Only slight pressure dependences were found for SiH2 + MeSiH3 ( 485 and 602 K) and for SiH2 + Me2SiH2 ( 600 K). The high pressure rate constants gave the following Arrhenius parameters: [GRAPHICS] These are consistent with fast, near to collision-controlled, association processes. RRKM modelling calculations are consistent with the observed pressure dependences ( and also the lack of them for SiH2 + Me3SiH). Ab initio calculations at both second order perturbation theory (MP2) and coupled cluster (CCSD(T)) levels, showed the presence of weakly-bound complexes along the reaction pathways. In the case of SiH2 + MeSiH3 two complexes, with different geometries, were obtained consistent with earlier studies of SiH2 + SiH4. These complexes were stabilised by methyl substitution in the substrate silane, but all had exceedingly low barriers to rearrangement to product disilanes. Although methyl groups in the substrate silane enhance the intrinsic SiH2 insertion rates, it is doubtful whether the intermediate complexes have a significant effect on the kinetics. A further calculation on the reaction MeSiH + SiH4 shows that the methyl substitution in the silylene should have a much more significant kinetic effect ( as observed in other studies).

A pressure-jump time-resolved x-ray diffraction study of cubic-cubic transition kinetics in monoolein

In the past two decades, the geometric pathways involved in the transformations between inverse bicontinuous cubic phases in amphiphilic systems have been extensively theoretically modeled. However, little experimental data exists on the cubic-cubic transformation in pure lipid systems. We have used pressure-jump time-resolved X-ray diffraction to investigate the transition between the gyroid Q(II)(G) and double-diamond Q(II)(D) phases in mixtures of 1-monoolein in 30 wt% water. We find for this system that the cubic-cubic transition occurs without any detectable intermediate structures. In addition, we have determined the kinetics of the transition, in both the forward and reverse directions, as a function of pressure-jump amplitude, temperature, and water content. A recently developed model allows (at least in principle) the calculation of the activation energy for lipid phase transitions from such data. The analysis is applicable only if kinetic reproducibility is achieved, at least within one sample, and achievement of such kinetic reproducibility is shown here, by carrying out prolonged pressure-cycling. The rate of transformation shows clear and consistent trends with pressure-jump amplitude, temperature, and water content, all of which are shown to be in agreement with the effect of the shift in the position of the cubic-cubic phase boundary following a change in the thermodynamic parameters.