We present high-resolution I-band imaging of the core of the globular cluster M15 obtained at the 2.5 m Nordic Optical Telescope with FastCam, a low readout noise L3CCD based instrument. Short exposure times (30 ms) were used to record 200000 images (512 x 512 pixels each) over a period of 2 hours 43 min. The lucky imaging technique was then applied to generate a final image of the cluster centre with FWHM ~ 0″.1 and 13″ x 13″ FoV. We obtained a catalogue of objects in this region with a limiting magnitude of I=19.5. I-band photometry and astrometry are reported for 1181 stars. This is the deepest I-band observation of the M15 core at this spatial resolution. Simulations show that crowding is limiting the completeness of the catalogue. At shorter wavelengths, a similar number of objects has been reported using HST/WFPC observations of the same field. The cross-match with the available HST catalogues allowed us to produce colour-magnitude diagrams where we identify new Blue Straggler star candidates and previously known stars of this class.

This document describes the exposure time calculator for the Wide-Field Infrared Survey Telescope (WFIRST) high-latitude survey. The calculator works in both imaging and spectroscopic modes. In addition to the standard ETC functions (e.g. background and S/N determination), the calculator integrates over the galaxy population and forecasts the density and redshift distribution of galaxy shapes usable for weak lensing (in imaging mode) and the detected emission lines (in spectroscopic mode). The source code is made available for public use.

Weak-lensing shear estimates show a troublesome dependence on the apparent brightness of the galaxies used to measure the ellipticity: In several studies, the amplitude of the inferred shear falls sharply with decreasing source significance. This dependence limits the overall ability of upcoming large weak-lensing surveys to constrain cosmological parameters. We seek to provide a concise overview of the impact of pixel noise on weak-lensing measurements, covering the entire path from noisy images to shear estimates. We show that there are at least three distinct layers, where pixel noise not only obscures but biases the outcome of the measurements: 1) the propagation of pixel noise to the non-linear observable ellipticity; 2) the response of the shape-measurement methods to limited amount of information extractable from noisy images; and 3) the reaction of shear estimation statistics to the presence of noise and outliers in the measured ellipticities. We identify and discuss several fundamental problems and show that each of them is able to introduce biases in the range of a few tenths to a few percent for galaxies with typical significance levels. Furthermore, all of these biases do not only depend on the brightness of galaxies but also on their ellipticity, with more elliptical galaxies often being harder to measure correctly. We also discuss existing possibilities to mitigate and novel ideas to avoid the biases induced by pixel noise. We present a new shear estimator that shows a more robust performance for noisy ellipticity samples. Finally, we release the open-source python code to predict and efficiently sample from the noisy ellipticity distribution and the shear estimators used in this work at https://github.com/pmelchior/epsnoise.

We describe an FPGA based Front-End Data Acquisition Module (FEDAM) for implementing Time-over-Threshold (ToT) Time-to-Digital conversion (TDC) of pulses obtained from the COSMICi astroparticle telescope detector system photomultiplier tubes. The telescope system consists of a minimum of three scintillation detectors configured to detect particle airshowers likely initiated by Ultra High Energy Cosmic Rays (UHECR). The relative time delay of detection events between the detectors is used to estimate the angle of incidence of the shower. The FEDAM provides time-over-threshold measurements with a resolution of 2 ns. This allows determination of shower direction to an error of 0.035 (cos {\theta})-1 radians where {\theta} is the angle between the baseline axis through a pair of detectors and the plane representing the shower front.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

The ANTARES detector is the largest deep sea underwater neutrino telescope in operation. The apparatus comprises a matrix of 885 photomultiplier tubes (PMTs) which detect the Cherenkov light emitted by the charged leptons produced in the charged current interactions of high energy neutrinos with the matter inside or near the detector. Reconstruction of the muon track and energy can be achieved using the time, position and charge information of the hits arriving to the PMTs. A good calibration of the detector is necessary in order to ensure its optimal performance. This contribution reviews the different calibration systems and methods developed by the ANTARES Collaboration.

LWA1 is a new radio telescope operating in the frequency range 10-88 MHz, located in central New Mexico. The telescope consists of 258 pairs of dipole-type antennas whose outputs are individually digitized and formed into beams. Simultaneously, signals from all dipoles can be recorded using one of the instrument’s “all dipoles” modes, facilitating all-sky imaging. Notable features of the instrument include high intrinsic sensitivity (about 6 kJy zenith system equivalent flux density), large instantaneous bandwidth (up to 78 MHz), and 4 independently-steerable beams utilizing digital “true time delay” beamforming. This paper summarizes the design of LWA1 and its performance as determined in commissioning experiments. We describe the method currently in use for array calibration, and report on measurements of sensitivity and beamwidth.

LWA1 is a new radio telescope operating in the frequency range 10-88 MHz, located in central New Mexico. The telescope consists of 258 pairs of dipole-type antennas whose outputs are individually digitized and formed into beams. Simultaneously, signals from all dipoles can be recorded using one of the instrument’s “all dipoles” modes, facilitating all-sky imaging. Notable features of the instrument include high intrinsic sensitivity (about 6 kJy zenith system equivalent flux density), large instantaneous bandwidth (up to 78 MHz), and 4 independently-steerable beams utilizing digital “true time delay” beamforming. This paper summarizes the design of LWA1 and its performance as determined in commissioning experiments. We describe the method currently in use for array calibration, and report on measurements of sensitivity and beamwidth.

We propose Paraiso, a domain specific language embedded in functional programming language Haskell, for automated tuning of explicit solvers of partial differential equations (PDEs) on GPUs as well as multicore CPUs. In Paraiso, one can describe PDE solving algorithms succinctly using tensor equations notation. Hydrodynamic properties, interpolation methods and other building blocks are described in abstract, modular, re-usable and combinable forms, which lets us generate versatile solvers from little set of Paraiso source codes. We demonstrate Paraiso by implementing a compressive hydrodynamics solver. A single source code less than 500 lines can be used to generate solvers of arbitrary dimensions, for both multicore CPUs and GPUs. We demonstrate both manual annotation based tuning and evolutionary computing based automated tuning of the program.

We propose Paraiso, a domain specific language embedded in functional programming language Haskell, for automated tuning of explicit solvers of partial differential equations (PDEs) on GPUs as well as multicore CPUs. In Paraiso, one can describe PDE solving algorithms succinctly using tensor equations notation. Hydrodynamic properties, interpolation methods and other building blocks are described in abstract, modular, re-usable and combinable forms, which lets us generate versatile solvers from little set of Paraiso source codes. We demonstrate Paraiso by implementing a compressive hydrodynamics solver. A single source code less than 500 lines can be used to generate solvers of arbitrary dimensions, for both multicore CPUs and GPUs. We demonstrate both manual annotation based tuning and automated tuning of the program.

A critical challenge in measuring the power spectrum of 21cm emission from cosmic reionization is compensating for the frequency dependence of an interferometer’s sampling pattern, which can cause smooth-spectrum foregrounds to appear unsmooth and degrade the separation between foregrounds and the target signal. In this paper, we present an approach to foreground removal that explicitly accounts for this frequency dependence. We apply the delay transformation introduced in Parsons & Backer (2009) to each baseline of an interferometer to concentrate smooth-spectrum foregrounds within the bounds of the maximum geometric delays physically realizable on that baseline. By focusing on delay-modes that correspond to image-domain regions beyond the horizon, we show that it is possible to avoid the bulk of smooth-spectrum foregrounds. We show that delay-modes that are uncorrupted by foregrounds also represent samples of the three-dimensional power spectrum, and can be used to constrain cosmic reionization. Because it uses only spectral smoothness to differentiate foregrounds from the targeted 21cm signature, this per-baseline analysis approach relies on spectrally- and spatially-smooth instrumental responses for foreground removal. For sufficient levels of instrumental smoothness relative to the brightness of interfering foregrounds, this technique substantially reduces the level of calibration previously thought necessary to detect 21cm reionization. As a result, this approach places fewer constraints on antenna configuration within an array, and in particular, facilitates the adoption of configurations that are optimized for power-spectrum sensitivity. Under these assumptions, we demonstrate the potential for the PAPER array to detect 21cm reionization at an amplitude of 10 mK^2 near k~0.2h Mpc^-1 with 128 dipoles in 7 months of observing.

There is a growing trend toward using high-level tools for design and implementation of radio astronomy digital signal processing (DSP) systems. Such tools, for example, those from the Collaboration for Astronomy Signal Processing and Electronics Research (CASPER), are usually platform-specific, and lack high-level, platform-independent, portable, scalable application specifications. This limits the designer’s ability to experiment with designs at a high-level of abstraction and early in the development cycle. We address some of these issues using a model-based design approach employing dataflow models. We demonstrate this approach by applying it to the design of a tunable digital downconverter (TDD) used for narrow-bandwidth spectroscopy. Our design is targeted toward an FPGA platform, called the Interconnect Break-out Board (IBOB), that is available from the CASPER. We use the term TDD to refer to a digital downconverter for which the decmation factor and center frequency can be reconfigured without the need for regenerating the hardware code. Such a design is currently not available in the CASPER DSP library. The work presented in this paper focuses on two aspects. Firstly, we introduce and demonstrate a dataflow-based design approach using the dataflow interchange format (DIF) tool for high-level application specification, and we integrate this approach with the CASPER tool flow. Secondly, we explore the trade-off between the flexibility of TDD designs and the low hardware cost of fixed-configuration digital downconverter (FDD) designs that use the available CASPER DSP library. We further explore this trade-off in the context of a two-stage downconversion scheme employing a combination of TDD or FDD designs.

A number of physical processes show some form of bifurcation or splintering around a given point. Recently, it has been noted that cavity searches for interactions between photons and exotic fields may also show a bifurcation[1]. This paper describes the simulation of bifurcation of an optical beam in the presence of periodic focusing. The work is applicable to searches for exotic particles that employ cavities to extend the optical path of a laser beam through an interacting field. Significantly, simulations of cavity bifurcation reveal both a redistribution of the beam’s energy density and a shifting of the beam’s profile.

We present measurements of the total radio flux density as well as very-long-baseline interferometry (VLBI) images of the star, IM Pegasi, which was used as the guide star for the NASA/Stanford relativity mission Gravity Probe B. We obtained flux densities and images from 35 sessions of observations at 8.4 GHz (wavelength = 3.6 cm) between 1997 January and 2005 July. The observations were accurately phase-referenced to several extragalactic reference sources, and we present the images in a star-centered frame, aligned by the position of the star as derived from our fits to its orbital motion, parallax, and proper motion. Both the flux density and the morphology of IM Peg are variable. For most sessions, the emission region has a single-peaked structure, but 25% of the time, we observed a two-peaked (and on one occasion perhaps a three-peaked) structure. On average, the emission region is elongated by 1.4 +- 0.4 mas (FWHM), with the average direction of elongation being close to that of the sky projection of the orbit normal. The average length of the emission region is approximately equal to the diameter of the primary star. No significant correlation with the orbital phase is found for either the flux density or the direction of elongation, and no preference for any particular longitude on the star is shown by the emission region.

