Monte Carlo investigation of the critical behavior of Stavskaya's probabilistic cellular automaton

Stavskaya's model is a one-dimensional probabilistic cellular automaton (PCA) introduced in the end of the 1960s as an example of a model displaying a nonequilibrium phase transition. Although its absorbing state phase transition is well understood nowadays, the model never received a full numerical treatment to investigate its critical behavior. In this Brief Report we characterize the critical behavior of Stavskaya's PCA by means of Monte Carlo simulations and finite-size scaling analysis. The critical exponents of the model are calculated and indicate that its phase transition belongs to the directed percolation universality class of critical behavior, as would be expected on the basis of the directed percolation conjecture. We also explicitly establish the relationship of the model with the Domany-Kinzel PCA on its directed site percolation line, a connection that seems to have gone unnoticed in the literature so far.

Comparing performance of instrumental drift correction by linear and quadratic adjusting in inductive electromagnetic data

Electromagnetic induction (EMI) method results are shown for vertical magnetic dipole (VMD) configuration by using the EM38 equipment. Performance in the location of metallic pipes and electrical cables is compared as a function of instrumental drift correction by linear and quadratic adjusting under controlled conditions. Metallic pipes and electrical cables are buried at the IAG/USP shallow geophysical test site in Sao Paulo City. Brazil. Results show that apparent electrical conductivity and magnetic susceptibility data were affected by ambient temperature variation. In order to obtain better contrast between background and metallic targets it was necessary to correct the drift. This correction was accomplished by using linear and quadratic relation between conductivity/susceptibility and temperature intending comparative studies. The correction of temperature drift by using a quadratic relation was effective, showing that all metallic targets were located as well deeper targets were also improved. (C) 2010 Elsevier B.V. All rights reserved.

An extinction-survival-type phase transition in the probabilistic cellular automaton p182-q200

