Showing papers in "arXiv: Instrumentation and Methods for Astrophysics in 2014"

The K2 Mission: Characterization and Early results

Abstract: The K2 mission will make use of the Kepler spacecraft and its assets to expand upon Kepler's groundbreaking discoveries in the fields of exoplanets and astrophysics through new and exciting observations. K2 will use an innovative way of operating the spacecraft to observe target fields along the ecliptic for the next 2-3 years. Early science commissioning observations have shown an estimated photometric precision near 400 ppm in a single 30 minute observation, and a 6-hour photometric precision of 80 ppm (both at V=12). The K2 mission offers long-term, simultaneous optical observation of thousands of objects at a precision far better than is achievable from ground-based telescopes. Ecliptic fields will be observed for approximately 75-days enabling a unique exoplanet survey which fills the gaps in duration and sensitivity between the Kepler and TESS missions, and offers pre-launch exoplanet target identification for JWST transit spectroscopy. Astrophysics observations with K2 will include studies of young open clusters, bright stars, galaxies, supernovae, and asteroseismology.

A Technique for Extracting Highly Precise Photometry for the Two-Wheeled Kepler Mission

Abstract: The original Kepler mission achieved high photometric precision thanks to ultra-stable pointing enabled by use of four reaction wheels. The loss of two of these reaction wheels reduced the telescope's ability to point precisely for extended periods of time, and as a result, the photometric precision has suffered. We present a technique for generating photometric light curves from pixel-level data obtained with the two-wheeled extended Kepler mission, K2. Our photometric technique accounts for the non-uniform pixel response function of the Kepler detectors by correlating flux measurements with the spacecraft's pointing and removing the dependence. When we apply our technique to the ensemble of stars observed during the Kepler Two-Wheel Concept Engineering Test, we find improvements over raw K2 photometry by factors of 2-5, with noise properties qualitatively similar to Kepler targets at the same magnitudes. We find evidence that the improvement in photometric precision depends on each target's position in the Kepler field of view, with worst precision near the edges of the field. Overall, this technique restores the median attainable photometric precision to within a factor of two of the original Kepler photometric precision for targets ranging from 10$^{th}$ to 15$^{th}$ magnitude in the Kepler bandpass, peaking with a median precision within 35% that of Kepler for stars between 12$^{th}$ and 13$^{th}$ magnitude in the Kepler bandpass.

SPT-3G: A Next-Generation Cosmic Microwave Background Polarization Experiment on the South Pole Telescope

Abstract: We describe the design of a new polarization sensitive receiver, SPT-3G, for the 10-meter South Pole Telescope (SPT). The SPT-3G receiver will deliver a factor of ~20 improvement in mapping speed over the current receiver, SPTpol. The sensitivity of the SPT-3G receiver will enable the advance from statistical detection of B-mode polarization anisotropy power to high signal-to-noise measurements of the individual modes, i.e., maps. This will lead to precise (~0.06 eV) constraints on the sum of neutrino masses with the potential to directly address the neutrino mass hierarchy. It will allow a separation of the lensing and inflationary B-mode power spectra, improving constraints on the amplitude and shape of the primordial signal, either through SPT-3G data alone or in combination with BICEP-2/KECK, which is observing the same area of sky. The measurement of small-scale temperature anisotropy will provide new constraints on the epoch of reionization. Additional science from the SPT-3G survey will be significantly enhanced by the synergy with the ongoing optical Dark Energy Survey (DES), including: a 1% constraint on the bias of optical tracers of large-scale structure, a measurement of the differential Doppler signal from pairs of galaxy clusters that will test General Relativity on ~200 Mpc scales, and improved cosmological constraints from the abundance of clusters of galaxies.

The Tycho-Gaia astrometric solution. How to get 2.5 million parallaxes with less than one year of Gaia data

