Bayesian L-optimal exact design of experiments for biological kinetic models

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

A joint theoretical and kinetic investigation on the fragmentation of (N-halo)-2-amino cycloalkanecarboxylates

A combined theoretical and experimental study to elucidate the molecular mechanism for the Grob fragmentation of different (N-halo)-2-amino cyclocarboxylates with the nitrogen atom in exocyclic position: (N-Cl)-2-amino cyclopropanecarboxylate (1), (N-Cl)-2-amino cyclobutanecarboxylate (2), (N-Cl)-2-amino cyclopentanecarboxylate (3) and (N-Cl)-2-amino cyclohexanecarboxylate (4), and the corresponding acyclic compounds, (N-Cl)-2-amino isobutyric acid (A), (N-Cl)-2-amino butyric acid (B), has been carried out. The kinetics of decomposition for these compounds and related bromine derivatives were experimentally determined by conventional and stopped-flow UV spectrophotometry. The reaction products have been analyzed by GC and spectrophotometry. Theoretical analysis is based in the localization of stationary points (reactants and transition structures) on the potential energy surface. Calculations were carried out at B3LYP/6-31+G* and MP2/6-31+G* computing methods in the gas phase, while solvent effects have been included by means the self-consistent reaction field theory, PCM continuum model, at MP2/6-31+G* and MP4/6-31+G*//MP2/6-31+G* calculation levels. Based on both experimental and theoretical results, the different Grob fragmentation processes show a global synchronicity index close to 0.9, corresponding to a nearly concerted process. At the TSs, the N-Cl bond breaking is more advanced than the C-C cleavage process. An antiperiplanar configuration of these bonds is reached at the TSs, and this geometrical arrangement is the key factor governing the decomposition. In the case of 1 and 2 the ring strain prevents this spatial disposition, leading to a larger value of the activation barrier. Natural population analysis shows that the polarization of the N-Cl and C-C bonds along the bond-breaking process can be considered the driving force for the decomposition and that a negative charge flows from the carboxylate group to the chlorine atom to assist the reaction pathway. A comparison of theoretical and experimental results shows the relevance of calculation level and the inclusion of solvent effects for determining accurate unimolecular rate coefficients for the decomposition process. © 2002 Published by Elsevier Science B.V.

Compression waves and kinetic energy losses in collisions between balls and rods of different lengths

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

On the Limit of Frequency of Electrochemical Oscillators and Its Relationship to Kinetic Parameters

The class of electrochemical oscillators characterized by a partially hidden negative differential resistance in an N-shaped current potential curve encompasses a myriad of experimental examples. We present a comprehensive methodological analysis of the oscillation frequency of this class of systems and discuss its dependence on electrical and kinetic parameters. The analysis is developed from a skeleton ordinary differential equation model, and an equation for the oscillation frequency is obtained. Simulations are carried out for a model system, namely, the nickel electrodissolution, and the numerical results are confirmed by experimental data on this system. In addition, the treatment is further applied to the electro-oxidation of ethylene glycol where unusually large oscillation frequencies have been reported. Despite the distinct chemistry underlying the oscillatory dynamics of these systems, a very good agreement between experiments and theoretical predictions is observed. The application of the developed theory is suggested as an important step for primary kinetic characterization.

Redox front geochemistry and weathering: theory with application to the Osamu Utsumi uranium mine, Poços de Caldas, Brazil

Numerical calculations describing weathering of the Poços de Caldas alkaline complex (Minas Gerais, Brazil) by infiltrating groundwater are carried out for time spans up to two million years in the absence of pyrite, and up to 500,000 years with pyrite present. Deposition of uranium resulting from infiltration of oxygenated, uranium bearing groundwater through the hydrothermally altered phonolitic host rock at the Osamu Utsumi uranium mine is also included in the latter calculation. The calculations are based on the quasi-stationary state approximation to mass conservation equations for pure advective transport. This approximation enables the prediction of solute concentrations, mineral abundances and porosity as functions of time and distance over geologic time spans. Mineral reactions are described by kinetic rate laws for both precipitation and dissolution. Homogeneous equilibrium is assumed to be maintained within the aqueous phase. No other constraints are imposed on the calculations other than the initial composition of the unaltered host rock and the composition of the inlet fluid, taken as rainwater modified by percolation through a soil zone. The results are in qualitative agreement with field observations at the Osamu Utsumi uranium mine. They predict a lateritic cover followed by a highly porous saprolitic zone, a zone of oxidized rock with pyrite replaced by iron-hydroxide, a sharp redox front at which uranium is deposited, and the reduced unweathered host rock. Uranium is deposited in a narrow zone located on the reduced side of the redox front in association with pyrite, in agreement with field observations. The calculations predict the formation of a broad dissolution front of primary kaolinite that penetrates deep into the host rock accompanied by the precipitation of secondary illite. Secondary kaolinite occurs in a saprolitic zone near the surface and in the vicinity of the redox front. Gibbsite forms a bi-modal distribution consisting of a maximum near the surface followed by a thin tongue extending downward into the weathered profile in agreement with field observations. The results are found to be insensitive to the kinetic rate constants used to describe mineral reactions.

