Stability and ergodicity of piecewise deterministic Markov processes

The main goal of this paper is to establish some equivalence results on stability, recurrence, and ergodicity between a piecewise deterministic Markov process ( PDMP) {X( t)} and an embedded discrete-time Markov chain {Theta(n)} generated by a Markov kernel G that can be explicitly characterized in terms of the three local characteristics of the PDMP, leading to tractable criterion results. First we establish some important results characterizing {Theta(n)} as a sampling of the PDMP {X( t)} and deriving a connection between the probability of the first return time to a set for the discrete-time Markov chains generated by G and the resolvent kernel R of the PDMP. From these results we obtain equivalence results regarding irreducibility, existence of sigma-finite invariant measures, and ( positive) recurrence and ( positive) Harris recurrence between {X( t)} and {Theta(n)}, generalizing the results of [ F. Dufour and O. L. V. Costa, SIAM J. Control Optim., 37 ( 1999), pp. 1483-1502] in several directions. Sufficient conditions in terms of a modified Foster-Lyapunov criterion are also presented to ensure positive Harris recurrence and ergodicity of the PDMP. We illustrate the use of these conditions by showing the ergodicity of a capacity expansion model.

The non-equilibrium nature of culinary evolution

Food is an essential part of civilization, with a scope that ranges from the biological to the economic and cultural levels. Here, we study the statistics of ingredients and recipes taken from Brazilian, British, French and Medieval cookery books. We find universal distributions with scale invariant behaviour. We propose a copy-mutate process to model culinary evolution that fits our empirical data very well. We find a cultural 'founder effect' produced by the non-equilibrium dynamics of the model. Both the invariant and idiosyncratic aspects of culture are accounted for by our model, which may have applications in other kinds of evolutionary processes.

Trapping of strangelets in the geomagnetic field

Strangelets arriving from the interstellar medium are an interesting target for experiments searching for evidence of this hypothetical state of hadronic matter. We entertain the possibility of a trapped strangelet population, quite analogous to ordinary nuclei and electron belts. For a population of strangelets to be trapped by the geomagnetic field, these incoming particles would have to fulfill certain conditions, namely, having magnetic rigidities above the geomagnetic cutoff and below a certain threshold for adiabatic motion to hold. We show in this work that, for fully ionized strangelets, there is a narrow window for stable trapping. An estimate of the stationary population is presented and the dominant loss mechanisms discussed. It is shown that the population would be substantially enhanced with respect to the interstellar medium flux (up to 2 orders of magnitude) due to quasistable trapping.

Confronting dark energy models using galaxy cluster number counts

The mass function of cluster-size halos and their redshift distribution are computed for 12 distinct accelerating cosmological scenarios and confronted to the predictions of the conventional flat Lambda CDM model. The comparison with Lambda CDM is performed by a two-step process. First, we determine the free parameters of all models through a joint analysis involving the latest cosmological data, using supernovae type Ia, the cosmic microwave background shift parameter, and baryon acoustic oscillations. Apart from a braneworld inspired cosmology, it is found that the derived Hubble relation of the remaining models reproduces the Lambda CDM results approximately with the same degree of statistical confidence. Second, in order to attempt to distinguish the different dark energy models from the expectations of Lambda CDM, we analyze the predicted cluster-size halo redshift distribution on the basis of two future cluster surveys: (i) an X-ray survey based on the eROSITA satellite, and (ii) a Sunayev-Zeldovich survey based on the South Pole Telescope. As a result, we find that the predictions of 8 out of 12 dark energy models can be clearly distinguished from the Lambda CDM cosmology, while the predictions of 4 models are statistically equivalent to those of the Lambda CDM model, as far as the expected cluster mass function and redshift distribution are concerned. The present analysis suggests that such a technique appears to be very competitive to independent tests probing the late time evolution of the Universe and the associated dark energy effects.

Comparative Analysis of Activation Phenotype, Proliferation, and IFN-gamma Production by Spleen NK1.1(+) and NK1.1(-) T Cells During Plasmodium chabaudi AS Malaria

