Dynamical features of reaction-diffusion fronts in fractals

The speed of front propagation in fractals is studied by using (i) the reduction of the reaction-transport equation into a Hamilton-Jacobi equation and (ii) the local-equilibrium approach. Different equations proposed for describing transport in fractal media, together with logistic reaction kinetics, are considered. Finally, we analyze the main features of wave fronts resulting from this dynamic process, i.e., why they are accelerated and what is the exact form of this acceleration

Metamorphism and kinetics in the Torres del Paine contact aureole

Contact aureoles provide an excellent geologic environment to study the mechanisms of metamorphic reactions in a natural system. The Torres del Paine (TP) intrusion is one of the most spectacular natural laboratories because of its excellent outcrop conditions. It formed in a period from 12.59 to 12.43 Ma and consists of three large granite and four smaller mafic batches. The oldest granite is on top, the youngest at the bottom of the granitic complex, and the granites overly the mafic laccolith. The TP intruded at a depth of 2-3 km into regional metamorphic anchizone to greenschist facies pelites, sandstones, and conglomerates of the Cerro Toro and Punta Barrosa formations. It formed a thin contact aureole of 150-400 m width. This thesis focuses on the reaction kinetics of the mineral cordierite in the contact aureole using quantitative textural analysis methods. First cordierite was formed from chlorite break¬down (zone I, ca. 480 °C, 750 bar). The second cordierite forming reaction was the muscovite break-down, which is accompanied by a modal decrease in biotite and the appearance of k- feldspar (zone II, 540-550 °C, 750 bar). Crystal sizes of the roundish, poikiloblastic cordierites were determined from microscope thin section images by manually marking each crystal. Images were then automatically processed with Matlab. The correction for the intersection probability of each crystal radius yields the crystal size distribution in the rock. Samples from zone I below the laccolith have the largest crystals (0.09 mm). Cordierites from zone II are smaller, with a maximum crystal radius of 0.057 mm. Rocks from zone II have a larger number of small cordierite crystals than rocks from zone I. A combination of these quantitative analysis with numerical modeling of nucleation and growth, is used to infer nucleation and growth parameters which are responsible for the observed mineral textures. For this, the temperature-time paths of the samples need to be known. The thermal history is complex because the main body of the intrusion was formed by several intrusive batches. The emplacement mechanism and duration of each batch can influence the thermal structure in the aureole. A possible subdivision of batches in smaller increments, so called pulses, will focus heat at the side of the intrusion. Focusing all pulses on one side increases the contact aureole size on that side, but decreases it on the other side. It forms a strongly asymmetric contact aureole. Detailed modeling shows that the relative thicknesses of the TP contact aureole above and below the intrusion (150 and 400 m) are best explained by a rapid emplacement of at least the oldest granite batch. Nevertheless, temperatures are significantly too low in all models, compared to observed mineral assemblages in the hornfelses. Hence, an other important thermal mechanisms needs to take place in the host rock. Clastic minerals in the immature sediments outside the contact aureole are hydrated due to small amounts of expelled fluids during contact metamorphism. This leads to a temperature increase of up to 50 °C. The origin of fluids can be traced by stable isotopes. Whole rock stable isotope data (6D and δ180) and chlorine concentrations in biotite document that the TP intrusion induced only very small amounts of fluid flow. Oxygen whole rock data show δ180 values between 9.0 and 10.0 %o within the first 5 m of the contact. Values increase to 13.0 - 15.0 %o further away from the intrusion. Whole rock 6D values display a more complex zoning. First, host rock values (-90 to -70 %o) smoothly decrease towards the contact by ca. 20 %o, up to a distance of ca. 150 m. This is followed by an increase of ca. 20 %o within the innermost 150 m of the aureole (-97.0 to -78 %o at the contact). The initial decrease in 6D values is interpreted to be due to Rayleigh fractionation accompanying the dehydration reactions forming cordierite, while the final increase reflects infiltration of water-rich fluids from the intrusion. An over-estimate on the quantity and the corresponding thermal effect yields a temperature increase of less than 30 °C. This suggests that fluid flow might have contributed only for a small amount to the thermal evolution of the system. A combination of the numerical growth model with the thermal model, including the hydration reaction enthalpies but neglecting fluid flow and incremental growth, can be used to numerically reproduce the observed cordierite textures in the contact aureole. This yields kinetic parameters which indicate fast cordierite crystallization before the thermal peak in the inner aureole, and continued reaction after the thermal peak in the outermost aureole. Only small temperature dependencies of the kinetic parameters seem to be needed to explain the obtained crystal size data. - Les auréoles de contact offrent un cadre géologique privilégié pour l'étude des mécanismes de réactions métamorphiques associés à la mise en place de magmas dans la croûte terrestre. Par ses conditions d'affleurements excellentes, l'intrusion de Torres del Paine représente un site exceptionnel pour améliorer nos connaissances de ces processus. La formation de cette intrusion composée de trois injections granitiques principales et de quatre injections mafiques de volume inférieur couvre une période allant de 12.50 à 12.43 Ma. Le plus vieux granite forme la partie sommitale de l'intrusion alors que l'injection la plus jeune s'observe à la base du complexe granitique; les granites recouvrent la partie mafique du laccolite. L'intrusion du Torres del Paine s'est mise en place a 2-3 km de profondeur dans un encaissant métamorphique. Cet encaissant est caractérisé par un métamorphisme régional de faciès anchizonal à schiste vert et est composé de pélites, de grès, et des conglomérats des formations du Cerro Toro et Punta Barrosa. La mise en place des différentes injections granitiques a généré une auréole de contact de 150-400 m d'épaisseur autour de l'intrusion. Cette thèse se concentre sur la cinétique de réaction associée à la formation de la cordiérite dans les auréoles de contact en utilisant des méthodes quantitatives d'analyses de texture. On observe plusieurs générations de cordiérite dans l'auréole de contact. La première cordiérite est formée par la décomposition de la chlorite (zone I, environ 480 °C, 750 bar), alors qu'une seconde génération de cordiérite est associée à la décomposition de la muscovite, laquelle est accompagnée par une diminution modale de la teneur en biotite et l'apparition de feldspath potassique (zone II, 540-550 °C, 750 bar). Les tailles des cristaux de cordiérites arrondies et blastic ont été déterminées en utilisant des images digitalisées des lames minces et en marquant individuellement chaque cristal. Les images sont ensuite traitées automatiquement à l'aide du programme Matlab. La correction de la probabilité d'intersection en fonction du rayon des cristaux permet de déterminer la distribution de la taille des cristaux dans la roche. Les échantillons de la zone I, en dessous du lacolite, sont caractérisés par de relativement grands cristaux (0.09 mm). Les cristaux de cordiérite de la zone II sont plus petits, avec un rayon maximal de 0.057 mm. Les roches de la zone II présentent un plus grand nombre de petits cristaux de cordiérite que les roches de la zone I. Une combinaison de ces analyses quantitatives avec un modèle numérique de nucléation et croissance a été utilisée pour déduire les paramètres de nucléation et croissance contrôlant les différentes textures minérales observées. Pour développer le modèle de nucléation et de croissance, il est nécessaire de connaître le chemin température - temps des échantillons. L'histoire thermique est complexe parce que l'intrusion est produite par plusieurs injections successives. En effet, le mécanisme d'emplace¬ment et la durée de chaque injection peuvent influencer la structure thermique dans l'auréole. Une subdivision des injections en plus petits incréments, appelés puises, permet de concentrer la chaleur dans les bords de l'intrusion. Une mise en place préférentielle de ces puises sur un côté de l'intrusion modifie l'apport thermique et influence la taille de l'auréole de contact produite, auréole qui devient asymétrique. Dans le cas de la première injection de granite, une modélisation détaillée montre que l'épaisseur relative de l'auréole de contact de Torres del Paine au-dessus et en dessous de l'intrusion (150 et 400 m) est mieux expliquée par un emplacement rapide du granite. Néanmoins, les températures calculées dans l'auréole de con¬tact sont trop basses pour que les modèles thermiques soient cohérants par rapport à la taille de cette auréole. Ainsi, un autre mecanisme exothermique est nécessaire pour permettre à la roche encais¬sante de produire les assemblages observés. L'observation des roches encaissantes entourant les granites montre que les minéraux clastiques dans les sédiments immatures au-dehors de l'auréole sont hydratés suite à la petite quantité de fluide expulsée durant le métamorphisme de contact et/ou la mise en place des granites. Les réactions d'hydratation peuvent permettre une augmentation de la température jusqu'à 50 °C. Afin de déterminer l'origine des fluides, une étude isotopique de roches de l'auréole de contact a été entreprise. Les isotopes stables d'oxygène et d'hydrogène sur la roche totale ainsi que la concentration en chlore dans la biotite indiquent que la mise en place des granites du Torres del Paine n'induit qu'une circulation de fluide limitée. Les données d'oxygène sur roche totale montrent des valeurs δ180 entre 9.0 et 10.0%o au sein des cinq premiers mètres du contact. Les valeurs augmentent jusqu'à 13.0 - 15.0 plus on s'éloigne de l'intrusion. Les valeurs 5D sur roche totale montrent une zonation plus complexe. Les valeurs de la roche encaissante (-90 à -70%o) diminuent progressivement d'environ 20%o depuis l'extérieur de l'auréole jusqu'à une distance d'environ 150 m du granite. Cette diminution est suivie par une augmentation d'environ 20%o au sein des 150 mètres les plus proches du contact (-97.0 à -78%o au contact). La diminution initiale des valeurs de 6D est interprétée comme la conséquence du fractionnement de Rayleigh qui accompagne les réactions de déshydratation formant la cordiérite, alors que l'augmentation finale reflète l'infiltration de fluide riche en eau venant de l'intrusion. A partir de ces résultats, le volume du fluide issu du granite ainsi que son effet thermique a pu être estimé. Ces résultats montrent que l'augmentation de température associée à ces fluides est limitée à un maximum de 30 °C. La contribution de ces fluides dans le bilan thermique est donc faible. Ces différents résultats nous ont permis de créer un modèle thermique associé à la for¬mation de l'auréole de contact qui intègre la mise en place rapide du granite et les réactions d'hydratation lors du métamorphisme. L'intégration de ce modèle thermique dans le modèle numérique de croissance minérale nous permet de calculer les textures des cordiérites. Cepen¬dant, ce modèle est dépendant de la vitesse de croissance et de nucléation de ces cordiérites. Nous avons obtenu ces paramètres en comparant les textures prédites par le modèle et les textures observées dans les roches de l'auréole de contact du Torres del Paine. Les paramètres cinétiques extraits du modèle optimisé indiquent une cristallisation rapide de la cordiérite avant le pic thermique dans la partie interne de l'auréole, et une réaction continue après le pic thermique dans la partie la plus externe de l'auréole. Seules de petites dépendances de température des paramètres de cinétique semblent être nécessaires pour expliquer les don¬nées obtenues sur la distribution des tailles de cristaux. Ces résultats apportent un éclairage nouveau sur la cinétique qui contrôle les réactions métamorphiques.