Object Kinetic Monte Carlo calculations of irradiated Fe-Cr dilute alloys: The effect of the interaction radius between substitutional Cr and self-interstitial Fe

ObjectKineticMonteCarlo models allow for the study of the evolution of the damage created by irradiation to time scales that are comparable to those achieved experimentally. Therefore, the essential ObjectKineticMonteCarlo parameters can be validated through comparison with experiments. However, this validation is not trivial since a large number of parameters is necessary, including migration energies of point defects and their clusters, binding energies of point defects in clusters, as well as the interactionradii. This is particularly cumbersome when describing an alloy, such as the Fe–Cr system, which is of interest for fusion energy applications. In this work we describe an ObjectKineticMonteCarlo model for Fe–Cr alloys in the dilute limit. The parameters used in the model come either from density functional theory calculations or from empirical interatomic potentials. This model is used to reproduce isochronal resistivity recovery experiments of electron irradiateddiluteFe–Cr alloys performed by Abe and Kuramoto. The comparison between the calculated results and the experiments reveal that an important parameter is the capture radius between substitutionalCr and self-interstitialFe atoms. A parametric study is presented on the effect of the capture radius on the simulated recovery curves.

MMonCa: A flexible and powerful new Kinetic Monte Carlo simulator

Kinetic Monte Carlo (KMC) is a widely used technique to simulate the evolution of radiation damage inside solids. Despite de fact that this technique was developed several decades ago, there is not an established and easy to access simulating tool for researchers interested in this field, unlike in the case of molecular dynamics or density functional theory calculations. In fact, scientists must develop their own tools or use unmaintained ones in order to perform these types of simulations. To fulfil this need, we have developed MMonCa, the Modular Monte Carlo simulator. MMonCa has been developed using professional C++ programming techniques and has been built on top of an interpreted language to allow having a powerful yet flexible, robust but customizable and easy to access modern simulator. Both non lattice and Lattice KMC modules have been developed. We will present in this conference, for the first time, the MMonCa simulator. Along with other (more detailed) contributions in this meeting, the versatility of MMonCa to study a number of problems in different materials (particularly, Fe and W) subject to a wide range of conditions will be shown. Regarding KMC simulations, we have studied neutron-generated cascade evolution in Fe (as a model material). Starting with a Frenkel pair distribution we have followed the defect evolution up to 450 K. Comparison with previous simulations and experiments shows excellent agreement. Furthermore, we have studied a more complex system (He-irradiated W:C) using a previous parametrization [1]. He-irradiation at 4 K followed by isochronal annealing steps up to 500 K has been simulated with MMonCa. The He energy was 400 eV or 3 keV. In the first case, no damage is associated to the He implantation, whereas in the second one, a significant Frenkel pair concentration (evolving into complex clusters) is associated to the He ions. We have been able to explain He desorption both in the absence and in the presence of Frenkel pairs and we have also applied MMonCa to high He doses and fluxes at elevated temperatures. He migration and trapping dominate the kinetics of He desorption. These processes will be discussed and compared to experimental results. [1] C.S. Becquart et al. J. Nucl. Mater. 403 (2010) 75

The theory of electrostatic probes in strong magnetic fields : Thesis submitted to the Faculty of The Graduate School of the University of Colorado for the Degree Doctor of Philosophy Department of Aerospace Engineering Science

A theory is developed of an electrostatic probe in a fully-ionized plasma in the presence of a strong magnetic field. The ratio of electron Larmor radius to probe transverse dimension is assumed to be small. Poisson's equation, together with kinetic equations for ions and electrons are considered. An asymptotic perturbation method of multiple scales is used by considering the characteristic lengths appearing in the problem. The leading behavior of the solution is found. The results obtained appear to apply to weaker fields also, agreeing with the solutions known in the limit of no magnetic field. The range of potentials for wich results are presented is limited. The basic effects produced by the field are a depletion of the plasma near the probe and a non-monotonic potential surrounding the probe. The ion saturation current is not changed but changes appear in both the floating potential Vf and the slope of the current-voltage diagram at Vf. The transition region extends beyond the space potential Vs,at wich point the current is largely reduced. The diagram does not have an exponential form in this region as commonly assumed. There exists saturation in electron collection. The extent to which the plasma is disturbed is determined. A cylindrical probe has no solution because of a logarithmic singularity at infinity. Extensions of the theory are considered.

Theory of a Probe in a Strong Magnetic Field. AT(30-1)-1238-MATT 599

