Using non-homogeneous Poisson models with multiple change-points to estimate the number of ozone exceedances in Mexico City

In this paper, we consider some non-homogeneous Poisson models to estimate the probability that an air quality standard is exceeded a given number of times in a time interval of interest. We assume that the number of exceedances occurs according to a non-homogeneous Poisson process (NHPP). This Poisson process has rate function lambda(t), t >= 0, which depends on some parameters that must be estimated. We take into account two cases of rate functions: the Weibull and the Goel-Okumoto. We consider models with and without change-points. When the presence of change-points is assumed, we may have the presence of either one, two or three change-points, depending of the data set. The parameters of the rate functions are estimated using a Gibbs sampling algorithm. Results are applied to ozone data provided by the Mexico City monitoring network. In a first instance, we assume that there are no change-points present. Depending on the adjustment of the model, we assume the presence of either one, two or three change-points. Copyright (C) 2009 John Wiley & Sons, Ltd.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Iterative strong coupling of dimensionally heterogeneous models

In this article we address decomposition strategies especially tailored to perform strong coupling of dimensionally heterogeneous models, under the hypothesis that one wants to solve each submodel separately and implement the interaction between subdomains by boundary conditions alone. The novel methodology takes full advantage of the small number of interface unknowns in this kind of problems. Existing algorithms can be viewed as variants of the `natural` staggered algorithm in which each domain transfers function values to the other, and receives fluxes (or forces), and vice versa. This natural algorithm is known as Dirichlet-to-Neumann in the Domain Decomposition literature. Essentially, we propose a framework in which this algorithm is equivalent to applying Gauss-Seidel iterations to a suitably defined (linear or nonlinear) system of equations. It is then immediate to switch to other iterative solvers such as GMRES or other Krylov-based method. which we assess through numerical experiments showing the significant gain that can be achieved. indeed. the benefit is that an extremely flexible, automatic coupling strategy can be developed, which in addition leads to iterative procedures that are parameter-free and rapidly converging. Further, in linear problems they have the finite termination property. Copyright (C) 2009 John Wiley & Sons, Ltd.

Transformed symmetric models

For the first time, we introduce a class of transformed symmetric models to extend the Box and Cox models to more general symmetric models. The new class of models includes all symmetric continuous distributions with a possible non-linear structure for the mean and enables the fitting of a wide range of models to several data types. The proposed methods offer more flexible alternatives to Box-Cox or other existing procedures. We derive a very simple iterative process for fitting these models by maximum likelihood, whereas a direct unconditional maximization would be more difficult. We give simple formulae to estimate the parameter that indexes the transformation of the response variable and the moments of the original dependent variable which generalize previous published results. We discuss inference on the model parameters. The usefulness of the new class of models is illustrated in one application to a real dataset.

On the unification of long-term survival models

In this paper we extend the long-term survival model proposed by Chen et al. [Chen, M.-H., Ibrahim, J.G., Sinha, D., 1999. A new Bayesian model for survival data with a surviving fraction. journal of the American Statistical Association 94, 909-919] via the generating function of a real sequence introduced by Feller [Feller, W., 1968. An Introduction to Probability Theory and its Applications, third ed., vol. 1, Wiley, New York]. A direct consequence of this new formulation is the unification of the long-term survival models proposed by Berkson and Gage [Berkson, J., Gage, R.P., 1952. Survival cure for cancer patients following treatment. journal of the American Statistical Association 47, 501-515] and Chen et al. (see citation above). Also, we show that the long-term survival function formulated in this paper satisfies the proportional hazards property if, and only if, the number of competing causes related to the occurrence of an event of interest follows a Poisson distribution. Furthermore, a more flexible model than the one proposed by Yin and Ibrahim [Yin, G., Ibrahim, J.G., 2005. Cure rate models: A unified approach. The Canadian journal of Statistics 33, 559-570] is introduced and, motivated by Feller`s results, a very useful competing index is defined. (c) 2008 Elsevier B.V. All rights reserved.

Dielectric Dispersion of Water in the Frequency Range from 10 mHz to 30 MHz

We investigate the dielectric dispersion of water, specially in the low-frequency range, by using the impedance spectroscopy technique. The frequency dependencies of the real R and imaginary Z parts of the impedance Could not be explained by means of the Usual description of the dielectric properties of the water as all insulating liquid containing ions. This is due to the incomplete knowledge of the parameters entering in the fundamental equations describing the evolution of the system, and oil the mechanisms regulating the exchange of charge of the cell with the external circuit. We propose a simple description of our experimental data based on the model of Debye, by invoking a dc conductivity of the cell, related to the nonblocking character of the electrodes. A discussion on the electric Circuits able to simulate the cell under investigation, based oil bulk and Surface elements, is also reported. We find that the simple circuit formed by a series of two parallels of resistance and capacitance is able to reproduce the experimental data concerning the real and imaginary part of the electrical impedance of the cell for frequency larger than 1 Hz. According to this description, one of the parallels takes into account the electrical properties of interface between the electrode and water, and the other of the bulk. For frequency lower than 1 Hz, a good agreement with the experimental data is obtained by simulating the electrical properties of the interface by means of the constant phase element.

Topological interactions in warped extra dimensions

Topological interactions will be generated in theories with compact extra dimensions where fermionic chiral zero modes have different localizations. This is the case in many warped extra dimension models where the right-handed top quark is typically localized away from the left-handed one. Using deconstruction techniques, we study the topological interactions in these models. These interactions appear as trilinear and quadrilinear gauge boson couplings in low energy effective theories with three or more sites, as well as in the continuum limit. We derive the form of these interactions for various cases, including examples of Abelian, non-Abelian and product gauge groups of phenomenological interest. The topological interactions provide a window into the more fundamental aspects of these theories and could result in unique signatures at the Large Hadron Collider, some of which we explore.

Stability of the curvature perturbation in dark sectors` mutual interacting models

