Data collapse, scaling functions, and analytical solutions of generalized growth models

We consider a nontrivial one-species population dynamics model with finite and infinite carrying capacities. Time-dependent intrinsic and extrinsic growth rates are considered in these models. Through the model per capita growth rate we obtain a heuristic general procedure to generate scaling functions to collapse data into a simple linear behavior even if an extrinsic growth rate is included. With this data collapse, all the models studied become independent from the parameters and initial condition. Analytical solutions are found when time-dependent coefficients are considered. These solutions allow us to perceive nontrivial transitions between species extinction and survival and to calculate the transition's critical exponents. Considering an extrinsic growth rate as a cancer treatment, we show that the relevant quantity depends not only on the intensity of the treatment, but also on when the cancerous cell growth is maximum.

Transmission electron microscopy for characterization of acrosomal damage after Percoll gradient centrifugation of cryopreserved bovine spermatozoa

The objective of this study was to characterize acrosomal ultrastructure following discontinuous Percoll gradient centrifugation of cryopreserved bovine sperm. Semen was collected from six bulls of different breeds and three ejaculates per bull were evaluated. Frozen semen samples were thawed and the acrosomal region of sperm cells was evaluated by transmission electron microscopy (TEM) before (n = 18) and after (n = 18) Percoll centrifugation. The evaluation of 20 sperm heads from each of the 36 samples analyzed ensured that a large number of cells were investigated. The data were subjected to analysis of variance at a level of significance of 5%. Percoll centrifugation reduced the percentage of sperm exhibiting normal acrosomes (from 61.77 to 30.24%), reduced the percentage of sperm presenting atypical acrosome reactions (from 28.38 to 4.84%) and increased the percentage of sperm exhibiting damage in the acrosome (from 6.14 to 64.26%). The percentage of sperm with typical acrosome reactions was not significantly different before (3.70%) and after (0.67%) centrifugation. TEM distinguished four different types of acrosomal status and enabled ultrastructural characterization of acrosomal injuries. The percentage of sperm exhibiting normal acrosomes decreased and damage in the acrosome was the most frequent acrosomal injury with the Percoll gradient centrifugation protocol utilized.

Generalized connectivity between any two nodes in a complex network

This article focuses on the identification of the number of paths with different lengths between pairs of nodes in complex networks and how these paths can be used for characterization of topological properties of theoretical and real-world complex networks. This analysis revealed that the number of paths can provide a better discrimination of network models than traditional network measurements. In addition, the analysis of real-world networks suggests that the long-range connectivity tends to be limited in these networks and may be strongly related to network growth and organization.

Formation of vortices in a dense Bose-Einstein condensate

A relaxation method is employed to study a rotating dense Bose-Einstein condensate beyond the Thomas-Fermi approximation. We use a slave-boson model to describe the strongly interacting condensate and derive a generalized nonlinear Schrodinger equation with a kinetic term for the rotating condensate. In comparison with previous calculations, based on the Thomas-Fermi approximation, significant improvements are found in regions where the condensate in a trap potential is not smooth. The critical angular velocity of the vortex formation is higher than in the Thomas-Fermi prediction.

Generalized Sznajd model for opinion propagation

In the last decade the Sznajd model has been successfully employed in modeling some properties and scale features of both proportional and majority elections. We propose a version of the Sznajd model with a generalized bounded confidence rule-a rule that limits the convincing capability of agents and that is essential to allow coexistence of opinions in the stationary state. With an appropriate choice of parameters it can be reduced to previous models. We solved this model both in a mean-field approach (for an arbitrary number of opinions) and numerically in a Barabaacutesi-Albert network (for three and four opinions), studying the transient and the possible stationary states. We built the phase portrait for the special cases of three and four opinions, defining the attractors and their basins of attraction. Through this analysis, we were able to understand and explain discrepancies between mean-field and simulation results obtained in previous works for the usual Sznajd model with bounded confidence and three opinions. Both the dynamical system approach and our generalized bounded confidence rule are quite general and we think it can be useful to the understanding of other similar models.

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.

Hydro-kinetic approach to relativistic heavy ion collisions