The NK1.1 molecule participates in NK, NKT, and T-cell activation, contributing to IFN-gamma production and cytotoxicity. To characterize the early immune response to Plasmodium chabaudi AS, spleen NK1.1(+) and NK1.1(-) T cells were compared in acutely infected C57BL/6 mice. The first parasitemia peak in C57BL/6 mice correlated with increase in CD4(+)NK1.1(+)TCR-alpha beta(+), CD8(+)NK1.1(+)TCR-alpha beta(+), and CD4(+)NK1.1(-)TCR-alpha beta(+) cell numbers per spleen, where a higher increment was observed for NK1.1(+) T cells compared to NK1.1(-) T cells. According to the ability to recognize the CD1d-alpha-GalCer tetramer, CD4(+)NK1.1(+) cells in 7-day infected mice were not predominantly invariant NKT cells. At that time, nearly all NK1.1(+) T cells and around 30% of NK1.1(-) T cells showed an experienced/activated (CD44(HI)CD69(HI)CD122(HI)) cell phenotype, with high expression of Fas and PD-L1 correlating with their low proliferative capacity. Moreover, whereas IFN-gamma production by CD4(+)NK1.1(+) cells peaked at day 4 p.i., the IFN-gamma response of CD4(+)NK1.1(-) cells continued to increase at day 5 of infection. We also observed, at day 7 p.i., 2-fold higher percentages of perforin(+) cells in CD8(+)NK1.1(+) cells compared to CD8(+)NK1.1(-) cells. These results indicate that spleen NK1.1(+) and NK1.1(-) T cells respond to acute P. chabaudi malaria with different kinetics in terms of activation, proliferation, and IFN-gamma production.

FLOWS WITHOUT WANDERING POINTS ON COMPACT CONNECTED SURFACES

Given a compact 2 dimensional manifold M we classify all continuous flows phi without wandering points on M. This classification is performed by finding finitely many pairwise disjoint open phi-invariant subsets {U(1), U(2), ..., U(n)} of M such that U(i=1)(n) (U(i)) over bar = M and each U(i) is either a suspension of an interval exchange transformation, or a maximal open cylinder made up of closed trajectories of phi.

SOME PROPERTIES OF THE MULTIPLICITY SEQUENCE FOR ARBITRARY IDEALS

In this work we prove that the Achilles-Manaresi multiplicity sequence, like the classical Hilbert-Samuel multiplicity, is additive with respect to the exact sequence of modules. We also prove the associativity formula for his mulitplicity sequence. As a consequence, we give new proofs for two results already known. First, the Achilles-Manaresi multiplicity sequence is an invariant up to reduction, a result first proved by Ciuperca. Second, I subset of J is a reduction of (J,M) if and only if c(0)(I(p), M(p)) = c(0)(J(p), M(p)) for all p is an element of Spec(A), a result first proved by Flenner and Manaresi.

AFFINE INTERVAL EXCHANGE TRANSFORMATIONS WITH FLIPS AND WANDERING INTERVALS

There exist uniquely ergodic affine interval exchange transformations of [0,1] with flips which have wandering intervals and are such that the support of the invariant measure is a Cantor set.

CHARACTERIZATION OF EXPONENTIAL DIVERGENCE OF THE KALMAN FILTER FOR TIME-VARYING SYSTEMS

This paper studies semistability of the recursive Kalman filter in the context of linear time-varying (LTV), possibly nondetectable systems with incorrect noise information. Semistability is a key property, as it ensures that the actual estimation error does not diverge exponentially. We explore structural properties of the filter to obtain a necessary and sufficient condition for the filter to be semistable. The condition does not involve limiting gains nor the solution of Riccati equations, as they can be difficult to obtain numerically and may not exist. We also compare semistability with the notions of stability and stability w.r.t. the initial error covariance, and we show that semistability in a sense makes no distinction between persistent and nonpersistent incorrect noise models, as opposed to stability. In the linear time invariant scenario we obtain algebraic, easy to test conditions for semistability and stability, which complement results available in the context of detectable systems. Illustrative examples are included.

Reduction of chaotic particle transport driven by drift waves in sheared flows

Investigations of chaotic particle transport by drift waves propagating in the edge plasma of tokamaks with poloidal zonal flow are described. For large aspect ratio tokamaks, the influence of radial electric field profiles on convective cells and transport barriers, created by the nonlinear interaction between the poloidal flow and resonant waves, is investigated. For equilibria with edge shear flow, particle transport is seen to be reduced when the electric field shear is reversed. The transport reduction is attributed to the robust invariant tori that occur in nontwist Hamiltonian systems. This mechanism is proposed as an explanation for the transport reduction in Tokamak Chauffage Alfven Bresilien [R. M. O. Galvao , Plasma Phys. Controlled Fusion 43, 1181 (2001)] for discharges with a biased electrode at the plasma edge.

