A master equation model for bimolecular reaction via multi-well isomerization intermediates

A reversible linear master equation model is presented for pressure- and temperature-dependent bimolecular reactions proceeding via multiple long-lived intermediates. This kinetic treatment, which applies when the reactions are measured under pseudo-first-order conditions, facilitates accurate and efficient simulation of the time dependence of the populations of reactants, intermediate species and products. Detailed exploratory calculations have been carried out to demonstrate the capabilities of the approach, with applications to the bimolecular association reaction C3H6 + H reversible arrow C3H7 and the bimolecular chemical activation reaction C2H2 +(CH2)-C-1--> C3H3+H. The efficiency of the method can be dramatically enhanced through use of a diffusion approximation to the master equation, and a methodology for exploiting the sparse structure of the resulting rate matrix is established.

The influence of collision energy and strain accumulation on the kinetics of mechanical alloying

The kinetics of mechanical alloying have been investigated by examining the effect that ball mass has on the rate at which titanium carbide forms from the elements. By varying the ball density while keeping the ball diameter and the charge ratio constant, the collision energy was independently controlled. Grinding media with a density from 3.8 g cm(-3) (agate) to 16.4 g cm(-3) (tungsten carbide) were used. The reaction rate increases exponentially with ball mass until a critical level is reached, which is determined by the induced temperature rise. Above this level, collisions of higher energy have no advantage. It is also shown that the reaction rate increases exponentially with the rate at which strain accumulates in the reactants. It is suggested that the strain accumulation rate in mechanically induced reactions is analogous to temperature in thermally induced chemical reactions.

Quantitative affinity chromatography

The objective of this review is to summarize developments in the use of quantitative affinity chromatography to determine equilibrium constants for solute interactions of biological interest. Affinity chromatography is an extremely versatile method for characterizing interactions between dissimilar reactants because the biospecificity incorporated into the design of the affinity matrix ensures applicability of the method regardless of the relative sizes of the two reacting solutes. Adoption of different experimental strategies, such as column chromatography, simple partition equilibrium experiments, solid-phase immunoassay, and biosensor technology, has led to a situation whereby affinity chromatography affords a means of characterizing interactions governed by an extremely broad range of binding affinities-relatively weak interactions (binding constants below 10(3) M-1) through to interactions with binding constants in excess of 10(9) M-1. In addition to its important role in solute separation and purification, affinity chromatography thus also possesses considerable potential for investigating the functional roles of the reactants thereby purified. (C) 2001 Elsevier Science B.V. All rights reserved.

Interpreting the effects of small uncharged solutes on protein-folding equilibria

Proteins are designed to function in environments crowded by cosolutes, but most studies of protein equilibria are conducted in dilute solution. While there is no doubt that crowding changes protein equilibria, interpretations of the changes remain controversial. This review combines experimental observations on the effect of small uncharged cosolutes (mostly sugars) on protein stability with a discussion of the thermodynamics of cosolute-induced nonideality and critical assessments of the most commonly applied interpretations. Despite the controversy surrounding the most appropriate manner for interpreting these effects of thermodynamic nonideality arising from the presence of small cosolutes, experimental advantage may still be taken of the ability of the cosolute effect to promote not only protein stabilization but also protein self-association and complex formation between dissimilar reactants. This phenomenon clearly has potential ramifications in the cell, where the crowded environment could well induce the same effects.

The mechanisms of combustion and continuous reactions during mechanical alloying

Some materials exhibit a combustion event during mechanical alloying, which results in the rapid transformation of reactants into products, while others show a slow transformation of reactants into products, In this paper, the continuous W + C --> WC reaction is compared to the Ti + C --> TiC combustion reaction. Rietveld refinement of X-ray diffraction patterns is used to show that these particular reactions proceed through different pathways, determined by crystallographic factors of the reactants. When a crystallographic relationship exists between the reactants and the products, such as that between W and WC, the product forms slowly over a period of time. In contrast, insertion of C into the Ti structure is associated with atomic rearrangements within the crowded lattice planes and the subsequent catastrophic failure of the reactant lattices results in combustion to form TiC. (C) 2001 Academic Press.

Automatic learning for the classification of chemical reactions and in statistical thermodynamics