Abstract: Context. The first release of astrometric data from Gaia will contain the mean stellar positions and magnitudes from the first year of observations, and proper motions from the combination of Gaia data with Hipparcos prior information (HTPM). Aims. We study the potential of using the positions from the Tycho-2 Catalogue as additional information for a joint solution with early Gaia data. We call this the Tycho-Gaia astrometric solution (TGAS). Methods. We adapt Gaia's Astrometric Global Iterative Solution (AGIS) to incorporate Tycho information, and use simulated Gaia observations to demonstrate the feasibility of TGAS and to estimate its performance. Results. Using six to twelve months of Gaia data, TGAS could deliver positions, parallaxes and annual proper motions for the 2.5 million Tycho-2 stars, with sub-milliarcsecond accuracy. TGAS overcomes some of the limitations of the HTPM project and allows its execution half a year earlier. Furthermore, if the parallaxes from Hipparcos are not incorporated in the solution, they can be used as a consistency check of the TGAS/HTPM solution.

Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder

Abstract: A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.

The MaNGA Integral Field Unit Fiber Feed System for the Sloan 2.5 m Telescope

Abstract: We describe the design, manufacture, and performance of bare-fiber integral field units (IFUs) for the SDSS-IV survey MaNGA (Mapping Nearby Galaxies at APO) on the the Sloan 2.5 m telescope at Apache Point Observatory (APO). MaNGA is a luminosity-selected integral-field spectroscopic survey of 10,000 local galaxies covering 360-1030 nm at R ~ 2200. The IFUs have hexagonal dense packing of fibers with packing regularity of 3 um (RMS), and throughput of 96+/-0.5% from 350 nm to 1 um in the lab. Their sizes range from 19 to 127 fibers (3-7 hexagonal layers) using Polymicro FBP 120:132:150 um core:clad:buffer fibers to reach a fill fraction of 56%. High throughput (and low focal-ratio degradation) is achieved by maintaining the fiber cladding and buffer intact, ensuring excellent surface polish, and applying a multi-layer AR coating of the input and output surfaces. In operations on-sky, the IFUs show only an additional 2.3% FRD-related variability in throughput despite repeated mechanical stressing during plate plugging (however other losses are present). The IFUs achieve on-sky throughput 5% above the single-fiber feeds used in SDSS-III/BOSS, attributable to equivalent performance compared to single fibers and additional gains from the AR coating. The manufacturing process is geared toward mass-production of high-multiplex systems. The low-stress process involves a precision ferrule with hexagonal inner shape designed to lead inserted fibers to settle in a dense hexagonal pattern. The ferrule inner diameter is tapered at progressively shallower angles toward its tip and the final 2 mm are straight and only a few um larger than necessary to hold the desired number of fibers. This process scales to accommodate other fiber sizes and to IFUs with substantially larger fiber count. (Abridged)

BRITE-Constellation: nanosatellites for precision photometry of bright stars

Abstract: BRITE-Constellation (where BRITE stands for BRIght Target Explorer) is an international nanosatellite mission to monitor photometrically, in two colours, the brightness and temperature variations of stars generally brighter than mag(V) ~ 4, with precision and time coverage not possible from the ground. The current mission design consists of six nanosats (hence Constellation): two from Austria, two from Canada, and two from Poland. Each 7 kg nanosat carries an optical telescope of aperture 3 cm feeding an uncooled CCD. One instrument in each pair is equipped with a blue filter, the other with a red filter. Each BRITE instrument has a wide field of view (~24 degrees), so up to about 15 bright stars can be observed simultaneously, sampled in 32 pixel x 32 pixel sub-rasters. Photometry of additional fainter targets, with reduced precision but thorough time sampling, will be possible through onboard data processing. The BRITE sample is dominated by the most intrinsically luminous stars: massive stars seen at all evolutionary stages, and evolved medium-class stars at the very end of their nuclear burning phases. The goals of BRITE-Constellation are to (1) measure p- and g-mode pulsations to probe the interiors and ages of stars through asteroseismology; (2) look for varying spots on the stars surfaces carried across the stellar disks by rotation, which are the sources of co-rotating interaction regions in the winds of the most luminous stars, probably arising from magnetic subsurface convection; and (3) search for planetary transits.

Observations of Transiting Exoplanets with the James Webb Space Telescope (JWST), Publications of the Astronomical Society of the Pacific (PASP), December 2014

Abstract: This article summarizes a workshop held on March, 2014, on the potential of the James Webb Space Telescope (JWST) to revolutionize our knowledge of the physical properties of exoplanets through transit observations. JWST's unique combination of high sensitivity and broad wavelength coverage will enable the accurate measurement of transits with high signal-to-noise. Most importantly, JWST spectroscopy will investigate planetary atmospheres to determine atomic and molecular compositions, to probe vertical and horizontal structure, and to follow dynamical evolution, i.e. exoplanet weather. JWST will sample a diverse population of planets of varying masses and densities in a wide variety of environments characterized by a range of host star masses and metallicities, orbital semi-major axes and eccentricities. A broad program of exoplanet science could use a substantial fraction of the overall JWST mission.