Kinetics of selective oxidation of lactose to lactobionic acid over Pd-active carbon catalyst

Laktoosi eli maitosokeri on tärkein ainesosa useimpien nisäkkäiden tuottamassa maidossa. Sitä erotetaan herasta, juustosta ja maidosta. Laktoosia käytetään elintarvike- ja lääketeollisuuden raaka-aineena monissaeri tuotteissa. Lääketeollisuudessa laktoosia käytetään esimerkiksi tablettien täyteaineena. Hapettamalla laktoosia voidaan valmistaa laktobionihappoa, 2-keto-laktobionihappoa ja laktuloosia. Laktobionihappoa käytetään biohajoavien pintojen ja kosmetiikkatuotteiden valmistuksessa, sekä sisäelinten säilöntäliuoksissa, joissa laktobionihappo estää happiradikaalien aiheuttamien kudosvaurioiden syntymistä. Tässä työssä laktoosia hapetettiin laktobionihapoksi sekoittimella varustetussa laboratoriomittakaavaisessa panosreaktorissa käyttäenkatalyyttinä palladiumia aktiivihiilellä. Muutamissa kokeissa katalyytin promoottorina käytettiin vismuttia, joka hidastaa katalyytin deaktivoitumista. Työn tarkoituksena oli saada lisää tietoa laktoosin hapettamisen kinetiikasta. Laktoosin hapettumisessa laktobionihapoksi havaittiin selektiivisyyteen vaikuttavan muunmuassa reaktiolämpötila, paine, pH ja käytetyn katalyytin määrä. Katalyyttiä kierrättämällä eri kokeiden välillä saatiin paremmat konversiot, selektiivisyydet ja saannot. Parhaat koetulokset saatiin hapetettaessa synteettisellä ilmalla 60 oC lämpötilassa ja 1 bar paineessa. Tehdyissä kokeissa pH:n säätö tehtiin manuaalisesti, joten pH ei pysynyt koko ajan haluttuna. Laktoosin konversio oli parhaimmillaan 95 %. Laktobionihapon suhteellinen selektiivisyys oli 100% ja suhteellinen saanto 100 %. Kinetiikan matemaattinen mallinnus tehtiin Modest-ohjelmalla käyttäen kokeista saatuja mittaustuloksia.Ohjelman avulla estimoitiin parametreja ja saatiin matemaattinen malli reaktorille. Tässä työssä tehtiin kineettinen mallinnus myös ravistelureaktorissa tehdyille laktoosin hapetuskokeille, missä pH pysyi koko ajan haluttuna 'in-situ' titrauksen avulla. Työn yhteydessä selvitettiin myös mahdollisuutta käyttää monoliittikatalyyttejä laktoosin hapetusreaktiossa.

