Parallelization of kinetic Monte Carlo algorithm to simulate AL3Sc precipitation

The present paper reports the precipitation process of Al3Sc structures in an aluminum scandium alloy, which has been simulated with a synchronous parallel kinetic Monte Carlo (spkMC) algorithm. The spkMC implementation is based on the vacancy diffusion mechanism. To filter the raw data generated by the spkMC simulations, the density-based clustering with noise (DBSCAN) method has been employed. spkMC and DBSCAN algorithms were implemented in the C language and using MPI library. The simulations were conducted in the SeARCH cluster located at the University of Minho. The Al3Sc precipitation was successfully simulated at the atomistic scale with the spkMC. DBSCAN proved to be a valuable aid to identify the precipitates by performing a cluster analysis of the simulation results. The achieved simulations results are in good agreement with those reported in the literature under sequential kinetic Monte Carlo simulations (kMC). The parallel implementation of kMC has provided a 4x speedup over the sequential version.

Extracting kinetic information from literature with KineticRE

To better understand the dynamic behavior of metabolic networks in a wide variety of conditions, the field of Systems Biology has increased its interest in the use of kinetic models. The different databases, available these days, do not contain enough data regarding this topic. Given that a significant part of the relevant information for the development of such models is still wide spread in the literature, it becomes essential to develop specific and powerful text mining tools to collect these data. In this context, this work has as main objective the development of a text mining tool to extract, from scientific literature, kinetic parameters, their respective values and their relations with enzymes and metabolites. The approach proposed integrates the development of a novel plug-in over the text mining framework @Note2. In the end, the pipeline developed was validated with a case study on Kluyveromyces lactis, spanning the analysis and results of 20 full text documents.

Arrested kinetic Li isotope fractionation at the margin of the llimaussaq complex, South Greenland: Evidence for open-system processes during final cooling of peralkaline igneous rocks

Li contents [Li] and isotopic composition (delta Li-7) of mafic minerals (mainly amphibole and clinopyroxene) from the alkaline to peralkaline Ilimaussaq plutonic complex, South Greenland, track the behavior of Li and its isotopes during magmatic differentiation and final cooling of an alkaline igneous system. [Li] in amphibole increase from < 10 ppm in Caamphiboles of the least differentiated unit to >3000 ppm in Na-amphiboles of the highly evolved units. In contrast, [Li] in clinopyroxene are comparatively low (<85 ppm) and do not vary systematically with differentiation. The distribution of Li between amphibole and pyroxene is controlled by the major element composition of the minerals (Ca-rich and Na-rich, respectively) and changes in oxygen fugacity (due to Li incorporation via coupled substitution with ferric iron) during magmatic differentiation. delta(7) Li values of all minerals span a wide range from + 17 to - 8 parts per thousand, with the different intrusive units of the complex having distinct Li isotopic systematics. Amphiboles, which dominate the Li budget of whole-rocks from the inner part of the complex, have constant delta Li-7 of + 1.8 +/- 2.2 parts per thousand (2 sigma, n = 15). This value reflects a homogeneous melt reservoir and is consistent with their mantle derivation, in agreement with published O and Nd isotopic data. Clinopyroxenes of these samples are consistently lighter, with Delta Li-7(amph-cpx). as large as 8 parts per thousand and are thus not in Li isotope equilibrium. These low values probably reflect late-stage diffusion of Li into clinopyroxene during final cooling of the rocks, thus enriching the clinopyroxene in 6 Li. At the margin of the complex delta(7) Li in the syenites increases systematically, from +2 to high values of + 14 parts per thousand. This, coupled with the observed Li isotope systematics of the granitic country rocks, reflects post-magmatic open-system processes occurring during final cooling of the intrusion. Although the shape and magnitude of the Li isotope and elemental profiles through syenite and country rock are suggestive of diffusion-driven isotope fractionation, they cannot be modeled by one-dimensional diffusive transport and point to circulation of a fluid having a high 67 Li value (possibly seawater) along the chilled contact. In all, this study demonstrates that Li isotopes can be used to identify complex fluid- and diffusion-governed processes taking place during the final cooling of such rocks. (c) 2007 Elsevier B.V All rights reserved.

