The log-bimodal-skew-normal model. A geochemical application

The main objective of this paper is to study a logarithm extension of the bimodal skew normal model introduced by Elal-Olivero et al. [1]. The model can then be seen as an alternative to the log-normal model typically used for fitting positive data. We study some basic properties such as the distribution function and moments, and discuss maximum likelihood for parameter estimation. We report results of an application to a real data set related to nickel concentration in soil samples. Model fitting comparison with several alternative models indicates that the model proposed presents the best fit and so it can be quite useful in real applications for chemical data on substance concentration. Copyright (C) 2011 John Wiley & Sons, Ltd.

Bayesian density estimation using Skew Student-t-Normal mixtures

We present a Bayesian approach for modeling heterogeneous data and estimate multimodal densities using mixtures of Skew Student-t-Normal distributions [Gomez, H.W., Venegas, O., Bolfarine, H., 2007. Skew-symmetric distributions generated by the distribution function of the normal distribution. Environmetrics 18, 395-407]. A stochastic representation that is useful for implementing a MCMC-type algorithm and results about existence of posterior moments are obtained. Marginal likelihood approximations are obtained, in order to compare mixture models with different number of component densities. Data sets concerning the Gross Domestic Product per capita (Human Development Report) and body mass index (National Health and Nutrition Examination Survey), previously studied in the related literature, are analyzed. (c) 2008 Elsevier B.V. All rights reserved.

Robustness and stabilization properties of monetary policy rules in Brazil

Based on three versions of a small macroeconomic model for Brazil, this paper presents empirical evidence on the effects of parameter uncertainty on monetary policy rules and on the robustness of optimal and simple rules over different model specifications. By comparing the optimal policy rule under parameter uncertainty with the rule calculated under purely additive uncertainty, we find that parameter uncertainty should make policymakers react less aggressively to the economy's state variables, as suggested by Brainard's "conservatism principIe", although this effect seems to be relatively small. We then informally investigate each rule's robustness by analyzing the performance of policy rules derived from each model under each one of the alternative models. We find that optimal rules derived from each model perform very poorly under alternative models, whereas a simple Taylor rule is relatively robusto We also fmd that even within a specific model, the Taylor rule may perform better than the optimal rule under particularly unfavorable realizations from the policymaker' s loss distribution function.

Statistical modelling of extreme rainfall in Madeira Island

Extreme rainfall events have triggered a significant number of flash floods in Madeira Island along its past and recent history. Madeira is a volcanic island where the spatial rainfall distribution is strongly affected by its rugged topography. In this thesis, annual maximum of daily rainfall data from 25 rain gauge stations located in Madeira Island were modelled by the generalised extreme value distribution. Also, the hypothesis of a Gumbel distribution was tested by two methods and the existence of a linear trend in both distributions parameters was analysed. Estimates for the 50– and 100–year return levels were also obtained. Still in an univariate context, the assumption that a distribution function belongs to the domain of attraction of an extreme value distribution for monthly maximum rainfall data was tested for the rainy season. The available data was then analysed in order to find the most suitable domain of attraction for the sampled distribution. In a different approach, a search for thresholds was also performed for daily rainfall values through a graphical analysis. In a multivariate context, a study was made on the dependence between extreme rainfall values from the considered stations based on Kendall’s τ measure. This study suggests the influence of factors such as altitude, slope orientation, distance between stations and their proximity of the sea on the spatial distribution of extreme rainfall. Groups of three pairwise associated stations were also obtained and an adjustment was made to a family of extreme value copulas involving the Marshall–Olkin family, whose parameters can be written as a function of Kendall’s τ association measures of the obtained pairs.

Análise por meio de espectroscopia de emissão óptica das espécies ativas em nitretação iônica e gaiola catódica

Plasma process like ionic nitriding and cathodic cage plasma nitriding are utilized in order to become hard surface of steels. The ionic nitriding is already accepted in the industry while cathodic cage plasma nitriding process is in industrial implementation stage. Those process depend of plasma parameters like electronic and ionic temperature (Te, Ti), species density (ne, ni) and of distribution function of these species. In the present work, the plasma used to those two processes has been observed through Optical Emission Spectroscopy OES technique in order to identify presents species in the treatment ambient and relatively quantify them. So plasma of typical mixtures like N2 H2 has been monitored through in order to study evolution of those species during the process. Moreover, it has been realized a systematic study about leaks, also thought OES, that accomplish the evolution of contaminant species arising because there is flux of atmosphere to inside nitriding chamber and in what conditions the species are sufficiently reduced. Finally, to describe the physic mechanism that acts on both coating techniques ionic nitriding and cathodic cage plasma nitriding

Simulação de fluxo de fluidos em meios porosos desordenados uma análise de efeito de escala na estimativa da permeabilidade e do coeficiente de arrasto