Photoinhibition of Photosystem II. Kinetics, Photoprotection and Mechanism

Photosystem II (PSII) is susceptible to light-induced damage defined as photoinhibition. In natural conditions, plants are capable of repairing the photoinhibited PSII by on-going degradation and re-synthesis of the D1 reaction centre protein of PSII. Photoinhibition is induced by both visible and ultraviolet light and photoinhibition occurs under all light intensities with the same efficiency per photon. In my thesis work, I studied the reaction kinetics and mechanism of photoinhibition of PSII, as well as photoprotection in leaves of higher plants. Action spectroscopy was used to identify photoreceptors of photoinhibition. I found that the action spectrum of photoinhibition in vivo shows resemblance to the absorption spectra of manganese model compounds of the oxygen evolving complex (OEC) suggesting a role for manganese as a photoreceptor of photoinhibition under UV and visible light. In order to study the protective effect of non-photochemical quenching, the action spectrum was measured from leaves of wild type Arabidopsis thaliana and two mutants impaired in nonphotochemical quenching of chlorophyll a excitations. The findings of action spectroscopy and simulations of chlorophyll-based photoinhibition mechanisms suggested that quenching of antenna excitations protects less efficiently than would be expected if antenna chlorophylls were the only photoreceptors of photoinhibition. The reaction kinetics of prolonged photoinhibition was studied in leaves of Cucurbita maxima and Capsicum annuum. The results indicated that photoinhibitory decrease in both the oxygen evolution activity and ratio of variable to maximum fluorescence follows firstorder kinetics in vivo. The persistence of first-order kinetics suggests that already photoinhibited reaction centres do not protect against photoinhibition and that the mechanism of photoinhibition does not have a reversible intermediate. When Cucurbita maxima leaves were photoinhibited with saturating single-turnover flashes and continuous light, the light response curve of photoinhibition was found to be essentially a straight line with both types of illumination, suggesting that similar photoinhibition mechanisms might function during illumination with continuous light and during illumination with short flashes.