Barycentric Corrections at 1 cm/s for precise Doppler velocities

Abstract: The goal of this paper is to establish the requirements of a barycentric correction with an RMS of $\lesssim 1$ cm/s, which is an order of magnitude better than necessary for the Doppler detection of true Earth analogs ($\sim9$ cm/s) We describe the theory and implementation of accounting for the effects on precise Doppler measurements of motion of the telescope through space, primarily from rotational and orbital motion of the Earth, and the motion of the solar system with respect to target star (ie the "barycentric correction") We describe the minimal algorithm necessary to accomplish this and how it differs from a na\"ive subtraction of velocities (ie a Galilean transformation) We demonstrate the validity of code we have developed from the California Planet Survey code via comparison with the pulsar timing package, TEMPO2 We estimate the magnitude of various terms and effects, including relativistic effects, and the errors associated with incomplete knowledge of telescope position, timing, and stellar position and motion We note that chromatic aberration will create uncertainties in the time of observation, which will complicate efforts to detect true Earth analogs Our code is available for public use and validation

25 Years of Self-Organized Criticality: Solar and Astrophysics

Abstract: Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC concept was initially applied to solar flares, stellar flares, and magnetospheric substorms, and later extended to the radiation belt, the heliosphere, lunar craters, the asteroid belt, the Saturn ring, pulsar glitches, soft X-ray repeaters, blazars, black-hole objects, cosmic rays, and boson clouds. The application of SOC concepts has been performed by numerical cellular automaton simulations, by analytical calculations of statistical (powerlaw-like) distributions based on physical scaling laws, and by observational tests of theoretically predicted size distributions and waiting time distributions. Attempts have been undertaken to import physical models into the numerical SOC toy models, such as the discretization of magneto-hydrodynamics (MHD) processes. The novel applications stimulated also vigorous debates about the discrimination between SOC models, SOC-like, and non-SOC processes, such as phase transitions, turbulence, random-walk diffusion, percolation, branching processes, network theory, chaos theory, fractality, multi-scale, and other complexity phenomena. We review SOC studies from the last 25 years and highlight new trends, open questions, and future challenges, as discussed during two recent ISSI workshops on this theme.

SKIRT: an Advanced Dust Radiative Transfer Code with a User-Friendly Architecture

Abstract: We discuss the architecture and design principles that underpin the latest version of SKIRT, a state-of-the-art open source code for simulating continuum radiation transfer in dusty astrophysical systems, such as spiral galaxies and accretion disks. SKIRT employs the Monte Carlo technique to emulate the relevant physical processes including scattering, absorption and emission by the dust. The code features a wealth of built-in geometries, radiation source spectra, dust characterizations, dust grids, and detectors, in addition to various mechanisms for importing snapshots generated by hydrodynamical simulations. The configuration for a particular simulation is defined at run-time through a user-friendly interface suitable for both occasional and power users. These capabilities are enabled by careful C++ code design. The programming interfaces between components are well defined and narrow. Adding a new feature is usually as simple as adding another class; the user interface automatically adjusts to allow configuring the new options. We argue that many scientific codes, like SKIRT, can benefit from careful object-oriented design and from a friendly user interface, even if it is not a graphical user interface.

Introduction to astroML: Machine Learning for Astrophysics

Abstract: Astronomy and astrophysics are witnessing dramatic increases in data volume as detectors, telescopes and computers become ever more powerful. During the last decade, sky surveys across the electromagnetic spectrum have collected hundreds of terabytes of astronomical data for hundreds of millions of sources. Over the next decade, the data volume will enter the petabyte domain, and provide accurate measurements for billions of sources. Astronomy and physics students are not traditionally trained to handle such voluminous and complex data sets. In this paper we describe astroML; an initiative, based on Python and scikit-learn, to develop a compendium of machine learning tools designed to address the statistical needs of the next generation of students and astronomical surveys. We introduce astroML and present a number of example applications that are enabled by this package.

