The Kinetics of Dissolution Revisited

NASA Astrophysics Data System (ADS)

Antonel, Paula S.; Hoijemberg, Pablo A.; Maiante, Leandro M.; Lagorio, M. Gabriela

2003-09-01

An experiment analyzing the kinetics of dissolution of a solid with cylindrical geometry in water is presented. The dissolution process is followed by measuring the solid mass and its size parameters (thickness and diameter) as a function of time. It is verified that the dissolution rate follows the Nernst model. Data treatment is compared with the dissolution of a spherical solid previously described. Kinetics, diffusion concepts, and polynomial fitting of experimental data are combined in this simple experiment.

A-Priori Tuning of Modified Magnussen Combustion Model

NASA Technical Reports Server (NTRS)

Norris, A. T.

2016-01-01

In the application of CFD to turbulent reacting flows, one of the main limitations to predictive accuracy is the chemistry model. Using a full or skeletal kinetics model may provide good predictive ability, however, at considerable computational cost. Adding the ability to account for the interaction between turbulence and chemistry improves the overall fidelity of a simulation but adds to this cost. An alternative is the use of simple models, such as the Magnussen model, which has negligible computational overhead, but lacks general predictive ability except for cases that can be tuned to the flow being solved. In this paper, a technique will be described that allows the tuning of the Magnussen model for an arbitrary fuel and flow geometry without the need to have experimental data for that particular case. The tuning is based on comparing the results of the Magnussen model and full finite-rate chemistry when applied to perfectly and partially stirred reactor simulations. In addition, a modification to the Magnussen model is proposed that allows the upper kinetic limit for the reaction rate to be set, giving better physical agreement with full kinetic mechanisms. This procedure allows a simple reacting model to be used in a predictive manner, and affords significant savings in computational costs for simulations.

A solvable model of Vlasov-kinetic plasma turbulence in Fourier-Hermite phase space

NASA Astrophysics Data System (ADS)

Adkins, T.; Schekochihin, A. A.

2018-02-01

A class of simple kinetic systems is considered, described by the one-dimensional Vlasov-Landau equation with Poisson or Boltzmann electrostatic response and an energy source. Assuming a stochastic electric field, a solvable model is constructed for the phase-space turbulence of the particle distribution. The model is a kinetic analogue of the Kraichnan-Batchelor model of chaotic advection. The solution of the model is found in Fourier-Hermite space and shows that the free-energy flux from low to high Hermite moments is suppressed, with phase mixing cancelled on average by anti-phase-mixing (stochastic plasma echo). This implies that Landau damping is an ineffective route to dissipation (i.e. to thermalisation of electric energy via velocity space). The full Fourier-Hermite spectrum is derived. Its asymptotics are -3/2$ at low wavenumbers and high Hermite moments ( ) and -1/2k-2$ at low Hermite moments and high wavenumbers ( ). These conclusions hold at wavenumbers below a certain cutoff (analogue of Kolmogorov scale), which increases with the amplitude of the stochastic electric field and scales as inverse square of the collision rate. The energy distribution and flows in phase space are a simple and, therefore, useful example of competition between phase mixing and nonlinear dynamics in kinetic turbulence, reminiscent of more realistic but more complicated multi-dimensional systems that have not so far been amenable to complete analytical solution.

Exponential growth kinetics for Polyporus versicolor and Pleurotus ostreatus in submerged culture

DOE Office of Scientific and Technical Information (OSTI.GOV)

Carroad, P.A.; Wilke, C.R.

1977-04-01

Simple mathematical models for a batch culture of pellet-forming fungi in submerged culture were tested on growth data for Polyporus versicolor (ATCC 12679) and Pleurotus ostreatus (ATCC 9415). A kinetic model based on a growth rate proportional to the two-thirds power of the cell mass was shown to be satisfactory. A model based on a growth rate directly proportional to the cell mass fitted the data equally well, however, and may be preferable because of mathematical simplicity.

Modeling of Non-isothermal Austenite Formation in Spring Steel

NASA Astrophysics Data System (ADS)

Huang, He; Wang, Baoyu; Tang, Xuefeng; Li, Junling

2017-12-01

The austenitization kinetics description of spring steel 60Si2CrA plays an important role in providing guidelines for industrial production. The dilatometric curves of 60Si2CrA steel were measured using a dilatometer DIL805A at heating rates of 0.3 K to 50 K/s (0.3 °C/s to 50 °C/s). Based on the dilatometric curves, a unified kinetics model using the internal state variable (ISV) method was derived to describe the non-isothermal austenitization kinetics of 60Si2CrA, and the abovementioned model models the incubation and transition periods. The material constants in the model were determined using a genetic algorithm-based optimization technique. Additionally, good agreement between predicted and experimental volume fractions of transformed austenite was obtained, indicating that the model is effective for describing the austenitization kinetics of 60Si2CrA steel. Compared with other modeling methods of austenitization kinetics, this model, which uses the ISV method, has some advantages, such as a simple formula and explicit physics meaning, and can be probably used in engineering practice.

Principles of protein folding--a perspective from simple exact models.

PubMed Central

Dill, K. A.; Bromberg, S.; Yue, K.; Fiebig, K. M.; Yee, D. P.; Thomas, P. D.; Chan, H. S.

1995-01-01

General principles of protein structure, stability, and folding kinetics have recently been explored in computer simulations of simple exact lattice models. These models represent protein chains at a rudimentary level, but they involve few parameters, approximations, or implicit biases, and they allow complete explorations of conformational and sequence spaces. Such simulations have resulted in testable predictions that are sometimes unanticipated: The folding code is mainly binary and delocalized throughout the amino acid sequence. The secondary and tertiary structures of a protein are specified mainly by the sequence of polar and nonpolar monomers. More specific interactions may refine the structure, rather than dominate the folding code. Simple exact models can account for the properties that characterize protein folding: two-state cooperativity, secondary and tertiary structures, and multistage folding kinetics--fast hydrophobic collapse followed by slower annealing. These studies suggest the possibility of creating "foldable" chain molecules other than proteins. The encoding of a unique compact chain conformation may not require amino acids; it may require only the ability to synthesize specific monomer sequences in which at least one monomer type is solvent-averse. PMID:7613459

Uncovering Oscillations, Complexity, and Chaos in Chemical Kinetics Using Mathematica

NASA Astrophysics Data System (ADS)

Ferreira, M. M. C.; Ferreira, W. C., Jr.; Lino, A. C. S.; Porto, M. E. G.

1999-06-01

Unlike reactions with no peculiar temporal behavior, in oscillatory reactions concentrations can rise and fall spontaneously in a cyclic or disorganized fashion. In this article, the software Mathematica is used for a theoretical study of kinetic mechanisms of oscillating and chaotic reactions. A first simple example is introduced through a three-step reaction, called the Lotka model, which exhibits a temporal behavior characterized by damped oscillations. The phase plane method of dynamic systems theory is introduced for a geometric interpretation of the reaction kinetics without solving the differential rate equations. The equations are later numerically solved using the built-in routine NDSolve and the results are plotted. The next example, still with a very simple mechanism, is the Lotka-Volterra model reaction, which oscillates indefinitely. The kinetic process and rate equations are also represented by a three-step reaction mechanism. The most important difference between this and the former reaction is that the undamped oscillation has two autocatalytic steps instead of one. The periods of oscillations are obtained by using the discrete Fourier transform (DFT)-a well-known tool in spectroscopy, although not so common in this context. In the last section, it is shown how a simple model of biochemical interactions can be useful to understand the complex behavior of important biological systems. The model consists of two allosteric enzymes coupled in series and activated by its own products. This reaction scheme is important for explaining many metabolic mechanisms, such as the glycolytic oscillations in muscles, yeast glycolysis, and the periodic synthesis of cyclic AMP. A few of many possible dynamic behaviors are exemplified through a prototype glycolytic enzymatic reaction proposed by Decroly and Goldbeter. By simply modifying the initial concentrations, limit cycles, chaos, and birhythmicity are computationally obtained and visualized.

Integral Design Methodology of Photocatalytic Reactors for Air Pollution Remediation.

PubMed

Passalía, Claudio; Alfano, Orlando M; Brandi, Rodolfo J

2017-06-07

An integral reactor design methodology was developed to address the optimal design of photocatalytic wall reactors to be used in air pollution control. For a target pollutant to be eliminated from an air stream, the proposed methodology is initiated with a mechanistic derived reaction rate. The determination of intrinsic kinetic parameters is associated with the use of a simple geometry laboratory scale reactor, operation under kinetic control and a uniform incident radiation flux, which allows computing the local superficial rate of photon absorption. Thus, a simple model can describe the mass balance and a solution may be obtained. The kinetic parameters may be estimated by the combination of the mathematical model and the experimental results. The validated intrinsic kinetics obtained may be directly used in the scaling-up of any reactor configuration and size. The bench scale reactor may require the use of complex computational software to obtain the fields of velocity, radiation absorption and species concentration. The complete methodology was successfully applied to the elimination of airborne formaldehyde. The kinetic parameters were determined in a flat plate reactor, whilst a bench scale corrugated wall reactor was used to illustrate the scaling-up methodology. In addition, an optimal folding angle of the corrugated reactor was found using computational fluid dynamics tools.

Predicting perturbation patterns from the topology of biological networks.

PubMed

Santolini, Marc; Barabási, Albert-László

2018-06-20

High-throughput technologies, offering an unprecedented wealth of quantitative data underlying the makeup of living systems, are changing biology. Notably, the systematic mapping of the relationships between biochemical entities has fueled the rapid development of network biology, offering a suitable framework to describe disease phenotypes and predict potential drug targets. However, our ability to develop accurate dynamical models remains limited, due in part to the limited knowledge of the kinetic parameters underlying these interactions. Here, we explore the degree to which we can make reasonably accurate predictions in the absence of the kinetic parameters. We find that simple dynamically agnostic models are sufficient to recover the strength and sign of the biochemical perturbation patterns observed in 87 biological models for which the underlying kinetics are known. Surprisingly, a simple distance-based model achieves 65% accuracy. We show that this predictive power is robust to topological and kinetic parameter perturbations, and we identify key network properties that can increase up to 80% the recovery rate of the true perturbation patterns. We validate our approach using experimental data on the chemotactic pathway in bacteria, finding that a network model of perturbation spreading predicts with ∼80% accuracy the directionality of gene expression and phenotype changes in knock-out and overproduction experiments. These findings show that the steady advances in mapping out the topology of biochemical interaction networks opens avenues for accurate perturbation spread modeling, with direct implications for medicine and drug development.

A comparison of simple global kinetic models for coal devolatilization with the CPD model

DOE PAGES

Richards, Andrew P.; Fletcher, Thomas H.

2016-08-01

Simulations of coal combustors and gasifiers generally cannot incorporate the complexities of advanced pyrolysis models, and hence there is interest in evaluating simpler models over ranges of temperature and heating rate that are applicable to the furnace of interest. In this paper, six different simple model forms are compared to predictions made by the Chemical Percolation Devolatilization (CPD) model. The model forms included three modified one-step models, a simple two-step model, and two new modified two-step models. These simple model forms were compared over a wide range of heating rates (5 × 10 3 to 10 6 K/s) at finalmore » temperatures up to 1600 K. Comparisons were made of total volatiles yield as a function of temperature, as well as the ultimate volatiles yield. Advantages and disadvantages for each simple model form are discussed. In conclusion, a modified two-step model with distributed activation energies seems to give the best agreement with CPD model predictions (with the fewest tunable parameters).« less

DOE Office of Scientific and Technical Information (OSTI.GOV)

Seapan, M.; Crynes, B.L.; Dale, S.

The objectives of this study were to analyze alternate crudes kinetic hydrotreatment data in the literature, develop a mathematical model for interpretation of these data, develop an experimental procedure and apparatus to collect accurate kinetic data, and finally, to combine the model and experimental data to develop a general model which, with a few experimental parameters, could be used in design of future hydrotreatment processes. These objectives were to cover a four year program (1980 to 1984) and were subjective to sufficient funding. Only partial funding has been available thus far to cover activities for two years. A hydrotreatment datamore » base is developed which contains over 2000 citations, stored in a microcomputer. About 50% of these are reviewed, classified and can be identified by feedstock, catalyst, reactor type and other process characteristics. Tests of published hydrodesulfurization data indicate the problems with simple n-th order, global kinetic models, and point to the value of developing intrinsic reaction kinetic models to describe the reaction processes. A Langmuir-Hinshelwood kinetic model coupled with a plug flow reactor design equation has been developed and used for published data evaluation. An experimental system and procedure have been designed and constructed, which can be used for kinetic studies. 30 references, 4 tables.« less

Receptor binding kinetics equations: Derivation using the Laplace transform method.

PubMed

Hoare, Sam R J

Measuring unlabeled ligand receptor binding kinetics is valuable in optimizing and understanding drug action. Unfortunately, deriving equations for estimating kinetic parameters is challenging because it involves calculus; integration can be a frustrating barrier to the pharmacologist seeking to measure simple rate parameters. Here, a well-known tool for simplifying the derivation, the Laplace transform, is applied to models of receptor-ligand interaction. The method transforms differential equations to a form in which simple algebra can be applied to solve for the variable of interest, for example the concentration of ligand-bound receptor. The goal is to provide instruction using familiar examples, to enable investigators familiar with handling equilibrium binding equations to derive kinetic equations for receptor-ligand interaction. First, the Laplace transform is used to derive the equations for association and dissociation of labeled ligand binding. Next, its use for unlabeled ligand kinetic equations is exemplified by a full derivation of the kinetics of competitive binding equation. Finally, new unlabeled ligand equations are derived using the Laplace transform. These equations incorporate a pre-incubation step with unlabeled or labeled ligand. Four equations for measuring unlabeled ligand kinetics were compared and the two new equations verified by comparison with numerical solution. Importantly, the equations have not been verified with experimental data because no such experiments are evident in the literature. Equations were formatted for use in the curve-fitting program GraphPad Prism 6.0 and fitted to simulated data. This description of the Laplace transform method will enable pharmacologists to derive kinetic equations for their model or experimental paradigm under study. Application of the transform will expand the set of equations available for the pharmacologist to measure unlabeled ligand binding kinetics, and for other time-dependent pharmacological activities. Copyright © 2017 Elsevier Inc. All rights reserved.

Relativistic Corrections to the Bohr Model of the Atom

ERIC Educational Resources Information Center

Kraft, David W.

1974-01-01

Presents a simple means for extending the Bohr model to include relativistic corrections using a derivation similar to that for the non-relativistic case, except that the relativistic expressions for mass and kinetic energy are employed. (Author/GS)

Kinetic electron model for plasma thruster plumes

NASA Astrophysics Data System (ADS)

Merino, Mario; Mauriño, Javier; Ahedo, Eduardo

2018-03-01

A paraxial model of an unmagnetized, collisionless plasma plume expanding into vacuum is presented. Electrons are treated kinetically, relying on the adiabatic invariance of their radial action integral for the integration of Vlasov's equation, whereas ions are treated as a cold species. The quasi-2D plasma density, self-consistent electric potential, and electron pressure, temperature, and heat fluxes are analyzed. In particular, the model yields the collisionless cooling of electrons, which differs from the Boltzmann relation and the simple polytropic laws usually employed in fluid and hybrid PIC/fluid plume codes.

A model for tracking concentration of chemical compounds within a tank of an automatic film processor.

PubMed

Sobol, Wlad T

2002-01-01

A simple kinetic model that describes the time evolution of the chemical concentration of an arbitrary compound within the tank of an automatic film processor is presented. It provides insights into the kinetics of chemistry concentration inside the processor's tank; the results facilitate the tasks of processor tuning and quality control (QC). The model has successfully been used in several troubleshooting sessions of low-volume mammography processors for which maintaining consistent QC tracking was difficult due to fluctuations of bromide levels in the developer tank.

Glassy Behavior due to Kinetic Constraints: from Topological Foam to Backgammon

NASA Astrophysics Data System (ADS)

Sherrington, David

A study is reported of a series of simple model systems with only non-interacting Hamiltonians, and hence simple equilibrium thermodynamics, but with constrained kinetics of a type initially suggested by topological considerations of foams and two-dimensional covalent glasses. It is demonstrated that oscopic dynamical features characteristic of real glasses, such as two-time decays in energy and auto-correlation functions, arise and may be understood in terms of annihilation-diffusion concepts and theory. This recognition leads to a sequence of further models which (i) encapsulate the essense but are more readily simulated and open to easier analytic study, and (ii) allow generalization and extension to higher dimension. Fluctuation-dissipation relations are also considered and show novel aspects. The comparison is with strong glasses.

KINETIC MODELING OF NOX FORMATION AND DESTRUCTION AND COMBUSTIBLES BURNOUT

EPA Science Inventory

The report describes a model of the gas-phase chemistry involved in the combustion of simple hydrocarbon fuels and the interconversion of fixed nitrogen species. One focus of the work was on modeling the chemistry involved in reburning and other advanced NOx control strategies. A...

Shear-induced reaction-limited aggregation kinetics of Brownian particles at arbitrary concentrations

NASA Astrophysics Data System (ADS)

Zaccone, Alessio; Gentili, Daniele; Wu, Hua; Morbidelli, Massimo

2010-04-01

The aggregation of interacting Brownian particles in sheared concentrated suspensions is an important issue in colloid and soft matter science per se. Also, it serves as a model to understand biochemical reactions occurring in vivo where both crowding and shear play an important role. We present an effective medium approach within the Smoluchowski equation with shear which allows one to calculate the encounter kinetics through a potential barrier under shear at arbitrary colloid concentrations. Experiments on a model colloidal system in simple shear flow support the validity of the model in the concentration range considered. By generalizing Kramers' rate theory to the presence of shear and collective hydrodynamics, our model explains the significant increase in the shear-induced reaction-limited aggregation kinetics upon increasing the colloid concentration.

Topological and kinetic determinants of the modal matrices of dynamic models of metabolism

PubMed Central

2017-01-01

Large-scale kinetic models of metabolism are becoming increasingly comprehensive and accurate. A key challenge is to understand the biochemical basis of the dynamic properties of these models. Linear analysis methods are well-established as useful tools for characterizing the dynamic response of metabolic networks. Central to linear analysis methods are two key matrices: the Jacobian matrix (J) and the modal matrix (M-1) arising from its eigendecomposition. The modal matrix M-1 contains dynamically independent motions of the kinetic model near a reference state, and it is sparse in practice for metabolic networks. However, connecting the structure of M-1 to the kinetic properties of the underlying reactions is non-trivial. In this study, we analyze the relationship between J, M-1, and the kinetic properties of the underlying network for kinetic models of metabolism. Specifically, we describe the origin of mode sparsity structure based on features of the network stoichiometric matrix S and the reaction kinetic gradient matrix G. First, we show that due to the scaling of kinetic parameters in real networks, diagonal dominance occurs in a substantial fraction of the rows of J, resulting in simple modal structures with clear biological interpretations. Then, we show that more complicated modes originate from topologically-connected reactions that have similar reaction elasticities in G. These elasticities represent dynamic equilibrium balances within reactions and are key determinants of modal structure. The work presented should prove useful towards obtaining an understanding of the dynamics of kinetic models of metabolism, which are rooted in the network structure and the kinetic properties of reactions. PMID:29267329

Simulations of thermodynamics and kinetics on rough energy landscapes with milestoning.

PubMed

Bello-Rivas, Juan M; Elber, Ron

2016-03-05

We investigated by computational means the kinetics and stationary behavior of stochastic dynamics on an ensemble of rough two-dimensional energy landscapes. There are no obvious separations of temporal scales in these systems, which constitute a simple model for the behavior of glasses and some biomaterials. Even though there are significant computational challenges present in these systems due to the large number of metastable states, the Milestoning method is able to compute their kinetic and thermodynamic properties exactly. We observe two clearly distinguished regimes in the overall kinetics: one in which diffusive behavior dominates and another that follows an Arrhenius law (despite the absence of a dominant barrier). We compare our results with those obtained with an exactly-solvable one-dimensional model, and with the results from the rough one-dimensional energy model introduced by Zwanzig. © 2015 Wiley Periodicals, Inc. © 2015 Wiley Periodicals, Inc.

A comparison of quantitative methods for clinical imaging with hyperpolarized (13)C-pyruvate.

PubMed

Daniels, Charlie J; McLean, Mary A; Schulte, Rolf F; Robb, Fraser J; Gill, Andrew B; McGlashan, Nicholas; Graves, Martin J; Schwaiger, Markus; Lomas, David J; Brindle, Kevin M; Gallagher, Ferdia A

2016-04-01

Dissolution dynamic nuclear polarization (DNP) enables the metabolism of hyperpolarized (13)C-labelled molecules, such as the conversion of [1-(13)C]pyruvate to [1-(13)C]lactate, to be dynamically and non-invasively imaged in tissue. Imaging of this exchange reaction in animal models has been shown to detect early treatment response and correlate with tumour grade. The first human DNP study has recently been completed, and, for widespread clinical translation, simple and reliable methods are necessary to accurately probe the reaction in patients. However, there is currently no consensus on the most appropriate method to quantify this exchange reaction. In this study, an in vitro system was used to compare several kinetic models, as well as simple model-free methods. Experiments were performed using a clinical hyperpolarizer, a human 3 T MR system, and spectroscopic imaging sequences. The quantitative methods were compared in vivo by using subcutaneous breast tumours in rats to examine the effect of pyruvate inflow. The two-way kinetic model was the most accurate method for characterizing the exchange reaction in vitro, and the incorporation of a Heaviside step inflow profile was best able to describe the in vivo data. The lactate time-to-peak and the lactate-to-pyruvate area under the curve ratio were simple model-free approaches that accurately represented the full reaction, with the time-to-peak method performing indistinguishably from the best kinetic model. Finally, extracting data from a single pixel was a robust and reliable surrogate of the whole region of interest. This work has identified appropriate quantitative methods for future work in the analysis of human hyperpolarized (13)C data. © 2016 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

Kinetics of steel slag leaching: Batch tests and modeling

DOE Office of Scientific and Technical Information (OSTI.GOV)

De Windt, Laurent, E-mail: laurent.dewindt@mines-paristech.fr; Chaurand, Perrine; Rose, Jerome

2011-02-15

Reusing steel slag as an aggregate for road construction requires to characterize the leaching kinetics and metal releases. In this study, basic oxygen furnace (BOF) steel slag were subjected to batch leaching tests at liquid to solid ratios (L/S) of 10 and 100 over 30 days; the leachate chemistry being regularly sampled in time. A geochemical model of the steel slag is developed and validated from experimental data, particularly the evolution with leaching of mineralogical composition of the slag and trace element speciation. Kinetics is necessary for modeling the primary phase leaching, whereas a simple thermodynamic equilibrium approach can bemore » used for secondary phase precipitation. The proposed model simulates the kinetically-controlled dissolution (hydrolysis) of primary phases, the precipitation of secondary phases (C-S-H, hydroxide and spinel), the pH and redox conditions, and the progressive release of major elements as well as the metals Cr and V. Modeling indicates that the dilution effect of the L/S ratio is often coupled to solubility-controlled processes, which are sensitive to both the pH and the redox potential. A sensitivity analysis of kinetic uncertainties on the modeling of element releases is performed.« less

Comparison of ACCENT 2000 Shuttle Plume Data with SIMPLE Model Predictions

NASA Astrophysics Data System (ADS)

Swaminathan, P. K.; Taylor, J. C.; Ross, M. N.; Zittel, P. F.; Lloyd, S. A.

2001-12-01

The JHU/APL Stratospheric IMpact of PLume Effluents (SIMPLE)model was employed to analyze the trace species in situ composition data collected during the ACCENT 2000 intercepts of the space shuttle Space Transportation Launch System (STS) rocket plume as a function of time and radial location within the cold plume. The SIMPLE model is initialized using predictions for species depositions calculated using an afterburning model based on standard TDK/SPP nozzle and SPF plume flowfield codes with an expanded chemical kinetic scheme. The time dependent ambient stratospheric chemistry is fully coupled to the plume species evolution whose transport is based on empirically derived diffusion. Model/data comparisons are encouraging through capturing observed local ozone recovery times as well as overall morphology of chlorine chemistry.

Modified kinetic theory applied to the shear flows of granular materials

DOE PAGES

Duan, Yifei; Feng, Zhi -Gang; Michaelides, Efstathios E.; ...

2017-04-11

Here, granular materials are characterized by large collections of discrete particles, where the particle-particle interactions are significantly more important than the particle-fluid interactions. The current kinetic theory captures fairly accurately the granular flow behavior in the dilute case, when only binary interactions are significant, but is not accurate at all in the dense flow regime, where multi-particle interactions and contacts must be modeled. To improve the kinetic theory results for granular flows in the dense flow regime, we propose a Modified Kinetic Theory (MKT) model that utilizes the contact duration or cut-off time to account for the complex particle-particle interactionsmore » in the dense regime. The contact duration model, also called TC model, is originally proposed by Luding and McNamara to solve the inelastic collapse issue existing in the Inelastic Hard Sphere (IHS) model. This model defines a cut-off time t c such that dissipation is not counted if the time between two consecutive contacts is less than t c. As shown in their study, the use of a cut-off time t c can also reduce the dissipation during multi-particle contacts. In this paper we relate the TC model with the Discrete Element Method (DEM) by choosing the cut-off time t c to be the duration of contact calculated from the linear-spring-dashpot soft-sphere model of the DEM. We examine two types of granular flows: simple shear flow and the plane shear flow, and compare the results of the classical Kinetic Theory (KT) model, the present MKT model, and the DEM model. Here, we show that the MKT model entails a significant improvement over the KT model for simple shear flows at inertial regimes. With the MKT model the calculations are close to the DEM results at solid fractions as high as 0.57. Even for the plane shear flows, where shear rate and solid fraction are inhomogeneous, the results of the MKT model agree very well with the DEM results.« less

Modified kinetic theory applied to the shear flows of granular materials

DOE Office of Scientific and Technical Information (OSTI.GOV)

Duan, Yifei; Feng, Zhi -Gang; Michaelides, Efstathios E.

Here, granular materials are characterized by large collections of discrete particles, where the particle-particle interactions are significantly more important than the particle-fluid interactions. The current kinetic theory captures fairly accurately the granular flow behavior in the dilute case, when only binary interactions are significant, but is not accurate at all in the dense flow regime, where multi-particle interactions and contacts must be modeled. To improve the kinetic theory results for granular flows in the dense flow regime, we propose a Modified Kinetic Theory (MKT) model that utilizes the contact duration or cut-off time to account for the complex particle-particle interactionsmore » in the dense regime. The contact duration model, also called TC model, is originally proposed by Luding and McNamara to solve the inelastic collapse issue existing in the Inelastic Hard Sphere (IHS) model. This model defines a cut-off time t c such that dissipation is not counted if the time between two consecutive contacts is less than t c. As shown in their study, the use of a cut-off time t c can also reduce the dissipation during multi-particle contacts. In this paper we relate the TC model with the Discrete Element Method (DEM) by choosing the cut-off time t c to be the duration of contact calculated from the linear-spring-dashpot soft-sphere model of the DEM. We examine two types of granular flows: simple shear flow and the plane shear flow, and compare the results of the classical Kinetic Theory (KT) model, the present MKT model, and the DEM model. Here, we show that the MKT model entails a significant improvement over the KT model for simple shear flows at inertial regimes. With the MKT model the calculations are close to the DEM results at solid fractions as high as 0.57. Even for the plane shear flows, where shear rate and solid fraction are inhomogeneous, the results of the MKT model agree very well with the DEM results.« less

Automatic network coupling analysis for dynamical systems based on detailed kinetic models.

PubMed

Lebiedz, Dirk; Kammerer, Julia; Brandt-Pollmann, Ulrich

2005-10-01

We introduce a numerical complexity reduction method for the automatic identification and analysis of dynamic network decompositions in (bio)chemical kinetics based on error-controlled computation of a minimal model dimension represented by the number of (locally) active dynamical modes. Our algorithm exploits a generalized sensitivity analysis along state trajectories and subsequent singular value decomposition of sensitivity matrices for the identification of these dominant dynamical modes. It allows for a dynamic coupling analysis of (bio)chemical species in kinetic models that can be exploited for the piecewise computation of a minimal model on small time intervals and offers valuable functional insight into highly nonlinear reaction mechanisms and network dynamics. We present results for the identification of network decompositions in a simple oscillatory chemical reaction, time scale separation based model reduction in a Michaelis-Menten enzyme system and network decomposition of a detailed model for the oscillatory peroxidase-oxidase enzyme system.

A new approach for development of kinetics of wastewater treatment in aerobic biofilm reactor

NASA Astrophysics Data System (ADS)

Goswami, S.; Sarkar, S.; Mazumder, D.

2017-09-01

Biofilm process is widely used for the treatment of a variety of wastewater especially containing slowly biodegradable substances. It provides resistance against toxic environment and is capable of retaining biomass under continuous operation. Development of kinetics is very much pertinent for rational design of a biofilm process for the treatment of wastewater with or without inhibitory substances. A simple approach for development of such kinetics for an aerobic biofilm reactor has been presented using a novel biofilm model. The said biofilm model is formulated from the correlations between substrate concentrations in the influent/effluent and at biofilm liquid interface along with substrate flux and biofilm thickness complying Monod's growth kinetics. The methodology for determining the kinetic coefficients for substrate removal and biomass growth has been demonstrated stepwise along with graphical representations. Kinetic coefficients like K, k, Y, b t, b s, and b d are determined either from the intercepts of X- and Y-axis or from the slope of the graphical plots.

Rapid photogeneration of silver nanoparticles in ethanolic solution: a kinetic study.

PubMed

Yahyaei, Bahareh; Azizian, Saeid

2013-01-15

Ag nanoparticles have been synthesized via UV irradiation of ethanolic solution of AgNO3 in presence of pluronic F127 surfactant. This study is aimed at developing a rapid, simple and green method to prepare Ag nanoparticles and understanding its generation kinetics. The formation dependency of silver nanoparticles on the concentration of reactants, UV exposure time and temperature has been investigated by using UV-vis spectroscopy. The 2D map technique has been used for the first time to estimate the switching time between the nucleation and growth of Ag nanoparticles. Appropriate kinetic models were used for modelling of both stages. Copyright © 2012 Elsevier B.V. All rights reserved.

Stressed Oxidation Life Prediction for C/SiC Composites

NASA Technical Reports Server (NTRS)

Levine, Stanley R.

2004-01-01

The residual strength and life of C/SiC is dominated by carbon interface and fiber oxidation if seal coat and matrix cracks are open to allow oxygen ingress. Crack opening is determined by the combination of thermal, mechanical and thermal expansion mismatch induced stresses. When cracks are open, life can be predicted by simple oxidation based models with reaction controlled kinetics at low temperature, and by gas phase diffusion controlled kinetics at high temperatures. Key life governing variables in these models include temperature, stress, initial strength, oxygen partial pressure, and total pressure. These models are described in this paper.

Stability of Retained Austenite in High-Al, Low-Si TRIP-Assisted Steels Processed via Continuous Galvanizing Heat Treatments

NASA Astrophysics Data System (ADS)

McDermid, J. R.; Zurob, H. S.; Bian, Y.

2011-12-01

Two galvanizable high-Al, low-Si transformation-induced plasticity (TRIP)-assisted steels were subjected to isothermal bainitic transformation (IBT) temperatures compatible with the continuous galvanizing (CGL) process and the kinetics of the retained austenite (RA) to martensite transformation during room temperature deformation studied as a function of heat treatment parameters. It was determined that there was a direct relationship between the rate of strain-induced transformation and optimal mechanical properties, with more gradual transformation rates being favored. The RA to martensite transformation kinetics were successfully modeled using two methodologies: (1) the strain-based model of Olsen and Cohen and (2) a simple relationship with the normalized flow stress, ( {{{σ_{{flow}} - σ_{YS} }/{σ_{YS }}}} ) . For the strain-based model, it was determined that the model parameters were a strong function of strain and alloy thermal processing history and a weak function of alloy chemistry. It was verified that the strain-based model in the present work agrees well with those derived by previous workers using TRIP-assisted steels of similar composition. It was further determined that the RA to martensite transformation kinetics for all alloys and heat treatments could be described using a simple model vs the normalized flow stress, indicating that the RA to martensite transformation is stress-induced rather than strain-induced for temperatures above the Ms^{σ }.

Generalized first-order kinetic model for biosolids decomposition and oxidation during hydrothermal treatment.

PubMed

Shanableh, A

2005-01-01

The main objective of this study was to develop generalized first-order kinetic models to represent hydrothermal decomposition and oxidation of biosolids within a wide range of temperatures (200-450 degrees C). A lumping approach was used in which oxidation of the various organic ingredients was characterized by the chemical oxygen demand (COD), and decomposition was characterized by the particulate (i.e., nonfilterable) chemical oxygen demand (PCOD). Using the Arrhenius equation (k = k(o)e(-Ea/RT)), activation energy (Ea) levels were derived from 42 continuous-flow hydrothermal treatment experiments conducted at temperatures in the range of 200-450 degrees C. Using predetermined values for k(o) in the Arrhenius equation, the activation energies of the various organic ingredients were separated into 42 values for oxidation and a similar number for decomposition. The activation energy values were then classified into levels representing the relative ease at which the organic ingredients of the biosolids were oxidized or decomposed. The resulting simple first-order kinetic models adequately represented, within the experimental data range, hydrothermal decomposition of the organic particles as measured by PCOD and oxidation of the organic content as measured by COD. The modeling approach presented in the paper provide a simple and general framework suitable for assessing the relative reaction rates of the various organic ingredients of biosolids.

Influence of temperature on the hydrolysis, acidogenesis and methanogenesis in mesophilic anaerobic digestion: parameter identification and modeling application.

PubMed

Donoso-Bravo, A; Retamal, C; Carballa, M; Ruiz-Filippi, G; Chamy, R

2009-01-01

The effect of temperature on the kinetic parameters involved in the main reactions of the anaerobic digestion process was studied. Batch tests with starch, glucose and acetic acid as substrates for hydrolysis, acidogenesis and methanogenesis, respectively, were performed in a temperature range between 15 and 45 degrees C. First order kinetics was assumed to determine the hydrolysis rate constant, while Monod and Haldane kinetics were considered for acidogenesis and methanogenesis, respectively. The results obtained showed that the anaerobic process is strongly influenced by temperature, with acidogenesis exerting the highest effect. The Cardinal Temperature Model 1 with an inflection point (CTM1) fitted properly the experimental data in the whole temperature range, except for the maximum degradation rate of acidogenesis. A simple case-study assessing the effect of temperature on an anaerobic CSTR performance indicated that with relatively simple substrates, like starch, the limiting reaction would change depending on temperature. However, when more complex substrates are used (e.g. sewage sludge), the hydrolysis might become more quickly into the limiting step.

Smouldering (thermal) remediation of soil contaminated with industrial organic liquids: novel insights into heat transfer and kinetics uncovered by integrating experiments and modelling

NASA Astrophysics Data System (ADS)

Gerhard, J.; Zanoni, M. A. B.; Torero, J. L.

2017-12-01

Smouldering (i.e., flameless combustion) underpins the technology Self-sustaining Treatment for Active Remediation (STAR). STAR achieves the in situ destruction of nonaqueous phase liquids (NAPLs) by generating a self-sustained smouldering reaction that propagates through the source zone. This research explores the nature of the travelling reaction and the influence of key in situ and engineered characteristics. A novel one-dimensional numerical model was developed (in COMSOL) to simulate the smouldering remediation of bitumen-contaminated sand. This model was validated against laboratory column experiments. Achieving model validation depended on correctly simulating the energy balance at the reaction front, including properly accounting for heat transfer, smouldering kinetics, and heat losses. Heat transfer between soil and air was demonstrated to be generally not at equilibrium. Moreover, existing heat transfer correlations were found to be inappropriate for the low air flow Reynold's numbers (Re < 30) relevant in this and similar thermal remediation systems. Therefore, a suite of experiments were conducted to generate a new heat transfer correlation, which generated correct simulations of convective heat flow through soil. Moreover, it was found that, for most cases of interest, a simple two-step pyrolysis/oxidation set of kinetic reactions was sufficient. Arrhenius parameters, calculated independently from thermogravimetric experiments, allowed the reaction kinetics to be validated in the smouldering model. Furthermore, a simple heat loss term sufficiently accounted for radial heat losses from the column. Altogether, these advances allow this simple model to reasonably predict the self-sustaining process including the peak reaction temperature, the reaction velocity, and the complete destruction of bitumen behind the front. Simulations with the validated model revealed numerous unique insights, including how the system inherently recycles energy, how air flow rate and NAPL saturation dictate contaminant destruction rates, and the extremes that lead to extinction. Overall, this research provides unique insights into the complex interplay of thermochemical processes that govern the success of smouldering as well as other thermal remediation approaches.

Hindered rotor models with variable kinetic functions for accurate thermodynamic and kinetic predictions

NASA Astrophysics Data System (ADS)

Reinisch, Guillaume; Leyssale, Jean-Marc; Vignoles, Gérard L.

2010-10-01

We present an extension of some popular hindered rotor (HR) models, namely, the one-dimensional HR (1DHR) and the degenerated two-dimensional HR (d2DHR) models, allowing for a simple and accurate treatment of internal rotations. This extension, based on the use of a variable kinetic function in the Hamiltonian instead of a constant reduced moment of inertia, is extremely suitable in the case of rocking/wagging motions involved in dissociation or atom transfer reactions. The variable kinetic function is first introduced in the framework of a classical 1DHR model. Then, an effective temperature and potential dependent constant is proposed in the cases of quantum 1DHR and classical d2DHR models. These methods are finally applied to the atom transfer reaction SiCl3+BCl3→SiCl4+BCl2. We show, for this particular case, that a proper accounting of internal rotations greatly improves the accuracy of thermodynamic and kinetic predictions. Moreover, our results confirm (i) that using a suitably defined kinetic function appears to be very adapted to such problems; (ii) that the separability assumption of independent rotations seems justified; and (iii) that a quantum mechanical treatment is not a substantial improvement with respect to a classical one.

The problem of the second wind turbine - a note on a common but flawed wind power estimation method

NASA Astrophysics Data System (ADS)

Gans, F.; Miller, L. M.; Kleidon, A.

2012-06-01

Several recent wind power estimates suggest that this renewable energy resource can meet all of the current and future global energy demand with little impact on the atmosphere. These estimates are calculated using observed wind speeds in combination with specifications of wind turbine size and density to quantify the extractable wind power. However, this approach neglects the effects of momentum extraction by the turbines on the atmospheric flow that would have effects outside the turbine wake. Here we show with a simple momentum balance model of the atmospheric boundary layer that this common methodology to derive wind power potentials requires unrealistically high increases in the generation of kinetic energy by the atmosphere. This increase by an order of magnitude is needed to ensure momentum conservation in the atmospheric boundary layer. In the context of this simple model, we then compare the effect of three different assumptions regarding the boundary conditions at the top of the boundary layer, with prescribed hub height velocity, momentum transport, or kinetic energy transfer into the boundary layer. We then use simulations with an atmospheric general circulation model that explicitly simulate generation of kinetic energy with momentum conservation. These simulations show that the assumption of prescribed momentum import into the atmospheric boundary layer yields the most realistic behavior of the simple model, while the assumption of prescribed hub height velocity can clearly be disregarded. We also show that the assumptions yield similar estimates for extracted wind power when less than 10% of the kinetic energy flux in the boundary layer is extracted by the turbines. We conclude that the common method significantly overestimates wind power potentials by an order of magnitude in the limit of high wind power extraction. Ultimately, environmental constraints set the upper limit on wind power potential at larger scales rather than detailed engineering specifications of wind turbine design and placement.

