Identifying frequencies with low signal-to-noise ratios in time series of stellar photometry and spectroscopy, and measuring their amplitude ratios and peak widths accurately, are critical goals for asteroseismology. These are also challenges for time series with gaps or whose data are not sampled at a constant rate, even with modern Discrete Fourier Transform (DFT) software. Also the False-Alarm Probability introduced by Lomb and Scargle is an approximation which becomes less reliable in time series with longer data gaps. A rigorous statistical treatment of how to determine the significance of a peak in a DFT, called SigSpec, is presented here. SigSpec is based on an analytical solution of the probability that a DFT peak of a given amplitude does not arise from white noise in a non-equally spaced data set. The underlying Probability Density Function (PDF) of the amplitude spectrum generated by white noise can be derived explicitly if both frequency and phase are incorporated into the solution. In this paper, I define and evaluate an unbiased statistical estimator, the "spectral significance", which depends on frequency, amplitude, and phase in the DFT, and which takes into account the time-domain sampling. I also compare this estimator to results from other well established techniques and demonstrate the effectiveness of SigSpec with a few examples of ground- and space-based photometric data, illustratring how SigSpec deals with the effects of noise and time-domain sampling in determining significant frequencies.

Identifying frequencies with low signal-to-noise ratios in time series of stellar photometry and spectroscopy, and measuring their amplitude ratios and peak widths accurately, are critical goals for asteroseismology. These are also challenges for time series with gaps or whose data are not sampled at a constant rate, even with modern Discrete Fourier Transform (DFT) software. Also the False-Alarm Probability introduced by Lomb and Scargle is an approximation which becomes less reliable in time series with longer data gaps. A rigorous statistical treatment of how to determine the significance of a peak in a DFT, called SigSpec, is presented here. SigSpec is based on an analytical solution of the probability that a DFT peak of a given amplitude does not arise from white noise in a non-equally spaced data set. The underlying Probability Density Function (PDF) of the amplitude spectrum generated by white noise can be derived explicitly if both frequency and phase are incorporated into the solution. In this paper, I define and evaluate an unbiased statistical estimator, the "spectral significance", which depends on frequency, amplitude, and phase in the DFT, and which takes into account the time-domain sampling. I also compare this estimator to results from other well established techniques and demonstrate the effectiveness of SigSpec with a few examples of ground- and space-based photometric data, illustratring how SigSpec deals with the effects of noise and time-domain sampling in determining significant frequencies.

This paper deals with the theoretical principle and optical design of a phase-shifting telescope-interferometer. What is called a “Telescope-Interferometer” (T-I) is indeed a novel, recently proposed Wavefront Error (WFE) sensing technique, whose basic idea consists in combining the main pupil of a telescope with a second, off-axis reference arm. Then a weak modulation of the Point Spread Function (PSF) is generated at the focal plane, allowing for direct phase measurements. We propose a notable improvement of the method, inspired from classical principles of phase shifting interferometry. Herein are presented the alternative principle and its achievable measurement accuracy. The technique shows high performance excepted on narrow areas located near the pupil boundary. It is applicable to both ground or space telescopes and is suitable for the co-phasing of segmented mirrors, which is of prime importance in view of future giant telescope projects.

This paper presents the basic principle and theoretical relationships of an original method allowing to retrieve the Wavefront Errors (WFE) of a ground or space-borne telescope when combining its main pupil with a second, decentered reference optical arm. The measurement accuracy of such a “telescope-interferometer” is then estimated by means of various numerical simulations, demonstrating a high performance excepted on limited areas near the telescope pupil rim. In particular, it allows direct phase evaluation (thus avoiding the use of first or second-order derivatives), which is of special interest for the co-phasing of segmented mirrors in future giant telescopes projects. We finally define the useful practical domain of the method, which seems to be better suited for periodical diagnostics of space or ground based telescopes, or to real-time scientific observations in some very specific cases (e.g. the central star in extrasolar planets searching instruments).

The use of Wavefront Sensors (WFS) is nowadays fundamental in the field of instrumental optics. This paper discusses the principle of an original and recently proposed new class of WFS. Their principle consists in evaluating the slopes of the wavefront errors by means of varying density filters placed into the image plane of the tested optical system. The device, sometimes called ‘optical differentiation WFS’ is completed by a digital data-processing system reconstructing the wavefront from the obtained slopes. Various luminous sources of different wavelengths and spectral widths can be employed. The capacities of the method are discussed from the geometrical and Fourier optics points of view, then by means of numerical simulations showing that the ultimate accuracy can be well below lambda/10 and lambda/100 Peak-to-Valley (PTV) and RMS respectively, provided that certain precautions are taken.

We have developed a characterization of the geological evolution of the Earths atmosphere and surface in order to model the observable spectra of an Earth-like planet through its geological history. These calculations are designed to guide the interpretation of an observed spectrum of such a planet by future instruments that will characterize exoplanets. Our models focus on spectral features that either imply habitability or are required for habitability. These features are generated by H2O, CO2, CH4, O2, O3, N2O, and vegetation-like surface albedos. We chose six geological epochs to characterize. These epochs exhibit a wide range in abundance for these molecules, ranging from a CO2 rich early atmosphere, to a CO2/CH4-rich atmosphere around 2 billion years ago to a present-day atmosphere. We analyzed the spectra to quantify the strength of each important spectral feature in both the visible and thermal infrared spectral regions, and the resolutions required to unambiguously observe the features for each epoch. We find a wide range of spectral resolutions required for observing the different features. For example, H2O and O3 can be observed with relatively low resolution, while O2 and N2O require higher resolution. We also find that the inclusion of clouds in our models significantly affects both the strengths and resolutions required to observe all spectral features.

The requirements on space missions designed to study Terrestrial exoplanets are discussed. We then investigate whether the design of such a mission, specifically the Darwin nulling interferometer, can be carried out in a simplified scenario. The key element here is accepting somewhat higher levels of stellar leakage. We establish detailed requirements resulting from the scientific rationale for the mission, and calculate detailed parameters for the stellar suppression required to achieve those requirements. We do this utilizing the Darwin input catalogue. The dominating noise source for most targets in this sample is essentially constant for all targets, while the leakage diminishes with the square of the distance. This means that the stellar leakage has an effect on the integration time only for the nearby stars, while for the more distant targets its influence decreases significantly. We assess the impact of different array configurations and nulling profiles and identify the stars for which the detection efficiency can be maximized.

The Darwin and TPF-I missions are Infrared free flying interferometer missions based on nulling interferometry. Their main objective is to detect and characterize other Earth-like planets, analyze the composition of their atmospheres and their capability to sustain life, as we know it. Darwin and TPF-I are currently in study phase. A number of mission architectures of 3 and 4 free flying telescopes are evaluated on the basis of the interferometer’s response, ability to distinguish multiple planet signatures and starlight rejection capabilities. The characteristics of the new configurations are compared also to the former, more complex Bowtie baseline architectures as well as evaluated on base of their science capability.

The search for life on extrasolar planets is based on the assumption that one can screen extrasolar planets for habitability spectroscopically. The first space born instruments able to detect as well as characterize extrasolar planets, Darwin and terrestrial planet finder (TPF-I and TPF-C) are scheduled to launch before the end of the next decade. The composition of the planetary surface, atmosphere, and its temperature-pressure profile influence a detectable spectroscopic signal considerably. For future space-based missions it will be crucial to know this influence to interpret the observed signals and detect signatures of life in remotely observed atmospheres. We give an overview of biomarkers in the visible and IR range, corresponding to the TPF-C and TPF-I/DARWIN concepts, respectively. We also give an overview of the evolution of biomarkers over time and its implication for the search for life on extrasolar Earth-like planets. We show that atmospheric features on Earth can provide clues of biological activities for at least 2 billion years.

The DARWIN mission is an Infrared free flying interferometer mission based on the new technique of nulling interferometry. Its main objective is to detect and characterize other Earth-like planets, analyze the composition of their atmospheres and their capability to sustain life, as we know it. DARWIN is currently in definition phase. This PhD work that has been undertaken within the DARWIN team at the European Space Agency (ESA) addresses two crucial aspects of the mission. Firstly, a DARWIN target star list has been established that includes characteristics of the target star sample that will be critical for final mission design, such as, luminosity, distance, spectral classification, stellar variability, multiplicity, location and radius of the star. Constrains were applied as set by planet evolution theory and mission architecture. Secondly, a number of alternative mission architectures have been evaluated on the basis of interferometer response as a function of wavelength, achievable modulation efficiency, number of telescopes and starlight rejection capabilities. The study has shown that the core mission goals should be achievable with a lower level of complexity as compared to the current baseline configuration.

We report on the feasibility study of a Laue lens for hard X-rays (>60 keV)based on mosaic crystals, for astrophysical applications. In particular we discuss the scientific motivations, its functioning principle, the procedure followed to select the suitable crystal materials, the criteria adopted to establish crystal dimensions and their distribution on the lens in order to obtain the best lens focusing capabilities, and the criteria for optimizing the lens effective area in a given passband. We also discuss the effects of misalignments of the crystal tiles due to unavoidable mechanical errors in assembling the lens. A software was developed to face all these topics and to evaluate the expected lens performance.

The Antarctic Muon And Neutrino Detector Array (AMANDA) is a high-energy neutrino telescope operating at the geographic South Pole. It is a lattice of photo-multiplier tubes buried deep in the polar ice between 1500m and 2000m. The primary goal of this detector is to discover astrophysical sources of high energy neutrinos. A high-energy muon neutrino coming through the earth from the Northern Hemisphere can be identified by the secondary muon moving upward through the detector. The muon tracks are reconstructed with a maximum likelihood method. It models the arrival times and amplitudes of Cherenkov photons registered by the photo-multipliers. This paper describes the different methods of reconstruction, which have been successfully implemented within AMANDA. Strategies for optimizing the reconstruction performance and rejecting background are presented. For a typical analysis procedure the direction of tracks are reconstructed with about 2 degree accuracy.

The Stratospheric Observatory for Infrared Astronomy (SOFIA) completed its first light flight in May of 2010 using the facility mid-infrared instrument FORCAST. Since then, FORCAST has successfully completed thirteen science flights on SOFIA. In this paper we describe the design, operation and performance of FORCAST as it relates to the initial three Short Science flights. FORCAST was able to achieve near diffraction-limited images for lambda > 30 microns allowing unique science results from the start with SOFIA. We also describe ongoing and future modifications that will improve overall capabilities and performance of FORCAST.