PyNeb: a new tool for analyzing emission lines. I. Code description and validation of results

Abstract: Analysis of emission lines in gaseous nebulae yields direct measures of physical conditions and chemical abundances and is the cornerstone of nebular astrophysics. Although the physical problem is conceptually simple, its practical complexity can be overwhelming since the amount of data to be analyzed steadily increases; furthermore, results depend crucially on the input atomic data, whose determination also improves each year. To address these challenges we created PyNeb, an innovative code for analyzing emission lines. PyNeb computes physical conditions and ionic and elemental abundances, and produces both theoretical and observational diagnostic plots. It is designed to be portable, modular, and largely customizable in aspects such as the atomic data used, the format of the observational data to be analyzed, and the graphical output. It gives full access to the intermediate quantities of the calculation, making it possible to write scripts tailored to the specific type of analysis one wants to carry out. In the case of collisionally excited lines, PyNeb works by solving the equilibrium equations for an n-level atom; in the case of recombination lines, it works by interpolation in emissivity tables. The code offers a choice of extinction laws and ionization correction factors, which can be complemented by user-provided recipes. It is entirely written in the python programming language and uses standard python libraries. It is fully vectorized, making it apt for analyzing huge amounts of data. The code is stable and has been benchmarked against IRAF/NEBULAR. It is public, fully documented, and has already been satisfactorily used in a number of published papers.

Probing for Exoplanets Hiding in Dusty Debris Disks: Disk Imaging, Characterization, and Exploration with HST/STIS Multi-Roll Coronagraphy

Abstract: Spatially resolved scattered-light images of circumstellar (CS) debris in exoplanetary systems constrain the physical properties and orbits of the dust particles in these systems. They also inform on co-orbiting (but unseen) planets, systemic architectures, and forces perturbing starlight-scattering CS material. Using HST/STIS optical coronagraphy, we have completed the observational phase of a program to study the spatial distribution of dust in ten CS debris systems, and one "mature" protoplanetrary disk all with HST pedigree, using PSF-subtracted multi-roll coronagraphy. These observations probe stellocentric distances > 5 AU for the nearest stars, and simultaneously resolve disk substructures well beyond, corresponding to the giant planet and Kuiper belt regions in our Solar System. They also disclose diffuse very low-surface brightness dust at larger stellocentric distances. We present new results inclusive of fainter disks such as HD92945 confirming, and better revealing, the existence of a narrow inner debris ring within a larger diffuse dust disk. Other disks with ring-like sub-structures, significant asymmetries and complex morphologies include: HD181327 with a posited spray of ejecta from a recent massive collision in an exo-Kuiper belt; HD61005 suggested interacting with the local ISM; HD15115 & HD32297, discussed also in the context of environmental interactions. These disks, and HD15745, suggest debris system evolution cannot be treated in isolation. For AU Mic's edge-on disk, out-of-plane surface brightness asymmetries at > 5 AU may implicate one or more planetary perturbers. Time resolved images of the MP Mus proto-planetary disk provide spatially resolved temporal variability in the disk illumination. These and other new images from our program enable direct inter-comparison of the architectures of these exoplanetary debris systems in the context of our own Solar System.

WISE photometry for 400 million SDSS sources

Abstract: We present photometry of images from the Wide-Field Infrared Survey Explorer (WISE; Wright et al. 2010) of over 400 million sources detected by the Sloan Digital Sky Survey (SDSS; York et al. 2000). We use a "forced photometry" technique, using measured SDSS source positions, star-galaxy separation and galaxy profiles to define the sources whose fluxes are to be measured in the WISE images. We perform photometry with The Tractor image modeling code, working on our "unWISE" coaddds and taking account of the WISE point-spread function and a noise model. The result is a measurement of the flux of each SDSS source in each WISE band. Many sources have little flux in the WISE bands, so often the measurements we report are consistent with zero. However, for many sources we get three- or four-sigma measurements; these sources would not be reported by the WISE pipeline and will not appear in the WISE catalog, yet they can be highly informative for some scientific questions. In addition, these small-signal measurements can be used in stacking analyses at catalog level. The forced photometry approach has the advantage that we measure a consistent set of sources between SDSS and WISE, taking advantage of the resolution and depth of the SDSS images to interpret the WISE images; objects that are resolved in SDSS but blended together in WISE still have accurate measurements in our photometry. Our results, and the code used to produce them, are publicly available at this http URL