Glassy behaviour in simple kinetically constrained models: topological networks, lattice analogues and annihilation-diffusion

NASA Astrophysics Data System (ADS)

Sherrington, David; Davison, Lexie; Buhot, Arnaud; Garrahan, Juan P.

2002-02-01

We report a study of a series of simple model systems with only non-interacting Hamiltonians, and hence simple equilibrium thermodynamics, but with constrained dynamics of a type initially suggested by foams and idealized covalent glasses. We demonstrate that macroscopic dynamical features characteristic of real and more complex model glasses, such as two-time decays in energy and auto-correlation functions, arise from the dynamics and we explain them qualitatively and quantitatively in terms of annihilation-diffusion concepts and theory. The comparison is with strong glasses. We also consider fluctuation-dissipation relations and demonstrate subtleties of interpretation. We find no FDT breakdown when the correct normalization is chosen.

Relationship between a solar drying model of red pepper and the kinetics of pure water evaporation (2)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Passamai, V.; Saravia, L.

1997-05-01

In part one, a simple drying model of red pepper related to water evaporation was developed. In this second part the drying model is applied by means of related experiments. Both laboratory and open air drying experiments were carried out to validate the model and simulation results are presented.

Computer-Aided Construction of Chemical Kinetic Models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Green, William H.

2014-12-31

The combustion chemistry of even simple fuels can be extremely complex, involving hundreds or thousands of kinetically significant species. The most reasonable way to deal with this complexity is to use a computer not only to numerically solve the kinetic model, but also to construct the kinetic model in the first place. Because these large models contain so many numerical parameters (e.g. rate coefficients, thermochemistry) one never has sufficient data to uniquely determine them all experimentally. Instead one must work in “predictive” mode, using theoretical rather than experimental values for many of the numbers in the model, and as appropriatemore » refining the most sensitive numbers through experiments. Predictive chemical kinetics is exactly what is needed for computer-aided design of combustion systems based on proposed alternative fuels, particularly for early assessment of the value and viability of proposed new fuels before those fuels are commercially available. This project was aimed at making accurate predictive chemical kinetics practical; this is a challenging goal which requires a range of science advances. The project spanned a wide range from quantum chemical calculations on individual molecules and elementary-step reactions, through the development of improved rate/thermo calculation procedures, the creation of algorithms and software for constructing and solving kinetic simulations, the invention of methods for model-reduction while maintaining error control, and finally comparisons with experiment. Many of the parameters in the models were derived from quantum chemistry calculations, and the models were compared with experimental data measured in our lab or in collaboration with others.« less

Kinetic modeling and fitting software for interconnected reaction schemes: VisKin.

PubMed

Zhang, Xuan; Andrews, Jared N; Pedersen, Steen E

2007-02-15

Reaction kinetics for complex, highly interconnected kinetic schemes are modeled using analytical solutions to a system of ordinary differential equations. The algorithm employs standard linear algebra methods that are implemented using MatLab functions in a Visual Basic interface. A graphical user interface for simple entry of reaction schemes facilitates comparison of a variety of reaction schemes. To ensure microscopic balance, graph theory algorithms are used to determine violations of thermodynamic cycle constraints. Analytical solutions based on linear differential equations result in fast comparisons of first order kinetic rates and amplitudes as a function of changing ligand concentrations. For analysis of higher order kinetics, we also implemented a solution using numerical integration. To determine rate constants from experimental data, fitting algorithms that adjust rate constants to fit the model to imported data were implemented using the Levenberg-Marquardt algorithm or using Broyden-Fletcher-Goldfarb-Shanno methods. We have included the ability to carry out global fitting of data sets obtained at varying ligand concentrations. These tools are combined in a single package, which we have dubbed VisKin, to guide and analyze kinetic experiments. The software is available online for use on PCs.

Mathematical modeling of enzyme production using Trichoderma harzianum P49P11 and sugarcane bagasse as carbon source.

PubMed

Gelain, Lucas; da Cruz Pradella, José Geraldo; da Costa, Aline Carvalho

2015-12-01

A mathematical model to describe the kinetics of enzyme production by the filamentous fungus Trichoderma harzianum P49P11 was developed using a low cost substrate as main carbon source (pretreated sugarcane bagasse). The model describes the cell growth, variation of substrate concentration and production of three kinds of enzymes (cellulases, beta-glucosidase and xylanase) in different sugarcane bagasse concentrations (5; 10; 20; 30; 40 gL(-1)). The 10 gL(-1) concentration was used to validate the model and the other to parameter estimation. The model for enzyme production has terms implicitly representing induction and repression. Substrate variation was represented by a simple degradation rate. The models seem to represent well the kinetics with a good fit for the majority of the assays. Validation results indicate that the models are adequate to represent the kinetics for a biotechnological process. Copyright © 2015 Elsevier Ltd. All rights reserved.

Boltzmann's "H"-Theorem and the Assumption of Molecular Chaos

ERIC Educational Resources Information Center

Boozer, A. D.

2011-01-01

We describe a simple dynamical model of a one-dimensional ideal gas and use computer simulations of the model to illustrate two fundamental results of kinetic theory: the Boltzmann transport equation and the Boltzmann "H"-theorem. Although the model is time-reversal invariant, both results predict that the behaviour of the gas is time-asymmetric.…

Determining Kinetic Parameters for Isothermal Crystallization of Glasses

NASA Technical Reports Server (NTRS)

Ray, C. S.; Zhang, T.; Reis, S. T.; Brow, R. K.

2006-01-01

Non-isothermal crystallization techniques are frequently used to determine the kinetic parameters for crystallization in glasses. These techniques are experimentally simple and quick compared to the isothermal techniques. However, the analytical models used for non-isothermal data analysis, originally developed for describing isothermal transformation kinetics, are fundamentally flawed. The present paper describes a technique for determining the kinetic parameters for isothermal crystallization in glasses, which eliminates most of the common problems that generally make the studies of isothermal crystallization laborious and time consuming. In this technique, the volume fraction of glass that is crystallized as a function of time during an isothermal hold was determined using differential thermal analysis (DTA). The crystallization parameters for the lithium-disilicate (Li2O.2SiO2) model glass were first determined and compared to the same parameters determined by other techniques to establish the accuracy and usefulness of the present technique. This technique was then used to describe the crystallization kinetics of a complex Ca-Sr-Zn-silicate glass developed for sealing solid oxide fuel cells.

Predicting the impact of biocorona formation kinetics on interspecies extrapolations of nanoparticle biodistribution modeling.

PubMed

Sahneh, Faryad Darabi; Scoglio, Caterina M; Monteiro-Riviere, Nancy A; Riviere, Jim E

2015-01-01

To assess the impact of biocorona kinetics on expected tissue distribution of nanoparticles (NPs) across species. The potential fate of NPs in vivo is described through a simple and descriptive pharmacokinetic model using rate processes dependent upon basal metabolic rate coupled to dynamics of protein corona. Mismatch of time scales between interspecies allometric scaling and the kinetics of corona formation is potentially a fundamental issue with interspecies extrapolations of NP biodistribution. The impact of corona evolution on NP biodistribution across two species is maximal when corona transition half-life is close to the geometric mean of NP half-lives of the two species. While engineered NPs can successfully reach target cells in rodent models, the results may be different in humans due to the fact that the longer circulation time allows for further biocorona evolution.

Kinetics of hydrogen peroxide decomposition by catalase: hydroxylic solvent effects.

PubMed

Raducan, Adina; Cantemir, Anca Ruxandra; Puiu, Mihaela; Oancea, Dumitru

2012-11-01

The effect of water-alcohol (methanol, ethanol, propan-1-ol, propan-2-ol, ethane-1,2-diol and propane-1,2,3-triol) binary mixtures on the kinetics of hydrogen peroxide decomposition in the presence of bovine liver catalase is investigated. In all solvents, the activity of catalase is smaller than in water. The results are discussed on the basis of a simple kinetic model. The kinetic constants for product formation through enzyme-substrate complex decomposition and for inactivation of catalase are estimated. The organic solvents are characterized by several physical properties: dielectric constant (D), hydrophobicity (log P), concentration of hydroxyl groups ([OH]), polarizability (α), Kamlet-Taft parameter (β) and Kosower parameter (Z). The relationships between the initial rate, kinetic constants and medium properties are analyzed by linear and multiple linear regression.

Solute transport with multiple equilibrium-controlled or kinetically controlled chemical reactions

USGS Publications Warehouse

Friedly, John C.; Rubin, Jacob

1992-01-01

A new approach is applied to the problem of modeling solute transport accompanied by many chemical reactions. The approach, based on concepts of the concentration space and its stoichiometric subspaces, uses elements of the subspaces as primary dependent variables. It is shown that the resulting model equations are compact in form, isolate the chemical reaction expressions from flow expressions, and can be used for either equilibrium or kinetically controlled reactions. The implications of the results on numerical algorithms for solving the equations are discussed. The application of the theory is illustrated throughout with examples involving a simple but broadly representative set of reactions previously considered in the literature. Numerical results are presented for four interconnected reactions: a homogeneous complexation reaction, two sorption reactions, and a dissolution/precipitation reaction. Three cases are considered: (1) four kinetically controlled reactions, (2) four equilibrium-controlled reactions, and (3) a system with two kinetically controlled reactions and two equilibrium-controlled reactions.

Measurement techniques for analysis of fission fragment excited gases

NASA Technical Reports Server (NTRS)

Schneider, R. T.; Carroll, E. E.; Davis, J. F.; Davie, R. N.; Maguire, T. C.; Shipman, R. G.

1976-01-01

Spectroscopic analysis of fission fragment excited He, Ar, Xe, N2, Ne, Ar-N2, and Ne-N2 have been conducted. Boltzmann plot analysis of He, Ar and Xe have indicated a nonequilibrium, recombining plasma, and population inversions have been found in these gases. The observed radiating species in helium have been adequately described by a simple kinetic model. A more extensive model for argon, nitrogen and Ar-N2 mixtures was developed which adequately describes the energy flow in the system and compares favorably with experimental measurements. The kinetic processes involved in these systems are discussed.

Jordan frame supergravity and inflation in the NMSSM

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ferrara, Sergio; INFN - Laboratori Nazionali di Frascati, Via Enrico Fermi 40, 00044 Frascati; Kallosh, Renata

2010-08-15

We present a complete explicit N=1, d=4 supergravity action in an arbitrary Jordan frame with nonminimal scalar-curvature coupling of the form {Phi}(z,z)R. The action is derived by suitably gauge fixing the superconformal action. The theory has a modified Kaehler geometry, and it exhibits a significant dependence on the frame function {Phi}(z,z) and its derivatives over scalars, in the bosonic as well as in the fermionic part of the action. Under certain simple conditions, the scalar kinetic terms in the Jordan frame have a canonical form. We consider an embedding of the next-to-minimal supersymmetric standard model (NMSSM) gauge theory into supergravity,more » clarifying the Higgs inflation model recently proposed by Einhorn and Jones. We find that the conditions for canonical kinetic terms are satisfied for the NMSSM scalars in the Jordan frame, which leads to a simple action. However, we find that the gauge singlet field experiences a strong tachyonic instability during inflation in this model. Thus, a modification of the model is required to support the Higgs-type inflation.« less

Limiting Energy Dissipation Induces Glassy Kinetics in Single-Cell High-Precision Responses

PubMed Central

Das, Jayajit

2016-01-01

Single cells often generate precise responses by involving dissipative out-of-thermodynamic-equilibrium processes in signaling networks. The available free energy to fuel these processes could become limited depending on the metabolic state of an individual cell. How does limiting dissipation affect the kinetics of high-precision responses in single cells? I address this question in the context of a kinetic proofreading scheme used in a simple model of early-time T cell signaling. Using exact analytical calculations and numerical simulations, I show that limiting dissipation qualitatively changes the kinetics in single cells marked by emergence of slow kinetics, large cell-to-cell variations of copy numbers, temporally correlated stochastic events (dynamic facilitation), and ergodicity breaking. Thus, constraints in energy dissipation, in addition to negatively affecting ligand discrimination in T cells, can create a fundamental difficulty in determining single-cell kinetics from cell-population results. PMID:26958894

Characterizing acid diffusion lengths in chemically amplified resists from measurements of deprotection kinetics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Patil, Abhijit A.; Pandey, Yogendra Narayan; Doxastakis, Manolis

2014-10-01

The acid-catalyzed deprotection of glassy poly(4-hydroxystyrene-co-tertbutyl acrylate) films was studied with infrared absorbance spectroscopy and stochastic simulations. Experimental data were interpreted with a simple description of subdiffusive acid transport coupled to second-order acid loss. This model predicts key attributes of observed deprotection rates, such as fast reaction at short times, slow reaction at long times, and a nonlinear dependence on acid loading. Fickian diffusion is approached by increasing the post-exposure bake temperature or adding plasticizing agents to the polymer resin. These findings demonstrate that acid mobility and overall deprotection kinetics are coupled to glassy matrix dynamics. To complement the analysismore » of bulk kinetics, acid diffusion lengths were calculated from the anomalous transport model and compared with nanopattern line widths. The consistent scaling between experiments and simulations suggests that the anomalous diffusion model could be further developed into a predictive lithography tool.« less

Kinetic energy of Ps formed by Ore mechanism in Ar gas

NASA Astrophysics Data System (ADS)

Sano, Yosuke; Kino, Yasushi; Oka, Toshitaka; Sekine, Tsutomu

2015-06-01

In order to investigate kinetic energy of positronium(Ps) formed by Ore mechanism, we performed positron annihilation age-momentum correlation (AMOC) measurements in Argas for 5.0 MPa and 7.5 MPa at room temperature. From the time dependence of Doppler broadening of para-Ps (p-Ps) self-annihilation gramma-ray component, we observed Ps slowing down process. Using a simple slowing down model, we obtained the initial kinetic energy of Ps formed by Ore mechanism and Ps-Armomentum transfer cross section. The initial kinetic energy was 3.9 eV which was higher than the kinetic energy of Ps formed at the upper limit of Ore gap. The momentum transfer cross section was 0.019 ± 0.010 nm2 in between 1 eV and 3.9 eV, and was close to the theoretical calculation.

Bringing metabolic networks to life: convenience rate law and thermodynamic constraints

PubMed Central

Liebermeister, Wolfram; Klipp, Edda

2006-01-01

Background Translating a known metabolic network into a dynamic model requires rate laws for all chemical reactions. The mathematical expressions depend on the underlying enzymatic mechanism; they can become quite involved and may contain a large number of parameters. Rate laws and enzyme parameters are still unknown for most enzymes. Results We introduce a simple and general rate law called "convenience kinetics". It can be derived from a simple random-order enzyme mechanism. Thermodynamic laws can impose dependencies on the kinetic parameters. Hence, to facilitate model fitting and parameter optimisation for large networks, we introduce thermodynamically independent system parameters: their values can be varied independently, without violating thermodynamical constraints. We achieve this by expressing the equilibrium constants either by Gibbs free energies of formation or by a set of independent equilibrium constants. The remaining system parameters are mean turnover rates, generalised Michaelis-Menten constants, and constants for inhibition and activation. All parameters correspond to molecular energies, for instance, binding energies between reactants and enzyme. Conclusion Convenience kinetics can be used to translate a biochemical network – manually or automatically - into a dynamical model with plausible biological properties. It implements enzyme saturation and regulation by activators and inhibitors, covers all possible reaction stoichiometries, and can be specified by a small number of parameters. Its mathematical form makes it especially suitable for parameter estimation and optimisation. Parameter estimates can be easily computed from a least-squares fit to Michaelis-Menten values, turnover rates, equilibrium constants, and other quantities that are routinely measured in enzyme assays and stored in kinetic databases. PMID:17173669

Modelling of Tc migration in an un-oxidized fractured drill core from Äspö, Sweden

NASA Astrophysics Data System (ADS)

Huber, F. M.; Totskiy, Y.; Montoya Garcia, V.; Enzmann, F.; Trumm, M.; Wenka, A.; Geckeis, H.; Schaefer, T.

2015-12-01

The radionuclide retention of redox sensitive radionuclides (e.g. Pu, Np, U, Tc) in crystalline host rock greatly depends on the rock matrix and the rock redox capacity. Preservation of drill cores concerning oxidation is therefore of paramount importance to reliably predict the near-natural radionuclide retention properties. Here, experimental results of HTO and Tc laboratory migration experiments in a naturally single fractured Äspö un-oxidized drill core are modelled using two different 2D models. Both models employ geometrical information obtained by μ-computed tomography (μCT) scanning of the drill core. The models differ in geometrical complexity meaning the first model (PPM-MD) consists of a simple parallel plate with a porous matrix adjacent to the fracture whereas the second model (MPM) uses the mid-plane of the 3D fracture only (no porous matrix). Simulation results show that for higher flow rates (Peclet number > 1), the MPM satisfactorily describes the HTO breakthrough curves (BTC) whereas the PPM-MD model nicely reproduces the HTO BTC for small Pe numbers (<1). These findings clearly highlight the influence of fracture geometry/flow field complexity on solute transport for Pe numbers > 1 and the dominating effect of matrix diffusion for Peclet numbers < 1. Retention of Tc is modelled using a simple Kd-approach in case of the PPM-MD and including 1st order sorptive reduction/desorption kinetics in case of the MPM. Batch determined sorptive reduction/desorption kinetic rates and Kd values for Tc on non-oxidized Äspö diorite are used in the model and compared to best fit values. By this approach, the transferability of kinetic data concerning sorptive reduction determined in static batch experiments to dynamic transport experiments is examined.

Casimir rack and pinion as a miniaturized kinetic energy harvester

NASA Astrophysics Data System (ADS)

Miri, MirFaez; Etesami, Zahra

2016-08-01

We study a nanoscale machine composed of a rack and a pinion with no contact, but intermeshed via the lateral Casimir force. We adopt a simple model for the random velocity of the rack subject to external random forces, namely, a dichotomous noise with zero mean value. We show that the pinion, even when it experiences random thermal torque, can do work against a load. The device thus converts the kinetic energy of the random motions of the rack into useful work.

Kinetic study of Chromium VI adsorption onto palm kernel shell activated carbon

NASA Astrophysics Data System (ADS)

Mohammad, Masita; Sadeghi Louyeh, Shiva; Yaakob, Zahira

2018-04-01

Heavy metal contamination of industrial effluent is one of the significant environmental problems due to their toxicity and its accumulation throughout the food chain. Adsorption is one of the promising methods for removal of heavy metals from aqua solution because of its simple technique, efficient, reliable and low-cost due to the utilization of residue from the agricultural industry. In this study, activated carbon from palm kernel shells has been produced through chemical activation process using zinc chloride as an activating agent and carbonized at 800 °C. Palm kernel shell activated carbon, PAC was assessed for its efficiency to remove Chromium (VI) ions from aqueous solutions through a batch adsorption process. The kinetic mechanisms have been analysed using Lagergren first-order kinetics model, second-order kinetics model and intra-particle diffusion model. The characterizations such as BET surface area, surface morphology, SEM-EDX have been done. The result shows that the activation process by ZnCl2 was successfully improved the porosity and modified the functional group of palm kernel shell. The result shows that the maximum adsorption capacity of Cr is 11.40mg/g at 30ppm initial metal ion concentration and 0.1g/50mL of adsorbent concentration. The adsorption process followed the pseudo second orders kinetic model.

Hydrogen peroxide and caustic soda: Dancing with a dragon while bleaching

Treesearch

Peter W. Hart; Carl Houtman; Kolby Hirth

2013-01-01

When hydrogen peroxide is mixed with caustic soda, an auto-accelerating reaction can lead to generation of significant amounts of heat and oxygen. On the basis of experiments using typical pulp mill process concentration and temperatures, a relatively simple kinetic model has been developed. Evaluation of these model results reveals that hydrogen peroxide-caustic soda...

qPIPSA: Relating enzymatic kinetic parameters and interaction fields

PubMed Central

Gabdoulline, Razif R; Stein, Matthias; Wade, Rebecca C

2007-01-01

Background The simulation of metabolic networks in quantitative systems biology requires the assignment of enzymatic kinetic parameters. Experimentally determined values are often not available and therefore computational methods to estimate these parameters are needed. It is possible to use the three-dimensional structure of an enzyme to perform simulations of a reaction and derive kinetic parameters. However, this is computationally demanding and requires detailed knowledge of the enzyme mechanism. We have therefore sought to develop a general, simple and computationally efficient procedure to relate protein structural information to enzymatic kinetic parameters that allows consistency between the kinetic and structural information to be checked and estimation of kinetic constants for structurally and mechanistically similar enzymes. Results We describe qPIPSA: quantitative Protein Interaction Property Similarity Analysis. In this analysis, molecular interaction fields, for example, electrostatic potentials, are computed from the enzyme structures. Differences in molecular interaction fields between enzymes are then related to the ratios of their kinetic parameters. This procedure can be used to estimate unknown kinetic parameters when enzyme structural information is available and kinetic parameters have been measured for related enzymes or were obtained under different conditions. The detailed interaction of the enzyme with substrate or cofactors is not modeled and is assumed to be similar for all the proteins compared. The protein structure modeling protocol employed ensures that differences between models reflect genuine differences between the protein sequences, rather than random fluctuations in protein structure. Conclusion Provided that the experimental conditions and the protein structural models refer to the same protein state or conformation, correlations between interaction fields and kinetic parameters can be established for sets of related enzymes. Outliers may arise due to variation in the importance of different contributions to the kinetic parameters, such as protein stability and conformational changes. The qPIPSA approach can assist in the validation as well as estimation of kinetic parameters, and provide insights into enzyme mechanism. PMID:17919319

Modeling of Fusarium redolens Dzf2 mycelial growth kinetics and optimal fed-batch fermentation for beauvericin production.

PubMed

Xu, Li-Jian; Liu, Yuan-Shuai; Zhou, Li-Gang; Wu, Jian-Yong

2011-09-01

Beauvericin (BEA) is a cyclic hexadepsipeptide mycotoxin with notable phytotoxic and insecticidal activities. Fusarium redolens Dzf2 is a highly BEA-producing fungus isolated from a medicinal plant. The aim of the current study was to develop a simple and valid kinetic model for F. redolens Dzf2 mycelial growth and the optimal fed-batch operation for efficient BEA production. A modified Monod model with substrate (glucose) and product (BEA) inhibition was constructed based on the culture characteristics of F. redolens Dzf2 mycelia in a liquid medium. Model parameters were derived by simulation of the experimental data from batch culture. The model fitted closely with the experimental data over 20-50 g l(-1) glucose concentration range in batch fermentation. The kinetic model together with the stoichiometric relationships for biomass, substrate and product was applied to predict the optimal feeding scheme for fed-batch fermentation, leading to 54% higher BEA yield (299 mg l(-1)) than in the batch culture (194 mg l(-1)). The modified Monod model incorporating substrate and product inhibition was proven adequate for describing the growth kinetics of F. redolens Dzf2 mycelial culture at suitable but not excessive initial glucose levels in batch and fed-batch cultures.

Multiplexed immunosensing and kinetics monitoring in nanofluidic devices with highly enhanced target capture efficiency

PubMed Central

Lin, Yii-Lih; Huang, Yen-Jun; Teerapanich, Pattamon; Leïchlé, Thierry

2016-01-01

Nanofluidic devices promise high reaction efficiency and fast kinetic responses due to the spatial constriction of transported biomolecules with confined molecular diffusion. However, parallel detection of multiple biomolecules, particularly proteins, in highly confined space remains challenging. This study integrates extended nanofluidics with embedded protein microarray to achieve multiplexed real-time biosensing and kinetics monitoring. Implementation of embedded standard-sized antibody microarray is attained by epoxy-silane surface modification and a room-temperature low-aspect-ratio bonding technique. An effective sample transport is achieved by electrokinetic pumping via electroosmotic flow. Through the nanoslit-based spatial confinement, the antigen-antibody binding reaction is enhanced with ∼100% efficiency and may be directly observed with fluorescence microscopy without the requirement of intermediate washing steps. The image-based data provide numerous spatially distributed reaction kinetic curves and are collectively modeled using a simple one-dimensional convection-reaction model. This study represents an integrated nanofluidic solution for real-time multiplexed immunosensing and kinetics monitoring, starting from device fabrication, protein immobilization, device bonding, sample transport, to data analysis at Péclet number less than 1. PMID:27375819

The study of zinc ions binding to casein.

PubMed

Pomastowski, P; Sprynskyy, M; Buszewski, B

2014-08-01

The presented research was focused on physicochemical study of casein properties and the kinetics of zinc ions binding to the protein. Moreover, a fast and simple method of casein extraction from cow's milk has been proposed. Casein isoforms, zeta potential (ζ) and particle size of the separated caseins were characterized with the use of capillary electrophoresis, zeta potential analysis and field flow fractionation (FFF) technique, respectively. The kinetics of the metal-binding process was investigated in batch adsorption experiments. Intraparticle diffusion model, first-order and zero-order kinetic models were applied to test the kinetic experimental data. Analysis of changes in infrared bands registered for casein before and after zinc binding was also performed. The obtained results showed that the kinetic process of zinc binding to casein is not homogeneous but is expressed with an initial rapid stage with about 70% of zinc ions immobilized by casein and with a much slower second step. Maximum amount of bound zinc in the experimental conditions was 30.04mgZn/g casein. Copyright © 2014 Elsevier B.V. All rights reserved.

Comparisons of dense-plasma-focus kinetic simulations with experimental measurements.

PubMed

Schmidt, A; Link, A; Welch, D; Ellsworth, J; Falabella, S; Tang, V

2014-06-01

Dense-plasma-focus (DPF) Z-pinch devices are sources of copious high-energy electrons and ions, x rays, and neutrons. The mechanisms through which these physically simple devices generate such high-energy beams in a relatively short distance are not fully understood and past optimization efforts of these devices have been largely empirical. Previously we reported on fully kinetic simulations of a DPF and compared them with hybrid and fluid simulations of the same device. Here we present detailed comparisons between fully kinetic simulations and experimental data on a 1.2 kJ DPF with two electrode geometries, including neutron yield and ion beam energy distributions. A more intensive third calculation is presented which examines the effects of a fully detailed pulsed power driver model. We also compare simulated electromagnetic fluctuations with direct measurement of radiofrequency electromagnetic fluctuations in a DPF plasma. These comparisons indicate that the fully kinetic model captures the essential physics of these plasmas with high fidelity, and provide further evidence that anomalous resistivity in the plasma arises due to a kinetic instability near the lower hybrid frequency.

Relaxation processes in a low-order three-dimensional magnetohydrodynamics model

NASA Technical Reports Server (NTRS)

Stribling, Troy; Matthaeus, William H.

1991-01-01

The time asymptotic behavior of a Galerkin model of 3D magnetohydrodynamics (MHD) has been interpreted using the selective decay and dynamic alignment relaxation theories. A large number of simulations has been performed that scan a parameter space defined by the rugged ideal invariants, including energy, cross helicity, and magnetic helicity. It is concluded that time asymptotic state can be interpreted as a relaxation to minimum energy. A simple decay model, based on absolute equilibrium theory, is found to predict a mapping of initial onto time asymptotic states, and to accurately describe the long time behavior of the runs when magnetic helicity is present. Attention is also given to two processes, operating on time scales shorter than selective decay and dynamic alignment, in which the ratio of kinetic to magnetic energy relaxes to values 0(1). The faster of the two processes takes states initially dominant in magnetic energy to a state of near-equipartition between kinetic and magnetic energy through power law growth of kinetic energy. The other process takes states initially dominant in kinetic energy to the near-equipartitioned state through exponential growth of magnetic energy.

Cylindrically symmetric Green's function approach for modeling the crystal growth morphology of ice.

PubMed

Libbrecht, K G

1999-08-01

We describe a front-tracking Green's function approach to modeling cylindrically symmetric crystal growth. This method is simple to implement, and with little computer power can adequately model a wide range of physical situations. We apply the method to modeling the hexagonal prism growth of ice crystals, which is governed primarily by diffusion along with anisotropic surface kinetic processes. From ice crystal growth observations in air, we derive measurements of the kinetic growth coefficients for the basal and prism faces as a function of temperature, for supersaturations near the water saturation level. These measurements are interpreted in the context of a model for the nucleation and growth of ice, in which the growth dynamics are dominated by the structure of a disordered layer on the ice surfaces.

Kinetic model of water disinfection using peracetic acid including synergistic effects.

PubMed

Flores, Marina J; Brandi, Rodolfo J; Cassano, Alberto E; Labas, Marisol D

2016-01-01

The disinfection efficiencies of a commercial mixture of peracetic acid against Escherichia coli were studied in laboratory scale experiments. The joint and separate action of two disinfectant agents, hydrogen peroxide and peracetic acid, were evaluated in order to observe synergistic effects. A kinetic model for each component of the mixture and for the commercial mixture was proposed. Through simple mathematical equations, the model describes different stages of attack by disinfectants during the inactivation process. Based on the experiments and the kinetic parameters obtained, it could be established that the efficiency of hydrogen peroxide was much lower than that of peracetic acid alone. However, the contribution of hydrogen peroxide was very important in the commercial mixture. It should be noted that this improvement occurred only after peracetic acid had initiated the attack on the cell. This synergistic effect was successfully explained by the proposed scheme and was verified by experimental results. Besides providing a clearer mechanistic understanding of water disinfection, such models may improve our ability to design reactors.

Enhanced axion-photon coupling in GUT with hidden photon

NASA Astrophysics Data System (ADS)

Daido, Ryuji; Takahashi, Fuminobu; Yokozaki, Norimi

2018-05-01

We show that the axion coupling to photons can be enhanced in simple models with a single Peccei-Quinn field, if the gauge coupling unification is realized by a large kinetic mixing χ = O (0.1) between hypercharge and unbroken hidden U(1)H. The key observation is that the U(1)H gauge coupling should be rather strong to induce such large kinetic mixing, leading to enhanced contributions of hidden matter fields to the electromagnetic anomaly. We find that the axion-photon coupling is enhanced by about a factor of 10-100 with respect to the GUT-axion models with E / N = 8 / 3.

A Simple Analytical Model for Predicting the Detectable Ion Current in Ion Mobility Spectrometry Using Corona Discharge Ionization Sources

NASA Astrophysics Data System (ADS)

Kirk, Ansgar Thomas; Kobelt, Tim; Spehlbrink, Hauke; Zimmermann, Stefan

2018-05-01

Corona discharge ionization sources are often used in ion mobility spectrometers (IMS) when a non-radioactive ion source with high ion currents is required. Typically, the corona discharge is followed by a reaction region where analyte ions are formed from the reactant ions. In this work, we present a simple yet sufficiently accurate model for predicting the ion current available at the end of this reaction region when operating at reduced pressure as in High Kinetic Energy Ion Mobility Spectrometers (HiKE-IMS) or most IMS-MS instruments. It yields excellent qualitative agreement with measurement results and is even able to calculate the ion current within an error of 15%. Additional interesting findings of this model are the ion current at the end of the reaction region being independent from the ion current generated by the corona discharge and the ion current in High Kinetic Energy Ion Mobility Spectrometers (HiKE-IMS) growing quadratically when scaling down the length of the reaction region. [Figure not available: see fulltext.

A Simple Analytical Model for Predicting the Detectable Ion Current in Ion Mobility Spectrometry Using Corona Discharge Ionization Sources.

PubMed

Kirk, Ansgar Thomas; Kobelt, Tim; Spehlbrink, Hauke; Zimmermann, Stefan

2018-05-08

Corona discharge ionization sources are often used in ion mobility spectrometers (IMS) when a non-radioactive ion source with high ion currents is required. Typically, the corona discharge is followed by a reaction region where analyte ions are formed from the reactant ions. In this work, we present a simple yet sufficiently accurate model for predicting the ion current available at the end of this reaction region when operating at reduced pressure as in High Kinetic Energy Ion Mobility Spectrometers (HiKE-IMS) or most IMS-MS instruments. It yields excellent qualitative agreement with measurement results and is even able to calculate the ion current within an error of 15%. Additional interesting findings of this model are the ion current at the end of the reaction region being independent from the ion current generated by the corona discharge and the ion current in High Kinetic Energy Ion Mobility Spectrometers (HiKE-IMS) growing quadratically when scaling down the length of the reaction region. Graphical Abstract ᅟ.

A Piagetian Learning Cycle for Introductory Chemical Kinetics.

ERIC Educational Resources Information Center

Batt, Russell H.

1980-01-01

Described is a Piagetian learning cycle based on Monte Carlo modeling of several simple reaction mechanisms. Included are descriptions of learning cycle phases (exploration, invention, and discovery) and four BASIC-PLUS computer programs to be used in the explanation of chemical reacting systems. (Author/DS)

Photocatalytic mineralization of commercial herbicides in a pilot-scale solar CPC reactor: photoreactor modeling and reaction kinetics constants independent of radiation field.

PubMed

Colina-Márquez, Jose; Machuca-Martínez, Fiderman; Li Puma, Gianluca

2009-12-01

The six-flux absorption-scattering model (SFM) of the radiation field in the photoreactor, combined with reaction kinetics and fluid-dynamic models, has proved to be suitable to describe the degradation of water pollutants in heterogeneous photocatalytic reactors, combining simplicity and accuracy. In this study, the above approach was extended to model the photocatalytic mineralization of a commercial herbicides mixture (2,4-D, diuron, and ametryne used in Colombian sugar cane crops) in a solar, pilot-scale, compound parabolic collector (CPC) photoreactor using a slurry suspension of TiO(2). The ray-tracing technique was used jointly with the SFM to determine the direction of both the direct and diffuse solar photon fluxes and the spatial profile of the local volumetric rate of photon absorption (LVRPA) in the CPC reactor. Herbicides mineralization kinetics with explicit photon absorption effects were utilized to remove the dependence of the observed rate constants from the reactor geometry and radiation field in the photoreactor. The results showed that the overall model fitted the experimental data of herbicides mineralization in the solar CPC reactor satisfactorily for both cloudy and sunny days. Using the above approach kinetic parameters independent of the radiation field in the reactor can be estimated directly from the results of experiments carried out in a solar CPC reactor. The SFM combined with reaction kinetics and fluid-dynamic models proved to be a simple, but reliable model, for solar photocatalytic applications.

Imperfect dark energy from kinetic gravity braiding

DOE Office of Scientific and Technical Information (OSTI.GOV)

Deffayet, Cédric; Pujolàs, Oriol; Sawicki, Ignacy

2010-10-01

We introduce a large class of scalar-tensor models with interactions containing the second derivatives of the scalar field but not leading to additional degrees of freedom. These models exhibit peculiar features, such as an essential mixing of scalar and tensor kinetic terms, which we have named kinetic braiding. This braiding causes the scalar stress tensor to deviate from the perfect-fluid form. Cosmology in these models possesses a rich phenomenology, even in the limit where the scalar is an exact Goldstone boson. Generically, there are attractor solutions where the scalar monitors the behaviour of external matter. Because of the kinetic braiding,more » the position of the attractor depends both on the form of the Lagrangian and on the external energy density. The late-time asymptotic of these cosmologies is a de Sitter state. The scalar can exhibit phantom behaviour and is able to cross the phantom divide with neither ghosts nor gradient instabilities. These features provide a new class of models for Dark Energy. As an example, we study in detail a simple one-parameter model. The possible observational signatures of this model include a sizeable Early Dark Energy and a specific equation of state evolving into the final de-Sitter state from a healthy phantom regime.« less

Local rules simulation of the kinetics of virus capsid self-assembly.

PubMed

Schwartz, R; Shor, P W; Prevelige, P E; Berger, B

1998-12-01

A computer model is described for studying the kinetics of the self-assembly of icosahedral viral capsids. Solution of this problem is crucial to an understanding of the viral life cycle, which currently cannot be adequately addressed through laboratory techniques. The abstract simulation model employed to address this is based on the local rules theory of. Proc. Natl. Acad. Sci. USA. 91:7732-7736). It is shown that the principle of local rules, generalized with a model of kinetics and other extensions, can be used to simulate complicated problems in self-assembly. This approach allows for a computationally tractable molecular dynamics-like simulation of coat protein interactions while retaining many relevant features of capsid self-assembly. Three simple simulation experiments are presented to illustrate the use of this model. These show the dependence of growth and malformation rates on the energetics of binding interactions, the tolerance of errors in binding positions, and the concentration of subunits in the examples. These experiments demonstrate a tradeoff within the model between growth rate and fidelity of assembly for the three parameters. A detailed discussion of the computational model is also provided.

Molecular-dynamics simulation of mutual diffusion in nonideal liquid mixtures

NASA Astrophysics Data System (ADS)

Rowley, R. L.; Stoker, J. M.; Giles, N. F.

1991-05-01

The mutual-diffusion coefficients, D 12, of n-hexane, n-heptane, and n-octane in chloroform were modeled using equilibrium molecular-dynamics (MD) simulations of simple Lennard-Jones (LJ) fluids. Pure-component LJ parameters were obtained by comparison of simulations to experimental self-diffusion coefficients. While values of “effective” LJ parameters are not expected to simulate accurately diverse thermophysical properties over a wide range of conditions, it was recently shown that effective parameters obtained from pure self-diffusion coefficients can accurately model mutual diffusion in ideal, liquid mixtures. In this work, similar simulations are used to model diffusion in nonideal mixtures. The same combining rules used in the previous study for the cross-interaction parameters were found to be adequate to represent the composition dependence of D 12. The effect of alkane chain length on D 12 is also correctly predicted by the simulations. A commonly used assumption in empirical correlations of D 12, that its kinetic portion is a simple, compositional average of the intradiffusion coefficients, is inconsistent with the simulation results. In fact, the value of the kinetic portion of D 12 was often outside the range of values bracketed by the two intradiffusion coefficients for the nonideal system modeled here.

A surprisingly complex aqueous chemistry of the simplest amino acid. A pulse radiolysis and theoretical study on H/D kinetic isotope effects in the reaction of glycine anions with hydroxyl radicals.

PubMed

Stefanić, I; Ljubić, I; Bonifacić, M; Sabljić, A; Asmus, K-D; Armstrong, D A

2009-04-07