The historical development of ground based astronomical telescopes leads us to expect that space-based astronomical telescopes will need to be operational for many decades. The exchange of scientific instruments in space will be a prerequisite for the long lasting scientific success of such missions. Operationally, the possibility to repair or replace key spacecraft components in space will be mandatory. We argue that these requirements can be fulfilled with robotic missions and see the development of the required engineering as the main challenge. Ground based operations, scientifically and technically, will require a low operational budget of the running costs. These can be achieved through enhanced autonomy of the spacecraft and mission independent concepts for the support of the software. This concept can be applied to areas where the mirror capabilities do not constrain the lifetime of the mission.

We have investigated the possibility of measuring the electron neutrino mass with sub-eV sensitivity by studying the electron capture decay of 163-Ho with cryogenic microcalorimeters. In this paper we will introduce an experiment’s concept, discuss the technical requirements, and identify a roadmap to reach a sensitivity of 0.1 eV and beyond.

Southern Africa has some of the world’s best sites for air Cherenkov telescopes. South Africa has only one viable site, which is south of Sutherland and also close to the Southern African Large Telescope (SALT). This site has very good infrastructure and is easy to access, but only 47% of the night-time has a cloudless sky usable for observations. Namibia, which already hosts the H.E.S.S telescope, has a number of potential sites with much less cloud coverage. The H.E.S.S. site is one of the highest of these sites at 1840 m a.s.l. with about 64% of the night-time cloudless. It also has very low night sky background levels and is relatively close (about 100 km) to Windhoek. Moving further away from Windhoek to the south, the cloud coverage and artificial night sky brightness becomes even less, with the site at Kuibis (between Keetmanshoop and Luderitz) at 1640 m a.s.l. having clear night skies 73% of the time. Even though this site seems remote (being 660 km from Windhoek by road), it is close to the national B4 highway, a railway line, a power line and an optical fiber line. It is also less than two hours drive away from a harbour and national airports. The Namibian sites also receive very little snow, if any, and the wind speeds are less than 50 km/h for more than 90% of the time with maximum wind speeds of around 100 km/h. Seismically the whole Southern African region is very stable.

Accretion onto supermassive black holes (SMBHs) in galaxy formation simulations is frequently modelled by the Bondi-Hoyle formalism. Here we examine the validity of this approach analytically and numerically. We argue that the character of the flow where one evaluates the gas properties is unlikely to satisfy the simple Bondi-Hoyle model. Only in the specific case of hot virialised gas with zero angular momentum and negligible radiative cooling is the Bondi-Hoyle solution relevant. In the opposite extreme, where the gas is in a state of free-fall at the evaluation radius due to efficient cooling and the dominant gravity of the surrounding halo, the Bondi-Hoyle formalism can be erroneous by orders of magnitude in either direction. This may impose artificial trends with halo mass in cosmological simulations by being wrong by different factors for different halo masses. We propose an expression for the sub-grid accretion rate which interpolates between the free-fall regime and the Bondi-Hoyle regime, therefore taking account of the contribution of the halo to the gas dynamics.

In this short communication I highlight how the number of collaborators on papers in the main astronomy journals has evolved over time. We see a trend of moving away from single-author papers. This communication is based on data in the holdings of the SAO/NASA Astrophysics Data System (ADS). The ADS is funded by NASA Grant NNX09AB39G.

We consider the task of distinguishing between two different alternative models that can roughly equally explain observed time series data, mainly focusing on the period ambiguity case (aliasing). We propose a test for checking whether the rival models are observationally equivalent for now or they are already distinguishable. It is the Vuong closeness test, which is based on the Kullback-Leibler Information Criterion. It is asymptotically normal and can work (in certain sense) even in the misspecified case, when the both proposed alternatives are actually wrong. This test is also very simple for practical use. We apply it to several known extrasolar planetary systems and find that our method often helps to resolve various model ambiguities emerging in astronomical practice, but preventing us from hasty conclusions in other cases.

We present a detailed analysis of four of the most widely used radio source finding packages in radio astronomy, and a program being developed for the Australian Square Kilometer Array Pathfinder (ASKAP) telescope. The four packages; SExtractor, SFind, IMSAD and Selavy are shown to produce source catalogues with high completeness and reliability. In this paper we analyse the small fraction (~1%) of cases in which these packages do not perform well. This small fraction of sources will be of concern for the next generation of radio surveys which will produce many thousands of sources on a daily basis, in particular for blind radio transients surveys. From our analysis we identify the ways in which the underlying source finding algorithms fail. We demonstrate a new source finding algorithm Aegean, based on the application of a Laplacian kernel, which can avoid these problems and can produce complete and reliable source catalogues for the next generation of radio surveys.

We present new mid-infrared images of the central region of the Orion Nebula using the newly commissioned SOFIA airborne telescope and its 5 — 40 micron camera FORCAST. The 37.1 micron images represent the highest resolution observations (<4") ever obtained of this region at these wavelengths. After BN/KL (which is described in a separate letter in this issue), the dominant source at all wavelengths except 37.1 micron is the Ney-Allen Nebula, a crescent-shaped extended source associated with theta 1D. The morphology of the Ney-Allen nebula in our images is consistent with the interpretation that it is ambient dust swept up by the stellar wind from theta 1D, as suggested by Smith et al. (2005). Our observations also reveal emission from two "proplyds" (proto-planetary disks), and a few embedded young stellar objects (YSOs; IRc9, and OMC1S IRS1, 2, and 10). The spectral energy distribution for IRc9 is presented and fitted with standard YSO models from Robitaille et al. (2007) to constrain the total luminosity, disk size, and envelope size. The diffuse, nebular emission we observe at all FORCAST wavelengths is most likely from the background photodissociation region (PDR) and shows structure that coincides roughly with H_alpha and [N II] emission. We conclude that the spatial variations in the diffuse emission are likely due to undulations in the surface of the background PDR.

Pioneering efforts aiming at the development of multi-messenger gravitational wave and electromagnetic astronomy have been made. An electromagnetic observation follow-up program of candidate gravitational wave events has been performed (Dec 17 2009 to Jan 8 2010 and Sep 4 to Oct 20 2010) during the recent runs of the LIGO and Virgo gravitational wave detectors. It involved ground-based and space electromagnetic facilities observing the sky at optical, X-ray and radio wavelengths. The joint gravitational wave and electromagnetic observation study requires the development of specific image analysis procedures able to discriminate the possible electromagnetic counterpart of gravitational wave triggers from contaminant/background events. The paper presents an overview of the electromagnetic follow-up program and the image analysis procedures.

2011 was a successful year for the GREGOR project. The telescope was finally completed in May with the installation of the 1.5-meter primary mirror. The installation of the first-light focal plane instruments was completed by the end of the year. At the same time, the preparations for the installation of the high-order adaptive optics were finished, its integration to the telescope is scheduled for early 2012. This paper describes the telescope and its instrumentation in their present first light configuration, and provides a brief overview of the science goals of GREGOR.

In the standard formulation of the smoothed particle hydrodynamics (SPH), it is assumed that the local density distribution is differentiable. This assumption is used to derive the spatial derivatives of other quantities. However, this assumption breaks down at the contact discontinuity, which appears often in simulations of astronomical objects. At the contact discontinuity, the density of the low-density side is overestimated while that of the high-density side is underestimated. As a result, the pressure of the low (high) density side is over (under) estimated. Thus, unphysical repulsive force appears at the contact discontinuity, resulting in the effective surface tension. This effective surface tension suppresses instabilities such as the Kelvin-Helmholtz and Rayleigh-Taylor instabilities. In this paper, we present a new formulation of SPH, which does not require the differentiability of density and thus can handle contact discontinuity without numerical problems. The results of standard tests such as the shock tube, Kelvin-Helmholtz and Rayleigh-Taylor instabilities, and the blob tests are all very favorable to our new formulation. We conclude that our new formulation solved practically all known difficulties of the standard SPH, without introducing additional numerical diffusion or breaking the exact force symmetry or energy conservation.

In the standard formulation of the smoothed particle hydrodynamics (SPH), it is assumed that the local density distribution is differentiable. This assumption is used to derive the spatial derivatives of other quantities. However, this assumption breaks down at the contact discontinuity, which appears often in simulations of astronomical objects. At the contact discontinuity, the density of the low-density side is overestimated while that of the high-density side is underestimated. As a result, the pressure of the low (high) density side is over (under) estimated. Thus, unphysical repulsive force appears at the contact discontinuity, resulting in the effective surface tension. This effective surface tension suppresses instabilities such as the Kelvin-Helmholtz and Rayleigh-Taylor instabilities. In this paper, we present a new formulation of SPH, which does not require the differentiability of density and thus can handle contact discontinuity without numerical problems. The results of standard tests such as the shock tube, Kelvin-Helmholtz and Rayleigh-Taylor instabilities, and the blob tests are all very favorable to our new formulation. We conclude that our new formulation solved practically all known difficulties of the standard SPH, without introducing additional numerical diffusion or breaking the exact force symmetry or energy conservation.

Since long time it has been pointed out by the DAMA collaboration that muons surviving the Gran Sasso mountain cannot mimic the Dark Matter annual modulation signature exploited by the DAMA/NaI and DAMA/LIBRA experiments. In the present paper, these and further arguments are gathered together in order to enable a wider community to suitably approach this point.

Expanding upon Pimbblet’s informative 2011 analysis of career h-indices for members of the Astronomical Society of Australia, we provide additional citation metrics which are geared to a) quantifying the current performance of b) all professional astronomers in Australia. We have trawled the staff web-pages of Australian Universities, Observatories and Research Organisations hosting professional astronomers, and identified 383 PhD-qualified, research-active, astronomers in the nation – 131 of these are not members of the Astronomical Society of Australia. Using the SAO/NASA Astrophysics Data System, we provide the three following common metrics based on publications in the first decade of the 21st century (2001-2010): h-index, author-normalised citation count and lead-author citation count. We additionally present a somewhat more inclusive analysis, applicable for many early-career researchers, that is based on publications from 2006-2010. Histograms and percentiles, plus top-performer lists, are presented for each category. Finally, building on Hirsch’s empirical equation, we find that the (10-year) h-index and (10-year) total citation count T can be approximated by the relation h = (0.5+sqrt{T})/sqrt{5} for h > 5.