We present a physical interpretation for the locations of the sources of radio emission in IM Pegasi (IM Peg, HR 8703), the guide star for the NASA/Stanford relativity mission Gravity Probe B. This emission is seen in each of our 35 epochs of 8.4-GHz VLBI observations taken from 1997 to 2005. We found that the mean position of the radio emission is at or near the projected center of the primary to within about 27% of its radius, identifying this active star as the radio emitter. The positions of the radio brightness peaks are scattered across the disk of the primary and slightly beyond, preferentially along an axis with position angle, p.a. = (-38 +- 8) deg, which is closely aligned with the sky projections of the orbit normal (p.a. = -49.5 +- 8.6 deg) and the expected spin axis of the primary. Comparison with simulations suggests that brightness peaks are 3.6 (+0.4,-0.7) times more likely to occur (per unit surface area) near the pole regions of the primary (|latitude| >= 70 deg) than near the equator (|latitude| <= 20 deg), and to also occur close to the surface with ~2/3 of them at altitudes not higher than 25% of the radius of the primary.

We used 8.4 GHz VLBI images obtained at up to 35 epochs between 1997 and 2005 to examine the radio structures of the main reference source, 3C 454.3, and two secondary reference sources, B2250+194 and B2252+172, for the guide star for the NASA/Stanford relativity mission Gravity Probe B (GP-B). For one epoch in 2004 May, we also obtained images at 5.0 and 15.4 GHz. The 35 8.4 GHz images for quasar 3C 454.3 confirm a complex, evolving, core-jet structure. We identified at each epoch a component, C1, near the easternmost edge of the core region. Simulations of the core region showed that C1 is located, on average, 0.18 +- 0.06 mas west of the unresolved “core” identified in 43 GHz images. We also identified in 3C 454.3 at 8.4 GHz several additional components which moved away from C1 with proper motions ranging in magnitude between 0.9c and 5c. The detailed motions of the components exhibit two distinct bends in the jet axis located ~3 and ~5.5 mas west of C1. The spectra between 5.0 and 15.4 GHz for the “moving” components are steeper than that for C1. The 8.4 GHz images of B2250+194 and B2252+172, in contrast to those of 3C 454.3, reveal compact structures. The spectrum between 5.0 and 15.4 GHz for B2250+194 is inverted while that for B2252+172 is flat. Based on its position near the easternmost edge of the 8.4 GHz radio structure, close spatial association with the 43 GHz core, and relatively flat spectrum, we believe 3C 454.3 component C1 to be the best choice for the ultimate reference point for the GP-B guide star. The compact structures and inverted to flat spectra of B2250+194 and B2252+172 make these objects valuable secondary reference sources

When VLBI observations are used to determine the position or motion of a radio source relative to reference sources nearby on the sky, the astrometric information is usually obtained via: (i) phase-referenced maps; or (ii) parametric model fits to measured fringe phases or multiband delays. In this paper we describe a “merged” analysis technique which combines some of the most important advantages of these other two approaches. In particular, our merged technique combines the superior model-correction capabilities of parametric model fits with the ability of phase-referenced maps to yield astrometric measurements of sources that are too weak to be used in parametric model fits. We compare the results from this merged technique with the results from phase-referenced maps and from parametric model fits in the analysis of astrometric VLBI observations of the radio-bright star IM Pegasi (HR 8703) and the radio source B2252+172 nearby on the sky. In these studies we use central-core components of radio sources 3C 454.3 and B2250+194 as our positional references. We obtain astrometric results for IM Peg with our merged technique even when the source is too weak to be used in parametric model fits, and we find that our merged technique yields superior astrometric results to the phase-referenced mapping technique. We used our merged technique to estimate the proper motion and other astrometric parameters of IM Peg in support of the NASA/Stanford Gravity Probe B mission.

We present the principal astrometric results of the very-long-baseline interferometry (VLBI) program undertaken in support of the Gravity Probe B (GP-B) relativity mission. VLBI observations of the GP-B guide star, the RS CVn binary IM Pegasi (HR 8703), yielded positions at 35 epochs between 1997 and 2005. We discuss the statistical assumptions behind these results and our methods for estimating the systematic errors. We find the proper motion of IM Peg in an extragalactic reference frame closely related to the International Celestial Reference Frame 2 (ICRF2) to be -20.83 +- 0.03 +- 0.09 mas/yr in right ascension and -27.27 +- 0.03 +- 0.09 mas/yr in declination. For each component the first uncertainty is the statistical standard error and the second is the total standard error (SE) including plausible systematic errors. We also obtain a parallax of 10.37 +- 0.07 mas (distance: 96.4 +- 0.7 pc), for which there is no evidence of any significant contribution of systematic error. Our parameter estimates for the ~25-day-period orbital motion of the stellar radio emission have SEs corresponding to ~0.10 mas on the sky in each coordinate. The total SE of our estimate of IM Peg’s proper motion is ~30% smaller than the accuracy goal set by the GP-B project before launch: 0.14 mas/yr for each coordinate of IM Peg’s proper motion. Our results ensure that the uncertainty in IM Peg’s proper motion makes only a very small contribution to the uncertainty of the GP-B relativity tests.

We made VLBI observations at 8.4 GHz between 1997 and 2005 to estimate the coordinates of the “core” component of the superluminal quasar, 3C 454.3, the ultimate reference point in the distant universe for the NASA/Stanford Gyroscope Relativity Mission, Gravity Probe B. These coordinates are determined relative to those of the brightness peaks of two other compact extragalactic sources, B2250+194 and B2252+172, nearby on the sky, and within a celestial reference frame (CRF), defined by a large suite of compact extragalactic radio sources, and nearly identical to the International Celestial Reference Frame 2 (ICRF2). We find that B2250+194 and B2252+172 are stationary relative to each other, and also in the CRF, to within 1-sigma upper limits of 15 and 30 micro-arcsec/yr in RA and decl., respectively. The core of 3C 454.3 appears to jitter in its position along the jet direction over ~0.2 mas, likely due to activity close to the putative supermassive black hole nearby, but on average is stationary in the CRF within 1-sigma upper limits on its proper motion of 39 micro-arcsec/yr (1.0c) and 30 micro-arcsec/yr (0.8c) in RA and decl., respectively, for the period 2002 – 2005. Our corresponding limit over the longer interval, 1998 – 2005, of more importance to GP-B, is 46 and 56 micro-arcsec/yr in RA and decl., respectively. Some of 3C 454.3’s jet components show significantly superluminal motion with speeds of up to ~200 micro-arcsec/yr or 5c in the CRF. The core of 3C 454.3 thus provides for Gravity Probe B a sufficiently stable reference in the distant universe.

We describe the NASA/Stanford gyroscope relativity mission, Gravity Probe B (GP-B), and provide an overview of the following series of six astrometric and astrophysical papers that report on our radio observations and analyses made in support of this mission. The main goal of this 8.5 year program of differential VLBI astrometry was to determine the proper motion of the guide star of the GP-B mission, the RS CVn binary IM Pegasi (IM Peg; HR 8703). This proper motion is determined with respect to compact, extragalactic reference sources. The results are: -20.833 +- 0.090 mas/yr and -27.267 +- 0.095 mas/yr for, respectively, the right ascension and declination, in local Cartesian coordinates, of IM Peg’s proper motion, and 10.370 +- 0.074 mas (i.e., 96.43 +- 0.69 pc) for its parallax (and distance). Each quoted uncertainty is meant to represent an ~70% confidence interval that includes the estimated contribution from systematic error. These results are accurate enough not to discernibly degrade the GP-B estimates of its gyroscopes’ relativistic precessions: the frame-dragging and geodetic effects.

Observations with interferometric gravitational-wave detectors result in probability sky maps that are multimodal and spread over 10-100 deg^2. We present a scheme for maximizing the probability of imaging optical counterparts to gravitational-wave transients given limited observing resources. Our framework is capable of coordinating many telescopes with different fields of view and limiting magnitudes. We present a case study comparing three different planning algorithms. We find that, with the network of telescopes that was used in the most recent joint LIGO-Virgo science run, a relatively straightforward coordinated approach doubles the detection efficiency relative to each telescope observing independently.

We present a multi-scale, multi-wavelength source extraction algorithm called getsources. Although it has been designed primarily for use in the far-infrared surveys of Galactic star-forming regions with Herschel, the method can be applied to many other astronomical images. Instead of the traditional approach of extracting sources in the observed images, the new method analyzes fine spatial decompositions of original images across a wide range of scales and across all wavebands. It cleans those single-scale images of noise and background, and constructs wavelength-independent single-scale detection images that preserve information in both spatial and wavelength dimensions. Sources are detected in the combined detection images by following the evolution of their segmentation masks across all spatial scales. Measurements of the source properties are done in the original background-subtracted images at each wavelength; the background is estimated by interpolation under the source footprints and overlapping sources are deblended in an iterative procedure. In addition to the main catalog of sources, various catalogs and images are produced that aid scientific exploitation of the extraction results. We illustrate the performance of getsources on Herschel images by extracting sources in sub-fields of the Aquila and Rosette star-forming regions. The source extraction code and validation images with a reference extraction catalog are freely available.

The new CRISP filter at the Swedish 1-m Solar Telescope provides opportunities for observing the solar atmosphere with unprecedented spatial resolution and cadence. In order to benefit from the high quality of observational data from this instrument, we have developed methods for calibrating and restoring polarized Stokes images, obtained at optical and near infrared wavelengths, taking into account field-of-view variations of the filter properties. In order to facilitate velocity measurements, a time series from a 3D hydrodynamical granulation simulation is used to compute quiet Sun spectral line profiles at different heliocentric angles. The synthetic line profiles, with their convective blueshifts, can be used as absolute references for line-of-sight velocities. Observations of the Ca II 8542 {\AA} line are used to study magnetic fields in chromospheric fibrils. The line wings show the granulation pattern at mid-photospheric heights whereas the overlying chromosphere is seen in the core of the line. Using full Stokes data, we have attempted to observationally verify the alignment of chromospheric fibrils with the magnetic field. Our results suggest that in most cases fibrils are aligned along the magnetic field direction, but we also find examples where this is not the case. Detailed interpretation of Stokes data from spectral lines formed in the chromospheric data can be made using non-LTE inversion codes. For the first time, we use a realistic 3D MHD chromospheric simulation of the quiet Sun to assess how well NLTE inversions recover physical quantities from spectropolarimetric observations of Ca II 8542 {\AA}. We demonstrate that inversions provide realistic estimates of depth-averaged quantities in the chromosphere, although high spectral resolution and high sensitivity are needed to measure quiet Sun chromospheric magnetic fields.

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: 1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; 2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image reconstruction process was “blind” in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with 6 to 10 telescopes, should be well suited to provide images of complex sceneries.