We investigate the critical behaviour of a probabilistic mixture of cellular automata (CA) rules 182 and 200 (in Wolfram`s enumeration scheme) by mean-field analysis and Monte Carlo simulations. We found that as we switch off one CA and switch on the other by the variation of the single parameter of the model, the probabilistic CA (PCA) goes through an extinction-survival-type phase transition, and the numerical data indicate that it belongs to the directed percolation universality class of critical behaviour. The PCA displays a characteristic stationary density profile and a slow, diffusive dynamics close to the pure CA 200 point that we discuss briefly. Remarks on an interesting related stochastic lattice gas are addressed in the conclusions.

Cellular automaton model for molecular traffic jams

We consider the time evolution of an exactly solvable cellular automaton with random initial conditions both in the large-scale hydrodynamic limit and on the microscopic level. This model is a version of the totally asymmetric simple exclusion process with sublattice parallel update and thus may serve as a model for studying traffic jams in systems of self-driven particles. We study the emergence of shocks from the microscopic dynamics of the model. In particular, we introduce shock measures whose time evolution we can compute explicitly, both in the thermodynamic limit and for open boundaries where a boundary-induced phase transition driven by the motion of a shock occurs. The motion of the shock, which results from the collective dynamics of the exclusion particles, is a random walk with an internal degree of freedom that determines the jump direction. This type of hopping dynamics is reminiscent of some transport phenomena in biological systems.

Alveolar osseous defect in rat for cell therapy: preliminary report

PURPOSE: To study were to reproduce an alveolar bone defect model in Wistar rats to be used for testing the efficacy of stem cell therapies. Additionally, we also aimed to determine the osteogenesis process of this osseous defect in the 1 month period post-surgery. METHODS: The animals were randomly divided into two groups of 7 animals each. A gingivobuccal incision was made, and a bone defect of 28 mm² of area was performed in the alveolar region. Animals were killed at 2 weeks after surgery (n=7) and 4 weeks after surgery (n=7). RESULTS: The average area of the alveolar defect at time point of 2 weeks was 22.27 ± 1.31 mm² and the average area of alveolar defect at time point of 4 weeks was 9.03 ± 1.17 mm². The average amount of bone formation at time point of 2 weeks was 5.73 ± 1.31 mm² and the average amount of bone formation at time point of 4 weeks was 19 ± 1.17 mm². Statistically significant differences between the amount of bone formation at 2 weeks and 4 weeks after surgery were seen (p=0.003). CONCLUSION: The highest rate of ossification occurred mostly from 2 to 4 weeks after surgery. This observation suggests that 4 weeks after the bone defect creation should be a satisfactory timing to assess the potential of bone inductive stem cells to accelerate bone regeneration in Wistar rats.

Avaliação da memória de trabalho em crianças com transtorno do déficit de atenção e hiperatividade

Esse projeto investigou as habilidades cognitivas da Memória de Trabalho de crianças e adolescentes com o “Transtorno do Déficit de Atenção e Hiperatividade” (TDAH). Foi aplicado teste TIHC - Teste Infantil de Habilidades Cognitivas informatizado, avaliando principalmente cinco habilidades cognitivas Raciocínio Indutivas (RI), Memória de Armazenamento Auditiva (MAA) e Visual e Memória de Trabalho Auditiva (MTA) e Visual (MTV) usando o modelo de Cattell-Horn-Carroll. O teste foi aplicado em 62 sujeitos sendo 32 com TDAH do Hospital das Clínicas-SEPIA-USP e 30 sujeitos controle “normal” do Ensino Fundamental da rede municipal de 1ª à 8ª séries com idades entre 7 e 15 anos, de ambos os sexos. As análises executadas indicaram a existência de diferenças em diversos aspectos relacionados à memória nos tipos peculiares de crianças com TDAH. Verificou-se que as crianças com TDAH obtiveram bom desempenho nas provas de memória visual em detrimento as provas de memória auditiva. A correlação dos resultados do TIHC verificara que o tempo de reação da prova de memória de armazenamento visual manteve uma correlação significativa com a prova de cálculo.

The Effect of Spatial Scale on Predicting Time Series: A Study on Epidemiological System Identification

A susceptible-infective-recovered (SIR) epidemiological model based on probabilistic cellular automaton (PCA) is employed for simulating the temporal evolution of the registered cases of chickenpox in Arizona, USA, between 1994 and 2004. At each time step, every individual is in one of the states S, I, or R. The parameters of this model are the probabilities of each individual (each cell forming the PCA lattice ) passing from a state to another state. Here, the values of these probabilities are identified by using a genetic algorithm. If nonrealistic values are allowed to the parameters, the predictions present better agreement with the historical series than if they are forced to present realistic values. A discussion about how the size of the PCA lattice affects the quality of the model predictions is presented. Copyright (C) 2009 L. H. A. Monteiro et al.

System Identification and Prediction of Dengue Fever Incidence in Rio de Janeiro

Identification, prediction, and control of a system are engineering subjects, regardless of the nature of the system. Here, the temporal evolution of the number of individuals with dengue fever weekly recorded in the city of Rio de Janeiro, Brazil, during 2007, is used to identify SIS (susceptible-infective-susceptible) and SIR (susceptible-infective-removed) models formulated in terms of cellular automaton (CA). In the identification process, a genetic algorithm (GA) is utilized to find the probabilities of the state transition S -> I able of reproducing in the CA lattice the historical series of 2007. These probabilities depend on the number of infective neighbors. Time-varying and non-time-varying probabilities, three different sizes of lattices, and two kinds of coupling topology among the cells are taken into consideration. Then, these epidemiological models built by combining CA and GA are employed for predicting the cases of sick persons in 2008. Such models can be useful for forecasting and controlling the spreading of this infectious disease.

Critical behavior of the susceptible-infected-recovered model on a square lattice

By means of numerical simulations and epidemic analysis, the transition point of the stochastic asynchronous susceptible-infected-recovered model on a square lattice is found to be c(0)=0.176 500 5(10), where c is the probability a chosen infected site spontaneously recovers rather than tries to infect one neighbor. This point corresponds to an infection/recovery rate of lambda(c)=(1-c(0))/c(0)=4.665 71(3) and a net transmissibility of (1-c(0))/(1+3c(0))=0.538 410(2), which falls between the rigorous bounds of the site and bond thresholds. The critical behavior of the model is consistent with the two-dimensional percolation universality class, but local growth probabilities differ from those of dynamic percolation cluster growth, as is demonstrated explicitly.

Role of noise in population dynamics cycles

Noise is an intrinsic feature of population dynamics and plays a crucial role in oscillations called phase-forgetting quasicycles by converting damped into sustained oscillations. This function of noise becomes evident when considering Langevin equations whose deterministic part yields only damped oscillations. We formulate here a consistent and systematic approach to population dynamics, leading to a Fokker-Planck equation and the associate Langevin equations in accordance with this conceptual framework, founded on stochastic lattice-gas models that describe spatially structured predator-prey systems. Langevin equations in the population densities and predator-prey pair density are derived in two stages. First, a birth-and-death stochastic process in the space of prey and predator numbers and predator-prey pair number is obtained by a contraction method that reduces the degrees of freedom. Second, a van Kampen expansion in the inverse of system size is then performed to get the Fokker-Planck equation. We also study the time correlation function, the asymptotic behavior of which is used to characterize the transition from the cyclic coexistence of species to the ordinary coexistence.

Time correlation function in systems with two coexisting biological species

We study a stochastic lattice model describing the dynamics of coexistence of two interacting biological species. The model comprehends the local processes of birth, death, and diffusion of individuals of each species and is grounded on interaction of the predator-prey type. The species coexistence can be of two types: With self-sustained coupled time oscillations of population densities and without oscillations. We perform numerical simulations of the model on a square lattice and analyze the temporal behavior of each species by computing the time correlation functions as well as the spectral densities. This analysis provides an appropriate characterization of the different types of coexistence. It is also used to examine linked population cycles in nature and in experiment.

Sensitivity to noise and ergodicity of an assembly line of cellular automata that classifies density

We investigate the sensitivity of the composite cellular automaton of H. Fuks [Phys. Rev. E 55, R2081 (1997)] to noise and assess the density classification performance of the resulting probabilistic cellular automaton (PCA) numerically. We conclude that the composite PCA performs the density classification task reliably only up to very small levels of noise. In particular, it cannot outperform the noisy Gacs-Kurdyumov-Levin automaton, an imperfect classifier, for any level of noise. While the original composite CA is nonergodic, analyses of relaxation times indicate that its noisy version is an ergodic automaton, with the relaxation times decaying algebraically over an extended range of parameters with an exponent very close (possibly equal) to the mean-field value.

Characterization and application of the solid waste generated in the production of calcium carbonate precipitate as a corrective for the soil acidity

The calcium carbonate industry generates solid waste products which, because of their high alkaline content (CaO, CaCO(3) and Ca (OH)(2)), have a substantial impact on the environment. The objectives of this study are to characterize and classify the solid waste products, which are generated during the hydration process of the calcium carbonate industry, according to ABNT`s NBR 10.000 series, and to determine the potential and efficiency of using these solid residues to correct soil acidity. Initially, the studied residue was submitted to gross mass, leaching, solubility, pH. X-ray Diffractometry, Inductive Coupled Plasma - Atomic Emission Spectrometry (ICP-AES), granularity and humidity analyses. The potential and efficiency of the residue for correcting soil acidity was determined by analysis of the quality attributes for soil correctives (PN, PRNT, Ca and Mg contents, granularity). Consequently, the results show that the studied residue may be used as a soil acidity corrective, considering that a typical corrective compound is recommended for each different type of soil. Additionally, the product must be further treated (dried and ground) to suit the specific requirements of the consumer market.

Embrittlement of case hardened steel chain link

Various steel chain links presented cracking during their manufacturing process, which includes induction case hardening and electrogalvanizing steps. Fractographic examination of the exposed crack surfaces revealed intergranular cracking with some areas featuring a thin layer of iron oxide, indicating that the cracking took place after the electrogalvanizing step. The location of the cracks coincided with the position of the deepest case hardened layer, suggesting the occurrence of localized overheating during the induction case hardening step. Inductive heating finite element analysis (COSMOS Designstar Software) confirmed that during the case hardening the austenitising temperature reached in the crack region values of approximately 1050 degrees C. The results indicated that intergranular cracking was caused by hydrogen embrittlement. (C) 2009 Elsevier Ltd. All rights reserved.

Recombinant rabies virus glycoprotein synthesis in bioreactor by transfected Drosophila melanogaster S2 cells carrying a constitutive or an inducible promoter

S2 cell populations (S2AcRVGP2K and S2MtRVGP-Hy) were selected after transfection of gene expression vectors carrying the cDNA encoding the rabies virus glycoprotein (RVGP) gene under the control of the constitutive (actin) or inductive (metallothionein) promoters. These cell populations were cultivated in a 1 L bioreactor mimicking a large scale bioprocess. Cell cultures were carried out at 90 rpm and monitored/controlled for temperature (28 degrees C) and dissolved oxygen (10 or 50% air saturation). Cell growth attained similar to 1.5-3 x 10(7) cells/mL after 3-4 clays of cultivation. The constitutive synthesis of RVGP in S2AcRVGP2K cells led to values of 0.76 mu g/10(7) cells at day 4 of culture. The RVGP synthesis in S2MtRVGP-Hy cell fraction increased upon CuSO(4) induction attaining specific productivities of 1.5-2 mu g/10(7) cells at clays 4-5. RVGP values in supernatant as a result of cell lysis were always very low (<0.2 mu g/mL) indicating good integrity of cells in culture. Overall the RVGP productivity was of 1.5-3 mg/L. Our data showed an important influence of dissolved oxygen on RVGP synthesis allowing a higher and sustained productivity by S2MtRVGP-Hy cells when cultivated with a DO of 10% air saturation. The RVGP productivity in bioreactors shown here mirrors those previously observed for T-flasks and shaker bottles and allow the preparation of the large RVGP quantities required for studies of structure and function. (C) 2010 Elsevier B.V. All rights reserved.