Expanding upon Pimbblet’s informative 2011 analysis of career h-indices for members of the Astronomical Society of Australia, we provide additional citation metrics which are geared to a) quantifying the current performance of b) all professional astronomers in Australia. We have trawled the staff web-pages of Australian Universities, Observatories and Research Organisations hosting professional astronomers, and identified 383 PhD-qualified, research-active, astronomers in the nation – 131 of these are not members of the Astronomical Society of Australia. Using the SAO/NASA Astrophysics Data System, we provide the three following common metrics based on publications in the first decade of the 21st century (2001-2010): h-index, author-normalised citation count and lead-author citation count. We additionally present a somewhat more inclusive analysis, applicable for many early-career researchers, that is based on publications from 2006-2010. Histograms and percentiles, plus top-performer lists, are presented for each category. Finally, building on Hirsch’s empirical equation, we find that the (10-year) h-index and (10-year) total citation count T can be approximated by the relation h = (0.5+sqrt{T})/sqrt{5} for h > 5.

Expanding upon Pimbblet’s informative 2011 analysis of career h-indices for members of the Astronomical Society of Australia, we provide additional citation metrics which are geared to a) quantifying the current performance of b) all professional astronomers in Australia. We have trawled the staff web-pages of Australian Universities, Observatories and Research Organisations hosting professional astronomers, and identified 383 PhD-qualified, research-active, astronomers in the nation – 131 of these are not members of the Astronomical Society of Australia. Using the SAO/NASA Astrophysics Data System, we provide the three following common metrics based on publications in the first decade of the 21st century (2001-2010): h-index, author-normalised citation count and lead-author citation count. We additionally present a somewhat more inclusive analysis, applicable for many early-career researchers, that is based on publications from 2006–2010. Histograms and percentiles, plus top-performer lists, are presented for each category. Finally, building on Hirsch’s empirical equation, we find that the (10-year) h-index and total citation count T can be approximated by the relation h = (0.5+\sqrt{T})/\sqrt{5} for h > 5.

A Bayesian approach to calibrating period-luminosity (PL) relations has substantial benefits over generic least-squares fits. In particular, the Bayesian approach takes into account the full prior distribution of the model parameters, such as the a priori distances, and refits these parameters as part of the process of settling on the most highly-constrained final fit. Additionally, the Bayesian approach can naturally ingest data from multiple wavebands and simultaneously fit the parameters of PL relations for each waveband in a procedure that constrains the parameter posterior distributions so as to minimize the scatter of the final fits appropriately in all wavebands. Here we describe the generalized approach to Bayesian model fitting and then specialize to a detailed description of applying Bayesian linear model fitting to the mid-infrared PL relations of RR Lyrae variable stars. For this example application we quantify the improvement afforded by using a Bayesian model fit. We also compare distances previously predicted in our example application to recently published parallax distances measured with the Hubble Space Telescope and find their agreement to be a vindication of our methodology. Our intent with this article is to spread awareness of the benefits and applicability of this Bayesian approach and encourage future PL relation investigations to consider employing this powerful analysis method.

The ANTARES neutrino telescope, located 40km off the coast of Toulon in the Mediterranean Sea at a mooring depth of about 2475m, consists of twelve detection lines equipped typically with 25 storeys. Every storey carries three optical modules that detect Cherenkov light induced by charged secondary particles (typically muons) coming from neutrino interactions. As these lines are flexible structures fixed to the sea bed and held taut by a buoy, sea currents cause the lines to move and the storeys to rotate. The knowledge of the position of the optical modules with a precision better than 10cm is essential for a good reconstruction of particle tracks. In this paper the ANTARES positioning system is described. It consists of an acoustic positioning system, for distance triangulation, and a compass-tiltmeter system, for the measurement of the orientation and inclination of the storeys. Necessary corrections are discussed and the results of the detector alignment procedure are described.

We propose a revolutionary way of studying the surface of Mars using a wind-driven network of mobile sensors- Gone with the Wind ON_Mars (GOWON). GOWON is a scalable architecture that will allow in-situ mapping of a wide range of phenomena, exploiting existing capabilities, but radically improving our ability to study Mars. GOWON has the following characteristics: 1.it consists of a dynamic wireless network of many compact mobile sensors. 2.the mobile sensors (called moballs) are spherically-shaped and wind-driven; they are lightweight and bouncy. 3. moballs communicate with each other and earth through a satellite system orbiting Mars. There is also peer-to-peer communication between the moballs, creating a network of shared data, computing, and tasks. Motivation and Rationale Thanks to earlier exploration missions to Mars we now have a much better understanding of many of the natural characteristics of the red planet. We now know that there is an abundance of wind (with average speeds of 10 m/s and much higher maximum speeds [1]), dust storms, high levels of saltation [2], crustal magnetic fields generated by an abundance of ferromagnetic minerals and rocks [3], etc. Future Mars missions must therefore attempt to leverage these characteristics, and must do so by exploiting recent advances in low power micro-devices using MEMS (Micro-Electro-Mechanical Systems that harvest vibration and other movement into energy) technologies and others, miniature cameras, miniature wet chemistry labs, integrated circuits, low power wireless devices, etc. We believe the system proposed here addresses this opportunity heads on. Our proposed system is much more than a collection of sensors; the system is larger than the sum of its parts. In addition to communicating with the satellite, the moballs can communicate with each other and therefore distribute tasks, data, computations, etc.

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.

The EMCCD is a type of CCD that delivers fast readout times and negligible readout noise, making it an ideal detector for high frame rate applications which improve resolution, like lucky imaging or shift-and-add. This improvement in resolution can potentially improve the photometry of faint stars in extremely crowded fields significantly by alleviating crowding. This improvement in resolution is a prerequisite for observing gravitational microlensing in main sequence stars towards the galactic bulge. However, the photometric stability of this device has not been assessed. The EMCCD has sources of noise not found in conventional CCDs, and new methods for photometric reduction must be developed. We aim to investigate how the normal photometric reduction steps from conventional CCDs should be adjusted to be applicable to EMCCD data. One complication is that a bias frame cannot be obtained conventionally, as the output from an EMCCD is not normally distributed. Also, the readout process generates spurious charges in any CCD, but in EMCCD data, these charges are visible as opposed to the conventional CCD. A simple probabilistic model for the dark output of an EMCCD is developed. Fitting this model with the expectation-maximization algorithm allows us to estimate the bias, readout noise, amplification, and spurious charge rate per pixel and thus correct for these phenomena. To investigate the stability of the photometry, corrected frames of a crowded field are reduced with a PSF fitting photometry package. We find that it is possible to develop an algorithm that elegantly reduces EMCCD data and produces stable photometry at the 1% level in an extremely crowded field.

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 http://danfm.ca/emcee under the GNU General Public License v2.

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 http://danfm.ca/emcee under the GNU General Public License v2.

Hipparcos, the first ever experiment of global astrometry, was launched by ESA in 1989 and its results published in 1997 (Perryman et al., Astron. Astrophys. 323, L49, 1997; Perryman & ESA (eds), The Hipparcos and Tycho catalogues, ESA SP-1200, 1997). A new reduction was later performed using an improved satellite attitude reconstruction leading to an improved accuracy for stars brighter than 9th magnitude (van Leeuwen & Fantino, Astron. Astrophys. 439, 791, 2005; van Leeuwen, Astron. Astrophys. 474, 653, 2007). The Hipparcos Catalogue provided an extended dataset of very accurate astrometric data (positions, trigonometric parallaxes and proper motions), enlarging by two orders of magnitude the quantity and quality of distance determinations and luminosity calibrations. The availability of more than 20000 stars with a trigonometric parallax known to better than 10% opened the way to a drastic revision of our 3-D knowledge of the solar neighbourhood and to a renewal of the calibration of many distance indicators and age estimations. The prospects opened by Gaia, the next ESA cornerstone, planned for launch in June 2013 (Perryman et al., Astron. Astrophys. 369, 339, 2001), are still much more dramatic: a billion objects with systematic and quasi simultaneous astrometric, spectrophotometric and spectroscopic observations, about 150 million stars with expected distances to better than 10%, all over the Galaxy. All stellar distance indicators, in very large numbers, will be directly measured, providing a direct calibration of their luminosity and making possible detailed studies of the impacts of various effects linked to chemical element abundances, age or cluster membership. With the help of simulations of the data expected from Gaia, obtained from the mission simulator developed by DPAC, we will illustrate what Gaia can provide with some selected examples.

The massive star forming region W3 was observed with the faint object infrared camera for the SOFIA telescope (FORCAST) as part of the Short Science program. The 6.4, 6.6, 7.7, 19.7, 24.2, 31.5 and 37.1 \um bandpasses were used to observe the emission of Polycyclic Aromatic Hydrocarbon (PAH) molecules, Very Small Grains and Big Grains. Optical depth and color temperature maps of W3A show that IRS2 has blown a bubble devoid of gas and dust of $\sim$0.05 pc radius. It is embedded in a dusty shell of ionized gas that contributes 40% of the total 24 \um emission of W3A. This dust component is mostly heated by far ultraviolet, rather than trapped Ly$\alpha$ photons. This shell is itself surrounded by a thin ($\sim$0.01 pc) photodissociation region where PAHs show intense emission. The infrared spectral energy distribution (SED) of three different zones located at 8, 20 and 25\arcsec from IRS2, show that the peak of the SED shifts towards longer wavelengths, when moving away from the star. Adopting the stellar radiation field for these three positions, DUSTEM model fits to these SEDs yield a dust-to-gas mass ratio in the ionized gas similar to that in the diffuse ISM. However, the ratio of the IR-to-UV opacity of the dust in the ionized shell is increased by a factor $\simeq$3 compared to the diffuse ISM.

We present the last results of our PSF reconstruction (PSF-R) project for the Keck-II and Gemini-North AO systems in natural guide star mode. Our initial tests have shown that the most critical aspects of PSF-R are the determination of the system static aberrations and the optical turbulence parameters, and we have set up a specific observation campaign on the two systems to explore this. We demonstrate that deformable mirror based seeing monitor works well, and 10% accuracy is easily obtained. Phase diversity has been demonstrated to work on sky sources. Besides, residual phase stationarity is an important assumption in PSF-R, and we demonstrate here that it is basically true. As a result of these tests and verifications, we have been able for the first time to obtain a very good PSF reconstruction for the Keck-II system, in bright natural guide star mode.