In order to maximize the sensitivity of pulsar timing arrays to a stochastic gravitational wave background, we present computational techniques to optimize observing schedules. The techniques are applicable to both single and multi-telescope experiments. The observing schedule is optimized for each telescope by adjusting the observing time allocated to each pulsar while keeping the total amount of observing time constant. The optimized schedule depends on the timing noise characteristics of each individual pulsar as well as the performance of instrumentation. Several examples are given to illustrate the effects of different types of noise. A method to select the most suitable pulsars to be included in a pulsar timing array project is also presented.

Supermassive black hole binaries (SMBHBs) are expected to emit continuous gravitational waves in the pulsar timing array (PTA) frequency band ($10^{-9}$–$10^{-7}$ Hz). The development of data analysis techniques aimed at efficient detection and characterization of these signals is critical to the gravitational wave detection effort. In this paper we leverage methods developed for LIGO continuous wave gravitational searches, and explore the use of the $\mathcal{F}$-statistic for such searches in pulsar timing data. Babak & Sesana 2012 have already used this approach in the context of PTAs to show that one can resolve multiple SMBHB sources in the sky. Our work improves on several aspects of prior continuous wave search methods developed for PTA data analysis. The algorithm is implemented fully in the time domain, which naturally deals with the irregular sampling typical of PTA data and avoids spectral leakage problems associated with frequency domain methods. We take into account the fitting of the timing model, and have generalized our approach to deal with both correlated and uncorrelated colored noise sources. We also develop an incoherent detection statistic that maximizes over all pulsar dependent contributions to the likelihood. To test the effectiveness and sensitivity of our detection statistics, we perform a number of monte-carlo simulations. We produce sensitivity curves for PTAs of various configurations, and outline an implementation of a fully functional data analysis pipeline. Finally, we present a derivation of the likelihood maximized over the gravitational wave phases at the pulsar locations, which results in a vast reduction of the search parameter space.

We present the results of a machine-learning (ML) based search for new R Coronae Borealis (RCB) stars and DY Persei-like stars (DYPers) in the Galaxy using cataloged light curves obtained by the All-Sky Automated Survey (ASAS). RCB stars – a rare class of hydrogen-deficient carbon-rich supergiants – are of great interest owing to the insights they can provide on the late stages of stellar evolution. DYPers are possibly the low-temperature, low-luminosity analogs to the RCB phenomenon, though additional examples are needed to fully establish this connection. While RCB stars and DYPers are traditionally identified by epochs of extreme dimming that occur without regularity, the ML search framework more fully captures the richness and diversity of their photometric behavior. We demonstrate that our ML method recovers ASAS candidates that would have been missed by traditional search methods employing hard cuts on amplitude and periodicity. Our search yields 13 candidates that we consider likely RCB stars/DYPers: new and archival spectroscopic observations confirm that four of these candidates are RCB stars and four are DYPers. Our discovery of four new DYPers increases the number of known Galactic DYPers from two to six; noteworthy is that one of the new DYPers has a measured parallax and is m ~ 7 mag, making it the brightest known DYPer to date. Future observations of these new DYPers should prove instrumental in establishing the RCB connection. We consider these results, derived from a machine-learned probabilistic classification catalog, as an important proof-of-concept for the efficient discovery of rare sources with time-domain surveys.

With growing data volumes from synoptic surveys, astronomers must become more abstracted from the discovery and introspection processes. Given the scarcity of follow-up resources, there is a particularly sharp onus on the frameworks that replace these human roles to provide accurate and well-calibrated probabilistic classification catalogs. Such catalogs inform the subsequent follow-up, allowing consumers to optimize the selection of specific sources for further study and permitting rigorous treatment of purities and efficiencies for population studies. Here, we describe a process to produce a probabilistic classification catalog of variability with machine learning from a multi-epoch photometric survey. In addition to producing accurate classifications, we show how to estimate calibrated class probabilities, and motivate the importance of probability calibration. We also introduce a methodology for feature-based anomaly detection, which allows discovery of objects in the survey that do not fit within the predefined class taxonomy. Finally, we apply these methods to sources observed by the All Sky Automated Survey (ASAS), and unveil the Machine-learned ASAS Classification Catalog (MACC), which is a 28-class probabilistic classification catalog of 50,124 ASAS sources. We estimate that MACC achieves a sub-20% classification error rate, and demonstrate that the class posterior probabilities are reasonably calibrated. MACC classifications compare favorably to the classifications of several previous domain-specific ASAS papers and to the ASAS Catalog of Variable Stars, which had classified only 24% of those sources into one of 12 science classes. The MACC is publicly available at http://www.bigmacc.info.

With growing data volumes from synoptic surveys, astronomers must become more abstracted from the discovery and introspection processes. Given the scarcity of follow-up resources, there is a particularly sharp onus on the frameworks that replace these human roles to provide accurate and well-calibrated probabilistic classification catalogs. Such catalogs inform the subsequent follow-up, allowing consumers to optimize the selection of specific sources for further study and permitting rigorous treatment of purities and efficiencies for population studies. Here, we describe a process to produce a probabilistic classification catalog of variability with machine learning from a multi-epoch photometric survey. In addition to producing accurate classifications, we show how to estimate calibrated class probabilities, and motivate the importance of probability calibration. We also introduce a methodology for feature-based anomaly detection, which allows discovery of objects in the survey that do not fit within the predefined class taxonomy. Finally, we apply these methods to sources observed by the All Sky Automated Survey (ASAS), and unveil the Machine-learned ASAS Classification Catalog (MACC), which is a 28-class probabilistic classification catalog of 50,124 ASAS sources. We estimate that MACC achieves a sub-20% classification error rate, and demonstrate that the class posterior probabilities are reasonably calibrated. MACC classifications compare favorably to the classifications of several previous domain-specific ASAS papers and to the ASAS Catalog of Variable Stars, which had classified only 24% of those sources into one of 11 science classes. The MACC is publicly available at http://www.bigmacc.info.

In this paper we present results from the Mapping Dark Matter competition that expressed the weak lensing shape measurement task in its simplest form and as a result attracted over 700 submissions in 2 months and a factor of 3 improvement in shape measurement accuracy on high signal to noise galaxies, over previously published results, and a factor 10 improvement over methods tested on constant shear blind simulations. We also review weak lensing shape measurement challenges, including the Shear TEsting Programmes (STEP1 and STEP2) and the GRavitational lEnsing Accuracy Testing competitions (GREAT08 and GREAT10).

The Light-Time Effect (LTE) is observed whenever the distance between the observer and any kind of periodic event changes in time. The usual cause of this distance change is the reflex motion about the system’s barycenter due to the gravitational influence of one or more additional bodies. We analyze 5032 eclipsing contact (EC) and detached (ED) binaries from the All Sky Automated Survey (ASAS) catalogue to detect variations in the times of eclipses which possible can be due to the LTE effect. To this end we use an approach known from the radio pulsar timing where a template radio pulse of a pulsar is used as a reference to measure the times of arrivals of the collected pulses. In our analysis as a template for a photometric time series from ASAS, we use a best-fitting trigonometric series representing the light curve of a given EC or ED. Subsequently, an O-C diagram is built by comparing the template light curve with light curves obtained from subsets of a given time series. Most of the variations we detected in O-Cs correspond to a linear period change. Three show evidence of more than one complete LTE-orbit. For these objects we obtained preliminary orbital solutions. Our results demonstrate that the timing analysis employed in radio pulsar timing can be effectively used to study large data sets from photometric surveys.

In this paper we describe the state-of-the art status of multi-frequency detection techniques for compact sources in microwave astronomy. From the simplest cases where the spectral behaviour is well-known (i.e. thermal SZ clusters) to the more complex cases where there is little a priori information (i.e. polarized radio sources) we will review the main advances and the most recent results in the detection problem.

We explore the sensitivity and completeness of long baseline interferometric observations for detecting unknown, faint companions around bright unresolved stars. We derive a linear expression for the closure phase signature of a faint companion in the high contrast regime (<0.1), and provide a quantitative estimation of the detection efficiency for the currently offered four-telescope configurations at the Very Large Telescope Interferometer. The results are compared to the performances provided by linear and Y-shaped interferometric configurations in order to identify the ideal array. We find that all configurations have a similar efficiency in discovering companions wider than 10mas. Assuming a closure phase accuracy of 0.25deg, that is typical of state-of-the-art instruments, we predict a median dynamic range of up to six magnitudes when stacking observations obtained at five different hour angles. Surveying bright stars to search for faint companions can be considered as an ideal filler programme for modern interferometric facilities because that places few constraints on the choice of the interferometric configuration.

We examine stochastic variability in the dynamics of X-ray emission from the black hole system GRS 1915+105, a strongly variable microquasar commonly used for studying relativistic jets and the physics of black hole accretion. The analysis of sample observations for 13 different states in both soft (low) and hard (high) energy bands is performed by flicker-noise spectroscopy (FNS), a phenomenological time series analysis method operating on structure functions and power spectrum estimates. We find the values of FNS parameters, including the Hurst exponent, flicker-noise parameter, and characteristic time scales, for each observation based on multiple 2,500-second continuous data segments. We identify four modes of stochastic variability driven by dissipative processes that may be related to viscosity fluctuations in the accretion disk around the black hole: random (RN), power-law (1F), one-scale (1S), and two-scale (2S). The variability modes are generally the same in soft and hard energy bands of the same observation. We discuss the potential for future FNS studies of accreting black holes.

With its unprecedented light-collecting area for night-sky observations, the Cherenkov Telescope Array (CTA) holds great potential for also optical stellar astronomy, in particular as a multi-element intensity interferometer for realizing imaging with sub-milliarcsecond angular resolution. Such an order-of-magnitude increase of the spatial resolution achieved in optical astronomy will reveal the surfaces of rotationally flattened stars with structures in their circumstellar disks and winds, or the gas flows between close binaries. Image reconstruction is feasible from the second-order coherence of light, measured as the temporal correlations of arrival times between photons recorded in different telescopes. This technique (once pioneered by Hanbury Brown and Twiss) connects telescopes only with electronic signals and is practically insensitive to atmospheric turbulence and to imperfections in telescope optics. Detector and telescope requirements are very similar to those for imaging air Cherenkov observatories, the main difference being the signal processing (calculating cross correlations between single camera pixels in pairs of telescopes). Observations of brighter stars are not limited by sky brightness, permitting efficient CTA use during also bright-Moon periods. While other concepts have been proposed to realize kilometer-scale optical interferometers of conventional amplitude (phase-) type, both in space and on the ground, their complexity places them much further into the future than CTA, which thus could become the first kilometer-scale optical imager in astronomy.

In 2011, the AAVSO conducted a survey of 615 people who are or were recently active in the 101-year old organization. The survey included questions about their demographic background and variable star interests. Data are descriptively analyzed and compared with prior surveys. Results show an organization of very highly educated, largely male amateur and professional astronomers distributed across 108 countries. Participants tend to be loyal, with the average time of involvement in the AAVSO reported as 14 years. Most major demographic factors have not changed much over time. However, the average age of new members is increasing. Also, a significant portion of the respondents report being strictly active in a non-observing capacity, reflecting the growing mission of the organization. Motivations of participants are more aligned with scientific contribution than with that reported by other citizen science projects. This may help explain why a third of all respondents are an author or co-author of a paper in an astronomical journal. Finally, there is some evidence that participation in the AAVSO has a greater impact on the respondents’ view of their role in astronomy compared to that expected through increasing amateur astronomy experience alone. Results paint a picture of participants in a modern, advanced citizen science organization.