Robust relations between CCN and the vertical evolution of cloud drop size distribution in deep convective clouds

In-situ measurements in convective clouds (up to the freezing level) over the Amazon basin show that smoke from deforestation fires prevents clouds from precipitating until they acquire a vertical development of at least 4 km, compared to only 1-2 km in clean clouds. The average cloud depth required for the onset of warm rain increased by similar to 350 m for each additional 100 cloud condensation nuclei per cm(3) at a super-saturation of 0.5% (CCN0.5%). In polluted clouds, the diameter of modal liquid water content grows much slower with cloud depth (at least by a factor of similar to 2), due to the large number of droplets that compete for available water and to the suppressed coalescence processes. Contrary to what other studies have suggested, we did not observe this effect to reach saturation at 3000 or more accumulation mode particles per cm(3). The CCN0.5% concentration was found to be a very good predictor for the cloud depth required for the onset of warm precipitation and other microphysical factors, leaving only a secondary role for the updraft velocities in determining the cloud drop size distributions. The effective radius of the cloud droplets (r(e)) was found to be a quite robust parameter for a given environment and cloud depth, showing only a small effect of partial droplet evaporation from the cloud's mixing with its drier environment. This supports one of the basic assumptions of satellite analysis of cloud microphysical processes: the ability to look at different cloud top heights in the same region and regard their r(e) as if they had been measured inside one well developed cloud. The dependence of r(e) on the adiabatic fraction decreased higher in the clouds, especially for cleaner conditions, and disappeared at r(e)>=similar to 10 mu m. We propose that droplet coalescence, which is at its peak when warm rain is formed in the cloud at r(e)=similar to 10 mu m, continues to be significant during the cloud's mixing with the entrained air, cancelling out the decrease in r(e) due to evaporation.

Measurement of neutral mesons in p plus p collisions at root s=200 GeV and scaling properties of hadron production

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of K(S)(0), omega, eta', and phi mesons in p + p collisions at root s 200 GeV. Measurements of omega and phi production in different decay channels give consistent results. New results for the omega are in agreement with previously published data and extend the measured p(T) coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, n and T, determining the high-p(T) and characterizing the low-p(T) regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.

Nuclear modification factors of phi mesons in d plus Au, Cu plus Cu, and Au plus Au collisions at root s(NN)=200 GeV

The PHENIX experiment at the Relativistic Heavy Ion Collider has performed systematic measurements of phi meson production in the K(+)K(-) decay channel at midrapidity in p + p, d + Au, Cu + Cu, and Au + Au collisions at root s(NN) = 200 GeV. Results are presented on the phi invariant yield and the nuclear modification factor R(AA) for Au + Au and Cu + Cu, and R(dA) for d + Au collisions, studied as a function of transverse momentum (1 < p(T) < 7 GeV/c) and centrality. In central and midcentral Au + Au collisions, the R(AA) of phi exhibits a suppression relative to expectations from binary scaled p + p results. The amount of suppression is smaller than that of the pi(0) and the. in the intermediate p(T) range (2-5 GeV/c), whereas, at higher p(T), the phi, pi(0), and. show similar suppression. The baryon (proton and antiproton) excess observed in central Au + Au collisions at intermediate p(T) is not observed for the phi meson despite the similar masses of the proton and the phi. This suggests that the excess is linked to the number of valence quarks in the hadron rather than its mass. The difference gradually disappears with decreasing centrality, and, for peripheral collisions, the R(AA) values for both particle species are consistent with binary scaling. Cu + Cu collisions show the same yield and suppression as Au + Au collisions for the same number of N(part). The R(dA) of phi shows no evidence for cold nuclear effects within uncertainties.

Enhanced Production of Direct Photons in Au plus Au Collisions at root s(NN)=200 GeV and Implications for the Initial Temperature

The production of e(+)e(-) pairs for m(e+e-) < 0.3 GeV/c(2) and 1< p(T) < 5 GeV/c is measured in p + p and Au + Au collisions at root s(NN) = 200 GeV. An enhanced yield above hadronic sources is observed. Treating the excess as photon internal conversions, the invariant yield of direct photons is deduced. In central Au + Au collisions, the excess of the direct photon yield over p + p is exponential in transverse momentum, with an inverse slope T = 221 +/- 19(stat) +/- 19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T(init) similar to 300-600 MeV at times of similar to 0.6-0.15 fm/c after the collision are in qualitative agreement with the data. Lattice QCD predicts a phase transition to quark gluon plasma at similar to 170 MeV.

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.