A pulse radiolysis study was carried out of the reaction rate constants and kinetic isotope effects of hydroxyl-radical-induced H/D abstraction from the most-simple alpha-amino acid glycine in its anionic form in water. The rate constants and yields of three predominantly formed radical products, glycyl (NH2-*CH-CO2-), aminomethyl (NH2-*CH2), and aminyl (*NH-CH2-CO2-) radicals, as well as of their partially or fully deuterated analogs, were found to be of comparable magnitude. The primary, secondary, and primary/secondary H/D kinetic isotope effects on the rate constants were determined with respect to each of the three radicals. The unusual variety of products for such an elementary reaction between two small and simple species indicates a complex mechanism with several reactions taking place simultaneously. Thus, a theoretical modeling of the reaction mechanism and kinetics in the gas- and aqueous phase was performed by using the unrestricted density functional theory with the BB1K functional (employing the polarizable continuum model for the aqueous phase), unrestricted coupled cluster UCCSD(T) method, and improved canonical variational theory. Several hydrogen-bonded prereaction complexes and transition states were detected. In particular, the calculations pointed to a significant mechanistic role of the three-electron two-orbital (sigma/sigma* N therefore O) hemibonded prereaction complexes in the aqueous phase. A good agreement with the experimental rate constants and kinetic isotope effects was achieved by downshifting the calculated reaction barriers by 3 kcal mol(-1) and damping the NH(D) stretching frequency by a factor of 0.86.

Simulations of Flame Acceleration and Deflagration-to-Detonation Transitions in Methane-Air Systems

DTIC Science & Technology

2010-03-17

are neglected. 3. Model parameter calibration The one-step Arrhenius kinetics used in this model cannot ex- actly reproduce all properties of laminar...with obstacles are compared to previ- ously reported experimental data. The results obtained using the simple reaction model qualitatively, and in...have taken in developing a multidimensional numerical model to study explosions in large-scale systems containing mixtures of nat- ural gas and air

Modeling the Mousetrap Car

ERIC Educational Resources Information Center

Jumper, William D.

2012-01-01

Many high school and introductory college physics courses make use of mousetrap car projects and competitions as a way of providing an engaging hands-on learning experience incorporating Newton's laws, conversion of potential to kinetic energy, dissipative forces, and rotational mechanics. Presented here is a simple analytical and finite element…

Cluster dynamics and cluster size distributions in systems of self-propelled particles

NASA Astrophysics Data System (ADS)

Peruani, F.; Schimansky-Geier, L.; Bär, M.

2010-12-01

Systems of self-propelled particles (SPP) interacting by a velocity alignment mechanism in the presence of noise exhibit rich clustering dynamics. Often, clusters are responsible for the distribution of (local) information in these systems. Here, we investigate the properties of individual clusters in SPP systems, in particular the asymmetric spreading behavior of clusters with respect to their direction of motion. In addition, we formulate a Smoluchowski-type kinetic model to describe the evolution of the cluster size distribution (CSD). This model predicts the emergence of steady-state CSDs in SPP systems. We test our theoretical predictions in simulations of SPP with nematic interactions and find that our simple kinetic model reproduces qualitatively the transition to aggregation observed in simulations.

Intercomparison of granular stress and turbulence models for unidirectional sheet flow applications

NASA Astrophysics Data System (ADS)

Chauchat, J.; Cheng, Z.; Hsu, T. J.

2016-12-01

The intergranular stresses are one of the key elements in two-phase sediment transport models. There are two main existing approaches, the kinetic theory of granular flows (Jenkins and Hanes, 1998; Hsu et al., 2004) and the phenomenological rheology such as the one proposed by Bagnold (Hanes and Bowen, 1985) or the μ(I) dense granular flow rheology (Revil-Baudard and Chauchat, 2013). Concerning the turbulent Reynolds stress, mixing length and k-ɛ turbulence models have been validated by previous studies (Revil-Baudard and Chauchat, 2013; Hsu et al., 2004). Recently, sedFoam was developed based on kinetic theory of granular flows and k-ɛ turbulence models (Cheng and Hsu, 2014). In this study, we further extended sedFoam by implementing the mixing length and the dense granular flow rheology by following Revil-Baudard and Chauchat (2013). This allows us to objectively compare the different combinations of intergranular stresses (kinetic theory or the dense granular flow rheology) and turbulence models (mixing length or k-ɛ) under unidirectional sheet flow conditions. We found that the calibrated mixing length and k-ɛ models predicts similar velocity and concentration profiles. The differences observed between the kinetic theory and the dense granular flow rheology requires further investigation. In particular, we hypothesize that the extended kinetic theory proposed by Berzi (2011) would probably improve the existing combination of the kinetic theory with a simple Coulomb frictional model in sedFoam. A semi-analytical solution proposed by Berzi and Fraccarollo(2013) for sediment transport rate and sheet layer thickness versus the Shields number is compared with the results obtained by using the dense granular flow rheology and the mixing length model. The results are similar which demonstrate that both the extended kinetic theory and the dense granular flow rheology can be used to model intergranular stresses under sheet flow conditions.

Limiting Energy Dissipation Induces Glassy Kinetics in Single-Cell High-Precision Responses.

PubMed

Das, Jayajit

2016-03-08

Single cells often generate precise responses by involving dissipative out-of-thermodynamic-equilibrium processes in signaling networks. The available free energy to fuel these processes could become limited depending on the metabolic state of an individual cell. How does limiting dissipation affect the kinetics of high-precision responses in single cells? I address this question in the context of a kinetic proofreading scheme used in a simple model of early-time T cell signaling. Using exact analytical calculations and numerical simulations, I show that limiting dissipation qualitatively changes the kinetics in single cells marked by emergence of slow kinetics, large cell-to-cell variations of copy numbers, temporally correlated stochastic events (dynamic facilitation), and ergodicity breaking. Thus, constraints in energy dissipation, in addition to negatively affecting ligand discrimination in T cells, can create a fundamental difficulty in determining single-cell kinetics from cell-population results. Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.

Approximate models for the ion-kinetic regime in inertial-confinement-fusion capsule implosions

DOE PAGES

Hoffman, Nelson M.; Zimmerman, George B.; Molvig, Kim; ...

2015-05-19

“Reduced” (i.e., simplified or approximate) ion-kinetic (RIK) models in radiation-hydrodynamic simulations permit a useful description of inertial-confinement-fusion (ICF) implosions where kinetic deviations from hydrodynamic behavior are important. For implosions in or near the kinetic regime (i.e., when ion mean free paths are comparable to the capsule size), simulations using a RIK model give a detailed picture of the time- and space-dependent structure of imploding capsules, allow an assessment of the relative importance of various kinetic processes during the implosion, enable explanations of past and current observations, and permit predictions of the results of future experiments. The RIK simulation method describedmore » here uses moment-based reduced kinetic models for transport of mass, momentum, and energy by long-mean-free-path ions, a model for the decrease of fusion reactivity owing to the associated modification of the ion distribution function, and a model of hydrodynamic turbulent mixing. Transport models are based on local gradient-diffusion approximations for the transport of moments of the ion distribution functions, with coefficients to impose flux limiting or account for transport modification. After calibration against a reference set of ICF implosions spanning the hydrodynamic-to-kinetic transition, the method has useful, quantifiable predictive ability over a broad range of capsule parameter space. Calibrated RIK simulations show that an important contributor to ion species separation in ICF capsule implosions is the preferential flux of longer-mean-free-path species out of the fuel and into the shell, leaving the fuel relatively enriched in species with shorter mean free paths. Also, the transport of ion thermal energy is enhanced in the kinetic regime, causing the fuel region to have a more uniform, lower ion temperature, extending over a larger volume, than implied by clean simulations. Furthermore, we expect that the success of our simple approach will motivate continued theoretical research into the development of first-principles-based, comprehensive, self-consistent, yet useable models of kinetic multispecies ion behavior in ICF plasmas.« less

Role of spatial inhomogenity in GPCR dimerisation predicted by receptor association-diffusion models

NASA Astrophysics Data System (ADS)

Deshpande, Sneha A.; Pawar, Aiswarya B.; Dighe, Anish; Athale, Chaitanya A.; Sengupta, Durba

2017-06-01

G protein-coupled receptor (GPCR) association is an emerging paradigm with far reaching implications in the regulation of signalling pathways and therapeutic interventions. Recent super resolution microscopy studies have revealed that receptor dimer steady state exhibits sub-second dynamics. In particular the GPCRs, muscarinic acetylcholine receptor M1 (M1MR) and formyl peptide receptor (FPR), have been demonstrated to exhibit a fast association/dissociation kinetics, independent of ligand binding. In this work, we have developed a spatial kinetic Monte Carlo model to investigate receptor homo-dimerisation at a single receptor resolution. Experimentally measured association/dissociation kinetic parameters and diffusion coefficients were used as inputs to the model. To test the effect of membrane spatial heterogeneity on the simulated steady state, simulations were compared to experimental statistics of dimerisation. In the simplest case the receptors are assumed to be diffusing in a spatially homogeneous environment, while spatial heterogeneity is modelled to result from crowding, membrane micro-domains and cytoskeletal compartmentalisation or ‘corrals’. We show that a simple association-diffusion model is sufficient to reproduce M1MR association statistics, but fails to reproduce FPR statistics despite comparable kinetic constants. A parameter sensitivity analysis is required to reproduce the association statistics of FPR. The model reveals the complex interplay between cytoskeletal components and their influence on receptor association kinetics within the features of the membrane landscape. These results constitute an important step towards understanding the factors modulating GPCR organisation.

Kinetic Analysis of Drug Release from Compounded Slow-release Capsules of Liothyronine Sodium (T3).

PubMed

Bakhteyar, Hamid; Cassone, Clayton; Kohan, Hamed Gilzad; Sani, Shabnam N

2017-01-01

The purpose of this study was to formulate extemporaneously compounded Liothyronine Sodium (T3) slow-release capsules and to evaluate their in vitro drug release performance. Twenty-one formulations containing T3 (7.5 µg) with various compositions of two different grades of Methocel E4M and K100M premium (30% to 90%), and/or SimpleCap/Lactose (10% to 70%) were examined. Quality assessment of the capsules was conducted by standard quality control criteria of the United States Pharmacopeia (i.e., weight variation, content uniformity) to ensure their compliance. The dissolution release profile of the formulations was evaluated using United States Pharmacopeia Apparatus type II (paddle method) at a speed of 50 rpm and temperature of 37°C in phosphate buffered saline media ( pH = 7.2 to 7.4). Aliquots from the media were taken periodically up to 24 hours and analyzed using a validated enzyme-linked immunosorbent assay method. The cumulative percentage of drug release for each formulation was fitted to eleven major release kinetic equations to determine the best-fit model of drug release, as well as the mechanism of release. Assay sensitivity was as low as 1 ng/mL and the optimal calibration range was found to be between 0 ng/mL and 7.5 ng/mL, which corresponded well with the average physiological plasma concentrations of T3. Liothyronine sodium with either SimpleCap (100%) or Methocel E4M (100%) exhibited slowrelease kinetic patterns of Peppas and Zero Order, respectively. The formulation with SimpleCap (100%) had a higher percentage of drug release (as compared to 100% Methocel E4M) within the first four hours; this formulation released 80% of the drug within 12 hours when the release was plateaued thereafter. The formulation with 30% Methocel E4M and 70% SimpleCap released 100% of the drug within the initial 12 hours and exhibited a Zero Order slow-release kinetic pattern. In general, the release kinetic rate of the formulations containing Methocel K100M appeared to be slower than Methocel E4M. This alteration may be due to a higher molecular weight and apparent viscosity of Methocel K100M. While most of the formulations were fitted to a slow-release kinetic pattern, several others including Methocel E4M 100%, 30% Methocel E4M+ 70% Simple Cap, 40% Methocel K100M+ 60% SimpleCap, 50% Methocel K100M+ 50% SimpleCap, 30% Methocel E4M+ 70% Lactose, 90% Methocel E4M+ 10% Lactose, 40% Methocel K100M+ 60% Lactose, and 50% Methocel K100M+ 50% Lactose followed an ideal slow-release kinetic pattern of Zero Order or Higuchi. The results of this study successfully demonstrated the optiomal composition of slow-release compounded capsules of T3. Future studies are warranted to evaluate the in vivo performance of the optimal formulations and to establish an in vitro-in vivo correlation. Copyright© by International Journal of Pharmaceutical Compounding, Inc.

Smouldering Combustion of Soil Organic Matter: Inverse Modelling of the Thermal and Oxidative Degradation Kinetics

NASA Astrophysics Data System (ADS)

Huang, Xinyan; Rein, Guillermo

2013-04-01

Smouldering combustion of soil organic matter (SOM) such as peatlands leads to the largest fires on Earth and posses a possible positive feedback mechanism to climate change. In this work, a kinetic model, including 3-step chemical reactions and 1-step water evaporation is proposed to describe drying, pyrolysis and oxidation behaviour of peat. Peat is chosen as the most important type of SOM susceptible to smoudering, and a Chinese boreal peat sample is selected from the literature. A lumped model of mass loss based on four Arrhenius-type reactions is developed to predict its thermal and oxidative degradation under a range of heating rates. A genetic algorithm is used to solve the inverse problem, and find a group of kinetic and stoichiometric parameters for this peat that provides the best match to the thermogravimetric (TG) data from literature. A multi-objective fitness function is defined using the measurements of both mass loss and mass-loss rate in inert and normal atmospheres under a range of heating rates. Piece-wise optimization is conducted to separate the low temperature drying (<450 K) from the higher temperature pyrolysis and oxidation reaction (>450 K). Modelling results shows the proposed 3-step chemistry is the unique simplest scheme to satisfy all given TG data of this particular peat type. Afterward, this kinetic model and its kinetic parameters are incorporated into a simple one-dimensional species model to study the relative position of each reaction inside a smoulder front. Computational results show that the species model agrees with experimental observations. This is the first time that the smouldering kinetics of SOM is explained and predicted, thus helping to understanding this important natural and widespread phenomenon.

λ-Repressor Oligomerization Kinetics at High Concentrations Using Fluorescence Correlation Spectroscopy in Zero-Mode Waveguides

PubMed Central

Samiee, K. T.; Foquet, M.; Guo, L.; Cox, E. C.; Craighead, H. G.

2005-01-01

Fluorescence correlation spectroscopy (FCS) has demonstrated its utility for measuring transport properties and kinetics at low fluorophore concentrations. In this article, we demonstrate that simple optical nanostructures, known as zero-mode waveguides, can be used to significantly reduce the FCS observation volume. This, in turn, allows FCS to be applied to solutions with significantly higher fluorophore concentrations. We derive an empirical FCS model accounting for one-dimensional diffusion in a finite tube with a simple exponential observation profile. This technique is used to measure the oligomerization of the bacteriophage λ repressor protein at micromolar concentrations. The results agree with previous studies utilizing conventional techniques. Additionally, we demonstrate that the zero-mode waveguides can be used to assay biological activity by measuring changes in diffusion constant as a result of ligand binding. PMID:15613638

A simple theory of motor protein kinetics and energetics. II.

PubMed

Qian, H

2000-01-10

A three-state stochastic model of motor protein [Qian, Biophys. Chem. 67 (1997) pp. 263-267] is further developed to illustrate the relationship between the external load on an individual motor protein in aqueous solution with various ATP concentrations and its steady-state velocity. A wide variety of dynamic motor behavior are obtained from this simple model. For the particular case of free-load translocation being the most unfavorable step within the hydrolysis cycle, the load-velocity curve is quasi-linear, V/Vmax = (cF/Fmax-c)/(1-c), in contrast to the hyperbolic relationship proposed by A.V. Hill for macroscopic muscle. Significant deviation from the linearity is expected when the velocity is less than 10% of its maximal (free-load) value--a situation under which the processivity of motor diminishes and experimental observations are less certain. We then investigate the dependence of load-velocity curve on ATP (ADP) concentration. It is shown that the free load Vmax exhibits a Michaelis-Menten like behavior, and the isometric Fmax increases linearly with ln([ATP]/[ADP]). However, the quasi-linear region is independent of the ATP concentration, yielding an apparently ATP-independent maximal force below the true isometric force. Finally, the heat production as a function of ATP concentration and external load are calculated. In simple terms and solved with elementary algebra, the present model provides an integrated picture of biochemical kinetics and mechanical energetics of motor proteins.

Study of the ion kinetic effects in ICF run-away burn using a quasi-1D hybrid model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huang, Cheng -Kun; Molvig, Kim; Albright, Brian James

Here, the loss of fuel ions in the Gamow peak and other kinetic effects related to the α particles during ignition, run-away burn, and disassembly stages of an inertial confinement fusion D-T capsule are investigated with a quasi-1D hybrid volume ignition model that includes kinetic ions, fluid electrons, Planckian radiation photons, and a metallic pusher. The fuel ion loss due to the Knudsen effect at the fuel-pusher interface is accounted for by a local-loss model by with an albedo model for ions returning from the pusher wall. The tail refilling and relaxation of the fuel ion distribution are captured withmore » a nonlinear Fokker-Planck solver. Alpha heating of the fuel ions is modeled kinetically while simple models for finite alpha range and electron heating are used. This dynamical model is benchmarked with a 3 T hydrodynamic burn model employing similar assumptions. For an energetic pusher (~40 kJ) that compresses the fuel to an areal density of ~1.07g/cm 2 at ignition, the simulation shows that the Knudsen effect can substantially limit ion temperature rise in runaway burn. While the final yield decreases modestly from kinetic effects of the α particles, large reduction of the fuel reactivity during ignition and runaway burn may require a higher Knudsen loss rate compared to the rise time of the temperatures above ~25 keV when the broad D-T Gamow peak merges into the bulk Maxwellian distribution.« less

Study of the ion kinetic effects in ICF run-away burn using a quasi-1D hybrid model

DOE PAGES

Huang, Cheng -Kun; Molvig, Kim; Albright, Brian James; ...

2017-02-21

Here, the loss of fuel ions in the Gamow peak and other kinetic effects related to the α particles during ignition, run-away burn, and disassembly stages of an inertial confinement fusion D-T capsule are investigated with a quasi-1D hybrid volume ignition model that includes kinetic ions, fluid electrons, Planckian radiation photons, and a metallic pusher. The fuel ion loss due to the Knudsen effect at the fuel-pusher interface is accounted for by a local-loss model by with an albedo model for ions returning from the pusher wall. The tail refilling and relaxation of the fuel ion distribution are captured withmore » a nonlinear Fokker-Planck solver. Alpha heating of the fuel ions is modeled kinetically while simple models for finite alpha range and electron heating are used. This dynamical model is benchmarked with a 3 T hydrodynamic burn model employing similar assumptions. For an energetic pusher (~40 kJ) that compresses the fuel to an areal density of ~1.07g/cm 2 at ignition, the simulation shows that the Knudsen effect can substantially limit ion temperature rise in runaway burn. While the final yield decreases modestly from kinetic effects of the α particles, large reduction of the fuel reactivity during ignition and runaway burn may require a higher Knudsen loss rate compared to the rise time of the temperatures above ~25 keV when the broad D-T Gamow peak merges into the bulk Maxwellian distribution.« less

EXPOSURE RECONSTRUCTION FOR REDUCING UNCERTAINTY IN RISK ASSESSMENT: EXAMPLE USING MTBE BIOMARKERS AND SIMPLE PHARMACOKINETIC MODEL

EPA Science Inventory

Adverse health risks from environmental agents are generally related to average (long term) exposures. We used results from a series of controlled human exposure tests and classical first order rate kinetics calculations to estimate how well spot measurements of methyl tertiary ...

Alternative Analysis of the Michaelis-Menten Equations

ERIC Educational Resources Information Center

Krogstad, Harald E.; Dawed, Mohammed Yiha; Tegegne, Tadele Tesfa

2011-01-01

Courses in mathematical modelling are always in need of simple, illustrative examples. The Michaelis-Menten reaction kinetics equations have been considered to be a basic example of scaling and singular perturbation. However, the leading order approximations do not easily show the expected behaviour, and this note proposes a different perturbation…

Probabilistic parameter estimation in a 2-step chemical kinetics model for n-dodecane jet autoignition

NASA Astrophysics Data System (ADS)

Hakim, Layal; Lacaze, Guilhem; Khalil, Mohammad; Sargsyan, Khachik; Najm, Habib; Oefelein, Joseph

2018-05-01

This paper demonstrates the development of a simple chemical kinetics model designed for autoignition of n-dodecane in air using Bayesian inference with a model-error representation. The model error, i.e. intrinsic discrepancy from a high-fidelity benchmark model, is represented by allowing additional variability in selected parameters. Subsequently, we quantify predictive uncertainties in the results of autoignition simulations of homogeneous reactors at realistic diesel engine conditions. We demonstrate that these predictive error bars capture model error as well. The uncertainty propagation is performed using non-intrusive spectral projection that can also be used in principle with larger scale computations, such as large eddy simulation. While the present calibration is performed to match a skeletal mechanism, it can be done with equal success using experimental data only (e.g. shock-tube measurements). Since our method captures the error associated with structural model simplifications, we believe that the optimised model could then lead to better qualified predictions of autoignition delay time in high-fidelity large eddy simulations than the existing detailed mechanisms. This methodology provides a way to reduce the cost of reaction kinetics in simulations systematically, while quantifying the accuracy of predictions of important target quantities.

Dynamic modeling of parallel robots for computed-torque control implementation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Codourey, A.

1998-12-01

In recent years, increased interest in parallel robots has been observed. Their control with modern theory, such as the computed-torque method, has, however, been restrained, essentially due to the difficulty in establishing a simple dynamic model that can be calculated in real time. In this paper, a simple method based on the virtual work principle is proposed for modeling parallel robots. The mass matrix of the robot, needed for decoupling control strategies, does not explicitly appear in the formulation; however, it can be computed separately, based on kinetic energy considerations. The method is applied to the DELTA parallel robot, leadingmore » to a very efficient model that has been implemented in a real-time computed-torque control algorithm.« less

Sensitivity of Polar Stratospheric Ozone Loss to Uncertainties in Chemical Reaction Kinetics

NASA Technical Reports Server (NTRS)

Kawa, S. Randolph; Stolarski, Richard S.; Douglass, Anne R.; Newman, Paul A.

2008-01-01

Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspect of out understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to to simulate these process in numerical models of chemistry and transport. These models depend on laboratory-measured kinetic reaction rates and photlysis cross section to simulate molecular interactions. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluation are applied in random combinations. We determine the key reaction and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.

Calculation of Rate Spectra from Noisy Time Series Data

PubMed Central

Voelz, Vincent A.; Pande, Vijay S.

2011-01-01

As the resolution of experiments to measure folding kinetics continues to improve, it has become imperative to avoid bias that may come with fitting data to a predetermined mechanistic model. Towards this end, we present a rate spectrum approach to analyze timescales present in kinetic data. Computing rate spectra of noisy time series data via numerical discrete inverse Laplace transform is an ill-conditioned inverse problem, so a regularization procedure must be used to perform the calculation. Here, we show the results of different regularization procedures applied to noisy multi-exponential and stretched exponential time series, as well as data from time-resolved folding kinetics experiments. In each case, the rate spectrum method recapitulates the relevant distribution of timescales present in the data, with different priors on the rate amplitudes naturally corresponding to common biases toward simple phenomenological models. These results suggest an attractive alternative to the “Occam’s razor” philosophy of simply choosing models with the fewest number of relaxation rates. PMID:22095854

Fully kinetic simulations of dense plasma focus Z-pinch devices.

PubMed

Schmidt, A; Tang, V; Welch, D

2012-11-16

Dense plasma focus Z-pinch devices are sources of copious high energy electrons and ions, x rays, and neutrons. The mechanisms through which these physically simple devices generate such high-energy beams in a relatively short distance are not fully understood. We now have, for the first time, demonstrated a capability to model these plasmas fully kinetically, allowing us to simulate the pinch process at the particle scale. We present here the results of the initial kinetic simulations, which reproduce experimental neutron yields (~10(7)) and high-energy (MeV) beams for the first time. We compare our fluid, hybrid (kinetic ions and fluid electrons), and fully kinetic simulations. Fluid simulations predict no neutrons and do not allow for nonthermal ions, while hybrid simulations underpredict neutron yield by ~100x and exhibit an ion tail that does not exceed 200 keV. Only fully kinetic simulations predict MeV-energy ions and experimental neutron yields. A frequency analysis in a fully kinetic simulation shows plasma fluctuations near the lower hybrid frequency, possibly implicating lower hybrid drift instability as a contributor to anomalous resistivity in the plasma.

Solving Simple Kinetics without Integrals

ERIC Educational Resources Information Center

de la Pen~a, Lisandro Herna´ndez

2016-01-01

The solution of simple kinetic equations is analyzed without referencing any topic from differential equations or integral calculus. Guided by the physical meaning of the rate equation, a systematic procedure is used to generate an approximate solution that converges uniformly to the exact solution in the case of zero, first, and second order…

Kinetic Theory and Simulation of Single-Channel Water Transport

NASA Astrophysics Data System (ADS)

Tajkhorshid, Emad; Zhu, Fangqiang; Schulten, Klaus

Water translocation between various compartments of a system is a fundamental process in biology of all living cells and in a wide variety of technological problems. The process is of interest in different fields of physiology, physical chemistry, and physics, and many scientists have tried to describe the process through physical models. Owing to advances in computer simulation of molecular processes at an atomic level, water transport has been studied in a variety of molecular systems ranging from biological water channels to artificial nanotubes. While simulations have successfully described various kinetic aspects of water transport, offering a simple, unified model to describe trans-channel translocation of water turned out to be a nontrivial task.

The use of dwell time cross-correlation functions to study single-ion channel gating kinetics.

PubMed Central

Ball, F G; Kerry, C J; Ramsey, R L; Sansom, M S; Usherwood, P N

1988-01-01

The derivation of cross-correlation functions from single-channel dwell (open and closed) times is described. Simulation of single-channel data for simple gating models, alongside theoretical treatment, is used to demonstrate the relationship of cross-correlation functions to underlying gating mechanisms. It is shown that time irreversibility of gating kinetics may be revealed in cross-correlation functions. Application of cross-correlation function analysis to data derived from the locust muscle glutamate receptor-channel provides evidence for multiple gateway states and time reversibility of gating. A model for the gating of this channel is used to show the effect of omission of brief channel events on cross-correlation functions. PMID:2462924

Electron heating in a Monte Carlo model of a high Mach number, supercritical, collisionless shock

NASA Technical Reports Server (NTRS)

Ellison, Donald C.; Jones, Frank C.

1987-01-01

Preliminary work in the investigation of electron injection and acceleration at parallel shocks is presented. A simple model of electron heating that is derived from a unified shock model which includes the effects of an electrostatic potential jump is described. The unified shock model provides a kinetic description of the injection and acceleration of ions and a fluid description of electron heating at high Mach number, supercritical, and parallel shocks.

Toxicokinetics of perfluorooctane sulfonate in birds under environmentally realistic exposure conditions and development of a kinetic predictive model.

PubMed

Tarazona, J V; Rodríguez, C; Alonso, E; Sáez, M; González, F; San Andrés, M D; Jiménez, B; San Andrés, M I

2015-01-22

This article describes the toxicokinetics of perfluorooctane sulfonate (PFOS) in birds under low repeated dosing, equivalent to 0.085 μg/kg per day, representing environmentally realistic exposure conditions. The best fitting was provided by a simple pseudo monocompartmental first-order kinetics model, regulated by two rates, with a pseudo first-order dissipation half-life of 230 days, accounting for real elimination as well as binding of PFOS to non-exchangeable structures. The calculated assimilation efficiency was 0.66 with confidence intervals of 0.64 and 0.68. The model calculations confirmed that the measured maximum concentrations were still far from the steady state situation, which for this dose regime, was estimated at a value of about 65 μg PFOS/L serum achieved after a theoretical 210 weeks continuous exposure. The results confirm a very different kinetics than that observed in single-dose experiments confirming clear dose-related differences in apparent elimination rates in birds, as described for humans and monkeys; suggesting that a capacity-limited saturable process should also be considered in the kinetic behavior of PFOS in birds. Pseudo first-order kinetic models are highly convenient and frequently used for predicting bioaccumulation of chemicals in livestock and wildlife; the study suggests that previous bioaccumulation models using half-lives obtained at high doses are expected to underestimate the biomagnification potential of PFOS. The toxicokinetic parameters presented here can be used for higher-tier bioaccumulation estimations of PFOS in chickens and as surrogate values for modeling PFOS kinetics in wild bird species. Copyright © 2014 Elsevier Ireland Ltd. All rights reserved.

An extensible framework for capturing solvent effects in computer generated kinetic models.

PubMed

Jalan, Amrit; West, Richard H; Green, William H

2013-03-14

Detailed kinetic models provide useful mechanistic insight into a chemical system. Manual construction of such models is laborious and error-prone, which has led to the development of automated methods for exploring chemical pathways. These methods rely on fast, high-throughput estimation of species thermochemistry and kinetic parameters. In this paper, we present a methodology for extending automatic mechanism generation to solution phase systems which requires estimation of solvent effects on reaction rates and equilibria. The linear solvation energy relationship (LSER) method of Abraham and co-workers is combined with Mintz correlations to estimate ΔG(solv)°(T) in over 30 solvents using solute descriptors estimated from group additivity. Simple corrections are found to be adequate for the treatment of radical sites, as suggested by comparison with known experimental data. The performance of scaled particle theory expressions for enthalpic-entropic decomposition of ΔG(solv)°(T) is also presented along with the associated computational issues. Similar high-throughput methods for solvent effects on free-radical kinetics are only available for a handful of reactions due to lack of reliable experimental data, and continuum dielectric calculations offer an alternative method for their estimation. For illustration, we model liquid phase oxidation of tetralin in different solvents computing the solvent dependence for ROO• + ROO• and ROO• + solvent reactions using polarizable continuum quantum chemistry methods. The resulting kinetic models show an increase in oxidation rate with solvent polarity, consistent with experiment. Further work needed to make this approach more generally useful is outlined.

Mechanistic study of manganese-substituted glycerol dehydrogenase using a kinetic and thermodynamic analysis.

PubMed

Fang, Baishan; Niu, Jin; Ren, Hong; Guo, Yingxia; Wang, Shizhen

2014-01-01

Mechanistic insights regarding the activity enhancement of dehydrogenase by metal ion substitution were investigated by a simple method using a kinetic and thermodynamic analysis. By profiling the binding energy of both the substrate and product, the metal ion's role in catalysis enhancement was revealed. Glycerol dehydrogenase (GDH) from Klebsiella pneumoniae sp., which demonstrated an improvement in activity by the substitution of a zinc ion with a manganese ion, was used as a model for the mechanistic study of metal ion substitution. A kinetic model based on an ordered Bi-Bi mechanism was proposed considering the noncompetitive product inhibition of dihydroxyacetone (DHA) and the competitive product inhibition of NADH. By obtaining preliminary kinetic parameters of substrate and product inhibition, the number of estimated parameters was reduced from 10 to 4 for a nonlinear regression-based kinetic parameter estimation. The simulated values of time-concentration curves fit the experimental values well, with an average relative error of 11.5% and 12.7% for Mn-GDH and GDH, respectively. A comparison of the binding energy of enzyme ternary complex for Mn-GDH and GDH derived from kinetic parameters indicated that metal ion substitution accelerated the release of dioxyacetone. The metal ion's role in catalysis enhancement was explicated.

Ultrafiltration membrane reactors for enzymatic resolution of amino acids: design model and optimization.

PubMed

Bódalo, A; Gómez, J L.; Gómez, E; Bastida, J; Máximo, M F.; Montiel, M C.

2001-03-08

In this paper the possibility of continuous resolution of DL-phenylalanine, catalyzed by L-aminoacylase in a ultrafiltration membrane reactor (UFMR) is presented. A simple design model, based on previous kinetic studies, has been demonstrated to be capable of describing the behavior of the experimental system. The model has been used to determine the optimal experimental conditions to carry out the asymmetrical hydrolysis of N-acetyl-DL-phenylalanine.

Thermodynamics of information processing based on enzyme kinetics: An exactly solvable model of an information pump.

PubMed

Cao, Yuansheng; Gong, Zongping; Quan, H T

2015-06-01

Motivated by the recent proposed models of the information engine [Proc. Natl. Acad. Sci. USA 109, 11641 (2012)] and the information refrigerator [Phys. Rev. Lett. 111, 030602 (2013)], we propose a minimal model of the information pump and the information eraser based on enzyme kinetics. This device can either pump molecules against the chemical potential gradient by consuming the information to be encoded in the bit stream or (partially) erase the information initially encoded in the bit stream by consuming the Gibbs free energy. The dynamics of this model is solved exactly, and the "phase diagram" of the operation regimes is determined. The efficiency and the power of the information machine is analyzed. The validity of the second law of thermodynamics within our model is clarified. Our model offers a simple paradigm for the investigating of the thermodynamics of information processing involving the chemical potential in small systems.

PyFolding: Open-Source Graphing, Simulation, and Analysis of the Biophysical Properties of Proteins.

PubMed

Lowe, Alan R; Perez-Riba, Albert; Itzhaki, Laura S; Main, Ewan R G

2018-02-06

For many years, curve-fitting software has been heavily utilized to fit simple models to various types of biophysical data. Although such software packages are easy to use for simple functions, they are often expensive and present substantial impediments to applying more complex models or for the analysis of large data sets. One field that is reliant on such data analysis is the thermodynamics and kinetics of protein folding. Over the past decade, increasingly sophisticated analytical models have been generated, but without simple tools to enable routine analysis. Consequently, users have needed to generate their own tools or otherwise find willing collaborators. Here we present PyFolding, a free, open-source, and extensible Python framework for graphing, analysis, and simulation of the biophysical properties of proteins. To demonstrate the utility of PyFolding, we have used it to analyze and model experimental protein folding and thermodynamic data. Examples include: 1) multiphase kinetic folding fitted to linked equations, 2) global fitting of multiple data sets, and 3) analysis of repeat protein thermodynamics with Ising model variants. Moreover, we demonstrate how PyFolding is easily extensible to novel functionality beyond applications in protein folding via the addition of new models. Example scripts to perform these and other operations are supplied with the software, and we encourage users to contribute notebooks and models to create a community resource. Finally, we show that PyFolding can be used in conjunction with Jupyter notebooks as an easy way to share methods and analysis for publication and among research teams. Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.

A simple hydrodynamic model of tornado-like vortices

NASA Astrophysics Data System (ADS)

Kurgansky, M. V.

2015-05-01

Based on similarity arguments, a simple fluid dynamic model of tornado-like vortices is offered that, with account for "vortex breakdown" at a certain height above the ground, relates the maximal azimuthal velocity in the vortex, reachable near the ground surface, to the convective available potential energy (CAPE) stored in the environmental atmosphere under pre-tornado conditions. The relative proportion of the helicity (kinetic energy) destruction (dissipation) in the "vortex breakdown" zone and, accordingly, within the surface boundary layer beneath the vortex is evaluated. These considerations form the basis of the dynamic-statistical analysis of the relationship between the tornado intensity and the CAPE budget in the surrounding atmosphere.

Development of simple-to-apply biogas kinetic models for the co-digestion of food waste and maize husk.

PubMed

Owamah, H I; Izinyon, O C

2015-10-01

Biogas kinetic models are often used to characterize substrate degradation and prediction of biogas production potential. Most of these existing models are however difficult to apply to substrates they were not developed for since their applications are usually substrate specific. Biodegradability kinetic (BIK) model and maximum biogas production potential and stability assessment (MBPPSA) model were therefore developed in this study for better understanding of the anaerobic co-digestion of food waste and maize husk for biogas production. Biodegradability constant (k) was estimated as 0.11 d(-1) using the BIK model. The results of maximum biogas production potential (A) obtained using the MBPPSA model were found to be in good correspondence, both in value and trend with the results obtained using the popular but complex modified Gompertz model for digesters B-1, B-2, B-3, B-4, and B-5. The (If) value of MBPPSA model also showed that digesters B-3, B-4, and B-5 were stable, while B-1 and B-2 were inhibited/unstable. Similar stability observation was also obtained using the modified Gompertz model. The MBPPSA model can therefore be used as an alternative model for anaerobic digestion feasibility studies and plant design. Copyright © 2015 Elsevier Ltd. All rights reserved.

Reverse Stability Kinetics of Meat Pigment Oxidation in Aqueous Extract from Fresh Beef.

PubMed

Frelka, John C; Phinney, David M; Wick, Macdonald P; Heldman, Dennis R

2017-12-01

The use of kinetic models is an evolving approach to describing quality changes in foods during processes, including storage. Previous studies indicate that the oxidation rate of myoglobin is accelerated under frozen storage conditions, a phenomenon termed reverse stability. The goal of this study was to develop a model for meat pigment oxidation to incorporate the phenomenon of reverse stability. In this investigation, the model system was an aqueous extract from beef which was stored under a range of temperatures, both unfrozen and frozen. The kinetic analysis showed that in unfrozen solutions, the temperature dependence of oxidation rate followed Arrhenius kinetics. However, under in frozen solutions the rate of oxidation increased with decreasing temperature until reaching a local maximum around -20 °C. The addition of NaCl to the model system increased oxidation rates at all temperatures, even above the initial freezing temperature. This observation suggests that this reaction is dependent on the ionic strength of the solution as well as temperature. The mechanism of this deviant kinetic behavior is not fully understood, but this study shows that the interplay of temperature and composition on the rate of oxidation of meat pigments is complicated and may involve multiple mechanisms. A better understanding of the kinetics of quality loss in a meat system allows for a re-examination of the current recommendations for frozen storage. The deviant kinetic behavior observed in this study indicates that the relationship between quality loss and temperature in a frozen food is not as simple as once thought. Product-specific recommendations could be implemented in the future that would allow for a decrease in energy consumption without a significant loss of quality. © 2017 Institute of Food Technologists®.

A consistent hierarchy of generalized kinetic equation approximations to the master equation applied to surface catalysis.

PubMed

Herschlag, Gregory J; Mitran, Sorin; Lin, Guang

2015-06-21

We develop a hierarchy of approximations to the master equation for systems that exhibit translational invariance and finite-range spatial correlation. Each approximation within the hierarchy is a set of ordinary differential equations that considers spatial correlations of varying lattice distance; the assumption is that the full system will have finite spatial correlations and thus the behavior of the models within the hierarchy will approach that of the full system. We provide evidence of this convergence in the context of one- and two-dimensional numerical examples. Lower levels within the hierarchy that consider shorter spatial correlations are shown to be up to three orders of magnitude faster than traditional kinetic Monte Carlo methods (KMC) for one-dimensional systems, while predicting similar system dynamics and steady states as KMC methods. We then test the hierarchy on a two-dimensional model for the oxidation of CO on RuO2(110), showing that low-order truncations of the hierarchy efficiently capture the essential system dynamics. By considering sequences of models in the hierarchy that account for longer spatial correlations, successive model predictions may be used to establish empirical approximation of error estimates. The hierarchy may be thought of as a class of generalized phenomenological kinetic models since each element of the hierarchy approximates the master equation and the lowest level in the hierarchy is identical to a simple existing phenomenological kinetic models.

Kinetic analysis of the effects of target structure on siRNA efficiency

NASA Astrophysics Data System (ADS)

Chen, Jiawen; Zhang, Wenbing

2012-12-01