The Lunar University Network for Astrophysics Research (LUNAR) is a team of researchers and students at leading universities, NASA centers, and federal research laboratories undertaking investigations aimed at using the Moon as a platform for space science. LUNAR research includes Lunar Interior Physics & Gravitation using Lunar Laser Ranging (LLR), Low Frequency Cosmology and Astrophysics (LFCA), Planetary Science and the Lunar Ionosphere, Radio Heliophysics, and Exploration Science. The LUNAR team is exploring technologies that are likely to have a dual purpose, serving both exploration and science. There is a certain degree of commonality in much of LUNAR’s research. Specifically, the technology development for a lunar radio telescope involves elements from LFCA, Heliophysics, Exploration Science, and Planetary Science; similarly the drilling technology developed for LLR applies broadly to both Exploration and Lunar Science.

The Lunar University Network for Astrophysics Research (LUNAR) is a team of researchers and students at leading universities, NASA centers, and federal research laboratories undertaking investigations aimed at using the Moon as a platform for space science. LUNAR research includes Lunar Interior Physics & Gravitation using Lunar Laser Ranging (LLR), Low Frequency Cosmology and Astrophysics (LFCA), Planetary Science and the Lunar Ionosphere, Radio Heliophysics, and Exploration Science. The LUNAR team is exploring technologies that are likely to have a dual purpose, serving both exploration and science. There is a certain degree of commonality in much of LUNAR’s research. Specifically, the technology development for a lunar radio telescope involves elements from LFCA, Heliophysics, Exploration Science, and Planetary Science; similarly the drilling technology developed for LLR applies broadly to both Exploration and Lunar Science.

SaVi, a program for visualizing satellite orbits, movement, and coverage, is maintained at the University of Surrey. This tool has been used for research in academic papers, and by industry companies designing and intending to deploy satellite constellations. It has also proven useful for demonstrating aspects of satellite constellations and their geometry, coverage and movement for educational and teaching purposes. SaVi is introduced and described briefly here.

A networking transport protocol, named Saratoga, has been developed at the University of Surrey for efficient delivery of imagery from Internet-Protocol-based remote-sensing satellites. Saratoga is now being implemented and evaluated for use for the high-end data-delivery needs of astronomers using large, advanced, radio telescopes. These telescopes are expected to take advantage of Internet technologies. This brief paper outlines the reasons for the creation and adoption of this protocol, discusses how it differs from and complements other protocols, and summarises the worldwide collaboration that is making this development possible.

The Cisco router in Low Earth Orbit (CLEO) was launched into space as an experimental secondary payload onboard the UK Disaster Monitoring Constellation (UK-DMC) satellite in September 2003. The UK-DMC satellite is one of an increasing number of DMC satellites in orbit that rely on the Internet Protocol (IP) for command and control and for delivery of data from payloads. The DMC satellites, built by Surrey Satellite Technology Ltd (SSTL), have imaged the effects of Hurricane Katrina, the Indian Ocean Tsunami, and other events for disaster relief under the International Space and Major Disasters Charter. It was possible to integrate the Cisco mobile access router into the UK-DMC satellite as a result of the DMC satellites’ adoption of existing commercial networking standards, using IP over Frame Relay over standard High-Level Data Link Control, or HDLC (ISO 13239) on standard serial interfaces. This approach came from work onboard SSTL’s earlier UoSAT-12 satellite

The realization of low-cost instruments with high technical performance is a goal which deserves some efforts in an epoch of fast technological developments: indeed such instruments can be easily reproduced and therefore allow to open new research programs in several Observatories. We realized a fast optical photometer based on the SiPM technology, using commercially available modules. Using low-cost components we have developed a custom electronic chain to extract the signal produced by a commercial MPPC module produced by Hamamatsu, in order to obtain sub millisecond sampling of the light curve of astronomical sources, typically pulsars. In the early February 2011 we observed the Crab Pulsar at the Cassini telescope with our prototype photometer, deriving its period, power spectrum and shape of its light curve in very good agreement with the results obtained in the past with other instruments.

We report on the development of W-band (75-110 GHz) heterodyne receiver technology for large-format astronomical arrays. The receiver system is designed to be both mass-producible, so that the designs could be scaled to thousands of receiver elements, and modular. Most of the receiver function- ality is integrated into compact Monolithic Microwave Integrated Circuit (MMIC) amplifier-based multichip modules. The MMIC modules include a chain of InP MMIC low-noise amplifiers, coupled-line bandpass filters and sub-harmonic Schottky diode mixers. The receiver signals will be routed to and from the MMIC modules on a multilayer high frequency laminate, which includes splitters, amplifiers, and frequency doublers. A prototype MMIC module has exhibited a band-averaged noise temperature of 41 K from 82-100 GHz and a gain of 29 dB at 15 K, which is the state- of-the-art for heterodyne multi-chip modules.

Using the planned antenna locations of the 128 antenna buildout of the Murchison Widefield Array (MWA), we accurately calculate its sensitivity to the Epoch of Reionization (EoR) power spectrum of redshifted 21 cm emission. Our calculation takes into account synthesis rotation, chromatic and asymmetrical baseline effects, and excludes modes that will be contaminated by foreground subtraction. With one full season of observation on two fields (900 and 700 hours), the MWA will be capable of a 14$\sigma$ detection of the EoR signal along with slope constraints.

Using the planned antenna locations of the 128 antenna buildout of the Murchison Widefield Array (MWA), we accurately calculate its sensitivity to the Epoch of Reionization (EoR) power spectrum of redshifted 21 cm emission. Our calculation takes into account synthesis rotation, chromatic and asymmetrical baseline effects, and excludes modes that will be contaminated by foreground subtraction. With one full season of observation on two fields (900 and 700 hours), the MWA will be capable of a 14$\sigma$ detection of the EoR signal along with slope constraints.

In this short review we present the developments over the last 5 decades that have led to the use of Graphics Processing Units (GPUs) for astrophysical simulations. Since the introduction of NVIDIA’s Compute Unified Device Architecture (CUDA) in 2007 the GPU has become a valuable tool for N-body simulations and is so popular these days that almost all papers about high precision N-body simulations use methods that are accelerated by GPUs. With the GPU hardware becoming more advanced and being used for more advanced algorithms like gravitational tree-codes we see a bright future for GPU like hardware in computational astrophysics.

We present a data model describing the structure of spectrophotometric datasets with spectral and temporal coordinates and associated metadata. This data model may be used to represent spectra, time series data, segments of SED (Spectral Energy Distributions) and other spectral or temporal associations.

The newsletter of the Australian Astronomical Observatory. In this issue: Using 2dF and AAOmega to Harness the Full Power of the Supernova Legacy Survey; Emission Lines in the Near Infrared: Tracing the Violent ISM; Dancing Starbugs: vacuum adhesion, field rotation and other progress; A message of progress from HERMES; Imaging with the 2dF Focal Plane Imager; and all the usual columns and news from the Observatory.

Intermediate-mass stars end their lives by ejecting the bulk of their envelope via a slow dense wind back into the interstellar medium, to form the next generation of stars and planets. Stellar pulsations are thought to elevate gas to an altitude cool enough for the condensation of dust, which is then accelerated by radiation pressure from starlight, entraining the gas and driving the wind. However accounting for the mass loss has been a problem due to the difficulty in observing tenuous gas and dust tens of milliarcseconds from the star, and there is accordingly no consensus on the way sufficient momentum is transferred from the starlight to the outflow. Here, we present spatially-resolved, multi-wavelength observations of circumstellar dust shells of three stars on the asymptotic giant branch of the HR diagram. When imaged in scattered light, dust shells were found at remarkably small radii (<~ 2 stellar radii) and with unexpectedly large grains (~300 nm radius). This proximity to the photosphere argues for dust species that are transparent to starlight and therefore resistant to sublimation by the intense radiation field. While transparency usually implies insufficient radiative pressure to drive a wind, the radiation field can accelerate these large grains via photon scattering rather than absorption – a plausible mass-loss mechanism for lower-amplitude pulsating stars.

Optical interferometry provides us with a unique opportunity to improve our understanding of stellar structure and evolution. Through direct observation of rotationally distorted photospheres at sub-milliarcsecond scales, we are now able to characterize latitude dependencies of stellar radius, temperature structure, and even energy transport. These detailed new views of stars are leading to revised thinking in a broad array of associated topics, such as spectroscopy, stellar evolution, and exoplanet detection. As newly advanced techniques and instrumentation mature, this topic in astronomy is poised to greatly expand in depth and influence.

We investigate the possibility that the anomalous acceleration of the Pioneer 10 and 11 spacecraft is due to the recoil force associated with an anisotropic emission of thermal radiation off the vehicles. To this end, relying on the project and spacecraft design documentation, we constructed a comprehensive finite-element thermal model of the two spacecraft. Then, we numerically solve thermal conduction and radiation equations using the actual flight telemetry as boundary conditions. We use the results of this model to evaluate the effect of the thermal recoil force on the Pioneer 10 spacecraft at various heliocentric distances. We found that the magnitude, temporal behavior, and direction of the resulting thermal acceleration are all similar to the properties of the observed anomaly. As a novel element of our investigation, we develop a parameterized model for the thermal recoil force and estimate the coefficients of this model independently from navigational Doppler data. We find no statistically significant difference between the two estimates and conclude that once the thermal recoil force is properly accounted for, no anomalous acceleration remains.

We investigate the suitability of the Wendland functions as smoothing kernels for smoothed particle hydrodynamics (SPH) and compare them with the traditional B-splines. Linear stability analysis in three dimensions and test simulations demonstrate that the Wendland kernels avoid the clumping (or pairing) instability for all neighbour numbers NH, despite having vanishing derivative at the origin. This disproves traditional ideas about the origin of this instability. Instead, we give an explanation based on the kernel Fourier transform, but also an interpretation in terms of the SPH density estimator. The Wendland kernels are computationally more convenient than the higher-order B-splines and thus allow large NH, which we show are required to obtain decent numerical accuracy for strongly shearing flows (note that computational costs rise sub-linear with NH). At low NH the quartic B-spline kernel with NH = 60 obtains much better convergence then the standard cubic B-spline with NH<=57.