A reaction engineering approach to homogeneously catalyzed glycerol hydrochlorination

The production of biodiesel through transesterification has created a surplus of glycerol on the international market. In few years, glycerol has become an inexpensive and abundant raw material, subject to numerous plausible valorisation strategies. Glycerol hydrochlorination stands out as an economically attractive alternative to the production of biobased epichlorohydrin, an important raw material for the manufacturing of epoxy resins and plasticizers. Glycerol hydrochlorination using gaseous hydrogen chloride (HCl) was studied from a reaction engineering viewpoint. Firstly, a more general and rigorous kinetic model was derived based on a consistent reaction mechanism proposed in the literature. The model was validated with experimental data reported in the literature as well as with new data of our own. Semi-batch experiments were conducted in which the influence of the stirring speed, HCl partial pressure, catalyst concentration and temperature were thoroughly analysed and discussed. Acetic acid was used as a homogeneous catalyst for the experiments. For the first time, it was demonstrated that the liquid-phase volume undergoes a significant increase due to the accumulation of HCl in the liquid phase. Novel and relevant features concerning hydrochlorination kinetics, HCl solubility and mass transfer were investigated. An extended reaction mechanism was proposed and a new kinetic model was derived. The model was tested with the experimental data by means of regression analysis, in which kinetic and mass transfer parameters were successfully estimated. A dimensionless number, called Catalyst Modulus, was proposed as a tool for corroborating the kinetic model. Reactive flash distillation experiments were conducted to check the commonly accepted hypothesis that removal of water should enhance the glycerol hydrochlorination kinetics. The performance of the reactive flash distillation experiments were compared to the semi-batch data previously obtained. An unforeseen effect was observed once the water was let to be stripped out from the liquid phase, exposing a strong correlation between the HCl liquid uptake and the presence of water in the system. Water has revealed to play an important role also in the HCl dissociation: as water was removed, the dissociation of HCl was diminished, which had a retarding effect on the reaction kinetics. In order to obtain a further insight on the influence of water on the hydrochlorination reaction, extra semi-batch experiments were conducted in which initial amounts of water and the desired product were added. This study revealed the possibility to use the desired product as an ideal “solvent” for the glycerol hydrochlorination process. A co-current bubble column was used to investigate the glycerol hydrochlorination process under continuous operation. The influence of liquid flow rate, gas flow rate, temperature and catalyst concentration on the glycerol conversion and product distribution was studied. The fluid dynamics of the system showed a remarkable behaviour, which was carefully investigated and described. Highspeed camera images and residence time distribution experiments were conducted to collect relevant information about the flow conditions inside the tube. A model based on the axial dispersion concept was proposed and confronted with the experimental data. The kinetic and solubility parameters estimated from the semi-batch experiments were successfully used in the description of mass transfer and fluid dynamics of the bubble column reactor. In light of the results brought by the present work, the glycerol hydrochlorination reaction mechanism has been finally clarified. It has been demonstrated that the reactive distillation technology may cause drawbacks to the glycerol hydrochlorination reaction rate under certain conditions. Furthermore, continuous reactor technology showed a high selectivity towards monochlorohydrins, whilst semibatch technology was demonstrated to be more efficient towards the production of dichlorohydrins. Based on the novel and revealing discoveries brought by the present work, many insightful suggestions are made towards the improvement of the production of αγ-dichlorohydrin on an industrial scale.