We consider perturbations in a cosmological model with a small coupling between dark energy and dark matter. We prove that the stability of the curvature perturbation depends on the type of coupling between dark sectors. When the dark energy is of quintessence type, if the coupling is proportional to the dark matter energy density, it will drive the instability in the curvature perturbations: however if the coupling is proportional to the energy density of dark energy, there is room for the stability in the curvature perturbations. When the dark energy is of phantom type, the perturbations are always stable, no matter whether the coupling is proportional to the one or the other energy density. (C) 2008 Elsevier B.V. All rights reserved.

Renormalization of the S parameter in holographic models of electroweak symmetry breaking

We show that the S parameter is not finite in theories of electroweak symmetry breaking in a slice of anti-de Sitter five-dimensional space, with the light fermions localized in the ultraviolet. We compute the one-loop contributions to S from the Higgs sector and show that they are logarithmically dependent on the cutoff of the theory. We discuss the renormalization of S, as well as the implications for bounds from electroweak precision measurements on these models. We argue that, although in principle the choice of renormalization condition could eliminate the S parameter constraint, a more consistent condition would still result in a large and positive S. On the other hand, we show that the dependence on the Higgs mass in S can be entirely eliminated by the renormalization procedure, making it impossible in these theories to extract a Higgs mass bound from electroweak precision constraints.

Magnetic interaction in exchange-biased bilayers: a first-order reversal curve analysis

The distributions of coercivities and magnetic interactions in a set of polycrystalline Ni(0.8)Fe(0.2)/FeMn bilayers have been determined using the first-order reversal curve (FORC) formalism. The thickness of the permalloy (Py) film was fixed at 10 nm (nominal), while that of the FeMn film varied within the range 0-20 nm. The FORC diagrams of each bilayer displayed two clearly distinguishable regions. The main region was generated by Py particles whose coercivities were enhanced in comparison with those in which the FeMn film was absent (sample O). The minor region was produced by Py particles with coercivities similar to or slightly higher than those of particles in the Py film of sample O. Each sample presented two distributions of interaction fields, one for each region, and both were centred slightly below the exchange-bias field, thus indicating a prevalence of magnetizing interactions. These results are consistent with a grain size distribution in the Py layer and the presence of uncompensated antiferromagnetic moments.

First-order-reversal-curve analysis of Pr-Fe-B-based exchange spring magnets

Ribbons of nominal composition (Pr(9.5)Fe(84.5)B(6))(0.96)Cr(0.01)(TiC)(0.03) were produced by arc-melting and melt-spinning the alloys on a Cu wheel. X-ray diffraction (XRD) reveals two main phases, one based upon alpha-Fe and the other upon Pr(2)Fe(14)B. The ribbons show exchange spring behavior with H (c) = 12.5 kOe and (BH)(max) = 13.6 MGOe when these two phases are well coupled. Transmission electron microscopy revealed the coupled behavior is observed when the microstructure consists predominantly of alpha-Fe grains (diameter similar to 100 nm.) surrounded by hard material containing Pr(2)Fe(14)B. The microstructure is discussed in terms of a calculation by Skomski and Coey. A first-order-reversal-curve (FORC) analysis was performed for both a well-coupled sample and a poorly coupled sample. The FORC diagrams show two strong peaks for both the poorly coupled sample and for the well-coupled material. In both cases, the localization of the FORC probability suggests magnetizing interactions between particles. Switching field distributions were calculated and are consistent with the sample microstructure.