RNAi efficiency for target cleavage and protein expression is related to the target structure. Considering the RNA-induced silencing complex (RISC) as a multiple turnover enzyme, we investigated the effect of target mRNA structure on siRNA efficiency with kinetic analysis. The 4-step model was used to study the target cleavage kinetic process: hybridization nucleation at an accessible target site, RISC-mRNA hybrid elongation along with mRNA target structure melting, target cleavage, and enzyme reactivation. At this model, the terms accounting for the target accessibility, stability, and the seed and the nucleation site effects are all included. The results are in good agreement with that of experiments which show different arguments about the structure effects on siRNA efficiency. It shows that the siRNA efficiency is influenced by the integrated factors of target's accessibility, stability, and the seed effects. To study the off-target effects, a simple model of one siRNA binding to two mRNA targets was designed. By using this model, the possibility for diminishing the off-target effects by the concentration of siRNA was discussed.

A Simple Classroom Teaching Technique to Help Students Understand Michaelis-Menten Kinetics

ERIC Educational Resources Information Center

Runge, Steven W.; Hill, Brent J. F.; Moran, William M.

2006-01-01

A new, simple classroom technique helps cell biology students understand principles of Michaelis-Menten enzyme kinetics. A student mimics the enzyme and the student's hand represents the enzyme's active site. The catalytic event is the transfer of marbles (substrate molecules) by hand from one plastic container to another. As predicted, increases…

Kinetics of rapid covalent bond formation of aniline with humic acid: ESR investigations with nitroxide spin labels

NASA Astrophysics Data System (ADS)

Glinka, Kevin; Matthies, Michael; Theiling, Marius; Hideg, Kalman; Steinhoff, Heinz-Jürgen

2016-04-01

Sulfonamide antibiotics used in livestock farming are distributed to farmland by application of slurry as fertilizer. Previous work suggests rapid covalent binding of the aniline moiety to humic acids found in soil. In the current work, kinetics of this binding were measured in X-band EPR spectroscopy by incubating Leonardite humic acid (LHA) with a paramagnetic aniline spin label (anilino-NO (2,5,5-Trimethyl-2-(3-aminophenyl)pyrrolidin-1-oxyl)). Binding was detected by a pronounced broadening of the spectral lines after incubation of LHA with anilino-NO. The time evolution of the amplitude of this feature was used for determining the reaction kinetics. Single- and double-exponential models were fitted to the data obtained for modelling one or two first-order reactions. Reaction rates of 0.16 min-1 and 0.012 min-1, were found respectively. Addition of laccase peroxidase did not change the kinetics but significantly enhanced the reacting fraction of anilino-NO. This EPR-based method provides a technically simple and effective method for following rapid binding processes of a xenobiotic substance to humic acids.

Kinetic models of gene expression including non-coding RNAs

NASA Astrophysics Data System (ADS)

Zhdanov, Vladimir P.

2011-03-01

In cells, genes are transcribed into mRNAs, and the latter are translated into proteins. Due to the feedbacks between these processes, the kinetics of gene expression may be complex even in the simplest genetic networks. The corresponding models have already been reviewed in the literature. A new avenue in this field is related to the recognition that the conventional scenario of gene expression is fully applicable only to prokaryotes whose genomes consist of tightly packed protein-coding sequences. In eukaryotic cells, in contrast, such sequences are relatively rare, and the rest of the genome includes numerous transcript units representing non-coding RNAs (ncRNAs). During the past decade, it has become clear that such RNAs play a crucial role in gene expression and accordingly influence a multitude of cellular processes both in the normal state and during diseases. The numerous biological functions of ncRNAs are based primarily on their abilities to silence genes via pairing with a target mRNA and subsequently preventing its translation or facilitating degradation of the mRNA-ncRNA complex. Many other abilities of ncRNAs have been discovered as well. Our review is focused on the available kinetic models describing the mRNA, ncRNA and protein interplay. In particular, we systematically present the simplest models without kinetic feedbacks, models containing feedbacks and predicting bistability and oscillations in simple genetic networks, and models describing the effect of ncRNAs on complex genetic networks. Mathematically, the presentation is based primarily on temporal mean-field kinetic equations. The stochastic and spatio-temporal effects are also briefly discussed.

Generalized time-dependent model of radiation-induced chromosomal aberrations in normal and repair-deficient human cells.

PubMed

Ponomarev, Artem L; George, Kerry; Cucinotta, Francis A

2014-03-01

We have developed a model that can simulate the yield of radiation-induced chromosomal aberrations (CAs) and unrejoined chromosome breaks in normal and repair-deficient cells. The model predicts the kinetics of chromosomal aberration formation after exposure in the G₀/G₁ phase of the cell cycle to either low- or high-LET radiation. A previously formulated model based on a stochastic Monte Carlo approach was updated to consider the time dependence of DNA double-strand break (DSB) repair (proper or improper), and different cell types were assigned different kinetics of DSB repair. The distribution of the DSB free ends was derived from a mechanistic model that takes into account the structure of chromatin and DSB clustering from high-LET radiation. The kinetics of chromosomal aberration formation were derived from experimental data on DSB repair kinetics in normal and repair-deficient cell lines. We assessed different types of chromosomal aberrations with the focus on simple and complex exchanges, and predicted the DSB rejoining kinetics and misrepair probabilities for different cell types. The results identify major cell-dependent factors, such as a greater yield of chromosome misrepair in ataxia telangiectasia (AT) cells and slower rejoining in Nijmegen (NBS) cells relative to the wild-type. The model's predictions suggest that two mechanisms could exist for the inefficiency of DSB repair in AT and NBS cells, one that depends on the overall speed of joining (either proper or improper) of DNA broken ends, and another that depends on geometric factors, such as the Euclidian distance between DNA broken ends, which influences the relative frequency of misrepair.

Electrodiffusion kinetics of ionic transport in a simple membrane channel.

PubMed

Valent, Ivan; Petrovič, Pavol; Neogrády, Pavel; Schreiber, Igor; Marek, Miloš

2013-11-21

We employ numerical techniques for solving time-dependent full Poisson-Nernst-Planck (PNP) equations in 2D to analyze transient behavior of a simple ion channel subject to a sudden electric potential jump across the membrane (voltage clamp). Calculated spatiotemporal profiles of the ionic concentrations and electric potential show that two principal exponential processes can be distinguished in the electrodiffusion kinetics, in agreement with original Planck's predictions. The initial fast process corresponds to the dielectric relaxation, while the steady state is approached in a second slower exponential process attributed to the nonlinear ionic redistribution. Effects of the model parameters such as the channel length, height of the potential step, boundary concentrations, permittivity of the channel interior, and ionic mobilities on electrodiffusion kinetics are studied. Numerical solutions are used to determine spatiotemporal profiles of the electric field, ionic fluxes, and both the conductive and displacement currents. We demonstrate that the displacement current is a significant transient component of the total electric current through the channel. The presented results provide additional information about the classical voltage-clamp problem and offer further physical insights into the mechanism of electrodiffusion. The used numerical approach can be readily extended to multi-ionic models with a more structured domain geometry in 2D or 3D, and it is directly applicable to other systems, such as synthetic nanopores, nanofluidic channels, and nanopipettes.

A simple analytical model for signal amplification by reversible exchange (SABRE) process.

PubMed

Barskiy, Danila A; Pravdivtsev, Andrey N; Ivanov, Konstantin L; Kovtunov, Kirill V; Koptyug, Igor V

2016-01-07

We demonstrate an analytical model for the description of the signal amplification by reversible exchange (SABRE) process. The model relies on a combined analysis of chemical kinetics and the evolution of the nuclear spin system during the hyperpolarization process. The presented model for the first time provides rationale for deciding which system parameters (i.e. J-couplings, relaxation rates, reaction rate constants) have to be optimized in order to achieve higher signal enhancement for a substrate of interest in SABRE experiments.

DNA denaturation through a model of the partition points on a one-dimensional lattice

NASA Astrophysics Data System (ADS)

Mejdani, R.; Huseini, H.

1994-08-01

We have shown that by using a model of the partition points gas on a one-dimensional lattice, we can study, besides the saturation curves obtained before for the enzyme kinetics, also the denaturation process, i.e. the breaking of the hydrogen bonds connecting the two strands, under treatment by heat of DNA. We think that this model, as a very simple model and mathematically transparent, can be advantageous for pedagogic goals or other theoretical investigations in chemistry or modern biology.

The physiological kinetics of nitrogen and the prevention of decompression sickness.

PubMed

Doolette, D J; Mitchell, S J

2001-01-01

Decompression sickness (DCS) is a potentially crippling disease caused by intracorporeal bubble formation during or after decompression from a compressed gas underwater dive. Bubbles most commonly evolve from dissolved inert gas accumulated during the exposure to increased ambient pressure. Most diving is performed breathing air, and the inert gas of interest is nitrogen. Divers use algorithms based on nitrogen kinetic models to plan the duration and degree of exposure to increased ambient pressure and to control their ascent rate. However, even correct execution of dives planned using such algorithms often results in bubble formation and may result in DCS. This reflects the importance of idiosyncratic host factors that are difficult to model, and deficiencies in current nitrogen kinetic models. Models describing the exchange of nitrogen between tissues and blood may be based on distributed capillary units or lumped compartments, either of which may be perfusion- or diffusion-limited. However, such simplistic models are usually poor predictors of experimental nitrogen kinetics at the organ or tissue level, probably because they fail to account for factors such as heterogeneity in both tissue composition and blood perfusion and non-capillary exchange mechanisms. The modelling of safe decompression procedures is further complicated by incomplete understanding of the processes that determine bubble formation. Moreover, any formation of bubbles during decompression alters subsequent nitrogen kinetics. Although these factors mandate complex resolutions to account for the interaction between dissolved nitrogen kinetics and bubble formation and growth, most decompression schedules are based on relatively simple perfusion-limited lumped compartment models of blood: tissue nitrogen exchange. Not surprisingly, all models inevitably require empirical adjustment based on outcomes in the field. Improvements in the predictive power of decompression calculations are being achieved using probabilistic bubble models, but divers will always be subject to the possibility of developing DCS despite adherence to prescribed limits.

Compartmental analysis of [11C]flumazenil kinetics for the estimation of ligand transport rate and receptor distribution using positron emission tomography.

PubMed

Koeppe, R A; Holthoff, V A; Frey, K A; Kilbourn, M R; Kuhl, D E

1991-09-01

The in vivo kinetic behavior of [11C]flumazenil ([11C]FMZ), a non-subtype-specific central benzodiazepine antagonist, is characterized using compartmental analysis with the aim of producing an optimized data acquisition protocol and tracer kinetic model configuration for the assessment of [11C]FMZ binding to benzodiazepine receptors (BZRs) in human brain. The approach presented is simple, requiring only a single radioligand injection. Dynamic positron emission tomography data were acquired on 18 normal volunteers using a 60- to 90-min sequence of scans and were analyzed with model configurations that included a three-compartment, four-parameter model, a three-compartment, three-parameter model, with a fixed value for free plus nonspecific binding; and a two-compartment, two-parameter model. Statistical analysis indicated that a four-parameter model did not yield significantly better fits than a three-parameter model. Goodness of fit was improved for three- versus two-parameter configurations in regions with low receptor density, but not in regions with moderate to high receptor density. Thus, a two-compartment, two-parameter configuration was found to adequately describe the kinetic behavior of [11C]FMZ in human brain, with stable estimates of the model parameters obtainable from as little as 20-30 min of data. Pixel-by-pixel analysis yields functional images of transport rate (K1) and ligand distribution volume (DV"), and thus provides independent estimates of ligand delivery and BZR binding.

Accounting for Space—Quantification of Cell-To-Cell Transmission Kinetics Using Virus Dynamics Models.

PubMed

Kumberger, Peter; Durso-Cain, Karina; Uprichard, Susan L; Dahari, Harel; Graw, Frederik

2018-04-17

Mathematical models based on ordinary differential equations (ODE) that describe the population dynamics of viruses and infected cells have been an essential tool to characterize and quantify viral infection dynamics. Although an important aspect of viral infection is the dynamics of viral spread, which includes transmission by cell-free virions and direct cell-to-cell transmission, models used so far ignored cell-to-cell transmission completely, or accounted for this process by simple mass-action kinetics between infected and uninfected cells. In this study, we show that the simple mass-action approach falls short when describing viral spread in a spatially-defined environment. Using simulated data, we present a model extension that allows correct quantification of cell-to-cell transmission dynamics within a monolayer of cells. By considering the decreasing proportion of cells that can contribute to cell-to-cell spread with progressing infection, our extension accounts for the transmission dynamics on a single cell level while still remaining applicable to standard population-based experimental measurements. While the ability to infer the proportion of cells infected by either of the transmission modes depends on the viral diffusion rate, the improved estimates obtained using our novel approach emphasize the need to correctly account for spatial aspects when analyzing viral spread.

Simulations of Quantum Dot Growth on Semiconductor Surfaces: Morphological Design of Sensor Concepts

DTIC Science & Technology

2008-12-01

size equalization can be clearly illustrated during the growth process. In this work we develop a fast multiscale 3D kinetic Monte Carlo ( KMC ) QD...model will provide an attractive means for producing predictably ordered nanostructures. MODEL DESCRIPTION The 3D layer-by-layer KMC growth model...Voter, 2001) and KMC simulation experience (Pan et al., 2004; Pan et al., 2006; Meixner et al, 2003) in 2D, we therefore propose the following simple

Modified ADM1 for modeling free ammonia inhibition in anaerobic acidogenic fermentation with high-solid sludge.

PubMed

Bai, Jie; Liu, He; Yin, Bo; Ma, Huijun; Chen, Xinchun

2017-02-01

Anaerobic acidogenic fermentation with high-solid sludge is a promising method for volatile fatty acid (VFA) production to realize resource recovery. In this study, to model inhibition by free ammonia in high-solid sludge fermentation, the anaerobic digestion model No. 1 (ADM1) was modified to simulate the VFA generation in batch, semi-continuous and full scale sludge. The ADM1 was operated on the platform AQUASIM 2.0. Three kinds of inhibition forms, e.g., simple inhibition, Monod and non-inhibition forms, were integrated into the ADM1 and tested with the real experimental data for batch and semi-continuous fermentation, respectively. The improved particle swarm optimization technique was used for kinetic parameter estimation using the software MATLAB 7.0. In the modified ADM1, the K s of acetate is 0.025, the k m,ac is 12.51, and the K I_NH3 is 0.02, respectively. The results showed that the simple inhibition model could simulate the VFA generation accurately while the Monod model was the better inhibition kinetics form in semi-continuous fermentation at pH10.0. Finally, the modified ADM1 could successfully describe the VFA generation and ammonia accumulation in a 30m 3 full-scale sludge fermentation reactor, indicating that the developed model can be applicable in high-solid sludge anaerobic fermentation. Copyright © 2016. Published by Elsevier B.V.

Kinetics and thermodynamics associated with Bi adsorption transitions at Cu and Ni grain boundaries

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tai, Kaiping; Feng, Lin; Dillon, Shen J.

The grain boundary diffusivity of Au in Cu and Cu-Bi, and Cu in Ni and Ni-Bi are characterized by secondary ion mass spectroscopy depth profiling. Samples are equilibrated in a Bi containing atmosphere at temperatures above and below the onset of grain boundary adsorption transitions, sometimes called complexion transitions. A simple thermo-kinetic model is used to estimate the relative entropic contributions to the grain boundary energies. The results indicate that the entropy term plays a major role in promoting thermally and chemically induced grain boundary complexion transition.

Evolution of large amplitude Alfven waves in solar wind plasmas: Kinetic-fluid models

NASA Astrophysics Data System (ADS)

Nariyuki, Y.

2014-12-01

Large amplitude Alfven waves are ubiquitously observed in solar wind plasmas. Mjolhus(JPP, 1976) and Mio et al(JPSJ, 1976) found that nonlinear evolution of the uni-directional, parallel propagating Alfven waves can be described by the derivative nonlinear Schrodinger equation (DNLS). Later, the multi-dimensional extension (Mjolhus and Wyller, JPP, 1988; Passot and Sulem, POP, 1993; Gazol et al, POP, 1999) and ion kinetic modification (Mjolhus and Wyller, JPP, 1988; Spangler, POP, 1989; Medvedev and Diamond, POP, 1996; Nariyuki et al, POP, 2013) of DNLS have been reported. Recently, Nariyuki derived multi-dimensional DNLS from an expanding box model of the Hall-MHD system (Nariyuki, submitted). The set of equations including the nonlinear evolution of compressional wave modes (TDNLS) was derived by Hada(GRL, 1993). DNLS can be derived from TDNLS by rescaling of the variables (Mjolhus, Phys. Scr., 2006). Nariyuki and Hada(JPSJ, 2007) derived a kinetically modified TDNLS by using a simple Landau closure (Hammet and Perkins, PRL, 1990; Medvedev and Diamond, POP, 1996). In the present study, we revisit the ion kinetic modification of multi-dimensional TDNLS through more rigorous derivations, which is consistent with the past kinetic modification of DNLS. Although the original TDNLS was derived in the multi-dimensional form, the evolution of waves with finite propagation angles in TDNLS has not been paid much attention. Applicability of the resultant models to solar wind turbulence is discussed.

Modeling the internal combustion engine

NASA Technical Reports Server (NTRS)

Zeleznik, F. J.; Mcbride, B. J.

1985-01-01

A flexible and computationally economical model of the internal combustion engine was developed for use on large digital computer systems. It is based on a system of ordinary differential equations for cylinder-averaged properties. The computer program is capable of multicycle calculations, with some parameters varying from cycle to cycle, and has restart capabilities. It can accommodate a broad spectrum of reactants, permits changes in physical properties, and offers a wide selection of alternative modeling functions without any reprogramming. It readily adapts to the amount of information available in a particular case because the model is in fact a hierarchy of five models. The models range from a simple model requiring only thermodynamic properties to a complex model demanding full combustion kinetics, transport properties, and poppet valve flow characteristics. Among its many features the model includes heat transfer, valve timing, supercharging, motoring, finite burning rates, cycle-to-cycle variations in air-fuel ratio, humid air, residual and recirculated exhaust gas, and full combustion kinetics.

Solid phase extraction of copper(II) by fixed bed procedure on cation exchange complexing resins.

PubMed

Pesavento, Maria; Sturini, Michela; D'Agostino, Girolamo; Biesuz, Raffaela

2010-02-19

The efficiency of the metal ion recovery by solid phase extraction (SPE) in complexing resins columns is predicted by a simple model based on two parameters reflecting the sorption equilibria and kinetics of the metal ion on the considered resin. The parameter related to the adsorption equilibria was evaluated by the Gibbs-Donnan model, and that related to the kinetics by assuming that the ion exchange is the adsorption rate determining step. The predicted parameters make it possible to evaluate the breakthrough volume of the considered metal ion, Cu(II), from different kinds of complexing resins, and at different conditions, such as acidity and ionic composition. Copyright 2009. Published by Elsevier B.V.

The Biotin/Avidin complex adhesion force

NASA Astrophysics Data System (ADS)

Balsera, Manel A.; Izrailev, Sergei; Stepaniants, Sergey; Oono, Yoshitsugu; Schulten, Klaus

1997-03-01

The vitamin Biotin and the protein avidin form one of the strongest non-covalent bonds between biological molecules. We have performed molecular and stochastic dynamic modeling of the unbinding of this complex(Izrailev et al., Biophysical Journal, In press). These simulations provide insight into the effect of particular residues and water on the tight binding of the system. With the aid of simple phenomenological models we have related qualitatively our results to Atomic Force Microscopy adhesion force measurements (E.-L. Florin, V. T. Moy and H. E. Gaub Science) 264:415-417 and kinetic dissociation experiments( A. Chilcotti and P. S. Stayton, J. Am. Chem. Soc.) 117:10622-10628. We will discuss the difficulties preventing a more quantitative understanding of the unbinding force and kinetics.

Estimating kinetic mechanisms with prior knowledge I: Linear parameter constraints.

PubMed

Salari, Autoosa; Navarro, Marco A; Milescu, Mirela; Milescu, Lorin S

2018-02-05

To understand how ion channels and other proteins function at the molecular and cellular levels, one must decrypt their kinetic mechanisms. Sophisticated algorithms have been developed that can be used to extract kinetic parameters from a variety of experimental data types. However, formulating models that not only explain new data, but are also consistent with existing knowledge, remains a challenge. Here, we present a two-part study describing a mathematical and computational formalism that can be used to enforce prior knowledge into the model using constraints. In this first part, we focus on constraints that enforce explicit linear relationships involving rate constants or other model parameters. We develop a simple, linear algebra-based transformation that can be applied to enforce many types of model properties and assumptions, such as microscopic reversibility, allosteric gating, and equality and inequality parameter relationships. This transformation converts the set of linearly interdependent model parameters into a reduced set of independent parameters, which can be passed to an automated search engine for model optimization. In the companion article, we introduce a complementary method that can be used to enforce arbitrary parameter relationships and any constraints that quantify the behavior of the model under certain conditions. The procedures described in this study can, in principle, be coupled to any of the existing methods for solving molecular kinetics for ion channels or other proteins. These concepts can be used not only to enforce existing knowledge but also to formulate and test new hypotheses. © 2018 Salari et al.

Removal of copper by oxygenated pyrolytic tire char: kinetics and mechanistic insights.

PubMed

Quek, Augustine; Balasubramanian, Rajashekhar

2011-04-01

The kinetics of copper ion (Cu(II)) removal from aqueous solution by pyrolytic tire char was modeled using five different conventional models. A modification to these models was also developed through a modified equation that accounts for precipitation. Conventional first- and second-order reaction models did not fit the copper sorption kinetics well, indicating a lack of simple rate-order dependency on solute concentration. Instead, a reversible first-order rate reaction showed the best fit to the data, indicating a dependence on surface functional groups. Due to the varying solution pH during the sorption process, modified external and internal mass transfer models were employed. Results showed that the sorption of copper onto oxygenated chars was limited by external mass transfer and internal resistance with and without the modification. However, the modification of the sorption process produced very different results for unoxygenated chars, which showed neither internal nor external limitation to sorption. Instead, its slow sorption rate indicates a lack of surface functional groups. The sorption of Cu(II) by oxygenated and unoxygenated chars was also found to occur via three and two distinct stages, respectively. Copyright © 2010 Elsevier Inc. All rights reserved.

The effect of time-dependent macromolecular crowding on the kinetics of protein aggregation: a simple model for the onset of age-related neurodegenerative disease

NASA Astrophysics Data System (ADS)

Minton, Allen

2014-08-01

A linear increase in the concentration of "inert" macromolecules with time is incorporated into simple excluded volume models for protein condensation or fibrillation. Such models predict a long latent period during which no significant amount of protein aggregates, followed by a steep increase in the total amount of aggregate. The elapsed time at which these models predict half-conversion of model protein to aggregate varies by less than a factor of two when the intrinsic rate constant for condensation or fibril growth of the protein is varied over many orders of magnitude. It is suggested that this concept can explain why the symptoms of neurodegenerative diseases associated with the aggregation of very different proteins and peptides appear at approximately the same advanced age in humans.

Cellular level models as tools for cytokine design.

PubMed

Radhakrishnan, Mala L; Tidor, Bruce

2010-01-01

Cytokines and growth factors are critical regulators that connect intracellular and extracellular environments through binding to specific cell-surface receptors. They regulate a wide variety of immunological, growth, and inflammatory response processes. The overall signal initiated by a population of cytokine molecules over long time periods is controlled by the subtle interplay of binding, signaling, and trafficking kinetics. Building on the work of others, we abstract a simple kinetic model that captures relevant features from cytokine systems as well as related growth factor systems. We explore a large range of potential biochemical behaviors, through systematic examination of the model's parameter space. Different rates for the same reaction topology lead to a dramatic range of biochemical network properties and outcomes. Evolution might productively explore varied and different portions of parameter space to create beneficial behaviors, and effective human therapeutic intervention might be achieved through altering network kinetic properties. Quantitative analysis of the results reveals the basis for tensions among a number of different network characteristics. For example, strong binding of cytokine to receptor can increase short-term receptor activation and signal initiation but decrease long-term signaling due to internalization and degradation. Further analysis reveals the role of specific biochemical processes in modulating such tensions. For instance, the kinetics of cytokine binding and receptor activation modulate whether ligand-receptor dissociation can generally occur before signal initiation or receptor internalization. Beyond analysis, the same models and model behaviors provide an important basis for the design of more potent cytokine therapeutics by providing insight into how binding kinetics affect ligand potency. (c) 2010 American Institute of Chemical Engineers

Unfolding and melting of DNA (RNA) hairpins: the concept of structure-specific 2D dynamic landscapes.

PubMed

Lin, Milo M; Meinhold, Lars; Shorokhov, Dmitry; Zewail, Ahmed H

2008-08-07

A 2D free-energy landscape model is presented to describe the (un)folding transition of DNA/RNA hairpins, together with molecular dynamics simulations and experimental findings. The dependence of the (un)folding transition on the stem sequence and the loop length is shown in the enthalpic and entropic contributions to the free energy. Intermediate structures are well defined by the two coordinates of the landscape during (un)zipping. Both the free-energy landscape model and the extensive molecular dynamics simulations totaling over 10 mus predict the existence of temperature-dependent kinetic intermediate states during hairpin (un)zipping and provide the theoretical description of recent ultrafast temperature-jump studies which indicate that hairpin (un)zipping is, in general, not a two-state process. The model allows for lucid prediction of the collapsed state(s) in simple 2D space and we term it the kinetic intermediate structure (KIS) model.

General mechanism of two-state protein folding kinetics.

PubMed

Rollins, Geoffrey C; Dill, Ken A

2014-08-13

We describe here a general model of the kinetic mechanism of protein folding. In the Foldon Funnel Model, proteins fold in units of secondary structures, which form sequentially along the folding pathway, stabilized by tertiary interactions. The model predicts that the free energy landscape has a volcano shape, rather than a simple funnel, that folding is two-state (single-exponential) when secondary structures are intrinsically unstable, and that each structure along the folding path is a transition state for the previous structure. It shows how sequential pathways are consistent with multiple stochastic routes on funnel landscapes, and it gives good agreement with the 9 order of magnitude dependence of folding rates on protein size for a set of 93 proteins, at the same time it is consistent with the near independence of folding equilibrium constant on size. This model gives estimates of folding rates of proteomes, leading to a median folding time in Escherichia coli of about 5 s.

Optical Measurements in Non-Equilibrium Plasmas and Flows

DTIC Science & Technology

2009-09-01

collision model, the exponent x is equal to 0.5, from simple kinetic theory. For most realistic inter-molecular potentials, the exponent x is in the range...Chemical Physics, Vol. 89, p. 5568 (1988). 9. Rosasco, G.J., Lempert, W., Hurst , W.S., and Fein, A., in “Spectral Line Shapes, Vol 2, Walter de Gruyter

Simple stochastic simulation.

PubMed

Schilstra, Maria J; Martin, Stephen R

2009-01-01

Stochastic simulations may be used to describe changes with time of a reaction system in a way that explicitly accounts for the fact that molecules show a significant degree of randomness in their dynamic behavior. The stochastic approach is almost invariably used when small numbers of molecules or molecular assemblies are involved because this randomness leads to significant deviations from the predictions of the conventional deterministic (or continuous) approach to the simulation of biochemical kinetics. Advances in computational methods over the three decades that have elapsed since the publication of Daniel Gillespie's seminal paper in 1977 (J. Phys. Chem. 81, 2340-2361) have allowed researchers to produce highly sophisticated models of complex biological systems. However, these models are frequently highly specific for the particular application and their description often involves mathematical treatments inaccessible to the nonspecialist. For anyone completely new to the field to apply such techniques in their own work might seem at first sight to be a rather intimidating prospect. However, the fundamental principles underlying the approach are in essence rather simple, and the aim of this article is to provide an entry point to the field for a newcomer. It focuses mainly on these general principles, both kinetic and computational, which tend to be not particularly well covered in specialist literature, and shows that interesting information may even be obtained using very simple operations in a conventional spreadsheet.

Growth morphologies of wax in the presence of kinetic inhibitors

NASA Astrophysics Data System (ADS)

Tetervak, Alexander A.

Driven by the need to prevent crystallization of normal alkanes from diesel fuels in cold climates, the petroleum industry has developed additives to slow the growth of these crystals and alter their morphologies. Although the utility of these kinetic inhibitors has been well demonstrated in the field, few studies have directly monitored their effect at microscopic morphology, and the mechanisms by which they act remain poorly understood. Here we present a study of the effects of such additives on the crystallization of long-chain n-alkanes from solution. The additives change the growth morphology from plate-like crystals to a microcrystalline mesh. When we impose a front velocity by moving the sample through a temperature gradient, the mesh growth may form a macroscopic banded pattern and also exhibit a burst-crystallization behavior. In this study, we characterize these crystallization phenomena and also two growth models: a continuum model that demonstrates the essential behavior of the banded crystallization, and a simple qualitative cellular automata model that captures basics of the burst-crystallization process. Keywords: solidification; mesh crystallization; kinetic inhibitor; burst growth.

Understanding system dynamics of an adaptive enzyme network from globally profiled kinetic parameters.

PubMed

Chiang, Austin W T; Liu, Wei-Chung; Charusanti, Pep; Hwang, Ming-Jing

2014-01-15

A major challenge in mathematical modeling of biological systems is to determine how model parameters contribute to systems dynamics. As biological processes are often complex in nature, it is desirable to address this issue using a systematic approach. Here, we propose a simple methodology that first performs an enrichment test to find patterns in the values of globally profiled kinetic parameters with which a model can produce the required system dynamics; this is then followed by a statistical test to elucidate the association between individual parameters and different parts of the system's dynamics. We demonstrate our methodology on a prototype biological system of perfect adaptation dynamics, namely the chemotaxis model for Escherichia coli. Our results agreed well with those derived from experimental data and theoretical studies in the literature. Using this model system, we showed that there are motifs in kinetic parameters and that these motifs are governed by constraints of the specified system dynamics. A systematic approach based on enrichment statistical tests has been developed to elucidate the relationships between model parameters and the roles they play in affecting system dynamics of a prototype biological network. The proposed approach is generally applicable and therefore can find wide use in systems biology modeling research.

Mass accommodation of water: bridging the gap between molecular dynamics simulations and kinetic condensation models.

PubMed

Julin, Jan; Shiraiwa, Manabu; Miles, Rachael E H; Reid, Jonathan P; Pöschl, Ulrich; Riipinen, Ilona

2013-01-17

The condensational growth of submicrometer aerosol particles to climate relevant sizes is sensitive to their ability to accommodate vapor molecules, which is described by the mass accommodation coefficient. However, the underlying processes are not yet fully understood. We have simulated the mass accommodation and evaporation processes of water using molecular dynamics, and the results are compared to the condensation equations derived from the kinetic gas theory to shed light on the compatibility of the two. Molecular dynamics simulations were performed for a planar TIP4P-Ew water surface at four temperatures in the range 268-300 K as well as two droplets, with radii of 1.92 and 4.14 nm at T = 273.15 K. The evaporation flux from molecular dynamics was found to be in good qualitative agreement with that predicted by the simple kinetic condensation equations. Water droplet growth was also modeled with the kinetic multilayer model KM-GAP of Shiraiwa et al. [Atmos. Chem. Phys. 2012, 12, 2777]. It was found that, due to the fast transport across the interface, the growth of a pure water droplet is controlled by gas phase diffusion. These facts indicate that the simple kinetic treatment is sufficient in describing pure water condensation and evaporation. The droplet size was found to have minimal effect on the value of the mass accommodation coefficient. The mass accommodation coefficient was found to be unity (within 0.004) for all studied surfaces, which is in agreement with previous simulation work. Additionally, the simulated evaporation fluxes imply that the evaporation coefficient is also unity. Comparing the evaporation rates of the mass accommodation and evaporation simulations indicated that the high collision flux, corresponding to high supersaturation, present in typical molecular dynamics mass accommodation simulations can under certain conditions lead to an increase in the evaporation rate. Consequently, in such situations the mass accommodation coefficient can be overestimated, but in the present cases the corrected values were still close to unity with the lowest value at ≈0.99.

Mass Accommodation of Water: Bridging the Gap Between Molecular Dynamics Simulations and Kinetic Condensation Models

PubMed Central

2012-01-01

The condensational growth of submicrometer aerosol particles to climate relevant sizes is sensitive to their ability to accommodate vapor molecules, which is described by the mass accommodation coefficient. However, the underlying processes are not yet fully understood. We have simulated the mass accommodation and evaporation processes of water using molecular dynamics, and the results are compared to the condensation equations derived from the kinetic gas theory to shed light on the compatibility of the two. Molecular dynamics simulations were performed for a planar TIP4P-Ew water surface at four temperatures in the range 268–300 K as well as two droplets, with radii of 1.92 and 4.14 nm at T = 273.15 K. The evaporation flux from molecular dynamics was found to be in good qualitative agreement with that predicted by the simple kinetic condensation equations. Water droplet growth was also modeled with the kinetic multilayer model KM-GAP of Shiraiwa et al. [Atmos. Chem. Phys.2012, 117, 2777]. It was found that, due to the fast transport across the interface, the growth of a pure water droplet is controlled by gas phase diffusion. These facts indicate that the simple kinetic treatment is sufficient in describing pure water condensation and evaporation. The droplet size was found to have minimal effect on the value of the mass accommodation coefficient. The mass accommodation coefficient was found to be unity (within 0.004) for all studied surfaces, which is in agreement with previous simulation work. Additionally, the simulated evaporation fluxes imply that the evaporation coefficient is also unity. Comparing the evaporation rates of the mass accommodation and evaporation simulations indicated that the high collision flux, corresponding to high supersaturation, present in typical molecular dynamics mass accommodation simulations can under certain conditions lead to an increase in the evaporation rate. Consequently, in such situations the mass accommodation coefficient can be overestimated, but in the present cases the corrected values were still close to unity with the lowest value at ≈0.99. PMID:23253100

Predicting protein decomposition: the case of aspartic-acid racemization kinetics.

PubMed Central

Collins, M J; Waite, E R; van Duin, A C

1999-01-01

The increase in proportion of the non-biological (D-) isomer of aspartic acid (Asp) relative to the L-isomer has been widely used in archaeology and geochemistry as a tool for dating. the method has proved controversial, particularly when used for bones. The non-linear kinetics of Asp racemization have prompted a number of suggestions as to the underlying mechanism(s) and have led to the use of mathematical transformations which linearize the increase in D-Asp with respect to time. Using one example, a suggestion that the initial rapid phase of Asp racemization is due to a contribution from asparagine (Asn), we demonstrate how a simple model of the degradation and racemization of Asn can be used to predict the observed kinetics. A more complex model of peptide bound Asx (Asn + Asp) racemization, which occurs via the formation of a cyclic succinimide (Asu), can be used to correctly predict Asx racemization kinetics in proteins at high temperatures (95-140 degrees C). The model fails to predict racemization kinetics in dentine collagen at 37 degrees C. The reason for this is that Asu formation is highly conformation dependent and is predicted to occur extremely slowly in triple helical collagen. As conformation strongly influences the rate of Asu formation and hence Asx racemization, the use of extrapolation from high temperatures to estimate racemization kinetics of Asx in proteins below their denaturation temperature is called into question. In the case of archaeological bone, we argue that the D:L ratio of Asx reflects the proportion of non-helical to helical collagen, overlain by the effects of leaching of more soluble (and conformationally unconstrained) peptides. Thus, racemization kinetics in bone are potentially unpredictable, and the proposed use of Asx racemization to estimate the extent of DNA depurination in archaeological bones is challenged. PMID:10091247

2D Kinetic Particle in Cell Simulations of a Shear-Flow Stabilized Z-Pinch

NASA Astrophysics Data System (ADS)

Tummel, Kurt; Higginson, Drew; Schmidt, Andrea; Link, Anthony; McLean, Harry; Shumlak, Uri; Nelson, Brian; Golingo, Raymond; Claveau, Elliot; Lawrence Livermore National Lab Team; University of Washington Team

2016-10-01

The Z-pinch is a relatively simple and attractive potential fusion reactor design, but attempts to develop such a reactor have consistently struggled to overcome Z-pinch instabilities. The ``sausage'' and ``kink'' modes are among the most robust and prevalent Z-pinch instabilities, but theory and simulations suggest that axial flow-shear, dvz / dr ≠ 0 , can suppress these modes. Experiments have confirmed that Z-pinch plasmas with embedded axial flow-shear display a significantly enhanced resilience to the sausage and kink modes at a demonstration current of 50kAmps. A new experiment is under way to test the concept at higher current, and efforts to model these plasmas are being expanded. The performance and stability of these devices will depend on features like the plasma viscosity, anomalous resistivity, and finite Larmor radius effects, which are most accurately characterized in kinetic models. To predict these features, kinetic simulations using the particle in cell code LSP are now in development, and initial benchmarking and 2D stability analyses of the sausage mode are presented here. These results represent the first kinetic modeling of the flow-shear stabilized Z-pinch. This work is funded by the USDOE/ARPAe Alpha Program. Prepared by LLNL under Contract DE-AC52-07NA27344.

Combination of Markov state models and kinetic networks for the analysis of molecular dynamics simulations of peptide folding.

PubMed

Radford, Isolde H; Fersht, Alan R; Settanni, Giovanni

2011-06-09

Atomistic molecular dynamics simulations of the TZ1 beta-hairpin peptide have been carried out using an implicit model for the solvent. The trajectories have been analyzed using a Markov state model defined on the projections along two significant observables and a kinetic network approach. The Markov state model allowed for an unbiased identification of the metastable states of the system, and provided the basis for commitment probability calculations performed on the kinetic network. The kinetic network analysis served to extract the main transition state for folding of the peptide and to validate the results from the Markov state analysis. The combination of the two techniques allowed for a consistent and concise characterization of the dynamics of the peptide. The slowest relaxation process identified is the exchange between variably folded and denatured species, and the second slowest process is the exchange between two different subsets of the denatured state which could not be otherwise identified by simple inspection of the projected trajectory. The third slowest process is the exchange between a fully native and a partially folded intermediate state characterized by a native turn with a proximal backbone H-bond, and frayed side-chain packing and termini. The transition state for the main folding reaction is similar to the intermediate state, although a more native like side-chain packing is observed.

The apparent solubility of aluminum (III) in Hanford high-level waste.

PubMed

Reynolds, Jacob G

2012-01-01

The solubility of aluminum in Hanford nuclear waste impacts on the processability of the waste by a number of proposed treatment options. For many years, Hanford staff has anecdotally noted that aluminum appears to be considerably more soluble in Hanford waste than the simpler electrolyte solutions used as analogues. There has been minimal scientific study to confirm these anecdotal observations, however. The present study determines the apparent solubility product for gibbsite in 50 tank samples. The ratio of hydroxide to aluminum in the liquid phase for the samples is calculated and plotted as a function of total sodium molarity. Total sodium molarity is used as a surrogate for ionic strength, because the relative ratios of mono-, di- and trivalent anions are not available for all of the samples. These results were compared to the simple NaOH-NaAl(OH)(4)-H(2)O system, and the NaOH-NaAl(OH)(4)-NaCl-H(2)O system data retrieved from the literature. The results show that gibbsite is apparently more soluble in the samples than in the simple systems whenever the sodium molarity is greater than 2M. This apparent enhanced solubility cannot be explained solely by differences in ionic strength. The change in solubility with ionic strength in simple systems is small compared to the difference between aluminum solubility in Hanford waste and the simple systems. The reason for the apparent enhanced solubility is unknown, but could include kinetic or thermodynamic factors that are not present in the simple electrolyte systems. Any kinetic explanation would have to explain why the samples are always supersaturated whenever the sodium molarity is above 2M. Real waste characterization data should not be used to validate thermodynamic solubility models until it can be confirmed that the apparent enhanced gibbsite solubility is a thermodynamic effect and not a kinetic effect.

