A general threshold stress hybrid hazard model for lifetime data

In this paper we propose a hybrid hazard regression model with threshold stress which includes the proportional hazards and the accelerated failure time models as particular cases. To express the behavior of lifetimes the generalized-gamma distribution is assumed and an inverse power law model with a threshold stress is considered. For parameter estimation we develop a sampling-based posterior inference procedure based on Markov Chain Monte Carlo techniques. We assume proper but vague priors for the parameters of interest. A simulation study investigates the frequentist properties of the proposed estimators obtained under the assumption of vague priors. Further, some discussions on model selection criteria are given. The methodology is illustrated on simulated and real lifetime data set.

Critical behavior in a cross-situational lexicon learning scenario

The associationist account for early word learning is based on the co-occurrence between referents and words. Here we introduce a noisy cross-situational learning scenario in which the referent of the uttered word is eliminated from the context with probability gamma, thus modeling the noise produced by out-of-context words. We examine the performance of a simple associative learning algorithm and find a critical value of the noise parameter gamma(c) above which learning is impossible. We use finite-size scaling to show that the sharpness of the transition persists across a region of order tau(-1/2) about gamma(c), where tau is the number of learning trials, as well as to obtain the learning error (scaling function) in the critical region. In addition, we show that the distribution of durations of periods when the learning error is zero is a power law with exponent -3/2 at the critical point. Copyright (C) EPLA, 2012

Quantum critical scaling at a Bose-glass/superfluid transition: Theory and experiment for a model quantum magnet

In this paper we investigate the quantum phase transition from magnetic Bose Glass to magnetic Bose-Einstein condensation induced by amagnetic field in NiCl2 center dot 4SC(NH2)(2) (dichloro-tetrakis-thiourea-nickel, or DTN), doped with Br (Br-DTN) or site diluted. Quantum Monte Carlo simulations for the quantum phase transition of the model Hamiltonian for Br-DTN, as well as for site-diluted DTN, are consistent with conventional scaling at the quantum critical point and with a critical exponent z verifying the prediction z = d; moreover the correlation length exponent is found to be nu = 0.75(10), and the order parameter exponent to be beta = 0.95(10). We investigate the low-temperature thermodynamics at the quantum critical field of Br-DTN both numerically and experimentally, and extract the power-law behavior of the magnetization and of the specific heat. Our results for the exponents of the power laws, as well as previous results for the scaling of the critical temperature to magnetic ordering with the applied field, are incompatible with the conventional crossover-scaling Ansatz proposed by Fisher et al. [Phys. Rev. B 40, 546 (1989)]. However they can all be reconciled within a phenomenological Ansatz in the presence of a dangerously irrelevant operator.

Unraveling the drivers of community dissimilarity and species extinction in fragmented landscapes

Communities in fragmented landscapes are often assumed to be structured by species extinction due to habitat loss, which has led to extensive use of the species-area relationship (SAR) in fragmentation studies. However, the use of the SAR presupposes that habitat loss leads species to extinction but does not allow for extinction to be offset by colonization of disturbed-habitat specialists. Moreover, the use of SAR assumes that species richness is a good proxy of community changes in fragmented landscapes. Here, we assessed how communities dwelling in fragmented landscapes are influenced by habitat loss at multiple scales; then we estimated the ability of models ruled by SAR and by species turnover in successfully predicting changes in community composition, and asked whether species richness is indeed an informative community metric. To address these issues, we used a data set consisting of 140 bird species sampled in 65 patches, from six landscapes with different proportions of forest cover in the Atlantic Forest of Brazil. We compared empirical patterns against simulations of over 8 million communities structured by different magnitudes of the power-law SAR and with species-specific rules to assign species to sites. Empirical results showed that, while bird community composition was strongly influenced by habitat loss at the patch and landscape scale, species richness remained largely unaffected. Modeling results revealed that the compositional changes observed in the Atlantic Forest bird metacommunity were only matched by models with either unrealistic magnitudes of the SAR or by models ruled by species turnover, akin to what would be observed along natural gradients. We show that, in the presence of such compositional turnover, species richness is poorly correlated with species extinction, and z values of the SAR strongly underestimate the effects of habitat loss. We suggest that the observed compositional changes are driven by each species reaching its individual extinction threshold: either a threshold of forest cover for species that disappear with habitat loss, or of matrix cover for species that benefit from habitat loss.

The influence of connectivity on the firing rate in a neuronal network with electrical and chemical synapses

We study the firing rate properties of a cellular automaton model for a neuronal network with chemical synapses. We propose a simple mechanism in which the nonlocal connections are included, through electrical and chemical synapses. In the latter case, we introduce a time delay which produces self-sustained activity. Nonlocal connections, or shortcuts, are randomly introduced according to a specified connection probability. There is a range of connection probabilities for which neuron firing occurs, as well as a critical probability for which the firing ceases in the absence of time delay. The critical probability for nonlocal shortcuts depends on the network size according to a power-law. We also compute the firing rate amplification factor by varying both the connection probability and the time delay for different network sizes. (C) 2011 Elsevier B.V. All rights reserved.

