A fast Bayesian reconstruction algorithm for emission tomography with entropy prior converging to feasible images

The development and tests of an iterative reconstruction algorithm for emission tomography based on Bayesian statistical concepts are described. The algorithm uses the entropy of the generated image as a prior distribution, can be accelerated by the choice of an exponent, and converges uniformly to feasible images by the choice of one adjustable parameter. A feasible image has been defined as one that is consistent with the initial data (i.e. it is an image that, if truly a source of radiation in a patient, could have generated the initial data by the Poisson process that governs radioactive disintegration). The fundamental ideas of Bayesian reconstruction are discussed, along with the use of an entropy prior with an adjustable contrast parameter, the use of likelihood with data increment parameters as conditional probability, and the development of the new fast maximum a posteriori with entropy (FMAPE) Algorithm by the successive substitution method. It is shown that in the maximum likelihood estimator (MLE) and FMAPE algorithms, the only correct choice of initial image for the iterative procedure in the absence of a priori knowledge about the image configuration is a uniform field.

A numerical model for the gamma-ray emission of the microquasar LS 5039

The possible association between the microquasar LS 5039 and the EGRET source 3EG J1824-1514 suggests that microquasars could also be sources of high energy gamma-rays. In this paper, we explore, with a detailed numerical model, if this system can produce the emission detected by EGRET (>100 MeV) through inverse Compton (IC) scattering. Our numerical approach considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting with both the radiation and the magnetic fields, taking into account the Thomson and Klein-Nishina regimes of interaction. The computed spectrum reproduces the observed spectral characteristics at very high energy.

High-energy gamma-ray emission from microquasars: LS 5039 and LS I+61303

The possible associations between the microquasars LS 5039 and LS I +61 303 and the EGRET sources 3EG J1824-1514 and 3EG J0241+6103 suggest that microquasars could also be sources of high-energy gamma-rays. In this work, we present a detailed numerical inverse Compton (IC) model, based on a microquasar scenario, that reproduces the high-energy gamma-ray spectra and variability observed by EGRET for the mentioned sources. Our model considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet that interact through IC scattering with both the radiation and the magnetic fields.

Exploring the high frequency emission of radio loud X-ray binaries

We report millimetre-wave continuum observations of the X-ray binaries Cygnus X-3, SS 433, LSI+61 303, Cygnus X-1 and GRS 1915+105. The observations were carried out with the IRAM 30 m-antenna at 250 GHz (1.25 mm) from 1998 March 14 to March 20. These millimetre measurements are complemented with centimetre observations from the Ryle Telescope, at 15 GHz (2.0 cm) and from the Green Bank Interferometer at 2.25 and 8.3 GHz (13 and 3.6 cm). Both Cygnus X-3 and SS 433 underwent moderate flaring events during our observations, whose main spectral evolution properties are described and interpreted. A significant spectral steepening was observed in both sources during the flare decay, that is likely to be caused by adiabatic expansion, inverse Compton and synchrotron losses. Finally, we also report 250 GHz upper limits for three additional undetected X-ray binary stars: LSI+65 010, LSI+61 235 and X Per.

Search for very high energy gamma-ray emission from pulsar-pulsar wind nebula systems with the MAGIC telescope

The MAGIC collaboration has searched for high-energy gamma-ray emission of some of the most promising pulsar candidates above an energy threshold of 50 GeV, an energy not reachable up to now by other ground-based instruments. Neither pulsed nor steady gamma-ray emission has been observed at energies of 100 GeV from the classical radio pulsars PSR J0205+6449 and PSR J2229+6114 (and their nebulae 3C58 and Boomerang, respectively) and the millisecond pulsar PSR J0218+4232. Here, we present the flux upper limits for these sources and discuss their implications in the context of current model predictions.

Statistical analysis of maximum likelihood estimator images of human brain FDG PET studies

The work presented evaluates the statistical characteristics of regional bias and expected error in reconstructions of real positron emission tomography (PET) data of human brain fluoro-deoxiglucose (FDG) studies carried out by the maximum likelihood estimator (MLE) method with a robust stopping rule, and compares them with the results of filtered backprojection (FBP) reconstructions and with the method of sieves. The task of evaluating radioisotope uptake in regions-of-interest (ROIs) is investigated. An assessment of bias and variance in uptake measurements is carried out with simulated data. Then, by using three different transition matrices with different degrees of accuracy and a components of variance model for statistical analysis, it is shown that the characteristics obtained from real human FDG brain data are consistent with the results of the simulation studies.

Faint arc-minute extended radio emission around Cygnus X-3