Ultra Low Surface Brightness Imaging with the Dragonfly Telephoto Array

Abstract: We describe the Dragonfly Telephoto Array, a robotic imaging system optimized for the detection of extended ultra low surface brightness structures. The array consists of eight Canon 400mm f/2.8 telephoto lenses coupled to eight science-grade commercial CCD cameras. The lenses are mounted on a common framework and are co-aligned to image simultaneously the same position on the sky. The system provides an imaging capability equivalent to a 0.4m aperture f/1.0 refractor with a 2.6 deg X 1.9 deg field of view. The system has no obstructions in the light path, optimized baffling, and internal optical surfaces coated with a new generation of anti-reflection coatings based on sub-wavelength nanostructures. As a result, the array's point spread function has a factor of ~10 less scattered light at large radii than well-baffled reflecting telescopes. The Dragonfly Telephoto Array is capable of imaging extended structures to surface brightness levels below 30 mag/arcsec^2 in 10h integrations (without binning or foreground star removal). This is considerably deeper than the surface brightness limit of any existing wide-field telescope. At present no systematic errors limiting the usefulness of much longer integration times has been identified. With longer integrations (50-100h), foreground star removal and modest binning the Dragonfly Telephoto Array is capable of probing structures with surface brightnesses below 32 mag/arcsec^2. Detection of structures at these surface brightness levels may hold the key to solving the "missing substructure" and "missing satellite" problems of conventional hierarchical galaxy formation models. The Dragonfly Telephoto Array is therefore executing a fully-automated multi-year imaging survey of a complete sample of nearby galaxies in order to undertake the first census of ultra-faint substructures in the nearby Universe.

Estimating extragalactic Faraday rotation

Abstract: (abridged) Observations of Faraday rotation for extragalactic sources probe magnetic fields both inside and outside the Milky Way. Building on our earlier estimate of the Galactic contribution, we set out to estimate the extragalactic contributions. We discuss the problems involved; in particular, we point out that taking the difference between the observed values and the Galactic foreground reconstruction is not a good estimate for the extragalactic contributions. We point out a degeneracy between the contributions to the observed values due to extragalactic magnetic fields and observational noise and comment on the dangers of over-interpreting an estimate without taking into account its uncertainty information. To overcome these difficulties, we develop an extended reconstruction algorithm based on the assumption that the observational uncertainties are accurately described for a subset of the data, which can overcome the degeneracy with the extragalactic contributions. We present a probabilistic derivation of the algorithm and demonstrate its performance using a simulation, yielding a high quality reconstruction of the Galactic Faraday rotation foreground, a precise estimate of the typical extragalactic contribution, and a well-defined probabilistic description of the extragalactic contribution for each data point. We then apply this reconstruction technique to a catalog of Faraday rotation observations. We vary our assumptions about the data, showing that the dispersion of extragalactic contributions to observed Faraday depths is most likely lower than 7 rad/m^2, in agreement with earlier results, and that the extragalactic contribution to an individual data point is poorly constrained by the data in most cases.

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

Abstract: The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1 0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

ESPRESSO: The next European exoplanet hunter

Abstract: The acronym ESPRESSO stems for Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations; this instrument will be the next VLT high resolution spectrograph. The spectrograph will be installed at the Combined-Coud\'e Laboratory of the VLT and linked to the four 8.2 m Unit Telescopes (UT) through four optical Coud\'e trains. ESPRESSO will combine efficiency and extreme spectroscopic precision. ESPRESSO is foreseen to achieve a gain of two magnitudes with respect to its predecessor HARPS, and to improve the instrumental radial-velocity precision to reach the 10 cm/s level. It can be operated either with a single UT or with up to four UTs, enabling an additional gain in the latter mode. The incoherent combination of four telescopes and the extreme precision requirements called for many innovative design solutions while ensuring the technical heritage of the successful HARPS experience. ESPRESSO will allow to explore new frontiers in most domains of astrophysics that require precision and sensitivity. The main scientific drivers are the search and characterization of rocky exoplanets in the habitable zone of quiet, nearby G to M-dwarfs and the analysis of the variability of fundamental physical constants. The project passed the final design review in May 2013 and entered the manufacturing phase. ESPRESSO will be installed at the Paranal Observatory in 2016 and its operation is planned to start by the end of the same year.

