Coexistence of interacting opinions in a generalized Sznajd model

The Sznajd model is a sociophysics model that mimics the propagation of opinions in a closed society, where the interactions favor groups of agreeing people. It is based in the Ising and Potts ferromagnetic models and, although the original model used only linear chains, it has since been adapted to general networks. This model has a very rich transient, which has been used to model several aspects of elections, but its stationary states are always consensus states. In order to model more complex behaviors, we have, in a recent work, introduced the idea of biases and prejudices to the Sznajd model by generalizing the bounded confidence rule, which is common to many continuous opinion models, to what we called confidence rules. In that work we have found that the mean field version of this model (corresponding to a complete network) allows for stationary states where noninteracting opinions survive, but never for the coexistence of interacting opinions. In the present work, we provide networks that allow for the coexistence of interacting opinions for certain confidence rules. Moreover, we show that the model does not become inactive; that is, the opinions keep changing, even in the stationary regime. This is an important result in the context of understanding how a rule that breeds local conformity is still able to sustain global diversity while avoiding a frozen stationary state. We also provide results that give some insights on how this behavior approaches the mean field behavior as the networks are changed.

Power-law creep behavior of a semiflexible chain

Rheological properties of adherent cells are essential for their physiological functions, and microrheological measurements on living cells have shown that their viscoelastic responses follow a weak power law over a wide range of time scales. This power law is also influenced by mechanical prestress borne by the cytoskeleton, suggesting that cytoskeletal prestress determines the cell's viscoelasticity, but the biophysical origins of this behavior are largely unknown. We have recently developed a stochastic two-dimensional model of an elastically joined chain that links the power-law rheology to the prestress. Here we use a similar approach to study the creep response of a prestressed three-dimensional elastically jointed chain as a viscoelastic model of semiflexible polymers that comprise the prestressed cytoskeletal lattice. Using a Monte Carlo based algorithm, we show that numerical simulations of the chain's creep behavior closely correspond to the behavior observed experimentally in living cells. The power-law creep behavior results from a finite-speed propagation of free energy from the chain's end points toward the center of the chain in response to an externally applied stretching force. The property that links the power law to the prestress is the chain's stiffening with increasing prestress, which originates from entropic and enthalpic contributions. These results indicate that the essential features of cellular rheology can be explained by the viscoelastic behaviors of individual semiflexible polymers of the cytoskeleton.

Connectivity and dynamics of neuronal networks as defined by the shape of individual neurons

Biological neuronal networks constitute a special class of dynamical systems, as they are formed by individual geometrical components, namely the neurons. In the existing literature, relatively little attention has been given to the influence of neuron shape on the overall connectivity and dynamics of the emerging networks. The current work addresses this issue by considering simplified neuronal shapes consisting of circular regions (soma/axons) with spokes (dendrites). Networks are grown by placing these patterns randomly in the two-dimensional (2D) plane and establishing connections whenever a piece of dendrite falls inside an axon. Several topological and dynamical properties of the resulting graph are measured, including the degree distribution, clustering coefficients, symmetry of connections, size of the largest connected component, as well as three hierarchical measurements of the local topology. By varying the number of processes of the individual basic patterns, we can quantify relationships between the individual neuronal shape and the topological and dynamical features of the networks. Integrate-and-fire dynamics on these networks is also investigated with respect to transient activation from a source node, indicating that long-range connections play an important role in the propagation of avalanches.

Timoshenko versus Euler beam theory: Pitfalls of a deterministic approach