We present the last results of our PSF reconstruction (PSF-R) project for the Keck-II and Gemini-North AO systems in natural guide star mode. Our initial tests have shown that the most critical aspects of PSF-R are the determination of the system static aberrations and the optical turbulence parameters, and we have set up a specific observation campaign on the two systems to explore this. We demonstrate that deformable mirror based seeing monitor works well, and 10% accuracy is easily obtained. Phase diversity has been demonstrated to work on sky sources. Besides, residual phase stationarity is an important assumption in PSF-R, and we demonstrate here that it is basically true. As a result of these tests and verifications, we have been able for the first time to obtain a very good PSF reconstruction for the Keck-II system, in bright natural guide star mode.

Directional detection of Dark Matter is a promising search strategy. However, to perform such detection, a given set of parameters has to be retrieved from the recoiling tracks : direction, sense and position in the detector volume. In order to optimize the track reconstruction and to fully exploit the data of forthcoming directional detectors, we present a likelihood method dedicated to 3D track reconstruction. This new analysis method is applied to the MIMAC detector. It requires a full simulation of track measurements in order to compare real tracks to simulated ones. We conclude that a good spatial resolution can be achieved, i.e. sub-mm in the anode plane and cm along the drift axis. This opens the possibility to perform a fiducialization of directional detectors. The angular resolution is shown to range between 20$^\circ$ to 80$^\circ$, depending on the recoil energy, which is however enough to achieve a high significance discovery of Dark Matter. On the contrary, we show that sense recognition capability of directional detectors depends strongly on the recoil energy and the drift distance, with small efficiency values (50%-70%). We suggest not to consider this information either for exclusion or discovery of Dark Matter for recoils below 100 keV and then to focus on axial directional data.

The Very Large Telescope (VLT) of the European Southern Observatory (ESO) on Cerro Paranal is one of the most influential observing complexes in the world. In this overview lecture, an introduction to the ESO organization and VLT telescopes is presented, along with a small selection of scientific works. Particular attention is given to the importance of the VLT for the Russian scientific community, as well as to present and future perspectives for making use of ESO and VLT data.

The Very Large Telescope (VLT) of the European Southern Observatory (ESO) on Cerro Paranal is one of the most influential observing complexes in the world. In this overview lecture, an introduction to the ESO organization and VLT telescopes is presented, along with a small selection of scientific works. Particular attention is given to the importance of the VLT for the Russian scientific community, as well as to present and future perspectives for making use of ESO and VLT data.

In this paper we demonstrate a methodology to remove the power of the drift induced from random acceleration on LISA proof mass in the frequency domain. The drift must be cleaned from LISA time series data in advance of any further analysis. The cleaning is usually performed in the time domain by using a quadratic function to fit the time series data, and then removing the fitted part from the data. Having Fourier transformed the residuals, and then convolved with LISA transfer function, LISA sensitivity curve can be obtained. However, cosmic gravitational-wave background cannot be retrieved with this approach due to its random nature. Here we provide a new representation of power spectrum given by discrete Fourier transform, which is applied to find the function of the drift power for the cleaning in the frequency domain. We also give the probability distribution used to analyze the data in the frequency domain. We combine several techniques, including Markov Chain Monte Carlo method, simulated annealing, and Gelman & Rubin’s method, with Baye’s theorem to build the algorithm. The algorithm is utilized to analyze 24 simulations of LISA instrumental noise. We prove that the LISA sensitivity can be recovered through this approach. It can help us to build algorithms for some tasks which are must accomplished in the frequency domain for LISA data analysis. This method can be applied to other space-borne interferometers if charges on their proof masses cannot be perfectly cancelled.

We review current methods for building PSF-matching kernels for the purposes of image subtraction or coaddition. Such methods use a linear decomposition of the kernel on a series of basis functions. The correct choice of these basis functions is fundamental to the efficiency and effectiveness of the matching – the chosen bases should represent the underlying signal using a reasonably small number of shapes, and/or have a minimum number of user-adjustable tuning parameters. We examine methods whose bases comprise multiple Gauss-Hermite polynomials, as well as a form free basis composed of delta-functions. Kernels derived from delta-functions are unsurprisingly shown to be more expressive; they are able to take more general shapes and perform better in situations where sum-of-Gaussian methods are known to fail. However, due to its many degrees of freedom (the maximum number allowed by the kernel size) this basis tends to overfit the problem, and yields noisy kernels having large variance. We introduce a new technique to regularize these delta-function kernel solutions, which bridges the gap between the generality of delta-function kernels, and the compactness of sum-of-Gaussian kernels. Through this regularization we are able to create general kernel solutions that represent the intrinsic shape of the PSF-matching kernel with only one degree of freedom, the strength of the regularization lambda. The role of lambda is effectively to exchange variance in the resulting difference image with variance in the kernel itself. We examine considerations in choosing the value of lambda, including statistical risk estimators and the ability of the solution to predict solutions for adjacent areas. Both of these suggest moderate strengths of lambda between 0.1 and 1.0, although this optimization is likely dataset dependent.

A number of well-motivated extensions of the LCDM concordance cosmological model postulate the existence of a population of sources embedded in the cosmic microwave background (CMB). One such example is the signature of cosmic bubble collisions which arise in models of eternal inflation. The most unambiguous way to test these scenarios is to evaluate the full posterior probability distribution of the global parameters defining the theory; however, a direct evaluation is computationally impractical on large datasets, such as those obtained by the Wilkinson Microwave Anisotropy Probe (WMAP) and Planck. A method to approximate the full posterior has been developed recently, which requires as an input a set of candidate sources which are most likely to give the largest contribution to the likelihood. In this article, we present an improved algorithm for detecting candidate sources using optimal filters, and apply it to detect candidate bubble collision signatures in WMAP 7-year observations. We show both theoretically and through simulations that this algorithm provides an enhancement in sensitivity over previous methods by a factor of approximately two. Moreover, no other filter-based approach can provide a superior enhancement of these signatures. Applying our algorithm to WMAP 7-year observations, we detect eight new candidate bubble collision signatures for follow-up analysis.

The IceTop air shower array is the surface component of the IceCube Neutrino Observatory at the geographic South Pole. The combination of IceTop and IceCube provides a new and powerful tool to measure cosmic ray composition in the energy range between about 300 TeV and 1 EeV by detecting the electromagnetic component at the surface in coincidence with the muon bundle in the deep underground detector. The paper will give an overview of the current status of the detector and the first physics results will be presented.

Within the efforts to bring frontline interactive astrophysics and astronomy to the classroom, the Hands on Universe (HOU) developed a set of exercises and platform using real data obtained by some of the most advanced ground and space observatories. The backbone of this endeavour is a new free software Web tool – Such a Lovely Software for Astronomy based on Image J (Salsa J). It is student-friendly and developed specifically for the HOU project and targets middle and high schools. It allows students to display, analyze, and explore professionally obtained astronomical images, while learning concepts on gravitational dynamics, kinematics, nuclear fusion, electromagnetism. The continuous evolving set of exercises and tutorials is being completed with real (professionally obtained) data to download and detailed tutorials. The flexibility of the Salsa J platform tool enables students and teachers to extend the exercises with their own observations. The software developed for the HOU program has been designed to be a multi-platform, multi-lingual experience for image manipulation and analysis in the classroom. Its design enables easy implementation of new facilities (extensions and plugins), minimal in-situ maintenance and flexibility for exercise plugin. Here, we describe some of the most advanced exercises about astrophysics in the classroom, addressing particular examples on gravitational dynamics, concepts currently introduced in most sciences curricula in middle and high schools.

In the United States, the National Science Foundation (NSF), anticipating the no growth in funding, has commissioned a review of NSF-funded astronomy assets with the goal of determining how to best allocate funding for this decade. Inputs from members of the US community were sought. It is a matter of simple arithmetic that for a fixed level of funding many significant aspirations of Astro2010 cannot be met. Here, accepting the boundary conditions posed above, I have focused on fields centered on optical astronomy which offer the best opportunity for progress in this decade and thus offer the highest cost-benefit ratio. Readers may profit from reading the first seven sections.

The Hamamatsu R11410-10 PMT has been tested extensively at UCLA and Hamamatsu for use in future ton-scale dual-phase xenon dark matter detectors. The R11410-10 PMT has been shown to have a quantum efficiency of greater than 30% at xenon scintillation wavelengths, a radioactivity of 20.4 mBq/piece and a dark count rate of 50 Hz at 0.3 photoelectron threshold. These results in conjunction with the high gain (up to 10^7), good single photoelectron resolution and stability at low temperature provide a satisfactory photosensor for the future of low background event searches.

The contributions of everyday individuals to significant research has grown dramatically beyond the early days of classical birdwatching and endeavors of amateurs of the 19th century. Now people who are casually interested in science can participate directly in research covering diverse scientific fields. Regarding astronomy, volunteers, either as individuals or as networks of people, are involved in a variety of types of studies. Citizen Science is intuitive, engaging, yet necessarily robust in its adoption of sci-entific principles and methods. Herein, we discuss Citizen Science, focusing on fully participatory projects such as Zooniverse (by several of the au-thors CL, AS, LF, SB), with mention of other programs. In particular, we make the case that citizen science (CS) can be an important aspect of the scientific data analysis pipelines provided to scientists by observatories.

Doppler spectroscopy has uncovered or confirmed all the known planets orbiting nearby stars. Two different approaches are used to obtain precision Doppler measurements at optical wavelengths. The first approach is the gas cell method, which is based on the least-squares matching of the absorption spectrum of Iodine over-imposed to the spectrum of the star. The second method relies on the construction of a stabilized spectrograph calibrated in wavelength with an externally fed calibration source. The most precise stabilized spectrometer in operation is HARPS, operated by ESO in La Silla/Chile. In the case of HARPS, the Doppler measurement is obtained using the so–called Cross-Correlation Function technique (CCF). It consists of multiplying the stellar spectrum with a binary mask and finding the minimum of such product as a function of the stellar Doppler shift. Such mask is weighted to account for the different depths of the stellar lines. It is known that CCF is suboptimal in exploiting the Doppler information in the stellar spectrum. Here, we describe an algorithm to obtain precision RV measurements based on least squares matching of each observation to a high signal-to-noise ratio template. Such algorithm is implemented in our software called HARPS-TERRA (Template Enhanced Radial velocity Re-analysis Application). We show that, compared to CCF, template matching provides a significant improvement in accuracy, specially when applied to M dwarfs. We conclude that other stabilized spectrographs should use a similar approach to achieve the sub \ms precision required to detect potentially habitable worlds around nearby stars.