We develop a combined hydro-kinetic approach which incorporates a hydrodynamical expansion of the systems formed in A + A collisions and their dynamical decoupling described by escape probabilities. The method corresponds to a generalized relaxation time (tau(rel)) approximation for the Boltzmann equation applied to inhomogeneous expanding systems; at small tau(rel) it also allows one to catch the viscous effects in hadronic component-hadron-resonance gas. We demonstrate how the approximation of sudden freeze-out can be obtained within this dynamical picture of continuous emission and find that hypersurfaces, corresponding to a sharp freeze-out limit, are momentum dependent. The pion m(T) spectra are computed in the developed hydro-kinetic model, and compared with those obtained from ideal hydrodynamics with the Cooper-Frye isothermal prescription. Our results indicate that there does not exist a universal freeze-out temperature for pions with different momenta, and support an earlier decoupling of higher p(T) particles. By performing numerical simulations for various initial conditions and equations of state we identify several characteristic features of the bulk QCD matter evolution preferred in view of the current analysis of heavy ion collisions at RHIC energies.

Approximation to density functional theory for the calculation of band gaps of semiconductors

The local-density approximation (LDA) together with the half occupation (transitionstate) is notoriously successful in the calculation of atomic ionization potentials. When it comes to extended systems, such as a semiconductor infinite system, it has been very difficult to find a way to half ionize because the hole tends to be infinitely extended (a Bloch wave). The answer to this problem lies in the LDA formalism itself. One proves that the half occupation is equivalent to introducing the hole self-energy (electrostatic and exchange correlation) into the Schrodinger equation. The argument then becomes simple: The eigenvalue minus the self-energy has to be minimized because the atom has a minimal energy. Then one simply proves that the hole is localized, not infinitely extended, because it must have maximal self-energy. Then one also arrives at an equation similar to the self- interaction correction equation, but corrected for the removal of just 1/2 electron. Applied to the calculation of band gaps and effective masses, we use the self- energy calculated in atoms and attain a precision similar to that of GW, but with the great advantage that it requires no more computational effort than standard LDA.

Nonuniversal behavior for aperiodic interactions within a mean-field approximation

We study the spin-1/2 Ising model on a Bethe lattice in the mean-field limit, with the interaction constants following one of two deterministic aperiodic sequences, the Fibonacci or period-doubling one. New algorithms of sequence generation were implemented, which were fundamental in obtaining long sequences and, therefore, precise results. We calculate the exact critical temperature for both sequences, as well as the critical exponents beta, gamma, and delta. For the Fibonacci sequence, the exponents are classical, while for the period-doubling one they depend on the ratio between the two exchange constants. The usual relations between critical exponents are satisfied, within error bars, for the period-doubling sequence. Therefore, we show that mean-field-like procedures may lead to nonclassical critical exponents.

Abundance of sedentary consumers and sessile organisms along the wave exposure gradient of subtropical rocky shores of the south-west Atlantic

Sedentary consumers play an important role on populations of prey and, hence, their patterns of abundance, distribution and coexistence on shores are important to evaluate their potential influence on ecosystem dynamics. Here, we aimed to describe their spatio-temporal distribution and abundance in relation to wave exposure in the intertidal rocky shores of the south-west Atlantic to provide a basis for further understanding of ecological processes in this system. The abundance and composition of the functional groups of sessile organisms and sedentary consumers were taken by sampling the intertidal of sheltered and moderately exposed shores during a period of one year. The sublittoral fringe of sheltered areas was dominated by macroalgae, while the low midlittoral was dominated by bare rock and barnacles. In contrast, filter-feeding animals prevailed at exposed shores, probably explaining the higher abundance of the predator Stramonita haemastoma at these locations. Limpets were more abundant at the midlittoral zone of all shores while sea urchins were exclusively found at the sublittoral fringe of moderately exposed shores, therefore, adding grazing pressure on these areas. The results showed patterns of coexistence, distribution and abundance of those organisms in this subtropical area, presumably as a result of wave action, competition and prey availability. It also brought insights on the influence of top-down and bottom-up processes in this area.

Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets

Outgassing of carbon dioxide (CO(2)) from rivers and streams to the atmosphere is a major loss term in the coupled terrestrial-aquatic carbon cycle of major low-gradient river systems (the term ""river system"" encompasses the rivers and streams of all sizes that compose the drainage network in a river basin). However, the magnitude and controls on this important carbon flux are not well quantified. We measured carbon dioxide flux rates (F(CO2)), gas transfer velocity (k), and partial pressures (p(CO2)) in rivers and streams of the Amazon and Mekong river systems in South America and Southeast Asia, respectively. F(CO2) and k values were significantly higher in small rivers and streams (channels <100 m wide) than in large rivers (channels >100 m wide). Small rivers and streams also had substantially higher variability in k values than large rivers. Observed F(CO2) and k values suggest that previous estimates of basinwide CO(2) evasion from tropical rivers and wetlands have been conservative and are likely to be revised upward substantially in the future. Data from the present study combined with data compiled from the literature collectively suggest that the physical control of gas exchange velocities and fluxes in low-gradient river systems makes a transition from the dominance of wind control at the largest spatial scales (in estuaries and river mainstems) toward increasing importance of water current velocity and depth at progressively smaller channel dimensions upstream. These results highlight the importance of incorporating scale-appropriate k values into basinwide models of whole ecosystem carbon balance.