This Thesis describes the application of automatic learning methods for a) the classification of organic and metabolic reactions, and b) the mapping of Potential Energy Surfaces(PES). The classification of reactions was approached with two distinct methodologies: a representation of chemical reactions based on NMR data, and a representation of chemical reactions from the reaction equation based on the physico-chemical and topological features of chemical bonds. NMR-based classification of photochemical and enzymatic reactions. Photochemical and metabolic reactions were classified by Kohonen Self-Organizing Maps (Kohonen SOMs) and Random Forests (RFs) taking as input the difference between the 1H NMR spectra of the products and the reactants. The development of such a representation can be applied in automatic analysis of changes in the 1H NMR spectrum of a mixture and their interpretation in terms of the chemical reactions taking place. Examples of possible applications are the monitoring of reaction processes, evaluation of the stability of chemicals, or even the interpretation of metabonomic data. A Kohonen SOM trained with a data set of metabolic reactions catalysed by transferases was able to correctly classify 75% of an independent test set in terms of the EC number subclass. Random Forests improved the correct predictions to 79%. With photochemical reactions classified into 7 groups, an independent test set was classified with 86-93% accuracy. The data set of photochemical reactions was also used to simulate mixtures with two reactions occurring simultaneously. Kohonen SOMs and Feed-Forward Neural Networks (FFNNs) were trained to classify the reactions occurring in a mixture based on the 1H NMR spectra of the products and reactants. Kohonen SOMs allowed the correct assignment of 53-63% of the mixtures (in a test set). Counter-Propagation Neural Networks (CPNNs) gave origin to similar results. The use of supervised learning techniques allowed an improvement in the results. They were improved to 77% of correct assignments when an ensemble of ten FFNNs were used and to 80% when Random Forests were used. This study was performed with NMR data simulated from the molecular structure by the SPINUS program. In the design of one test set, simulated data was combined with experimental data. The results support the proposal of linking databases of chemical reactions to experimental or simulated NMR data for automatic classification of reactions and mixtures of reactions. Genome-scale classification of enzymatic reactions from their reaction equation. The MOLMAP descriptor relies on a Kohonen SOM that defines types of bonds on the basis of their physico-chemical and topological properties. The MOLMAP descriptor of a molecule represents the types of bonds available in that molecule. The MOLMAP descriptor of a reaction is defined as the difference between the MOLMAPs of the products and the reactants, and numerically encodes the pattern of bonds that are broken, changed, and made during a chemical reaction. The automatic perception of chemical similarities between metabolic reactions is required for a variety of applications ranging from the computer validation of classification systems, genome-scale reconstruction (or comparison) of metabolic pathways, to the classification of enzymatic mechanisms. Catalytic functions of proteins are generally described by the EC numbers that are simultaneously employed as identifiers of reactions, enzymes, and enzyme genes, thus linking metabolic and genomic information. Different methods should be available to automatically compare metabolic reactions and for the automatic assignment of EC numbers to reactions still not officially classified. In this study, the genome-scale data set of enzymatic reactions available in the KEGG database was encoded by the MOLMAP descriptors, and was submitted to Kohonen SOMs to compare the resulting map with the official EC number classification, to explore the possibility of predicting EC numbers from the reaction equation, and to assess the internal consistency of the EC classification at the class level. A general agreement with the EC classification was observed, i.e. a relationship between the similarity of MOLMAPs and the similarity of EC numbers. At the same time, MOLMAPs were able to discriminate between EC sub-subclasses. EC numbers could be assigned at the class, subclass, and sub-subclass levels with accuracies up to 92%, 80%, and 70% for independent test sets. The correspondence between chemical similarity of metabolic reactions and their MOLMAP descriptors was applied to the identification of a number of reactions mapped into the same neuron but belonging to different EC classes, which demonstrated the ability of the MOLMAP/SOM approach to verify the internal consistency of classifications in databases of metabolic reactions. RFs were also used to assign the four levels of the EC hierarchy from the reaction equation. EC numbers were correctly assigned in 95%, 90%, 85% and 86% of the cases (for independent test sets) at the class, subclass, sub-subclass and full EC number level,respectively. Experiments for the classification of reactions from the main reactants and products were performed with RFs - EC numbers were assigned at the class, subclass and sub-subclass level with accuracies of 78%, 74% and 63%, respectively. In the course of the experiments with metabolic reactions we suggested that the MOLMAP / SOM concept could be extended to the representation of other levels of metabolic information such as metabolic pathways. Following the MOLMAP idea, the pattern of neurons activated by the reactions of a metabolic pathway is a representation of the reactions involved in that pathway - a descriptor of the metabolic pathway. This reasoning enabled the comparison of different pathways, the automatic classification of pathways, and a classification of organisms based on their biochemical machinery. The three levels of classification (from bonds to metabolic pathways) allowed to map and perceive chemical similarities between metabolic pathways even for pathways of different types of metabolism and pathways that do not share similarities in terms of EC numbers. Mapping of PES by neural networks (NNs). In a first series of experiments, ensembles of Feed-Forward NNs (EnsFFNNs) and Associative Neural Networks (ASNNs) were trained to reproduce PES represented by the Lennard-Jones (LJ) analytical potential function. The accuracy of the method was assessed by comparing the results of molecular dynamics simulations (thermal, structural, and dynamic properties) obtained from the NNs-PES and from the LJ function. The results indicated that for LJ-type potentials, NNs can be trained to generate accurate PES to be used in molecular simulations. EnsFFNNs and ASNNs gave better results than single FFNNs. A remarkable ability of the NNs models to interpolate between distant curves and accurately reproduce potentials to be used in molecular simulations is shown. The purpose of the first study was to systematically analyse the accuracy of different NNs. Our main motivation, however, is reflected in the next study: the mapping of multidimensional PES by NNs to simulate, by Molecular Dynamics or Monte Carlo, the adsorption and self-assembly of solvated organic molecules on noble-metal electrodes. Indeed, for such complex and heterogeneous systems the development of suitable analytical functions that fit quantum mechanical interaction energies is a non-trivial or even impossible task. The data consisted of energy values, from Density Functional Theory (DFT) calculations, at different distances, for several molecular orientations and three electrode adsorption sites. The results indicate that NNs require a data set large enough to cover well the diversity of possible interaction sites, distances, and orientations. NNs trained with such data sets can perform equally well or even better than analytical functions. Therefore, they can be used in molecular simulations, particularly for the ethanol/Au (111) interface which is the case studied in the present Thesis. Once properly trained, the networks are able to produce, as output, any required number of energy points for accurate interpolations.