Interaction of oligonucleotide-based amphiphilic block copolymers with cell membrane models

Oligonucleotides have unique molecular recognition properties, being involved in biological mechanisms such as cell-surface receptor recognition or gene silencing. For their use in human therapy for drug or gene delivery, the cell membrane remains a barrier, but this can be obviated by grafting a hydrophobic tail to the oligonucleotide. Here we demonstrate that two oligonucleotides, one consisting of 12 guanosine units (G(12)), and the other one consisting of five adenosine and seven guanosine (A(5)G(7)) units, when functionalized with poly(butadiene), namely PB-G(12) and PB-A(5)G(7), can be inserted into Langmuir monolayers of dipalmitoyl phosphatidyl choline (DPPC), which served as a cell membrane model. PB-G(12) and PB-A(5)G(7) were found to affect the DPPC monolayer even at high surface pressures. The effects from PB-G(12) were consistently stronger, particularly in reducing the elasticity of the DPPC monolayers, which may have important biological implications. Multilayers of DPPC and nucleotide-based copolymers could be adsorbed onto solid supports, in the form of Y-type LB films, in which the molecular-level interaction led to lower energies in the vibrational spectra of the nucleotide-based copolymers. This successful deposition of solid films opens the way for devices to be produced which exploit the molecular recognition properties of the nucleotides. (C) 2010 Elsevier Inc. All rights reserved.

Spin-density functional for exchange anisotropic Heisenberg model

Ground-state energies for anti ferromagnetic Heisenberg models with exchange anisotropy are estimated by means of a local-spin approximation made in the context of the density functional theory. Correlation energy is obtained using the non-linear spin-wave theory for homogeneous systems from which the spin functional is built. Although applicable to chains of any size, the results are shown for small number of sites, to exhibit finite-size effects and allow comparison with exact-numerical data from direct diagonalization of small chains. (C) 2009 Elsevier B.V. All rights reserved.

The specificity of frutalin lectin using biomembrane models

Frutalin is a homotetrameric alpha-D-galactose (D-Gal)-binding lectin that activates natural killer cells in vitro and promotes leukocyte migration in vivo. Because lectins are potent lymphocyte stimulators, understanding the interactions that occur between them and cell surfaces can help to the action mechanisms involved in this process. In this paper, we present a detailed investigation of the interactions of frutalin with phospho- and glycolipids using Langmuir monolayers as biomembrane models. The results confirm the specificity of frutalin for D-Gal attached to a biomembrane. Adsorption of frutalin was more efficient for the galactose polar head lipids, in contrast to the one for sulfated galactose, in which a lag time is observed, indicating a rearrangement of the monolayer to incorporate the protein. Regarding ganglioside GM1 monolayers, lower quantities of the protein were adsorbed, probably due to the farther apart position of D-galactose from the interface. Binary mixtures containing galactocerebroside revealed small domains formed at high lipid packing in the presence of frutalin, suggesting that lectin induces the clusterization and the forming of domains in vitro, which may be a form of receptor internalization. This is the first experimental evidence of such lectin effect, and it may be useful to understand the mechanism of action of lectins at the molecular level. (C) 2010 Elsevier B.V. All rights reserved.

Comparative study of liponucleosides in Langmuir monolayers as cell membrane models

Liponucleosides may assist the anchoring of nucleic acid nitrogen bases into biological membranes for tailored nanobiotechnological applications. To this end precise knowledge about the biophysical and chemical details at the membrane surface is required. In this paper, we used Langmuir monolayers as simplified cell membrane models and studied the insertion of five lipidated nucleosides. These molecules varied in the type of the covalently attached lipid group, the nucleobase, and the number of hydrophobic moieties attached to the nucleoside. All five lipidated nucleosides were found to be surface-active and capable of forming stable monolayers. They could also be incorporated into dipalmitoylphosphatidylcholine (DPPC) monolayers, four of which induced expansion in the surface pressure isotherm and a decrease in the surface compression modulus of DPPC. In contrast, one nucleoside possessing three alkyl chain modifications formed very condensed monolayers and induced film condensation and an increase in the compression modulus for the DPPC monolayer, thus reflecting the importance of the ability of the nucleoside molecules to be arranged in a closely packed manner. The implications of these results lie on the possibility of tuning nucleic acid pairing by modifying structural characteristics of the liponucleosides. (C) 2010 Elsevier B.V. All rights reserved.