CLASS: The Cosmology Large Angular Scale Surveyor

Abstract: The Cosmology Large Angular Scale Surveyor (CLASS) is an experiment to measure the signature of a gravita-tional-wave background from inflation in the polarization of the cosmic microwave background (CMB). CLASS is a multi-frequency array of four telescopes operating from a high-altitude site in the Atacama Desert in Chile. CLASS will survey 70\% of the sky in four frequency bands centered at 38, 93, 148, and 217 GHz, which are chosen to straddle the Galactic-foreground minimum while avoiding strong atmospheric emission lines. This broad frequency coverage ensures that CLASS can distinguish Galactic emission from the CMB. The sky fraction of the CLASS survey will allow the full shape of the primordial B-mode power spectrum to be characterized, including the signal from reionization at low $\ell$. Its unique combination of large sky coverage, control of systematic errors, and high sensitivity will allow CLASS to measure or place upper limits on the tensor-to-scalar ratio at a level of $r=0.01$ and make a cosmic-variance-limited measurement of the optical depth to the surface of last scattering, $\tau$.

Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets

Abstract: We present the use of continuous-time autoregressive moving average (CARMA) models as a method for estimating the variability features of a light curve, and in particular its power spectral density (PSD). CARMA models fully account for irregular sampling and measurement errors, making them valuable for quantifying variability, forecasting and interpolating light curves, and for variability-based classification. We show that the PSD of a CARMA model can be expressed as a sum of Lorentzian functions, which makes them extremely flexible and able to model a broad range of PSDs. We present the likelihood function for light curves sampled from CARMA processes, placing them on a statistically rigorous foundation, and we present a Bayesian method to infer the probability distribution of the PSD given the measured lightcurve. Because calculation of the likelihood function scales linearly with the number of data points, CARMA modeling scales to current and future massive time-domain data sets. We conclude by applying our CARMA modeling approach to light curves for an X-ray binary, two AGN, a long-period variable star, and an RR-Lyrae star, in order to illustrate their use, applicability, and interpretation.

DIAMONDS: a new Bayesian Nested Sampling tool. Application to Peak Bagging of solar-like oscillations

Abstract: To exploit the full potential of Kepler light curves, sophisticated and robust analysis tools are now required more than ever. Characterizing single stars with an unprecedented level of accuracy and subsequently analyzing stellar populations in detail are fundamental to further constrain stellar structure and evolutionary models. We developed a new code, termed Diamonds, for Bayesian parameter estimation and model comparison by means of the nested sampling Monte Carlo (NSMC) algorithm, an efficient and powerful method very suitable for high-dimensional and multi-modal problems. A detailed description of the features implemented in the code is given with a focus on the novelties and differences with respect to other existing methods based on NSMC. Diamonds is then tested on the bright F8 V star KIC~9139163, a challenging target for peak-bagging analysis due to its large number of oscillation peaks observed, which are coupled to the blending that occurs between $\ell=2,0$ peaks, and the strong stellar background signal. We further strain the performance of the approach by adopting a 1147.5 days-long Kepler light curve. The Diamonds code is able to provide robust results for the peak-bagging analysis of KIC~9139163. We test the detection of different astrophysical backgrounds in the star and provide a criterion based on the Bayesian evidence for assessing the peak significance of the detected oscillations in detail. We present results for 59 individual oscillation frequencies, amplitudes and linewidths and provide a detailed comparison to the existing values in the literature. Lastly, we successfully demonstrate an innovative approach to peak bagging that exploits the capability of Diamonds to sample multi-modal distributions, which is of great potential for possible future automatization of the analysis technique.

Direct detection of exoplanets in the 3 -- 10 micron range with E-ELT/METIS