Forest structure and live aboveground biomass variation along an elevational gradient of tropical Atlantic moist forest (Brazil)

Live aboveground biomass (AGB) is an important source of uncertainty in the carbon balance from the tropical regions in part due scarcity of reliable estimates of live AGB and its variation across landscapes and forest types. Studies of forest structure and biomass stocks of Neotropical forests are biased toward Amazonian and Central American sites. In particular, standardized estimates of aboveground biomass stocks for the Brazilian Atlantic forest are rarely available. Notwithstanding the role of environmental variables that control the distribution and abundance of biomass in tropical lowland forests has been the subject of considerable research, the effect of short, steep elevational gradients on tropical forest structure and carbon dynamics is not well known. In order to evaluate forest structure and live AGB variation along an elevational gradient (0-1100 m a.s.l.) of coastal Atlantic Forest in SE Brazil, we carried out a standard census of woody stems >= 4.8 cm dbh in 13 1-ha permanent plots established on four different sites in 2006-2007. Live AGB ranged from 166.3 Mg ha(-1) (bootstrapped 95% CI: 1444,187.0) to 283.2 Mg ha(-1) (bootstrapped 95% CI: 253.0,325.2) and increased with elevation. We found that local-scale topographic variation associated with elevation influences the distribution of trees >50 cm dbh and total live AGB. Across all elevations, we found more stems (64-75%) with limited crown illumination but the largest proportion of the live AGB (68-85%) was stored in stems with highly illuminated or fully exposed crowns. Topography, disturbance and associated changes in light and nutrient supply probably control biomass distribution along this short but representative elevational gradient. Our findings also showed that intact Atlantic forest sites stored substantial amounts of carbon aboveground. The live tree AGB of the stands was found to be lower than Central Amazonian forests, but within the range of Neotropical forests, in particular when compared to Central American forests. Our comparative data suggests that differences in live tree AGB among Neotropical forests are probably related to the heterogeneous distribution of large and medium-sized diameter trees within forests and how the live biomass is partitioned among those size classes, in accordance with general trends found by previous studies. In addition, the elevational variation in live AGB stocks suggests a large spatial variability over coastal Atlantic forests in Brazil, clearly indicating that it is important to consider regional differences in biomass stocks for evaluating the role of this threatened tropical biome in the global carbon cycle. (C) 2010 Elsevier B.V. All rights reserved.

Sequential injections as an alternative to gradient exploitation for implementing differential kinetic analysis in a flow injection system

A novel flow-based strategy for implementing simultaneous determinations of different chemical species reacting with the same reagent(s) at different rates is proposed and applied to the spectrophotometric catalytic determination of iron and vanadium in Fe-V alloys. The method relies on the influence of Fe(II) and V(IV) on the rate of the iodide oxidation by Cr(VI) under acidic conditions, the Jones reducing agent is then needed Three different plugs of the sample are sequentially inserted into an acidic KI reagent carrier stream, and a confluent Cr(VI) solution is added downstream Overlap between the inserted plugs leads to a complex sample zone with several regions of maximal and minimal absorbance values. Measurements performed on these regions reveal the different degrees of reaction development and tend to be more precise Data are treated by multivariate calibration involving the PLS algorithm The proposed system is very simple and rugged Two latent variables carried out ca 95% of the analytical information and the results are in agreement with ICP-OES. (C) 2010 Elsevier B V. All rights reserved.

A hybrid-Trefftz formulation for plane elasticity with selective enrichment of the approximations

This work is related to the so-called non-conventional finite element formulations. Essentially, a methodology for the enrichment of the initial approximation which is typical of the meshless methods and based on the clouds concept is introduced in the hybrid-Trefftz formulation for plane elasticity. The formulation presented allows for the approximation and direct enrichment of two independent fields: stresses in the domains and displacements on the boundaries of the elements. Defined by a set of elements and interior boundaries sharing a common node, the cloud notion is employed to select the enrichment support for the approximation fields. The numerical analysis performed reveals an excellent performance of the resulting formulation, characterized by the good approximation ability and a reduced computational effort. Copyright (C) 2009 John Wiley & Sons, Ltd.

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.