Residence Time Distribution Models Derived from Non-Ideal Laminar Velocity Profiles in Tubes

The theoretical E-curve for the laminar flow of non-Newtonian fluids in circular tubes may not be accurate for real tubular systems with diffusion, mechanical vibration, wall roughness, pipe fittings, curves, coils, or corrugated walls. Deviations from the idealized laminar flow reactor (LFR) cannot be well represented using the axial dispersion or the tanks-in-series models of residence time distribution (RTD). In this work, four RTD models derived from non-ideal velocity profiles in segregated tube flow are proposed. They were used to represent the RTD of three tubular systems working with Newtonian and pseudoplastic fluids. Other RTD models were considered for comparison. The proposed models provided good adjustments, and it was possible to determine the active volumes. It is expected that these models can be useful for the analysis of LFR or for the evaluation of continuous thermal processing of viscous foods.

Charge dynamics in half-filled Hubbard chains with finite on-site interaction

We study the charge dynamic structure factor of the one-dimensional Hubbard model with finite on-site repulsion U at half-filling. Numerical results from the time-dependent density matrix renormalization group are analyzed by comparison with the exact spectrum of the model. The evolution of the line shape as a function of U is explained in terms of a relative transfer of spectral weight between the two-holon continuum that dominates in the limit U -> infinity and a subset of the two-holon-two-spinon continuum that reconstructs the electron-hole continuum in the limit U -> 0. Power-law singularities along boundary lines of the spectrum are described by effective impurity models that are explicitly invariant under spin and eta-spin SU(2) rotations. The Mott-Hubbard metal-insulator transition is reflected in a discontinuous change of the exponents of edge singularities at U = 0. The sharp feature observed in the spectrum for momenta near the zone boundary is attributed to a van Hove singularity that persists as a consequence of integrability.

Low-temperature thermodynamic properties near the field-induced quantum critical point in NiCl2-4SC(NH2)(2)

We present a comprehensive experimental and theoretical investigation of the thermodynamic properties: specific heat, magnetization, and thermal expansion in the vicinity of the field-induced quantum critical point (QCP) around the lower critical field H-c1 approximate to 2 T in NiCl2-4SC(NH2)(2). A T-3/2 behavior in the specific heat and magnetization is observed at very low temperatures at H = H-c1, which is consistent with the universality class of Bose-Einstein condensation of magnons. The temperature dependence of the thermal expansion coefficient at H-c1 shows minor deviations from the expected T-1/2 behavior. Our experimental study is complemented by analytical calculations and quantum Monte Carlo simulations, which reproduce nicely the measured quantities. We analyze the thermal and the magnetic Gruneisen parameters, which are ideal quantities to identify QCPs. Both parameters diverge at H-c1 with the expected T-1 power law. By using the Ehrenfest relations at the second-order phase transition, we are able to estimate the pressure dependencies of the characteristic temperature and field scales.

Identification of literary movements using complex networks to represent texts

The use of statistical methods to analyze large databases of text has been useful in unveiling patterns of human behavior and establishing historical links between cultures and languages. In this study, we identified literary movements by treating books published from 1590 to 1922 as complex networks, whose metrics were analyzed with multivariate techniques to generate six clusters of books. The latter correspond to time periods coinciding with relevant literary movements over the last five centuries. The most important factor contributing to the distinctions between different literary styles was the average shortest path length, in particular the asymmetry of its distribution. Furthermore, over time there has emerged a trend toward larger average shortest path lengths, which is correlated with increased syntactic complexity, and a more uniform use of the words reflected in a smaller power-law coefficient for the distribution of word frequency. Changes in literary style were also found to be driven by opposition to earlier writing styles, as revealed by the analysis performed with geometrical concepts. The approaches adopted here are generic and may be extended to analyze a number of features of languages and cultures.

Evolution of pi(0) Suppression in Au plus Au Collisions from root s(NN)=39 to 200 GeV

Neutral-pion pi(0) spectra were measured at midrapidity (vertical bar y vertical bar < 0.35) in Au + Au collisions at root s(NN) = 39 and 62.4 GeV and compared with earlier measurements at 200 GeV in a transverse-momentum range of 1 < p(T) < 10 GeV/c. The high-p(T) tail is well described by a power law in all cases, and the powers decrease significantly with decreasing center-of-mass energy. The change of powers is very similar to that observed in the corresponding spectra for p + p collisions. The nuclear modification factors (RAA) show significant suppression, with a distinct energy, centrality, and p(T) dependence. Above p(T) = 7 GeV/c, R-AA is similar for root sNN = 62.4 and 200 GeV at all centralities. Perturbative-quantum-chromodynamics calculations that describe R-AA well at 200 GeV fail to describe the 39 GeV data, raising the possibility that, for the same p(T) region, the relative importance of initial-state effects and soft processes increases at lower energies. The p(T) range where pi(0) spectra in central Au + Au collisions have the same power as in p + p collisions is approximate to 5 and 7 GeV/c for root sNN = 200 and 62.4 GeV, respectively. For the root sNN = 39 GeV data, it is not clear whether such a region is reached, and the x(T) dependence of the x(T)-scaling power-law exponent is very different from that observed in the root sNN = 62 and 200 GeV data, providing further evidence that initial-state effects and soft processes mask the in-medium suppression of hardscattered partons to higher p(T) as the collision energy decreases.