Abstract: We quantify the scientific potential for exoplanet imaging with the Mid-infrared E-ELT Imager and Spectrograph (METIS) foreseen as one of the instruments of the European Extremely Large Telescope (E-ELT). We focus on two main science cases: (1) the direct detection of known gas giant planets found by radial velocity (RV) searches; and (2) the direct detection of small (1 - 4 R_earth) planets around the nearest stars. Under the assumptions made in our modeling, in particular on the achievable inner working angle and sensitivity, our analyses reveal that within a reasonable amount of observing time METIS is able to image >20 already known, RV-detected planets in at least one filter. Many more suitable planets with dynamically determined masses are expected to be found in the coming years with the continuation of RV-surveys and the results from the GAIA astrometry mission. In addition, by extrapolating the statistics for close-in planets found by \emph{Kepler}, we expect METIS might detect ~10 small planets with equilibrium temperatures between 200 - 500 K around the nearest stars. This means that (1) METIS will help constrain atmospheric models for gas giant planets by determining for a sizable sample their luminosity, temperature and orbital inclination; and (2) METIS might be the first instrument to image a nearby (super-)Earth-sized planet with an equilibrium temperature near that expected to enable liquid water on a planet surface.

UVMULTIFIT: A versatile tool for fitting astronomical radio interferometric data

Abstract: The analysis of astronomical interferometric data is often performed on the images obtained after deconvolution of the interferometer's point spread function (PSF). This strategy can be understood (especially for cases of sparse arrays) as fitting models to models, since the deconvolved images are already non-unique model representations of the actual data (i.e., the visibilities). Indeed, the interferometric images may be affected by visibility gridding, weighting schemes (e.g., natural vs. uniform), and the particulars of the (non-linear) deconvolution algorithms. Fitting models to the direct interferometric observables (i.e., the visibilities) is preferable in the cases of simple (analytical) sky intensity distributions. In this paper, we present UVMULTIFIT, a versatile library for fitting visibility data, implemented in a Python-based framework. Our software is currently based on the CASA package, but can be easily adapted to other analysis packages, provided they have a Python API. We have tested the software with synthetic data, as well as with real observations. In some cases (e.g., sources with sizes smaller than the diffraction limit of the interferometer), the results from the fit to the visibilities (e.g., spectra of close by sources) are far superior to the output obtained from the mere analysis of the deconvolved images. UVMULTIFIT is a powerful improvement of existing tasks to extract the maximum amount of information from visibility data, especially in cases close to the sensitivity/resolution limits of interferometric observations.

An improved source-subtracted and destriped 408 MHz all-sky map

Abstract: The all-sky 408 MHz map of Haslam et al. is one the most important total-power radio surveys. It has been widely used to study diffuse synchrotron radiation from our Galaxy and as a template to remove foregrounds in cosmic microwave background data. However, there are a number of issues associated with it that must be dealt with, including large-scale striations and contamination from extragalactic radio sources. We have re-evaluated and re-processed the rawest data available to produce a new and improved 408 MHz all-sky map. We first quantify the positional accuracy ($\approx 7$ arcmin) and effective beam ($56.0\pm1.0$ arcmin) of the four individual surveys from which it was assembled. Large-scale striations associated with $1/f$ noise in the scan direction are reduced to a level $\ll 1$ K using a Fourier-based filtering technique. The most important improvement results from the removal of extragalactic sources. We have used an iterative combination of two techniques -- two-dimensional Gaussian fitting and minimum curvature spline surface inpainting -- to remove the brightest sources ($\gtrsim 2$ Jy), which provides a significant improvement over previous versions of the map. We quantify the impact with power spectra and a template fitting analysis of foregrounds to the WMAP data. The new map is publicly available and is recommended as the template of choice for large-scale diffuse Galactic synchrotron emission. We also provide a higher resolution map with small-scale fluctuations added, assuming a power-law angular power spectrum down to the pixel scale (1.7 arcmin). This should prove useful in simulations used for studying the feasibility of detecting HI fluctuations from the Epoch of Reionization.

Laboratory formation of fullerenes from PAHs: Top-down interstellar chemistry

Abstract: Interstellar molecules are thought to build up in the shielded environment of molecular clouds or in the envelope of evolved stars. This follows many sequential reaction steps of atoms and simple molecules in the gas phase and/or on (icy) grain surfaces. However, these chemical routes are highly inefficient for larger species in the tenuous environment of space as many steps are involved and, indeed, models fail to explain the observed high abundances. This is definitely the case for the C$_{60}$ fullerene, recently identified as one of the most complex molecules in the interstellar medium. Observations have shown that, in some PDRs, its abundance increases close to strong UV-sources. In this letter we report laboratory findings in which C$_{60}$ formation can be explained by characterizing the photochemical evolution of large PAHs. Sequential H losses lead to fully dehydrogenated PAHs and subsequent losses of C$_{2}$ units convert graphene into cages. Our results present for the first time experimental evidence that PAHs in excess of 60 C-atoms efficiently photo-isomerize to Buckminsterfullerene, C$_{60}$. These laboratory studies also attest to the importance of top-down synthesis routes for chemical complexity in space.