Detecting the first electron pairs with nuclear emulsion allows a precise measurement of the direction of incident gamma-rays as well as their polarization. With recent innovations in emulsion scanning, emulsion analyzing capability is becoming increasingly powerful. Presently, we are developing a balloon-borne gamma-ray telescope using nuclear emulsion. An overview and a status of our telescope is given.

In general Evolutionary Computation (EC) includes a number of optimization methods inspired by biological mechanisms of evolution. The methods catalogued in this area use the Darwinian principles of life evolution to produce algorithms that returns high quality solutions to hard-to-solve optimization problems. The main strength of EC is precisely that they provide good solutions even if the computational resources (e.g., running time) are limited. Astronomy and Astrophysics are two fields that often require optimizing problems of high complexity or analyzing a huge amount of data and the so-called complete optimization methods are inherently limited by the size of the problem/data. For instance, reliable analysis of large amounts of data is central to modern astrophysics and astronomical sciences in general. EC techniques perform well where other optimization methods are inherently limited (as complete methods applied to NP-hard problems), and in the last ten years, numerous proposals have come up that apply with greater or lesser success methodologies of evolutional computation to common engineering problems. Some of these problems, such as the estimation of non-lineal parameters, the development of automatic learning techniques, the implementation of control systems, or the resolution of multi-objective optimization problems, have had (and have) a special repercussion in the fields. For these reasons EC emerges as a feasible alternative for traditional methods. In this paper, we discuss some promising applications in this direction and a number of recent works in this area; the paper also includes a general description of EC to provide a global perspective to the reader and gives some guidelines of application of EC techniques for future research

Sticking of H and D atoms on interstellar dust grains is the first step in molecular hydrogen formation, which is a key reaction in the InterStellar Medium (ISM). After studying the sticking coefficients of H2 and D2 molecules on amorphous silicate surfaces experimentally and theoretically, we extrapolate the results to the sticking coefficient of atoms and propose a formulae that gives the sticking coefficients of H and D on both silicates and icy dust grains. In our experiments, we used the King and Wells method for measuring the sticking coefficients of H2 and D2 molecules on a silicate surface held at 10 K. It consists of measuring with a QMS (quadrupole mass spectrometer) the signals of H2 and D2 molecules reflected by the surface during the exposure of the sample to the molecular beam at a temperature ranging from 20 K to 340 K. We tested the efficiency of a physical model, developed previously for sticking on water-ice surfaces. We applied this model to our experimental results for the sticking coefficients of H2 and D2 molecules on a silicate surface and estimated the sticking coefficient of atoms by a single measurement of atomic recombination and propose an extrapolation. Sticking of H, D, HD, H2, and D2 on silicates grains behaves the same as on icy dust grains. The sticking decreases with the gas temperature, and is dependent on the mass of the impactor. The sticking coefficient for both surfaces and impactors can be modeled by an analytical formulae S(T), which describes both the experiments and the thermal distribution expected in an astrophysical context. The parameters S0 and T0 are summarized in a table. Previous estimates for the sticking coefficient of H atoms are close to the new estimation; however, we find that, when isotopic effects are taken into account, the sticking coefficient variations can be as much as a factor of 2 at T=100 K.

Made-to-measure methods such as the parallel code NMAGIC are powerful tools to build galaxy models reproducing observational data. They work by adapting the particle weights in an N-body system until the target observables are well matched. Here we introduce a moving prior regularization (MPR) method for such particle models. It is based on determining from the particles a distribution of priors in phase-space, which are updated in parallel with the weight adaptation. This method allows one to construct smooth models from noisy data without erasing global phase-space gradients. We first apply MPR to a spherical system for which the distribution function can in theory be uniquely recovered from idealized data. We show that NMAGIC with MPR indeed converges to the true solution with very good accuracy, independent of the initial particle model. Compared to the standard weight entropy regularization, biases in the anisotropy structure are removed and local fluctuations in the intrinsic distribution function are reduced. We then investigate how the uncertainties in the inferred dynamical structure increase with less complete and noisier kinematic data, and how the dependence on the initial particle model also increases. Finally, we apply the MPR technique to the two intermediate-luminosity elliptical galaxies NGC 4697 and NGC 3379, obtaining smoother dynamical models in luminous and dark matter potentials.

Recent developments in instrumentation (e.g., in particular the Kepler & CoRoT satellites) provide a new opportunity to improve the models of stellar pulsations. Surface layers, rotation, and magnetic fields imprint erratic frequency shifts, trends, and other non-random behaviour in the frequency spectra. As our observational uncertainties become smaller, these are increasingly important and difficult to deal with using standard fitting techniques. To improve the models, new ways to compare their predictions with observations need to be conceived. In this paper we present a completely probabilistic (Bayesian) approach to asteroseismic model fitting. It allows for varying degrees of prior mode identification, corrections for the discrete nature of the grid, and most importantly implements a treatment of systematic errors, such as the “surface effects.” It removes the need to apply semi- empirical corrections to the observations prior to fitting them to the models and results in a consistent set of probabilities with which the model physics can be probed and compared. As an example, we show a detailed asteroseismic analysis of the Sun. We find a most probable solar age, including a 35 +- 5 million year pre-main sequence phase, of 4.591 billion years, and initial element mass fractions of X_0 = 0.72, Y_0 = 0.264, Z_0 = 0.016, consistent with recent asteroseismic and non-asteroseismic studies.

This paper describes an applicative software tool, named IONORT (IONOspheric Ray Tracing), for calculating a three-dimensional ray tracing of high frequency waves in the ionospheric medium. This tool runs under Windows operating systems and its friendly graphical user interface facilitates both the numerical data input/output and the two/three-dimensional visualization of the ray path. In order to calculate the coordinates of the ray and the three components of the wave vector along the path as dependent variables, the core of the program solves a system of six first order differential equations, the group path being the independent variable of integration. IONORT uses a three-dimensional electron density specification of the ionosphere, as well as by geomagnetic field and neutral particles-electrons collision frequency models having validity in the area of interest.

The measurement of the positions, distances, motions and luminosities of stars represents the foundations of modern astronomical knowledge. Launched at the end of the eighties, the ESA Hipparcos satellite was the first space mission dedicated to such measurements. Hipparcos improved position accuracies by a factor of 100 compared to typical ground-based results and provided astrometric and photometric multi-epoch observations of 118,000 stars over the entire sky. The impact of Hipparcos on astrophysics has been extremely valuable and diverse. Building on this important European success, the ESA Gaia cornerstone mission promises an even more impressive advance. Compared to Hipparcos, it will bring a gain of a factor 50 to 100 in position accuracy and of a factor of 10,000 in star number, collecting photometric, spectrophotometric and spectroscopic data for one billion celestial objects. During its 5-year flight, Gaia will measure objects repeatedly, up to a few hundred times, providing an unprecedented database to study the variability of all types of celestial objects. Gaia will bring outstanding contributions, directly or indirectly, to most fields of research in astrophysics, such as the study of our Galaxy and of its stellar constituents, the search for planets outside the solar system.

This paper presents a hands-on introduction to the medieval astrolabe, based around a working model which can be constructed from photocopies of the supplied figures. As well as describing how to assemble the model, I also provide a brief explanation of how each of its various parts might be used. The printed version of this paper includes only the parts needed to build a single model prepared for use at latitudes around 52{\deg}N, but an accompanying electronic file archive includes equivalent images which can be used to build models prepared for use at any other latitude. The vector graphics scripts used to generate the models are also available for download, allowing customised astrolabes to be made.

In this paper we investigate the performance of the likelihood ratio method as a tool for identifying optical and infrared counterparts to proposed radio continuum surveys with SKA precursor and pathfinder telescopes. We present a comparison of the infrared counterparts identified by the likelihood ratio in the VISTA Deep Extragalactic Observations (VIDEO) survey to radio observations with 6, 10 and 15 arcsec resolution. We cross-match a deep radio catalogue consisting of radio sources with peak flux density $>$ 60 $\mu$Jy with deep near-infrared data limited to $K_{\mathrm{s}}\lesssim$ 22.6. Comparing the infrared counterparts from this procedure to those obtained when cross-matching a set of simulated lower resolution radio catalogues indicates that degrading the resolution from 6 arcsec to 10 and 15 arcsec decreases the completeness of the cross-matched catalogue by approximately 3 and 7 percent respectively. When matching against shallower infrared data, comparable to that achieved by the VISTA Hemisphere Survey, the fraction of radio sources with reliably identified counterparts drops from $\sim$89%, at $K_{\mathrm{s}}\lesssim$22.6, to 47% with $K_{\mathrm{s}}\lesssim$20.0. Decreasing the resolution at this shallower infrared limit does not result in any further decrease in the completeness produced by the likelihood ratio matching procedure. However, we note that radio continuum surveys with the MeerKAT and eventually the SKA, will require long baselines in order to ensure that the resulting maps are not limited by instrumental confusion noise.

We present a new algorithm called ‘Fast Integrated Spectra Analyzer” (FISA) that permits fast and reasonably accurate age and reddening determinations for small angular diameter open clusters by using their integrated spectra in the (3600-7400) \AA \ range and currently available template spectrum libraries. This algorithm and its implementation help to achieve astrophysical results in shorter times than from other methods. A brief review is given of the integrated spectroscopic technique applied to the study of open clusters as well as the basic assumptions that justify its use. We describe the numerical algorithm employed in detail, show examples of its application, and provide a link to the code. Our method has successfully been applied to integrated spectroscopy of open clusters, both in the Galaxy and in the Magellanic Clouds, to determine ages and reddenings.

We performed photometric calibration of the PhotoMultiplier Tube (PMT) and readout electronics used for the new fluorescence detectors of the Telescope Array (TA) experiment using Rayleigh scattered photons from a pulsed nitrogen laser beam. The experimental setup, measurement procedure, and results of calibration are described. The total systematic uncertainty of the calibration is estimated to be 7.5%. An additional uncertainty of 3.7% is introduced by the transport of the calibrated PMTs from the laboratory to the TA experimental site.

