Whole-History Rating: A Bayesian Rating System for Players of Time-Varying Strength

International audienceWhole-History Rating (WHR) is a new method to estimate the time-varying strengths of players involved in paired comparisons. Like many variations of the Elo rating system, the whole-history approach is based on the dynamic Bradley-Terry model. But, instead of using incremental approximations, WHR directly computes the exact maximum a posteriori over the whole rating history of all players. This additional accuracy comes at a higher computational cost than traditional methods, but computation is still fast enough to be easily applied in real time to large-scale game servers (a new game is added in less than 0.001 second). Experiments demonstrate that, in comparison to Elo, Glicko, TrueSkill, and decayed-history algorithms, WHR produces better predictions

A perturbation particle method for stability studies of stellar systems

In order to investigate the stability of stellar systems, we discuss a numerical method that uses an analytic distribution function to describe a stellar system in equilibrium and "perturbation particles" to represent departures from that equilibrium state. Thus, all the particles are used only to represent the perturbation, and statistical fluctuations due to the finite number of particles are much less severe than in full N-body codes. We provide a general description of the method, recipes for particular aspects of its implementation, and an example of its application to a simple model with known analytical solution.Asociación Argentina de Astronomí

Calculation of Spectral Darkening and Visibility Functions for Solar Oscillations

Calculations of spectral darkening and visibility functions for the brightness oscillations of the Sun resulting from global solar oscillations are presented. This has been done for a broad range of the visible and infrared continuum spectrum. The procedure for the calculations of these functions includes the numerical computation of depth-dependent derivatives of the opacity caused by p modes in the photosphere. A radiative-transport code was used for this purpose to get the disturbances of the opacities from temperature and density fluctuations. The visibility and darkening functions are obtained for adiabatic oscillations under the assumption that the temperature disturbances are proportional to the undisturbed temperature of the photosphere. The latter assumption is the only way to explore any opacity effects since the eigenfunctions of p-mode oscillations have not been obtained so far. This investigation reveals that opacity effects have to be taken into account because they dominate the violet and infrared part of the spectrum. Because of this dominance, the visibility functions are negative for those parts of the spectrum. Furthermore, the darkening functions show a wavelength-dependent change of sign for some wavelengths owing to these opacity effects. However, the visibility and darkening functions under the assumptions used contradict the observations of global p-mode oscillations, but it is beyond doubt that the opacity effects influence the brightness fluctuations of the Sun resulting from global oscillations

Multidimensional relativistic MHD simulations of Pulsar Wind Nebulae: dynamics and emission

Pulsar Wind Nebulae, and the Crab nebula in particular, are the best cosmic laboratories to investigate the dynamics of magnetized relativistic outflows and particle acceleration up to PeV energies. Multidimensional MHD modeling by means of numerical simulations has been very successful at reproducing, to the very finest details, the innermost structure of these synchrotron emitting nebulae, as observed in the X-rays. Therefore, the comparison between the simulated source and observations can be used as a powerful diagnostic tool to probe the physical conditions in pulsar winds, like their composition, magnetization, and degree of anisotropy. However, in spite of the wealth of observations and of the accuracy of current MHD models, the precise mechanisms for magnetic field dissipation and for the acceleration of the non-thermal emitting particles are mysteries still puzzling theorists to date. Here we review the methodologies of the computational approach to the modeling of Pulsar Wind Nebulae, discussing the most relevant results and the recent progresses achieved in this fascinating field of high-energy astrophysics.Comment: 29 pages review, preliminary version. To appear in the book "Modelling Nebulae" edited by D. Torres for Springer, based on the invited contributions to the workshop held in Sant Cugat (Barcelona), June 14-17, 201

The thermal and kinematic Sunyaev-Zel'dovich effects revisited

This paper shows that a simple convolution integral expression based on the mean value of the isotropic frequency distribution corresponding to photon scattering off electrons leads to useful analytical expressions describing the thermal Sunyaev-Zel'dovich effect. The approach, to first order in the Compton parameter is able to reproduce the Kompaneets equation describing the effect. Second order effects in the parameter $z=\frac{kT_{e}}{mc^{2}}$ induce a slight increase in the crossover frequency.Comment: 7 pages, 2 figure

Axion-like particle effects on the polarization of cosmic high-energy gamma sources

Various satellite-borne missions are being planned whose goal is to measure the polarization of a large number of gamma-ray bursts (GRBs). We show that the polarization pattern predicted by current models of GRB emission can be drastically modified by the existence of very light axion-like particles (ALPs), which are present in many extensions of the Standard Model of particle physics. Basically, the propagation of photons emitted by a GRB through cosmic magnetic fields with a domain-like structure induces photon-ALP mixing, which is expected to produce a strong modification of the original photon polarization. Because of the random orientation of the magnetic field in each domain, this effect strongly depends on the orientation of the photon line of sight. As a consequence, photon-ALP conversion considerably broadens the original polarization distribution. Searching for such a peculiar feature through future high-statistics polarimetric measurements is therefore a new opportunity to discover very light ALPs.Comment: Final version (21 pages, 8 eps figures). Matches the version published on JCAP. Added a Section on the effects of cosmic expansion on photon-ALP conversions. Figures modified to take into account this effect. References updated. Conclusions unchanged

Signatures of photon and axion-like particle mixing in the gamma-ray burst jet