Subcritical fatigue in fuse networks

We obtain the Paris law of fatigue crack propagation in a fuse network model where the accumulated damage in each resistor increases with time as a power law of the local current amplitude. When a resistor reaches its fatigue threshold, it burns irreversibly. Over time, this drives cracks to grow until the system is fractured into two parts. We study the relation between the macroscopic exponent of the crack-growth rate -entering the phenomenological Paris law-and the microscopic damage accumulation exponent, gamma, under the influence of disorder. The way the jumps of the growing crack, Delta a, and the waiting time between successive breaks, Delta t, depend on the type of material, via gamma, are also investigated. We find that the averages of these quantities, <Delta a > and <Delta t >/< t(r)>, scale as power laws of the crack length a, <Delta a > proportional to a(alpha) and <Delta t >/< t(r)> proportional to a(-beta), where < t(r)> is the average rupture time. Strikingly, our results show, for small values of gamma, a decrease in the exponent of the Paris law in comparison with the homogeneous case, leading to an increase in the lifetime of breaking materials. For the particular case of gamma = 0, when fatigue is exclusively ruled by disorder, an analytical treatment confirms the results obtained by simulation. Copyright (C) EPLA, 2012

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.

Pancreatic iron stores assessed by magnetic resonance imaging (MRI) in beta thalassemic patients

Purpose: To assess the correlation between MRI findings of the pancreas with those of the heart and liver in patients with beta thalassemia; to compare the pancreas T2* MRI results with glucose and ferritin levels and labile plasma iron (LPI). Materials and methods: We retrospectively evaluated chronically transfused patients, testing glucose with enzymatic tests, serum ferritin with chemiluminescence, LPI with cellular fluorescence, and T2* MRI to assess iron content in the heart, liver, and pancreas. MRI results were compared with one another and with serum glucose, ferritin, and LPI. Liver iron concentration (LIC) was determined in 11 patients' liver biopsies by atomic absorption spectrometry. Results: 289 MRI studies were available from 115 patients during the period studied. 9.4% of patients had overt diabetes and an additional 16% of patients had impaired fasting glucose. Both pancreatic and cardiac R2* had predictive power (p < 0.0001) for identifying diabetes. Cardiac and pancreatic R2* were modestly correlated with one another (r(2) = 0.20, p < 0.0001). Both were weakly correlated with LIC (r(2) = 0.09, p < 0.0001 for both) and serum ferritin (r(2) = 0.14, p < 0.0001 and r(2) = 0.03, p < 0.02, respectively). None of the three served as a screening tool for single observations. There is a strong log-log, or power-law, relationship between ratio of signal intensity (SIR) values and pancreas R2* with an r(2) of 0.91. Conclusions: Pancreatic iron overload can be assessed by MRI, but siderosis in other organs did not correlate significantly with pancreatic hemosiderosis. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

The spectrum of high-energy cosmic rays measured with KASCADE-Grande

The energy spectrum of cosmic rays between 10(16) eV and 10(18) eV, derived from measurements of the shower size (total number of charged particles) and the total muon number of extensive air showers by the KASCADE-Grande experiment, is described. The resulting all-particle energy spectrum exhibits strong hints for a hardening of the spectrum at approximately 2 . 10(16) eV and a significant steepening at approximate to 8 . 10(16) eV. These observations challenge the view that the spectrum is a single power law between knee and ankle. Possible scenarios generating such features are discussed in terms of astrophysical processes that may explain the transition region from galactic to extragalactic origin of cosmic rays. (C) 2012 Elsevier B.V. All rights reserved.

Density and rheological parameters of goat milk

The rheological behavior and density of goat milk was studied as a function of solids concentration (10.5 to 50.0%) and temperature (273 to 331 k). Newtonian behavior was observed for values of total solids (TS) between 10.5 and 22.0% and temperatures from 276 to 331 k changing to pseudoplastic behavior without yield stress for TS from 25.0 to 39.4% at the same range of temperature. Goat milk with TS between 44.3 to 50.0% and temperatures of 273 to 296 k showed yield stress in addition to pseudoplastic behavior. At 303 to 331 k the power law model was observed again, without yield stress. The density of goat milk ranged from 991.7 to 1232.4 kg.m-3.