Dynamical behaviour on the parameter space: new populational growth models proportional to beta densities

We present new populational growth models, generalized logistic models which are proportional to beta densities with shape parameters p and 2, where p > 1, with Malthusian parameter r. The complex dynamical behaviour of these models is investigated in the parameter space (r, p), in terms of topological entropy, using explicit methods, when the Malthusian parameter r increases. This parameter space is split into different regions, according to the chaotic behaviour of the models.

Modeling Allee Effect from Beta(p, 2) Densities

In this work we develop and investigate generalized populational growth models, adjusted from Beta(p, 2) densities, with Allee effect. The use of a positive parameter leads the presented generalization, which yields some more flexible models with variable extinction rates. An Allee limit is incorporated so that the models under study have strong Allee effect.

Densities and refractive indices for the ternary mixture methanol/propan-1-Ol/acetonitrile

Refractive indices, n(D), and densities, rho, at 298.15 K were measured for the ternary mixture methanol (MeOH)/propan-1-ol (1-PrOH)/acetonitrile (MeCN) for a total of 22 mole fractions, along with 18 mole fractions of each of the corresponding binary mixtures, methanol/propan-1-ol, propan-1-ol/acetonitrile and methanol/acetonitrile. The variation of excess refractive indices and excess molar volumes with composition was modeled by the Redlich-Kister polynomial function in the case of binary mixtures and by the Cibulka equation for the ternary mixture. A thermodynamic approach to excess refractive indices, recently proposed by other authors, was applied for the first time to ternary liquid mixtures. Structural effects were identified and interpreted both in the binary and ternary systems. A complex relationship between excess refractive indices and excess molar volumes was identified, revealing all four possible sign combinations between these two properties. Structuring of the mixtures was also discussed on the basis of partial molar volumes of the binary and ternary mixtures.

Generalized models from Beta(p,2) densities with strong allee effect: dynamical approach

A dynamical approach to study the behaviour of generalized populational growth models from Bets(p, 2) densities, with strong Allee effect, is presented. The dynamical analysis of the respective unimodal maps is performed using symbolic dynamics techniques. The complexity of the correspondent discrete dynamical systems is measured in terms of topological entropy. Different populational dynamics regimes are obtained when the intrinsic growth rates are modified: extinction, bistability, chaotic semistability and essential extinction.

Equilibrium self-assembly of colloids with distinct interaction sites: Thermodynamics, percolation, and cluster distribution functions

We calculate the equilibrium thermodynamic properties, percolation threshold, and cluster distribution functions for a model of associating colloids, which consists of hard spherical particles having on their surfaces three short-ranged attractive sites (sticky spots) of two different types, A and B. The thermodynamic properties are calculated using Wertheim's perturbation theory of associating fluids. This also allows us to find the onset of self-assembly, which can be quantified by the maxima of the specific heat at constant volume. The percolation threshold is derived, under the no-loop assumption, for the correlated bond model: In all cases it is two percolated phases that become identical at a critical point, when one exists. Finally, the cluster size distributions are calculated by mapping the model onto an effective model, characterized by a-state-dependent-functionality (f) over bar and unique bonding probability (p) over bar. The mapping is based on the asymptotic limit of the cluster distributions functions of the generic model and the effective parameters are defined through the requirement that the equilibrium cluster distributions of the true and effective models have the same number-averaged and weight-averaged sizes at all densities and temperatures. We also study the model numerically in the case where BB interactions are missing. In this limit, AB bonds either provide branching between A-chains (Y-junctions) if epsilon(AB)/epsilon(AA) is small, or drive the formation of a hyperbranched polymer if epsilon(AB)/epsilon(AA) is large. We find that the theoretical predictions describe quite accurately the numerical data, especially in the region where Y-junctions are present. There is fairly good agreement between theoretical and numerical results both for the thermodynamic (number of bonds and phase coexistence) and the connectivity properties of the model (cluster size distributions and percolation locus).

Learning dependent sources using mixtures of Dirichlet: applications on hyperspectral unmixing

This paper is an elaboration of the DECA algorithm [1] to blindly unmix hyperspectral data. The underlying mixing model is linear, meaning that each pixel is a linear mixture of the endmembers signatures weighted by the correspondent abundance fractions. The proposed method, as DECA, is tailored to highly mixed mixtures in which the geometric based approaches fail to identify the simplex of minimum volume enclosing the observed spectral vectors. We resort then to a statitistical framework, where the abundance fractions are modeled as mixtures of Dirichlet densities, thus enforcing the constraints on abundance fractions imposed by the acquisition process, namely non-negativity and constant sum. With respect to DECA, we introduce two improvements: 1) the number of Dirichlet modes are inferred based on the minimum description length (MDL) principle; 2) The generalized expectation maximization (GEM) algorithm we adopt to infer the model parameters is improved by using alternating minimization and augmented Lagrangian methods to compute the mixing matrix. The effectiveness of the proposed algorithm is illustrated with simulated and read data.