Supplementing the publications based on the first-light observations with the German Receiver for Astronomy at Terahertz frequencies (GREAT) on SOFIA, we present background information on the underlying heterodyne detector technology. We describe the superconducting hot electron bolometer (HEB) detectors that are used as frequency mixers in the L1 (1400 GHz), L2 (1900 GHz), and M (2500 GHz) channels of GREAT. Measured performance of the detectors is presented and background information on their operation in GREAT is given. Our mixer units are waveguide-based and couple to free-space radiation via a feedhorn antenna. The HEB mixers are designed, fabricated, characterized, and flight-qualified in-house. We are able to use the full intermediate frequency bandwidth of the mixers using silicon-germanium multi-octave cryogenic low-noise amplifiers with very low input return loss. Superconducting HEB mixers have proven to be practical and sensitive detectors for high-resolution THz frequency spectroscopy on SOFIA. We show that our niobium-titanium-nitride (NbTiN) material HEBs on silicon nitride (SiN) membrane substrates have an intermediate frequency (IF) noise roll-off frequency above 2.8 GHz, which does not limit the current receiver IF bandwidth. Our mixer technology development efforts culminate in the first successful operation of a waveguide-based HEB mixer at 2.5 THz and deployment for radioastronomy. A significant contribution to the success of GREAT is made by technological development, thorough characterization and performance optimization of the mixer and its IF interface for receiver operation on SOFIA. In particular, the development of an optimized mixer IF interface contributes to the low passband ripple and excellent stability, which GREAT demonstrated during its initial successful astronomical observation runs.

We present a gravitational hierarchical N-body code that is designed to run efficiently on Graphics Processing Units (GPUs). All parts of the algorithm are executed on the GPU which eliminates the need for data transfer between the Central Processing Unit (CPU) and the GPU. Our tests indicate that the gravitational tree-code outperforms tuned CPU code for all parts of the algorithm and show an overall performance improvement of more than a factor 20, resulting in a processing rate of more than 2.8 million particles per second.

In this paper we discuss a simple method of testing for the presence of energy-dependent dispersion in high energy data-sets. It uses the minimisation of the Kolmogorov distance between the cumulative distribution of two probability functions as the statistical metric to estimate the magnitude of any spectral dispersion within transient features in a light-curve and we also show that it performs well in the presence of modest energy resolutions (~20%) typical of gamma-ray observations. After presenting the method in detail we apply it to a parameterised simulated lightcurve based on the extreme VHE gamma-ray flare of PKS 2155-304 observed with H.E.S.S. in 2006, in order to illustrate its potential through the concrete example of setting constraints on quantum-gravity induced Lorentz invariance violation (LIV) effects. We obtain comparable limits to those of the most advanced techniques used in LIV searches applied to similar datasets, but the present method has the advantage of being particularly straightforward to use. Whilst the development of the method was motivated by LIV searches, it is also applicable to other astrophysical situations where energy-dependent dispersion is expected, such as spectral lags from the acceleration and cooling of particles in relativistic outflows.