The apparent solubility of aluminum (III) in Hanford high-level waste

DOE Office of Scientific and Technical Information (OSTI.GOV)

Reynolds, Jacob G.

2012-12-01

The solubility of aluminum in Hanford nuclear waste impacts on the process ability of the waste by a number of proposed treatment options. For many years, Hanford staff has anecdotally noted that aluminum appears to be considerably more soluble in Hanford waste than the simpler electrolyte solutions used as analogues. There has been minimal scientific study to confirm these anecdotal observations, however. The present study determines the apparent solubility product for gibbsite in 50 tank samples. The ratio of hydroxide to aluminum in the liquid phase for the samples is calculated and plotted as a function of total sodium molarity.more » Total sodium molarity is used as a surrogate for ionic strength, because the relative ratios of mono, di and trivalent anions are not available for all of the samples. These results were compared to the simple NaOH-NaAl(OH{sub 4})H{sub 2}O system, and the NaOH-NaAl(OH{sub 4})NaCl-H{sub 2}O system data retrieved from the literature. The results show that gibbsite is apparently more soluble in the samples than in the simple systems whenever the sodium molarity is greater than two. This apparent enhanced solubility cannot be explained solely by differences in ionic strength. The change in solubility with ionic strength in simple systems is small compared to the difference between aluminum solubility in Hanford waste and the simple systems. The reason for the apparent enhanced solubility is unknown, but could include. kinetic or thermodynamic factors that are not present in the simple electrolyte systems. Any kinetic explanation would have to explain why the samples are always supersaturated whenever the sodium molarity is above two. Real waste characterization data should not be used to validate thermodynamic solubility models until it can be confirmed that the apparent enhanced gibbsite solubility is a thermodynamic effect and not a kinetic effect.« less

Thermodynamics of information processing based on enzyme kinetics: An exactly solvable model of an information pump

NASA Astrophysics Data System (ADS)

Cao, Yuansheng; Gong, Zongping; Quan, H. T.

2015-06-01

Motivated by the recent proposed models of the information engine [Proc. Natl. Acad. Sci. USA 109, 11641 (2012), 10.1073/pnas.1204263109] and the information refrigerator [Phys. Rev. Lett. 111, 030602 (2013), 10.1103/PhysRevLett.111.030602], we propose a minimal model of the information pump and the information eraser based on enzyme kinetics. This device can either pump molecules against the chemical potential gradient by consuming the information to be encoded in the bit stream or (partially) erase the information initially encoded in the bit stream by consuming the Gibbs free energy. The dynamics of this model is solved exactly, and the "phase diagram" of the operation regimes is determined. The efficiency and the power of the information machine is analyzed. The validity of the second law of thermodynamics within our model is clarified. Our model offers a simple paradigm for the investigating of the thermodynamics of information processing involving the chemical potential in small systems.

Modeling the Free Carrier Recombination Kinetics in PTB7:PCBM Organic Photovoltaics

DOE PAGES

Oosterhout, Stefan D.; Ferguson, Andrew J.; Larson, Bryon W.; ...

2016-10-03

Currently the exact recombination mechanism of free carriers in organic photovoltaic (OPV) devices is poorly understood. Often a reduced Langevin model is used to describe the decay behavior of electrons and holes. Here we propose a novel, simple kinetic model that accurately describes the decay behavior of free carriers in the PTB7:PCBM organic photovoltaic blend. In order to accurately describe the recombination behavior of free carriers as measured by time-resolved microwave conductivity (TRMC), this model needs to only take into account free and trapped holes in the polymer, and free electrons in the fullerene. The model is consistent for differentmore » PTB7:PCBM blend ratios and spans a light intensity range of over 3 orders of magnitude. Furthermore, the model demonstrates that dark carriers exist in the polymer and interact with photoinduced charge carriers, and that the trapping and detrapping rates of the holes are of high importance to the overall carrier lifetime.« less

A Simple Laboratory Experiment to Determine the Kinetics of Mutarotation of D-Glucose Using a Blood Glucose Meter

ERIC Educational Resources Information Center

Perles, Carlos E.; Volpe, Pedro L. O.

2008-01-01

A simple commercial blood glucose meter is used to follow the kinetics of mutarotation of D-glucose in aqueous solution. The results may be compared with those obtained using an automatic polarimeter, if this is available This experiment is proposed for use by students in a general chemistry, biology, organic chemistry, and physical chemistry…

The Distribution of Chromosomal Aberrations in Human Cells Predicted by a Generalized Time-Dependent Model of Radiation-Induced Formation of Aberrations

NASA Technical Reports Server (NTRS)

Ponomarev, Artem L.; George, K.; Cucinotta, F. A.

2011-01-01

New experimental data show how chromosomal aberrations for low- and high-LET radiation are dependent on DSB repair deficiencies in wild-type, AT and NBS cells. We simulated the development of chromosomal aberrations in these cells lines in a stochastic track-structure-dependent model, in which different cells have different kinetics of DSB repair. We updated a previously formulated model of chromosomal aberrations, which was based on a stochastic Monte Carlo approach, to consider the time-dependence of DSB rejoining. The previous version of the model had an assumption that all DSBs would rejoin, and therefore we called it a time-independent model. The chromosomal-aberrations model takes into account the DNA and track structure for low- and high-LET radiations, and provides an explanation and prediction of the statistics of rare and more complex aberrations. We compared the program-simulated kinetics of DSB rejoining to the experimentally-derived bimodal exponential curves of the DSB kinetics. We scored the formation of translocations, dicentrics, acentric and centric rings, deletions, and inversions. The fraction of DSBs participating in aberrations was studied in relation to the rejoining time. Comparisons of simulated dose dependence for simple aberrations to the experimental dose-dependence for HF19, AT and NBS cells will be made.

How Kinetics within the Unfolded State Affects Protein Folding: an Analysis Based on Markov State Models and an Ultra-Long MD Trajectory

PubMed Central

Deng, Nan-jie; Dai, Wei

2013-01-01

Understanding how kinetics in the unfolded state affects protein folding is a fundamentally important yet less well-understood issue. Here we employ three different models to analyze the unfolded landscape and folding kinetics of the miniprotein Trp-cage. The first is a 208 μs explicit solvent molecular dynamics (MD) simulation from D. E. Shaw Research containing tens of folding events. The second is a Markov state model (MSM-MD) constructed from the same ultra-long MD simulation; MSM-MD can be used to generate thousands of folding events. The third is a Markov state model built from temperature replica exchange MD simulations in implicit solvent (MSM-REMD). All the models exhibit multiple folding pathways, and there is a good correspondence between the folding pathways from direct MD and those computed from the MSMs. The unfolded populations interconvert rapidly between extended and collapsed conformations on time scales ≤ 40 ns, compared with the folding time of ≈ 5 μs. The folding rates are independent of where the folding is initiated from within the unfolded ensemble. About 90 % of the unfolded states are sampled within the first 40 μs of the ultra-long MD trajectory, which on average explores ~27 % of the unfolded state ensemble between consecutive folding events. We clustered the folding pathways according to structural similarity into “tubes”, and kinetically partitioned the unfolded state into populations that fold along different tubes. From our analysis of the simulations and a simple kinetic model, we find that when the mixing within the unfolded state is comparable to or faster than folding, the folding waiting times for all the folding tubes are similar and the folding kinetics is essentially single exponential despite the presence of heterogeneous folding paths with non-uniform barriers. When the mixing is much slower than folding, different unfolded populations fold independently leading to non-exponential kinetics. A kinetic partition of the Trp-cage unfolded state is constructed which reveals that different unfolded populations have almost the same probability to fold along any of the multiple folding paths. We are investigating whether the results for the kinetics in the unfolded state of the twenty-residue Trp-cage is representative of larger single domain proteins. PMID:23705683

Thermodynamic Analysis of Chemically Reacting Mixtures-Comparison of First and Second Order Models.

PubMed

Pekař, Miloslav

2018-01-01

Recently, a method based on non-equilibrium continuum thermodynamics which derives thermodynamically consistent reaction rate models together with thermodynamic constraints on their parameters was analyzed using a triangular reaction scheme. The scheme was kinetically of the first order. Here, the analysis is further developed for several first and second order schemes to gain a deeper insight into the thermodynamic consistency of rate equations and relationships between chemical thermodynamic and kinetics. It is shown that the thermodynamic constraints on the so-called proper rate coefficient are usually simple sign restrictions consistent with the supposed reaction directions. Constraints on the so-called coupling rate coefficients are more complex and weaker. This means more freedom in kinetic coupling between reaction steps in a scheme, i.e., in the kinetic effects of other reactions on the rate of some reaction in a reacting system. When compared with traditional mass-action rate equations, the method allows a reduction in the number of traditional rate constants to be evaluated from data, i.e., a reduction in the dimensionality of the parameter estimation problem. This is due to identifying relationships between mass-action rate constants (relationships which also include thermodynamic equilibrium constants) which have so far been unknown.

Simple, Fast and Accurate Implementation of the Diffusion Approximation Algorithm for Stochastic Ion Channels with Multiple States

PubMed Central

Orio, Patricio; Soudry, Daniel

2012-01-01

Background The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled gating particles, while the DA was modeled using uncoupled gating particles. Implementations of DA with coupled particles, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. Main Contributions We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable – allowing an easy, transparent and efficient DA implementation, avoiding unnecessary approximations. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods, except when short time steps or low channel numbers were used. PMID:22629320

Modelling of individual subject ozone exposure response kinetics.

PubMed

Schelegle, Edward S; Adams, William C; Walby, William F; Marion, M Susan

2012-06-01

A better understanding of individual subject ozone (O(3)) exposure response kinetics will provide insight into how to improve models used in the risk assessment of ambient ozone exposure. To develop a simple two compartment exposure-response model that describes individual subject decrements in forced expiratory volume in one second (FEV(1)) induced by the acute inhalation of O(3) lasting up to 8 h. FEV(1) measurements of 220 subjects who participated in 14 previously completed studies were fit to the model using both particle swarm and nonlinear least squares optimization techniques to identify three subject-specific coefficients producing minimum "global" and local errors, respectively. Observed and predicted decrements in FEV(1) of the 220 subjects were used for validation of the model. Further validation was provided by comparing the observed O(3)-induced FEV(1) decrements in an additional eight studies with predicted values obtained using model coefficients estimated from the 220 subjects used in cross validation. Overall the individual subject measured and modeled FEV(1) decrements were highly correlated (mean R(2) of 0.69 ± 0.24). In addition, it was shown that a matrix of individual subject model coefficients can be used to predict the mean and variance of group decrements in FEV(1). This modeling approach provides insight into individual subject O(3) exposure response kinetics and provides a potential starting point for improving the risk assessment of environmental O(3) exposure.

Film growth kinetics and electric field patterning during electrospray deposition of block copolymer thin films

NASA Astrophysics Data System (ADS)

Toth, Kristof; Hu, Hanqiong; Choo, Youngwoo; Loewenberg, Michael; Osuji, Chinedum

The delivery of sub-micron droplets of dilute polymer solutions to a heated substrate by electrospray deposition (ESD) enables precisely controlled and continuous growth of block copolymer (BCP) thin films. Here we explore patterned deposition of BCP films by spatially varying the electric field at the substrate using an underlying charged grid, as well as film growth kinetics. Numerical analysis was performed to examine pattern fidelity by considering the trajectories of charged droplets during flight through imposed periodic field variations in the vicinity of the substrate. Our work uncovered an unexpected modality for improving the resolution of the patterning process via stronger field focusing through the use of a second oppositely charged grid beneath a primary focusing array, with an increase in highly localized droplet deposition on the intersecting nodes of the grid. Substrate coverage kinetics are considered for homopolymer deposition in the context of simple kinetic models incorporating temperature and molecular weight dependence of diffusivity. By contrast, film coverage kinetics for block copolymer depositions are additionally convoluted with preferential wetting and thickness-periodicity commensurability effects. NSF GRFP.

A new methodology to determine kinetic parameters for one- and two-step chemical models

NASA Technical Reports Server (NTRS)

Mantel, T.; Egolfopoulos, F. N.; Bowman, C. T.

1996-01-01

In this paper, a new methodology to determine kinetic parameters for simple chemical models and simple transport properties classically used in DNS of premixed combustion is presented. First, a one-dimensional code is utilized to performed steady unstrained laminar methane-air flame in order to verify intrinsic features of laminar flames such as burning velocity and temperature and concentration profiles. Second, the flame response to steady and unsteady strain in the opposed jet configuration is numerically investigated. It appears that for a well determined set of parameters, one- and two-step mechanisms reproduce the extinction limit of a laminar flame submitted to a steady strain. Computations with the GRI-mech mechanism (177 reactions, 39 species) and multicomponent transport properties are used to validate these simplified models. A sensitivity analysis of the preferential diffusion of heat and reactants when the Lewis number is close to unity indicates that the response of the flame to an oscillating strain is very sensitive to this number. As an application of this methodology, the interaction between a two-dimensional vortex pair and a premixed laminar flame is performed by Direct Numerical Simulation (DNS) using the one- and two-step mechanisms. Comparison with the experimental results of Samaniego et al. (1994) shows a significant improvement in the description of the interaction when the two-step model is used.

Simulations of 4D edge transport and dynamics using the TEMPEST gyro-kinetic code

NASA Astrophysics Data System (ADS)

Rognlien, T. D.; Cohen, B. I.; Cohen, R. H.; Dorr, M. R.; Hittinger, J. A. F.; Kerbel, G. D.; Nevins, W. M.; Xiong, Z.; Xu, X. Q.

2006-10-01

Simulation results are presented for tokamak edge plasmas with a focus on the 4D (2r,2v) option of the TEMPEST continuum gyro-kinetic code. A detailed description of a variety of kinetic simulations is reported, including neoclassical radial transport from Coulomb collisions, electric field generation, dynamic response to perturbations by geodesic acoustic modes, and parallel transport on open magnetic-field lines. Comparison is made between the characteristics of the plasma solutions on closed and open magnetic-field line regions separated by a magnetic separatrix, and simple physical models are used to qualitatively explain the differences observed in mean flow and electric-field generation. The status of extending the simulations to 5D turbulence will be summarized. The code structure used in this ongoing project is also briefly described, together with future plans.

Ultrasensitive dual phosphorylation dephosphorylation cycle kinetics exhibits canonical competition behavior

NASA Astrophysics Data System (ADS)

Huang, Qingdao; Qian, Hong

2009-09-01

We establish a mathematical model for a cellular biochemical signaling module in terms of a planar differential equation system. The signaling process is carried out by two phosphorylation-dephosphorylation reaction steps that share common kinase and phosphatase with saturated enzyme kinetics. The pair of equations is particularly simple in the present mathematical formulation, but they are singular. A complete mathematical analysis is developed based on an elementary perturbation theory. The dynamics exhibits the canonical competition behavior in addition to bistability. Although widely understood in ecological context, we are not aware of a full range of biochemical competition in a simple signaling network. The competition dynamics has broad implications to cellular processes such as cell differentiation and cancer immunoediting. The concepts of homogeneous and heterogeneous multisite phosphorylation are introduced and their corresponding dynamics are compared: there is no bistability in a heterogeneous dual phosphorylation system. A stochastic interpretation is also provided that further gives intuitive understanding of the bistable behavior inside the cells.

Analysis of Mathematical Modelling on Potentiometric Biosensors

PubMed Central

Mehala, N.; Rajendran, L.

2014-01-01

A mathematical model of potentiometric enzyme electrodes for a nonsteady condition has been developed. The model is based on the system of two coupled nonlinear time-dependent reaction diffusion equations for Michaelis-Menten formalism that describes the concentrations of substrate and product within the enzymatic layer. Analytical expressions for the concentration of substrate and product and the corresponding flux response have been derived for all values of parameters using the new homotopy perturbation method. Furthermore, the complex inversion formula is employed in this work to solve the boundary value problem. The analytical solutions obtained allow a full description of the response curves for only two kinetic parameters (unsaturation/saturation parameter and reaction/diffusion parameter). Theoretical descriptions are given for the two limiting cases (zero and first order kinetics) and relatively simple approaches for general cases are presented. All the analytical results are compared with simulation results using Scilab/Matlab program. The numerical results agree with the appropriate theories. PMID:25969765

Analysis of mathematical modelling on potentiometric biosensors.

PubMed

Mehala, N; Rajendran, L

2014-01-01

A mathematical model of potentiometric enzyme electrodes for a nonsteady condition has been developed. The model is based on the system of two coupled nonlinear time-dependent reaction diffusion equations for Michaelis-Menten formalism that describes the concentrations of substrate and product within the enzymatic layer. Analytical expressions for the concentration of substrate and product and the corresponding flux response have been derived for all values of parameters using the new homotopy perturbation method. Furthermore, the complex inversion formula is employed in this work to solve the boundary value problem. The analytical solutions obtained allow a full description of the response curves for only two kinetic parameters (unsaturation/saturation parameter and reaction/diffusion parameter). Theoretical descriptions are given for the two limiting cases (zero and first order kinetics) and relatively simple approaches for general cases are presented. All the analytical results are compared with simulation results using Scilab/Matlab program. The numerical results agree with the appropriate theories.

Statistical validity of using ratio variables in human kinetics research.

PubMed

Liu, Yuanlong; Schutz, Robert W

2003-09-01

The purposes of this study were to investigate the validity of the simple ratio and three alternative deflation models and examine how the variation of the numerator and denominator variables affects the reliability of a ratio variable. A simple ratio and three alternative deflation models were fitted to four empirical data sets, and common criteria were applied to determine the best model for deflation. Intraclass correlation was used to examine the component effect on the reliability of a ratio variable. The results indicate that the validity, of a deflation model depends on the statistical characteristics of the particular component variables used, and an optimal deflation model for all ratio variables may not exist. Therefore, it is recommended that different models be fitted to each empirical data set to determine the best deflation model. It was found that the reliability of a simple ratio is affected by the coefficients of variation and the within- and between-trial correlations between the numerator and denominator variables. It was recommended that researchers should compute the reliability of the derived ratio scores and not assume that strong reliabilities in the numerator and denominator measures automatically lead to high reliability in the ratio measures.

Experimental determination of gap flow-conditioned forces at turbine stages and their effect on the running stability of simple rotors. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Wohlrab, R.

1983-01-01

Instabilities in turbine operation can be caused by forces which are produced in connection with motions involving the oil film in the bearings. An experimental investigation regarding the characteristics of such forces in the case of three typical steam turbine stages is conducted, taking into account the effect of various parameters. Supplementary kinetic tests are carried out to obtain an estimate of the flow forces which are proportional to the velocity. The measurements are based on the theoretical study of the damping characteristics of a vibrational model. A computational analysis of the effect of the measured fluid forces on the stability characteristics of simple rotor model is also conducted.

Charge transfer kinetics at the solid-solid interface in porous electrodes

NASA Astrophysics Data System (ADS)

Bai, Peng; Bazant, Martin Z.

2014-04-01

Interfacial charge transfer is widely assumed to obey the Butler-Volmer kinetics. For certain liquid-solid interfaces, the Marcus-Hush-Chidsey theory is more accurate and predictive, but it has not been applied to porous electrodes. Here we report a simple method to extract the charge transfer rates in carbon-coated LiFePO4 porous electrodes from chronoamperometry experiments, obtaining curved Tafel plots that contradict the Butler-Volmer equation but fit the Marcus-Hush-Chidsey prediction over a range of temperatures. The fitted reorganization energy matches the Born solvation energy for electron transfer from carbon to the iron redox site. The kinetics are thus limited by electron transfer at the solid-solid (carbon-LixFePO4) interface rather than by ion transfer at the liquid-solid interface, as previously assumed. The proposed experimental method generalizes Chidsey’s method for phase-transforming particles and porous electrodes, and the results show the need to incorporate Marcus kinetics in modelling batteries and other electrochemical systems.

Effect of nonlinearity in hybrid kinetic Monte Carlo-continuum models.

PubMed

Balter, Ariel; Lin, Guang; Tartakovsky, Alexandre M

2012-01-01

Recently there has been interest in developing efficient ways to model heterogeneous surface reactions with hybrid computational models that couple a kinetic Monte Carlo (KMC) model for a surface to a finite-difference model for bulk diffusion in a continuous domain. We consider two representative problems that validate a hybrid method and show that this method captures the combined effects of nonlinearity and stochasticity. We first validate a simple deposition-dissolution model with a linear rate showing that the KMC-continuum hybrid agrees with both a fully deterministic model and its analytical solution. We then study a deposition-dissolution model including competitive adsorption, which leads to a nonlinear rate, and show that in this case the KMC-continuum hybrid and fully deterministic simulations do not agree. However, we are able to identify the difference as a natural result of the stochasticity coming from the KMC surface process. Because KMC captures inherent fluctuations, we consider it to be more realistic than a purely deterministic model. Therefore, we consider the KMC-continuum hybrid to be more representative of a real system.

Net growth rate of continuum heterogeneous biofilms with inhibition kinetics.

PubMed

Gonzo, Elio Emilio; Wuertz, Stefan; Rajal, Veronica B

2018-01-01

Biofilm systems can be modeled using a variety of analytical and numerical approaches, usually by making simplifying assumptions regarding biofilm heterogeneity and activity as well as effective diffusivity. Inhibition kinetics, albeit common in experimental systems, are rarely considered and analytical approaches are either lacking or consider effective diffusivity of the substrate and the biofilm density to remain constant. To address this obvious knowledge gap an analytical procedure to estimate the effectiveness factor (dimensionless substrate mass flux at the biofilm-fluid interface) was developed for a continuum heterogeneous biofilm with multiple limiting-substrate Monod kinetics to different types of inhibition kinetics. The simple perturbation technique, previously validated to quantify biofilm activity, was applied to systems where either the substrate or the inhibitor is the limiting component, and cases where the inhibitor is a reaction product or the substrate also acts as the inhibitor. Explicit analytical equations are presented for the effectiveness factor estimation and, therefore, the calculation of biomass growth rate or limiting substrate/inhibitor consumption rate, for a given biofilm thickness. The robustness of the new biofilm model was tested using kinetic parameters experimentally determined for the growth of Pseudomonas putida CCRC 14365 on phenol. Several additional cases have been analyzed, including examples where the effectiveness factor can reach values greater than unity, characteristic of systems with inhibition kinetics. Criteria to establish when the effectiveness factor can reach values greater than unity in each of the cases studied are also presented.

Phenomenology of wall-bounded Newtonian turbulence.

PubMed

L'vov, Victor S; Pomyalov, Anna; Procaccia, Itamar; Zilitinkevich, Sergej S

2006-01-01

We construct a simple analytic model for wall-bounded turbulence, containing only four adjustable parameters. Two of these parameters are responsible for the viscous dissipation of the components of the Reynolds stress tensor. The other two parameters control the nonlinear relaxation of these objects. The model offers an analytic description of the profiles of the mean velocity and the correlation functions of velocity fluctuations in the entire boundary region, from the viscous sublayer, through the buffer layer, and further into the log-law turbulent region. In particular, the model predicts a very simple distribution of the turbulent kinetic energy in the log-law region between the velocity components: the streamwise component contains a half of the total energy whereas the wall-normal and cross-stream components contain a quarter each. In addition, the model predicts a very simple relation between the von Kármán slope k and the turbulent velocity in the log-law region v+ (in wall units): v+=6k. These predictions are in excellent agreement with direct numerical simulation data and with recent laboratory experiments.

Kinetic Damage from Meteorites

NASA Technical Reports Server (NTRS)

Cooke, William; Brown, Peter; Matney, Mark

2017-01-01

A Near Earth object impacting into Earth's atmosphere may produce damaging effects at the surface due to airblast, thermal pulse, or kinetic impact in the form of meteorites. At large sizes (>many tens of meters), the damage is amplified by the hypersonic impact of these large projectiles moving with cosmic velocity, leaving explosively produced craters. However, much more common is simple "kinetic" damage caused by the impact of smaller meteorites moving at terminal speeds. As of this date a handful of instances are definitively known of people or structures being directly hit and/or damaged by the kinetic impact of meteorites. Meteorites known to have struck humans include the Sylacauga, Alabama fall (1954) and the Mbale meteorite fall (1992). Much more common is kinetic meteorite damage to cars, buildings, and even a post box (Claxton, Georgia - 1984). Historical accounts indicate that direct kinetic damage by meteorites may be more common than recent accounts suggest (Yau et al., 1994). In this talk we will examine the contemporary meteorite flux and estimate the frequency of kinetic damage to various structures, as well as how the meteorite flux might affect the rate of human casualties. This will update an earlier study by Halliday et al (1985), adding variations expected in meteorite flux with latitude (Le Feuvre and Wieczorek, 2008) and validating these model predictions of speed and entry angle with observations from the NASA and SOMN fireball networks. In particular, we explore the physical characteristics of bright meteors which may be used as a diagnostic for estimating which fireballs produce meteorites and hence how early warning of such kinetic damage may be estimated in advance through observations and modelling.

Kinetic Damage from Meteorites

NASA Technical Reports Server (NTRS)

Cooke, William; Brown, Peter; Matney, Mark

2017-01-01

A Near Earth object impacting into Earth's atmosphere may produce damaging effects at the surface due to airblast, thermal pulse, or kinetic impact in the form of meteorites. At large sizes (greater than many tens of meters), the damage is amplified by the hypersonic impact of these large projectiles moving with cosmic velocity, leaving explosively produced craters. However, much more common is simple "kinetic" damage caused by the impact of smaller meteorites moving at terminal speeds. As of this date a handful of instances are definitively known of people or structures being directly hit and/or damaged by the kinetic impact of meteorites. Meteorites known to have struck humans include the Sylacauga, Alabama fall (1954) and the Mbale meteorite fall (1992). Much more common is kinetic meteorite damage to cars, buildings, and even a post box (Claxton, Georgia - 1984). Historical accounts indicate that direct kinetic damage by meteorites may be more common than recent accounts suggest (Yau et al., 1994). In this talk we will examine the contemporary meteorite flux and estimate the frequency of kinetic damage to various structures, as well as how the meteorite flux might affect the rate of human casualties. This will update an earlier study by Halliday et al (1985), adding variations expected in meteorite flux with latitude (Le Feuvre and Wieczorek, 2008) and validating these model predictions of speed and entry angle with observations from the NASA and SOMN fireball networks. In particular, we explore the physical characteristics of bright meteors which may be used as a diagnostic for estimating which fireballs produce meteorites and hence how early warning of such kinetic damage may be estimated in advance through observations and modeling.

The Study of a Simple Redox Reaction as an Experimental Approach to Chemical Kinetics.

ERIC Educational Resources Information Center

Elias, Horst; Zipp, Arden P.

1988-01-01

Recommends using iodide ions and peroxodisulfate ions for studying rate laws instead of the standard iodine clock for kinetic study. Presents the methodology and a discussion of the kinetics involved for a laboratory experiment for a high school or introductory college course. (ML)

Hidden Quantum Processes, Quantum Ion Channels, and 1/ fθ-Type Noise.

PubMed

Paris, Alan; Vosoughi, Azadeh; Berman, Stephen A; Atia, George

2018-07-01

In this letter, we perform a complete and in-depth analysis of Lorentzian noises, such as those arising from [Formula: see text] and [Formula: see text] channel kinetics, in order to identify the source of [Formula: see text]-type noise in neurological membranes. We prove that the autocovariance of Lorentzian noise depends solely on the eigenvalues (time constants) of the kinetic matrix but that the Lorentzian weighting coefficients depend entirely on the eigenvectors of this matrix. We then show that there are rotations of the kinetic eigenvectors that send any initial weights to any target weights without altering the time constants. In particular, we show there are target weights for which the resulting Lorenztian noise has an approximately [Formula: see text]-type spectrum. We justify these kinetic rotations by introducing a quantum mechanical formulation of membrane stochastics, called hidden quantum activated-measurement models, and prove that these quantum models are probabilistically indistinguishable from the classical hidden Markov models typically used for ion channel stochastics. The quantum dividend obtained by replacing classical with quantum membranes is that rotations of the Lorentzian weights become simple readjustments of the quantum state without any change to the laboratory-determined kinetic and conductance parameters. Moreover, the quantum formalism allows us to model the activation energy of a membrane, and we show that maximizing entropy under constrained activation energy yields the previous [Formula: see text]-type Lorentzian weights, in which the spectral exponent [Formula: see text] is a Lagrange multiplier for the energy constraint. Thus, we provide a plausible neurophysical mechanism by which channel and membrane kinetics can give rise to [Formula: see text]-type noise (something that has been occasionally denied in the literature), as well as a realistic and experimentally testable explanation for the numerical values of the spectral exponents. We also discuss applications of quantum membranes beyond [Formula: see text]-type -noise, including applications to animal models and possible impact on quantum foundations.

A study of the effect of space-dependent neutronics on stochastically-induced bifurcations in BWR dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Analytis, G.T.

1995-09-01

A non-linear one-group space-dependent neutronic model for a finite one-dimensional core is coupled with a simple BWR feed-back model. In agreement with results obtained by the authors who originally developed the point-kinetics version of this model, we shall show numerically that stochastic reactivity excitations may result in limit-cycles and eventually in a chaotic behaviour, depending on the magnitude of the feed-back coefficient K. In the framework of this simple space-dependent model, the effect of the non-linearities on the different spatial harmonics is studied and the importance of the space-dependent effects is exemplified and assessed in terms of the importance ofmore » the higher harmonics. It is shown that under certain conditions, when the limit-cycle-type develop, the neutron spectra may exhibit strong space-dependent effects.« less

Study of drug concentration effects on in vitro lipolysis kinetics in medium-chain triglycerides by considering oil viscosity and surface tension.

PubMed

Arnold, Yvonne Elisabeth; Imanidis, Georgios; Kuentz, Martin

2011-10-09

Simple oil formulations are widely used in oral drug delivery and the fate of these systems is governed mainly by the dispersion and digestion process. The current work aimed to study concentration effects of six poorly water-soluble drugs on the in vitro lipolysis rate of medium-chain triglycerides. The results were compared with drug effects on oil viscosity and surface tension. First the different drugs were characterized by molecular modeling and their influence on physical oil properties was assessed. Herein capillary viscosimetry was employed as well as dynamic surface tensiometry. Subsequently, an apparent in vitro lipolysis rate was determined in biorelevant medium using an automated pH stat titrator linked to a thermo-controlled vessel. The different drugs exhibited varying effects on oil viscosity and surface tension. However, all drugs significantly lowered the apparent lipolysis rate of the oil. This effect was very similar among the different compounds with exception of orlistat, which practically blocked lipolysis because of a potent direct inhibition. The other drugs affected lipolysis kinetics most likely by different mechanism(s). In light of the obtained results, a drug effect on oil viscosity or surface tension appeared to play a minor role in reducing the lipolysis rate. The lipolysis kinetics was further not affected by the drug load, which was deemed advantageous from a pharmaceutical viewpoint. Different dose strengths are therefore not assumed to alter lipolysis kinetics, which is beneficial for limiting the variability of in vivo drug release. Further studies of drug solubility kinetics in the evolving digestion phases are, however, needed to finally assess potential effects of dosage strength in simple oil formulations. Copyright © 2011. Published by Elsevier B.V.

Polymer Fluid Dynamics: Continuum and Molecular Approaches.

PubMed

Bird, R B; Giacomin, A J

2016-06-07

To solve problems in polymer fluid dynamics, one needs the equations of continuity, motion, and energy. The last two equations contain the stress tensor and the heat-flux vector for the material. There are two ways to formulate the stress tensor: (a) One can write a continuum expression for the stress tensor in terms of kinematic tensors, or (b) one can select a molecular model that represents the polymer molecule and then develop an expression for the stress tensor from kinetic theory. The advantage of the kinetic theory approach is that one gets information about the relation between the molecular structure of the polymers and the rheological properties. We restrict the discussion primarily to the simplest stress tensor expressions or constitutive equations containing from two to four adjustable parameters, although we do indicate how these formulations may be extended to give more complicated expressions. We also explore how these simplest expressions are recovered as special cases of a more general framework, the Oldroyd 8-constant model. Studying the simplest models allows us to discover which types of empiricisms or molecular models seem to be worth investigating further. We also explore equivalences between continuum and molecular approaches. We restrict the discussion to several types of simple flows, such as shearing flows and extensional flows, which are of greatest importance in industrial operations. Furthermore, if these simple flows cannot be well described by continuum or molecular models, then it is not necessary to lavish time and energy to apply them to more complex flow problems.

High resolution approach to the native state ensemble kinetics and thermodynamics.

PubMed

Wu, Sangwook; Zhuravlev, Pavel I; Papoian, Garegin A

2008-12-15

Many biologically interesting functions such as allosteric switching or protein-ligand binding are determined by the kinetics and mechanisms of transitions between various conformational substates of the native basin of globular proteins. To advance our understanding of these processes, we constructed a two-dimensional free energy surface (FES) of the native basin of a small globular protein, Trp-cage. The corresponding order parameters were defined using two native substructures of Trp-cage. These calculations were based on extensive explicit water all-atom molecular dynamics simulations. Using the obtained two-dimensional FES, we studied the transition kinetics between two Trp-cage conformations, finding that switching process shows a borderline behavior between diffusive and weakly-activated dynamics. The transition is well-characterized kinetically as a biexponential process. We also introduced a new one-dimensional reaction coordinate for the conformational transition, finding reasonable qualitative agreement with the two-dimensional kinetics results. We investigated the distribution of all the 38 native nuclear magnetic resonance structures on the obtained FES, analyzing interactions that stabilize specific low-energy conformations. Finally, we constructed a FES for the same system but with simple dielectric model of water instead of explicit water, finding that the results were surprisingly similar in a small region centered on the native conformations. The dissimilarities between the explicit and implicit model on the larger-scale point to the important role of water in mediating interactions between amino acid residues.

Parameter estimation in tree graph metabolic networks.

PubMed

Astola, Laura; Stigter, Hans; Gomez Roldan, Maria Victoria; van Eeuwijk, Fred; Hall, Robert D; Groenenboom, Marian; Molenaar, Jaap J

2016-01-01

We study the glycosylation processes that convert initially toxic substrates to nutritionally valuable metabolites in the flavonoid biosynthesis pathway of tomato (Solanum lycopersicum) seedlings. To estimate the reaction rates we use ordinary differential equations (ODEs) to model the enzyme kinetics. A popular choice is to use a system of linear ODEs with constant kinetic rates or to use Michaelis-Menten kinetics. In reality, the catalytic rates, which are affected among other factors by kinetic constants and enzyme concentrations, are changing in time and with the approaches just mentioned, this phenomenon cannot be described. Another problem is that, in general these kinetic coefficients are not always identifiable. A third problem is that, it is not precisely known which enzymes are catalyzing the observed glycosylation processes. With several hundred potential gene candidates, experimental validation using purified target proteins is expensive and time consuming. We aim at reducing this task via mathematical modeling to allow for the pre-selection of most potential gene candidates. In this article we discuss a fast and relatively simple approach to estimate time varying kinetic rates, with three favorable properties: firstly, it allows for identifiable estimation of time dependent parameters in networks with a tree-like structure. Secondly, it is relatively fast compared to usually applied methods that estimate the model derivatives together with the network parameters. Thirdly, by combining the metabolite concentration data with a corresponding microarray data, it can help in detecting the genes related to the enzymatic processes. By comparing the estimated time dynamics of the catalytic rates with time series gene expression data we may assess potential candidate genes behind enzymatic reactions. As an example, we show how to apply this method to select prominent glycosyltransferase genes in tomato seedlings.

Nuclear-coupled thermal-hydraulic stability analysis of boiling water reactors

NASA Astrophysics Data System (ADS)

Karve, Atul A.

We have studied the nuclear-coupled thermal-hydraulic stability of boiling water reactors (BWRs) using a model we developed from: the space-time modal neutron kinetics equations based on spatial omega-modes, the equations for two-phase flow in parallel boiling channels, the fuel rod heat conduction equations, and a simple model for the recirculation loop. The model is represented as a dynamical system comprised of time-dependent nonlinear ordinary differential equations, and it is studied using stability analysis, modern bifurcation theory, and numerical simulations. We first determine the stability boundary (SB) in the most relevant parameter plane, the inlet-subcooling-number/external-pressure-drop plane, for a fixed control rod induced external reactivity equal to the 100% rod line value and then transform the SB to the practical power-flow map. Using this SB, we show that the normal operating point at 100% power is very stable, stability of points on the 100% rod line decreases as the flow rate is reduced, and that points are least stable in the low-flow/high-power region. We also determine the SB when the modal kinetics is replaced by simple point reactor kinetics and show that the first harmonic mode has no significant effect on the SB. Later we carry out the relevant numerical simulations where we first show that the Hopf bifurcation, that occurs as a parameter is varied across the SB is subcritical, and that, in the important low-flow/high-power region, growing oscillations can result following small finite perturbations of stable steady-states on the 100% rod line. Hence, a point on the 100% rod line in the low-flow/high-power region, although stable, may nevertheless be a point at which a BWR should not be operated. Numerical simulations are then done to calculate the decay ratios (DRs) and frequencies of oscillations for various points on the 100% rod line. It is determined that the NRC requirement of DR < 0.75-0.8 is not rigorously satisfied in the low-flow/high-power region and hence these points should be avoided during normal startup and shutdown operations. The frequency of oscillation is shown to decrease as the flow rate is reduced and the frequency of 0.5Hz observed in the low-flow/high-power region is consistent with those observed during actual instability incidents. Additional numerical simulations show that in the low-flow/high-power region, for the same initial conditions, the use of point kinetics leads to damped oscillations, whereas the model that includes the modal kinetics equations results in growing nonlinear oscillations. Thus, we show that side-by-side out-of-phase growing power oscillations result due to the very important first harmonic mode effect and that the use of point kinetics, which fails to predict these growing oscillations, leads to dramatically nonconservative results. Finally, the effect of a simple recirculation loop model that we develop is studied by carrying out additional stability analyses and additional numerical simulations. It is shown that the loop has a stabilizing effect on certain points on the 100% rod line for time delays equal to integer multiples of the natural period of oscillation, whereas it has a destabilizing effect for half-integer multiples. However, for more practical time delays, it is determined that the overall effect generally is destabilizing.

Historical perspective on lead biokinetic models.

PubMed Central

Rabinowitz, M

1998-01-01