Statistics and kinetics of single-molecule electron transfer dynamics in complex environments: A simulation model study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Effects of temperature and of the addition of accelerating and retarding agents on the kinetics of hydration of tricalcium silicate

The formation of calcium silicate hydrates (C-S-H) during the hydration of tricalcium silicate (C3S) in pure water and in water solutions containing 1% CaCl2 (accelerator) and 0.01% saccharose (retarder) was studied by small-angle X-ray scattering (SAXS). SAXS measurements were performed under isothermal conditions within the temperature range 25 °C T < 52 °C. The experimental results indicate that the time variation of the mass fraction of the C-S-H product phase, α(f), can be fitted, under all conditions of paste setting, by Avrami equation, α(t) = 1 -exp(-(kt)′), k being a rate parameter and n an exponent depending on the characteristics of the transformation. The parameter n is approximately equal to 2 for hydration of C^S in pure water. Depending on temperature, n varies from 2 to 2.65 for hydration in the presence of CaC^ and saccharose. The value n = 2 is theoretically expected for lateral growth of thin C-S-H plates of constant thickness. The time dependence of SAXS intensity indicates that the transformed phase (C-S-H) consists of colloidal particles in early stages of hydration, evolving by two-dimensional growth toward a disordered lamellar structure composed of very thin plates. The activation energy ΔE for the growth of C-S-H phase was determined from the time dependence of X-ray scattering intensity. These data were obtained by in situ measurements at different temperatures of hydration. The values of ΔE are 37.7, 49.4, and 44.3 kJ/mol for hydration in pure water and in water solutions containing CaCl2 and saccharose, respectively. © 2000 American Chemical Society.

Thermogravimetric studies of the decomposition kinetics of four different Hevea rubber clones using Ozawa's approach

Thermal degradation of natural rubber extracted from four different rubber clones of the Hevea brasiliensis species was investigated by thermogravimetry using Ozawa's approach to assess the kinetic parameters of the decomposition process. The results are discussed in terms of the order of reaction, kinetic parameters such as activation energy as a function of conversion degree, and thermal stability for each of the clones.

Application of non-isothermal cure kinetics on the interaction of poly(ethylene terephthalate) - Alkyd resin paints