The present study provides a methodology that gives a predictive character the computer simulations based on detailed models of the geometry of a porous medium. We using the software FLUENT to investigate the flow of a viscous Newtonian fluid through a random fractal medium which simplifies a two-dimensional disordered porous medium representing a petroleum reservoir. This fractal model is formed by obstacles of various sizes, whose size distribution function follows a power law where exponent is defined as the fractal dimension of fractionation Dff of the model characterizing the process of fragmentation these obstacles. They are randomly disposed in a rectangular channel. The modeling process incorporates modern concepts, scaling laws, to analyze the influence of heterogeneity found in the fields of the porosity and of the permeability in such a way as to characterize the medium in terms of their fractal properties. This procedure allows numerically analyze the measurements of permeability k and the drag coefficient Cd proposed relationships, like power law, for these properties on various modeling schemes. The purpose of this research is to study the variability provided by these heterogeneities where the velocity field and other details of viscous fluid dynamics are obtained by solving numerically the continuity and Navier-Stokes equations at pore level and observe how the fractal dimension of fractionation of the model can affect their hydrodynamic properties. This study were considered two classes of models, models with constant porosity, MPC, and models with varying porosity, MPV. The results have allowed us to find numerical relationship between the permeability, drag coefficient and the fractal dimension of fractionation of the medium. Based on these numerical results we have proposed scaling relations and algebraic expressions involving the relevant parameters of the phenomenon. In this study analytical equations were determined for Dff depending on the geometrical parameters of the models. We also found a relation between the permeability and the drag coefficient which is inversely proportional to one another. As for the difference in behavior it is most striking in the classes of models MPV. That is, the fact that the porosity vary in these models is an additional factor that plays a significant role in flow analysis. Finally, the results proved satisfactory and consistent, which demonstrates the effectiveness of the referred methodology for all applications analyzed in this study.

Qualidade do gasto público em educação para os municípios da região metropolitana de Natal-2007

The aim of this thesis is to evaluate the quality of public spending on education for the municipalities of the Metropolitan Region of Natal (RMN) in 2009 by use of two theories: The Theory of Welfare (Welfare State) and the Public Choice Theory (TEP), both important to understand the relationship between education and economics. The study also uses principles of microeconomics and public sector economics to get a better idea of the role of education in economy and society. It describes the development of the educational policy in Brazil from 1988 to the Federal Constitution of 2010, following the major changes in basic education during each government. The characteristics of the RMN municipalities were illustrated with socioeconomic indicators, while educational indicators were used to characterize each municipality regarding education. The model used in this study was developed by Bertê, Brunet and Borges, the data was collected on the back of the School Census 2009 and the Brazil Exam 2009 and it was processed quantitavely in the Information System on Public Budgets in Education (SIOPE) by use of the statistical method called standardized score of the normal cumulative distribution function. The quality of public spending on education is the result of the relation between performance indicator ratio and expense ratio. For the qualitative analysis of results, the criteria of efficiency, efficacy and effectiveness were used. The study found that municipalities with higher expenses showed a worse quality of spending and failed to convert the expenditure incurred into performance, thus confirming ineffectiveness

Teoria cinética não extensiva: efeitos físicos em gases e plasmas

The standard kinetic theory for a nonrelativistic diluted gas is generalized in the spirit of the nonextensive statistic distribution introduced by Tsallis. The new formalism depends on an arbitrary q parameter measuring the degree of nonextensivity. In the limit q = 1, the extensive Maxwell-Boltzmann theory is recovered. Starting from a purely kinetic deduction of the velocity q-distribution function, the Boltzmann H-teorem is generalized for including the possibility of nonextensive out of equilibrium effects. Based on this investigation, it is proved that Tsallis' distribution is the necessary and sufficient condition defining a thermodynamic equilibrium state in the nonextensive context. This result follows naturally from the generalized transport equation and also from the extended H-theorem. Two physical applications of the nonextensive effects have been considered. Closed analytic expressions were obtained for the Doppler broadening of spectral lines from an excited gas, as well as, for the dispersion relations describing the eletrostatic oscillations in a diluted electronic plasma. In the later case, a comparison with the experimental results strongly suggests a Tsallis distribution with the q parameter smaller than unity. A complementary study is related to the thermodynamic behavior of a relativistic imperfect simple fluid. Using nonequilibrium thermodynamics, we show how the basic primary variables, namely: the energy momentum tensor, the particle and entropy fluxes depend on the several dissipative processes present in the fluid. The temperature variation law for this moving imperfect fluid is also obtained, and the Eckart and Landau-Lifshitz formulations are recovered as particular cases

Monte Carlo investigations of intermolecular interactions in water-amide mixtures