New RuII(arene) complexes with halogen-substituted bis- and tris(pyrazol-1-YL)borate ligands

[RuCl(arene)(-Cl)](2) dimers were treated in a 1:2 molar ratio with sodium or thallium salts of bis- and tris(pyrazolyl)borate ligands [Na(BpBr3)], [Tl(TpBr3)], and [Tl(Tp(iPr,4Br))]. Mononuclear neutral complexes [RuCl(arene)((2)-BpBr3)] (1: arene=p-cymene (cym); 2: arene=hexamethylbenzene (hmb); 3: arene=benzene (bz)), [RuCl(arene)((2)-TpBr3)] (4: arene=cym; 6: arene=bz), and [RuCl(arene)((2)-Tp(iPr,4Br))] (7: arene=cym, 8: arene=hmb, 9: arene=bz) have been always obtained with the exception of the ionic [Ru-2(hmb)(2)(-Cl)(3)][TpBr3] (5), which formed independently of the ratio of reactants and reaction conditions employed. The ionic [Ru(CH3OH)(cym)((2)-BpBr3)][X] (10: X=PF6, 12: X=O3SCF3) and the neutral [Ru(O2CCF3)(cym)((2)-BpBr3)] (11) have been obtained by a metathesis reaction with corresponding silver salts. All complexes 1-12 have been characterized by analytical and spectroscopic data (IR, ESI-MS, H-1 and (CNMR)-C-13 spectroscopy). The structures of the thallium and calcium derivatives of ligand TpBr3, [Tl(TpBr3)] and [Ca(dmso)(6)][TpBr3](2)2DMSO, of the complexes 1, 4, 5, 6, 11, and of the decomposition product [RuCl(cym)(Hpz(iPr,4Br))(2)][Cl] (7) have been confirmed by using single-crystal X-ray diffraction. Electrochemical studies showed that 1-9 and 11 undergo a single-electron (RuRuIII)-Ru-II oxidation at a potential, measured by cyclic voltammetry, which allows comparison of the electron-donor characters of the bis- and tris(pyrazol-1-yl)borate and arene ligands, and to estimate, for the first time, the values of the Lever E-L ligand parameter for BpBr3, TpBr3, and Tp(iPr,4Br). Theoretical calculations at the DFT level indicated that both oxidation and reduction of the Ru complexes under study are mostly metal-centered with some involvement of the chloride ligand in the former case, and also demonstrated that the experimental isolation of the (3)-binuclear complex 5 (instead of the mononuclear 5) is accounted for by the low thermodynamic stability of the latter species due to steric reasons.

Projeto de uma unidade de produção de biodiesel num reator ultrassónico pela via supercrítica