Samples of paint (P), reused PET (PET-R) and paint/PET-R mixtures (PPET-R) were evaluated using DSC to verify their physical-chemical properties and thermal behavior. Films from paints and PPET-R are visually similar. It was possible to establish that the maximum amount of PET-R that can be added to paint without significantly altering its filming properties is 2%. The cure process (80-203°C) was identified through DSC curves. The kinetic parameters, activation energy (E a) and Arrhenius parameters (A) for the samples containing 0.5 to 1% of PET-R, were calculated using the Flynn-Wall-Ozawa isoconversional method. It was observed that for greater amounts of PET-R added, there is a decrease in the E a values for the cure process. A Kinetic compensation effect (KCE), represented by the equation InA=-2.70+0.31E a was observed for all the samples. The most suitable kinetic model to describe this cure process is the autocatalytic Šesták-Berggreen, model applied to heterogeneous systems. © 2007 Springer Science+Business Media, LLC.

Reaction Control in Quiescent Systems of Free-Radical Retrograde-Precipitation Polymerization

Free-radical retrograde-precipitation polymerization, FRRPP in short, is a novel polymerization process discovered by Dr. Gerard Caneba in the late 1980s. The current study is aimed at gaining a better understanding of the reaction mechanism of the FRRPP and its thermodynamically-driven features that are predominant in controlling the chain reaction. A previously developed mathematical model to represent free radical polymerization kinetics was used to simulate a classic bulk polymerization system from the literature. Unlike other existing models, such a sparse-matrix-based representation allows one to explicitly accommodate the chain length dependent kinetic parameters. Extrapolating from the past results, mixing was experimentally shown to be exerting a significant influence on reaction control in FRRPP systems. Mixing alone drives the otherwise severely diffusion-controlled reaction propagation in phase-separated polymer domains. Therefore, in a quiescent system, in the absence of mixing, it is possible to retard the growth of phase-separated domains, thus producing isolated polymer nanoparticles (globules). Such a diffusion-controlled, self-limiting phenomenon of chain growth was also observed using time-resolved small angle x-ray scattering studies of reaction kinetics in quiescent systems of FRRPP. Combining the concept of self-limiting chain growth in quiescent FRRPP systems with spatioselective reaction initiation of lithography, microgel structures were synthesized in a single step, without the use of molds or additives. Hard x-rays from the bending magnet radiation of a synchrotron were used as an initiation source, instead of the more statistally-oriented chemical initiators. Such a spatially-defined reaction was shown to be self-limiting to the irradiated regions following a polymerization-induced self-assembly phenomenon. The pattern transfer aspects of this technique were, therefore, studied in the FRRP polymerization of N-isopropylacrylamide (NIPAm) and methacrylic acid (MAA), a thermoreversible and ionic hydrogel, respectively. Reaction temperature increases the contrast between the exposed and unexposed zones of the formed microgels, while the irradiation dose is directly proportional to the extent of phase separation. The response of Poly (NIPAm) microgels prepared from the technique described in this study was also characterized by small angle neutron scattering.

Reaction of the Escherichia coli quinol oxidase cytochrome bo3 with dioxygen: the role of a bound ubiquinone molecule.

We have studied the kinetics of the oxygen reaction of the fully reduced quinol oxidase, cytochrome bo3, using flow-flash and stopped flow techniques. This enzyme belongs to the heme-copper oxidase family but lacks the CuA center of the cytochrome c oxidases. Depending on the isolation procedure, the kinetics are found to be either nearly monophasic and very different from those of cytochrome c oxidase or multiphasic and quite similar to cytochrome c oxidase. The multiphasic kinetics in cytochrome c oxidase can largely be attributed to the presence Of CuA as the donor of a fourth electron, which rereduces the originally oxidized low-spin heme and completes the reduction of O2 to water. Monophasic kinetics would thus be expected, a priori, for cytochrome bo3 since it lacks the CuA center, and in this case we show that the oxygen reaction is incomplete and ends with the ferryl intermediate. Multiphasic kinetics thus suggest the presence of an extra electron donor (analogous to CuA). We observe such kinetics exclusively with cytochrome bo3 that contains a single equivalent of bound ubiquinone-8, whereas we find no bound ubiquinone in an enzyme exhibiting monophasic kinetics. Reconstitution with ubiquinone-8 converts the reaction kinetics from monophasic to multiphasic. We conclude that a single bound ubiquinone molecule in cytochrome bo3 is capable of fast rereduction of heme b and that the reaction with O2 is quite similar in quinol and cytochrome c oxidases.

