No Arabic abstract
We propose to model the image differentials of astrophysical source maps by Students t-distribution and to use them in the Bayesian source separation method as priors. We introduce an efficient Markov Chain Monte Carlo (MCMC) sampling scheme to unmix the astrophysical sources and describe the derivation details. In this scheme, we use the Langevin stochastic equation for transitions, which enables parallel drawing of random samples from the posterior, and reduces the computation time significantly (by two orders of magnitude). In addition, Students t-distribution parameters are updated throughout the iterations. The results on astrophysical source separation are assessed with two performance criteria defined in the pixel and the frequency domains.
A long-standing problem of astrophysical research is how to simultaneously obtain spectra of thousands of sources randomly positioned in the field of view of a telescope. Digital Micromirror Devices, used as optical switches, provide a most powerful solution allowing to design a new generation of instruments with unprecedented capabilities. We illustrate the key factors (opto-mechanical, cryo-thermal, cosmic radiation environment,...) that constrain the design of DMD-based multi-object spectrographs, with particular emphasis on the IR spectroscopic channel onboard the EUCLID mission, currently considered by the European Space Agency for a 2017 launch date.
Time domain astronomy has come of age with astronomers now able to monitor the sky at high cadence both across the electromagnetic spectrum and using neutrinos and gravitational waves. The advent of new observing facilities permits new science, but the ever increasing throughput of facilities demands efficient communication of coincident detections and better subsequent coordination among the scientific community so as to turn detections into scientific discoveries. To discuss the revolution occurring in our ability to monitor the Universe and the challenges it brings, on 2012 April 25-26 a group of scientists from observational and theoretical teams studying transients met with representatives of the major international transient observing facilities at the Kavli Royal Society International Centre, UK. This immediately followed the Royal Society Discussion meeting New windows on transients across the Universe held in London. Here we present a summary of the Kavli meeting at which the participants discussed the science goals common to the transient astronomy community and analysed how to better meet the challenges ahead as ever more powerful observational facilities come on stream.
We describe the current status of CATS (astrophysical CATalogs Support system), a publicly accessible tool maintained at Special Astrophysical Observatory of the Russian Academy of Sciences (SAO RAS) (http://cats.sao.ru) allowing one to search hundreds of catalogs of astronomical objects discovered all along the electromagnetic spectrum. Our emphasis is mainly on catalogs of radio continuum sources observed from 10 MHz to 245 GHz, and secondly on catalogs of objects such as radio and active stars, X-ray binaries, planetary nebulae, HII regions, supernova remnants, pulsars, nearby and radio galaxies, AGN and quasars. CATS also includes the catalogs from the largest extragalactic surveys with non-radio waves. In 2008 CATS comprised a total of about 10e9 records from over 400 catalogs in the radio, IR, optical and X-ray windows, including most source catalogs deriving from observations with the Russian radio telescope RATAN-600. CATS offers several search tools through different ways of access, e.g. via web interface and e-mail. Since its creation in 1997 CATS has managed about 10,000 requests. Currently CATS is used by external users about 1500 times per day and since its opening to the public in 1997 has received about 4000 requests for its selection and matching tasks.
Analysis of cosmic microwave background (CMB) datasets typically requires some filtering of the raw time-ordered data. Filtering is frequently used to minimize the impact of low frequency noise, atmospheric contributions and/or scan synchronous signals on the resulting maps. In this work we explicitly construct a general filtering operator, which can unambiguously remove any set of unwanted modes in the data, and then amend the map-making procedure in order to incorporate and correct for it. We show that such an approach is mathematically equivalent to the solution of a problem in which the sky signal and unwanted modes are estimated simultaneously and the latter are marginalized over. We investigate the conditions under which this amended map-making procedure can render an unbiased estimate of the sky signal in realistic circumstances. We then study the effects of time-domain filtering on the noise correlation structure in the map domain, as well as impact it may have on the performance of the popular pseudo-spectrum estimators. We conclude that although maps produced by the proposed estimators arguably provide the most faithful representation of the sky possible given the data, they may not straightforwardly lead to the best constraints on the power spectra of the underlying sky signal and special care may need to be taken to ensure this is the case. By contrast, simplified map-makers which do not explicitly correct for time-domain filtering, but leave it to subsequent steps in the data analysis, may perform equally well and be easier and faster to implement. We focus on polarization-sensitive measurements targeting the B-mode component of the CMB signal and apply the proposed methods to realistic simulations based on characteristics of an actual CMB polarization experiment, POLARBEAR.
For large scale on-line inference problems the update strategy is critical for performance. We derive an adaptive scan Gibbs sampler that optimizes the update frequency by selecting an optimum mini-batch size. We demonstrate performance of our adaptive batch-size Gibbs sampler by comparing it against the collapsed Gibbs sampler for Bayesian Lasso, Dirichlet Process Mixture Models (DPMM) and Latent Dirichlet Allocation (LDA) graphical models.