Aims.We revisit the vicinity of the microquasar Cygnus X-3 at radio wavelengths. We aim to improve our previous search for possible associated extended radio features/hot spots in the position angle of the Cygnus X-3 relativistic jets focusing on shorter angular scales than previously explored. Methods.Our work is mostly based on analyzing modern survey and archive radio data, mainly including observations carried out with the Very Large Array and the Ryle Telescopes. We also used deep near-infrared images that we obtained in 2005. Results.We present new radio maps of the Cygnus X-3 field computed after combining multi-configuration Very Large Array archive data at 6 cm and different observing runs at 2 cm with the Ryle Telescope. These are probably among the deepest radio images of Cygnus X-3 reported to date at cm wavelengths. Both interferometers reveal an extended radio feature within a few arc-minutes of the microquasar position, thus making our detection more credible. Moreover, this extended emission is possibly non-thermal, although this point still needs confirmation. Its physical connection with the microquasar is tentatively considered under different physical scenarios. We also report on the serendipitous discovery of a likely Fanaroff-Riley type II radio galaxy only away from Cygnus X-3.

Leptonic secondary emission in a hadronic microquasar model

Context.It has been proposed that the origin of the very high-energy photons emitted from high-mass X-ray binaries with jet-like features, so-called microquasars (MQs), is related to hadronic interactions between relativistic protons in the jet and cold protons of the stellar wind. Leptonic secondary emission should be calculated in a complete hadronic model that includes the effects of pairs from charged pion decays inside the jets and the emission from pairs generated by gamma-ray absorption in the photosphere of the system. Aims.We aim at predicting the broadband spectrum from a general hadronic microquasar model, taking into account the emission from secondaries created by charged pion decay inside the jet. Methods.The particle energy distribution for secondary leptons injected along the jets is consistently derived taking the energy losses into account. The spectral energy distribution resulting from these leptons is calculated after assuming different values of the magnetic field inside the jets. We also compute the spectrum of the gamma-rays produced by neutral pion-decay and processed by electromagnetic cascades under the stellar photon field. Results.We show that the secondary emission can dominate the spectral energy distribution at low energies (~1 MeV). At high energies, the production spectrum can be significantly distorted by the effect of electromagnetic cascades. These effects are phase-dependent, and some variability modulated by the orbital period is predicted.

Exploring the connection between the stellar wind and the non-thermal emission in LS 5039

Context.LS 5039 has been observed with several X-ray instruments so far showing quite steady emission in the long term and no signatures of accretion disk. The source also presents X-ray variability at orbital timescales in flux and photon index. The system harbors an O-type main sequence star with moderate mass-loss. At present, the link between the X-rays and the stellar wind is unclear. Aims.We study the X-ray fluxes, spectra, and absorption properties of LS 5039 at apastron and periastron passages during an epoch of enhanced stellar mass-loss, and the long term evolution of the latter in connection with the X-ray fluxes. Methods.New XMM-Newton observations were performed around periastron and apastron passages in September 2005, when the stellar wind activity was apparently higher. April 2005 Chandra observations on LS 5039 were revisited. Moreover, a compilation of H EW data obtained since 1992, from which the stellar mass-loss evolution can be approximately inferred, was carried out. Results.XMM-Newton observations show higher and harder emission around apastron than around periastron. No signatures of thermal emission or a reflection iron line indicating the presence of an accretion disk are found in the spectrum, and the hydrogen column density () is compatible with being the same in both observations and consistent with the interstellar value. 2005 Chandra observations show a hard X-ray spectrum, and possibly high fluxes, although pileup effects preclude conclusive results from being obtained. The H EW shows yearly variations of 10%, and does not seem to be correlated with X-ray fluxes obtained at similar phases, unlike what is expected in the wind accretion scenario. Conclusions.2005 XMM-Newton and Chandra observations are consistent with 2003 RXTE/PCA results, namely moderate flux and spectral variability at different orbital phases. The constancy of the seems to imply that either the X-ray emitter is located at 1012 cm from the compact object, or the density in the system is 3 to 27 times smaller than that predicted by a spherical symmetric wind model. We suggest that the multiwavelength non-thermal emission of LS 5039 is related to the observed extended radio jets and is unlikely to be produced inside the binary system.

Variable Very-High-Energy Gamma-Ray Emission from the Microquasar LS I +61 303

Microquasars are binary star systems with relativistic radio-emitting jets. They are potential sources of cosmic rays and can be used to elucidate the physics of relativistic jets. We report the detection of variable gamma-ray emission above 100 gigaelectron volts from the microquasar LS I 61 + 303. Six orbital cycles were recorded. Several detections occur at a similar orbital phase, which suggests that the emission is periodic. The strongest gamma-ray emission is not observed when the two stars are closest to one another, implying a strong orbital modulation of the emission or absorption processes.

Complementarity between human capital and trade in regional technological progress

The effect of openness and trade orientation on economic growth remains a highly contentious issue in the literature. Trade facilitates the spread of knowledge and the adoption of more advanced and efficient technologies, which hastens total factor productivity (TFP) growth and, hence, per capita income. New technologies that spread through trade require a sufficiently skilled labour force to adapt them to the domestic productive environment. Thus, openness and human capital accumulation will lead to TFP growth and the greater the complementarity between both variables, the higher the TFP growth. This paper discusses the implications of these assumptions and tests their empirical validity, using a pool of data for manufacturing industry in Spanish regions in a period in which both the stock of human capital and openness experienced a notable increase.