A concretização deste estágio na INCBIO teve como principal objetivo o projeto e dimensionamento de uma unidade de produção de biodiesel por transesterificação nãocatalítica num reator ultrassónico com condições supercríticas de pressão e temperatura. Com vista à concretização do trabalho, iniciou-se a realização do estudo do estado da arte relativo à produção de biodiesel, com particular enfoque na produção de biodiesel por via catalítica e por via supercrítica e na produção de biodiesel com uso de tecnologia ultrassónica. Conclui-se que nenhum estudo contempla a combinação simultânea da produção de biodiesel por via supercrítica através da tecnologia ultrassónica. Este estudo do estado da arte permitiu ainda definir as condições de temperatura, pressão e rácio mássico (250 °C, 95 bar e 1:1 respetivamente) a considerar no projeto da unidade de produção de biodiesel deste trabalho. Com base no estudo do estado da arte efetuado e com base nas características da matéria procedeu-se à definição do processo de produção de biodiesel. Para a definição do processo começou-se por elaborar o diagrama de blocos do processo (BFD) e o diagrama de fluxo do processo (PFD). Com base nos diagramas e na composição da matéria-prima, procedeu-se à quantificação dos reagentes (metanol) com base na estequiometria das reações envolvidas e ao cálculo do balanço de massa. O balanço de massa foi calculado com base na estequiometria das reações envolvidas e foi também calculado através do software de simulação ASPEN PLUS. Após o cálculo do balanço de massa elaborou-se o diagrama de tubulação e instrumentação (P&ID), que contém todos os equipamentos, válvulas, instrumentação e tubagens existentes na unidade. Após a definição do processo e cálculo do balanço de massa procedeu-se ao dimensionamento mecânico e cálculo hidráulico dos tanques, tubagem, bombas, permutador de calor, reator ultrassónico, válvulas de controlo e instrumentação de acordo com as normas ASME. Nesta fase do trabalho foram consultados diversos fornecedores possíveis para a compra de todo o material necessário. O dimensionamento mecânico e cálculo hidráulico efetuados permitiram, entre outras informações relevantes, obter as dimensões necessárias à construção do layout e à elaboração do desenho 3D. Com os resultados obtidos e desenhos elaborados, é possível avançar com a construção da unidade, pelo que pode-se inferir que o objetivo de projetar uma unidade de produção de biodiesel num reator ultrassónico com condições supercríticas foi alcançado. Por fim, efetuou-se uma análise económica detalhada que possibilita a comparação de uma unidade de produção de biodiesel por via catalítica (unidade de produção da INCBIO) com a unidade de produção de biodiesel num reator ultrassónico com condições supercríticas (unidade projetada neste trabalho). Por questões de confidencialidade, não foram revelados maior parte dos resultados da análise económica efetuada. No entanto, conclui-se que os custos de construção da unidade de produção de biodiesel num reator ultrassónico com condições supercríticas são mais baixos cerca de 35 a 40%, quando comparados com os custos de construção da unidade de produção de biodiesel por via catalítica, evidenciando assim que a combinação em simultâneo das condições supercríticas com a tecnologia ultrassónica possibilita a diminuição dos custos de produção.

Valorização da Borra de Café: Otimização da Produção de Biodiesel por Catálise Enzimática