The selection criteria for Euler-Bernoulli or Timoshenko beam theories are generally given by means of some deterministic rule involving beam dimensions. The Euler-Bernoulli beam theory is used to model the behavior of flexure-dominated (or ""long"") beams. The Timoshenko theory applies for shear-dominated (or ""short"") beams. In the mid-length range, both theories should be equivalent, and some agreement between them would be expected. Indeed, it is shown in the paper that, for some mid-length beams, the deterministic displacement responses for the two theories agrees very well. However, the article points out that the behavior of the two beam models is radically different in terms of uncertainty propagation. In the paper, some beam parameters are modeled as parameterized stochastic processes. The two formulations are implemented and solved via a Monte Carlo-Galerkin scheme. It is shown that, for uncertain elasticity modulus, propagation of uncertainty to the displacement response is much larger for Timoshenko beams than for Euler-Bernoulli beams. On the other hand, propagation of the uncertainty for random beam height is much larger for Euler beam displacements. Hence, any reliability or risk analysis becomes completely dependent on the beam theory employed. The authors believe this is not widely acknowledged by the structural safety or stochastic mechanics communities. (C) 2010 Elsevier Ltd. All rights reserved.

A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem.

Effect of bituminous impregnation on nonwoven geotextiles tensile and permeability properties

Geosynthetics interlayer systems are effective techniques to control reflective cracking in damaged pavements. It comprises the inclusion of nonwoven geotextiles between the damaged layer and the new overlay of the pavement to reduce the propagation of cracks and to extend pavement life. However, the success of this technique depends directly on the understanding of the geotextile`s behavior when impregnated with asphalt This paper evaluates different nonwoven geotextiles frequently used in anti-reflective cracking systems, focusing on initial stiffness gain and permeability reduction after asphalt impregnation. Fresh and impregnated samples of polyester and polypropylene nonwoven geotextiles were tested. Cationic rapid setting emulsified asphalt was used as asphalt binder. Wide-width tensile tests were carried out based on the specification of ABNT - NBR 12824 (1993). Water vapor transmission tests were conducted according to ASTM E 96M (2005). Results of tensile tests on impregnated geotextiles showed a significant increase on tensile strength values, probably due to the inter contact of the fibers. Results also showed high increase in strength values at strain levels less than 0.05% and decrease on stiffness gains with increase of strains. Water vapor transmission tests demonstrated that cationic asphalt emulsion applied on nonwoven geotextiles allows a drastic reduction in permeability values to turn nonwoven geotextiles into a low permeability barrier. (C) 2010 Elsevier Ltd. All rights reserved.

Influence of Microstructures on Polarization Resistance and Hydrogen-Induced Cracking of a Micro-Alloyed Steel Submitted to Different Cooling Rates

The effect of different microstructures on the polarization resistance (Rp) and the hydrogen-induced cracking (HIC) of a micro-alloyed steel austenitized and submitted to different cooling rates was studied. Samples 19.1 x 6 x 2 mm, containing the whole thickness of the plate were extracted from a 20 mm plate and heat treated on a quenching dilatometer, were submitted to Rp and HIC corrosion tests. Both Rp and HIC tests followed as close as possible ASTM G59 and NACE standard TM0284-2003, in this case, modified only with regard to the size of the samples. Steel samples transformed from austenite by a slow cooling (cooling rate of 0.5 degrees C.s(-1)) showed higher susceptibility to hydrogen-induced cracking, with large cracks in the middle of the sample propagating along segregation bands, corresponding to the centerline of the plate thickness. For cooling rates of 10 degrees C.s(-1), only small cracks were found in the matrix and micro cracks nucleated at non-metallic inclusions. For higher cooling rates (40 degrees C.s(-1)) very few small cracks were detected, linked to non-metallic inclusions. This result suggests that structures formed by polygonal structures and segregation bands (were cutectoid microconstituents predominate) have higher susceptibility to HIC. Structures predominantly formed by acicular ferrite make it difficult to propagate the cracks among non-oriented and interlaced acicular ferrite crystals. Smaller segregation bands containing eutectoid products also help inhibit cracking and crack propagation; segregation bands can function as pipelines for hydrogen diffusion and offer a path of stress concentration for the propagation of cracks, frequently associated to non-metallic inclusions. Polarization resistance essays performed on the steel in theas received condition, prior to any heat treatment, showed larger differences between the regions of the plate, with a considerably lower Rp in the centerline. The austenitization heat treatments followed by cooling rates of 0.5 e 10 degrees C.s(-1) made more uniform the corrosion resistance along the thickness of the plate. The effects of heat treatments on the corrosion resistance are probably related to the microconstituent formed, allied to the chemical homogenization of the impurities concentrated on the centerline of the plate.

Who should wear mask against airborne infections? Altering the contact network for controlling the spread of contagious diseases

There are several ways of controlling the propagation of a contagious disease. For instance, to reduce the spreading of an airborne infection, individuals can be encouraged to remain in their homes and/or to wear face masks outside their domiciles. However, when a limited amount of masks is available, who should use them: the susceptible subjects, the infective persons or both populations? Here we employ susceptible-infective-recovered (SIR) models described in terms of ordinary differential equations and probabilistic cellular automata in order to investigate how the deletion of links in the random complex network representing the social contacts among individuals affects the dynamics of a contagious disease. The inspiration for this study comes from recent discussions about the impact of measures usually recommended by health public organizations for preventing the propagation of the swine influenza A (H1N1) virus. Our answer to this question can be valid for other eco-epidemiological systems. (C) 2010 Elsevier BM. All rights reserved.

ESEM observations of SCC initiation for 4340 high strength steel in distilled water

The stress corrosion cracking (SCC) initiation process for 4340 high strength steel in distilled water at room temperature was studied using a new kind of instrument: an environmental scanning electron microscope (ESEM). It was found that the applied stress accelerated oxide film formation which has an important influence on the subsequent SCC initiation. SCC was observed to initiate in the following circumstances: (1) cracking of a thick oxide film leading to SCC initiation along metal grain boundaries, (2) the initiation of pits initiating SCC in the metal and (3) SCC initiating from the edge of the specimen. All these three SCC initiation circumstances are consistent with the following model which couples SCC initiation with cracking of a surface protective oxide. There is a dynamic interaction between oxide formation, the applied stress, oxide cracking, pitting and the initiation of SCC. An aspect of the dynamic interaction is cracks forming in a protective surface oxide because of the applied stress, exposing to the water bare metal at the oxide crack tip, and oxidation of the bare metal causing crack healing. Oxide crack healing would be competing with the initiation of intergranular SCC if an oxide crack meets the metal surface at a grain boundary. If the intergranular SCC penetration is sufficiently fast along the metal grain boundary, then the crack yaws open preventing healing of the oxide crack. If intergranular SCC penetration is not sufficiently fast, then the oxidation process could produce sufficient oxide to fill both the stress corrosion crack and the oxide crack; in this case there would be initiation of SCC but only limited propagation of SCC. Stress-induced cracks in very thin oxide can induce pits which initiate SCC, and under some conditions such stress induced cracks in a thin oxide can directly initiate SCC.

The crack tip strain field of AISI 4340 - Part II - Experimental measurements

Crack tip strain maps have been measured for AISI 4340 high strength steel. No significant creep was observed. The measured values of CTOD were greater than expected from the HRR model. Crack tip branching was observed in every experiment. The direction of crack branching was in the same direction as a major ridge'' of epsilon(yy) strain, which in turn was in the same direction as predicted by the HRR model. Furthermore, the measured magnitudes of the epsilon(y)y strain in this same direction were in general greater than the values predicted by the HRR model. This indicates more plasticity in the crack tip region than expected from the HRR model. This greater plasticity could be related to the larger than expected CTOD values. The following discrepancies between the measured strain fields for AISI 4340 and the HRR predictions are noteworthy: (1) The crack branching. (2) Values of CTOD significantly higher than predicted by HRR. (3) The major ridge'' of epsilon(yy) strain an angle of about 60 degrees with the direction of overall propagation of the fatigue precrack, in which the measured magnitudes of the epsilon(yy) strain were greater than the values predicted by the HRR model. (4) Asymmetric shape of the plastic zone as measured by the epsilon(yy) strain. (5) Values of shear strain gamma(xy) significantly higher than predicted by the HRR model. (C) 1999 Kluwer Academic Publishers.

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary's entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429-35) for two-dimensional non-interacting tidal waves. (C) 2000 Elsevier Science Ltd. All rights reserved.

Apoptotic and developmental effects of bovine Herpesvirus type-5 infection on in vitro-produced bovine embryos

Bovine Herpesvirus type-5 (BoHV-5), which is potentially neuropathogenic, was recently described to be related with reproductive disorders in cows. The objective was to elucidate mechanisms involved in propagation of BoHV-5 in embryonic cells. For this purpose, bovine embryos produced in vitro were assayed for apoptotic markers after experimental infection of oocytes, in vitro fertilization, and development. Host DNA fragmentation was detected with a TUNEL assay, expression of annexin-V was measured with indirect immunofluorescence, and viral DNA was detected with in situ hybridization. Infective BoHV-5 virus was recovered from embryos derived from exposed oocytes after two consecutive passages on Madin-Darby bovine kidney (MDBK) cells. The viral DNA corresponding to US9 gene, localized between nucleotides 126243 to 126493, was detected in situ and amplified. There was no significant difference between the ratio of TUNEL stained nuclei and total cells in good quality blastocysts (0.87 +/- 0.05, mean SD), but there were differences (P < 0.05) between infected (0.18 +/- 0.05) and uninfected blastocysts (0.73 +/- 0.07). The Annexin-V label was more intense in uninfected embryos (0.79 +/- 0.04; P < 0.05). The quality of infected and uninfected embryos was considered equal, with no significant effect on embryonic development. In conclusion, we inferred that BoHV-5 infected bovine oocytes, replicated, and suppressed some apoptotic pathways, without significantly affecting embryonic development. (C) 2010 Elsevier Inc. All rights reserved.

Pyramidal neurones in macaque visual cortex: Interareal phenotypic variation of dendritic branching patterns

The basal dendritic arbors of over 500-layer III pyramidal neurones of the macaque cortex were compared by fractal analyses, which provides a measure of the space filling (or branching pattern) of dendritic arbors. Fractal values (D) of individual cells were compared between the cytochrome oxidase (CO)-rich blobs and CO-poor interblobs, of middle and upper layer III, and between sublaminae, in the primary visual area (Vi). These data were compared with those in the CO compartments in the second visual area (V2), and seven other extrastriate cortical areas. (V4, MT, LIP, 7a, TEO, TE and STP). There were significant differences in the fractal dimensions, and therefore the dendritic branching patterns, of cells in striate and extrastriate areas. Of the 55 possible pairwise comparisons of fractal dimension of neurones in different cortical areas (or CO compartments), 39 proved to be significantly different. The markedly different morphologies of pyramidal cells in the different cortical areas may be one of the features that determine the functional signatures of these cells by influencing the number of inputs received by, and propagation of potentials through, their dendritic arbors.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Chebyshev real wave packet propagation: H+O-2 (J=0) state-to-state reactive scattering calculations

In this paper we explore the relative performance of two recently developed wave packet methodologies for reactive scattering, namely the real wave packet Chebyshev domain propagation of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)] and the Lanczos subspace wave packet approach of Smith [J. Chem. Phys. 116, 2354 (2002); Chem. Phys. Lett. 336, 149 (2001)]. In the former method, a modified Schrodinger equation is employed to propagate the real part of the wave packet via the well-known Chebyshev iteration. While the time-dependent wave packet from the modified Schrodinger equation is different from that obtained using the standard Schrodinger equation, time-to-energy Fourier transformation yields wave functions which differ only trivially by normalization. In the Lanczos subspace approach the linear system of equations defining the action of the Green operator may be solved via either time-dependent or time-independent methods, both of which are extremely efficient due to the simple tridiagonal structure of the Hamiltonian in the Lanczos representation. The two different wave packet methods are applied to three dimensional reactive scattering of H+O-2 (total J=0). State-to-state reaction probabilities, product state distributions, as well as initial-state-resolved cumulative reaction probabilities are examined. (C) 2002 American Institute of Physics.