A historical review of the development of biokinetic model of lead is presented. Biokinetics is interpreted narrowly to mean only physiologic processes happening within the body. Proceeding chronologically, for each epoch, the measurements of lead in the body are presented along with mathematical models in an attempt to trace the convergence of observations from two disparate fields--occupational medicine and radiologic health--into some unified models. Kehoe's early balance studies and the use of radioactive lead tracers are presented. The 1960s saw the joint application of radioactive lead techniques and simple compartmental kinetic models used to establish the exchange rates and residence times of lead in body pools. The applications of stable isotopes to questions of the magnitudes of respired and ingested inputs required the development of a simple three-pool model. During the 1980s more elaborate models were developed. One of their key goals was the establishment of the dose-response relationship between exposure to lead and biologic precursors of adverse health effects. PMID:9860905

A nanojet: propulsion of a molecular machine by an asymmetric distribution of reaction--products

NASA Astrophysics Data System (ADS)

Liverpool, Tanniemola; Golestanian, Ramin; Ajdari, Armand

2006-03-01

A simple model for the reaction-driven propulsion of a small device is proposed as a model for (part of) a molecular machine in aqueous media. Motion of the device is driven by an asymmetric distribution of reaction products. We calculate the propulsive velocity of the device as well as the scale of the velocity fluctuations. We also consider the effects of hydrodynamic flow as well as a number of different scenarios for the kinetics of the reaction.

Propulsion of a Molecular Machine by Asymmetric Distribution of Reaction Products

NASA Astrophysics Data System (ADS)

Golestanian, Ramin; Liverpool, Tanniemola B.; Ajdari, Armand

2005-06-01

A simple model for the reaction-driven propulsion of a small device is proposed as a model for (part of) a molecular machine in aqueous media. The motion of the device is driven by an asymmetric distribution of reaction products. The propulsive velocity of the device is calculated as well as the scale of the velocity fluctuations. The effects of hydrodynamic flow as well as a number of different scenarios for the kinetics of the reaction are addressed.

Propulsion of a molecular machine by asymmetric distribution of reaction products.

PubMed

Golestanian, Ramin; Liverpool, Tanniemola B; Ajdari, Armand

2005-06-10

A simple model for the reaction-driven propulsion of a small device is proposed as a model for (part of) a molecular machine in aqueous media. The motion of the device is driven by an asymmetric distribution of reaction products. The propulsive velocity of the device is calculated as well as the scale of the velocity fluctuations. The effects of hydrodynamic flow as well as a number of different scenarios for the kinetics of the reaction are addressed.

Current sheet in plasma as a system with a controlling parameter

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fridman, Yu. A., E-mail: yulya-fridman@yandex.ru; Chukbar, K. V., E-mail: Chukbar-KV@nrcki.ru

2015-08-15

A simple kinetic model describing stationary solutions with bifurcated and single-peaked current density profiles of a plane electron beam or current sheet in plasma is presented. A connection is established between the two-dimensional constructions arising in terms of the model and the one-dimensional considerations by Bernstein−Greene−Kruskal facilitating the reconstruction of the distribution function of trapped particles when both the profile of the electric potential and the free particles distribution function are known.

Refined elasticity sampling for Monte Carlo-based identification of stabilizing network patterns.

PubMed

Childs, Dorothee; Grimbs, Sergio; Selbig, Joachim

2015-06-15

Structural kinetic modelling (SKM) is a framework to analyse whether a metabolic steady state remains stable under perturbation, without requiring detailed knowledge about individual rate equations. It provides a representation of the system's Jacobian matrix that depends solely on the network structure, steady state measurements, and the elasticities at the steady state. For a measured steady state, stability criteria can be derived by generating a large number of SKMs with randomly sampled elasticities and evaluating the resulting Jacobian matrices. The elasticity space can be analysed statistically in order to detect network positions that contribute significantly to the perturbation response. Here, we extend this approach by examining the kinetic feasibility of the elasticity combinations created during Monte Carlo sampling. Using a set of small example systems, we show that the majority of sampled SKMs would yield negative kinetic parameters if they were translated back into kinetic models. To overcome this problem, a simple criterion is formulated that mitigates such infeasible models. After evaluating the small example pathways, the methodology was used to study two steady states of the neuronal TCA cycle and the intrinsic mechanisms responsible for their stability or instability. The findings of the statistical elasticity analysis confirm that several elasticities are jointly coordinated to control stability and that the main source for potential instabilities are mutations in the enzyme alpha-ketoglutarate dehydrogenase. © The Author 2015. Published by Oxford University Press.

Shear viscosity for dense plasmas by equilibrium molecular dynamics in asymmetric Yukawa ionic mixtures.

PubMed

Haxhimali, Tomorr; Rudd, Robert E; Cabot, William H; Graziani, Frank R

2015-11-01

We present molecular dynamics (MD) calculations of shear viscosity for asymmetric mixed plasma for thermodynamic conditions relevant to astrophysical and inertial confinement fusion plasmas. Specifically, we consider mixtures of deuterium and argon at temperatures of 100-500 eV and a number density of 10^{25} ions/cc. The motion of 30,000-120,000 ions is simulated in which the ions interact via the Yukawa (screened Coulomb) potential. The electric field of the electrons is included in this effective interaction; the electrons are not simulated explicitly. Shear viscosity is calculated using the Green-Kubo approach with an integral of the shear stress autocorrelation function, a quantity calculated in the equilibrium MD simulations. We systematically study different mixtures through a series of simulations with increasing fraction of the minority high-Z element (Ar) in the D-Ar plasma mixture. In the more weakly coupled plasmas, at 500 eV and low Ar fractions, results from MD compare very well with Chapman-Enskog kinetic results. In the more strongly coupled plasmas, the kinetic theory does not agree well with the MD results. We develop a simple model that interpolates between classical kinetic theories at weak coupling and the Murillo Yukawa viscosity model at higher coupling. This hybrid kinetics-MD viscosity model agrees well with the MD results over the conditions simulated, ranging from moderately weakly coupled to moderately strongly coupled asymmetric plasma mixtures.

Shear viscosity for dense plasmas by equilibrium molecular dynamics in asymmetric Yukawa ionic mixtures

NASA Astrophysics Data System (ADS)

Haxhimali, Tomorr; Rudd, Robert E.; Cabot, William H.; Graziani, Frank R.

2015-11-01

We present molecular dynamics (MD) calculations of shear viscosity for asymmetric mixed plasma for thermodynamic conditions relevant to astrophysical and inertial confinement fusion plasmas. Specifically, we consider mixtures of deuterium and argon at temperatures of 100-500 eV and a number density of 1025 ions/cc. The motion of 30 000-120 000 ions is simulated in which the ions interact via the Yukawa (screened Coulomb) potential. The electric field of the electrons is included in this effective interaction; the electrons are not simulated explicitly. Shear viscosity is calculated using the Green-Kubo approach with an integral of the shear stress autocorrelation function, a quantity calculated in the equilibrium MD simulations. We systematically study different mixtures through a series of simulations with increasing fraction of the minority high-Z element (Ar) in the D-Ar plasma mixture. In the more weakly coupled plasmas, at 500 eV and low Ar fractions, results from MD compare very well with Chapman-Enskog kinetic results. In the more strongly coupled plasmas, the kinetic theory does not agree well with the MD results. We develop a simple model that interpolates between classical kinetic theories at weak coupling and the Murillo Yukawa viscosity model at higher coupling. This hybrid kinetics-MD viscosity model agrees well with the MD results over the conditions simulated, ranging from moderately weakly coupled to moderately strongly coupled asymmetric plasma mixtures.

Genetic and molecular dosimetry of HZE radiation (7-IML-1)

NASA Technical Reports Server (NTRS)

Nelson, Gregory A.

1992-01-01

The objectives of the study are to determine the kinetics of production and to characterize the unique aspects of genetic and developmental lesion induced in animal cells by radiation present in the space environment. Special attention is given to heavy charged particles. The organism Caenorhabditis elegans, a simple nematode, is used as a model system for a coordinated set of ground-based and flight experiments.

Lumped mass formulations for modeling flexible body systems

NASA Technical Reports Server (NTRS)

Rampalli, Rajiv

1989-01-01

The efforts of Mechanical Dynamics, Inc. in obtaining a general formulation for flexible bodies in a multibody setting are discussed. The efforts being supported by MDI, both in house and externally are summarized. The feasibility of using lumped mass approaches to modeling flexibility in a multibody dynamics context is examined. The kinematics and kinetics for a simple system consisting of two rigid bodies connected together by an elastic beam are developed in detail. Accuracy, efficiency and ease of use using this approach are some of the issues that are then looked at. The formulation is then generalized to a superelement containing several nodes and connecting several bodies. Superelement kinematics and kinetics equations are developed. The feasibility and effectiveness of the method is illustrated by the use of some examples illustrating phenomena common in the context of spacecraft motions.

Coherent Rayleigh-Brillouin scattering measurements of bulk viscosity of polar and nonpolar gases, and kinetic theory.

PubMed

Meijer, A S; de Wijn, A S; Peters, M F E; Dam, N J; van de Water, W

2010-10-28

We investigate coherent Rayleigh-Brillouin spectroscopy as an efficient process to measure the bulk viscosity of gases at gigahertz frequencies. Scattered spectral distributions are measured using a Fizeau spectrometer. We discuss the statistical error due to the fluctuating mode structure of the used pump laser. Experiments were done for both polar and nonpolar gases and the bulk viscosity was obtained from the spectra using the Tenti S6 model. Results are compared to simple classical kinetic models of molecules with internal degrees of freedom. At the extremely high (gigahertz) frequencies of our experiment, most internal vibrational modes remain frozen and the bulk viscosity is dominated by the rotational degrees of freedom. Our measurements show that the molecular dipole moments have unexpectedly little influence on the bulk viscosity at room temperature and moderate pressure.

A Simple Kinetic Model for the Growth of Fe2B Layers on AISI 1026 Steel During the Powder-pack Boriding

NASA Astrophysics Data System (ADS)

Flores-Rentería, M. A.; Ortiz-Domínguez, M.; Keddam, M.; Damián-Mejía, O.; Elias-Espinosa, M.; Flores-González, M. A.; Medina-Moreno, S. A.; Cruz-Avilés, A.; Villanueva-Ibañez, M.

2015-02-01

This work focused on the determination of boron diffusion coefficient through the Fe2B layers on AISI 1026 steel using a mathematical model. The suggested model solves the mass balance equation at the (Fe2B/substrate) interface. This thermochemical treatment was carried out in the temperature range of 1123-1273 K for a treatment time ranging from 2 to 8 h. The generated boride layers were characterized by different experimental techniques such as light optical microscopy, scanning electron microscopy, XRD analysis and the Daimler-Benz Rockwell-C indentation technique. As a result, the boron activation energy for AISI 1026 steel was estimated as 178.4 kJ/mol. Furthermore, this kinetic model was validated by comparing the experimental Fe2B layer thickness with the predicted one at a temperature of 1253 K for 5 h of treatment. A contour diagram relating the layer thickness to the boriding parameters was proposed to be used in practical applications.

General Mechanism of Two-State Protein Folding Kinetics

PubMed Central

Rollins, Geoffrey C.; Dill, Ken A.

2016-01-01

We describe here a general model of the kinetic mechanism of protein folding. In the Foldon Funnel Model, proteins fold in units of secondary structures, which form sequentially along the folding pathway, stabilized by tertiary interactions. The model predicts that the free energy landscape has a volcano shape, rather than a simple funnel, that folding is two-state (single-exponential) when secondary structures are intrinsically unstable, and that each structure along the folding path is a transition state for the previous structure. It shows how sequential pathways are consistent with multiple stochastic routes on funnel landscapes, and it gives good agreement with the 9 order of magnitude dependence of folding rates on protein size for a set of 93 proteins, at the same time it is consistent with the near independence of folding equilibrium constant on size. This model gives estimates of folding rates of proteomes, leading to a median folding time in Escherichia coli of about 5 s. PMID:25056406

Kinetic model for dependence of thin film stress on growth rate, temperature, and microstructure

NASA Astrophysics Data System (ADS)

Chason, E.; Shin, J. W.; Hearne, S. J.; Freund, L. B.

2012-04-01

During deposition, many thin films go through a range of stress states, changing from compressive to tensile and back again. In addition, the stress depends strongly on the processing and material parameters. We have developed a simple analytical model to describe the stress evolution in terms of a kinetic competition between different mechanisms of stress generation and relaxation at the triple junction where the surface and grain boundary intersect. The model describes how the steady state stress scales with the dimensionless parameter D/LR where D is the diffusivity, R is the growth rate, and L is the grain size. It also explains the transition from tensile to compressive stress as the microstructure evolves from isolated islands to a continuous film. We compare calculations from the model with measurements of the stress dependence on grain size and growth rate in the steady state regime and of the evolution of stress with thickness for different temperatures.

QUENCHING OF CARBON MONOXIDE AND METHANE IN THE ATMOSPHERES OF COOL BROWN DWARFS AND HOT JUPITERS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Visscher, Channon; Moses, Julianne I., E-mail: visscher@lpi.usra.edu, E-mail: jmoses@spacescience.org

We explore CO{r_reversible}CH{sub 4} quench kinetics in the atmospheres of substellar objects using updated timescale arguments, as suggested by a thermochemical kinetics and diffusion model that transitions from the thermochemical-equilibrium regime in the deep atmosphere to a quench-chemical regime at higher altitudes. More specifically, we examine CO quench chemistry on the T dwarf Gliese 229B and CH{sub 4} quench chemistry on the hot-Jupiter HD 189733b. We describe a method for correctly calculating reverse rate coefficients for chemical reactions, discuss the predominant pathways for CO{r_reversible}CH{sub 4} interconversion as indicated by the model, and demonstrate that a simple timescale approach can bemore » used to accurately describe the behavior of quenched species when updated reaction kinetics and mixing-length-scale assumptions are used. Proper treatment of quench kinetics has important implications for estimates of molecular abundances and/or vertical mixing rates in the atmospheres of substellar objects. Our model results indicate significantly higher K{sub zz} values than previously estimated near the CO quench level on Gliese 229B, whereas current-model-data comparisons using CH{sub 4} permit a wide range of K{sub zz} values on HD 189733b. We also use updated reaction kinetics to revise previous estimates of the Jovian water abundance, based upon the observed abundance and chemical behavior of carbon monoxide. The CO chemical/observational constraint, along with Galileo entry probe data, suggests a water abundance of approximately 0.51-2.6 x solar (for a solar value of H{sub 2}O/H{sub 2} = 9.61 x 10{sup -4}) in Jupiter's troposphere, assuming vertical mixing from the deep atmosphere is the only source of tropospheric CO.« less

Direct Observations of a Dynamically Driven Phase Transition with in situ X-Ray Diffraction in a Simple Ionic Crystal

NASA Astrophysics Data System (ADS)

Kalita, Patricia; Specht, Paul; Root, Seth; Sinclair, Nicholas; Schuman, Adam; White, Melanie; Cornelius, Andrew L.; Smith, Jesse; Sinogeikin, Stanislav

2017-12-01

We report real-time observations of a phase transition in the ionic solid CaF2 , a model A B2 structure in high-pressure physics. Synchrotron x-ray diffraction coupled with dynamic loading to 27.7 GPa, and separately with static compression, follows, in situ, the fluorite to cotunnite structural phase transition, both on nanosecond and on minute time scales. Using Rietveld refinement techniques, we examine the kinetics and hysteresis of the transition. Our results give insight into the kinetic time scale of the fluorite-cotunnite phase transition under shock compression, which is relevant to a number of isomorphic compounds.

Kinetics of hairpin ribozyme cleavage in yeast.

PubMed Central

Donahue, C P; Fedor, M J

1997-01-01

Hairpin ribozymes catalyze a self-cleavage reaction that provides a simple model for quantitative analyses of intracellular mechanisms of RNA catalysis. Decay rates of chimeric mRNAs containing self-cleaving ribozymes give a direct measure of intracellular cleavage kinetics in yeast. Intracellular ribozyme-mediated cleavage occurs at similar rates and shows similar inhibition by ribozyme mutations as ribozyme-mediated reactions in vitro, but only when ribozymes are located in a favorable mRNA sequence context. The impact of cleavage on mRNA abundance is shown to depend directly on intrinsic mRNA stability. Surprisingly, cleavage products are no more labile than uncleaved mRNAs despite the loss of terminal cap structures or poly (A). PMID:9292496

Theoretical study of production of unique glasses in space

NASA Technical Reports Server (NTRS)

Larsen, D. C.

1974-01-01

Analytical functional relationships describing homogeneous nucleation and crystallization in various supercooled liquids were developed. The time and temperature dependent relationships of nucleation and crystallization (intrinsic properties) are being used to relate glass forming tendency to extrinsic parameters such as cooling rate through computer simulation. Single oxide systems are being studied initially to aid in developing workable kinetic models and to indicate the primary materials parameters affecting glass formation. The theory and analytical expressions developed for simple systems is then extended to complex oxide systems. A thorough understanding of nucleation and crystallization kinetics of glass forming systems provides a priori knowledge of the ability of a given system to form a glass.

On the origin of why static or breakloose friction is larger than kinetic friction, and how to reduce it: the role of aging, elasticity and sequential interfacial slip.

PubMed

Lorenz, B; Persson, B N J

2012-06-06

We discuss the origin of static friction and show how it can be reduced towards kinetic friction by the appropriate design of the sliding system. The basic idea is to use elastically soft solids and apply the external forces in such a way that different parts of the contacting interface start to slip at different times during the (tangential) loading process. In addition, the local slip must be large enough in order to result in a strong drop in the static friction force. We illustrate the theoretical predictions with the results of a simple model experiment.

Direct Observations of a Dynamically Driven Phase Transition with in situ X-Ray Diffraction in a Simple Ionic Crystal

DOE PAGES

Kalita, Patricia E.; Specht, Paul Elliot; Root, Seth; ...

2017-12-21

Here, we report real-time observations of a phase transition in the ionic solid CaF 2, a model AB 2 structure in high-pressure physics. Synchrotron x-ray diffraction coupled with dynamic loading to 27.7 GPa, and separately with static compression, follows, in situ, the fluorite to cotunnite structural phase transition, both on nanosecond and on minute time scales. Using Rietveld refinement techniques, we examine the kinetics and hysteresis of the transition. Our results give insight into the kinetic time scale of the fluorite-cotunnite phase transition under shock compression, which is relevant to a number of isomorphic compounds.

Describing the dynamics of processes consisting simultaneously of Poissonian and non-Poissonian kinetics

NASA Astrophysics Data System (ADS)

Eule, S.; Friedrich, R.

2013-03-01

Dynamical processes exhibiting non-Poissonian kinetics with nonexponential waiting times are frequently encountered in nature. Examples are biochemical processes like gene transcription which are known to involve multiple intermediate steps. However, often a second process, obeying Poissonian statistics, affects the first one simultaneously, such as the degradation of mRNA in the above example. The aim of the present article is to provide a concise treatment of such random systems which are affected by regular and non-Poissonian kinetics at the same time. We derive the governing master equation and provide a controlled approximation scheme for this equation. The simplest approximation leads to generalized reaction rate equations. For a simple model of gene transcription we solve the resulting equation and show how the time evolution is influenced significantly by the type of waiting time distribution assumed for the non-Poissonian process.

Non-Gaussianity in multi-sound-speed disformally coupled inflation

DOE Office of Scientific and Technical Information (OSTI.GOV)

De Bruck, Carsten van; Longden, Chris; Koivisto, Tomi, E-mail: C.vandeBruck@sheffield.ac.uk, E-mail: tomi.koivisto@nordita.org, E-mail: cjlongden1@sheffield.ac.uk

Most, if not all, scalar-tensor theories are equivalent to General Relativity with a disformally coupled matter sector. In extra-dimensional theories such a coupling can be understood as a result of induction of the metric on a brane that matter is confined to. This article presents a first look at the non-Gaussianities in disformally coupled inflation, a simple two-field model that features a novel kinetic interaction. Cases with both canonical and Dirac-Born-Infeld (DBI) kinetic terms are taken into account, the latter motivated by the possible extra-dimensional origin of the disformality. The computations are carried out for the equilateral configuration in themore » slow-roll regime, wherein it is found that the non-Gaussianity is typically rather small and negative. This is despite the fact that the new kinetic interaction causes the perturbation modes to propagate with different sounds speeds, which may both significantly deviate from unity during inflation.« less

Kinetics and Mechanisms of Chemical and Biological Agents Release from Biopolymeric Microcapsules.

PubMed

Vinceković, Marko; Jurić, Slaven; Đermić, Edyta; Topolovec-Pintarić, Snježana

2017-11-08

Kinetics and mechanisms of copper cations and Trichoderma viride spores release from uncoated and chitosan coated alginate microcapsules were investigated. The gelation of a fixed amount of sodium alginate at different concentrations of copper ion solutions resulted in distinct kinetics and release mechanisms. The increase in copper cation concentration promoted, but the presence of the chitosan layer on the microcapsule surface and the increase in microcapsule size reduced the rate of active agent release. Fitting to simple Korsmeyer-Peppas empirical model revealed that the underlying release mechanism (Fickian diffusion or a combination of the diffusion and erosion mechanisms) depends on the copper cation concentration and presence of T. viride spores. The investigation pointed out that the proper selection of formulation variables helps in designing microcapsules with the desirable release of copper ions and T. viride for plant protection and nutrition.

Kinetics of devolatilization and oxidation of a pulverized biomass in an entrained flow reactor under realistic combustion conditions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jimenez, Santiago; Remacha, Pilar; Ballester, Javier

2008-03-15

In this paper the results of a complete set of devolatilization and combustion experiments performed with pulverized ({proportional_to}500 {mu}m) biomass in an entrained flow reactor under realistic combustion conditions are presented. The data obtained are used to derive the kinetic parameters that best fit the observed behaviors, according to a simple model of particle combustion (one-step devolatilization, apparent oxidation kinetics, thermally thin particles). The model is found to adequately reproduce the experimental trends regarding both volatile release and char oxidation rates for the range of particle sizes and combustion conditions explored. The experimental and numerical procedures, similar to those recentlymore » proposed for the combustion of pulverized coal [J. Ballester, S. Jimenez, Combust. Flame 142 (2005) 210-222], have been designed to derive the parameters required for the analysis of biomass combustion in practical pulverized fuel configurations and allow a reliable characterization of any finely pulverized biomass. Additionally, the results of a limited study on the release rate of nitrogen from the biomass particle along combustion are shown. (author)« less

Discrete kinetic models from funneled energy landscape simulations.

PubMed

Schafer, Nicholas P; Hoffman, Ryan M B; Burger, Anat; Craig, Patricio O; Komives, Elizabeth A; Wolynes, Peter G

2012-01-01

A general method for facilitating the interpretation of computer simulations of protein folding with minimally frustrated energy landscapes is detailed and applied to a designed ankyrin repeat protein (4ANK). In the method, groups of residues are assigned to foldons and these foldons are used to map the conformational space of the protein onto a set of discrete macrobasins. The free energies of the individual macrobasins are then calculated, informing practical kinetic analysis. Two simple assumptions about the universality of the rate for downhill transitions between macrobasins and the natural local connectivity between macrobasins lead to a scheme for predicting overall folding and unfolding rates, generating chevron plots under varying thermodynamic conditions, and inferring dominant kinetic folding pathways. To illustrate the approach, free energies of macrobasins were calculated from biased simulations of a non-additive structure-based model using two structurally motivated foldon definitions at the full and half ankyrin repeat resolutions. The calculated chevrons have features consistent with those measured in stopped flow chemical denaturation experiments. The dominant inferred folding pathway has an "inside-out", nucleation-propagation like character.

Modifier mass transfer kinetic effect in the performance of solvent gradient simulated moving bed (SG-SMB) process

NASA Astrophysics Data System (ADS)

Câmara, L. D. T.

2015-09-01

The solvent-gradient simulated moving bed process (SG-SMB) is the new tendency in the performance improvement if compared to the traditional isocratic solvent conditions. In such SG-SMB separation process the modulation of the solvent strength leads to significant increase in the purities and productivity followed by reduction in the solvent consumption. A stepwise modelling approach was utilized in the representation of the interconnected chromatographic columns of the system combined with lumped mass transfer models between the solid and liquid phase. The influence of the solvent modifier was considered applying the Abel model which takes into account the effect of modifier volume fraction over the partition coefficient. The modelling and simulations were carried out and compared to the experimental SG-SMB separation of the amino acids phenylalanine and tryptophan. A lumped mass transfer kinetic model was applied for both the modifier (ethanol) as well as the solutes. The simulation results showed that such simple and global mass transfer models are enough to represent all the mass transfer effect between the solid adsorbent and the liquid phase. The separation performance can be improved reducing the interaction or the mass transfer kinetic effect between the solid adsorbent phase and the modifier. The simulations showed great agreement fitting the experimental data of the amino acids concentrations both at the extract as well as at the raffinate.

A review of some fish nutrition methodologies.

PubMed

Belal, Ibrahim E H

2005-03-01

Several classical warm blooded animal (poultry, sheep, cows, etc.) methods for dietary nutrients evaluation (digestibility, metabolizablity, and energy budget) are applied to fish, even though fish live in a different environment in addition to being cold blooded animals. These applications have caused significant errors that have made these methods non-additive and meaningless, as is explained in the text. In other words, dietary digestion and absorption could not adequately be measured due to the aquatic environment fish live in. Therefore, net nutrient deposition and/or growth are the only accurate measurement left to evaluate dietary nutrients intake in fish. In order to understand and predict dietary nutrient intake-growth response relationship, several mathematical models; (1) the simple linear equation, (2) the logarithmic equation, and (3) the quadratic equation are generally used. These models however, do not describe a full range of growth and have no biological meaning as explained in the text. On the other hand, a model called the saturation kinetic model. It has biological basis (the law of mass action and the enzyme kinetic) and it describes the full range of growth curve. Additionally, it has four parameters that summarize the growth curve and could also be used in comparing diets or nutrients effect on fish growth and/or net nutrient deposition. The saturation kinetic model is proposed to be adequate for dietary nutrient evaluation for fish. The theoretical derivation of this model is illustrated in the text.

Co-digestion of solid waste: Towards a simple model to predict methane production.

PubMed

Kouas, Mokhles; Torrijos, Michel; Schmitz, Sabine; Sousbie, Philippe; Sayadi, Sami; Harmand, Jérôme

2018-04-01

Modeling methane production is a key issue for solid waste co-digestion. Here, the effect of a step-wise increase in the organic loading rate (OLR) on reactor performance was investigated, and four new models were evaluated to predict methane yields using data acquired in batch mode. Four co-digestion experiments of mixtures of 2 solid substrates were conducted in semi-continuous mode. Experimental methane yields were always higher than the BMP values of mixtures calculated from the BMP of each substrate, highlighting the importance of endogenous production (methane produced from auto-degradation of microbial community and generated solids). The experimental methane productions under increasing OLRs corresponded well to the modeled data using the model with constant endogenous production and kinetics identified at 80% from total batch time. This model provides a simple and useful tool for technical design consultancies and plant operators to optimize the co-digestion and the choice of the OLRs. Copyright © 2018 Elsevier Ltd. All rights reserved.

Semiempirical equations for modeling solid-state kinetics based on a Maxwell-Boltzmann distribution of activation energies: applications to a polymorphic transformation under crystallization slurry conditions and to the thermal decomposition of AgMnO4 crystals.

PubMed

Skrdla, Peter J; Robertson, Rebecca T

2005-06-02

Many solid-state reactions and phase transformations performed under isothermal conditions give rise to asymmetric, sigmoidally shaped conversion-time (x-t) profiles. The mathematical treatment of such curves, as well as their physical interpretation, is often challenging. In this work, the functional form of a Maxwell-Boltzmann (M-B) distribution is used to describe the distribution of activation energies for the reagent solids, which, when coupled with an integrated first-order rate expression, yields a novel semiempirical equation that may offer better success in the modeling of solid-state kinetics. In this approach, the Arrhenius equation is used to relate the distribution of activation energies to a corresponding distribution of rate constants for the individual molecules in the reagent solids. This distribution of molecular rate constants is then correlated to the (observable) reaction time in the derivation of the model equation. In addition to providing a versatile treatment for asymmetric, sigmoidal reaction curves, another key advantage of our equation over other models is that the start time of conversion is uniquely defined at t = 0. We demonstrate the ability of our simple, two-parameter equation to successfully model the experimental x-t data for the polymorphic transformation of a pharmaceutical compound under crystallization slurry (i.e., heterogeneous) conditions. Additionally, we use a modification of this equation to model the kinetics of a historically significant, homogeneous solid-state reaction: the thermal decomposition of AgMnO4 crystals. The potential broad applicability of our statistical (i.e., dispersive) kinetic approach makes it a potentially attractive alternative to existing models/approaches.

A framework for expanding aqueous chemistry in the ...

EPA Pesticide Factsheets

This paper describes the development and implementation of an extendable aqueous-phase chemistry option (AQCHEM − KMT(I)) for the Community Multiscale Air Quality (CMAQ) modeling system, version 5.1. Here, the Kinetic PreProcessor (KPP), version 2.2.3, is used to generate a Rosenbrock solver (Rodas3) to integrate the stiff system of ordinary differential equations (ODEs) that describe the mass transfer, chemical kinetics, and scavenging processes of CMAQ clouds. CMAQ's standard cloud chemistry module (AQCHEM) is structurally limited to the treatment of a simple chemical mechanism. This work advances our ability to test and implement more sophisticated aqueous chemical mechanisms in CMAQ and further investigate the impacts of microphysical parameters on cloud chemistry. Box model cloud chemistry simulations were performed to choose efficient solver and tolerance settings, evaluate the implementation of the KPP solver, and assess the direct impacts of alternative solver and kinetic mass transfer on predicted concentrations for a range of scenarios. Month-long CMAQ simulations for winter and summer periods over the US reveal the changes in model predictions due to these cloud module updates within the full chemical transport model. While monthly average CMAQ predictions are not drastically altered between AQCHEM and AQCHEM − KMT, hourly concentration differences can be significant. With added in-cloud secondary organic aerosol (SOA) formation from bio

Nonlocal Models of Cosmic Acceleration

NASA Astrophysics Data System (ADS)

Woodard, R. P.

2014-02-01

I review a class of nonlocally modified gravity models which were proposed to explain the current phase of cosmic acceleration without dark energy. Among the topics considered are deriving causal and conserved field equations, adjusting the model to make it support a given expansion history, why these models do not require an elaborate screening mechanism to evade solar system tests, degrees of freedom and kinetic stability, and the negative verdict of structure formation. Although these simple models are not consistent with data on the growth of cosmic structures many of their features are likely to carry over to more complicated models which are in better agreement with the data.

A novel determination of calcite dissolution kinetics in seawater

NASA Astrophysics Data System (ADS)

Subhas, Adam V.; Rollins, Nick E.; Berelson, William M.; Dong, Sijia; Erez, Jonathan; Adkins, Jess F.

2015-12-01

We present a novel determination of the dissolution kinetics of inorganic calcite in seawater. We dissolved 13 C -labeled calcite in unlabeled seawater, and traced the evolving δ13 C composition of the fluid over time to establish dissolution rates. This method provides sensitive determinations of dissolution rate, which we couple with tight constraints on both seawater saturation state and surface area of the dissolving minerals. We have determined dissolution rates for two different abiotic calcite materials and three different grain sizes. Near-equilibrium dissolution rates are highly nonlinear, and are well normalized by geometric surface area, giving an empirical dissolution rate dependence on saturation state (Ω) of: This result substantiates the non-linear response of calcite dissolution to undersaturation. The bulk dissolution rate constant calculated here is in excellent agreement with those determined in far from equilibrium and dilute solution experiments. Plots of dissolution versus undersaturation indicates the presence of at least two dissolution mechanisms, implying a criticality in the calcite-seawater system. Finally, our new rate determination has implications for modeling of pelagic and seafloor dissolution. Nonlinear dissolution kinetics in a simple 1-D lysocline model indicate a possible transition from kinetic to diffusive control with increasing water depth, and also confirm the importance of respiration-driven dissolution in setting the shape of the calcite lysocline.

Ultrasonic-assisted synthesis of novel nanocomposite of poly(vinyl alcohol) and amino-modified MCM-41: A green adsorbent for Cd(II) removal.

PubMed

Soltani, Roozbeh; Dinari, Mohammad; Mohammadnezhad, Gholamhossein

2018-01-01

Amino-modified MCM-41/poly(vinyl alcohol) nanocomposite (M-MCM-41/PVOH NC) was developed for the adsorption of Cd(II) from aqueous media. M-MCM-41/PVOH NC was prepared through ultrasonic-assisted and simple blending procedure with economical and environmentally friendly polymer. The as-prepared adsorbent was characterized by FT-IR, TEM, FE-SEM and TGA. The contact time, solution pH and initial concentration of Cd(II) were found to affect the adsorption of Cd(II) from aqueous media. Kinetic studies were carried out and pseudo-first-order (PFO), pseudo-second-order (PSO), Elovich, and intra-particle diffusion (IPD) reaction kinetic models were examined. The kinetic results revealed that the adsorption of Cd(II) onto M-MCM-41/PVOH NC followed PSO kinetic model and is a chemical adsorption. The equilibrium adsorption data were evaluated by different isotherms viz. Langmuir, Freundlich, and Dubinin Radushkevich (D-R) equations. The equilibrium data fitted better with the Langmuir isotherm and the maximum adsorption capacity of M-MCM-41/PVOH NC at 298K was calculated to be 46.73mgg -1 for Cd(II) on a typical saturated monomolecular layer with a fixed number of localized adsorption sites. Copyright © 2017 Elsevier B.V. All rights reserved.

Kinetic and theoretical studies of novel biodegradable thermo-sensitive xerogels based on PEG/PVP/silica for sustained release of enrofloxacin

NASA Astrophysics Data System (ADS)

Ebadi, Azra; Rafati, Amir Abbas; Bavafa, Sadeghali; Mohammadi, Masoumah

2017-12-01

This study involves the synthesis of a new silica-based colloidal hybrid system. In this new hybrid system, poly (ethylene glycol) (PEG) and thermo-sensitive amphiphilic biocompatible poly (vinyl pyrrolidone) (PVP) were used to create suitable storage for hydrophobic drugs. The possibility of using variable PVP/PEG molar ratios to modulate drug release rate from silica nanoparticles was a primary goal of the current research. In addition, an investigation of the drug release kinetic was conducted. To achieve this, silica nanoparticles were synthesized in poly (ethylene glycol) (PEG) and poly (vinyl pyrrolidone) (PVP) solution incorporated with enrofloxacin (EFX) (as a model hydrophobic drug), using a simple synthetic strategy of hybrid materials which avoided waste and multi-step processes. The impacts of PVP/PEG molar ratios, temperature, and pH of the release medium on release kinetic were investigated. The physicochemical properties of the drug-loaded composites were studied by Fourier transform infrared (FT-IR) spectra, scanning electron microscopy (SEM), and thermogravimetric analysis (TGA). In vitro drug release studies demonstrated that the drug release rate, which was evaluated by analyzing the experimental data with seven kinetic models in a primarily non-Fickian diffusion-controlled process, aligned well with both Ritger-Peppas and Sahlin-Peppas equations.

A "Stationery" Kinetics Experiment.

ERIC Educational Resources Information Center

Hall, L.; Goberdhansingh, A.

1988-01-01

Describes a simple redox reaction that occurs between potassium permanganate and oxalic acid that can be used to prepare an interesting disappearing ink for demonstrating kinetics for introductory chemistry. Discusses laboratory procedures and factors that influence disappearance times. (CW)

Defining Nitrogen Kinetics for Air Break in Prebreath

NASA Technical Reports Server (NTRS)

Conkin, Johnny

2010-01-01

Actual tissue nitrogen (N2) kinetics are complex; the uptake and elimination is often approximated with a single half-time compartment in statistical descriptions of denitrogenation [prebreathe(PB)] protocols. Air breaks during PB complicate N2 kinetics. A comparison of symmetrical versus asymmetrical N2 kinetics was performed using the time to onset of hypobaric decompression sickness (DCS) as a surrogate for actual venous N2 tension. METHODS: Published results of 12 tests involving 179 hypobaric exposures in altitude chambers after PB, with and without airbreaks, provide the complex protocols from which to model N2 kinetics. DCS survival time for combined control and airbreaks were described with an accelerated log logistic model where N2 uptake and elimination before, during, and after the airbreak was computed with a simple exponential function or a function that changed half-time depending on ambient N2 partial pressure. P1N2-P2 = (Delta)P defined decompression dose for each altitude exposure, where P2 was the test altitude and P1N2 was computed N2 pressure at the beginning of the altitude exposure. RESULTS: The log likelihood (LL) without decompression dose (null model) was -155.6, and improved (best-fit) to -97.2 when dose was defined with a 240 min half-time for both N2 elimination and uptake during the PB. The description of DCS survival time was less precise with asymmetrical N2 kinetics, for example, LL was -98.9 with 240 min half-time elimination and 120 min half-time uptake. CONCLUSION: The statistical regression described survival time mechanistically linked to symmetrical N2 kinetics during PBs that also included airbreaks. The results are data-specific, and additional data may change the conclusion. The regression is useful to compute additional PB time to compensate for an airbreak in PB within the narrow range of tested conditions.

Defining Nitrogen Kinetics for Air Break in Prebreathe

NASA Technical Reports Server (NTRS)

Conkin, Johnny

2009-01-01

Actual tissue nitrogen (N2) kinetics are complex; the uptake and elimination is often approximated with a single half-time compartment in statistical descriptions of denitrogenation [prebreathe (PB)] protocols. Air breaks during PB complicate N2 kinetics. A comparison of symmetrical versus asymmetrical N2 kinetics was performed using the time to onset of hypobaric decompression sickness (DCS) as a surrogate for actual venous N2 tension. Published results of 12 tests involving 179 hypobaric exposures in altitude chambers after PB, with and without air breaks, provide the complex protocols from which to model N2 kinetics. DCS survival time for combined control and air breaks were described with an accelerated log logistic model where N2 uptake and elimination before, during, and after the air break was computed with a simple exponential function or a function that changed half-time depending on ambient N2 partial pressure. P1N2-P2 = delta P defined DCS dose for each altitude exposure, where P2 was the test altitude and P1N2 was computed N2 pressure at the beginning of the altitude exposure. The log likelihood (LL) without DCS dose (null model) was -155.6, and improved (best-fit) to -97.2 when dose was defined with a 240 min half-time for both N2 elimination and uptake during the PB. The description of DCS survival time was less precise with asymmetrical N2 kinetics, for example, LL was -98.9 with 240 min half-time elimination and 120 min half-time uptake. The statistical regression described survival time mechanistically linked to symmetrical N2 kinetics during PBs that also included air breaks. The results are data-specific, and additional data may change the conclusion. The regression is useful to compute additional PB time to compensate for an air break in PB within the narrow range of tested conditions.

Heat Pipe Reactor Dynamic Response Tests: SAFE-100 Reactor Core Prototype

NASA Technical Reports Server (NTRS)

Bragg-Sitton, Shannon M.

2005-01-01