Monte Carlo simulations of water-amides (amide=fonnamide-FOR, methylfonnamide-NMF and dimethylformamide-DMF) solutions have been carried out in the NpT ensemble at 308 K and 1 atm. The structure and excess enthalpy of the mixtures as a function of the composition have been investigated. The TIP4P model was used for simulating water and six-site models previously optimized in this laboratory were used for simulating the liquid amides. The intermolecular interaction energy was calculated using the classical 6-12 Lennard-Jones potential plus a Coulomb term. The interaction energy between solute and solvent has been partitioned what leads to a better understanding of the behavior of the enthalpy of mixture obtained for the three solutions experimentally. Radial distribution functions for the water-amides correlations permit to explore the intermolecular interactions between the molecules. The results show that three, two and one hydrogen bonds between the water and the amide molecules are formed in the FOR, NMF and DMF-water solutions, respectively. These H-bonds are, respectively, stronger for DMF-water, NMF-water and FOR-water. In the NMF-water solution, the interaction between the methyl group of the NMF and the oxygen of the water plays a role in the stabilization of the aqueous solution quite similar to that of an H-bond in the FOR-water solution. (c) 2005 Elsevier B.V. All rights reserved.

Modeling the dielectric response of lanthanum modified lead zirconate titanate ferroelectric ceramics-an approach to the phase transitions in relaxor ferroelectrics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Monte Carlo thermodynamic and structural properties of the TIP4P water model: dependence on the computational conditions

Classical Monte Carlo simulations were carried out on the NPT ensemble at 25°C and 1 atm, aiming to investigate the ability of the TIP4P water model [Jorgensen, Chandrasekhar, Madura, Impey and Klein; J. Chem. Phys., 79 (1983) 926] to reproduce the newest structural picture of liquid water. The results were compared with recent neutron diffraction data [Soper; Bruni and Ricci; J. Chem. Phys., 106 (1997) 247]. The influence of the computational conditions on the thermodynamic and structural results obtained with this model was also analyzed. The findings were compared with the original ones from Jorgensen et al [above-cited reference plus Mol. Phys., 56 (1985) 1381]. It is notice that the thermodynamic results are dependent on the boundary conditions used, whereas the usual radial distribution functions g(O/O(r)) and g(O/H(r)) do not depend on them.

Pion structure in the instanton liquid model

The covariant quark model of the pion based on the effective nonlocal quark-hadron Lagrangian involving nonlocality induced by instanton fluctuations of the QCD vacuum is reviewed. Explicit gauge invariant formalism allows us to construct the conserved vector and axial currents and to demonstrate their consistency with the Ward-Takahashi identities and low-energy theorems. The spontaneous breaking of chiral symmetry results in the dynamic quark mass and the vertex of the quark-pion interaction, both momentum-dependent. The parameters of the instanton vacuum, the average size of the instantons, and the effective quark mass are expressed in terms of the vacuum expectation values of the lowest dimension quark-gluon operators and low-energy pion observables. The transition pion form factor for the processes gamma*gamma --> pi (0) and gamma*gamma* --> pi (0) is analyzed in detail. The kinematic dependence of the transition form factor at high momentum transfers allows one to determine the relationship between the light-cone amplitude of the quark distribution in the pion and the quark-pion vertex function. Its dynamic dependence implies that the transition form factor gamma*gamma --> pi (0) at high momentum transfers is acutely sensitive to the size of the nonlocality of nonperturbative fluctuations in the QCD vacuum. In the leading twist, the distribution amplitude and the distribution function of the valence quarks in the pion are calculated at a low normalization point of the order of the inverse average instanton size rho (-1)(c). The QCD results are evolved to higher momentum transfers and are in reasonable agreement with available experimental data on the pion structure.

Measurement of the muon charge asymmetry from W boson decays

We present a measurement of the muon charge asymmetry from W boson decays using 0.3 fb(-1) of data collected at root s =1.96 GeV between 2002 and 2004 with the D0 detector at the Fermilab Tevatron (p) over bar Collider. We compare our findings with expectations from next-to-leading-order calculations performed using the CTEQ6.1M and MRST04 NLO parton distribution functions. Our findings can be used to constrain future parton distribution function fits.

Linking the hydrodynamic and kinetic description of a dissipative relativistic conformal theory

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Sol phase and sol-gel transition in SnO2 colloidal suspensions

The effect of concentration on the structure of SnO2 colloids in aqueous suspension, on their spatial correlation and on the gelation process was studied by small angle x-ray scattering (SAXS). The shape of the experimental SAXS curves varies with suspension concentration. For diluted suspensions ([SnO2] less than or equal to 0.13 mol L-1), SAXS results indicate the presence of colloidal fractal aggregates with an internal correlation length xi congruent to 20 Angstrom, without any noticeable spatial correlation between them. This suggests that the aggregates are spatially arranged without any significant interaction like in ideal gas structures. For higher concentrations ([SnO2] = 0.16, 0.32, and 0.64 mol L-1), the colloidal aggregates are larger (xi = 24 Angstrom) and exhibit a certain degree of spatial correlation between them. The pair correlation function corresponding to the sol with the highest concentration (0.92 mol L-1) reveals a rather strong short range order between aggregates, characteristic of a fluid-like structure, with an average nearest-neighbor distance between aggregates d(1) = 125 Angstrom and an average second-neighbor distance d(2) = 283 Angstrom. The pair distribution function remains essentially invariant during the sol-gel transition, suggesting that gelation involves the formation of a few points of connection between the aggregates resulting in a gel network constituted by essentially linear chains of clusters..