A new data analysis method based on physical observables for WIMP dark matter searches with noble liquid Xe dual-phase TPCs is presented. Traditionally, the nuclear recoil energy from a scatter in the liquid target has been estimated by means of the initial prompt scintillation light (S1) produced at the interaction vertex. The ionization charge (C2), or its secondary scintillation (S2), is combined with the primary scintillation in Log(S2/S1) vs. S1 only as a discrimination parameter against electron recoil background. Arguments in favor of C2 as the more reliable nuclear recoil energy estimator than S1 are presented. The new phase space of Log(S1/C2) vs. C2 is introduced as more efficient for nuclear recoil acceptance and exhibiting superior energy resolution. This is achieved without compromising the discrimination power of the LXe TPC, nor its 3D event reconstruction and fiducialization capability, as is the case for analyses that exploit only the ionization channel. Finally, the concept of two independent energy estimators for background rejection is presented: E2 as the primary (based on C2) and E1 as the secondary (based on S1). Log(E1/E2) vs. E2 is shown to be the most appropriate phase space in which to evaluate WIMP signal candidates.

A new data analysis method based on physical observables for WIMP dark matter searches with noble liquid Xe dual-phase TPCs is presented. Traditionally, the nuclear recoil energy from a scatter in the liquid target has been estimated by means of the initial prompt scintillation light (S1) produced at the interaction vertex. The ionization charge (C2), or its secondary scintillation (S2), is combined with the primary scintillation in Log(S2/S1) vs. S1 only as a discrimination parameter against electron recoil background. Arguments in favor of C2 as the more reliable nuclear recoil energy estimator than S1 are presented. The new phase space of Log(S1/C2) vs. C2 is introduced as more e?cient for nuclear recoil acceptance and exhibiting superior energy resolution. This is achieved without compromising the discrimination power of the LXe TPC, nor its 3D event reconstruction and ?ducialization capability, as is the case for analyses that exploit only the ionization channel. Finally, the concept of two independent energy estimators for background rejection is presented: E2 as the primary (based on C2) and E1 as the secondary (based on S1). Log(E1/E2) vs. E2 is shown to be the most appropriate phase space in which to evaluate WIMP signal candidates.

We present a laser frequency comb based upon a 250 MHz mode-locked erbium-doped fiber laser that spans more than 300 terahertz of bandwidth, from 660 nm to 2000 nm. The system generates 1.2 nJ, 70 fs pulses at 1050 nm by amplifying the 1580 nm laser light in Er:fiber, followed by nonlinear broadening to 1050 nm and amplification in Yb:fiber. Extension of the frequency comb into the visible is achieved by supercontinuum generation from the 1050 nm light. Comb coherence is verified with cascaded f-2f interferometry and comparison to a frequency stabilized laser.

The Kilodegree Extremely Little Telescope (KELT) project is a survey for new transiting planets around bright stars. KELT-South is a small-aperture, wide-field automated telescope located at Sutherland, South Africa. The telescope surveys a set of 26 degree by 26 degree fields around the southern sky, and targets stars in the range of 8 < V < 10 mag, searching for transits by Hot Jupiters. This paper describes the KELT-South system hardware and software and discusses the quality of the observations. We show that KELT-South is able to achieve the necessary photometric precision to detect transits of Hot Jupiters around solar-type main-sequence stars.

X-ray Multi-Mirror Mission (XMM-Newton) has been one of the most successful astronomy missions launched by the European Space Agency. It exploits innovative use of replication technology for the X-ray reflecting telescopes that has resulted in an unprecedented combination of effective area and resolution. Three telescopes are equipped with imaging cameras and spectrometers that operate simultaneously, together with a coaligned optical telescope. The key features of the payload are described, and the in-orbit performance and scientific achievements are summarised. Subject terms or keywords: XMM-Newton, X-ray astronomy, space telescopes

Genetic algorithms (GAs) emulate the process of biological evolution, in a computational setting, in order to generate good solutions to difficult search and optimisation problems. GA-based optimisers tend to be extremely robust and versatile compared to most traditional techniques used to solve optimisation problems. This review paper provides a very brief introduction to GAs and outlines their utility in astronomy and astrophysics.

Micro-electro-mechanical systems (MEMS) technology can provide for deformable mirrors (DMs) with excellent performance within a favorable economy of scale. Large MEMS-based astronomical adaptive optics (AO) systems such as the Gemini Planet Imager are coming on-line soon. As MEMS DM end-users, we discuss our decade of practice with the micromirrors, from inspecting and characterizing devices to evaluating their performance in the lab. We also show MEMS wavefront correction on-sky with the “Villages” AO system on a 1-m telescope, including open-loop control and visible-light imaging. Our work demonstrates the maturity of MEMS technology for astronomical adaptive optics.

Double-peaked line profiles are commonly considered a hallmark of rotating disks, with the distance between the peaks a measure of the rotation velocity. However, double-peaks can arise also from radiative transfer effects in optically thick non-rotating sources. Utilizing exact solutions of the line transfer problem we present a detailed study of line emission from geometrically thin Keplerian disks. We derive the conditions for emergence of kinematic double peaks in optically thin and thick disks, and find that it is generally impossible to disentangle the effects of kinematics and line opacity in observed double-peaked profiles. Unless supplemented by additional information, a double-peaked profile alone is not a reliable indicator of a rotating disk. In certain circumstances, triple and quadruple profiles might be better indicators of rotation in optically thick disks.

“Approximate Bayesian Computation” (ABC) represents a powerful methodology for the analysis of complex stochastic systems for which the likelihood of the observed data under an arbitrary set of input parameters may be entirely intractable-the latter condition rendering useless the standard machinery of tractable likelihood-based, Bayesian statistical inference (e.g. conventional Markov Chain Monte Carlo simulation; MCMC). In this article we demonstrate the potential of ABC for astronomical model analysis by application to a case study in the morphological transformation of high redshift galaxies. To this end we develop, first, a stochastic model for the competing processes of merging and secular evolution in the early Universe; and second, through an ABC-based comparison against the observed demographics of the first generation of massive (M_gal > 10^11 M_sun) galaxies (at 1.5 < z < 3) in the CANDELS/EGS dataset we derive posterior probability densities for the key parameters of this model. The "Sequential Monte Carlo" (SMC) implementation of ABC exhibited herein, featuring both a self-generating target sequence and self-refining MCMC kernel, is amongst the most efficient of contemporary approaches to this important statistical algorithm. We highlight as well through our chosen case study the value of careful summary statistic selection, and demonstrate two modern strategies for assessment and optimisation in this regard. Ultimately, our ABC analysis of the high redshift morphological mix returns tight constraints on the evolving merger rate in the early Universe and reveals merging, rather than secular evolution, as the most important mechanism for building up the first generation of bulges in early-type disks.

We study the asymptotic behaviour of needlets-based approximate maximum likelihood estimators for the spectral parameters of Gaussian and isotropic spherical random fields. We prove consistency and asymptotic Gaussianity, in the high-frequency limit, thus generalizing earlier results by Durastanti et al. (2011) based upon standard Fourier analysis on the sphere. The asymptotic results are then illustrated by an extensive Monte Carlo study.

Post-launch boresight of the IBEX-Lo instrument onboard the Interstellar Boundary Explorer (IBEX) is determined based on IBEX-Lo Star Sensor observations. Accurate information on the boresight of the neutral gas camera is essential for precise determination of interstellar gas flow parameters. Utilizing spin-phase information from the spacecraft attitude control system (ACS), positions of stars observed by the Star Sensor during two years of IBEX measurements were analyzed and compared with positions obtained from a star catalog. No statistically significant differences were observed beyond those expected from the pre-launch uncertainty in the Star Sensor mounting. Based on the star observations and their positions in the spacecraft reference system, pointing of the IBEX satellite spin axis was determined and compared with the pointing obtained from the ACS. Again, no statistically significant deviations were observed. We conclude that no systematic correction for boresight geometry is needed in the analysis of IBEX-Lo observations to determine neutral interstellar gas flow properties. A stack-up of uncertainties in attitude knowledge shows that the instantaneous IBEX-Lo pointing is determined to within $\sim 0.1\degr$ in both spin angle and elevation using either the Star Sensor or the ACS. Further, the Star Sensor can be used to independently determine the spacecraft spin axis. Thus, Star Sensor data can be used reliably to correct the spin phase when the Star Tracker (used by the ACS) is disabled by bright objects in its field-of-view. The Star Sensor can also determine the spin axis during most orbits and thus provides redundancy for the Star Tracker.

This analysis begins to explore the complex chromosphere-corona mass cycle using a blend of imaging and spectroscopic diagnostics. Single Gaussian fits to hot emission line profiles (formed above 1MK) at the base of coronal loop structures indicate material blue-shifts of 5-10km/s while cool emission line profiles (formed below 1MK) yield red-shifts of a similar magnitude – indicating, to zeroth order, that a temperature-dependent bifurcating flow exists on coronal structures. Image sequences of the same region reveal weakly emitting upward propagating disturbances in both hot and cool emission with apparent speeds of 50-150km/s. Spectroscopic observations indicate that these propagating disturbances produce a weak emission component in the blue wing at commensurate speed, but that they contribute only a few percent to the (ensemble) emission line profile in a single spatio-temporal resolution element. Subsequent analysis of imaging data shows material “draining” slowly (~10km/s) out of the corona, but only in the cooler passbands. We interpret the draining as the return-flow of coronal material at the end of the complex chromosphere-corona mass cycle. Further, we suggest that the efficient radiative cooling of the draining material produces a significant contribution to the red wing of cool emission lines that is ultimately responsible for their systematic red-shift as derived from a single Gaussian fit when compared to those formed in hotter (conductively dominated) domains. The presence of counter-streaming flows complicates the line profiles, their interpretation, and asymmetry diagnoses, but allows a different physical picture of the lower corona to develop.

The Clementine gnomon has been built in 1702 to measure the Earth’s obliquity variation. For this reason the pinhole was located in the walls of Diocletian’s times (305 a. D.) in order to remain stable along the centuries, but its original form and position have been modified. We used an astrometric method to recover the original position of the pinhole: reshaping the pinhole to a circle of 1.5 cm of diameter, the positions of the Northern and Southern limbs have been compared with the ephemerides. A sistematic shift of 4.5 mm Southward of the whole solar image shows that the original pinhole was 4.5 mm North of the actual position, as the images in the Bianchini’s book (1703) suggest. The oval shape of the actual pinhole is also wrong. Using a circle the larger solar spots are clearly visible. Some reference stars of the catalogue of Philippe de la Hire (1702), used originally for measuring the ecliptic latitude of the Sun, are written next to the meridian line, but after the last restauration (2000), four of them are wrongly located. Finally the deviation from the true North, of the meridian line’s azimuth confirms the value recovered in 1750. This, with the local deviations of a true line, will remain as systematic error, like for all these historical instruments.