Regional economic dynamics and convergence in the european union

Deepening in the European Union (EU) integration process has enhanced the question of economic disparities at a regional level. Theconvergence process observed until the late seventies was exhausted onwards incoincidence with important changes in the economic activity. The paper showshow these factors would have provoked a regional differenciated response that,despite being important, would have not strengthened the decrease in regionalinequalities. We use an alternative and (in our opinion) richer approach to thetraditional convergence analysis, where the evolution of the whole regionaldistribution is what matters and not that of a representative economy. Moreover,when analysing inequalities among regional economies, the geographical spaceacquire an outstanding role. Hence, we apply spatial association tests and relatethem to the convergence analysis

Regional structure of wages and external economies in Spain

Regional data on wages for the Spanish economy show that workers who live in developed regions earn more than workers in other regions.Literature on external economies provides a possible explanation of why firms do not move from these regions to others where wages are lower. Previous studies for the Spanish case use aggregated sectoral data to explain in terms of external economies why average wages are different across regions. The originalcontribution of this paper consists of using individual data to detect the existenceand nature of external economies as an explanatory cause of territorial wagedifferences. With this aim, we have used individual data from the EPF 1990-91(INE). This information permits us to control the influence of individual and jobcharacteristics on wages to, first, detect the existence of external economies and,second, to test alternative explanations of their presence. The empirical evidenceobtained confirms the relevance of territorial external economies and their influence on wages, as a result of improvements in the productive efficiency of the firm. In concrete terms, the more relevant external economies are associatedwith the regional human capital stock and geographical productive specialisation

Potencialidad de la modelización state-space y el filtro de Kalman para el anàlisis regional : una aplicación para el índice de actividad industrial

El déficit existente a nuestro país con respecto a la disponibilidad de indicadores cuantitativos con los que llevar a término un análisis coyuntural de la actividad industrial regional ha abierto un debate centrado en el estudio de cuál es la metodología más adecuada para elaborar indicadores de estas características. Dentro de este marco, en este trabajo se presentan las principales conclusiones obtenidas en anteriores estudios (Clar, et. al., 1997a, 1997b y 1998) sobre la idoneidad de extender las metodologías que actualmente se están aplicando a las regiones españolas para elaborar indicadores de la actividad industrial mediante métodos indirectos. Estas conclusiones llevan a plantear una estrategia distinta a las que actualmente se vienen aplicando. En concreto, se propone (siguiendo a Israilevich y Kuttner, 1993) un modelo de variables latentes para estimar el indicador de la producción industrial regional. Este tipo de modelo puede especificarse en términos de un modelo statespace y estimarse mediante el filtro de Kalman. Para validar la metodología propuesta se estiman unos indicadores de acuerdo con ella para tres de las cuatro regiones españolas que disponen d¿un Índice de Producción Industrial (IPI) elaborado mediante el método directo (Andalucía, Asturias y el País Vasco) y se comparan con los IPIs publicados (oficiales). Los resultados obtenidos muestran el buen comportamiento de l¿estrategia propuesta, abriendo así una línea de trabajo con la que subsanar el déficit al que se hacía referencia anteriormente

Algunes reflexions sobre la construcció d'indicadors indirectes pel seguiment de l'activitat industrial regional

En un gran nombre d'economies, l'evolució de la producció industrial s'analitza a partir de la informació sobre el Producte Industrial Brut i/o el Valor Afegit Brut que proporcionen les Comptabilitats Nacionals. A Espanya, la utilització d'aquestes dades presenta el problema que no estan disponibles tan ràpidament com seria desitjable. En conseqüència, no és possible realitzar un seguiment a curt termini de l'activitat industrial a partir dels mateixos. Per a solucionar aquest problema, l'Institut Nacional d'Estadística elabora un Índex de Producció Industrial mensual a partir de la informació obtinguda a través d'una enquesta dirigida a una mostra representativa de les empreses espanyoles. No obstant això, a nivell regional, les dificultats per a realitzar un seguiment de l'activitat industrial són majors a causa de l'escassesa d'informació estadística. Durant els últims anys, diferents institucions públiques i privades han començat a elaborar indicadors d'activitat per a algunes regions espanyoles, encara que a partir de metodologies no homogènies, de manera que aquests índexs no són directament comparables. Per a corregir aquesta situació, en diferents fòrums s'ha proposat emprar la metodologia utilitzada per l'Institut d'Estadística de Catalunya (IEC) per a la comunitat catalana com alternativa per a aquelles comunitats espanyoles que no disposen d'un indicador de l'activitat industrial, atès que per a Catalunya resulta una metodologia adequada. En aquest treball s'estudia la idoneïtat d'estendre aquesta metodologia a la resta de regions espanyoles. Per a això, es construeixen uns indicadors d'acord amb la metodologia del IEC i es comparen amb els índexs regionals obtinguts per mètodes directes per a tres de les quatre regions que existeixen: Andalusia, Astúries i Euskadi