Ventilation influence in occupational exposure to fungi and volatile organic compounds: poultry case

Introduction - In poultry houses, large-scale production has led to increased bird densities within buildings. Such high densities of animals kept within confined spaces are a source of human health problems related to occupational organic dust exposure. This organic dust is composed of both non-viable particles and viable particulate matter (also called bioaerosols). Bioaerosols are comprised by airborne bacteria, fungi, viruses and their by-products, endotoxins and mycotoxins. Exposure to fungi in broiler houses may vary depending upon the applied ventilation system. Ventilation can be an important resource in order to reduce air contamination in these type of settings. Nevertheless, some concerns regarding costs, sensitivity of the animal species to temperature differences, and also the type of building used define which type of ventilation is used. Aim of the study - A descriptive study was developed in one poultry unit aiming to assess occupational fungal and volatile organic compounds (VOCs) exposure.

Cape Verde hotspot from the upper crust to the top of the lower mantle

We investigate the crust, upper mantle and mantle transition zone of the Cape Verde hotspot by using seismic P and S receiver functions from several tens of local seismograph stations. We find a strong discontinuity at a depth of similar to 10 km underlain by a similar to 15-km thick layer with a high (similar to 1.9) Vp/Vs velocity ratio. We interpret this discontinuity and the underlying layer as the fossil Moho, inherited from the pre-hotspot era, and the plume-related magmatic underplate. Our uppermost-mantle models are very different from those previously obtained for this region: our S velocity is much lower and there are no indications of low densities. Contrary to previously published arguments for the standard transition zone thickness our data indicate that this thickness under the Cape Verde islands is up to similar to 30 km less than in the ambient mantle. This reduction is a combined effect of a depression of the 410-km discontinuity and an uplift of the 660-km discontinuity. The uplift is in contrast to laboratory data and some seismic data on a negligible dependence of depth of the 660-km discontinuity on temperature in hotspots. A large negative pressure-temperature slope which is suggested by our data implies that the 660-km discontinuity may resist passage of the plume. Our data reveal beneath the islands a reduction of S velocity of a few percent between 470-km and 510-km depths. The low velocity layer in the upper transition zone under the Cape Verde archipelago is very similar to that previously found under the Azores and a few other hotspots. In the literature there are reports on a regional 520-km discontinuity, the impedance of which is too large to be explained by the known phase transitions. Our observations suggest that the 520-km discontinuity may present the base of the low-velocity layer in the transition zone. (C) 2011 Elsevier B.V. All rights reserved.

Patching up dipoles: can dipolar particles be viewed as patchy colloids?

We investigate whether the liquid-vapour phase transition of strongly dipolar fluids can be understood using a model of patchy colloids. These consist of hard spherical particles with three short-ranged attractive sites (patches) on their surfaces. Two of the patches are of type A and one is of type B. Patches A on a particle may bond either to a patch A or to a patch B on another particle. Formation of an AA (AB) bond lowers the energy by epsilon AA (epsilon AB). In the limit [image omitted], this patchy model exhibits condensation driven by AB-bonds (Y-junctions). Y-junctions are also present in low-density, strongly dipolar fluids, and have been conjectured to play a key role in determining their critical behaviour. We map the dipolar Yukawa hard-sphere (DYHS) fluid onto this 2A + 1B patchy model by requiring that the latter reproduce the correct DYHS critical point as a function of the isotropic interaction strength epsilon Y. This is achieved for sensible values of epsilon AB and the bond volumes. Results for the internal energy and the particle coordination number are in qualitative agreement with simulations of DYHSs. Finally, by taking the limit [image omitted], we arrive at a new estimate for the critical point of the dipolar hard-sphere fluid, which agrees with extrapolations from simulation.

Re-entrant phase behaviour of network fluids: A patchy particle model with temperature-dependent valence

We study a model consisting of particles with dissimilar bonding sites ("patches"), which exhibits self-assembly into chains connected by Y-junctions, and investigate its phase behaviour by both simulations and theory. We show that, as the energy cost epsilon(j) of forming Y-junctions increases, the extent of the liquid-vapour coexistence region at lower temperatures and densities is reduced. The phase diagram thus acquires a characteristic "pinched" shape in which the liquid branch density decreases as the temperature is lowered. To our knowledge, this is the first model in which the predicted topological phase transition between a fluid composed of short chains and a fluid rich in Y-junctions is actually observed. Above a certain threshold for epsilon(j), condensation ceases to exist because the entropy gain of forming Y-junctions can no longer offset their energy cost. We also show that the properties of these phase diagrams can be understood in terms of a temperature-dependent effective valence of the patchy particles. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3605703]

Arrow Plot: a new graphical tool for selecting up and down regulated genes and genes differentially expressed on samples subgroups