A Lanczos-powered implementation of the Faber polynomial quantum time propagator for reaction probabilities

Complementing our recent work on subspace wavepacket propagation [Chem. Phys. Lett. 336 (2001) 149], we introduce a Lanczos-based implementation of the Faber polynomial quantum long-time propagator. The original version [J. Chem. Phys. 101 (1994) 10493] implicitly handles non-Hermitian Hamiltonians, that is, those perturbed by imaginary absorbing potentials to handle unwanted reflection effects. However, like many wavepacket propagation schemes, it encounters a bottleneck associated with dense matrix-vector multiplications. Our implementation seeks to reduce the quantity of such costly operations without sacrificing numerical accuracy. For some benchmark scattering problems, our approach compares favourably with the original. (C) 2004 Elsevier B.V. All rights reserved.

Structural evolution in a hydrothermal reaction between Nb2O5 and NaOH solution: From Nb2O5 grains to microporous Na2Nb2O6 center dot(2)/3H2O fibers and NaNbO3 cubes

Niobium pentoxide reacts actively with concentrate NaOH solution under hydrothermal conditions at as low as 120 degrees C. The reaction ruptures the corner-sharing of NbO7 decahedra and NbO6 octahedra in the reactant Nb2O5, yielding various niobates, and the structure and composition of the niobates depend on the reaction temperature and time. The morphological evolution of the solid products in the reaction at 180 degrees C is monitored via SEM: the fine Nb2O5 powder aggregates first to irregular bars, and then niobate fibers with an aspect ratio of hundreds form. The fibers are microporous molecular sieve with a monoclinic lattice, Na2Nb2O6 center dot(2)/3H2O. The fibers are a metastable intermediate of this reaction, and they completely convert to the final product NaNbO3 Cubes in the prolonged reaction of 1 h. This study demonstrates that by carefully optimizing the reaction condition, we can selectively fabricate niobate structures of high purity, including the delicate microporous fibers, through a direct reaction between concentrated NaOH solution and Nb2O5. This synthesis route is simple and suitable for the large-scale production of the fibers. The reaction first yields poorly crystallized niobates consisting of edge-sharing NbO6 octahedra, and then the microporous fibers crystallize and grow by assembling NbO6 octahedra or clusters of NbO6 octahedra and NaO6 units. Thus, the selection of the fibril or cubic product is achieved by control of reaction kinetics. Finally, niobates with different structures exhibit remarkable differences in light absorption and photoluminescence properties. Therefore, this study is of importance for developing new functional materials by the wet-chemistry process.

Mass transfer and chemical reaction on a distillation plate

A multistage distillation column in which mass transfer and a reversible chemical reaction occurred simultaneously, has been investigated to formulate a technique by which this process can be analysed or predicted. A transesterification reaction between ethyl alcohol and butyl acetate, catalysed by concentrated sulphuric acid, was selected for the investigation and all the components were analysed on a gas liquid chromatograph. The transesterification reaction kinetics have been studied in a batch reactor for catalyst concentrations of 0.1 - 1.0 weight percent and temperatures between 21.4 and 85.0 °C. The reaction was found to be second order and dependent on the catalyst concentration at a given temperature. The vapour liquid equilibrium data for six binary, four ternary and one quaternary systems are measured at atmospheric pressure using a modified Cathala dynamic equilibrium still. The systems with the exception of ethyl alcohol - butyl alcohol mixtures, were found to be non-ideal. Multicomponent vapour liquid equilibrium compositions were predicted by a computer programme which utilised the Van Laar constants obtained from the binary data sets. Good agreement was obtained between the predicted and experimental quaternary equilibrium vapour compositions. Continuous transesterification experiments were carried out in a six stage sieve plate distillation column. The column was 3" in internal diameter and of unit construction in glass. The plates were 8" apart and had a free area of 7.7%. Both the liquid and vapour streams were analysed. The component conversion was dependent on the boilup rate and the reflux ratio. Because of the presence of the reaction, the concentration of one of the lighter components increased below the feed plate. In the same region a highly developed foam was formed due to the presence of the catalyst. The experimental results were analysed by the solution of a series of simultaneous enthalpy and mass equations. Good agreement was obtained between the experimental and calculated results.