O decréscimo das reservas de petróleo e as consequências ambientais resultantes do recurso a combustíveis fósseis nos motores a diesel têm levado à procura de combustíveis alternativos. Esta pesquisa alicerçada nas fontes de energia renovável tornou-se essencial, face à crescente procura de energia e ao limitado fornecimento de combustíveis fósseis . Resíduos de óleo de cozinha, gordura animal, entre outros resíduos de origem biológica, tais como a borra de café, são exemplos de matérias-primas para a produção de biodiesel. A sua valorização tem interesse quer pela perspetiva ambiental, quer pela económica, pois aumenta não só a flexibilidade e diversificação das matérias-primas, mas também contribui para uma estabilidade de custos e alteração nas políticas agrícolas e de uso do solo. É neste contexto que se enquadra o biodiesel e a borra de café, pretendendo-se aqui efetuar o estudo da produção, à escala laboratorial, de biodiesel a partir da borra de café, por transesterificação enzimática, visando a procura das melhores condições reacionais. Iniciando-se com a caracterização da borra de café, foram avaliados antes e após a extração do óleo da borra de café, diversos parâmetros, de entre os quais se destacam: o teor de humidade (16,97% e 6,79%), teor de cinzas (1,91 e 1,57%), teor de azoto (1,71 e 2,30%), teor de proteínas (10,7 e 14,4%), teor de carbono (70,2 e 71,7%), teor de celulose bruta (14,77 e 18,48%), teor de lenhina (31,03% e 30,97%) e poder calorifico superior (19,5 MJ/kg e 19,9 MJ/kg). Sumariamente, constatou-se que os valores da maioria dos parâmetros não difere substancialmente dos valores encontrados na literatura, tendo sido evidenciado o potencial da utilização desta biomassa, como fonte calorifica para queima e geração de energia. Sendo a caracterização do óleo extraído da borra de café um dos objetivos antecedentes à produção do biodiesel, pretendeu-se avaliar os diferentes parâmetros mais significativos. No que diz respeito à caracterização do óleo extraído, distingue-se a sua viscosidade cinemática (38,04 mm2/s), densidade 0,9032 g/cm3, poder calorífico de 37,9 kcal/kg, índice de iodo igual a 63,0 gI2/ 100 g óleo, o teor de água do óleo foi de 0,15 %, o índice de acidez igual a 44,8 mg KOH/g óleo, ponto de inflamação superior a 120 ºC e teor em ácidos gordos de 82,8%. Inicialmente foram efetuados ensaios preliminares, a fim de selecionar a lipase (Lipase RMIM, TL 100L e CALB L) e álcool (metanol ou etanol puros) mais adequados à produção de biodiesel, pelo que o rendimento de 83,5% foi obtido através da transesterificação mediada pela lipase RMIM, utilizando como álcool o etanol. Sendo outro dos objetivos a otimização do processo de transesterificação enzimática, através de um desenho composto central a três variáveis (razão molar etanol: óleo, concentração de enzima e temperatura), recorrendo ao software JMP 8.0, determinou-se como melhores condições, uma razão molar etanol: óleo 5:1, adição de 4,5% (m/m) de enzima e uma temperatura de 45 ºC, que conduziram a um rendimento experimental equivalente a 96,7 % e teor de ésteres 87,6%. Nestas condições, o rendimento teórico foi de 99,98%. Procurou-se ainda estudar o efeito da adição de água ao etanol, isto é, o efeito da variação da concentração do etanol pela adição de água, para teores de etanol de 92%, 85% e 75%. Verificou-se que até 92% decorreu um aumento da transesterificação (97,2%) para um teor de ésteres de (92,2%), pelo que para teores superiores de água adicionada (75% e 85%) ocorreu um decréscimo no teor final em ésteres (77,2% e 89,9%) e no rendimento da reação (84,3% e 91,9%). Isto indica a ocorrência da reação de hidrólise em maior extensão, que leva ao desvio do equilíbrio no sentido contrário à reação de formação dos produtos, isto é, dos ésteres. Finalmente, relativamente aos custos associados ao processo de produção de biodiesel, foram estimados para o conjunto de 27 ensaios realizados neste trabalho, e que corresponderam a 767,4 g de biodiesel produzido, sendo o custo dos reagentes superior ao custo energético, de 156,16 € e 126,02 €, respetivamente. Naturalmente que não esperamos que, a nível industrial os custos sejam desta ordem de grandeza, tanto mais que há economia de escala e que as enzimas utilizadas no processo deveriam ser reutilizadas diversas vezes.

Development of integrated separation processes with green solvents

Dissertação para obtenção do Grau de Doutor em Química Sustentável

Optic perineuritis secondary to Wegener's granulomatosis.

BACKGROUND: Optic perineuritis (OPN) is an inflammatory condition involving the optic nerve sheath because of a variety of causes. We describe three patients in whom OPN was secondary to Wegener's granulomatosis (WG) and compare the clinical findings in these cases with those of idiopathic OPN. METHODS: This is a retrospective small case series derived from patients with OPN seen in an outpatient neuro-ophthalmology clinic. Medical records and imaging studies of these patients were reviewed. RESULTS: These patients shared clinical similarities with idiopathic OPN including age, sex, symptoms, radiographic findings and steroid responsiveness. However, recurrence of symptoms on lowering the prednisone dose below 40 mg distinguished these patients from those with idiopathic OPN. CONCLUSIONS: Steroid dependency in idiopathic OPN should raise suspicion of WG. Patients with OPN should be specifically questioned regarding pre-existing upper respiratory tract disorders and rheumatic symptoms and laboratory testing should include acute phase reactants, anti-neutrophil cytoplasmic antibodies and tests of renal function.