Photons couple to Axion-Like Particles (ALPs) or more generally to any pseudo Nambu-Goldstone boson in the presence of an external electromagnetic field. Mixing between photons and ALPs in the strong magnetic field of a Gamma-Ray Burst (GRB) jet during the prompt emission phase can leave observable imprints on the gamma-ray polarization and spectrum. Mixing in the intergalactic medium is not expected to modify these signatures for ALP mass > 10^(-14) eV and/or for < nG magnetic field. We show that the depletion of photons due to conversion to ALPs changes the linear degree of polarization from the values predicted by the synchrotron model of gamma ray emission. We also show that when the magnetic field orientation in the propagation region is perpendicular to the field orientation in the production region, the observed synchrotron spectrum becomes steeper than the theoretical prediction and as detected in a sizable fraction of GRB sample. Detection of the correlated polarization and spectral signatures from these steep-spectrum GRBs by gamma-ray polarimeters can be a very powerful probe to discover ALPs. Measurement of gamma-ray polarization from GRBs in general, with high statistics, can also be useful to search for ALPs.Comment: 17 pages, 3 figures. Accepted for publication in JCAP with minor change

The effect of two-temperature post-shock accretion flow on the linear polarization pulse in magnetic cataclysmic variables

The temperatures of electrons and ions in the post-shock accretion region of a magnetic cataclysmic variable (mCV) will be equal at sufficiently high mass flow rates or for sufficiently weak magnetic fields. At lower mass flow rates or in stronger magnetic fields, efficient cyclotron cooling will cool the electrons faster than the electrons can cool the ions and a two-temperature flow will result. Here we investigate the differences in polarized radiation expected from mCV post-shock accretion columns modeled with one- and two-temperature hydrodynamics. In an mCV model with one accretion region, a magnetic field >~30 MG and a specific mass flow rate of ~0.5 g/cm/cm/s, along with a relatively generic geometric orientation of the system, we find that in the ultraviolet either a single linear polarization pulse per binary orbit or two pulses per binary orbit can be expected, depending on the accretion column hydrodynamic structure (one- or two-temperature) modeled. Under conditions where the physical flow is two-temperature, one pulse per orbit is predicted from a single accretion region where a one-temperature model predicts two pulses. The intensity light curves show similar pulse behavior but there is very little difference between the circular polarization predictions of one- and two-temperature models. Such discrepancies indicate that it is important to model some aspect of two-temperature flow in indirect imaging procedures, like Stokes imaging, especially at the edges of extended accretion regions, were the specific mass flow is low, and especially for ultraviolet data.Comment: Accepted for publication in Astrophysics & Space Scienc

New Insights into White-Light Flare Emission from Radiative-Hydrodynamic Modeling of a Chromospheric Condensation

(abridged) The heating mechanism at high densities during M dwarf flares is poorly understood. Spectra of M dwarf flares in the optical and near-ultraviolet wavelength regimes have revealed three continuum components during the impulsive phase: 1) an energetically dominant blackbody component with a color temperature of T $\sim$ 10,000 K in the blue-optical, 2) a smaller amount of Balmer continuum emission in the near-ultraviolet at lambda $<$ 3646 Angstroms and 3) an apparent pseudo-continuum of blended high-order Balmer lines. These properties are not reproduced by models that employ a typical "solar-type" flare heating level in nonthermal electrons, and therefore our understanding of these spectra is limited to a phenomenological interpretation. We present a new 1D radiative-hydrodynamic model of an M dwarf flare from precipitating nonthermal electrons with a large energy flux of $10^{13}$ erg cm$^{-2}$ s$^{-1}$. The simulation produces bright continuum emission from a dense, hot chromospheric condensation. For the first time, the observed color temperature and Balmer jump ratio are produced self-consistently in a radiative-hydrodynamic flare model. We find that a T $\sim$ 10,000 K blackbody-like continuum component and a small Balmer jump ratio result from optically thick Balmer and Paschen recombination radiation, and thus the properties of the flux spectrum are caused by blue light escaping over a larger physical depth range compared to red and near-ultraviolet light. To model the near-ultraviolet pseudo-continuum previously attributed to overlapping Balmer lines, we include the extra Balmer continuum opacity from Landau-Zener transitions that result from merged, high order energy levels of hydrogen in a dense, partially ionized atmosphere. This reveals a new diagnostic of ambient charge density in the densest regions of the atmosphere that are heated during dMe and solar flares.Comment: 50 pages, 2 tables, 13 figures. Accepted for publication in the Solar Physics Topical Issue, "Solar and Stellar Flares". Version 2 (June 22, 2015): updated to include comments by Guest Editor. The final publication is available at Springer via http://dx.doi.org/10.1007/s11207-015-0708-

Statistics of the gravitational force in various dimensions of space: from Gaussian to Levy laws

We discuss the distribution of the gravitational force created by a Poissonian distribution of field sources (stars, galaxies,...) in different dimensions of space d. In d=3, it is given by a Levy law called the Holtsmark distribution. It presents an algebraic tail for large fluctuations due to the contribution of the nearest neighbor. In d=2, it is given by a marginal Gaussian distribution intermediate between Gaussian and Levy laws. In d=1, it is exactly given by the Bernouilli distribution (for any particle number N) which becomes Gaussian for N>>1. Therefore, the dimension d=2 is critical regarding the statistics of the gravitational force. We generalize these results for inhomogeneous systems with arbitrary power-law density profile and arbitrary power-law force in a d-dimensional universe