Showing papers in "arXiv: Instrumentation and Methods for Astrophysics in 2012"
TL;DR: This document introduces a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010).
Abstract: We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to $\sim N^2$ for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation and API. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at this http URL under the MIT License.
5,293 citations
Posted Content•
TL;DR: The Astrophysics Source Code Library (ASCL) as discussed by the authors is a free on-line registry for source codes of interest to astronomers and astrophysicists, which can be accessed at this http URL.
Abstract: The Astrophysics Source Code Library (ASCL), founded in 1999, is a free on-line registry for source codes of interest to astronomers and astrophysicists. The library is housed on the discussion forum for Astronomy Picture of the Day (APOD) and can be accessed at this http URL. The ASCL has a comprehensive listing that covers a significant number of the astrophysics source codes used to generate results published in or submitted to refereed journals and continues to grow. The ASCL currently has entries for over 500 codes; its records are citable and are indexed by ADS. The editors of the ASCL and members of its Advisory Committee were on hand at a demonstration table in the ADASS poster room to present the ASCL, accept code submissions, show how the ASCL is starting to be used by the astrophysics community, and take questions on and suggestions for improving the resource.
1,789 citations
TL;DR: In this paper, the Sloan Digital Sky Survey III (SDSS-III) presented the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS), along with the data presented in previous data releases.
Abstract: The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff -0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2).
The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.
1,443 citations
TL;DR: In this paper, the design and performance of the multi-object fiber spectrographs for the Sloan Digital Sky Survey (SDSS) and their upgrade for the Baryon Oscillation Spectroscopic Survey (BOSS) were presented.
Abstract: We present the design and performance of the multi-object fiber spectrographs for the Sloan Digital Sky Survey (SDSS) and their upgrade for the Baryon Oscillation Spectroscopic Survey (BOSS). Originally commissioned in Fall 1999 on the 2.5-m aperture Sloan Telescope at Apache Point Observatory, the spectrographs produced more than 1.5 million spectra for the SDSS and SDSS-II surveys, enabling a wide variety of Galactic and extra-galactic science including the first observation of baryon acoustic oscillations in 2005. The spectrographs were upgraded in 2009 and are currently in use for BOSS, the flagship survey of the third-generation SDSS-III project. BOSS will measure redshifts of 1.35 million massive galaxies to redshift 0.7 and Lyman-alpha absorption of 160,000 high redshift quasars over 10,000 square degrees of sky, making percent level measurements of the absolute cosmic distance scale of the Universe and placing tight constraints on the equation of state of dark energy.
The twin multi-object fiber spectrographs utilize a simple optical layout with reflective collimators, gratings, all-refractive cameras, and state-of-the-art CCD detectors to produce hundreds of spectra simultaneously in two channels over a bandpass covering the near ultraviolet to the near infrared, with a resolving power R = \lambda/FWHM ~ 2000. Building on proven heritage, the spectrographs were upgraded for BOSS with volume-phase holographic gratings and modern CCD detectors, improving the peak throughput by nearly a factor of two, extending the bandpass to cover 360 < \lambda < 1000 nm, and increasing the number of fibers from 640 to 1000 per exposure. In this paper we describe the original SDSS spectrograph design and the upgrades implemented for BOSS, and document the predicted and measured performances.
977 citations
TL;DR: The Murchison Widefield Array (MWA) as mentioned in this paper is one of three Square Kilometre Array Precursor telescopes and is located at the MUR-astronomy Observatory in Western Australia, a location chosen for its extremely low levels of radio frequency interference.
Abstract: The Murchison Widefield Array (MWA) is one of three Square Kilometre Array Precursor telescopes and is located at the Murchison Radio-astronomy Observatory in the Murchison Shire of the mid-west of Western Australia, a location chosen for its extremely low levels of radio frequency interference. The MWA operates at low radio frequencies, 80-300 MHz, with a processed bandwidth of 30.72 MHz for both linear polarisations, and consists of 128 aperture arrays (known as tiles) distributed over a ~3 km diameter area. Novel hybrid hardware/software correlation and a real-time imaging and calibration systems comprise the MWA signal processing backend. In this paper the as-built MWA is described both at a system and sub-system level, the expected performance of the array is presented, and the science goals of the instrument are summarised.
808 citations
Posted Content•
TL;DR: The LISA Consortium as discussed by the authors proposed a 4-year mission in response to ESA's call for missions for L3, which provides the closest ever view of the infant universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales.
Abstract: Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km.
LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.
582 citations
TL;DR: In this article, a Bayesian Maximum A Posteriori (MAP) approach is presented where a subset of highly correlated and quiet stars is used to generate a cotrending basis vector set which is in turn used to establish a range of "reasonable" robust fit parameters.
Abstract: With the unprecedented photometric precision of the Kepler Spacecraft, significant systematic and stochastic errors on transit signal levels are observable in the Kepler photometric data. These errors, which include discontinuities, outliers, systematic trends and other instrumental signatures, obscure astrophysical signals. The Presearch Data Conditioning (PDC) module of the Kepler data analysis pipeline tries to remove these errors while preserving planet transits and other astrophysically interesting signals. The completely new noise and stellar variability regime observed in Kepler data poses a significant problem to standard cotrending methods such as SYSREM and TFA. Variable stars are often of particular astrophysical interest so the preservation of their signals is of significant importance to the astrophysical community. We present a Bayesian Maximum A Posteriori (MAP) approach where a subset of highly correlated and quiet stars is used to generate a cotrending basis vector set which is in turn used to establish a range of "reasonable" robust fit parameters. These robust fit parameters are then used to generate a Bayesian Prior and a Bayesian Posterior Probability Distribution Function (PDF) which when maximized finds the best fit that simultaneously removes systematic effects while reducing the signal distortion and noise injection which commonly afflicts simple least-squares (LS) fitting. A numerical and empirical approach is taken where the Bayesian Prior PDFs are generated from fits to the light curve distributions themselves.
520 citations
TL;DR: In this article, a new method to achieve point spread function (PSF) subtractions for high contrast imaging using principal component analysis (PCA) that is applicable to both point sources or extended objects (disks) is described.
Abstract: We describe a new method to achieve point spread function (PSF) subtractions for high- contrast imaging using Principal Component Analysis (PCA) that is applicable to both point sources or extended objects (disks). Assuming a library of reference PSFs, a Karhunen-Lo`eve transform of theses references is used to create an orthogonal basis of eigenimages, on which the science target is projected. For detection this approach provides comparable suppression to the Locally Optimized Combination of Images (LOCI) algorithm, albeit with increased robustness to the algorithm parameters and speed enhancement. For characterization of detected sources the method enables forward modeling of astrophysical sources. This alleviates the biases in the astrometry and photometry of discovered faint sources, which are usually associated with LOCI- based PSF subtractions schemes. We illustrate the algorithm performance using archival Hubble Space Telescope (HST) images, but the approach may also be considered for ground-based data acquired with Angular Differential Imaging (ADI) or integral-field spectrographs (IFS).
515 citations
TL;DR: The Weizmann Interactive Supernova data REPository (WISeREP) as discussed by the authors is an SQL-based database with an interactive web-based graphical interface for supernova data.
Abstract: We have entered an era of massive data sets in astronomy. In particular, the number of supernova (SN) discoveries and classifications has substantially increased over the years from few tens to thousands per year. It is no longer the case that observations of a few prototypical events encapsulate most spectroscopic information about SNe, motivating the development of modern tools to collect, archive, organize and distribute spectra in general, and SN spectra in particular. For this reason we have developed the Weizmann Interactive Supernova data REPository - WISeREP - an SQL-based database (DB) with an interactive web-based graphical interface. The system serves as an archive of high quality SN spectra, including both historical (legacy) data as well as data that is accumulated by ongoing modern programs. The archive provides information about objects, their spectra, and related meta-data. Utilizing interactive plots, we provide a graphical interface to visualize data, perform line identification of the major relevant species, determine object redshifts, classify SNe and measure expansion velocities. Guest users may view and download spectra or other data that have been placed in the public domain. Registered users may also view and download data that are proprietary to specific programs with which they are associated. The DB currently holds >8000 spectra, of which >5000 are public; the latter include published spectra from the Palomar Transient Factory, all of the SUSPECT archive, the Caltech-Core-Collapse Program, the CfA SN spectra archive and published spectra from the UC Berkeley SNDB repository. It offers an efficient and convenient way to archive data and share it with colleagues, and we expect that data stored in this way will be easy to access, increasing its visibility, usefulness and scientific impact.
477 citations
TL;DR: ExOFAST as discussed by the authors is a suite of routines written in IDL which is designed to fit exoplanetary transits and radial velocity variations simultaneously or separately, and characterize the parameter uncertainties and covariances with a Differential Evolution Markov Chain Monte Carlo method.
Abstract: We present EXOFAST, a fast, robust suite of routines written in IDL which is designed to fit exoplanetary transits and radial velocity variations simultaneously or separately, and characterize the parameter uncertainties and covariances with a Differential Evolution Markov Chain Monte Carlo method. We describe how our code incorporates both data sets to simultaneously derive stellar parameters along with the transit and RV parameters, resulting in more self-consistent results on an example fit of the discovery data of HAT-P-3b that is well-mixed in under five minutes on a standard desktop computer. We describe in detail how our code works and outline ways in which the code can be extended to include additional effects or generalized for the characterization of other data sets -- including non-planetary data sets. We discuss the pros and cons of several common ways to parameterize eccentricity, highlight a subtle mistake in the implementation of MCMC that could bias the inferred eccentricity of intrinsically circular orbits to significantly non-zero results, discuss a problem with IDL's built-in random number generator in its application to large MCMC fits, and derive a method to analytically fit the linear and quadratic limb darkening coefficients of a planetary transit. Finally, we explain how we achieved improved accuracy and over a factor of 100 improvement in the execution time of the transit model calculation. Our entire source code, along with an easy-to-use online interface for several basic features of our transit and radial velocity fitting, are available online at this http URL .
383 citations
TL;DR: XENON1T as mentioned in this paper is the next detector at the ton scale, with an experimental aim of probing WIMP interaction cross-sections above the spin-independent cross-section.
Abstract: The worldwide race towards direct dark matter detection in the form of Weakly Interacting Massive Particles (WIMPs) has been dramatically accelerated by the remarkable progress and evolution of liquid xenon time projection chambers (LXeTPCs). With a realistic discovery potential, XENON100 has already reached a sensitivity of $7\times10^{-45}\,
{cm}^2$, and continues to accrue data at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy towards its ultimate sensitivity reach at the $\sigma_{
{SI}}\sim 2\times10^{-45}\,
{cm}^2$ level for the spin-independent WIMP-nucleon cross-section. To fully explore the favoured parameter space for WIMP dark matter in search of a first robust and statistically significant discovery, or to confirm any hint of a signal from \Xehund, the next phase of the XENON program will be a detector at the ton scale - XENON1T. The XENON1T detector, based on 2.2 ton of LXe viewed by low radioactivity photomultiplier tubes and housed in a water Cherenkov muon veto at LNGS, is presented. With an experimental aim of probing WIMP interaction cross-sections above of order $\sigma_{
{SI}}\sim 2\times10^{-47}\,
{cm}^2$ within 2 years of operation, XENON1T will provide the sensitivity to probe a particularly favourable region of electroweak physics on a timescale compatible with complementary ground and satellite based indirect searches and with accelerator dark matter searches at the LHC. Indeed, for a $\sigma_{
{SI}} \sim 10^{-45}\,
{cm}^2$ and $100 \,
{GeV/c^2}$ WIMP mass, XENON1T could detect of order 100 events in this exposure, providing statistics for placing significant constraints on the WIMP mass.
TL;DR: The Murchison Widefield Array (MWA) as discussed by the authors is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency.
Abstract: Significant new opportunities for astrophysics and cosmology have been identified at low radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency. The telescope will enable new advances along four key science themes, including searching for redshifted 21 cm emission from the epoch of reionisation in the early Universe; Galactic and extragalactic all-sky southern hemisphere surveys; time-domain astrophysics; and solar, heliospheric, and ionospheric science and space weather. The Murchison Widefield Array is located in Western Australia at the site of the planned Square Kilometre Array (SKA) low-band telescope and is the only low-frequency SKA precursor facility. In this paper, we review the performance properties of the Murchison Widefield Array and describe its primary scientific objectives.
TL;DR: In this paper, the authors developed an absolute flux density scale for cm-wavelength astronomy by combining accurate flux density ratios determined by the VLA between the planet Mars and a set of potential calibrators with the thermophysical emission model of Mars, adjusted to the absolute scale established by WMAP.
Abstract: We develop an absolute flux density scale for cm-wavelength astronomy by combining accurate flux density ratios determined by the VLA between the planet Mars and a set of potential calibrators with the Rudy thermophysical emission model of Mars, adjusted to the absolute scale established by WMAP. The radio sources 3C123, 3C196, 3C286 and 3C295 are found to be varying at a level of less than ~5% per century at all frequencies between 1 and 50 GHz, and hence are suitable as flux density standards. We present polynomial expressions for their spectral flux densities, valid from 1 to 50 GHz, with absolute accuracy estimated at 1-3% depending on frequency. Of the four sources, 3C286 is the most compact and has the flattest spectral index, making it the most suitable object on which to establish the spectral flux density scale. The sources 3C48, 3C138, 3C147, NGC7027, NGC6542, and MWC349 show significant variability on various timescales. Polynomial coefficients for the spectral flux density are developed for 3C48, 3C138, and 3C147 for each of the seventeen observation dates, spanning 1983 through 2012. The planets Venus, Uranus, and Neptune are included in our observations, and we derive their brightness temperatures over the same frequency range.
TL;DR: In this paper, a study of the correlation between extinction in the Milky Way and the equivalent width of the NaI D absorption doublet is presented, and the authors provide a brief guide on how to best estimate extinction to extragalactic sources such as supernovae, under a variety of circumstances.
Abstract: Dust extinction and reddening are ubiquitous in astronomical observations and are often a major source of systematic uncertainty. We present here a study of the correlation between extinction in the Milky Way and the equivalent width of the NaI D absorption doublet. Our sample includes more than 100 high resolution spectra from the KECK telescopes and nearly a million low resolution spectra from the Sloan Digital Sky Survey (SDSS). We measure the correlation to unprecedented precision, constrain its shape, and derive an empirical relation between these quantities with a dispersion of order 0.15 magnitude in E(B-V). From the shape of the curve of growth we further show that a typical sight line through the Galaxy, as seen within the SDSS footprint, crosses about three dust clouds. We provide a brief guide on how to best estimate extinction to extragalactic sources such as supernovae, using the NaI D absorption feature, under a variety of circumstances.
TL;DR: Gaia as mentioned in this paper is the next mission of the European Space Agency (ESA), following up on the success of the Hipparcos mission, and is the first mission to obtain absolute astrometry, broad-band photometry, and low-resolution spectro-photometry.
Abstract: Gaia is the next astrometry mission of the European Space Agency (ESA), following up on the success of the Hipparcos mission. With a focal plane containing 106 CCD detectors, Gaia will survey the entire sky and repeatedly observe the brightest 1,000 million objects, down to 20th magnitude, during its 5-year lifetime. Gaia's science data comprises absolute astrometry, broad-band photometry, and low-resolution spectro-photometry. Spectroscopic data with a resolving power of 11,500 will be obtained for the brightest 150 million sources, down to 17th magnitude. The thermo-mechanical stability of the spacecraft, combined with the selection of the L2 Lissajous point of the Sun-Earth/Moon system for operations, allows stellar parallaxes to be measured with standard errors less than 10 micro-arcsecond (muas) for stars brighter than 12th magnitude, 25 muas for stars at 15th magnitude, and 300 muas at magnitude 20. Photometric standard errors are in the milli-magnitude regime. The spectroscopic data allows the measurement of radial velocities with errors of 15 km/s at magnitude 17. Gaia's primary science goal is to unravel the kinematical, dynamical, and chemical structure and evolution of the Milky Way. In addition, Gaia's data will touch many other areas of science, e.g., stellar physics, solar-system bodies, fundamental physics, and exo-planets. The Gaia spacecraft is currently in the qualification and production phase. With a launch in 2013, the final catalogue is expected in 2021. The science community in Europe, organised in the Data Processing and Analysis Consortium (DPAC), is responsible for the processing of the data.
TL;DR: In this paper, the authors presented parameterized broadband spectral models valid at frequencies between 30-300 MHz for six bright radio sources selected from the 3C survey, spread in Right Ascension from 0-24 hours.
Abstract: We present parameterized broadband spectral models valid at frequencies between 30-300 MHz for six bright radio sources selected from the 3C survey, spread in Right Ascension from 0-24 hours. For each source, data from the literature are compiled and tied to a common flux density scale. These data are then used to parameterize an analytic polynomial spectral calibration model. The optimal polynomial order in each case is determined using the ratio of the Bayesian evidence for the candidate models. Maximum likelihood parameter values for each model are presented, with associated errors, and the percentage error in each model as a function of frequency is derived. These spectral models are intended as an initial reference for science from the new generation of low frequency telescopes now coming on line, with particular emphasis on the Low Frequency Array (LOFAR).
TL;DR: In this article, the authors apply the delay transformation to each baseline of an interferometer to concentrate smooth-spectrum foregrounds within the bounds of the maximum geometric delays physically realizable on that baseline.
Abstract: A critical challenge in measuring the power spectrum of 21cm emission from cosmic reionization is compensating for the frequency dependence of an interferometer's sampling pattern, which can cause smooth-spectrum foregrounds to appear unsmooth and degrade the separation between foregrounds and the target signal. In this paper, we present an approach to foreground removal that explicitly accounts for this frequency dependence. We apply the delay transformation introduced in Parsons & Backer (2009) to each baseline of an interferometer to concentrate smooth-spectrum foregrounds within the bounds of the maximum geometric delays physically realizable on that baseline. By focusing on delay-modes that correspond to image-domain regions beyond the horizon, we show that it is possible to avoid the bulk of smooth-spectrum foregrounds. We map the point-spread function of delay-modes to k-space, showing that delay-modes that are uncorrupted by foregrounds also represent samples of the three-dimensional power spectrum, and can be used to constrain cosmic reionization. Because it uses only spectral smoothness to differentiate foregrounds from the targeted 21cm signature, this per-baseline analysis approach relies on spectrally- and spatially-smooth instrumental responses for foreground removal. For sufficient levels of instrumental smoothness relative to interfering foregrounds, this technique substantially reduces the level of calibration previously thought necessary to detect 21cm reionization. As a result, this approach places fewer constraints on antenna configuration within an array, facilitating the adoption of configurations optimized for power-spectrum sensitivity. Under these assumptions, we demonstrate the potential for PAPER to detect 21cm reionization at an amplitude of 10 mK^2 near k~0.2h Mpc^-1 with 132 dipoles in 7 months of observing.
TL;DR: A new algorithm is described, that is shown to perform quite well, has linear time complexity and is fast enough to be applied in modern high resolution observatories.
Abstract: A technique is described that is used to improve the detection of radio-frequency interference in astronomical radio observatories. It is applied on a two-dimensional interference mask after regular detection in the time-frequency domain with existing techniques. The scale-invariant rank (SIR) operator is defined, which is a one-dimensional mathematical morphology technique that can be used to find adjacent intervals in the time or frequency domain that are likely to be affected by RFI. The technique might also be applicable in other areas in which morphological scale-invariant behaviour is desired, such as source detection. A new algorithm is described, that is shown to perform quite well, has linear time complexity and is fast enough to be applied in modern high resolution observatories. It is used in the default pipeline of the LOFAR observatory.
TL;DR: A core-collapse supernova will produce an enormous burst of neutrinos of all flavors in the few-tens-of-MeV range as mentioned in this paper, which will yield answers to many physics and astrophysics questions.
Abstract: A core-collapse supernova will produce an enormous burst of neutrinos of all flavors in the few-tens-of-MeV range. Measurement of the flavor, time and energy structure of a nearby core-collapse neutrino burst will yield answers to many physics and astrophysics questions. The neutrinos left over from past cosmic supernovae are also observable, and their detection will improve knowledge of core collapse rates and average neutrino emission. This review describes experimental techniques for detection of core-collapse neutrinos, as well as the sensitivities of current and future detectors.
Posted Content•
TL;DR: In 2011, the Science Definition Team (SDT) was chartered to work with NASA's Wide Field Infra-Red Survey Telescope (WFIRST) Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST as discussed by the authors.
Abstract: In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.
TL;DR: The 4MOST project as mentioned in this paper proposes a wide-field (3 square degree, goal >5 square degree), high-multiplex (1500 fibres, goal 3000 fibres) spectroscopic survey facility for an ESO 4m-class telescope.
Abstract: The 4MOST consortium is currently halfway through a Conceptual Design study for ESO with the aim to develop a wide-field (>3 square degree, goal >5 square degree), high-multiplex (>1500 fibres, goal 3000 fibres) spectroscopic survey facility for an ESO 4m-class telescope (VISTA). 4MOST will run permanently on the telescope to perform a 5 year public survey yielding more than 20 million spectra at resolution R~5000 ({\lambda}=390-1000 nm) and more than 2 million spectra at R~20,000 (395-456.5 nm & 587-673 nm). The 4MOST design is especially intended to complement three key all-sky, space-based observatories of prime European interest: Gaia, eROSITA and Euclid. Initial design and performance estimates for the wide-field corrector concepts are presented. We consider two fibre positioner concepts, a well-known Phi-Theta system and a new R-Theta concept with a large patrol area. The spectrographs are fixed configuration two-arm spectrographs, with dedicated spectrographs for the high- and low-resolution. A full facility simulator is being developed to guide trade-off decisions regarding the optimal field-of-view, number of fibres needed, and the relative fraction of high-to-low resolution fibres. Mock catalogues with template spectra from seven Design Reference Surveys are simulated to verify the science requirements of 4MOST. The 4MOST consortium aims to deliver the full 4MOST facility by the end of 2018 and start delivering high-level data products for both consortium and ESO community targets a year later with yearly increments.
TL;DR: In this article, an atmospheric model was developed as part of the Austrian ESO In-Kind contribution to increase the accuracy of the exposure time calculator of the European Southern Observatory's (ESO) Very Large Telescope (VLT) at Cerro Paranal.
Abstract: The Earth's atmosphere affects ground-based astronomical observations. Scattering, absorption, and radiation processes deteriorate the signal-to-noise ratio of the data received. For scheduling astronomical observations it is, therefore, important to accurately estimate the wavelength-dependent effect of the Earth's atmosphere on the observed flux. In order to increase the accuracy of the exposure time calculator of the European Southern Observatory's (ESO) Very Large Telescope (VLT) at Cerro Paranal, an atmospheric model was developed as part of the Austrian ESO In-Kind contribution. It includes all relevant components, such as scattered moonlight, scattered starlight, zodiacal light, atmospheric thermal radiation and absorption, and non-thermal airglow emission. This paper focuses on atmospheric scattering processes that mostly affect the blue ( 0.55 mum) wavelength regime. While the former is mainly investigated by means of radiative transfer models, the intensity and variability of the latter is studied with a sample of 1186 VLT FORS1 spectra. For a set of parameters such as the object altitude angle, Moon-object angular distance, ecliptic latitude, bimonthly period, and solar radio flux, our model predicts atmospheric radiation and transmission at a requested resolution. A comparison of our model with the FORS1 spectra and photometric data for the night-sky brightness from the literature, suggest a model accuracy of about 20%. This is a significant improvement with respect to existing predictive atmospheric models for astronomical exposure time calculators.
TL;DR: In this paper, a new and significantly improved version of the presearch data conditioning (PDC) module of the Kepler data analysis pipeline has been proposed, which utilizes a Bayesian approach for removing systematics, reliably corrects errors in the light curves while at the same time preserving planet transits and other astrophysically interesting signals.
Abstract: Kepler provides light curves of 156,000 stars with unprecedented precision. However, the raw data as they come from the spacecraft contain significant systematic and stochastic errors. These errors, which include discontinuities, systematic trends, and outliers, obscure the astrophysical signals in the light curves. To correct these errors is the task of the Presearch Data Conditioning (PDC) module of the Kepler data analysis pipeline. The original version of PDC in Kepler did not meet the extremely high performance requirements for the detection of miniscule planet transits or highly accurate analysis of stellar activity and rotation. One particular deficiency was that astrophysical features were often removed as a side-effect to removal of errors. In this paper we introduce the completely new and significantly improved version of PDC which was implemented in Kepler SOC 8.0. This new PDC version, which utilizes a Bayesian approach for removal of systematics, reliably corrects errors in the light curves while at the same time preserving planet transits and other astrophysically interesting signals. We describe the architecture and the algorithms of this new PDC module, show typical errors encountered in Kepler data, and illustrate the corrections using real light curve examples.
TL;DR: In this article, a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes is presented, which rigorously and self-consistently treats measurement uncertainties, including calibration and noise, resulting in more precise SED fits.
Abstract: We present a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes. Under the standard assumption of optically thin single temperature (T) sources the dust SED as represented by a power--law modified black body is subject to a strong degeneracy between T and the spectral index beta. The traditional non-hierarchical approaches, typically based on chi-square minimization, are severely limited by this degeneracy, as it produces an artificial anti-correlation between T and beta even with modest levels of observational noise. The hierarchical Bayesian method rigorously and self-consistently treats measurement uncertainties, including calibration and noise, resulting in more precise SED fits. As a result, the Bayesian fits do not produce any spurious anti-correlations between the SED parameters due to measurement uncertainty. We demonstrate that the Bayesian method is substantially more accurate than the chi-square fit in recovering the SED parameters, as well as the correlations between them. As an illustration, we apply our method to Herschel and sub millimeter ground-based observations of the star-forming Bok globule CB244. This source is a small, nearby molecular cloud containing a single low-mass protostar and a starless core. We find that T and beta are weakly positively correlated -- in contradiction with the chi-square fits, which indicate a T-beta anti-correlation from the same data-set. Additionally, in comparison to the chi-square fits the Bayesian SED parameter estimates exhibit a reduced range in values.
TL;DR: In this paper, the authors identify four fundamental power spectrum shapes produced by mode-mixing that will affect all upcoming observations and forecast the shape of an EoR window that is mostly free of contamination.
Abstract: Contamination from instrumental effects interacting with bright astrophysical sources is the primary impediment to measuring Epoch of Reionization and BAO 21 cm power spectra---an effect called mode-mixing. In this paper we identify four fundamental power spectrum shapes produced by mode-mixing that will affect all upcoming observations. We are able, for the first time, to explain the wedge-like structure seen in advanced simulations and to forecast the shape of an 'EoR window' that is mostly free of contamination. Understanding the origins of these contaminations also enables us to identify calibration and foreground subtraction errors below the imaging limit, providing a powerful new tool for precision observations.
TL;DR: In this article, the A-projection algorithm is applied to the Square Kilometer Array (LOFAR) to deal with non-unitary station beams and non-diagonal Mueller matrices.
Abstract: The aimed high sensitivities and large fields of view of the new generation of interferometers impose to reach high dynamic range of order $\sim$1:$10^6$ to 1:$10^8$ in the case of the Square Kilometer Array. The main problem is the calibration and correction of the Direction Dependent Effects (DDE) that can affect the electro-magnetic field (antenna beams, ionosphere, Faraday rotation, etc.). As shown earlier the A-Projection is a fast and accurate algorithm that can potentially correct for any given DDE in the imaging step. With its very wide field of view, low operating frequency ($\sim30-250$ MHz), long baselines, and complex station-dependent beam patterns, the Low Frequency Array (LOFAR) is certainly the most complex SKA precursor. In this paper we present a few implementations of A-Projection applied to LOFAR that can deal with non-unitary station beams and non-diagonal Mueller matrices. The algorithm is designed to correct for all the DDE, including individual antenna, projection of the dipoles on the sky, beam forming and ionospheric effects. We describe a few important algorithmic optimizations related to LOFAR's architecture allowing us to build a fast imager. Based on simulated datasets we show that A-Projection can give dramatic dynamic range improvement for both phased array beams and ionospheric effects. We will use this algorithm for the construction of the deepest extragalactic surveys, comprising hundreds of days of integration.
TL;DR: The DARWIN project as discussed by the authors is a design study for a next-generation, multi-ton dark matter detector in Europe, where liquid argon and liquid xenon are the target media for direct detection of dark matter candidates in the form of weakly interacting massive particles (WIMPs).
Abstract: DARWIN (dark matter wimp search with noble liquids) is a design study for a next-generation, multi-ton dark matter detector in Europe. Liquid argon and/or liquid xenon are the target media for the direct detection of dark matter candidates in the form of weakly interacting massive particles (WIMPs). Light and charge signals created by particle interactions in the active detector volume are observed via the time projection chamber technique. DARWIN is to probe the spin-independent, WIMP-nucleon cross section down 1e-48 cm2 and to measure WIMP-induced nuclear recoil spectra with high-statistics, should they be discovered by an existing or near-future experiment. After a brief introduction, I will describe the project, selected R&D topics, expected backgrounds and the physics reach.
TL;DR: In this article, the Parkes Pulsar Timing Array (PPTA) project has developed the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales.
Abstract: Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale.
TL;DR: In this article, the authors present the main techniques and R&D projects that will allow to build so-called ultimate WIMP detectors, capable of probing spin-independent interactions down to the unimaginably low cross section of 1e-48 cm2, before the irreducible neutrino background takes over.
Abstract: Direct dark matter searches are promising techniques to identify the nature of dark matter particles. I describe the future of this field of research, focussing on the question of what can be achieved in the next decade. I will present the main techniques and R&D projects that will allow to build so-called ultimate WIMP detectors, capable of probing spin-independent interactions down to the unimaginably low cross section of 1e-48 cm2, before the irreducible neutrino background takes over. If a discovery is within the reach of a near-future dark matter experiment, these detectors will be able to constrain WIMP properties such as its mass, scattering cross section and possibly spin. With input from the LHC and from indirect searches, direct detection experiments will hopefully allow to determine the local density and to constrain the local phase-space structure of our dark matter halo.
TL;DR: In this paper, the authors present the design and characterization of the POLARBEAR experiment, which utilizes a unique focal plane of 1,274 antenna-coupled, polarization sensitive TES bolometers cooled to 250 milliKelvin.
Abstract: We present the design and characterization of the POLARBEAR experiment. POLARBEAR will measure the polarization of the cosmic microwave background (CMB) on angular scales ranging from the experiment's 3.5 arcminute beam size to several degrees. The experiment utilizes a unique focal plane of 1,274 antenna-coupled, polarization sensitive TES bolometers cooled to 250 milliKelvin. Employing this focal plane along with stringent control over systematic errors, POLARBEAR has the sensitivity to detect the expected small scale B-mode signal due to gravitational lensing and search for the large scale B-mode signal from inflationary gravitational waves.
POLARBEAR was assembled for an engineering run in the Inyo Mountains of California in 2010 and was deployed in late 2011 to the Atacama Desert in Chile. An overview of the instrument is presented along with characterization results from observations in Chile.