The SAFE-I00a test article at the NASA Marshall Space Flight Center was used to simulate a variety of potential reactor transients; the SAFEl00a is a resistively heated, stainless-steel heat-pipe (HP)-reactor core segment, coupled to a gas-flow heat exchanger (HX). For these transients the core power was controlled by a point kinetics model with reactivity feedback based on core average temperature; the neutron generation time and the temperature feedback coefficient are provided as model inputs. This type of non-nuclear test is expected to provide reasonable approximation of reactor transient behavior because reactivity feedback is very simple in a compact fast reactor (simple, negative, and relatively monotonic temperature feedback, caused mostly by thermal expansion) and calculations show there are no significant reactivity effects associated with fluid in the HP (the worth of the entire inventory of Na in the core is .

Using Trypsin & Soybean Trypsin Inhibitor to Teach Principles of Enzyme Kinetics

ERIC Educational Resources Information Center

Howard, David R.; Herr, Julie; Hollister, Rhiannon

2006-01-01

Trypsin and soybean trypsin inhibitor (Kunitz inhibitor) can be used in a relatively simple and inexpensive student exercise to demonstrate the usefulness of enzyme kinetics. The study of enzyme kinetics is essential to biology because enzymes play such a crucial role in the biochemical pathways of all living organisms. The data from enzyme…

Detonation product EOS studies: Using ISLS to refine CHEETAH

NASA Astrophysics Data System (ADS)

Zaug, Joseph; Fried, Larry; Hansen, Donald

2001-06-01

Knowledge of an effective interatomic potential function underlies any effort to predict or rationalize the properties of solids and liquids. The experiments we undertake are directed towards determination of equilibrium and dynamic properties of simple fluids at densities sufficiently high that traditional computational methods and semi-empirical forms successful at ambient conditions may require reconsideration. In this paper we present high-pressure and temperature experimental sound speed data on a suite of non-ideal simple fluids and fluid mixtures. Impulsive Stimulated Light Scattering conducted in the diamond-anvil cell offers an experimental approach to determine cross-pair potential interactions through equation of state determinations. In addition the kinetics of structural relaxation in fluids can be studied. We compare our experimental results with our thermochemical computational model CHEETAH. Computational models are systematically improved with each addition of experimental data. Experimentally grounded computational models provide a good basis to confidently understand the chemical nature of reactions at extreme conditions.

Transporters, channels, or simple diffusion? Dogmas, atypical roles and complexity in transport systems.

PubMed

Conde, Artur; Diallinas, George; Chaumont, François; Chaves, Manuela; Gerós, Hernâni

2010-06-01

The recent breakthrough discoveries of transport systems assigned with atypical functions provide evidence for complexity in membrane transport biochemistry. Some channels are far from being simple pores creating hydrophilic passages for solutes and can, unexpectedly, act as enzymes, or mediate high-affinity uptake, and some transporters are surprisingly able to function as sensors, channels or even enzymes. Furthermore, numerous transport studies have demonstrated complex multiphasic uptake kinetics for organic and mineral nutrients. The biphasic kinetics of glucose uptake in Saccharomyces cerevisiae, a result of several genetically distinct uptake systems operating simultaneously, is a classical example that is a subject of continuous debate. In contrast, some transporters display biphasic kinetics, being bona fidae dual-affinity transporters, their kinetic properties often modulated by post-translational regulation. Also, aquaporins have recently been reported to exhibit diverse transport properties and can behave as highly adapted, multifunctional channels, transporting solutes such as CO(2), hydrogen peroxide, urea, ammonia, glycerol, polyols, carbamides, purines and pyrimidines, metalloids, glycine, and lactic acid, rather than being simple water pores. The present review provides an overview on some atypical functions displayed by transporter proteins and discusses how this novel knowledge on cellular uptake systems may be related to complex multiphasic uptake kinetics often seen in a wide variety of living organisms and the intriguing diffusive uptake of sugars and other solutes. Copyright 2009 Elsevier Ltd. All rights reserved.

Monte Carlo solution of Boltzmann equation for a simple model of highly nonequilibrium diatomic gases: Translational rotational energy relaxation

NASA Technical Reports Server (NTRS)

Yoshikawa, K. K.

1978-01-01

The semiclassical transition probability was incorporated in the simulation for energy exchange between rotational and translational energy. The results provide details on the fundamental mechanisms of gas kinetics where analytical methods were impractical. The validity of the local Maxwellian assumption and relaxation time, rotational-translational energy transition, and a velocity analysis of the inelastic collision were discussed in detail.

Single-particle trajectories reveal two-state diffusion-kinetics of hOGG1 proteins on DNA.

PubMed

Vestergaard, Christian L; Blainey, Paul C; Flyvbjerg, Henrik

2018-03-16

We reanalyze trajectories of hOGG1 repair proteins diffusing on DNA. A previous analysis of these trajectories with the popular mean-squared-displacement approach revealed only simple diffusion. Here, a new optimal estimator of diffusion coefficients reveals two-state kinetics of the protein. A simple, solvable model, in which the protein randomly switches between a loosely bound, highly mobile state and a tightly bound, less mobile state is the simplest possible dynamic model consistent with the data. It yields accurate estimates of hOGG1's (i) diffusivity in each state, uncorrupted by experimental errors arising from shot noise, motion blur and thermal fluctuations of the DNA; (ii) rates of switching between states and (iii) rate of detachment from the DNA. The protein spends roughly equal time in each state. It detaches only from the loosely bound state, with a rate that depends on pH and the salt concentration in solution, while its rates for switching between states are insensitive to both. The diffusivity in the loosely bound state depends primarily on pH and is three to ten times higher than in the tightly bound state. We propose and discuss some new experiments that take full advantage of the new tools of analysis presented here.

Facile fabrication of superparamagnetic graphene/polyaniline/Fe3O4 nanocomposites for fast magnetic separation and efficient removal of dye.

PubMed

Mu, Bin; Tang, Jie; Zhang, Long; Wang, Aiqin

2017-07-13

Using graphene as adsorbent for removal of pollutants from polluted water is commonly recognized to be costly because the graphene is usually produced by a very complex process. Herein, a simple and eco-friendly method was employed to fabricate efficient superparamagnetic graphene/polyaniline/Fe 3 O 4 nanocomposites for removal of dyes. The exfoliation of graphite as nanosheets and the functionalization of nanosheets with polyaniline and Fe 3 O 4 nanoparticles were simultaneously achieved via a one-pot reaction process combining the intercalation polymerization of aniline and the co-precipitation of the residual Fe 3+ and the generated Fe 2+ . The obtained graphene/polyaniline/Fe 3 O 4 nanocomposites exhibited excellent adsorption performance for Congo red, even in the presence of Brilliant green. The adsorption kinetics and adsorption isotherms were well fitted with pseudo second-order kinetic model and Langmuir isotherm model, respectively. In a word, this method is simple and industrially feasible, which provides a new approach to fabricate highly efficient graphene-based adsorbents on large scale for removal of dyes. In addition, it also can be used to exfoliate other two-dimensional materials, such as boron nitride, carbon nitride and MoS 2 for a range of possible applications.

An Example of Following Second-Order Kinetics by Simple Laboratory Means

ERIC Educational Resources Information Center

Schreiber, Gisela

1976-01-01

Describes a procedure for studying the kinetics of the second-order hydrolysis of ethylene bromohydrine in alkaline medium by incorporating a substance that changes color as one of the reacting components is depleted. (MLH)

An algorithm for the kinetics of tire pyrolysis under different heating rates.

PubMed

Quek, Augustine; Balasubramanian, Rajashekhar

2009-07-15

Tires exhibit different kinetic behaviors when pyrolyzed under different heating rates. A new algorithm has been developed to investigate pyrolysis behavior of scrap tires. The algorithm includes heat and mass transfer equations to account for the different extents of thermal lag as the tire is heated at different heating rates. The algorithm uses an iterative approach to fit model equations to experimental data to obtain quantitative values of kinetic parameters. These parameters describe the pyrolysis process well, with good agreement (r(2)>0.96) between the model and experimental data when the model is applied to three different brands of automobile tires heated under five different heating rates in a pure nitrogen atmosphere. The model agrees with other researchers' results that frequencies factors increased and time constants decreased with increasing heating rates. The model also shows the change in the behavior of individual tire components when the heating rates are increased above 30 K min(-1). This result indicates that heating rates, rather than temperature, can significantly affect pyrolysis reactions. This algorithm is simple in structure and yet accurate in describing tire pyrolysis under a wide range of heating rates (10-50 K min(-1)). It improves our understanding of the tire pyrolysis process by showing the relationship between the heating rate and the many components in a tire that depolymerize as parallel reactions.

Principal component analysis acceleration of rovibrational coarse-grain models for internal energy excitation and dissociation

NASA Astrophysics Data System (ADS)

Bellemans, Aurélie; Parente, Alessandro; Magin, Thierry

2018-04-01

The present work introduces a novel approach for obtaining reduced chemistry representations of large kinetic mechanisms in strong non-equilibrium conditions. The need for accurate reduced-order models arises from compression of large ab initio quantum chemistry databases for their use in fluid codes. The method presented in this paper builds on existing physics-based strategies and proposes a new approach based on the combination of a simple coarse grain model with Principal Component Analysis (PCA). The internal energy levels of the chemical species are regrouped in distinct energy groups with a uniform lumping technique. Following the philosophy of machine learning, PCA is applied on the training data provided by the coarse grain model to find an optimally reduced representation of the full kinetic mechanism. Compared to recently published complex lumping strategies, no expert judgment is required before the application of PCA. In this work, we will demonstrate the benefits of the combined approach, stressing its simplicity, reliability, and accuracy. The technique is demonstrated by reducing the complex quantum N2(g+1Σ) -N(S4u ) database for studying molecular dissociation and excitation in strong non-equilibrium. Starting from detailed kinetics, an accurate reduced model is developed and used to study non-equilibrium properties of the N2(g+1Σ) -N(S4u ) system in shock relaxation simulations.

Hydrogen Donor-Acceptor Fluctuations from Kinetic Isotope Effects: A Phenomenological Model

PubMed Central

Roston, Daniel; Cheatum, Christopher M.; Kohen, Amnon

2012-01-01

Kinetic isotope effects (KIEs) and their temperature dependence can probe the structural and dynamic nature of enzyme-catalyzed proton or hydride transfers. The molecular interpretation of their temperature dependence requires expensive and specialized QM/MM calculations to provide a quantitative molecular understanding. Currently available phenomenological models use a non-adiabatic assumption that is not appropriate for most hydride and proton-transfer reactions, while others require more parameters than the experimental data justify. Here we propose a phenomenological interpretation of KIEs based on a simple method to quantitatively link the size and temperature dependence of KIEs to a conformational distribution of the catalyzed reaction. The present model assumes adiabatic hydrogen tunneling, and by fitting experimental KIE data, the model yields a population distribution for fluctuations of the distance between donor and acceptor atoms. Fits to data from a variety of proton and hydride transfers catalyzed by enzymes and their mutants, as well as non-enzymatic reactions, reveal that steeply temperature-dependent KIEs indicate the presence of at least two distinct conformational populations, each with different kinetic behaviors. We present the results of these calculations for several published cases and discuss how the predictions of the calculations might be experimentally tested. The current analysis does not replace molecular quantum mechanics/molecular mechanics (QM/MM) investigations, but it provides a fast and accessible way to quantitatively interpret KIEs in the context of a Marcus-like model. PMID:22857146

Sensitivity of Polar Stratospheric Ozone Loss to Uncertainties in Chemical Reaction Kinetics

NASA Technical Reports Server (NTRS)

Kawa, S. Randolph; Stolarksi, Richard S.; Douglass, Anne R.; Newman, Paul A.

2008-01-01

Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspects of our understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to simulate these processes in numerical models of chemistry and transport. The fidelity of the models is assessed in comparison with a wide range of observations. These models depend on laboratory-measured kinetic reaction rates and photolysis cross sections to simulate molecular interactions. A typical stratospheric chemistry mechanism has on the order of 50- 100 species undergoing over a hundred intermolecular reactions and several tens of photolysis reactions. The rates of all of these reactions are subject to uncertainty, some substantial. Given the complexity of the models, however, it is difficult to quantify uncertainties in many aspects of system. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluations are applied in random combinations. We determine the key reactions and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.

Modeling a photoinduced planar-to-homeotropic anchoring transition triggered by surface azobenzene units in a nematic liquid crystal

NASA Astrophysics Data System (ADS)

Oswald, Patrick; Ignés-Mullol, Jordi

2017-09-01

The performance of light-controlled liquid crystal anchoring surfaces depends on the nature of the photosensitive moieties and on the concentration of spacer units. Here, we study the kinetics of photosensitive liquid crystal cells that incorporate an azobenzene-based self-assembled monolayer. We characterize the photoinduced homeotropic-to-planar transition and the subsequent reverse relaxation in terms of the underlying isomerization of the photosensitive layer. We show that the response time can be precisely adjusted by tuning the lateral packing of azobenzene units by means of inert spacer molecules. Using simple kinetic assumptions and a well-known model for the energetics of liquid crystal anchoring we are able to capture the details of the optical microscopy experimental observations. Our analysis provides fitted values for all the relevant material parameters, including the zenithal and the azimuthal anchoring strength.

Kinetic analysis of competition between aerosol particle removal and generation by ionization air purifiers.

PubMed

Alshawa, Ahmad; Russell, Ashley R; Nizkorodov, Sergey A

2007-04-01

Ionization air purifiers are increasingly used to remove aerosol particles from indoor air. However, certain ionization air purifiers also emit ozone. Reactions between the emitted ozone and unsaturated volatile organic compounds (VOC) commonly found in indoor air produce additional respirable aerosol particles in the ultrafine (<0.1 microm) and fine (<2.5 microm) size domains. A simple kinetic model is used to analyze the competition between the removal and generation of particulate matter by ionization air purifiers under conditions of a typical residential building. This model predicts that certain widely used ionization air purifiers may actually increase the mass concentration of fine and ultrafine particulates in the presence of common unsaturated VOC, such as limonene contained in many household cleaning products. This prediction is supported by an explicit observation of ultrafine particle nucleation events caused by the addition of D-limonene to a ventilated office room equipped with a common ionization air purifier.

Solute segregation kinetics and dislocation depinning in a binary alloy

NASA Astrophysics Data System (ADS)

Dontsova, E.; Rottler, J.; Sinclair, C. W.

2015-06-01

Static strain aging, a phenomenon caused by diffusion of solute atoms to dislocations, is an important contributor to the strength of substitutional alloys. Accurate modeling of this complex process requires both atomic spatial resolution and diffusional time scales, which is very challenging to achieve with commonly used atomistic computational methods. In this paper, we use the recently developed "diffusive molecular dynamics" (DMD) method that is capable of describing the kinetics of the solute segregation process at the atomic level while operating on diffusive time scales in a computationally efficient way. We study static strain aging in the Al-Mg system and calculate the depinning shear stress between edge and screw dislocations and their solute atmospheres formed for various waiting times with different solute content and for a range of temperatures. A simple phenomenological model is also proposed that describes the observed behavior of the critical shear stress as a function of segregation level.

Kinetics and modeling of hexavalent chromium reduction in Enterobacter cloacae

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yamamoto, Koji; Kato, Junichi; Yano, Takuo

1993-01-05

Kinetics of bacterial reduction of toxic hexavalent chromium (chromate: CrO[sub 4][sup [minus]2]) was investigated using batch and fed-batch cultures of Enterobacter cloacae strain HO1. In fed-batch cultures, the CrO[sub 4][sup [minus]2] feed was controlled on the basis of the rate of pH change. This control strategy has proven to be useful for avoiding toxic CrO[sub 3][sup [minus]2] overload. A simple mathematical model was developed to describe the bacterial process of CrO[sub 4][sup [minus]2] reduction. In this model, two types of bacterial cells were considered: induced, CrO[sub 4][sup [minus]2]-resistant cells and uninduced, sensitive ones. Only resistant cells were assumed to bemore » able to reduce CrO[sub 4][sup [minus]2]. These fundamental ideas were supported by the model predictions which well approximated all experimental data. In a simulation study, the model was also used to optimize fed-batch cultures, instead of lengthy and expensive laboratory experiments.« less

Modeling fixed and fluidized reactors for cassava starch Saccharification with immobilized enzyme

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zanin, G.M.; De Moraes, F.F.

1997-12-31

Cassava starch saccharification in fixed-and fluidized-bed reactors using immobilized enzyme was modeled in a previous paper using a simple model in which all dextrins were grouped in a single substrate. In that case, although good fit of the model to experimental data was obtained, physical inconsistency appeared as negative kinetic constants. In this work, a multisubstrate model, developed earlier for saccharification with free enzyme, is adapted for immobilized enzyme. This latter model takes into account the formation of intermediate substrates, which are dextrins competing for the catalytic site of the enzyme, reversibility of some reactions, inhibition by substrate and product,more » and the formation of isomaltose. Kinetic parameters to be used with this model were obtained from initial velocity saccharification tests using the immobilized enzyme and different liquefied starch concentrations. The new model was found to be valid for modeling both fixed- and fluidized-bed reactors. It did not present inconsistencies as the earlier one had and has shown that apparent glucose inhibition is about seven times higher in the fixed-bed than in fluidized-bed reactor. 13 refs., 5 figs., 1 tab.« less

Cluster kinetics model of particle separation in vibrated granular media.

PubMed

McCoy, Benjamin J; Madras, Giridhar

2006-01-01

We model the Brazil-nut effect (BNE) by hypothesizing that granules form clusters that fragment and aggregate. This provides a heterogeneous medium in which the immersed intruder particle rises (BNE) or sinks (reverse BNE) according to relative convection currents and buoyant and drag forces. A simple relationship proposed for viscous drag in terms of the vibrational intensity and the particle to grain density ratio allows simulation of published experimental data for rise and sink times as functions of particle radius, initial depth of the particle, and particle-grain density ratio. The proposed model correctly describes the experimentally observed maximum in risetime.

On some stochastic formulations and related statistical moments of pharmacokinetic models.

PubMed

Matis, J H; Wehrly, T E; Metzler, C M

1983-02-01

This paper presents the deterministic and stochastic model for a linear compartment system with constant coefficients, and it develops expressions for the mean residence times (MRT) and the variances of the residence times (VRT) for the stochastic model. The expressions are relatively simple computationally, involving primarily matrix inversion, and they are elegant mathematically, in avoiding eigenvalue analysis and the complex domain. The MRT and VRT provide a set of new meaningful response measures for pharmacokinetic analysis and they give added insight into the system kinetics. The new analysis is illustrated with an example involving the cholesterol turnover in rats.

Dependence of cancer cell adhesion kinetics on integrin ligand surface density measured by a high-throughput label-free resonant waveguide grating biosensor.

PubMed

Orgovan, Norbert; Peter, Beatrix; Bősze, Szilvia; Ramsden, Jeremy J; Szabó, Bálint; Horvath, Robert

2014-02-07

A novel high-throughput label-free resonant waveguide grating (RWG) imager biosensor, the Epic® BenchTop (BT), was utilized to determine the dependence of cell spreading kinetics on the average surface density (v(RGD)) of integrin ligand RGD-motifs. v(RGD) was tuned over four orders of magnitude by co-adsorbing the biologically inactive PLL-g-PEG and the RGD-functionalized PLL-g-PEG-RGD synthetic copolymers from their mixed solutions onto the sensor surface. Using highly adherent human cervical tumor (HeLa) cells as a model system, cell adhesion kinetic data of unprecedented quality were obtained. Spreading kinetics were fitted with the logistic equation to obtain the spreading rate constant (r) and the maximum biosensor response (Δλmax), which is assumed to be directly proportional to the maximum spread contact area (Amax). r was found to be independent of the surface density of integrin ligands. In contrast, Δλmax increased with increasing RGD surface density until saturation at high densities. Interpreting the latter behavior with a simple kinetic mass action model, a 2D dissociation constant of 1753 ± 243 μm(-2) (corresponding to a 3D dissociation constant of ~30 μM) was obtained for the binding between RGD-specific integrins embedded in the cell membrane and PLL-g-PEG-RGD. All of these results were obtained completely noninvasively without using any labels.

Adsorptive Removal of Toxic Chromium from Waste-Water Using Wheat Straw and Eupatorium adenophorum

PubMed Central

Song, Dagang; Pan, Kaiwen; Tariq, Akash; Azizullah, Azizullah; Sun, Feng; Li, Zilong; Xiong, Qinli

2016-01-01

Environmental pollution with heavy metals is a serious issue worldwide posing threats to humans, animals and plants and to the stability of overall ecosystem. Chromium (Cr) is one of most hazardous heavy metals with a high carcinogenic and recalcitrant nature. Aim of the present study was to select low-cost biosorbent using wheat straw and Eupatorium adenophorum through simple carbonization process, capable of removing Cr (VI) efficiently from wastewater. From studied plants a low cost adsorbent was prepared for removing Cr (VI) from aqueous solution following very simple carbonization method excluding activation process. Several factors such as pH, contact time, sorbent dosage and temperature were investigated for attaining ideal condition. For analysis of adsorption equilibrium isotherm data, Langmuir, Freundlich and Temkin models were used while pseudo-first-order, pseudo-second-order, external diffusion and intra-particle diffusion models were used for the analysis of kinetic data. The obtained results revealed that 99.9% of Cr (VI) removal was observed in the solution with a pH of 1.0. Among all the tested models Langmuir model fitted more closely according to the data obtained. Increase in adsorption capacity was observed with increasing temperature revealing endothermic nature of Cr (VI). The maximum Cr (VI) adsorption potential of E. adenophorum and wheat straw was 89.22 mg per 1 gram adsorbent at 308K. Kinetic data of absorption precisely followed pseudo-second-order model. Present study revealed highest potential of E. adenophorum and wheat straw for producing low cost adsorbent and to remove Cr (VI) from contaminated water. PMID:27911906

Numerical simulation by the molecular collision theory of two-phase mixture explosion characteristics in closed or vented vessels

NASA Astrophysics Data System (ADS)

Pascaud, J. M.; Brossard, J.; Lombard, J. M.

1999-09-01

The aim of this work consists in presenting a simple modelling (the molecular collision theory), easily usable in an industrial environment in order to predict the evolution of thermodynamical characteristics of the combustion of two-phase mixtures in a closed or a vented vessel. Basic characteristics of the modelling have been developed for ignition and combustion of propulsive powders and adapted with appropriate parameters linked to simplified kinetics. A simple representation of the combustion phenomena based on energy transfers and the action of specific molecules is presented. The model is generalized to various mixtures such as dust suspensions, liquid fuel drops and hybrid mixtures composed of dust and a gaseous supply such as methane or propane in the general case of vented explosions. The pressure venting due to the vent breaking is calculated from thermodynamical characteristics given by the model and taking into account, the mass rate of discharge of the different products deduced from the standard orifice equations. The application conditions determine the fuel ratio of the used mixtures, the nature of the chemical kinetics and the calculation of a universal set of parameters. The model allows to study the influence of the fuel concentration and the supply of gaseous additives, the influence of the vessel volume (2400ell leq V_bleq 250 000ell) and the influence of the venting pressure or the vent area. The first results have been compared with various experimental works available for two phase mixtures and indicate quite correct predictions.

Adsorptive Removal of Toxic Chromium from Waste-Water Using Wheat Straw and Eupatorium adenophorum.

PubMed

Song, Dagang; Pan, Kaiwen; Tariq, Akash; Azizullah, Azizullah; Sun, Feng; Li, Zilong; Xiong, Qinli

2016-01-01

Environmental pollution with heavy metals is a serious issue worldwide posing threats to humans, animals and plants and to the stability of overall ecosystem. Chromium (Cr) is one of most hazardous heavy metals with a high carcinogenic and recalcitrant nature. Aim of the present study was to select low-cost biosorbent using wheat straw and Eupatorium adenophorum through simple carbonization process, capable of removing Cr (VI) efficiently from wastewater. From studied plants a low cost adsorbent was prepared for removing Cr (VI) from aqueous solution following very simple carbonization method excluding activation process. Several factors such as pH, contact time, sorbent dosage and temperature were investigated for attaining ideal condition. For analysis of adsorption equilibrium isotherm data, Langmuir, Freundlich and Temkin models were used while pseudo-first-order, pseudo-second-order, external diffusion and intra-particle diffusion models were used for the analysis of kinetic data. The obtained results revealed that 99.9% of Cr (VI) removal was observed in the solution with a pH of 1.0. Among all the tested models Langmuir model fitted more closely according to the data obtained. Increase in adsorption capacity was observed with increasing temperature revealing endothermic nature of Cr (VI). The maximum Cr (VI) adsorption potential of E. adenophorum and wheat straw was 89.22 mg per 1 gram adsorbent at 308K. Kinetic data of absorption precisely followed pseudo-second-order model. Present study revealed highest potential of E. adenophorum and wheat straw for producing low cost adsorbent and to remove Cr (VI) from contaminated water.

Experimental validation of an analytical kinetic model for edge-localized modes in JET-ITER-like wall

NASA Astrophysics Data System (ADS)

Guillemaut, C.; Metzger, C.; Moulton, D.; Heinola, K.; O’Mullane, M.; Balboa, I.; Boom, J.; Matthews, G. F.; Silburn, S.; Solano, E. R.; contributors, JET

2018-06-01

The design and operation of future fusion devices relying on H-mode plasmas requires reliable modelling of edge-localized modes (ELMs) for precise prediction of divertor target conditions. An extensive experimental validation of simple analytical predictions of the time evolution of target plasma loads during ELMs has been carried out here in more than 70 JET-ITER-like wall H-mode experiments with a wide range of conditions. Comparisons of these analytical predictions with diagnostic measurements of target ion flux density, power density, impact energy and electron temperature during ELMs are presented in this paper and show excellent agreement. The analytical predictions tested here are made with the ‘free-streaming’ kinetic model (FSM) which describes ELMs as a quasi-neutral plasma bunch expanding along the magnetic field lines into the Scrape-Off Layer without collisions. Consequences of the FSM on energy reflection and deposition on divertor targets during ELMs are also discussed.

Kinetic Features in the Ion Flux Spectrum

NASA Astrophysics Data System (ADS)

Vafin, S.; Riazantseva, M.; Yoon, P. H.

2017-11-01

An interesting feature of solar wind fluctuations is the occasional presence of a well-pronounced peak near the spectral knee. These peaks are well investigated in the context of magnetic field fluctuations in the magnetosheath and they are typically related to kinetic plasma instabilities. Recently, similar peaks were observed in the spectrum of ion flux fluctuations of the solar wind and magnetosheath. In this paper, we propose a simple analytical model to describe such peaks in the ion flux spectrum based on the linear theory of plasma fluctuations. We compare our predictions with a sample observation in the solar wind. For the given observation, the peak requires ˜10 minutes to grow up to the observed level that agrees with the quasi-linear relaxation time. Moreover, our model well reproduces the form of the measured peak in the ion flux spectrum. The observed lifetime of the peak is about 50 minutes, which is relatively close to the nonlinear Landau damping time of 30-40 minutes. Overall, our model proposes a plausible scenario explaining the observation.

Cr(OH)3-NPs-CNC hybrid nanocomposite: a sorbent for adsorptive removal of methylene blue and malachite green from solutions.

PubMed

Nekouei, Farzin; Nekouei, Shahram; Keshtpour, Farzaneh; Noorizadeh, Hossein; Wang, Shaobin

2017-11-01

In this article, Cr(OH) 3 nanoparticle-modified cellulose nanocrystal (CNC) as a novel hybrid nanocomposite (Cr(OH) 3 -NPs-CNC) was prepared by a simple procedure and used as a sorbent for adsorptive removal of methylene blue (MB) and malachite green (MG) from aqueous solution. Different kinetic models were tested, and the pseudo-second-order kinetic model was found more suitable for the MB and MG adsorption processes. The BET and Langmuir models were more suitable for the adsorption processes of MB and MG. Thermodynamic studies suggested that the adsorption of MB and MG onto Cr(OH) 3 -NPs-CNC nanocomposite was a spontaneous and endothermic process. The maximum adsorption capacities for MB and MG were reached 106 and 104 mg/g, respectively, which were almost two times higher than unmodified CNC. The chemical stability and leaching tests of the Cr(OH) 3 -NPs-CNC hybrid nanocomposite showed that only small amounts of chromium were leached into the solution.

A unified model of Grignard reagent formation.

PubMed

Shao, Yunqi; Liu, Zhen; Huang, Pan; Liu, Boping

2018-04-25

Grignard reagents are among the most fundamental reagents in organic synthesis, yet studies have hitherto failed to fully explain the selectivity and kinetics of Grignard reagent formation (GRF). The present study provides new insights into the intermediates and pathways of GRF using density functional theory (DFT) calculations. Potential energy surfaces of RX dissociation along different directions reveal the origin of configuration retention of alkenyl and aromatic halides. Radical intermediates participate solely in the dissociation stage, and depend on the geometry of the reactant halide. Dissociation of organic halides yields stabilized surface anions, and the rest of the reaction is ionic in nature. MgX+/RMg+ were proposed as the key intermediates of Mg leaving from the surface in the self-activation of GRF, which explains the accelerated kinetics upon addition of RMgX or MgX2. The intermediacy of the cations was supported by a simple electrochemical experiment. To the best of our knowledge, this is the first unified ionic model (I-model) developed for resolving the controversial issues of GRF.

Nonlocal kinetic energy functional from the jellium-with-gap model: Applications to orbital-free density functional theory

NASA Astrophysics Data System (ADS)

Constantin, Lucian A.; Fabiano, Eduardo; Della Sala, Fabio

2018-05-01

Orbital-free density functional theory (OF-DFT) promises to describe the electronic structure of very large quantum systems, being its computational cost linear with the system size. However, the OF-DFT accuracy strongly depends on the approximation made for the kinetic energy (KE) functional. To date, the most accurate KE functionals are nonlocal functionals based on the linear-response kernel of the homogeneous electron gas, i.e., the jellium model. Here, we use the linear-response kernel of the jellium-with-gap model to construct a simple nonlocal KE functional (named KGAP) which depends on the band-gap energy. In the limit of vanishing energy gap (i.e., in the case of metals), the KGAP is equivalent to the Smargiassi-Madden (SM) functional, which is accurate for metals. For a series of semiconductors (with different energy gaps), the KGAP performs much better than SM, and results are close to the state-of-the-art functionals with sophisticated density-dependent kernels.

A Simple Parameterization of Mixing of Passive Scalars in Turbulent Flows

NASA Astrophysics Data System (ADS)

Nithianantham, Ajithshanthar; Venayagamoorthy, Karan

2015-11-01

A practical model for quantifying the turbulent diascalar diffusivity is proposed as Ks = 1 . 1γ'LTk 1 / 2 , where LT is defined as the Thorpe length scale, k is the turbulent kinetic energy and γ' is one-half of the mechanical to scalar time scale ratio, which was shown by previous researchers to be approximately 0 . 7 . The novelty of the proposed model lies in the use of LT, which is a widely used length scale in stably stratified flows (almost exclusively used in oceanography), for quantifying turbulent mixing in unstratified flows. LT can be readily obtained in the field using a Conductivity, Temperature and Depth (CTD) profiler. The turbulent kinetic energy is mostly contained in the large scales of the flow field and hence can be measured in the field or modeled in numerical simulations. Comparisons using DNS data show remarkably good agreement between the predicted and exact diffusivities. Office of Naval Research and National Science Foundation.

A radially resolved kinetic model for nonlocal electron ripple diffusion losses in tokamaks

DOE Office of Scientific and Technical Information (OSTI.GOV)

Robertson, Scott

A relatively simple radially resolved kinetic model is applied to the ripple diffusion problem for electrons in tokamaks. The distribution function f(r,v) is defined on a two-dimensional grid, where r is the radial coordinate and v is the velocity coordinate. Particle transport in the radial direction is from ripple and banana diffusion and transport in the velocity direction is described by the Fokker-Planck equation. Particles and energy are replaced by source functions that are adjusted to maintain a constant central density and temperature. The relaxed profiles of f(r,v) show that the electron distribution function at the wall contains suprathermal electronsmore » that have diffused from the interior that enhance ripple transport. The transport at the periphery is therefore nonlocal. The energy replacement times from the computational model are near to the experimental replacement times for tokamak discharges in the compilation by Pfeiffer and Waltz [Nucl. Fusion 19, 51 (1979)].« less

Energy Conservation and Conversion in NIMROD Computations of Magnetic Reconnection

NASA Astrophysics Data System (ADS)

Maddox, J. A.; Sovinec, C. R.

2017-10-01

Previous work modeling magnetic relaxation during non-inductive start-up at the Pegasus spherical tokamak indicates an order of magnitude gap between measured experimental temperature and simulated temperature in NIMROD. Potential causes of the plasma temperature gap include: insufficient transport modeling, too low modeled injector power input, and numerical loss of energy, as energy is not algorithmically conserved in NIMROD simulations. Simple 2D nonlinear MHD simulations explore numerical energy conservation discrepancies in NIMROD because understanding numerical loss of energy is fundamental to addressing the physical problems of the other potential causes of energy loss. Evolution of these configurations induces magnetic reconnection, which transfers magnetic energy to heat and kinetic energy. The kinetic energy is eventually damped so, magnetic energy loss must correspond to an increase in internal energy. Results in the 2D geometries indicate that numerical energy loss during reconnection depends on the temporal resolution of the dynamics. Work support from U.S. Department of Energy through a subcontract from the Plasma Science and Innovation Center.

Shear viscosity for dense plasmas by equilibrium molecular dynamics in asymmetric Yukawa ionic mixtures

DOE PAGES

Haxhimali, Tomorr; Rudd, Robert E.; Cabot, William H.; ...

2015-11-24

We present molecular dynamics (MD) calculations of shear viscosity for asymmetric mixed plasma for thermodynamic conditions relevant to astrophysical and inertial confinement fusion plasmas. Specifically, we consider mixtures of deuterium and argon at temperatures of 100–500 eV and a number density of 10 25 ions/cc. The motion of 30 000–120 000 ions is simulated in which the ions interact via the Yukawa (screened Coulomb) potential. The electric field of the electrons is included in this effective interaction; the electrons are not simulated explicitly. Shear viscosity is calculated using the Green-Kubo approach with an integral of the shear stress autocorrelation function,more » a quantity calculated in the equilibrium MD simulations. We systematically study different mixtures through a series of simulations with increasing fraction of the minority high- Z element (Ar) in the D-Ar plasma mixture. In the more weakly coupled plasmas, at 500 eV and low Ar fractions, results from MD compare very well with Chapman-Enskog kinetic results. In the more strongly coupled plasmas, the kinetic theory does not agree well with the MD results. Here, we develop a simple model that interpolates between classical kinetic theories at weak coupling and the Murillo Yukawa viscosity model at higher coupling. Finally, this hybrid kinetics-MD viscosity model agrees well with the MD results over the conditions simulated, ranging from moderately weakly coupled to moderately strongly coupled asymmetric plasma mixtures.« less

Relaxational effects in radiating stellar collapse

NASA Astrophysics Data System (ADS)

Govender, Megan; Maartens, Roy; Maharaj, Sunil D.

1999-12-01

Relaxational effects in stellar heat transport can in many cases be significant. Relativistic Fourier-Eckart theory is inherently quasi-stationary, and cannot incorporate these effects. The effects are naturally accounted for in causal relativistic thermodynamics, which provides an improved approximation to kinetic theory. Recent results, based on perturbations of a static star, show that relaxation effects can produce a significant increase in the central temperature and temperature gradient for a given luminosity. We use a simple stellar model that allows for non-perturbative deviations from staticity, and confirms qualitatively the predictions of the perturbative models.

Chemical control of the viscoelastic properties of vinylogous urethane vitrimers

PubMed Central

Denissen, Wim; Droesbeke, Martijn; Nicolaÿ, Renaud; Leibler, Ludwik; Winne, Johan M.; Du Prez, Filip E.

2017-01-01

Vinylogous urethane based vitrimers are polymer networks that have the intrinsic property to undergo network rearrangements, stress relaxation and viscoelastic flow, mediated by rapid addition/elimination reactions of free chain end amines. Here we show that the covalent exchange kinetics significantly can be influenced by combination with various simple additives. As anticipated, the exchange reactions on network level can be further accelerated using either Brønsted or Lewis acid additives. Remarkably, however, a strong inhibitory effect is observed when a base is added to the polymer matrix. These effects have been mechanistically rationalized, guided by low-molecular weight kinetic model experiments. Thus, vitrimer elastomer materials can be rationally designed to display a wide range of viscoelastic properties. PMID:28317893

Folding and unfolding single RNA molecules under tension

PubMed Central

Woodside, Michael T; García-García, Cuauhtémoc; Block, Steven M

2010-01-01

Single-molecule force spectroscopy constitutes a powerful method for probing RNA folding: it allows the kinetic, energetic, and structural properties of intermediate and transition states to be determined quantitatively, yielding new insights into folding pathways and energy landscapes. Recent advances in experimental and theoretical methods, including fluctuation theorems, kinetic theories, novel force clamps, and ultrastable instruments, have opened new avenues for study. These tools have been used to probe folding in simple model systems, for example, RNA and DNA hairpins. Knowledge gained from such systems is helping to build our understanding of more complex RNA structures composed of multiple elements, as well as how nucleic acids interact with proteins involved in key cellular activities, such as transcription and translation. PMID:18786653

Quantitative Detection of Nucleoside Analogues by Multi-enzyme Biosensors using Time-Resolved Kinetic Measurements.

PubMed

Muthu, Pravin; Lutz, Stefan

2016-04-05

Fast, simple and cost-effective methods for detecting and quantifying pharmaceutical agents in patients are highly sought after to replace equipment and labor-intensive analytical procedures. The development of new diagnostic technology including portable detection devices also enables point-of-care by non-specialists in resource-limited environments. We have focused on the detection and dose monitoring of nucleoside analogues used in viral and cancer therapies. Using deoxyribonucleoside kinases (dNKs) as biosensors, our chemometric model compares observed time-resolved kinetics of unknown analytes to known substrate interactions across multiple enzymes. The resulting dataset can simultaneously identify and quantify multiple nucleosides and nucleoside analogues in complex sample mixtures. © 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

A simple recipe for setting up the flux equations of cyclic and linear reaction schemes of ion transport with a high number of states: The arrow scheme.

PubMed

Hansen, Ulf-Peter; Rauh, Oliver; Schroeder, Indra

2016-01-01

The calculation of flux equations or current-voltage relationships in reaction kinetic models with a high number of states can be very cumbersome. Here, a recipe based on an arrow scheme is presented, which yields a straightforward access to the minimum form of the flux equations and the occupation probability of the involved states in cyclic and linear reaction schemes. This is extremely simple for cyclic schemes without branches. If branches are involved, the effort of setting up the equations is a little bit higher. However, also here a straightforward recipe making use of so-called reserve factors is provided for implementing the branches into the cyclic scheme, thus enabling also a simple treatment of such cases.

A simple recipe for setting up the flux equations of cyclic and linear reaction schemes of ion transport with a high number of states: The arrow scheme

PubMed Central

Hansen, Ulf-Peter; Rauh, Oliver; Schroeder, Indra

2016-01-01

abstract The calculation of flux equations or current-voltage relationships in reaction kinetic models with a high number of states can be very cumbersome. Here, a recipe based on an arrow scheme is presented, which yields a straightforward access to the minimum form of the flux equations and the occupation probability of the involved states in cyclic and linear reaction schemes. This is extremely simple for cyclic schemes without branches. If branches are involved, the effort of setting up the equations is a little bit higher. However, also here a straightforward recipe making use of so-called reserve factors is provided for implementing the branches into the cyclic scheme, thus enabling also a simple treatment of such cases. PMID:26646356

Development of the Concept of Energy Conservation using Simple Experiments for Grade 10 Students

NASA Astrophysics Data System (ADS)

Rachniyom, S.; Toedtanya, K.; Wuttiprom, S.

2017-09-01

The purpose of this research was to develop students’ concept of and retention rate in relation to energy conservation. Activities included simple and easy experiments that considered energy transformation from potential to kinetic energy. The participants were 30 purposively selected grade 10 students in the second semester of the 2016 academic year. The research tools consisted of learning lesson plans and a learning achievement test. Results showed that the experiments worked well and were appropriate as learning activities. The students’ achievement scores significantly increased at the statistical level of 05, the students’ retention rates were at a high level, and learning behaviour was at a good level. These simple experiments allowed students to learn to demonstrate to their peers and encouraged them to use familiar models to explain phenomena in daily life.

Degradation of organic pollutants by Vacuum-Ultraviolet (VUV): Kinetic model and efficiency.

PubMed

Xie, Pengchao; Yue, Siyang; Ding, Jiaqi; Wan, Ying; Li, Xuchun; Ma, Jun; Wang, Zongping

2018-04-15

Vacuum-Ultraviolet (VUV), an efficient and green method to produce hydroxyl radical (•OH), is effective in degrading numerous organic contaminants in aqueous solution. Here, we proposed an effective and simple kinetic model to describe the degradation of organic pollutants in VUV system, by taking the •OH scavenging effects of formed organic intermediates as co-existing organic matter in whole. Using benzoic acid (BA) as a •OH probe, •OH was regarded vital for pollutant degradation in VUV system, and the thus developed model successfully predicted its degradation kinetics under different conditions. Effects of typical influencing factors such as BA concentrations and UV intensity were investigated quantitatively by the model. Temperature was found to be an important influencing factor in the VUV system, and the quantum yield of •OH showed a positive linear dependence on temperature. Impacts of humic acid (HA), alkalinity, chloride, and water matrices (realistic waters) on the oxidation efficiency were also examined. BA degradation was significantly inhibited by HA due to its scavenging of •OH, but was influenced much less by the alkalinity and chloride; high oxidation efficiency was still obtained in the realistic water. The degradation kinetics of three other typical micropollutants including bisphenol A (BPA), nitrobenzene (NB) and dimethyl phthalate (DMP), and the mixture of co-existing BA, BPA and DMP were further studied, and the developed model predicted the experimental data well, especially in realistic water. It is expected that this study will provide an effective approach to predict the degradation of organic micropollutants by the promising VUV system, and broaden the application of VUV system in water treatment. Copyright © 2018 Elsevier Ltd. All rights reserved.

Detonation Product EOS Studies: Using ISLS to Refine Cheetah

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zaug, J M; Howard, W M; Fried, L E

2001-08-08

Knowledge of an effective interatomic potential function underlies any effort to predict or rationalize the properties of solids and liquids. The experiments we undertake are directed towards determination of equilibrium and dynamic properties of simple fluids at densities sufficiently high that traditional computational methods and semi-empirical forms successful at ambient conditions may require reconsideration. In this paper we present high-pressure and temperature experimental sound speed data on a simple fluid, methanol. Impulsive Stimulated Light Scattering (ISLS) conducted on diamond-anvil cell (DAC) encapsulated samples offers an experimental approach to determine cross-pair potential interactions through equation of state determinations. In addition themore » kinetics of structural relaxation in fluids can be studied. We compare our experimental results with our thermochemical computational model Cheetah. Computational models are systematically improved with each addition of experimental data.« less

Kinetic study of the thermal hydrolysis of Agave salmiana for mezcal production.

PubMed

Garcia-Soto, M J; Jimenez-Islas, H; Navarrete-Bolanos, J L; Rico-Martinez, R; Miranda-Lopez, R; Botello-Alvarez, J E

2011-07-13

The kinetics of the thermal hydrolysis of the fructans of Agave salmiana were determined during the cooking step of mezcal production in a pilot autoclave. Thermal hydrolysis was achieved at different temperatures and cooking times, ranging from 96 to 116 °C and from 20 to 80 h. A simple kinetic model of the depolymerization of fructans to monomers and other reducing sugars and of the degradation of reducing sugars to furans [principally 5-(hydroxymethyl)furfural, HMF] was developed. From this model, the rate constants of the reactions were calculated, as well as the pre-exponential factors and activation energies of the Arrhenius equation. The model was found to fit the experimental data well. The tradeoff between a maximum fructan hydrolysis and a critical furan concentration in allowing for the best ethanol yield during fermentation was investigated. The results indicated that the thermal hydrolysis of agave was optimal, from the point of view of ethanol yield in the ensuing fermentation, in the temperature range of 106-116 °C and the cooking range time of 6-14 h. The optimal conditions corresponded to a fructan hydrolysis of 80%, producing syrups with furan and reducing sugar concentrations of 1 ± 0.1 and 110 ± 10 g/L, respectively.

First Passage Times, Lifetimes, and Relaxation Times of Unfolded Proteins

NASA Astrophysics Data System (ADS)

Dai, Wei; Sengupta, Anirvan M.; Levy, Ronald M.

2015-07-01

The dynamics of proteins in the unfolded state can be quantified in computer simulations by calculating a spectrum of relaxation times which describes the time scales over which the population fluctuations decay to equilibrium. If the unfolded state space is discretized, we can evaluate the relaxation time of each state. We derive a simple relation that shows the mean first passage time to any state is equal to the relaxation time of that state divided by the equilibrium population. This explains why mean first passage times from state to state within the unfolded ensemble can be very long but the energy landscape can still be smooth (minimally frustrated). In fact, when the folding kinetics is two-state, all of the unfolded state relaxation times within the unfolded free energy basin are faster than the folding time. This result supports the well-established funnel energy landscape picture and resolves an apparent contradiction between this model and the recently proposed kinetic hub model of protein folding. We validate these concepts by analyzing a Markov state model of the kinetics in the unfolded state and folding of the miniprotein NTL9 (where NTL9 is the N -terminal domain of the ribosomal protein L9), constructed from a 2.9 ms simulation provided by D. E. Shaw Research.

First Passage Times, Lifetimes, and Relaxation Times of Unfolded Proteins.

PubMed

Dai, Wei; Sengupta, Anirvan M; Levy, Ronald M

2015-07-24

The dynamics of proteins in the unfolded state can be quantified in computer simulations by calculating a spectrum of relaxation times which describes the time scales over which the population fluctuations decay to equilibrium. If the unfolded state space is discretized, we can evaluate the relaxation time of each state. We derive a simple relation that shows the mean first passage time to any state is equal to the relaxation time of that state divided by the equilibrium population. This explains why mean first passage times from state to state within the unfolded ensemble can be very long but the energy landscape can still be smooth (minimally frustrated). In fact, when the folding kinetics is two-state, all of the unfolded state relaxation times within the unfolded free energy basin are faster than the folding time. This result supports the well-established funnel energy landscape picture and resolves an apparent contradiction between this model and the recently proposed kinetic hub model of protein folding. We validate these concepts by analyzing a Markov state model of the kinetics in the unfolded state and folding of the miniprotein NTL9 (where NTL9 is the N-terminal domain of the ribosomal protein L9), constructed from a 2.9 ms simulation provided by D. E. Shaw Research.

A Statistical Physicist's Approach to Biological Motion: From the the Kinesin Walk to Muscle Contraction

NASA Astrophysics Data System (ADS)

Vicsek, Tamas

1997-03-01

It is demonstrated that a wide range of experimental results on biological motion can be successfully interpreted in terms of statistical physics motivated models taking into account the relevant microscopic details of motor proteins and allowing analytic solutions. Two important examples are considered, i) the motion of a single kinesin molecule along microtubules inside individual cells and ii) muscle contraction which is a macroscopic phenomenon due to the collective action of a large number of myosin heads along actin filaments. i) Recently individual two-headed kinesin molecules have been studied in in vitro motility assays revealing a number of their peculiar transport properties. Here we propose a simple and robust model for the kinesin stepping process with elastically coupled Brownian heads showing all of these properties. The analytic treatment of our model results in a very good fit to the experimental data and practically has no free parameters. ii) Myosin is an ATPase enzyme that converts the chemical energy stored in ATP molecules into mechanical work. During muscle contraction, the myosin cross-bridges attach to the actin filaments and exert force on them yielding a relative sliding of the actin and myosin filaments. In this paper we present a simple mechanochemical model for the cross-bridge interaction involving the relevant kinetic data and providing simple analytic solutions for the mechanical properties of muscle contraction, such as the force-velocity relationship or the relative number of the attached cross-bridges. So far the only analytic formula which could be fitted to the measured force-velocity curves has been the well known Hill equation containing parameters lacking clear microscopic origin. The main advantages of our new approach are that it explicitly connects the mechanical data with the kinetic data and the concentration of the ATP and ATPase products and as such it leads to new analytic solutions which agree extremely well with a wide range of experimental curves, while the parameters of the corresponding expressions have well defined microscopic meaning.