Techniques and Review of Absolute Flux Calibration from the Ultraviolet to the Mid-Infrared

Abstract: The measurement of precise absolute fluxes for stellar sources has been pursued with increased vigor since the discovery of the dark energy and the realization that its detailed understanding requires accurate spectral energy distributions (SEDs) of redshifted Ia supernovae in the rest frame. The flux distributions of spectrophotometric standard stars were initially derived from the comparison of stars to laboratory sources of known flux but are now mostly based on calculated model atmospheres. For example, pure hydrogen white dwarf (WD) models provide the basis for the HST CALSPEC archive of flux standards. The basic equations for quantitative spectrophotometry and photometry are explained in detail. Several historical lab based flux calibrations are reviewed; and the SEDs of stars in the major on-line astronomical databases are compared to the CALSPEC reference standard spectrophotometry. There is good evidence that relative fluxes from the visible to the near-IR wavelength of ~2.5 micron are currently accurate to 1% for the primary reference standards; and new comparisons with lab flux standards show promise for improving that precision.

New Vacuum Solar Telescope and Observations with High Resolution

Abstract: The New Vacuum Solar Telescope (NVST) is a 1 meter vacuum solar telescope that aims to observe the fine structures on the Sun. The main tasks of NVST are high resolution imaging and spectral observations, including the measurements of solar magnetic field. NVST is the primary ground-based facility of Chinese solar community in this solar cycle. It is located by the Fuxian Lake of southwest China, where the seeing is good enough to perform high resolution observations. In this paper, we first introduce the general conditions of Fuxian Solar Observatory and the primary science cases of NVST. Then, the basic structures of this telescope and instruments are described in detail. Finally, some typical high resolution data of solar photosphere and chromosphere are also shown.

Determining stellar atmospheric parameters and chemical abundances of FGK stars with iSpec

Abstract: Context. An increasing number of high-resolution stellar spectra is available today thanks to many past and ongoing extensive spectroscopic surveys. Consequently, the scientific community needs automatic procedures to derive atmospheric parameters and individual element abundances. Aims. Based on the widely known SPECTRUM code by R. O. Gray, we developed an integrated spectroscopic software framework suitable for the determination of atmospheric parameters (i.e., effective temperature, surface gravity, metallicity) and individual chemical abundances. The code, named iSpec and freely distributed, is written mainly in Python and can be used on different platforms. Methods. iSpec can derive atmospheric parameters by using the synthetic spectral fitting technique and the equivalent width method. We validated the performance of both approaches by developing two different pipelines and analyzing the Gaia FGK benchmark stars spectral library. The analysis was complemented with several tests designed to assess other aspects, such as the interpolation of model atmospheres and the performance with lower quality spectra. Results. We provide a code ready to perform automatic stellar spectral analysis. We successfully assessed the results obtained for FGK stars with high-resolution and high signal-to-noise spectra.

Non-linear Kalman filters for calibration in radio interferometry

Abstract: We present a new calibration scheme based on a non-linear version of Kalman filter that aims at estimating the physical terms appearing in the Radio Interferometry Measurement Equation (RIME). We enrich the filter's structure with a tunable data representation model, together with an augmented measurement model for regularization. We show using simulations that it can properly estimate the physical effects appearing in the RIME. We found that this approach is particularly useful in the most extreme cases such as when ionospheric and clock effects are simultaneously present. Combined with the ability to provide prior knowledge on the expected structure of the physical instrumental effects (expected physical state and dynamics), we obtain a fairly cheap algorithm that we believe to be robust, especially in low signal-to-noise regime. Potentially the use of filters and other similar methods can represent an improvement for calibration in radio interferometry, under the condition that the effects corrupting visibilities are understood and analytically stable. Recursive algorithms are particularly well adapted for pre-calibration and sky model estimate in a streaming way. This may be useful for the SKA-type instruments that produce huge amounts of data that have to be calibrated before being averaged.