We develop a novel method to separate the components of a diffuse emission process based on an association with the energy spectra. Most of the existing methods use some information about the spatial distribution of components, e.g., closeness to an external template, independence of components etc., in order to separate them. In this paper we propose a method where one puts conditions on the spectra only. The advantages of our method are: 1) it is internal: the maps of the components are constructed as combinations of data in different energy bins, 2) the components may be correlated among each other, 3) the method is semi-blind: in many cases, it is sufficient to assume a functional form of the spectra and determine the parameters from a maximization of a likelihood function. As an example, we derive the CMB map and the foreground maps for seven yeas of WMAP data. In an Appendix, we present a generalization of the method, where one can also add a number of external templates.

Building on the “drPACS” contribution at ADASS XX of a simple Unix pipeline infrastructure, we implemented a pipeline toolkit using the package MIRIAD to combine Interferometric and Single Dish data (MIS). This was prompted by our observations made with the Combined Array For Research in Millimeter-wave Astronomy (CARMA) interferometer of the star-forming region NGC 1333, a large survey highlighting the new 23-element and singledish observing modes. The project consists of 20 CARMA datasets each containing interferometric as well as simultaneously obtained single dish data, for 3 molecular spectral lines and continuum, in 527 different pointings, covering an area of about 8 by 11 arcminutes. A small group of collaborators then shared this toolkit and their parameters via CVS, and scripts were developed to ensure uniform data reduction across the group. The pipeline was run end-to-end each night as new observations were obtained, producing maps that contained all the data to date. We will show examples of the scripts and data products. This approach could serve as a model for repeated calibration and mapping of large mixed-mode correlation datasets from ALMA.