NGA/eLISA is a new mission proposal with arm length 106 km and one interferometer down-scaled from LISA (http://elisa-ngo.org/). Just like LISA and ASTROD-GW, in order to attain the requisite sensitivity for NGO/eLISA, laser frequency noise must be suppressed below the secondary noises such as the optical path noise, acceleration noise etc. In previous papers, we have used the CGC 2.7 ephemeris to numerically simulate the time delay interferometry for LISA and ASTROD-GW with one arm dysfunctional and found that they are both well below their respective limits under which the laser frequency noise is required to be suppressed. In this paper, we follow the same procedure to simulate the time delay interferometry numerically. To do this, we work out a set of 1000-day optimized mission orbits of NGO/eLISA spacecraft starting at January 1st, 2021 using the CGC 2.7 ephemeris framework. We then use this numerical solution to calculate the residual optical path differences in the second-generation solutions of our previous papers. The maximum path length difference, for all configuration calculated, is below 12 mm (40 ps). This is well below the limit under which the laser frequency noise is required to be suppressed for NGO/eLISA. We compare and discuss the resulting differences due to different arm lengths for various mission proposals — NGO/eLISA, an NGO-LISA-type mission with a nominal arm length of 2 \times 10^6 km, LISA and ASTROD-GW.

V921 Scorpii is a close binary system (separation 0.025″) showing the B[e]-phenomenon. The system is surrounded by an enigmatic bipolar nebula, which might have been shaped by episodic mass-loss events, possibly triggered by dynamical interactions between the companion and the circumprimary disk (Kraus et al. 2012a). In this paper, we investigate the spatial structure and kinematics of the circumprimary disk, with the aim to obtain new insights into the still strongly debated evolutionary stage. For this purpose, we combine, for the first time, infrared spectro-interferometry (VLTI/AMBER, R=12,000) and spectro-astrometry (VLT/CRIRES, R=100,000), which allows us to study the AU-scale distribution of circumstellar gas and dust with an unprecedented velocity resolution of 3 km*s^-1. Using a model-independent photocenter analysis technique, we find that the Br-gamma-line emission rotates in the same plane as the dust disk. We can reproduce the wavelength-differential visibilities and phases and the double-peaked line profile using a Keplerian-rotating disk model. The derived mass of the central star is 5.4+/-0.4 M_sun*(d/1150 pc), which is considerably lower than expected from the spectral classification, suggesting that V921 Sco might be more distant (d approx 2kpc) than commonly assumed. Using the geometric information provided by our Br-gamma spectro-interferometric data and Paschen, Brackett, and Pfund line decrement measurements in 61 hydrogen recombination line transitions, we derive the density of the line-emitting gas (N_e=2…6*10^19 m^-3). Given that our measurements can be reproduced with a Keplerian velocity field without outflowing velocity component and the non-detection of age-indicating spectroscopic diagnostics, our study provides new evidence for the pre-main-sequence nature of V921 Sco.

V921 Scorpii is a close binary system (separation 0.025″) showing the B[e]-phenomenon. The system is surrounded by an enigmatic bipolar nebula, which might have been shaped by episodic mass-loss events, possibly triggered by dynamical interactions between the companion and the circumprimary disk (Kraus et al. 2012a). In this paper, we investigate the spatial structure and kinematics of the circumprimary disk, with the aim to obtain new insights into the still strongly debated evolutionary stage. For this purpose, we combine, for the first time, infrared spectro-interferometry (VLTI/AMBER, R=12,000) and spectro-astrometry (VLT/CRIRES, R=100,000), which allows us to study the AU-scale distribution of circumstellar gas and dust with an unprecedented velocity resolution of 3 km…^-1. Using a model-independent photocenter analysis technique, we find that the Br{\gamma}-line emission rotates in the same plane as the dust disk. We can reproduce the wavelength-differential visibilities and phases and the double-peaked line profile using a Keplerian-rotating disk model. The derived mass of the central star is 5.4+/-0.4 M_sun\cdot(d/1150 pc), which is considerably lower than expected from the spectral classification, suggesting that V921 Sco might be more distant (d approx 2kpc) than commonly assumed. Using the geometric information provided by our Br-gamma spectro-interferometric data and Paschen, Brackett, and Pfund line decrement measurements in 61 hydrogen recombination line transitions, we derive the density of the line-emitting gas (N_e=2…6\cdot10^19 m^-3). Given that our measurements can be reproduced with a Keplerian velocity field without outflowing velocity component and the non-detection of age-indicating spectroscopic diagnostics, our study provides new evidence for the pre-main-sequence nature of V921 Sco.

We have entered an era of massive data sets in astronomy. In particular, the number of supernova (SN) discoveries and classifications has substantially increased over the years from few tens to thousands per year. It is no longer the case that observations of a few prototypical events encapsulate most spectroscopic information about SNe, motivating the development of modern tools to collect, archive, organize and distribute spectra in general, and SN spectra in particular. For this reason we have developed the Weizmann Interactive Supernova data REPository – WISeREP – an SQL-based database (DB) with an interactive web-based graphical interface. The system serves as an archive of high quality SN spectra, including both historical (legacy) data as well as data that is accumulated by ongoing modern programs. The archive provides information about objects, their spectra, and related meta-data. Utilizing interactive plots, we provide a graphical interface to visualize data, perform line identification of the major relevant species, determine object redshifts, classify SNe and measure expansion velocities. Guest users may view and download spectra or other data that have been placed in the public domain. Registered users may also view and download data that are proprietary to specific programs with which they are associated. The DB currently holds >8000 spectra, of which >5000 are public; the latter include published spectra from the Palomar Transient Factory, all of the SUSPECT archive, the Caltech-Core-Collapse Program, the CfA SN spectra archive and published spectra from the UC Berkeley SNDB repository. It offers an efficient and convenient way to archive data and share it with colleagues, and we expect that data stored in this way will be easy to access, increasing its visibility, usefulness and scientific impact.

We have entered an era of massive data sets in astronomy. In particular, the number of supernova (SN) discoveries and classifications has substantially increased over the years from few tens to thousands per year. It is no longer the case that observations of a few prototypical events encapsulate most spectroscopic information about SNe, motivating the development of modern tools to collect, archive, organize and distribute spectra in general, and SN spectra in particular. For this reason we have developed the Weizmann Institute of Science Experimental Astrophysics Spectroscopy System – WISEASS — an SQL-based database (DB) with an interactive web-based graphical interface. The system serves as an archive of high quality SN spectra, including both historical (legacy) data as well as data that is accumulated by ongoing modern programs. The archive provides information about objects, their spectra, and related meta-data. Utilizing interactive plots, we provide a graphical interface to visualize data, perform line identification of the major relevant species, determine object redshifts, classify SNe and measure expansion velocities. Guest users may view and download spectra or other data that have been placed in the public domain. Registered users may also view and download data that are proprietary to specific programs with which they are associated. The DB currently holds >7700 spectra, of which >4600 are public; the latter include published spectra from the Palomar Transient Factory (PTF), the Caltech-Core-Collapse Program (CCCP), all of the SUSPECT (SUpernova SPECTrum) archive and the CfA Type Ia SN spectral archive. It offers an efficient and convenient way to archive data and share it with colleagues, and we expect that data stored in this way will be easy to access, increasing its visibility, usefulness and scientific impact.

Cosmic ray electrons and positrons constitute an important component of the background for imaging atmospheric Cherenkov Telescope Systems with very low energy thresholds. As the primary energy of electrons and positrons decreases, their contribution to the background trigger rate dominates over protons, at least in terms of differential rates against actual energies. After event reconstruction, this contribution might become comparable to the proton background at energies of the order of few GeV. It is well known that the flux of low energy charged particles is suppressed by the Earth’s magnetic field. This effect strongly depends on the geographical location, the direction of incidence of the charged particle and its mass. Therefore, the geomagnetic field can contribute to diminish the rate of the electrons and positrons detected by a given array of Cherenkov Telescopes. In this work we study the propagation of low energy primary electrons in the Earth’s magnetic field by using the backtracking technique. We use a more realistic geomagnetic field model than the one used in previous calculations. We consider some sites relevant for new generations of imaging atmospheric Cherenkov Telescopes. We also study in detail the case of 5@5, a proposed low energy Cherenkov Telescope array.

The Tianshan Radio Experiment for Neutrino Detection (TREND) is a sino-french collaboration (CNRS/IN2P3 and Chinese Academy of Science) developing an autonomous antenna array for the detection of high energy Extensive Air Showers (EAS) on the site of the 21CMA radio observatory. The autonomous detection and identification of EAS was achieved by TREND on a prototype array in 2009. This result was confirmed soon after when EAS radio-candidates could be tagged as cosmic ray events by an array of particle detectors running in parallel at the same location. This result is an important milestone for TREND, and more generally, for the maturation of the EAS radio-detection technique. The array is presently composed of 50 antennas covering a total area of ~1.2 km^2, running in steady conditions since March 2011. We are presently processing the data to identify EAS radio-candidates. In a long term perspective, TREND is intended to search for high energy tau neutrinos. Here we only report on the results achieved so far by TREND.

A prototype of digital frequency multiplexing electronics allowing the real time monitoring of microwave kinetic inductance detector (MKIDs) arrays for mm-wave astronomy has been developed. Thanks to the frequency multiplexing, it can monitor simultaneously 400 pixels over a 500 MHz bandwidth and requires only two coaxial cables for instrumenting such a large array. The chosen solution and the performance achieved are presented in this paper.

The wavelength dependence of atmospheric refraction causes elongation of finite-bandwidth images along the elevation vector, which produces spurious signals in weak gravitational lensing shear measurements unless this atmospheric dispersion is calibrated and removed to high precision. Because astrometric solutions and point spread function (PSF) characteristics are typically calibrated from stellar images, differences between the reference stars’ spectra and the galaxies’ spectra will leave residual errors in both the astrometric positions ($\Delta{\bar{R}}$) and in the second moment (width) of the wavelength-averaged PSF ($\Delta{v}$) for galaxies. We estimate the level of $\Delta{V}$ that will induce spurious weak lensing signals in PSF-corrected galaxy shapes that exceed the statistical errors of the {\em Dark Energy Survey (DES)} and the {\em Large Synoptic Survey Telescope (LSST)} cosmic-shear experiments. We also estimate the $\Delta{\bar{R}}$ signals that will produce unacceptable spurious distortions after stacking of exposures taken at different airmasses and hour angles. Using standard galaxy and stellar spectral templates we calculate the resultant errors in the $griz$ bands, and find that atmospheric dispersion differentials, left uncorrected, exceed the {\em DES} cosmic-shear requirements in the $g$ and $r$ bands, and exceed the stricter LSST requirements in $i$ band. We find that a simple correction linear in galaxy color is accurate enough to recover the use of $r$ band for DES and $i$ band for LSST. More complex approaches to correction of the atmospheric dispersion signal will be needed to use the $g$ band for DES cosmic-shear measurements or to use the $g$ or $r$ bands for LSST cosmic-shear measurements.

The field of TeV gamma-ray astronomy has produced many exciting results over the last decade. Both the source catalogue, and the range of astrophysical questions which can be addressed, continue to expand. This article presents a topical review of the field, with a focus on the observational results of the imaging atmospheric Cherenkov telescope arrays. The results encompass pulsars and their nebulae, supernova remnants, gamma-ray binary systems, star forming regions and starburst and active galaxies.

The technique of High Spectral Resolution Lidar (HSRL) for atmospheric monitoring allows the determination of the aerosol to molecular ratio and can be used in UHECR Observatories using air fluorescence telescopes. By this technique a more accurate estimate of the Cherenkov radiation superimposed to the fluorescence signal can be achieved. A laboratory setup was developed to determine the backscattering coefficients using microparticles diluted in water and diffusion interfaces. In this setup we used a CW SLM laser at 532 nm and a 250 mm Newtonian telescope. Simulations of the above experimental configuration have been made using Scatlab\c{opyright}, FINESSE\c{opyright} 0.99.8 and MATLAB\c{opyright} and are presented in this work. We compare the simulated 2-dimensional Fabry-Perot fringe images of the backscattering signal recorded in the CCD sensor with that of experimental ones. Additionally, we simulated the backscattering of the laser beam by the atmosphere at a height of 2000 m and we have studied the influence of the beam and its diameter on the fringe image.

The AMBRE Project is a collaboration between the European Southern Observatory (ESO) and the Observatoire de la Cote d’Azur (OCA) that has been established in order to carry out the determination of stellar atmospheric parameters for the archived spectra of four ESO spectrographs. The analysis of the FEROS archived spectra for their stellar parameters (effective temperatures, surface gravities, global metallicities, alpha element to iron ratios and radial velocities) has been completed in the first phase of the AMBRE Project. From the complete ESO:FEROS archive dataset that was received, a total of 21551 scientific spectra have been identified, covering the period 2005 to 2010. These spectra correspond to ~6285 stars. The determination of the stellar parameters was carried out using the stellar parameterisation algorithm, MATISSE (MATrix Inversion for Spectral SynthEsis), which has been developed at OCA to be used in the analysis of large scale spectroscopic studies in galactic archaeology. An analysis pipeline has been constructed that integrates spectral reduction and radial velocity correction procedures with MATISSE in order to automatically determine the stellar parameters of the FEROS spectra. Stellar atmospheric parameters (Teff, log g, [M/H] and [alpha/Fe]) were determined for 6508 (30.2%) of the FEROS archived spectra (~3087 stars). Radial velocities were determined for 11963 (56%) of the archived spectra. 2370 (11%) spectra could not be analysed within the pipeline. 12673 spectra (58.8%) were analysed in the pipeline but their parameters were discarded based on quality criteria and error analysis determined within the automated process. The majority of these rejected spectra were found to have broad spectral features indicating that they may be hot and/or fast rotating stars, which are not considered within the adopted reference synthetic spectra grid of FGKM stars.

Reliable timing calibration is essential for the accurate comparison of XMM-Newton light curves with those from other observatories, to ultimately use them to derive precise physical quantities. The XMM-Newton timing calibration is based on pulsar analysis. However, as pulsars show both timing noise and glitches, it is essential to monitor these calibration sources regularly. To this end, the XMM-Newton observatory performs observations twice a year of the Crab pulsar to monitor the absolute timing accuracy of the EPIC-pn camera in the fast Timing and Burst modes. We present the results of this monitoring campaign, comparing XMM-Newton data from the Crab pulsar (PSR B0531+21) with radio measurements. In addition, we use five pulsars (PSR J0537-69, PSR B0540-69, PSR B0833-45, PSR B1509-58 and PSR B1055-52) with periods ranging from 16 ms to 197 ms to verify the relative timing accuracy. We analysed 38 XMM-Newton observations (0.2-12.0 keV) of the Crab taken over the first ten years of the mission and 13 observations from the five complementary pulsars. All the data were processed with the SAS, the XMM-Newton Scientific Analysis Software, version 9.0. Epoch folding techniques coupled with \chi^{2} tests were used to derive relative timing accuracies. The absolute timing accuracy was determined using the Crab data and comparing the time shift between the main X-ray and radio peaks in the phase folded light curves. The relative timing accuracy of XMM-Newton is found to be better than 10^{-8}. The strongest X-ray pulse peak precedes the corresponding radio peak by 306\pm9 \mus, which is in agreement with other high energy observatories such as Chandra, INTEGRAL and RXTE. The derived absolute timing accuracy from our analysis is \pm48 \mus.

The reliability of astronomical observations at millimeter and sub-millimeter wavelengths closely depends on a low vertical content of water vapor as well as on high atmospheric emission stability. Although Concordia station at Dome C (Antarctica) enjoys good observing conditions in this atmospheric spectral windows, as shown by preliminary site-testing campaigns at different bands and in, not always, time overlapped periods, a dedicated instrument able to continuously determine atmospheric performance for a wide spectral range is not yet planned. In the absence of such measurements, in this paper we suggest a semi-empirical approach to perform an analysis of atmospheric transmission and emission at Dome C to compare the performance for 7 photometric bands ranging from 100 GHz to 2 THz. Radiosoundings data provided by the Routine Meteorological Observations (RMO) Research Project at Concordia station are corrected by temperature and humidity errors and dry biases and then employed to feed ATM (Atmospheric Transmission at Microwaves) code to generate synthetic spectra in the wide spectral range from 100 GHz to 2 THz. To quantify the atmospheric contribution in millimeter and sub-millimeter observations we are considering several photometric bands in which atmospheric quantities are integrated. The observational capabilities of this site at all the selected spectral bands are analyzed considering monthly averaged transmissions joined to the corresponding fluctuations. Transmission and pwv statistics at Dome C derived by our semi-empirical approach are consistent with previous works. It is evident the decreasing of the performance at high frequencies. We propose to introduce a new parameter to compare the quality of a site at different spectral bands, in terms of high transmission and emission stability, the Site Quality Factor.

A wavelength calibration system based on a laser frequency comb (LFC) was developed in a co-operation between the Kiepenheuer-Institut f\”ur Sonnenphysik, Freiburg, Germany and the Max-Planck-Institut f\”ur Quantenoptik, Garching, Germany for permanent installation at the German Vacuum Tower Telescope (VTT) on Tenerife, Canary Islands. The system was installed successfully in October 2011. By simultaneously recording the spectra from the Sun and the LFC, for each exposure a calibration curve can be derived from the known frequencies of the comb modes that is suitable for absolute calibration at the meters per second level. We briefly summarize some topics in solar physics that benefit from absolute spectroscopy and point out the advantages of LFC compared to traditional calibration techniques. We also sketch the basic setup of the VTT calibration system and its integration with the existing echelle spectrograph.

The Rosse Solar-Terrestrial Observatory (RSTO; www.rosseobservatory.ie) was established at Birr Castle, Co. Offaly, Ireland (53 05′38.9″, 7 55′12.7″) in 2010 to study solar radio bursts and the response of the Earth’s ionosphere and geomagnetic field. To date, three Compound Astronomical Low-cost Low-frequency Instrument for Spectroscopy and Transportable Observatory (CALLISTO) spectrometers have been installed, with the capability of observing in the frequency range 10-870 MHz. The receivers are fed simultaneously by biconical and log-periodic antennas. Nominally, frequency spectra in the range 10-400 MHz are obtained with 4 sweeps per second over 600 channels. Here, we describe the RSTO solar radio spectrometer set-up, and present dynamic spectra of a sample of Type II, III and IV radio bursts. In particular, we describe fine-scale structure observed in Type II bursts, including band splitting and rapidly varying herringbone features.

The Rosse Solar-Terrestrial Observatory (RSTO; www.rosseobservatory.ie) was established at Birr Castle, Co. Offaly, Ireland (53{\deg}05′38.9″, 7{\deg}55′12.7″) in 2010 to study solar radio bursts and the response of the Earth’s ionosphere and geomagnetic field. To date, three Compound Astronomical Low-cost Low-frequency Instrument for Spectroscopy and Transportable Observatory (CAL- LISTO) spectrometers have been installed, with the capability of observing in the frequency range 10-870 MHz. The receivers are fed simultaneously by biconical and log-periodic antennas. Nominally, frequency spectra in the range 10-400 MHz are obtained with 4 sweeps per second over 600 channels. Here, we describe the RSTO solar radio spectrometer set-up, and present dynamic spectra of a sample of Type II, III and IV radio bursts. In particular, we describe fine-scale structure observed in Type II bursts, including band splitting and rapidly varying herringbone features.

In this paper, we present the acoustic transceiver developed for the positioning system in underwater neutrino telescopes. These infrastructures are not completely rigid and need a positioning system in order to monitor the position of the optical sensors of the telescope which have some degree of motion due to sea currents. To have a highly reliable and versatile system in the infrastructure, the transceiver has the requirements of reduced cost, low power consumption, high intensity for emission, low intrinsic noise, arbitrary signals for emission and the capacity of acquiring and processing the received signal on the board. The solution proposed and presented here consists of an acoustic transducer that works in the 20-40 kHz region and withstands high pressures (up to 500 bars). The electronic-board can be configured from shore and is able to feed the transducer with arbitrary signals and to control the transmitted and received signals with very good timing precision. The results of the different tests done on the transceiver in the laboratory are described here, as well as the change implemented for its integration in the Instrumentation Line of ANTARES for the in situ tests. We consider the transceiver design is so versatile that it may be used in other kinds of marine positioning systems, alone or combined with other marine systems, or integrated in different Earth-Sea Observatories, where the localization of the sensors is an issue.

In this paper acoustic transmitters that were developed for use in underwater neutrino telescopes are presented. Firstly, an acoustic transceiver has been developed as part of the acoustic positioning system of neutrino telescopes. These infrastructures are not completely rigid and require a positioning system in order to monitor the position of the optical sensors which move due to sea currents. To guarantee a reliable and versatile system, the transceiver has the requirements of reduced cost, low power consumption, high pressure withstanding (up to 500 bars), high intensity for emission, low intrinsic noise, arbitrary signals for emission and the capacity of acquiring and processing received signals. Secondly, a compact acoustic transmitter array has been developed for the calibration of acoustic neutrino detection systems. The array is able to mimic the signature of ultra-high-energy neutrino interaction in emission directivity and signal shape. The technique of parametric acoustic sources has been used to achieve the proposed aim. The developed compact array has practical features such as easy manageability and operation. The prototype designs and the results of different tests are described. The techniques applied for these two acoustic systems are so powerful and versatile that may be of interest in other marine applications using acoustic transmitters.

The development of multi-layer optics which allow to focus photons up to 100 keV and more promises an enormous jump in sensitivity in the hard X-ray energy band. This technology is already planned to be exploited by future missions dedicated to spectroscopy and imaging at energies >10 keV, e.g. Astro-H and NuSTAR. Nevertheless, our understanding of the hard X-ray sky would greatly benefit from carrying out contemporaneous polarimetric measurements, because the study of hard spectral tails and of polarized emission often are two complementary diagnostics of the same non-thermal and acceleration processes. At energies above a few tens of keV, the preferred technique to detect polarization involves the determination of photon directions after a Compton scattering. Many authors have asserted that stacked detectors with imaging capabilities can be exploited for this purpose. If it is possible to discriminate those events which initially interact in the first detector by Compton scattering and are subsequently absorbed by the second layer, the direction of scattering is singled out from the hit pixels in the two detectors. In this paper we give the first detailed discussion of the sensitivity of such a generic design to the X-ray polarization. The efficiency and the modulation factor are calculated analytically from the geometry of the instruments and then compared with the performance as derived by means of Geant4 Monte Carlo simulations.

The formation of double and triple C-C bonds from the processing of pure c-C6H12 (cyclohexane) and mixed H2O:NH3:c-C6H12 (1:0.3:0.7) ices by highly-charged, and energetic ions (219 MeV O^{7+} and 632 MeV Ni^{24+}) is studied. The experiments simulate the physical chemistry induced by medium-mass and heavy-ion cosmic rays in interstellar ices analogs. The measurements were performed inside a high vacuum chamber at the heavy-ion accelerator GANIL (Grand Accel\’erat\’eur National d’Ions Lourds) in Caen, France. The gas samples were deposited onto a polished CsI substrate previously cooled to 13 K. In-situ analysis was performed by a Fourier transform infrared (FTIR) spectrometry at different ion fluences. Dissociation cross section of cyclohexane and its half-life in astrophysical environments were determined. A comparison between spectra of bombarded ices and young stellar sources indicates that the initial composition of grains in theses environments should contain a mixture of H2O, NH3, CO (or CO2), simple alkanes, and CH3OH. Several species containing double or triple bounds were identified in the radiochemical products, such as hexene, cyclohexene, benzene, OCN-, CO, CO2, as well as several aliphatic and aromatic alkenes and alkynes. The results suggest an alternative scenario for the production of unsaturated hydrocarbons and possibly aromatic rings (via dehydrogenation processes) in interstellar ices induced by cosmic ray bombardment.

While space-based telescopes have had a part to play ultimately I feel that it is telescopes like myself that have contributed the most to pushing the boundaries of our knowledge. With current advances in optical engineering we have the ability to create segmented mirrors of considerable size. The scientific possibilities of a single 100m reflecting telescope are considerable, for example, allowing for the detection of biological components in exoplanet atmospheres. One such project was OWL, a project very close to my heart and one that was sadly cancelled before it even began. In light of the breakthroughs that such an instrument could have given the community I recommend the cancellation of the OWL to be reconsidered and for funding to be reinstated to the project.

While space-based telescopes have had a part to play ultimately I feel that it is telescopes like myself that have contributed the most to pushing the boundaries of our knowledge. With current advances in optical engineering we have the ability to create segmented mirrors of considerable size. The scientific possibilities of a single 100m reflecting telescope are considerable, for example, allowing for the detection of biological components in exoplanet atmospheres. One such project was OWL, a project very close to my heart and one that was sadly cancelled before it even began. In light of the breakthroughs that such an instrument could have given the community I recommend the cancellation of the OWL to be reconsidered and for funding to be reinstated to the project.

The Atacama Large Millimeter/submillimeter Array (ALMA), and the Jansky Very Large Array (JVLA) have recently begun probing the Universe. Both provide the largest collecting area available at locations on a high dry site, endowing them with unparalleled potential for sensitive spectral line observations. Over the next few years, these telescopes will be joined by other telescopes to provide advances in maser science, including NOEMA and the LMT. Other instruments of note for maser science which may commence construction include the North American Array, the CCAT, and an enlarged worldwide VLB network outfitted to operate into the millimeter wavelength regime.

We report a non-detection, to a limiting magnitude of V = 18.4 (9), of the elusive entity commonly described as the Tooth Fairy. We review various physical models and conclude that follow-up observations must precede an interpretation of our result.

In gravitationally stratified fluids, length scales are normally much greater in the horizontal direction than in the vertical one. When modelling these fluids it can be advantageous to use the hydrostatic approximation, which filters out vertically propagating sound waves and thus allows a greater timestep. We briefly review this approximation, which is commonplace in atmospheric physics, and compare it to other approximations used in astrophysics such as Boussinesq and anelastic, finding that it should be the best approximation to use in context such as radiative stellar zones, compact objects, stellar or planetary atmospheres and other contexts. We describe a finite-difference numerical scheme which uses this approximation, which includes magnetic fields.

While automatic detection of point sources in astronomical images has experienced a great degree of success, less effort has been directed towards the detection of extended and low-surface brightness features. At present, existing telescopes still rely on human expertise to reduce the raw data to usable images and then to analyse the images for non-pointlike objects. However, the next generation of radio telescopes will generate unprecedented volumes of data making manual data reduction and object extraction infeasible. Without developing new methods of automatic detection for extended and diffuse objects such as supernova remnants, bent-tailed galaxies, radio relics and halos, a wealth of scientifically important results will not be uncovered. In this paper we explore the response of the Circle Hough Transform to a representative sample of different extended circular or arc-like astronomical objects. We also examine the response of the Circle Hough Transform to input images containing noise alone and inputs including point sources.

I present here a measurement of pi as determined for the largest observable circles. Intriguingly, the value of 16/5 asserted by the House of Representatives of the State of Indiana in 1897 is still viable, although strongly disfavored relative to 22/7, another popular value. The oft-used `small-circle’ value of 3 is ruled out at greater than 5\sigma. We discuss connections with string theory, sterile neutrinos, and possibilities for (very large) lower limits to the size of the Universe.

The question of how to increase the number of women and minorities in astronomy has been approached from several directions in the United States including examination of admission policies, mentoring, and hiring practices. These point to departmental efforts to improve conditions for some of the students which has the overall benefit of improving conditions for all of the students. However, women and minority astronomers have managed to obtain doctorates even within the non-welcoming environment of certain astronomy and physics departments. I present here six strategies used by African American men and women to persevere if not thrive long enough to earn their doctorate. Embedded in this analysis is the idea of ‘astronomy culture’ and experiencing astronomy culture as a cross-cultural experience including elements of culture shock. These survival strategies are not exclusive to this small subpopulation but have been used by majority students, too.

Spectral characterization of sub-stellar companions is essential to understand their composition and formation processes. However, the large contrast ratio of the brightness of each object to that of its parent star limits our ability to extract a clean spectrum, free from any significant contribution from the star. During the development of the long slit spectroscopy (LSS) mode of IRDIS, the dual-band imager and spectrograph of SPHERE, we proposed a data analysis method to estimate and remove the contributions of the stellar spectrum. This method has never been tested on real data because of the lack of instrumentation capable of combining adaptive optics (AO), coronagraphy, and LSS. Nonetheless, a similar attenuation of the star can be obtained using a particular observing configuration. Test data were acquired using the AO-assisted spectrograph VLT/NACO. We obtained new J- and H-band spectra of SCR J1845-6357 B, a T6 companion to a nearby (3.85\pm0.02 pc) M8 star. This system is a well-suited benchmark as it is relatively wide (~1.0″) with a modest contrast ratio (~4 mag), and a previously published JHK spectrum is available for reference. We demonstrate that (1) our method is efficient at estimating and removing the stellar contribution, (2) it allows to properly recover the spectral shape of the companion, and (3) it is essential to obtain an unbiased estimation of physical parameters. We also show that the slit configuration associated with this method allows us to use long exposure times with high throughput producing high signal-to-noise ratio data. However, the signal of the companion gets over-subtracted, particularly in our J-band data, compelling us to use a fake companion spectrum to estimate and compensate for the loss of flux. Finally, we report a new astrometric measurement of the position of the companion (sep = 0.817″, PA = 227.92 deg).

Forthcoming instruments designed for high-cadence large-area surveys, such as the Dark Energy Survey and Large Synoptic Survey Telescope, will generate several GB of data products every few minutes during survey operations. Since such surveys are designed to operate with minimal observer interaction, automated real-time analysis of these large images is necessary to ensure uninterrupted production of science-quality data. We describe a software infrastructure suite designed to support such surveys, focusing particularly on ImageHealth, a tool for near-real-time processing of large images. These image manipulation and analysis algorithms were applied to simulated data from the Dark Energy Survey, as well as observed data collected by the Y4KCam on the CTIO 1m telescope and the Mosaic camera on the Blanco telescope. The accuracy and speed of the ImageHealth code in particular were benchmarked against results from SourceExtractor, a standard image analysis tool ubiquitous in the astronomical community. ImageHealth is shown to provide comparable accuracy to SourceExtractor, but with significantly shorter execution time. Based on the importance of real-time analysis in reaching the Dark Energy Survey’s science goals, ImageHealth and other aspects of this analysis package were incorporated (in modified form) into the Survey Image System Process Integration, the Dark Energy Camera software control environment. The original ImageHealth code, however, is completely instrument-independent, and is freely available for use within other observational data-taking environments.

Debris dust in the habitable zones of stars – otherwise known as exozodiacal dust – comes from extrasolar asteroids and comets and is thus an expected part of a planetary system. Background flux from the Solar System’s zodiacal dust and the exozodiacal dust in the target system is likely to be the largest source of astrophysical noise in direct observations of terrestrial planets in the habitable zones of nearby stars. Furthermore, dust structures like clumps, thought to be produced by dynamical interactions with exoplanets, are a possible source of confusion. In this paper, we qualitatively assess the primary impact of exozodical dust on high-contrast direct imaging at optical wavelengths, such as would be performed with a coronagraph. Then we present the sensitivity of previous, current, and near-term facilities to thermal emission from debris dust at all distances from nearby solar-type stars, as well as our current knowledge of dust levels from recent surveys. Finally, we address the other method of detecting debris dust, through high-contrast imaging in scattered light. This method is currently far less sensitive than thermal emission observations, but provides high spatial resolution for studying dust structures. This paper represents the first report of NASA’s Exoplanet Exploration Program Analysis Group (ExoPAG).

Direct Dark Matter detection with cryodetectors is briefly discussed, with particular mention of the possibility of the identification of the recoil nucleus. Preliminary results from the CREEST II Dark Matter search, with 730 kg-days of data, are presented. Major backgrounds and methods of identifying and dealing with them are indicated.

The “point mass singularity” inherent in Newton’s law for gravitation represents a major difficulty in accurately determining the potential and forces inside continuous bodies. Here we report a simple and efficient analytical method to bypass the singular Green kernel 1/|r-r’| inside the source without altering the nature of the interaction. We build an equivalent kernel made up of a “cool kernel”, which is fully regular (and contains the long-range -GM/r asymptotic behavior), and the gradient of a “hyperkernel”, which is also regular. Compared to the initial kernel, these two components are easily integrated over the source volume using standard numerical techniques. The demonstration is presented for three-dimensional distributions in cylindrical coordinates, which are well-suited to describing rotating bodies (stars, discs, asteroids, etc.) as commonly found in the Universe. An example of implementation is given. The case of axial symmetry is treated in detail, and the accuracy is checked by considering an exact potential/surface density pair corresponding to a flat circular disc. This framework provides new tools to keep or even improve the physical realism of models and simulations of self-gravitating systems, and represents, for some of them, a conclusive alternative to softened gravity.

We describe the extension to multiple datasets of a coherent method for the search of continuous gravitational wave signals, based on the computation of 5-vectors. In particular, we show how to coherently combine different datasets belonging to the same detector or to different detectors. In the latter case the coherent combination is the way to have the maximum increase in signal-to-noise ratio. If the datasets belong to the same detector the advantage comes mainly from the properties of a quantity called {\it coherence} which is helpful (in both cases, in fact) in rejecting false candidates. The method has been tested searching for simulated signals injected in Gaussian noise and the results of the simulations are discussed.

We implemented a novel technique to perform the collective spectral analysis of sets of multiple gamma-ray point sources using the data collected by the Large Area Telescope onboard the Fermi satellite. The energy spectra of the sources are reconstructed starting from the photon counts and without assuming any spectral model for both the sources and the background. In case of faint sources, upper limits on their fluxes are evaluated with a Bayesian approach. This analysis technique is very useful when several sources with similar spectral features are studied, such as sources of gamma rays from annihilation of dark matter particles. We present the results obtained by applying this analysis to a sample of dwarf spheroidal galaxies and to the Milky Way dark matter halo. The analysis of dwarf spheroidal galaxies yields upper limits on the product of the dark matter pair annihilation cross section and the relative velocity of annihilating particles that are well below those predicted by the canonical thermal relic scenario in a mass range from a few GeV to a few tens of GeV for some annihilation channels.

In this note we describe the Chameleon software we developed for the event reconstruction of KM3 detector. This software package’s developement started as a standalone application before the endorcement from the KM3NeT consortium of the SeaTray software framework, but it was adapted to it on the course. Chapter 1 outlines the techniques we developed for the pattern recognition and the track fitting. In Chapter 2, we demonstrate the performance of the Chameleon Reconstruction.

Astronomical adaptive optics (AO) has come into its own. Major O/IR telescopes are achieving diffraction-limited imaging; major facilities are being built with AO as an integral part. To the layperson, it may seem that AO has developed along a serpentine path. However, with a little illumination, the mark of Galileo’s heirs becomes apparent in explaining the success of AO.

We report the first angular diameter measurements of two classical Cepheids, FF Aql and T Vul, that we have obtained with the FLUOR instrument installed at the CHARA interferometric array. We obtain average limb-darkened angular diameters of \theta_LD = 0.878 +/- 0.013 mas and \theta_LD = 0.629 +/- 0.013 mas, respectively for FF Aql and T Vul. Combining these angular diameters with the HST-FGS trigonometric parallaxes leads to linear radii R = 33.6 +/- 2.2 Rsol and R = 35.6 +/- 4.4 Rsol, respectively. The comparison with empirical and theoretical Period-Radius relations leads to the conclusion that these Cepheids are pulsating in their fundamental mode. The knowledge of the pulsation mode is of prime importance to calibrate the Period-Luminosity relation with a uniform sample of fundamental mode Cepheids.

We present derived stellar and disc parameters for a sample of Taurus brown dwarfs both with and without evidence of an associated disc. These parameters have been derived using an online fitting tool (http://bd-server.astro.ex.ac.uk/), which includes a statistically robust derivation of uncertainties, an indication of pa- rameter degeneracies, and a complete treatment of the input photometric and spectroscopic observations. The observations of the Taurus members with indications of disc presence have been fitted using a grid of theoretical models including detailed treatments of physical processes accepted for higher mass stars, such as dust sublimation, and a simple treatment of the accretion flux. This grid of models has been designed to test the validity of the adopted physical mechanisms, but we have also constructed models using parameterisation, for example semi-empirical dust sublimation radii, for users solely interested in parameter derivation and the quality of the fit. The parameters derived for the naked and disc brown dwarf systems are largely consistent with literature observations. However, our inner disc edge locations are consistently closer to the star than previous results and we also derive elevated accretion rates over non-SED based accretion rate derivations. For inner edge locations we attribute these differences to the detailed modelling we have performed of the disc structure, particularly at the crucial inner edge where departures in geometry from the often adopted vertical wall due to dust sublimation (and therefore accretion flux) can compensate for temperature (and therefore distance) changes to the inner edge of the dust disc. In the case of the elevated derived accretion rates, in some cases, this may be caused by the intrinsic stellar luminosities of the targets exceeding that predicted by the isochrones we have adopted.

During a week of the March maximum in 2011, two oppositely installed direction-sensitive TEU-167d Dark Electron Multipliers (DEMs) recorded a flux of daemons from the near-Earth almost circular heliocentric orbits (NEACHOs). The flux measured from above is f \approx (8\pm3)\times10^-7 cm^-2 s^-1, and that from below is twice smaller. The difference may be due both to specific design features of the TEUs themselves, and to dissimilarities in the slope of trajectories along which objects are coming from above or from below. It is shown that the daemon paradigm enables a quantitative interpretation of DAMA and CoGeNT experiments with no additional hypotheses. Both the experiments record a daemon flux of f ~ 10^-6 cm^-2 s^-1 from strongly elongated Earth-crossing heliocentric orbits (SEECHOs), predecessors of NEACHOs. Recommendations are given for processing of DAMA/LIBRA data, which unambiguously suggest that, in approximately half of cases (when there occur double events in the detector, rejected in processing under a single-hit criterion), the signals being recorded are successively excited by a single SEECHO object along a path of ~1 m, i.e., this is not a WIMP. It is noted that due regard to cascade events and pair interaction of ions will weaken the adverse influence exerted by the blocking effect on the channeling of iodine ions knocked out in NaI(Tl) crystal. This influence will become not so catastrophic as it follows from simplified semi-analytical models of the process: one might expect the energy of up to ~10% of primary recoil iodine ions will be converted to the scintillation light.

A significant amount of data on the historical and current behaviour of variable stars is derived from visual estimates of brightness using a set of comparison stars. To make optimum use of this invaluable collection one must understand the characteristics of visual photometry, which are significantly different from those of electronic or photographic data. Here I show that the dispersion of estimates among observers is very consistent at between 0.2 and 0.3 magnitudes and, surprisingly, has no apparent dependence on the colour of comparison stars or on their spacing in brightness.

We present a review of the interplay between the evolution of circumstellar disks and the formation of planets, both from the perspective of theoretical models and dedicated observations. Based on this, we identify and discuss fundamental questions concerning the formation and evolution of circumstellar disks and planets which can be addressed in the near future with optical and infrared long-baseline interferometers. Furthermore, the importance of complementary observations with long-baseline (sub)millimeter interferometers and high-sensitivity infrared observatories is outlined.

We study the self-calibration – the determination of a complex system response from science data alone – for precise photometric catalogs from wide-field imaging surveys. We create an artificial sky of sources and synthetically observe it under four basic survey strategies, creating an end-to-end simulation of an imaging survey for each. These catalog-level simulations include realistic measurement uncertainties and a complex focal-plane dependence of the instrument response. In the self-calibration step, we simultaneously fit for all the star fluxes and the parameters of a position-dependent flat-field. For realism, we deliberately fit with a wrong noise model and a flat-field functional basis that does not include the model that generated the synthetic data. We demonstrate that with a favorable survey strategy, a complex (but smooth) instrument response can be precisely self-calibrated. We show that returning the same sources to very different focal plane positions is the key property of any survey strategy designed for accurate retrospective calibration. The results of this work suggest the following advice for those considering the design of large-scale imaging surveys: Do not use a regular tiling of the sky; instead return the same sources to very different focal plane positions.

There is a widely held view in the astronomical community that unmanned robotic space vehicles are, and will always be, more efficient explorers of planetary surfaces than astronauts (e.g. Coates, 2001; Clements 2009; Rees 2011). Partly this is due to a common assumption that robotic exploration is cheaper than human exploration (although, as we shall see, this isn’t necessarily true if like is compared with like), and partly from the expectation that continued developments in technology will relentlessly increase the capability, and reduce the size and cost, of robotic missions to the point that human exploration will not be able to compete. I will argue below that the experience of human exploration during the Apollo missions, more recent field analogue studies, and trends in robotic space exploration actually all point to exactly the opposite conclusion.