Role of prohormone convertases in pro-neuropeptide Y processing: coexpression and in vitro kinetic investigations.

Proneuropeptide Y (ProNPY) undergoes cleavage at a single dibasic site Lys38-Arg39 resulting in the formation of 1-39 amino acid NPY which is further processed successively by carboxypeptidase-like and peptidylglycine alpha-amidating monooxygenase enzymes. To investigate whether prohormone convertases are involved in ProNPY processing, a vaccinia virus derived expression system was used to coexpress recombinant ProNPY with each of the prohormone convertases PC1/3, PC2, furin, and PACE4 in Neuro2A and NIH 3T3 cell lines as regulated neuroendocrine and constitutive prototype cell lines, respectively. The analysis of processed products shows that only PC1/3 generates NPY in NIH 3T3 cells while both PC1/3 and PC2 are able to generate NPY in Neuro2A cells. The convertases furin and PACE4 are unable to process ProNPY in either cell line. Moreover, comparative in vitro cleavage of recombinant NPY precursor by the enzymes PC1/3, PC2 and furin shows that only PC1/3 and PC2 are involved in specific cleavage of the dibasic site. Kinetic studies demonstrate that PC1/3 cleaves ProNPY more efficiently than PC2. The main difference between the cleavage efficiency is observed in the Vmax values whereas no major difference is observed in Km values. In addition the cleavage by PC1/3 and PC2 of two peptides reproducing the dibasic cleavage site with different amino acid sequence lengths namely (20-49)-ProNPY and (28-43)-ProNPY was studied. These shortened ProNPY substrates, when recognized by the enzymes, are more efficiently cleaved than ProNPY itself. The shortest peptide is not cleaved by PC2 while it is by PC1/3. On the basis of these observations it is proposed, first, that the constitutive secreted NPY does not result from the cleavage carried out by ubiquitously expressed enzymes furin and PACE4; second, that PC1/3 and PC2 are not equipotent in the cleavage of ProNPY; and third, substrate peptide length might discriminate PC1/3 and PC2 processing activity.

Nanoparticles in SiH4-Ar plasma: Modelling and comparison with experimental data

Experimental and theoretical investigations for growth of silicon nanoparticles (4 to 14 nm) in radio frequency discharge were carried out. Growth processes were performed with gas mixtures of SiH4 and Ar in a plasma chemical reactor at low pressure. A distinctive feature of presented kinetic model of generation and growth of nanoparticles (compared to our earlier model) is its ability to investigate small"critical" dimensions of clusters, determining the rate of particle production and taking into account the influence of SiH2 and Si2Hm dimer radicals. The experiments in the present study were extended to high pressure (≥20 Pa) and discharge power (≥40 W). Model calculations were compared to experimental measurements, investigating the dimension of silicon nanoparticles as a function of time, discharge power, gas mixture, total pressure, and gas flow.

Retrieval of kinetic rates in reactions with semi batch liquid phase using ill-posed inverse problem theory

A neural network procedure to solve inverse chemical kinetic problems is discussed in this work. Rate constants are calculated from the product concentration of an irreversible consecutive reaction: the hydrogenation of Citral molecule, a process with industrial interest. Simulated and experimental data are considered. Errors in the simulated data, up to 7% in the concentrations, were assumed to investigate the robustness of the inverse procedure. Also, the proposed method is compared with two common methods in nonlinear analysis; the Simplex and Levenberg-Marquardt approaches. In all situations investigated, the neural network approach was numerically stable and robust with respect to deviations in the initial conditions or experimental noises.

Methane combustion kinetic rate constants determination: an ill-posed inverse problem analysis