Background: A common task in analyzing microarray data is to determine which genes are differentially expressed across two (or more) kind of tissue samples or samples submitted under experimental conditions. Several statistical methods have been proposed to accomplish this goal, generally based on measures of distance between classes. It is well known that biological samples are heterogeneous because of factors such as molecular subtypes or genetic background that are often unknown to the experimenter. For instance, in experiments which involve molecular classification of tumors it is important to identify significant subtypes of cancer. Bimodal or multimodal distributions often reflect the presence of subsamples mixtures. Consequently, there can be genes differentially expressed on sample subgroups which are missed if usual statistical approaches are used. In this paper we propose a new graphical tool which not only identifies genes with up and down regulations, but also genes with differential expression in different subclasses, that are usually missed if current statistical methods are used. This tool is based on two measures of distance between samples, namely the overlapping coefficient (OVL) between two densities and the area under the receiver operating characteristic (ROC) curve. The methodology proposed here was implemented in the open-source R software. Results: This method was applied to a publicly available dataset, as well as to a simulated dataset. We compared our results with the ones obtained using some of the standard methods for detecting differentially expressed genes, namely Welch t-statistic, fold change (FC), rank products (RP), average difference (AD), weighted average difference (WAD), moderated t-statistic (modT), intensity-based moderated t-statistic (ibmT), significance analysis of microarrays (samT) and area under the ROC curve (AUC). On both datasets all differentially expressed genes with bimodal or multimodal distributions were not selected by all standard selection procedures. We also compared our results with (i) area between ROC curve and rising area (ABCR) and (ii) the test for not proper ROC curves (TNRC). We found our methodology more comprehensive, because it detects both bimodal and multimodal distributions and different variances can be considered on both samples. Another advantage of our method is that we can analyze graphically the behavior of different kinds of differentially expressed genes. Conclusion: Our results indicate that the arrow plot represents a new flexible and useful tool for the analysis of gene expression profiles from microarrays.

Comparação do cálculo de dose com os algoritmos Analytical Anisotropic Algorithm e Pencil Beam Convolution na patologia de pulmão

Mestrado em Radioterapia.

Dynamical analysis in growth models: Blumberg’s equation

We present a new dynamical approach to the Blumberg's equation, a family of unimodal maps. These maps are proportional to Beta(p, q) probability densities functions. Using the symmetry of the Beta(p, q) distribution and symbolic dynamics techniques, a new concept of mirror symmetry is defined for this family of maps. The kneading theory is used to analyze the effect of such symmetry in the presented models. The main result proves that two mirror symmetric unimodal maps have the same topological entropy. Different population dynamics regimes are identified, when the intrinsic growth rate is modified: extinctions, stabilities, bifurcations, chaos and Allee effect. To illustrate our results, we present a numerical analysis, where are demonstrated: monotonicity of the topological entropy with the variation of the intrinsic growth rate, existence of isentropic sets in the parameters space and mirror symmetry.

Three-dimensional patchy lattice model for empty fluids

The phase diagram of a simple model with two patches of type A and ten patches of type B (2A10B) on the face centred cubic lattice has been calculated by simulations and theory. Assuming that there is no interaction between the B patches the behavior of the system can be described in terms of the ratio of the AB and AA interactions, r. Our results show that, similarly to what happens for related off-lattice and two-dimensional lattice models, the liquid-vapor phase equilibria exhibit reentrant behavior for some values of the interaction parameters. However, for the model studied here the liquid-vapor phase equilibria occur for values of r lower than 1/3, a threshold value which was previously thought to be universal for 2AnB models. In addition, the theory predicts that below r = 1/3 (and above a new condensation threshold which is < 1/3) the reentrant liquid-vapor equilibria are so extreme that it exhibits a closed loop with a lower critical point, a very unusual behavior in single-component systems. An order-disorder transition is also observed at higher densities than the liquid-vapor equilibria, which shows that the liquid-vapor reentrancy occurs in an equilibrium region of the phase diagram. These findings may have implications in the understanding of the condensation of dipolar hard spheres given the analogy between that system and the 2AnB models considered here. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4771591]

An Extension of Gompertzian Growth Dynamics Weibull and Frechet Models

In this work a new probabilistic and dynamical approach to an extension of the Gompertz law is proposed. A generalized family of probability density functions, designated by Beta* (p, q), which is proportional to the right hand side of the Tsoularis-Wallace model, is studied. In particular, for p = 2, the investigation is extended to the extreme value models of Weibull and Frechet type. These models, described by differential equations, are proportional to the hyper-Gompertz growth model. It is proved that the Beta* (2, q) densities are a power of betas mixture, and that its dynamics are determined by a non-linear coupling of probabilities. The dynamical analysis is performed using techniques of symbolic dynamics and the system complexity is measured using topological entropy. Generally, the natural history of a malignant tumour is reflected through bifurcation diagrams, in which are identified regions of regression, stability, bifurcation, chaos and terminus.