Here we describe the Astrophysics Source Code Library (ASCL), which takes an active approach to sharing astrophysical source code. ASCL’s editor seeks out both new and old peer-reviewed papers that describe methods or experiments that involve the development or use of source code, and adds entries for the found codes to the library. This approach ensures that source codes are added without requiring authors to actively submit them, resulting in a comprehensive listing that covers a significant number of the astrophysics source codes used in peer-reviewed studies. The ASCL now has over 340 codes in it and continues to grow. In 2011, the ASCL (http://ascl.net) has on average added 19 new codes per month. An advisory committee has been established to provide input and guide the development and expansion of the new site, and a marketing plan has been developed and is being executed. All ASCL source codes have been used to generate results published in or submitted to a refereed journal and are freely available either via a download site or from an identified source. This paper provides the history and description of the ASCL. It lists the requirements for including codes, examines the benefits of the ASCL, and outlines some of its future plans.

Much of scientific progress now hinges on the reliability, falsifiability and reproducibility of computer source codes. Astrophysics in particular is a discipline that today leads other sciences in making useful scientific components freely available online, including data, abstracts, preprints, and fully published papers, yet even today many astrophysics source codes remain hidden from public view. We review the importance and history of source codes in astrophysics and previous efforts to develop ways in which information about astrophysics codes can be shared. We also discuss why some scientist coders resist sharing or publishing their codes, the reasons for and importance of overcoming this resistance, and alert the community to a reworking of one of the first attempts for sharing codes, the Astrophysics Source Code Library (ASCL). We discuss the implementation of the ASCL in an accompanying poster paper. We suggest that code could be given a similar level of referencing as data gets in repositories such as ADS.

A new spectrometer for the Green Bank Telescope (GBT) is being built jointly by the NRAO and the CASPER, University of California, Berkeley. The spectrometer uses 8 bit ADCs and will be capable of processing up to 1.25 GHz bandwidth from 8 dual polarized beams. This mode will be used to process data from focal plane arrays. The spectrometer supports observing mode with 8 tunable digital sub-bands within the 1.25 GHz bandwidth. The spectrometer can also be configured to process a bandwidth of up to 10 GHz with 64 tunable sub-bands from a dual polarized beam. The vastly enhanced backend capabilities will support several new science projects with the GBT.

Precise astronomical spectroscopic analyses routinely assume that individual pixels in charge-coupled devices (CCDs) have uniform sensitivity to photons. Intra-pixel sensitivity (IPS) variations may already cause small systematic errors in, for example, studies of extra-solar planets via stellar radial velocities and cosmological variability in fundamental constants via quasar spectroscopy, but future experiments requiring velocity precisions approaching ~1 cm/s will be more strongly affected. Laser frequency combs have been shown to provide highly precise wavelength calibration for astronomical spectrographs, but here we show that they can also be used to measure IPS variations in astronomical CCDs in situ. We successfully tested a laser frequency comb system on the Ultra-High Resolution Facility spectrograph at the Anglo-Australian Telescope. By modelling the 2-dimensional comb signal recorded in a single CCD exposure, we find that the average IPS deviates by <8 per cent if it is assumed to vary symmetrically about the pixel centre. We also demonstrate that series of comb exposures with absolutely known offsets between them can yield tighter constraints on symmetric IPS variations from ~100 pixels. We discuss measurement of asymmetric IPS variations and absolute wavelength calibration of astronomical spectrographs and CCDs using frequency combs.

Gravitational Waves (GWs) are tiny ripples in the fabric of space-time predicted by Einstein’s General Relativity. Pulsar timing arrays (PTAs) are well poised to detect low frequency ($10^{-9}$ — $10^{-7}$ Hz) GWs in the near future. There has been a significant amount of research into the detection of a stochastic background of GWs from supermassive black hole binaries (SMBHBs). Recent work has shown that single continuous sources standing out above the background may be detectable by PTAs operating at a sensitivity sufficient to detect the stochastic background. The most likely sources of continuous GWs in the pulsar timing frequency band are extremely massive and/or nearby SMBHBs. In this paper we present detection strategies including various forms of matched filtering and power spectral summing. We determine the efficacy and computational cost of such strategies. It is shown that it is computationally infeasible to use an optimal matched filter including the poorly constrained pulsar distances with a grid based method. We show that an Earth-term-matched filter constructed using only the correlated signal terms is both computationally viable and highly sensitive to GW signals. This technique is only a factor of two less sensitive than the computationally unrealizable optimal matched filter and a factor of two more sensitive than a power spectral summing technique. We further show that a pairwise matched filter, taking the pulsar distances into account is comparable to the optimal matched filter for the single template case and comparable to the Earth-term-matched filter for many search templates. Finally, using simulated data optimal quality, we place a theoretical minimum detectable strain amplitude of $h>2\times 10^{-15}$ from continuous GWs at frequencies on the order $\sim1/T_{\rm obs}$.

Stellar feedback drives the circulation of matter from the disk to the halo of galaxies. We perform three-dimensional magnetohydrodynamic simulations of a vertical column of the interstellar medium with initial conditions typical of the solar circle in which supernovae drive turbulence and determine the vertical stratification of the medium. The simulations were run using a stable, positivity-preserving scheme for ideal MHD implemented in the FLASH code. We find that the majority (\approx 90 %) of the mass is contained in thermally-stable temperature regimes of cold molecular and atomic gas at T < 200 K or warm atomic and ionized gas at 5000 K < T < 10^{4.2} K, with strong peaks in probability distribution functions of temperature in both the cold and warm regimes. The 200 – 10^{4.2} K gas fills 50-60 % of the volume near the plane, with hotter gas associated with supernova remnants (30-40 %) and cold clouds ( 3 kpc. The magnetic field in our models has no significant impact on the scale heights of gas in each temperature regime; the magnetic tension force is approximately equal to and opposite the magnetic pressure, so the addition of the field does not significantly affect the vertical support of the gas. The addition of a magnetic field does reduce the fraction of gas in the cold (< 200 K) regime with a corresponding increase in the fraction of warm (~ 10^4 K) gas. However, our models lack rotational shear and thus have no large-scale dynamo, which reduces the role of the field in the models compared to reality. The supernovae drive oscillations in the vertical distribution of halo gas, with the period of the oscillations ranging from ~ 30 Myr in the T < 200 K gas to ~ 100 Myr in the 10^6 K gas, in line with predictions by Walters & Cox.

Stellar feedback drives the circulation of matter from the disk to the halo of galaxies. We perform three-dimensional magnetohydrodynamic simulations of a vertical column of the interstellar medium with initial conditions typical of the solar circle in which supernovae drive turbulence and determine the vertical stratification of the medium. The simulations were run using a stable, positivity-preserving scheme for ideal MHD implemented in the FLASH code. We find that the majority (\approx 90 %) of the mass is contained in thermally-stable temperature regimes of cold molecular and atomic gas at T < 200 K or warm atomic and ionized gas at 5000 K < T < 10^{4.2} K, with strong peaks in probability distribution functions of temperature in both the cold and warm regimes. The 200 – 10^{4.2} K gas fills 50-60 % of the volume near the plane, with hotter gas associated with supernova remnants (30-40 %) and cold clouds ( 3 kpc. The magnetic field in our models has no significant impact on the scale heights of gas in each temperature regime; the magnetic tension force is approximately equal to and opposite the magnetic pressure, so the addition of the field does not significantly affect the vertical support of the gas. The addition of a magnetic field does reduce the fraction of gas in the cold (< 200 K) regime with a corresponding increase in the fraction of warm (~ 10^4 K) gas. However, our models lack rotational shear and thus have no large-scale dynamo, which reduces the role of the field in the models compared to reality. The supernovae drive oscillations in the vertical distribution of halo gas, with the period of the oscillations ranging from ~ 30 Myr in the T < 200 K gas to ~ 100 Myr in the 10^6 K gas, in line with predictions by Walters & Cox.

A probabilistic two-component mixture model allows one to separate the diffuse background from the celestial sources within a one-step algorithm without data censoring. The background is modeled with a thin-plate spline combined with the satellite’s exposure time. Source probability maps are created in a multi-resolution analysis for revealing faint and extended sources. All detected sources are automatically parametrized to produce a list of source positions, fluxes and morphological parameters. The present analysis is applied to the Chandra Deep Field South 2 Ms public released data. Within its 1.884 ks of exposure time and its angular resolution (0.984 arcsec), the Chandra Deep Field South data are particularly suited for testing the Background-Source separation algorithm.

Astronomical images in the Poisson regime are typically characterized by a spatially varying cosmic background, large variety of source morphologies and intensities, data incompleteness, steep gradients in the data, and few photon counts per pixel. The Background-Source separation technique is developed with the aim to detect faint and extended sources in astronomical images characterized by Poisson statistics. The technique employs Bayesian mixture models to reliably detect the background as well as the sources with their respective uncertainties. Background estimation and source detection is achieved in a single algorithm. A large variety of source morphologies is revealed. The technique is applied in the X-ray part of the electromagnetic spectrum on ROSAT and Chandra data sets and it is under a feasibility study for the forthcoming eROSITA mission.

We present a high-precision infrared radial velocity study of late-type stars using spectra obtained with NIRSPEC at the W. M. Keck Observatory. Radial velocity precisions of 50 m/s are achieved for old field mid-M dwarfs using telluric features for precise wavelength calibration. Using this technique, 20 young stars in the {\beta} Pic (age ~12 Myr) and TW Hya (age ~8 Myr) Associations were monitored over several years to search for low mass companions; we also included the chromospherically active field star GJ 873 (EV Lac) in this survey. Based on comparisons with previous optical observations of these young active stars, radial velocity measurements at infrared wavelengths mitigate the radial velocity noise caused by star spots by a factor of ~3. Nevertheless, star spot noise is still the dominant source of measurement error for young stars at 2.3 {\mu}m, and limits the precision to ~77 m/s for the slowest rotating stars (v sin i 12 km/s). The observations reveal both GJ 3305 and TWA 23 to be single-lined spectroscopic binaries; in the case of GJ 3305, the motion is likely caused by its 0.09″ companion, identified after this survey began. The large amplitude, short-timescale variations of TWA 13A are indicative of a hot Jupiter-like companion, but the available data are insufficient to confirm this. We label it as a candidate radial velocity variable. For the remainder of the sample, these observations exclude the presence of any ‘hot’ (P < 3 days) companions more massive than 8 MJup, and any 'warm' (P < 30 days) companions more massive than 17 MJup, on average. Assuming an edge-on orbit for the edge-on disk system AU Mic, these observations exclude the presence of any hot Jupiters more massive than 1.8 MJup or warm Jupiters more massive than 3.9 MJup.

A detailed and formal account of polarization measurements using Bayesian analysis is given based on the assumption of gaussian error for the Stokes parameters. This analysis is crucial for the measurement of the polarization degree and angle at very low (and very high) signal-to-noise. The treatment serves as a framework for customized analysis of data based on a particular prior suited to the experiment.

While fully automated methods for detecting faint moving objects in astronomical images – such as Kuiper belt objects (KBOs) – are constantly improving, visual detection still has a role to play especially when the fixed background is cluttered with stars. Color intensity projections (CIPs) using hue cycling – which combines a sequence of greyscale images into a single color image – aids in the visual detection of moving objects by highlighting them using color in an intuitive way. To demonstrate the usefulness of CIPs in detecting faint moving objects a sequence of 16 images from the SuprimeCam camera of the Subaru telescope were combined into a CIPs image. As well has making even faint moving objects easier to visually detect against a cluttered background, CCD artefacts were also more easily recognisable. The new Hyper SuprimeCam for the Subaru telescope – which will allow many short exposure images to be acquired with little dead time between images – should provide ideal data for use with the CIPs algorithm. In addition, the current search for KBOs to be targeted by the New Horizon’s spacecraft after its flyby of Pluto provides an excellent test case for the state of the art in faint moving object detection against a cluttered background.

The neutrino telescopes of the present generation, depending on their specific features, can reconstruct the neutrino spectra from a galactic burst. Since the optical counterpart could be not available, it is desirable to have at hand alternative methods to estimate the distance of the supernova explosion using only the neutrino data. In this work we present preliminary results on the method we are proposing to estimate the distance from a galactic supernova based only on the spectral shape of the neutrino burst and assumptions on the gravitational binding energy released an a typical supernova explosion due to stellar collapses.

Context. This study has been developed in the framework of the computational simulations executed for the preparation of the ESA Gaia astrometric mission. Aims. We focus on describing the objects and characteristics that Gaia will potentially observe without taking into consideration instrumental effects (detection efficiency, observing errors). Methods. The theoretical Universe Model prepared for the Gaia simulation has been statistically analyzed at a given time. Ingredients of the model are described, giving most attention to the stellar content, the double and multiple stars, and variability. Results. In this simulation the errors have not been included yet. Hence we estimate the number of objects and their theoretical photometric, astrometric and spectroscopic characteristics in the case that they are perfectly detected. We show that Gaia will be able to potentially observe 1.1 billion of stars (single or part of multiple star systems) of which about 2% are variable stars, 3% have one or two exoplanets. At the extragalactic level, observations will be potentially composed by several millions of galaxies, half million to 1 million of quasars and about 50,000 supernovas that will occur during the 5 years of mission. The simulated catalogue will be made publicly available by the DPAC on the Gaia portal of the ESA web site http://www.rssd.esa.int/gaia/.

Context. This study has been developed in the framework of the computational simulations executed for the preparation of the ESA Gaia astrometric mission. Aims. We focus on describing the objects and characteristics that Gaia will potentially observe without taking into consideration instrumental effects (detection efficiency, observing errors). Methods. The theoretical Universe Model prepared for the Gaia simulation has been statistically analyzed at a given time. Ingredients of the model are described, giving most attention to the stellar content, the double and multiple stars, and variability. Results. In this simulation the errors have not been included yet. Hence we estimate the number of objects and their theoretical photometric, astrometric and spectroscopic characteristics in the case that they are perfectly detected. We show that Gaia will be able to potentially observe 1.1 billion of stars (single or part of multiple star systems) of which about 2% are variable stars, 3% have one or two exoplanets. At the extragalactic level, observations will be potentially composed by several millions of galaxies, half million to 1 million of quasars and about 50,000 supernovas that will occur during the 5 years of mission. The simulated catalogue will be made publicly available by the DPAC on the Gaia portal of the ESA web site http://www.rssd.esa.int/gaia/.

The accuracy in the photometry of a point source depends on the point-spread function (PSF), detector pixelization, and observing strategy. The PSF and pixel response describe the spatial blurring of the source, the pixel scale describes the spatial sampling of a single exposure, and the observing strategy determines the set of dithered exposures with pointing offsets from which the source flux is inferred. In a wide-field imaging survey, sources of interest are randomly distributed within the field of view and hence are centered randomly within a pixel. A given hardware configuration and observing strategy therefore have a distribution of photometric uncertainty for sources of fixed flux that fall in the field. In this article we explore the ensemble behavior of photometric and position accuracies for different PSFs, pixel scales, and dithering patterns. We find that the average uncertainty in the flux determination depends slightly on dither strategy, whereas the position determination can be strongly dependent on the dithering. For cases with pixels much larger than the PSF, the uncertainty distributions can be non-Gaussian, with rms values that are particularly sensitive to the dither strategy. We also find that for these configurations with large pixels, pointings dithered by a fractional pixel amount do not always give minimal average uncertainties; this is in contrast to image reconstruction for which fractional dithers are optimal. When fractional pixel dithering is favored, a pointing accuracy of better than $\sim 0.15$ pixel width is required to maintain half the advantage over random dithers.

The late-type dwarf GJ 436 is known to host a transiting Neptune-mass planet in a 2.6-day orbit. We present results of our interferometric measurements to directly determine the stellar diameter ($R_{\star} = 0.455 \pm 0.018 R_{\odot}$) and effective temperature ($T_{\rm EFF} = 3416 \pm 54$ K). We combine our stellar parameters with literature time-series data, which allows us to calculate physical and orbital system parameters, including GJ 436’s stellar mass ($M_{\star} = 0.472^{+ 0.0636}_{- 0.0566} M_{\odot}$), planetary radius ($R_{p} = 0.370^{+ 0.0149}_{- 0.0145} R_{Jupiter}$), planetary mass ($M_{p} = 0.075^{+ 0.0076}_{- 0.0072} M_{Jupiter}$), implying a mean planetary density of $\rho_{p} = 1.48^{+ 0.116}_{- 0.103} \rho_{Jupiter}$. These values are generally in good agreement with previous literature estimates based on assumed stellar mass and photometric light curve fitting. Finally, we examine the expected phase curves of the hot Neptune GJ 436b, based on various assumptions concerning the efficiency of energy redistribution in the planetary atmosphere, and find that it could be constrained with {\it Spitzer} monitoring observations.

The problem of photometric identification will be extremely important for large surveys in the next decade. In this work, we propose the use of KPCA combined with k = 1 nearest neighbor algorithm (KPCA+1NN) as a framework for SNe photometric classification. The method does not rely on information about redshift or local enviromental variables, so it is less sensitive to bias than its template fitting counterparts. We applied the method to $\approx$ 20000 SNe light curve released after the \textit{Supernova Photometric Classification Challenge} (SNPCC). Results for the photometric sample achieved up to 89% efficiency (eff), 97% purity (pur), 96% successful classification (SC) rates and figure of merit (FoM) of 0.79 (SNR$\geq$5). If we impose no SNR cuts, we obtain up to 64% eff, 43% pur, 46% SC and FoM of 0.10. We also present the classification results using only pre-maximum epoches, obtaining 80% eff, 73% pur, 84% SC and FoM of 0.32 (SNR$\geq$5). Comparing the performance of our classifier with MLCS2k2 fit probability, we demonstrate that KPCA+1NN is able to improve the classification results of MLCS2k up to $\approx 15%$ without the need of redshift information. Results are sensitive to the information contained in each light curve, as a consequence, higher quality the data points lead to higher successfull classification rates. The method is flexible enough to be applied to other astrophysical transients, as long as a training and a template sample are provided.

We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7×10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.

The radiation pressure due to reflected chameleons is suggested for their detection. Sensitive photon radiation pressure or force measuring techniques could become the new antennas for chameleons or other particles with similar properties, with gravitational waves antennas being also of potential interest. This is eventually only at first sight a Gedankenexperiment. Also chameleon-through-a-wall experiments, in particular when equipped with a Fabry-Perot resonator, could apply the suggested techniques.

The action of the atmospheric seeing is blurring, image stretching and image motion. This happens even to the image of the Sun which is more than half degree wide. Low frequency seeing components affect the solar diameter values measured either throught the drift-scan or the heliometer methods. We present evidences of image motion and stretching down to 0.001 Hz.

The action of the atmospheric seeing is blurring, image stretching and image motion. This happens even to the image of the Sun which is more than half degree wide. Low frequency seeing components affect the solar diameter values measured either throught the drift-scan or the heliometer methods. We present evidences of image motion and stretching down to 0.01 Hz.