Bioconjugated lanthanide luminescent helicates as multilabels for lab-on-a-chip detection of cancer biomarkers.

The lanthanide binuclear helicate [Eu(2)(L(C2(CO(2)H)))(3)] is coupled to avidin to yield a luminescent bioconjugate EuB1 (Q = 9.3%, tau((5)D(0)) = 2.17 ms). MALDI/TOF mass spectrometry confirms the covalent binding of the Eu chelate and UV-visible spectroscopy allows one to determine a luminophore/protein ratio equal to 3.2. Bio-affinity assays involving the recognition of a mucin-like protein expressed on human breast cancer MCF-7 cells by a biotinylated monoclonal antibody 5D10 to which EuB1 is attached via avidin-biotin coupling demonstrate that (i) avidin activity is little affected by the coupling reaction and (ii) detection limits obtained by time-resolved (TR) luminescence with EuB1 and a commercial Eu-avidin conjugate are one order of magnitude lower than those of an organic conjugate (FITC-streptavidin). In the second part of the paper, conditions for growing MCF-7 cells in 100-200 microm wide microchannels engraved in PDMS are established; we demonstrate that EuB1 can be applied as effectively on this lab-on-a-chip device for the detection of tumour-associated antigens as on MCF-7 cells grown in normal culture vials. In order to exploit the versatility of the ligand used for self-assembling [Ln(2)(L(C2(CO(2)H)))(3)] helicates, which sensitizes the luminescence of both Eu(III) and Tb(III) ions, a dual on-chip assay is proposed in which estrogen receptors (ERs) and human epidermal growth factor receptors (Her2/neu) can be simultaneously detected on human breast cancer tissue sections. The Ln helicates are coupled to two secondary antibodies: ERs are visualized by red-emitting EuB4 using goat anti-mouse IgG and Her2/neu receptors by green-emitting TbB5 using goat anti-rabbit IgG. The fact that the assay is more than 6 times faster and requires 5 times less reactants than conventional immunohistochemical assays provides essential advantages over conventional immunohistochemistry for future clinical biomarker detection.

Modelling changes in stable isotope compositions of minerals during net transfer reactions in a contact aureole: Wollastonite growth at the northern Hunter Mountain Batholith (Death Valley National Park, USA)

One of the world's largest wollastonite deposits was formed at the contact of the northern Hunter Mountain Batholith (California, USA) in Paleozoic sediments. Wollastonite occurs as zones of variable thickness surrounding layers or nodules of quartzite in limestones. A minimum formation temperature of 650 degrees C is estimated from isolated periclase-bearing lenses in that area. Contact metamorphism of siliceous carbonates has produced mineral assemblages that are consistent with heterogeneous, and partly limited infiltration of water-rich fluids, compatible with O-18/O-16 and C-13/C-12 isotopic patterns recorded in carbonates. Oxygen isotope compositions of wollastonites in the study area may also not require infiltration of large quantities of externally-derived fluids that were out of equilibrium with the rocks. 8180 values of wollastonite are high (14.8 parts per thousand to 25.0 parts per thousand; median: 19.7 parts per thousand) and close to those of the host limestone (19.7 parts per thousand to 28 parts per thousand; median: 24.9 parts per thousand) and quartz (18.0 parts per thousand. to 29.1 parts per thousand; median: 22.6 parts per thousand). Isotopic disequilibrium exists at quartz/wollastonite and wollastonite/calcite boundaries. Therefore, classical batch/Rayleigh fractionation models based on reactant and product equilibrium are not applicable to the wollastonite rims. An approach that relies on local instantaneous mass balance for the reactants, based on the wollastonite-forming reaction is suggested as an alternative way to model wollastonite reaction rims. This model reproduces many of the measured delta O-18 values of wollastonite reaction rims of the current study to within +/- 1 parts per thousand, even though the wollastonite compositions vary by almost 10 parts per thousand. (C) 2011 Elsevier B.V. All rights reserved.

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Spatial correlations and cross sections of clusters in the A + B -> 0 reaction

We consider the distribution of cross sections of clusters and the density-density correlation functions for the A+B¿0 reaction. We solve the reaction-diffusion equations numerically for random initial distributions of reactants. When both reactant species have the same diffusion coefficients the distribution of cross sections and the correlation functions scale with the diffusion length and obey superuniversal laws (independent of dimension). For different diffusion coefficients the correlation functions still scale, but the scaling functions depend on the dimension and on the diffusion coefficients. Furthermore, we display explicitly the peculiarities of the cluster-size distribution in one dimension.