Modelling morphology evolution during solidification of IPP in processing conditions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pantani, R., E-mail: rpantani@unisa.it, E-mail: fedesantis@unisa.it, E-mail: vsperanza@unisa.it, E-mail: gtitomanlio@unisa.it; De Santis, F., E-mail: rpantani@unisa.it, E-mail: fedesantis@unisa.it, E-mail: vsperanza@unisa.it, E-mail: gtitomanlio@unisa.it; Speranza, V., E-mail: rpantani@unisa.it, E-mail: fedesantis@unisa.it, E-mail: vsperanza@unisa.it, E-mail: gtitomanlio@unisa.it

During polymer processing, crystallization takes place during or soon after flow. In most of cases, the flow field dramatically influences both the crystallization kinetics and the crystal morphology. On their turn, crystallinity and morphology affect product properties. Consequently, in the last decade, researchers tried to identify the main parameters determining crystallinity and morphology evolution during solidification In processing conditions. In this work, we present an approach to model flow-induced crystallization with the aim of predicting the morphology after processing. The approach is based on: interpretation of the FIC as the effect of molecular stretch on the thermodynamic crystallization temperature; modelingmore » the molecular stretch evolution by means of a model simple and easy to be implemented in polymer processing simulation codes; identification of the effect of flow on nucleation density and spherulites growth rate by means of simple experiments; determination of the condition under which fibers form instead of spherulites. Model predictions reproduce most of the features of final morphology observed in the samples after solidification.« less

On the derivation of a simple dynamic model of anaerobic digestion including the evolution of hydrogen.

PubMed

Giovannini, Giannina; Sbarciog, Mihaela; Steyer, Jean-Philippe; Chamy, Rolando; Vande Wouwer, Alain

2018-05-01

Hydrogen has been found to be an important intermediate during anaerobic digestion (AD) and a key variable for process monitoring as it gives valuable information about the stability of the reactor. However, simple dynamic models describing the evolution of hydrogen are not commonplace. In this work, such a dynamic model is derived using a systematic data driven-approach, which consists of a principal component analysis to deduce the dimension of the minimal reaction subspace explaining the data, followed by an identification of the kinetic parameters in the least-squares sense. The procedure requires the availability of informative data sets. When the available data does not fulfill this condition, the model can still be built from simulated data, obtained using a detailed model such as ADM1. This dynamic model could be exploited in monitoring and control applications after a re-identification of the parameters using actual process data. As an example, the model is used in the framework of a control strategy, and is also fitted to experimental data from raw industrial wine processing wastewater. Copyright © 2018 Elsevier Ltd. All rights reserved.

Autocatalytic polymerization generates persistent random walk of crawling cells.

PubMed

Sambeth, R; Baumgaertner, A

2001-05-28

The autocatalytic polymerization kinetics of the cytoskeletal actin network provides the basic mechanism for a persistent random walk of a crawling cell. It is shown that network remodeling by branching processes near the cell membrane is essential for the bimodal spatial stability of the network which induces a spontaneous breaking of isotropic cell motion. Details of the phenomena are analyzed using a simple polymerization model studied by analytical and simulation methods.

Unsteady Plasma Ejections from Hollow Accretion Columns of Galactic Neutron Stars as a Trigger for Gamma-Ray Bursts

NASA Astrophysics Data System (ADS)

Gvaramadze, V. V.

1995-09-01

We propose a model of gamma-ray bursts (GRBs) based on close Galactic neutron stars with accretion disks. We outline a simple mechanism of unsteady plasma ejections during episodic accretion events. The relative kinetic energy of ejected blobs can be converted into gamma-rays by internal shocks. The beaming of gamma-ray emission can be responsible for the observed isotropic angular distribution of GRBs.

A deterministic mathematical model for bidirectional excluded flow with Langmuir kinetics.

PubMed

Zarai, Yoram; Margaliot, Michael; Tuller, Tamir

2017-01-01

In many important cellular processes, including mRNA translation, gene transcription, phosphotransfer, and intracellular transport, biological "particles" move along some kind of "tracks". The motion of these particles can be modeled as a one-dimensional movement along an ordered sequence of sites. The biological particles (e.g., ribosomes or RNAPs) have volume and cannot surpass one another. In some cases, there is a preferred direction of movement along the track, but in general the movement may be bidirectional, and furthermore the particles may attach or detach from various regions along the tracks. We derive a new deterministic mathematical model for such transport phenomena that may be interpreted as a dynamic mean-field approximation of an important model from mechanical statistics called the asymmetric simple exclusion process (ASEP) with Langmuir kinetics. Using tools from the theory of monotone dynamical systems and contraction theory we show that the model admits a unique steady-state, and that every solution converges to this steady-state. Furthermore, we show that the model entrains (or phase locks) to periodic excitations in any of its forward, backward, attachment, or detachment rates. We demonstrate an application of this phenomenological transport model for analyzing ribosome drop off in mRNA translation.

Frontal Conversion and Uniformity in 3D Printing by Photopolymerisation

PubMed Central

Vitale, Alessandra; Cabral, João T.

2016-01-01

We investigate the impact of the non-uniform spatio-temporal conversion, intrinsic to photopolymerisation, in the context of light-driven 3D printing of polymers. The polymerisation kinetics of a series of model acrylate and thiol-ene systems, both neat and doped with a light-absorbing dye, is investigated experimentally and analysed according to a descriptive coarse-grained model for photopolymerisation. In particular, we focus on the relative kinetics of polymerisation with those of 3D printing, by comparing the evolution of the position of the conversion profile (zf) to the sequential displacement of the object stage (∆z). After quantifying the characteristic sigmoidal monomer-to-polymer conversion of the various systems, with a combination of patterning experiments, FT-IR mapping, and modelling, we compute representative regimes for which zf is smaller, commensurate with, or larger than ∆z. While non-monotonic conversion can be detrimental to 3D printing, for instance in causing differential shrinkage of inhomogeneity in material properties, we identify opportunities for facile fabrication of modulated materials in the z-direction (i.e., along the illuminated axis). Our simple framework and model, based on directly measured parameters, can thus be employed in photopolymerisation-based 3D printing, both in process optimisation and in the precise design of complex, internally stratified materials by coupling the z-stage displacement and frontal polymerisation kinetics. PMID:28773881

Removal of phosphate from aqueous solution using magnesium-alginate/chitosan modified biochar microspheres derived from Thalia dealbata.

PubMed

Cui, Xiaoqiang; Dai, Xi; Khan, Kiran Yasmin; Li, Tingqiang; Yang, Xiaoe; He, Zhenli

2016-10-01

The objective of this study was to determine the feasibility of using magnesium-alginate/chitosan modified biochar microspheres to enhance removal of phosphate from aqueous solution. The introduction of MgCl2 substantially increased surface area of biochar (116.2m(2)g(-1)), and both granulation with alginate/chitosan and modification with magnesium improved phosphate sorption on the biochars. Phosphate sorption on the biochars could be well described by a simple Langmuir model, and the MgCl2-alginate modified biochar microspheres exhibited the highest phosphate sorption capacity (up to 46.56mgg(-1)). The pseudo second order kinetic model better fitted the kinetic data, and both the Yoon-Nelson and Thomas models were superior to other models in describing phosphate dynamic sorption. Precipitation with minerals and ligand exchange were the possible mechanisms of phosphate sorption on the modified biochars. These results imply that MgCl2-alginate modified biochar microspheres have potential as a green cost-effective sorbent for remediating P contaminated water environment. Copyright © 2016 Elsevier Ltd. All rights reserved.

Acceleration of the direct reconstruction of linear parametric images using nested algorithms.

PubMed

Wang, Guobao; Qi, Jinyi

2010-03-07

Parametric imaging using dynamic positron emission tomography (PET) provides important information for biological research and clinical diagnosis. Indirect and direct methods have been developed for reconstructing linear parametric images from dynamic PET data. Indirect methods are relatively simple and easy to implement because the image reconstruction and kinetic modeling are performed in two separate steps. Direct methods estimate parametric images directly from raw PET data and are statistically more efficient. However, the convergence rate of direct algorithms can be slow due to the coupling between the reconstruction and kinetic modeling. Here we present two fast gradient-type algorithms for direct reconstruction of linear parametric images. The new algorithms decouple the reconstruction and linear parametric modeling at each iteration by employing the principle of optimization transfer. Convergence speed is accelerated by running more sub-iterations of linear parametric estimation because the computation cost of the linear parametric modeling is much less than that of the image reconstruction. Computer simulation studies demonstrated that the new algorithms converge much faster than the traditional expectation maximization (EM) and the preconditioned conjugate gradient algorithms for dynamic PET.

Kinetic rate constant prediction supports the conformational selection mechanism of protein binding.

PubMed

Moal, Iain H; Bates, Paul A

2012-01-01

The prediction of protein-protein kinetic rate constants provides a fundamental test of our understanding of molecular recognition, and will play an important role in the modeling of complex biological systems. In this paper, a feature selection and regression algorithm is applied to mine a large set of molecular descriptors and construct simple models for association and dissociation rate constants using empirical data. Using separate test data for validation, the predicted rate constants can be combined to calculate binding affinity with accuracy matching that of state of the art empirical free energy functions. The models show that the rate of association is linearly related to the proportion of unbound proteins in the bound conformational ensemble relative to the unbound conformational ensemble, indicating that the binding partners must adopt a geometry near to that of the bound prior to binding. Mirroring the conformational selection and population shift mechanism of protein binding, the models provide a strong separate line of evidence for the preponderance of this mechanism in protein-protein binding, complementing structural and theoretical studies.

Kinetics of SiHCl3 chemical vapor deposition and fluid dynamic simulations.

PubMed

Cavallotti, Carlo; Masi, Maurizio

2011-09-01

Though most of the current silicon photovoltaic technology relies on trichlorosilane (SiHCl3) as a precursor gas to deposit Si, only a few studies have been devoted to the investigation of its gas phase and surface kinetics. In the present work we propose a new kinetic mechanism apt to describe the gas phase and surface chemistry active during the deposition of Si from SiHCl3. Kinetic constants of key reactions were either taken from the literature or determined through ab initio calculations. The capability of the mechanism to reproduce experimental data was tested through the implementation of the kinetic scheme in a fluid dynamic model and in the simulation of both deposition and etching of Si in horizontal reactors. The results of the simulations show that the reactivity of HCl is of key importance in order to control the Si deposition rate. When HCl reaches a critical concentration in the gas phase it starts etching the Si surface, so that the net deposition rate is the net sum of the adsorption rate of the gas phase precursors and the etching rate due to HCl. In these conditions the possibility to further deposit Si is directly related to the rate of consumption of HCl through its reaction with SiHCl3 to give SiCl4. The proposed reaction mechanism was implemented in a 3D fluid dynamic model of a simple Siemens reactor. The simulation results indicate that the proposed interpretation of the growth process applies also to this class of reactors, which operate in what can be defined as a mixed kinetic-transport controlled regime.

Noninvasive in situ observation of the crystallization kinetics of biological macromolecules by confocal laser scanning microscopy.

PubMed

Mühlig, P; Klupsch, Th; Kaulmann, U; Hilgenfeld, R

2003-04-01

High-resolution confocal laser scanning microscopy (CLSM) is a powerful tool for in situ observation and analysis of protein crystal growth kinetics. Because the resolution of CLSM is not diffraction-limited by the object, it is possible to visualize, under certain conditions, objects in molecular dimensions. A modified batch technique is applied which allows the growth kinetics of sufficiently small crystallites fixed at the lower side of a cover glass, within a hanging drop, to be studied in reflected light near the total reflection angle. A gap, or cavity, filled with solution is formed between the cover glass and the upper crystal face, which acts to fix small crystallites by hydrodynamic friction forces. The cavity height enables the propagation of molecular steps across the upper crystal face without constraint, so that the propagation velocity and geometrical parameters can be measured by CLSM. The layer growth kinetics of monoclinic crystallites of a long-acting insulin derivative (Insulin Glargine) is investigated. For a twofold supersaturation of the solution, the growth is governed by 2D nucleation at the edges of the crystallites followed by a spreading of molecular steps. The layer growth kinetics are well fitted by the simple cubic kinetic lattice model. We find that only about one of a thousand solute (protein) molecules which push a kink place due to their Brownian motion becomes really incorporated into the growing crystal.

Supersymmetric quantum mechanics method for the Fokker-Planck equation with applications to protein folding dynamics

NASA Astrophysics Data System (ADS)

Polotto, Franciele; Drigo Filho, Elso; Chahine, Jorge; Oliveira, Ronaldo Junio de

2018-03-01

This work developed analytical methods to explore the kinetics of the time-dependent probability distributions over thermodynamic free energy profiles of protein folding and compared the results with simulation. The Fokker-Planck equation is mapped onto a Schrödinger-type equation due to the well-known solutions of the latter. Through a semi-analytical description, the supersymmetric quantum mechanics formalism is invoked and the time-dependent probability distributions are obtained with numerical calculations by using the variational method. A coarse-grained structure-based model of the two-state protein Tm CSP was simulated at a Cα level of resolution and the thermodynamics and kinetics were fully characterized. Analytical solutions from non-equilibrium conditions were obtained with the simulated double-well free energy potential and kinetic folding times were calculated. It was found that analytical folding time as a function of temperature agrees, quantitatively, with simulations and experiments from the literature of Tm CSP having the well-known 'U' shape of the Chevron Plots. The simple analytical model developed in this study has a potential to be used by theoreticians and experimentalists willing to explore, quantitatively, rates and the kinetic behavior of their system by informing the thermally activated barrier. The theory developed describes a stochastic process and, therefore, can be applied to a variety of biological as well as condensed-phase two-state systems.

Kinetic theory of shear thickening for a moderately dense gas-solid suspension: From discontinuous thickening to continuous thickening

NASA Astrophysics Data System (ADS)

Hayakawa, Hisao; Takada, Satoshi; Garzó, Vicente

2017-10-01

The Enskog kinetic theory for moderately dense gas-solid suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the environmental fluid on solid particles is modeled via a viscous drag force plus a stochastic Langevin-like term. The Enskog equation is solved by means of two independent but complementary routes: (i) Grad's moment method and (ii) event-driven Langevin simulation of hard spheres. Both approaches clearly show that the flow curve (stress-strain rate relation) depends significantly on the volume fraction of the solid particles. In particular, as the density increases, there is a transition from the discontinuous shear thickening (observed in dilute gases) to the continuous shear thickening for denser systems. The comparison between theory and simulations indicates that while the theoretical predictions for the kinetic temperature agree well with simulations for densities φ ≲0.5 , the agreement for the other rheological quantities (the viscosity, the stress ratio, and the normal stress differences) is limited to more moderate densities (φ ≲0.3 ) if the inelasticity during collisions between particles is not large.

Discrimination of ionic species from broad-beam ion sources

NASA Technical Reports Server (NTRS)

Anderson, J. R.

1993-01-01

The performance of a broad-beam, three-grid, ion extraction system incorporating radio frequency (RF) mass discrimination was investigated experimentally. This testing demonstrated that the system, based on a modified single-stage Bennett mass spectrometer, can discriminate between ionic species having about a 2-to-1 mass ratio while producing a broad-beam of ions with low kinetic energy (less than 15 eV). Testing was conducted using either argon and krypton ions or atomic and diatomic oxygen ions. A simple one-dimensional model, which ignores magnetic field and space-charge effects, was developed to predict the species separation capabilities as well as the kinetic energies of the extracted ions. The experimental results correlated well with the model predictions. This RF mass discrimination system can be used in applications where both atomic and diatomic ions are produced, but a beam of only one of the species is desired. An example of such an application is a 5 eV atomic oxygen source. This source would produce a beam of atomic oxygen with 5 eV kinetic energy, which would be directed onto a material specimen, to simulate the interaction between the surface of a satellite and the rarefied atmosphere encountered in low-Earth orbit.

Kinetic theory of shear thickening for a moderately dense gas-solid suspension: From discontinuous thickening to continuous thickening.

PubMed

Hayakawa, Hisao; Takada, Satoshi; Garzó, Vicente

2017-10-01

The Enskog kinetic theory for moderately dense gas-solid suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the environmental fluid on solid particles is modeled via a viscous drag force plus a stochastic Langevin-like term. The Enskog equation is solved by means of two independent but complementary routes: (i) Grad's moment method and (ii) event-driven Langevin simulation of hard spheres. Both approaches clearly show that the flow curve (stress-strain rate relation) depends significantly on the volume fraction of the solid particles. In particular, as the density increases, there is a transition from the discontinuous shear thickening (observed in dilute gases) to the continuous shear thickening for denser systems. The comparison between theory and simulations indicates that while the theoretical predictions for the kinetic temperature agree well with simulations for densities φ≲0.5, the agreement for the other rheological quantities (the viscosity, the stress ratio, and the normal stress differences) is limited to more moderate densities (φ≲0.3) if the inelasticity during collisions between particles is not large.

Understanding the effects of strain on morphological instabilities of a nanoscale island during heteroepitaxial growth

DOE Office of Scientific and Technical Information (OSTI.GOV)

Feng, Lu; Wang, Jing; Wang, Shibin

A comprehensive morphological stability analysis of a nanoscale circular island during heteroepitaxial growth is presented based on continuum elasticity theory. The interplay between kinetic and thermodynamic mechanisms is revealed by including strain-related kinetic processes. In the kinetic regime, the Burton-Cabrera-Frank model is adopted to describe the growth front of the island. Together with kinetic boundary conditions, various kinetic processes including deposition flow, adatom diffusion, attachment-detachment kinetics, and the Ehrlich-Schwoebel barrier can be taken into account at the same time. In the thermodynamic regime, line tension, surface energy, and elastic energy are considered. As the strain relief in the early stagesmore » of heteroepitaxy is more complicated than commonly suggested by simple consideration of lattice mismatch, we also investigate the effects of external applied strain and elastic response due to perturbations on the island shape evolution. The analytical expressions for elastic fields induced by mismatch strain, external applied strain, and relaxation strain are presented. A systematic approach is developed to solve the system via a perturbation analysis which yields the conditions of film morphological instabilities. Consistent with previous experimental and theoretical work, parametric studies show the kinetic evolution of elastic relaxation, island morphology, and film composition under various conditions. Our present work offers an effective theoretical approach to get a comprehensive understanding of the interplay between different growth mechanisms and how to tailor the growth mode by controlling the nature of the crucial factors.« less

Measuring the Acoustic Release of a Chemotherapeutic Agent from Folate-Targeted Polymeric Micelles.

PubMed

Abusara, Ayah; Abdel-Hafez, Mamoun; Husseini, Ghaleb

2018-08-01

In this paper, we compare the use of Bayesian filters for the estimation of release and re-encapsulation rates of a chemotherapeutic agent (namely Doxorubicin) from nanocarriers in an acoustically activated drug release system. The study is implemented using an advanced kinetic model that takes into account cavitation events causing the antineoplastic agent's release from polymeric micelles upon exposure to ultrasound. This model is an improvement over the previous representations of acoustic release that used simple zero-, first- and second-order release and re-encapsulation kinetics to study acoustically triggered drug release from polymeric micelles. The new model incorporates drug release and micellar reassembly events caused by cavitation allowing for the controlled release of chemotherapeutics specially and temporally. Different Bayesian estimators are tested for this purpose including Kalman filters (KF), Extended Kalman filters (EKF), Particle filters (PF), and multi-model KF and EKF. Simulated and experimental results are used to verify the performance of the above-mentioned estimators. The proposed methods demonstrate the utility and high-accuracy of using estimation methods in modeling this drug delivery technique. The results show that, in both cases (linear and non-linear dynamics), the modeling errors are expensive but can be minimized using a multi-model approach. In addition, particle filters are more flexible filters that perform reasonably well compared to the other two filters. The study improved the accuracy of the kinetic models used to capture acoustically activated drug release from polymeric micelles, which may in turn help in designing hardware and software capable of precisely controlling the delivered amount of chemotherapeutics to cancerous tissue.

Direct Monte Carlo simulation of chemical reaction systems: Simple bimolecular reactions

NASA Astrophysics Data System (ADS)

Piersall, Shannon D.; Anderson, James B.

1991-07-01

In applications to several simple reaction systems we have explored a ``direct simulation'' method for predicting and understanding the behavior of gas phase chemical reaction systems. This Monte Carlo method, originated by Bird, has been found remarkably successful in treating a number of difficult problems in rarefied dynamics. Extension to chemical reactions offers a powerful tool for treating reaction systems with nonthermal distributions, with coupled gas-dynamic and reaction effects, with emission and adsorption of radiation, and with many other effects difficult to treat in any other way. The usual differential equations of chemical kinetics are eliminated. For a bimolecular reaction of the type A+B→C+D with a rate sufficiently low to allow a continued thermal equilibrium of reactants we find that direct simulation reproduces the expected second order kinetics. Simulations for a range of temperatures yield the activation energies expected for the reaction models specified. For faster reactions under conditions leading to a depletion of energetic reactant species, the expected slowing of reaction rates and departures from equilibrium distributions are observed. The minimum sample sizes required for adequate simulations are as low as 1000 molecules for these cases. The calculations are found to be simple and straightforward for the homogeneous systems considered. Although computation requirements may be excessively high for very slow reactions, they are reasonably low for fast reactions, for which nonequilibrium effects are most important.

Nonlinear behavior during NO2 hydrogenation on a nanosized Pt-Rh catalyst sample

NASA Astrophysics Data System (ADS)

Barroo, Cédric; De Decker, Yannick; Jacobs, Luc; de Bocarmé, Thierry Visart

2017-08-01

Automotive pollution control crucially relies on the reactivity of metal alloy catalysts. Understanding how the chemistry of an alloy compares with that of pure metals forms a decisive step towards the rational development of applied formulations of such catalysts. In this context, we studied the hydrogenation of NO2 on Pt-Rh catalysts at the nanoscale with field emission microscopy (FEM). Previous studies have shown the presence of complex reaction kinetics at the surface of Pt for this reaction, including periodic oscillations at 390 K. As we briefly show here, similar kinetics can also be observed on Rh at higher temperatures. The alloy samples (Pt-17.4 at.%Rh) show signs of important reactivity and associated nonlinear dynamics in an intermediate temperature range. In particular, at 425 K isothermal oscillations are observed on this specific alloy catalyst. The role of the alloy composition on the window of reactivity is explained with a simple theoretical model for the kinetics of the reaction.

Optimization of Malachite Green Removal from Water by TiO₂ Nanoparticles under UV Irradiation.

PubMed

Ma, Yongmei; Ni, Maofei; Li, Siyue

2018-06-13

TiO₂ nanoparticles with surface porosity were prepared by a simple and efficient method and presented for the removal of malachite green (MG), a representative organic pollutant, from aqueous solution. Photocatalytic degradation experiments were systematically conducted to investigate the influence of TiO₂ dosage, pH value, and initial concentrations of MG. The kinetics of the reaction were monitored via UV spectroscopy and the kinetic process can be well predicted by the pseudo first-order model. The rate constants of the reaction kinetics were found to decrease as the initial MG concentration increased; increased via elevated pH value at a certain amount of TiO₂ dosage. The maximum efficiency of photocatalytic degradation was obtained when the TiO₂ dosage, pH value and initial concentrations of MG were 0.6 g/L, 8 and 10 −5 mol/L (M), respectively. Results from this study provide a novel optimization and an efficient strategy for water pollutant treatment.

α-Phase transformation kinetics of U – 8 wt% Mo established by in situ neutron diffraction

DOE PAGES

Garlea, Elena; Steiner, M. A.; Calhoun, C. A.; ...

2016-05-08

The α-phase transformation kinetics of as-cast U - 8 wt% Mo below the eutectoid temperature have been established by in situ neutron diffraction. α-phase weight fraction data acquired through Rietveld refinement at five different isothermal hold temperatures can be modeled accurately utilizing a simple Johnson-Mehl-Avrami-Kolmogorov impingement-based theory, and the results are validated by a corresponding evolution in the γ-phase lattice parameter during transformation that follows Vegard’s law. Neutron diffraction data is used to produce a detailed Time-Temperature-Transformation diagram that improves upon inconsistencies in the current literature, exhibiting a minimum transformation start time of 40 min at temperatures between 500 °Cmore » and 510 °C. Lastly, the transformation kinetics of U – 8 wt% Mo can vary significantly from as-cast conditions after extensive heat treatments, due to homogenization of the typical dendritic microstructure which possesses non-negligible solute segregation.« less

Tension-dependent structural deformation alters single-molecule transition kinetics.

PubMed

Sudhanshu, B; Mihardja, S; Koslover, E F; Mehraeen, S; Bustamante, C; Spakowitz, A J

2011-02-01

We analyze the response of a single nucleosome to tension, which serves as a prototypical biophysical measurement where tension-dependent deformation alters transition kinetics. We develop a statistical-mechanics model of a nucleosome as a wormlike chain bound to a spool, incorporating fluctuations in the number of bases bound, the spool orientation, and the conformations of the unbound polymer segments. With the resulting free-energy surface, we perform dynamic simulations that permit a direct comparison with experiments. This simple approach demonstrates that the experimentally observed structural states at nonzero tension are a consequence of the tension and that these tension-induced states cease to exist at zero tension. The transitions between states exhibit substantial deformation of the unbound polymer segments. The associated deformation energy increases with tension; thus, the application of tension alters the kinetics due to tension-induced deformation of the transition states. This mechanism would arise in any system where the tether molecule is deformed in the transition state under the influence of tension.

Tension-dependent structural deformation alters single-molecule transition kinetics

PubMed Central

Sudhanshu, B.; Mihardja, S.; Koslover, E. F.; Mehraeen, S.; Bustamante, C.; Spakowitz, A. J.

2011-01-01

We analyze the response of a single nucleosome to tension, which serves as a prototypical biophysical measurement where tension-dependent deformation alters transition kinetics. We develop a statistical-mechanics model of a nucleosome as a wormlike chain bound to a spool, incorporating fluctuations in the number of bases bound, the spool orientation, and the conformations of the unbound polymer segments. With the resulting free-energy surface, we perform dynamic simulations that permit a direct comparison with experiments. This simple approach demonstrates that the experimentally observed structural states at nonzero tension are a consequence of the tension and that these tension-induced states cease to exist at zero tension. The transitions between states exhibit substantial deformation of the unbound polymer segments. The associated deformation energy increases with tension; thus, the application of tension alters the kinetics due to tension-induced deformation of the transition states. This mechanism would arise in any system where the tether molecule is deformed in the transition state under the influence of tension. PMID:21245354

Sunflower shells utilization for energetic purposes in an integrated approach of energy crops: laboratory study pyrolysis and kinetics.

PubMed

Zabaniotou, A A; Kantarelis, E K; Theodoropoulos, D C

2008-05-01

Sunflower is a traditional crop which can be used for the production of bioenergy and liquid biofuels. A study of the pyrolytic behaviour of sunflower residues at temperatures from 300 to 600 degrees C has been carried out. The experiments were performed in a captive sample reactor under atmospheric pressure and helium as sweeping gas. The yields of the derived pyrolysis products were determined in relation to temperature, with constant sweeping gas flow of 50 cm3 min(-1) and heating rate of 40 degrees Cs(-1). The maximum gas yield of around 53 wt.% was obtained at 500 degrees C, whereas maximum oil yield of about 21 wt.% was obtained at 400 degrees C. A simple first order kinetic model has been applied for the devolatilization of biomass. Kinetic constants have been estimated: E=78.15 kJ mol(-1); k(0)=1.03 x 10(3)s(-1).

NIMROD modeling of poloidal flow damping in tokamaks using kinetic closures

NASA Astrophysics Data System (ADS)

Jepson, J. R.; Hegna, C. C.; Held, E. D.

2017-10-01

Calculations of poloidal flow damping in a tokamak are undertaken using two different implementations of the ion drift kinetic equation (DKE) in the extended MHD code NIMROD. The first approach is hybrid fluid/kinetic and uses a Chapman Enskog-like (CEL) Ansatz. Closure of the evolving lower-order fluid moment equations for n, V , and T is provided by solutions to the ion CEL-DKE written in the macroscopic flow reference frame. The second implementation solves the DKE using a delta-f approach. Here, the delta-f distribution describes all of the information beyond a static, lowest-order Maxwellian. We compare the efficiency and accuracy of these two approaches for a simple initial value problem that monitors the relaxation of the poloidal flow profile in high- and low-aspect-ratio tokamak geometry. The computation results are compared against analytic predictions of time dependent closures for the parallel viscous force. Supported by DoE Grants DE-FG02-86ER53218 and DE-FG02-04ER54746.

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

DOE PAGES

Tang, Jin-Yun; Riley, William J.

2017-09-05

Several land biogeochemical models used for studying carbon–climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate–consumer relationships for consumer-mediated redox reactions of the form A +  B E→  products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate–consumer relationships can be formulated as the computationally difficult full equilibrium chemistrymore » problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation – ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate–consumer interactions and can be applied to improve Earth system model (ESM) land models.« less

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tang, Jin-Yun; Riley, William J.

Several land biogeochemical models used for studying carbon–climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate–consumer relationships for consumer-mediated redox reactions of the form A +  B E→  products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate–consumer relationships can be formulated as the computationally difficult full equilibrium chemistrymore » problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation – ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate–consumer interactions and can be applied to improve Earth system model (ESM) land models.« less

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

NASA Astrophysics Data System (ADS)

Tang, Jin-Yun; Riley, William J.

2017-09-01

Several land biogeochemical models used for studying carbon-climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate-consumer relationships for consumer-mediated redox reactions of the form A + BE→ products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate-consumer relationships can be formulated as the computationally difficult full equilibrium chemistry problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation - ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate-consumer interactions and can be applied to improve Earth system model (ESM) land models.