A kinetic approach is used to develop a theory of electrostatic probes in a fully ionized plasma in the presence of a magnetic field. A consistent asymptotic expansion is obtained assuming that the electron Larmor radius is small compared to the radius of the probe. The order of magnitude of neglected terms is given. It is found that the electric potential within the tube of force defined by the cross section of the probe decays non-mono tonic ally from the probe; this bump disappears at a certain probe voltage and the theory is valid up to this voltage. The transition region, which extends beyond plasma potential, is not exponential. The possible saturation of the electron current is discussed. Restricted numerical results are given; they seem to be useful for weaker magnetic fields down to the zero-field limit. Extensions of the theory a r e considered.

MMonCa: An Object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion

In this work, we introduce the Object Kinetic Monte Carlo (OKMC) simulator MMonCa and simulate the defect evolution in three different materials. We start by explaining the theory of OKMC and showing some details of how such theory is implemented by creating generic structures and algorithms in the objects that we want to simulate. Then we successfully reproduce simulated results for defect evolution in iron, silicon and tungsten using our simulator and compare with available experimental data and similar simulations. The comparisons validate MMonCa showing that it is powerful and flexible enough to be customized and used to study the damage evolution of defects in a wide range of solid materials.

Force and kinetic barriers to unzipping of the DNA double helix

A theory of the unzipping of double-stranded DNA is presented and is compared to recent micromanipulation experiments. It is shown that the interactions that stabilize the double helix and the elastic rigidity of single strands simply determine the sequence-dependent ≈12-pN force threshold for DNA strand separation. Using a semimicroscopic model of the binding between nucleotide strands, we show that the greater rigidity of the strands when formed into double-stranded DNA, relative to that of isolated strands, gives rise to a potential barrier to unzipping. The effects of this barrier are derived analytically. The force to keep the extremities of the molecule at a fixed distance, the kinetic rates for strand unpairing at fixed applied force, and the rupture force as a function of loading rate are calculated. The dependence of the kinetics and of the rupture force on molecule length is also analyzed.

Well posedness of an integrodifferential kinetic model of Fokker-Planck type for angiogenesis.

Tumor induced angiogenesis processes including the effect of stochastic motion and branching of blood vessels can be described coupling a (nonlocal in time) integrodifferential kinetic equation of Fokker–Planck type with a diffusion equation for the tumor induced ingiogenic factor. The chemotactic force field depends on the flux of blood vessels through the angiogenic factor. We develop an existence and uniqueness theory for this system under natural assumptions on the initial data. The proof combines the construction of fundamental solutions for associated linearized problems with comparison principles, sharp estimates of the velocity integrals and compactness results for this type of kinetic and parabolic operators

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Mechanisms influencing levitation and the scaling laws in nanopores: Oscillator model theory

We provide here a detailed theoretical explanation of the floating molecule or levitation effect, for molecules diffusing through nanopores, using the oscillator model theory (Phys. Rev. Lett. 2003, 91, 126102) recently developed in this laboratory. It is shown that on reduction of pore size the effect occurs due to decrease in frequency of wall collision of diffusing particles at a critical pore size. This effect is, however, absent at high temperatures where the ratio of kinetic energy to the solid-fluid interaction strength is sufficiently large. It is shown that the transport diffusivities scale with this ratio. Scaling of transport diffusivities with respect to mass is also observed, even in the presence of interactions.

Kinetic studies of the metathesis of 1-hexene using Re2O7/γ-Al2O3 and the synthesis and metathesis of oxazolines

The kinetics of the metathesis of 1-hexene using Re2O7/-Al_2O_3 as the catalyst were investigated under a variety of conditions. The experiments were carried out under high vacuum conditions. The product solutions were characterised by gas liquid chromatography and mass spectroscopy. The initial kinetics of the metathesis of 1-hexene showed that the reaction was first order in the weight of the catalyst and second order in the concentration of 1-hexene. A kinetic scheme which correlated the experimental data with the metallocarbene chain mechanism postulated by Herisson and Chauvin and the kinetics of the reaction was explained using a model based on the Langmuir-Hinshelwood theory. The low conversion of 1-hexene to its products is due to termination reactions which most likely occur by the decomposition of the metallocyclobutane intermediate to produce a cyclopropane derivative and an inactive centre. The optimum temperature for the metathesis of 1-hexene over Re_2O_7/-Al2O3 is 45oC and above this temperature, the rate of metathesis decreases rapidly. Co-catalysts alter the active sites for metathesis so that the catalyst is more selective to the metathesis of 1-hexene. However, the regeneration of metathesis activity is much worse for promoted catalysts than for the unpromoted. The synthesis and metathesis of 4,4-dimethyl-2-allowbreak (9-decenyl)-1,3-oxazoline and 4,4-dimethyl-2-allowbreak (3-pentenyl)-1,3-oxazoline was attempted and the products were analysed by thin layer chromatography, infra-red, 13C and 1H nmr and mass spectroscopy. Obtaining the oxazolines in a good yield with high purity was difficult and consequently metathesis of the impure products did not occur.