Methane combustion was studied by the Westbrook and Dryer model. This well-established simplified mechanism is very useful in combustion science, for computational effort can be notably reduced. In the inversion procedure to be studied, rate constants are obtained from [CO] concentration data. However, when inherent experimental errors in chemical concentrations are considered, an ill-conditioned inverse problem must be solved for which appropriate mathematical algorithms are needed. A recurrent neural network was chosen due to its numerical stability and robustness. The proposed methodology was compared against Simplex and Levenberg-Marquardt, the most used methods for optimization problems.

Kinetic modeling of the biosorption of Cd2+ ions from aqueous solutions onto Eichhornia crassipes roots using potentiometry: low-cost alternative to conventional methods

This work presents the use of potentiometric measurements for kinetic studies of biosorption of Cd2+ ions from aqueous solutions on Eichhornia crassipes roots. The open circuit potential of the Cd/Cd2+ electrode of the first kind was measured during the bioadsorption process. The amount of Cd2+ ions accumulated was determined in real time. The data were fit to different models, with the pseudo-second-order model proving to be the best in describing the data. The advantages and limitations of the methodology proposed relative to the traditional method are discussed.

Wet oxidation of recalcitrant lignin waters: Experimental and kinetic studies

The dissertation is based on four articles dealing with recalcitrant lignin water purification. Lignin, a complicated substance and recalcitrant to most treatment technologies, inhibits seriously pulp and paper industry waste management. Therefore, lignin is studied, using WO as a process method for its degradation. A special attention is paid to the improvement in biodegradability and the reduction of lignin content, since they have special importance for any following biological treatment. In most cases wet oxidation is not used as a complete ' mineralization method but as a pre treatment in order to eliminate toxic components and to reduce the high level of organics produced. The combination of wet oxidation with a biological treatment can be a good option due to its effectiveness and its relatively low technology cost. The literature part gives an overview of Advanced Oxidation Processes (AOPs). A hot oxidation process, wet oxidation (WO), is investigated in detail and is the AOP process used in the research. The background and main principles of wet oxidation, its industrial applications, the combination of wet oxidation with other water treatment technologies, principal reactions in WO, and key aspects of modelling and reaction kinetics are presented. There is also given a wood composition and lignin characterization (chemical composition, structure and origin), lignin containing waters, lignin degradation and reuse possibilities, and purification practices for lignin containing waters. The aim of the research was to investigate the effect of the operating conditions of WO, such as temperature, partial pressure of oxygen, pH and initial concentration of wastewater, on the efficiency, and to enhance the process and estimate optimal conditions for WO of recalcitrant lignin waters. Two different waters are studied (a lignin water model solution and debarking water from paper industry) to give as appropriate conditions as possible. Due to the great importance of re using and minimizing the residues of industries, further research is carried out using residual ash of an Estonian power plant as a catalyst in wet oxidation of lignin-containing water. Developing a kinetic model that includes in the prediction such parameters as TOC gives the opportunity to estimate the amount of emerging inorganic substances (degradation rate of waste) and not only the decrease of COD and BOD. The degradation target compound, lignin is included into the model through its COD value (CODligning). Such a kinetic model can be valuable in developing WO treatment processes for lignin containing waters, or other wastewaters containing one or more target compounds. In the first article, wet oxidation of "pure" lignin water was investigated as a model case with the aim of degrading lignin and enhancing water biodegradability. The experiments were performed at various temperatures (110 -190°C), partial oxygen pressures (0.5 -1.5 MPa) and pH (5, 9 and 12). The experiments showed that increasing the temperature notably improved the processes efficiency. 75% lignin reduction was detected at the lowest temperature tested and lignin removal improved to 100% at 190°C. The effect of temperature on the COD removal rate was lower, but clearly detectable. 53% of organics were oxidized at 190°C. The effect of pH occurred mostly on lignin removal. Increasing the pH enhanced the lignin removal efficiency from 60% to nearly 100%. A good biodegradability ratio (over 0.5) was generally achieved. The aim of the second article was to develop a mathematical model for "pure" lignin wet oxidation using lumped characteristics of water (COD, BOD, TOC) and lignin concentration. The model agreed well with the experimental data (R2 = 0.93 at pH 5 and 12) and concentration changes during wet oxidation followed adequately the experimental results. The model also showed correctly the trend of biodegradability (BOD/COD) changes. In the third article, the purpose of the research was to estimate optimal conditions for wet oxidation (WO) of debarking water from the paper industry. The WO experiments were' performed at various temperatures, partial oxygen pressures and pH. The experiments showed that lignin degradation and organics removal are affected remarkably by temperature and pH. 78-97% lignin reduction was detected at different WO conditions. Initial pH 12 caused faster removal of tannins/lignin content; but initial pH 5 was more effective for removal of total organics, represented by COD and TOC. Most of the decrease in organic substances concentrations occurred in the first 60 minutes. The aim of the fourth article was to compare the behaviour of two reaction kinetic models, based on experiments of wet oxidation of industrial debarking water under different conditions. The simpler model took into account only the changes in COD, BOD and TOC; the advanced model was similar to the model used in the second article. Comparing the results of the models, the second model was found to be more suitable for describing the kinetics of wet oxidation of debarking water. The significance of the reactions involved was compared on the basis of the model: for instance, lignin degraded first to other chemically oxidizable compounds rather than directly to biodegradable products. Catalytic wet oxidation of lignin containing waters is briefly presented at the end of the dissertation. Two completely different catalysts were used: a commercial Pt catalyst and waste power plant ash. CWO showed good performance using 1 g/L of residual ash gave lignin removal of 86% and COD removal of 39% at 150°C (a lower temperature and pressure than with WO). It was noted that the ash catalyst caused a remarkable removal rate for lignin degradation already during the pre heating for `zero' time, 58% of lignin was degraded. In general, wet oxidation is not recommended for use as a complete mineralization method, but as a pre treatment phase to eliminate toxic or difficultly biodegradable components and to reduce the high level of organics. Biological treatment is an appropriate post treatment method since easily biodegradable organic matter remains after the WO process. The combination of wet oxidation with subsequent biological treatment can be an effective option for the treatment of lignin containing waters.

Emission, kinetic and magnetic phenomena in rare-earth and transition metal doped ZnSe single crystals

In this work emission, optical, electrical and magnetic properties of the d- and f- elements doped zinc selenide crystals were investigated within a wide temperature range. Doping was performed in various technological processes: during the growth by chemical vapor transport method; by thermal diffusion from the Bi or Zn melt. Concentration of the doping impurity in the crystals was controlled by amount of the dopant in the source material or by its concentration in the doping media. Special interest in the work was paid to the influence of the different concentrations of Cr and Yb impurities on ZnSe crystals’ properties, correlations between observed effects and similarities with the Ni, Mn and Gd dopants are analysed. Possibility of formation of the excitons bound to the doping d-ions was shown. In contrast to this, it was observed that f-elements do not bound excitons, but prevent formation of excitons bound to some uncontrolled impurities. A mechanism of Cr doping impurity interaction with background impurities and zinc selenide structural defects was proposed based on experimental data. An assumption about resonant energy transfer between double charged chromium ions and complexes based on crystals’ vacancy defects was made. A correlation between emission and magnetic properties of the d- ions doped samples was established. Based on this correlation a mechanism explaining the concentration quench of the emission was proposed. It was found that f-ions bind electrically active shallow and deep donor and acceptor states of background impurity to electrically neutral complexes. This may be observed as “purification” of ZnSe crystals by doping with the rare-earth elements, resulting i tendency of the properties of f-ion doped crystals to the properties of intrinsic crystals, but with smaller concentration of uncontrolled native and impurity defects. A possible interpretation of this effect was proposed. It was shown that selenium substituting impurities decrease efficiency of the Yb doping. Based on this experimental results an attempt to determine ytterbium ion surroundings in the crystal lattice was made. It was shown that co-doping of zinc selenide crystals with the d- and f- ions leads to the combination of the impurities influence on the material’s properties. On the basis of obtained data an interaction mechanism of the d- and f-elements co-dopants was proposed. Guided by the model of the ytterbium ion incorporation in the selenide sublattice of the ZnSe crystals, an assumption about stabilization of single charged chromium ions in the zinc sublattice crystal nodes, by means of formation of the local charge compensating clusters, was made.

Modeling and kinetic study of bio-ethanol production from soy protein concentrate by-product

Abstract The aim of this work was to evaluate a non-agitated process of bioethanol production from soybean molasses and the kinetic parameters of fermentation using a strain of Saccharomyces cerevisiae (ATCC® 2345). Kinetic experiment was conducted in medium with 30% (w v-1) of soluble solids without supplementation or pH adjustment. The maximum ethanol concentration was in 44 hours, the ethanol productivity was 0.946 g L-1 h-1, the yield over total initial sugars (Y1) was 47.87%, over consumed sugars (Y2) was 88.08% and specific cells production rate was 0.006 h-1. The mathematical polynomial was adjusted to the experimental data and provided very similar parameters of yield and productivity. Based in this study, for one ton of soybean molasses can be produced 103 kg of anhydrous bioethanol.

The regulation of oat coleoptile phosphoenolpyruvate carboxylase and malic enzyme by Hâ ½ and metabolites : kinetic evidence for and against a cytosolic pH-stat

Phosphoenolpyruvate carboxylase (PEPC) and malic enzyme activities in soluble protein extracts of Avena coleoptiles were investigated to determine whether their kinetics were consistent with a role in cytosol pH regulation. Malic enzyme activity was specific for NADP+ and Mn2+. Maximal labelled product formation from [14C]-substrates required the presence of all coenzymes, cofactors and substrates. Plots of rate versus malate concentration, and linear transformations there- 2 of, indicated typical Michaelis-Menten kinetics at non-saturating malate levels and substrate inhibition at higher malate levels. pH increases between 6.5 and 7.25 increased near-optimal activity, decreased the degree of substrate inhibition and the Kmapp(Mn2+) but did not affect the Vmax or Kmapp(malate). Transformed data of PEPC activity demonstrated non-linear plots indicative of non-Michaelian kinetics. pH increases between 7.0 and 7.6 increased the Vmax and decreased the Km app (Mg2+) but did not affect the Kmapp(PEP). Various carboxylic acids and phosphorylated sugars inhibited PEPC and malic enzyme activities, and these effects decreased with pH increases. Metabolite inhibited malic enzyme activity was non-competitive and resulted mainly from Mn2+ chelation. In contrast, metabolite inhibited PEPC activity was unique for each compound tested, being variously dependent on the PEP concentration and the pH employed. These results indicate that fluctuations in pH and metabolite levels affect PEPC and malic enzyme activities similarly and that 3 the in vitro properties of PEPC are consistent with its proposed role in a pH-stat, whereas the in vitro properties of the malic enzyme cannot be interpreted in terms of a role in pH regulation.

Synthesis, characterization and kinetic studies on complex formed between amantadine hydrochloride and sodium molybdate at physiological pH

The title reaction was undertaken to establish the interaction between amantadine and molybdate at physiological pH. Identical FTIR spectra, TG-DTA curves and CHN data of the complexes formed from three solutions at pH 1.5, 7.4 and 8.0 indicate that the same complex was formed at all the three pHs. The FTIR spectrum shows shift in peaks corresponding to primary amino group of the drug due to coordination to molybdate. An octahedral geometry is assigned to the complex. The kinetics of the complexation has been studied at low concentrations of the reactants using UV-visible spectrophotometry. At pH 7.4, the initial rate varies linearly with [molybdate]. A plot of initial rate versus [drug] is linear passing through origin. These results indicate that the drug and molybdate react at pH 7.4 even at low concentrations. At pH 1.5, the rate increases linearly with increase in [drug] but decreases with [molybdate]. The effect of pH and ionic strength on the rate of the reaction has also been studied. A suitable mechanism has been proposed for the reaction. Reaction between the drug and molybdate even at low concentrations and the fact that the amino group of amantadine required to be free for its function as antiviral, is bound to molybdate in the complex suggests that simultaneous administration of the drug and molybdate supplements should be avoided.

Drastic variation in the surface boundary layer parameters over Cochin during the annular solar eclipse: Analysis using sonic anemometer data

Atmospheric surface boundary layer parameters vary anomalously in response to the occurrence of annular solar eclipse on 15th January 2010 over Cochin. It was the longest annular solar eclipse occurred over South India with high intensity. As it occurred during the noon hours, it is considered to be much more significant because of its effects in all the regions of atmosphere including ionosphere. Since the insolation is the main driving factor responsible for the anomalous changes occurred in the surface layer due to annular solar eclipse, occurred on 15th January 2010, that played very important role in understanding dynamics of the atmosphere during the eclipse period because of its coincidence with the noon time. The Sonic anemometer is able to give data of zonal, meridional and vertical wind as well as the air temperature at a temporal resolution of 1 s. Different surface boundary layer parameters and turbulent fluxes were computed by the application of eddy correlation technique using the high resolution station data. The surface boundary layer parameters that are computed using the sonic anemometer data during the period are momentum flux, sensible heat flux, turbulent kinetic energy, frictional velocity (u*), variance of temperature, variances of u, v and w wind. In order to compare the results, a control run has been done using the data of previous day as well as next day. It is noted that over the specified time period of annular solar eclipse, all the above stated surface boundary layer parameters vary anomalously when compared with the control run. From the observations we could note that momentum flux was 0.1 Nm 2 instead of the mean value 0.2 Nm-2 when there was eclipse. Sensible heat flux anomalously decreases to 50 Nm 2 instead of the mean value 200 Nm 2 at the time of solar eclipse. The turbulent kinetic energy decreases to 0.2 m2s 2 from the mean value 1 m2s 2. The frictional velocity value decreases to 0.05 ms 1 instead of the mean value 0.2 ms 1. The present study aimed at understanding the dynamics of surface layer in response to the annular solar eclipse over a tropical coastal station, occurred during the noon hours. Key words: annular solar eclipse, surface boundary layer, sonic anemometer

Can climate trends be calculated from reanalysis data?

Several global quantities are computed from the ERA40 reanalysis for the period 1958-2001 and explored for trends. These are discussed in the context of changes to the global observing system. Temperature, integrated water vapor (IWV), and kinetic energy are considered. The ERA40 global mean temperature in the lower troposphere has a trend of +0.11 K per decade over the period of 1979-2001, which is slightly higher than the MSU measurements, but within the estimated error limit. For the period 1958 2001 the warming trend is 0.14 K per decade but this is likely to be an artifact of changes in the observing system. When this is corrected for, the warming trend is reduced to 0.10 K per decade. The global trend in IWV for the period 1979-2001 is +0.36 mm per decade. This is about twice as high as the trend determined from the Clausius-Clapeyron relation assuming conservation of relative humidity. It is also larger than results from free climate model integrations driven by the same observed sea surface temperature as used in ERA40. It is suggested that the large trend in IWV does not represent a genuine climate trend but an artifact caused by changes in the global observing system such as the use of SSM/I and more satellite soundings in later years. Recent results are in good agreement with GPS measurements. The IWV trend for the period 1958-2001 is still higher but reduced to +0.16 mm per decade when corrected for changes in the observing systems. Total kinetic energy shows an increasing global trend. Results from data assimilation experiments strongly suggest that this trend is also incorrect and mainly caused by the huge changes in the global observing system in 1979. When this is corrected for, no significant change in global kinetic energy from 1958 onward can be found.