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

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

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

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

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

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

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

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

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

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

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

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

We present details of improvements to data processing and analysis which were recently used for a re-reduction of the Very Large Array (VLA) Low-frequency Sky Survey (VLSS) data. Algorithms described are implemented in the data-reduction package Obit, and include smart-windowing to reduce clean bias, improved automatic radio frequency interference removal, improved bright-source peeling, and higher-order Zernike fits to model the ionospheric phase contributions. An additional, but less technical improvement was using the original VLSS catalog as a same-frequency/same-resolution reference for calculating ionospheric corrections, allowing more accuracy and a higher percentage of data for which solutions are found. We also discuss new algorithms for extracting a source catalog and analyzing ionospheric fluctuations present in the data. The improved reduction techniques led to substantial improvements including images of six previously unpublished fields (1% of the survey area) and reducing the clean bias by 50%. The largest angular size imaged has been roughly doubled, and the number of cataloged sources is increased by 35% to 95,000.

The Durham adaptive optics real-time controller is a generic, high performance real-time control system for astronomical adaptive optics systems. It has recently had new features added as well as performance improvements, and here we give details of these, as well as ways in which optimisations can be made for specific adaptive optics systems and hardware implementations. We also present new measurements that show how this real-time control system could be used with any existing adaptive optics system, and also show that when used with modern hardware, it has high enough performance to be used with most Extremely Large Telescope adaptive optics systems.

In this pedagogical text aimed at those wanting to start thinking about or brush up on probabilistic inference, I review the rules by which probability distribution functions can (and cannot) be combined. I connect these rules to the operations performed in probabilistic data analysis. Dimensional analysis is emphasized as a valuable tool for helping to construct non-wrong probabilistic statements. The applications of probability calculus in constructing likelihoods, marginalized likelihoods, posterior probabilities, and posterior predictions are all discussed.

In magnetohydrodynamics (MHD), the magnetic field is evolved by the induction equation and coupled to the gas dynamics by the Lorentz force. We perform numerical smoothed particle magnetohydrodynamics (Spmhd) simulations and study the influence of a numerical magnetic divergence. For instabilities arising from divergence B related errors, we find the hyperbolic/parabolic cleaning scheme suggested by Dedner et al. 2002 to give good results and prevent numerical artifacts from growing. Additionally, we demonstrate that certain current Spmhd implementations of magnetic field regularizations give rise to unphysical instabilities in long-time simulations. We also find this effect when employing Euler potentials (divergenceless by definition), which are not able to follow the winding-up process of magnetic field lines properly. Furthermore, we present cosmological simulations of galaxy cluster formation at extremely high resolution including the evolution of magnetic fields. We show synthetic Faraday rotation maps and derive structure functions to compare them with observations. Comparing all the simulations with and without divergence cleaning, we are able to confirm the results of previous simulations performed with the standard implementation of MHD in Spmhd at normal resolution. However, at extremely high resolution, a cleaning scheme is needed to prevent the growth of numerical errors at small scales.

We present the first attempts to build a user-friendly interface for the Virtual Observatory of the University of Guanajuato. The data tables will be accessible to the public through PHP scripts and SQL database managers, such as MySQL and PostgreSQL, all administrated through phpMyAdmin and pgMyAdmin. Although it is not made public yet, this interface will be the basis upon which the final front end for our VO will be built. Furthermore, we present a preliminary version of a web front end to the publicly available stellar population synthesis code STARLIGHT (starlight.ufsc.br) which will be made available with our VO. This front end aims to provide an easy and flexible access to the code itself, letting users fit their own observed spectra with their preferred combination of physical and technical parameters, rather than making available only the results of fitting a specific sample of spectra with predefined parameters.

Radio interferometers can measure the full polarization state of incoming waves by cross–correlating all combinations of two orthogonal polarizations at each antenna. The independent sets of electronics used to detect the two polarization states will introduce a differential instrumental delay between the two data streams. The usual technique of separate calibration of the parallel–hand sets of visibilities still allows for an arbitrary offset in group delay and phase between the two parallel systems. In order to use the cross–polarized visibilities, this instrumental offset must be determined and removed. This paper describes one such technique and explores its application in the Obit package. The technique is successfully applied to some EVLA data using both strongly and weakly polarized calibrators.

Reflective light collectors with hexagonal entrance and exit apertures are frequently used in front of the focal-plane camera of a very-high-energy gamma-ray telescope to increase the collection efficiency of atmospheric Cherenkov photons and reduce the night-sky background entering at large incident angles. The shape of a hexagonal light collector is usually based on Winston’s design, which is optimized for only two-dimensional optical systems. However, it is not known whether a hexagonal Winston cone is optimal for the real three-dimensional optical systems of gamma-ray telescopes. For the first time we optimize the shape of a hexagonal light collector using quadratic and cubic B\’ezier curves. We demonstrate that our optimized designs simultaneously achieve a higher collection efficiency and background reduction rate than traditional designs.

Context. The INTEGRAL observatory operating in a hard X-ray/gamma domain gathered a large observational data set over nine years since 2003. Dominant fraction of the observing time was dedicated to the Galactic source population study making the possibility of the deepest Galactic survey in hard X-rays ever compiled. Aims. The aim of the current Galactic survey is to make a basis for Galactic source population studies, and perform mapping of the Milky Way in hard X-rays over the maximum exposure available at |b|<17.5 deg. Methods. We used sky reconstruction algorithms specially developed for the high quality imaging of INTEGRAL/IBIS data. Results. We present sky images, sensitivity maps, and catalogs of detected sources in three energy bands: 17-60, 17-35, and 35-80 keV in the Galactic plane at |b|5 sigma has an identification completeness of ~91%, which is valuable for population studies.

We have developed the Smoothed Bandpass Calibration (SBC) method and the best suitable scan pattern to optimize radio spectroscopic observations. Adequate spectral smoothing is applied to the spectrum toward OFF-source blank sky adjacent to a target source direction for the purpose of bandpass correction. Because the smoothing process reduces noise, the integration time for OFF-source scans can be reduced keeping the signal-to-noise ratio. Since the smoothing is not applied to ON-source scans, the spectral resolution for line features is kept. An optimal smoothing window is determined by bandpass flatness evaluated by Spectral Allan Variance (SAV). An efficient scan pattern is designed to the OFF-source scans within the bandpass stability timescale estimated by Time-based Allan Variance (TAV). We have tested the SBC using the digital spectrometer, VESPA, on the VERA Iriki station. For the targeted noise level of 5e-4 as a ratio to the system noise, the optimal smoothing window was 32 – 60 ch in the whole bandwidth of 1024 ch, and the optimal scan pattern was designed as a sequence of 70-s ON + 10-s OFF scan pairs. The noise level with the SBC was reduced by a factor of 1.74 compared with the conventional method. The total telescope time to achieve the goal with the SBC was 400 s, which was 1/3 of 1200 s required by the conventional way. Improvement in telescope time efficiency with the SBC was calculated as 3x, 2x and 1.3x for single-beam, dual-beam, and on-the-fly (OTF) scans, respectively. The SBC works to optimize scan patterns for observations from now, and also works to improve signal-to-noise ratios of archival data if ON- and OFF-source spectra are individually recorded, though the efficiency depends on the spectral stability of the receiving system.

We have developed the Smoothed Bandpass Calibration (SBC) method and the best suitable scan pattern to optimize radio spectroscopic observations. Adequate spectral smoothing is applied to the spectrum toward OFF-source blank sky adjacent to a target source direction for the purpose of bandpass correction. Because the smoothing process reduces noise, the integration time for OFF-source scans can be reduced keeping the signal-to-noise ratio. Since the smoothing is not applied to ON-source scans, the spectral resolution for line features is kept. An optimal smoothing window is determined by bandpass flatness evaluated by Spectral Allan Variance (SAV). An efficient scan pattern is designed to the OFF-source scans within the bandpass stability timescale estimated by Time-based Allan Variance (TAV). We have tested the SBC using the digital spectrometer, VESPA, on the VERA Iriki station. For the targeted noise level of 5e-4 as a ratio to the system noise, the optimal smoothing window was 32 – 60 ch in the whole bandwidth of 1024 ch, and the optimal scan pattern was designed as a sequence of 70-s ON + 10-s OFF scan pairs. The noise level with the SBC was reduced by a factor of 1.74 compared with the conventional method. The total telescope time to achieve the goal with the SBC was 400 s, which was 1/3 of 1200 s required by the conventional way. Improvement in telescope time efficiency with the SBC was calculated as 3x, 2x and 1.3x for single-beam, dual-beam, and on-the-fly (OTF) scans, respectively. The SBC works to optimize scan patterns for observations from now, and also works to improve signal-to-noise ratios of archival data if ON- and OFF-source spectra are individually recorded, though the efficiency depends on the spectral stability of the receiving system.

We have developed a compact cryogenic system with a pulse tube refrigerator and a coaxial heat exchanger. This liquefaction-purification system not only saves the cooling power used to reach high gaseous recirculation rate, but also reduces the impurity level with high speed. The heat exchanger operates with an efficiency of 99%, which indicates the possibility for fast xenon gas recirculation in a highpressurized large-scale xenon storage with much less thermal losses.

In this white paper, we discuss the concept of a next-generation X-ray mission called BEST (Black hole Evolution and Space Time). The mission concept uses a 3000 square centimeter effective area mirror (at 6 keV) to achieve unprecedented sensitivities for hard X-ray imaging spectrometry (5-70 keV) and for broadband X-ray polarimetry (2-70 keV). BEST can make substantial contributions to our understanding of the inner workings of accreting black holes, our knowledge about the fabric of extremely curved spacetime, and the evolution of supermassive black holes. BEST will allow for time resolved studies of accretion disks. With a more than seven times larger mirror area and a seven times wider bandpass than GEMS, BEST will take X-ray polarimetry to a new level: it will probe the time variability of the X-ray polarization from stellar mass and supermassive black holes, and it will measure the polarization properties in 30 independent energy bins. These capabilities will allow BEST to conduct tests of accretion disk models and the underlying spacetimes. With three times larger mirror area and ten times better angular resolution than NuSTAR, BEST will be able to make deep field observations with a more than 15 times better sensitivity than NuSTAR. The mission will be able to trace the evolution of obscured and unobscured black holes in the redshift range from zero to six, covering the most important epoch of supermassive black hole growth. The hard X-ray sensitivity of BEST will enable a deep census of non-thermal particle populations. BEST will give us insights into AGN feedback by measuring the particle luminosity injected by AGNs into the interstellar medium (ISM) of their hosts, and will map the emission from particles accelerated at large scale structure shocks. Finally, BEST has the potential to constrain the equation of state of neutron stars (NS).

The advent of stable, highly squeezed states of light has generated great interest in the gravitational wave community as a means for improving the quantumnoise- limited performance of advanced interferometric detectors. To confidently measure these squeezed states, it is first necessary to measure the shot-noise across the frequency band of interest. Technical noise, such as non-stationary events, beam pointing, and parasitic interference, can corrupt shot-noise measurements at low Fourier frequencies, below tens of kilo-Hertz. In this paper we present a qualitative investigation into all of the relevant noise sources and the methods by which they can be identified and mitigated in order to achieve quantum noise limited balanced homodyne detection. Using these techniques, flat shot-noise down to Fourier frequencies below 0.5 Hz is produced. This enables the direct observation of large magnitudes of squeezing across the entire audio-band, of particular interest for ground-based interferometric gravitational wave detectors. 11.6 dB of shot-noise suppression is directly observed, with more than 10 dB down to 10 Hz.

The structure of the solar chromosphere is believed to be governed by magnetic fields, even in quiet-Sun regions that have a relatively weak photospheric field. During the past decade inversion methods have emerged as powerful tools for analyzing the chromosphere of active regions. The applicability of inversions to infer the stratification of the physical conditions in a dynamic 3D solar chromosphere has not yet been studied in detail. This study aims to establish the diagnostic capabilities of non-local thermodynamical equilibrium (NLTE) inversion techniques of Stokes profiles induced by the Zeeman effect in the Ca II 8542 line. We computed the Ca II atomic level populations in a snapshot from a 3D radiation-MHD simulation of the quiet solar atmosphere in non-LTE using the 3D radiative transfer code Multi3d. These populations were used to compute synthetic full-Stokes profiles in the Ca II 8542 line using 1.5D radiative transfer and the inversion code Nicole. The profiles were then spectrally degraded to account for finite filter width and Gaussian noise was added to account for finite photon flux. These profiles were inverted using Nicole and the results were compared with the original model atmosphere. Our NLTE inversions applied to quiet-Sun synthetic observations provide reasonably good estimates of the chromospheric magnetic field, line-of-sight velocities and somewhat less accurate, but still very useful, estimates of the temperature. Three dimensional scattering of photons cause cool pockets in the chromosphere to be invisible in the line profile and consequently they are also not recovered by the inversions. To successfully detect Stokes linear polarization in this quiet snapshot, a noise level below 10^{-3.5} is necessary.

We propose a novel algorithm for image reconstruction in radio interferometry. The ill-posed inverse problem associated with the incomplete Fourier sampling identified by the visibility measurements, is regularized by the assumption of average signal sparsity over representations in multiple wavelet bases. The algorithm, defined in the versatile framework of convex optimization, is dubbed Sparsity Averaging Reweighted Analysis (SARA). We show through simulations that the proposed approach largely outperforms state-of-the-art imaging methods in the field, which are based on the assumption of signal sparsity in a single basis only.

The changes in phase, time and frequency suffered by signals when retransmitted by a remote and inaccessible transponder and the propagation delays are major constraints to obtain accurate ranging measurements in various related applications. We present a new method and system to determine these delays for every single pulsed signal transmission. The process utilizes four ground-based reference stations, synchronized in time and installed at well known geodesic coordinates. The repeater station is located within the fields of view common to the four reference bases, such as in a platform transported by a satellite, balloon, aircraft, etc. Signal transmitted by one of the reference bases is retransmitted by the transponder, received back by the four bases, producing four ranging measurements which are processed to determine uniquely the time delays undergone in every retransmission process. The repeater’s positions with respect to each group of three out of four reference bases are given by a system of equations. A minimization function is derived comparing repeater’s positions referred to at least two groups of three reference bases. The minimum found by iterative methods provide the signal transit time at the repeater and propagation delays, providing the correct repeater position. The method is applicable to the transponder platform positioning and navigation, time synchronization of remote clocks, and location of targets. The algorithm has been fully demonstrated simulated for practical situation with the transponder carried by an aircraft moving over bases on the ground. The errors of the determinations have been evaluated for uncertainties in clock synchronization, in propagation time delays and other system parameters.

Residual uncertainty in the calibration of gravitational wave (GW) detector data leads to systematic errors which must be accounted for in setting limits on the strength of GW signals. When cross-correlation measurements are made using data from a pair of instruments, as in searches for a stochastic GW background, the calibration uncertainties associated with the two instruments can be combined into an uncertainty associated with the pair. With the advent of multi-baseline GW observation (e.g., networks consisting of multiple detectors such as the LIGO observatories and Virgo), a more sophisticated treatment is called for. We describe how the correlations between calibration factors associated with different pairs can be taken into account by marginalizing over the uncertainty associated with each instrument, defining two methods known as per-baseline and per-instrument marginalization.

The AARTFAAC project aims to implement an All-Sky Monitor (ASM), using the Low Frequency Array (LOFAR) telescope. It will enable real-time, 24×7 monitoring for low frequency radio transients over most of the sky locally visible to the LOFAR at timescales ranging from milliseconds to several days, and rapid triggering of follow-up observations with the full LOFAR on detection of potential transient candidates. These requirements pose several implementation challenges: imaging of an all-sky field of view, low latencies of processing, continuous availability and autonomous operation of the ASM. The first of these has already resulted in the correlator for the ASM being the largest in the world in terms of its number of input channels. It will generate $\sim 1.5 \cdot 10^5$ correlations per second per spectral channel when built. Test observations using existing LOFAR infrastructure were carried out to quantify and constrain crucial instrumental design criteria for the ASM. In this paper, we present an overview of the AARTFAAC data processing pipeline and illustrate some of the aforementioned challenges by showing all-sky images obtained from one of the test observations. These results provide quantitative estimates of the capabilities of the instrument.

We present the results of a study on the impact of particles of galactic (GCR) and solar origin for the X-ray Microcalorimeter Spectrometer (XMS) aboard an astronomical satellite flying in an orbit at the second Lagrangian point (L2). The detailed configuration presented in this paper is the one adopted for the International X-Ray Observatory (IXO) study, however the derived estimates can be considered a conservative limit for ATHENA, that is the IXO redefined mission proposed to ESA. This work is aimed at the estimate of the residual background level expected on the focal plane detector during the mission lifetime, a crucial information in the development of any instrumental configuration that optimizes the XMS scientific performances. We used the Geant4 toolkit, a Monte Carlo based simulator, to investigate the rejection efficiency of the anticoincidence system and assess the residual background on the detector.

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

In self-consistent N-body simulations of collisionless systems, gravitational interactions are modified on small scales to remove singularities and simplify the task of numerically integrating the equations of motion. This `gravitational softening’ is sometimes presented as an ad-hoc departure from Newtonian gravity. However, softening can also be described as a smoothing operation applied to the mass distribution; the gravitational potential and the smoothed density obey Poisson’s equation precisely. While `softening’ and `smoothing’ are mathematically equivalent descriptions, the latter has some advantages. For example, the smoothing description suggests a way to set up N-body initial conditions in almost perfect dynamical equilibrium.

We consider the problem of inferring constraints on a high-dimensional parameter space with a computationally expensive likelihood function. Markov chain Monte Carlo (MCMC) methods offer significant improvements in efficiency over grid-based searches and are easy to implement in a wide range of cases. However, MCMCs offer few guarantees that all of the interesting regions of parameter space are explored. We propose a machine learning algorithm that improves upon the performance of MCMC by intelligently targeting likelihood evaluations so as to quickly and accurately characterize the likelihood surface in both low- and high-likelihood regions. We compare our algorithm to MCMC on toy examples and the 7-year WMAP cosmic microwave background data release. Our algorithm finds comparable parameter constraints to MCMC in fewer calls to the likelihood function and with greater certainty that all of the interesting regions of parameter space have been explored.

Scanamorphos is one of the public softwares available to post-process scan observations performed with the Herschel photometer arrays. This post-processing mainly consists in subtracting the total low-frequency noise (both its thermal and non-thermal components), masking cosmic ray hit residuals, and projecting the data onto a map. Although it was developed for Herschel, it is also applicable with minimal adjustment to scan observations made with other bolometer arrays, provided they entail sufficient redundancy; it was successfully applied to P-Artemis, an instrument operating on the APEX telescope. Contrary to most other algorithms (first developed for microwave background experiments and later adapted to Herschel), Scanamorphos does not assume any particular noise model, and does not apply any Fourier-space filtering to the data, but is an empirical tool using purely the redundancy built in the observations — taking advantage of the fact that each portion of the sky is sampled at multiple times by multiple bolometers. It is an interactive software in the sense that the user is allowed to optionally visualize and control results at each intermediate step, but the processing is fully automated. This paper describes the principles and algorithm of Scanamorphos and presents several examples of application.

The Virtual Observatory (VO) is becoming the de-facto standard for astronomical data publication. However, the number of radio astronomical archives is still low in general, and even lower is the number of radio astronomical data available through the VO. In order to facilitate the building of new radio astronomical archives, easing at the same time their interoperability with VO framework, we have developed a VO-compliant data model which provides interoperable data semantics for radio data. That model, which we call the Radio Astronomical DAta Model for Single-dish (RADAMS) has been built using standards of (and recommendations from) the International Virtual Observatory Alliance (IVOA). This article describes the RADAMS and its components, including archived entities and their relationships to VO metadata. We show that by using IVOA principles and concepts, the effort needed for both the development of the archives and their VO compatibility has been lowered, and the joint development of two radio astronomical archives have been possible. We plan to adapt RADAMS to be able to deal with interferometry data in the future.

We present Iris, the VAO (Virtual Astronomical Observatory) application for analyzing SEDs (spectral energy distributions). Iris is the result of one of the major science initiatives of the VAO, and the first version was released in September 2011. Iris combines key features of several existing software applications to streamline and enhance SED analysis. With Iris, users may read and display SEDs, select data ranges for analysis, fit models to SEDs, and calculate confidence limits on best-fit parameters. SED data may be uploaded into the application from IVOA-compliant VOTable and FITS format files, or retrieved directly from NED. Data written in unsupported formats may be converted using SedImporter, a new application provided with Iris. The components of Iris have been contributed by members of the VAO. Specview, contributed by STScI, provides a GUI for reading, editing, and displaying SEDs, as well as defining models and parameter values. Sherpa, contributed by the Chandra project at SAO, provides a library of models, fit statistics, and optimization methods; the underlying I/O library, SEDLib, is a VAO product written by SAO to current IVOA (International Virtual Observatory Alliance) data model standards. NED is a service provided by IPAC for easy location of data for a given extragalactic source, including SEDs. SedImporter is a new tool for converting non-standard SED data files into a format supported by Iris. We demonstrate the use of SedImporter to retrieve SEDs from a variety of sources–from the NED SED service, from the user’s own data, and from other VO applications using SAMP (Simple Application Messaging Protocol). We also demonstrate the use of Iris to read, display, select ranges from, and fit models to SEDs. Finally, we discuss the architecture of Iris, and the use of IVOA standards so that Specview, Sherpa, SEDLib and SedImporter work together seamlessly.

We analyze the variability of the spectral solar irradiance during the period from 7 January, 2010 until 20 January, 2010 as measured by the Herzberg channel (190-222 nm) of the Large Yield RAdiometer (LYRA) onboard PROBA2. In this period of time observations by the LYRA nominal unit experienced degradation and the signal produced by the Herzberg channel frequently jumped from one level to another. Both these factors significantly complicates the analysis. We present the algorithm which allowed us to extract the solar variability from the LYRA data and compare the results with SORCE/SOLSTICE measurements and with modeling based on the Code for the Solar Irradiance (COSI).

The time-symmetric block time–step (TSBTS) algorithm is a newly developed efficient scheme for $N$–body integrations. It is constructed on an era-based iteration. In this work, we re-designed the TSBTS integration scheme with dynamically changing era size. A number of numerical tests were performed to show the importance of choosing the size of the era, especially for long time integrations. Our second aim was to show that the TSBTS scheme is as suitable as previously known schemes for developing parallel $N$–body codes. In this work, we relied on a parallel scheme using the copy algorithm for the time-symmetric scheme. We implemented a hybrid of data and task parallelization for force calculation to handle load balancing problems that can appear in practice. Using the Plummer model initial conditions for different numbers of particles, we obtained the expected efficiency and speedup for a small number of particles. Although parallelization of the direct $N$–body codes is negatively affected by the communication/calculation ratios, we obtained good load balance results. Moreover, we were able to conserve the advantages of the algorithm (e.g., energy conservation for long–term simulations).

Results are presented from radioactivity screening of two models of photomultiplier tubes designed for use in current and future liquid xenon experiments. The Hamamatsu 5.6 cm diameter R8778 PMT, used in the LUX dark matter experiment, has yielded a positive detection of four common radioactive isotopes: 238U, 232Th, 40K, and 60Co. Screening of LUX materials has rendered backgrounds from other detector materials subdominant to the R8778 contribution. A prototype Hamamatsu 7.6 cm diameter R11410 MOD PMT has also been screened, with benchmark isotope counts measured at <0.4 238 U / <0.3 232 Th / <8.3 40 K / 2.0+-0.2 60 Co mBq/PMT. This represents a large reduction, equal to a change of \times 1/24 238U / \times 1/9 232Th / \times 1/8 40K per PMT, between R8778 and R11410 MOD, concurrent with a doubling of the photocathode surface area (4.5 cm to 6.4 cm diameter). 60Co measurements are comparable between the PMTs, but can be significantly reduced in future R11410 MOD units through further material selection. Assuming PMT activity equal to the measured 90% upper limits, Monte Carlo estimates indicate that replacement of R8778 PMTs with R11410 MOD PMTs will change LUX PMT electron recoil background contributions by a factor of \times1/25 after further material selection for 60Co reduction, and nuclear recoil backgrounds by a factor of \times 1/36. The strong reduction in backgrounds below the measured R8778 levels makes the R11410 MOD a very competitive technology for use in large-scale liquid xenon detectors.

On the final day of the Stellar Polarimetry conference, participants split up into three “breakout sessions” to discuss the future of the field in the areas of instrumentation, upcoming opportunities, and community priorities. This contribution compiles the major recommendations arising from each breakout session. We hope that the polarimetric community will find these ideas useful as we consider how to maintain the vitality of polarimetry in the coming years.

The Auger Engineering Radio Array (AERA) is currently detecting cosmic rays of energies at and above 10^17 eV at the Pierre Auger Observatory, by triggering on the radio emission produced in the associated air showers. The radio-detection technique must cope with a significant background of man-made radio-frequency interference, but can provide information on shower development with a high duty cycle. We discuss our techniques to handle the challenges of self-triggered radio detection in a low-power autonomous array, including triggering and filtering algorithms, data acquisition design, and communication systems.

This review outlines concepts of mathematical statistics, elements of probability theory, hypothesis tests and point estimation for use in the analysis of modern astronomical data. Least squares, maximum likelihood, and Bayesian approaches to statistical inference are treated. Resampling methods, particularly the bootstrap, provide valuable procedures when distributions functions of statistics are not known. Several approaches to model selection and good- ness of fit are considered. Applied statistics relevant to astronomical research are briefly discussed: nonparametric methods for use when little is known about the behavior of the astronomical populations or processes; data smoothing with kernel density estimation and nonparametric regression; unsupervised clustering and supervised classification procedures for multivariate problems; survival analysis for astronomical datasets with nondetections; time- and frequency-domain times series analysis for light curves; and spatial statistics to interpret the spatial distributions of points in low dimensions. Two types of resources are presented: about 40 recommended texts and monographs in various fields of statistics, and the public domain R software system for statistical analysis. Together with its \sim 3500 (and growing) add-on CRAN packages, R implements a vast range of statistical procedures in a coherent high-level language with advanced graphics.

The Earth’s atmosphere affects ground-based astronomical observations. Scattering, absorption, and radiation processes deteriorate the signal-to-noise ratio of the data received. For scheduling astronomical observations it is, therefore, important to accurately estimate the wavelength-dependent effect of the Earth’s atmosphere on the observed flux. In order to increase the accuracy of the exposure time calculator of the European Southern Observatory’s (ESO) Very Large Telescope (VLT) at Cerro Paranal, an atmospheric model was developed as part of the Austrian ESO In-Kind contribution. It includes all relevant components, such as scattered moonlight, scattered starlight, zodiacal light, atmospheric thermal radiation and absorption, and non-thermal airglow emission. This paper focuses on atmospheric scattering processes that mostly affect the blue ( 0.55 mum) wavelength regime. While the former is mainly investigated by means of radiative transfer models, the intensity and variability of the latter is studied with a sample of 1186 VLT FORS1 spectra. For a set of parameters such as the object altitude angle, Moon-object angular distance, ecliptic latitude, bimonthly period, and solar radio flux, our model predicts atmospheric radiation and transmission at a requested resolution. A comparison of our model with the FORS1 spectra and photometric data for the night-sky brightness from the literature, suggest a model accuracy of about 20%. This is a significant improvement with respect to existing predictive atmospheric models for astronomical exposure time calculators.

The CoRoT field LRa02 has been observed with the Berlin Exoplanet Search Telescope II (BEST II) during the southern summer 2007/2008. A first analysis of stellar variability led to the publication of 345 newly discovered variable stars. Now, a deeper analysis of this data set was used to optimize the variability search procedure. Several methods and parameters have been tested in order to improve the selection process compared to the widely used J index for variability ranking. This paper describes an empirical approach to treat systematic trends in photometric data based upon the analysis of variance statistics that can significantly decrease the rate of false detections. Finally, the process of reanalysis and method improvement has virtually doubled the number of variable stars compared to the first analysis by Kabath et al. A supplementary catalog of 272 previously unknown periodic variables plus 52 stars with suspected variability is presented. Improved ephemerides are given for 19 known variables in the field. In addition, the BEST II results are compared with CoRoT data and its automatic variability classification.

The CoRoT field LRa02 has been observed with the Berlin Exoplanet Search Telescope II (BEST II) during the southern summer 2007/2008. A first analysis of stellar variability led to the publication of 345 newly discovered variable stars. Now, a deeper analysis of this data set was used to optimize the variability search procedure. Several methods and parameters have been tested in order to improve the selection process compared to the widely used J index for variability ranking. This paper describes an empirical approach to treat systematic trends in photometric data based upon the analysis of variance statistics that can significantly decrease the rate of false detections. Finally, the process of reanalysis and method improvement has virtually doubled the number of variable stars compared to the first analysis by Kabath et al. A supplementary catalog of 272 previously unknown periodic variables plus 52 stars with suspected variability is presented. Improved ephemerides are given for 19 known variables in the field. In addition, the BEST II results are compared with CoRoT data and its automatic variability classification.

This paper presents a conceptual design for a spectrometer designed specifically for characterizing transiting exoplanets with space-borne infrared telescopes. The design adopting cross-dispersion is intended to be simple, compact, highly stable, and has capability of simultaneous coverage over a wide wavelength region with high throughput. Typical wavelength coverage and spectral resolving power is 1-13 micron with a spectral resolving power of ~ a few hundred, respectively. The baseline design consists of two detectors, two prisms with a dichroic coating and microstructured grating surfaces, and three mirrors. Moving parts are not adopted. The effect of defocusing is evaluated for the case of a simple shift of the detector, and anisotropic defocusing to maintain the spectral resolving power. Variations in the design and its application to planned missions are also discussed.

Space-based gravitational-wave observatories such as the Laser Interferometer Space Antenna (LISA) use time-shifted and time-scaled linear combinations of differential laser-phase beat signals to cancel the otherwise overwhelming laser frequency noise. Nanosecond timing precision is needed to accurately form these Time-Delay Interferometry (TDI) combinations which defines a ~1 meter requirement on the inter-spacecraft ranging capability. The University of Florida Hardware-in-the-loop LISA Interferometry Simulator (UFLIS) has been used to test Time-Delay Interferometry in a configuration which incorporates variable delays, realistic Doppler shifts, and simulated gravitational-wave signals. The TDI 2.0 combinations are exploited to determine the time-changing delays with nanosecond accuracy using a TDI-ranging reference tone. These variable delays are used in forming the TDI combinations to achieve the LISA interferometry sensitivity resulting from 10 orders of magnitude laser frequency noise cancellation.

Nano and micro meter sized dust particles travelling through the heliosphere at several hundreds of km/s have been repeatedly detected by interplanetary spacecraft. When such fast moving dust particles hit a solid target in space, an expanding plasma cloud is formed through the vaporisation and ionisation of the dust particles itself and part of the target material at and near the impact point. Immediately after the impact the small and dense cloud is dominated by collisions and the expansion can be described by fluid equations. However, once the cloud has reached micro-m dimensions, the plasma may turn collisionless and a kinetic description is required to describe the subsequent expansion. In this paper we explore the late and possibly collisionless spherically symmetric unconstrained expansion of a single ionized ion-electron plasma using N-body simulations. Given the strong uncertainties concerning the early hydrodynamic expansion, we assume that at the time of the transition to the collisionless regime the cloud density and temperature are spatially uniform. We do also neglect the role of the ambient plasma. This is a reasonable assumption as long as the cloud density is substantially higher than the ambient plasma density. In the case of clouds generated by fast interplanetary dust grains hitting a solid target some 10^7 electrons and ions are liberated and the in vacuum approximation is acceptable up to meter order cloud dimensions. …

The Cherenkov Telescope Array (CTA), currently in its Preparatory Phase, will be the first open observatory for very high energy gamma-rays from galactic and extragalactic sources. The international consortium behind CTA is preparing the construction of two large arrays of Cherenkov telescopes in the Northern and Southern Hemispheres with a performance that will be significantly improved compared to the current generation of arrays. Its increased sensitivity and energy range will give CTA access to a large population of Active Galactic Nuclei (AGN) not yet detected at very high energies and provide much more details on known TeV sources. While the low end of the CTA energy coverage will close the current gap with the Fermi-LAT band, its high energy coverage will open a new window on the sky and help us understand the intrinsic shape of the hardest blazar spectra. We outline the current status of CTA and discuss the science case for AGN physics with the observatory. Predictions for source detections based on extrapolations of Fermi-LAT spectra are discussed. An overview is given of prospects for the detection of extended emission from radio galaxies, of rapid variability, and spectral features. The observation of AGN with CTA will also improve current constraints on the distribution of the extragalactic background light, the strength of the intergalactic magnetic field and Lorentz invariance violation.

The Cherenkov Telescope Array (CTA), currently in its Preparatory Phase, will be the first open observatory for very high energy gamma-rays from galactic and extragalactic sources. The international consortium behind CTA is preparing the construction of two large arrays of Cherenkov telescopes in the Northern and Southern Hemispheres with a performance that will be significantly improved compared to the current generation of arrays. Its increased sensitivity and energy range will give CTA access to a large population of Active Galactic Nuclei (AGN) not yet detected at very high energies and provide much more details on known TeV sources. While the low end of the CTA energy coverage will close the current gap with the Fermi-LAT band, its high energy coverage will open a new window on the sky and help us understand the intrinsic shape of the hardest blazar spectra. We outline the current status of CTA and discuss the science case for AGN physics with the observatory. Predictions for source detections based on extrapolations of Fermi-LAT spectra are discussed. An overview is given of prospects for the detection of extended emission from radio galaxies, of rapid variability, and spectral features. The observation of AGN with CTA will also improve current constraints on the distribution of the extragalactic background light, the strength of the intergalactic magnetic field and Lorentz invariance violation.

High rate sampling detectors measuring the potential difference between the main body and boom antennas of interplanetary spacecraft have been shown to be efficient means to measure the voltage pulses induced by nano dust impacts on the spacecraft body itself (see Meyer-Vernet et al, Solar Phys. 256, 463 (2009)). However, rough estimates of the free charge liberated in post impact expanding plasma cloud indicate that the cloud’s own internal electrostatic field is too weak to account for measured pulses as the ones from the TDS instrument on the STEREO spacecraft frequently exceeding 0.1 V/m. In this paper we argue that the detected pulses are not a direct measure of the potential structure of the plasma cloud, but are rather the consequence of a transitional interruption of the photoelectron return current towards the portion of the antenna located within the expanding cloud.

High rate sampling detectors measuring the potential difference between the main body and boom antennas of interplanetary spacecraft have been shown to be efficient means to measure the voltage pulses induced by nano dust impacts on the spacecraft body itself (see Meyer-Vernet et al, Solar Phys. 256, 463 (2009)). However, rough estimates of the free charge liberated in post impact expanding plasma cloud indicate that the cloud’s own internal electrostatic field is too weak to account for measured pulses as the ones from the TDS instrument on the STEREO spacecraft frequently exceeding 0.1 V/m$. In this paper we argue that the detected pulses are not a direct measure of the potential structure of the plasma cloud, but are rather the consequence of a transitional interruption of the photoelectron return current towards the portion of the antenna located within the expanding cloud.

High rate sampling detectors measuring the potential difference between the main body and boom antennas of interplanetary spacecraft have been shown to be efficient means to measure the voltage pulses induced by nano dust impacts on the spacecraft body itself (see Meyer-Vernet et al, Solar Phys. 256, 463 (2009)). However, rough estimates of the free charge liberated in post impact expanding plasma cloud indicate that the cloud’s own internal electrostatic field is too weak to account for measured pulses as the ones from the TDS instrument on the STEREO spacecraft frequently exceeding 0.1 V/m$. In this paper we argue that the detected pulses are not a direct measure of the potential structure of the plasma cloud, but are rather the consequence of a transitional interruption of the photoelectron return current towards the portion of the antenna located within the expanding cloud.

Our aim is to assess the benefits and limitations of using the redundant visibility information in regular phased array systems for improving the calibration. Regular arrays offer the possibility to use redundant visibility information to constrain the calibration of the array independent of a sky model and a beam models of the station elements. It requires a regular arrangement in the configuration of array elements and identical beam patterns. We revised a calibration method for phased array stations using the redundant visibility information in the system and applied it successfully to a LOFAR station. The performance and limitations of the method were demonstrated by comparing its use on real and simulated data. The main limitation is the mutual coupling between the station elements, which leads to non-identical beams and stronger baseline dependent noise. Comparing the variance of the estimated complex gains with the Cramer-Rao Bound (CRB) indicates that redundancy is a stable and optimum method for calibrating the complex gains of the system. Our study shows that the use of the redundant visibility does improve the quality of the calibration in phased array systems. In addition it provides a powerful tool for system diagnostics. Our results demonstrate that designing redundancy in both the station layout and the array configuration of future aperture arrays is strongly recommended. In particular in the case of the Square Kilometre Array with its dynamic range requirement which surpasses any existing array by an order of magnitude.

We present a method to optimize absorption cells for precise wavelength calibration in the near-infrared. We apply it to design and optimize methane isotopologue cells for precision radial velocity measurements in the K band. We also describe the construction and installation of two such cells for the CSHELL spectrograph at NASA’s IRTF. We have obtained their high-resolution laboratory spectra, which we can then use in precision radial velocity measurements and which can also have other applications. In terms of obtainable RV precision methane should out-perform other proposed cells, such as the ammonia cell ($^{14}$NH$_{3}$) recently demonstrated on CRIRES/VLT. The laboratory spectra of Ammonia and the Methane cells show strong absorption features in the H band that could also be exploited for precision Doppler measurements. We present spectra and preliminary radial velocity measurements obtained during our first-light run. These initial results show that a precision down to 20-30 m s$^{-1}$ can be obtained using a wavelength interval of only 5 nm in the K band and S/N$\sim$150. This supports the prediction that a precision down to a few m s$^{-1}$ can be achieved on late M dwarfs using the new generation of NIR spectrographs, thus enabling the detection of terrestrial planets in their habitable zones. Doppler measurements in the NIR can also be used to mitigate the radial velocity jitter due to stellar activity enabling more efficient surveys on young active stars.

We investigate the orbital stability of a putative Jovian planet reported by Ramm et. al in a compact binary \nu Octantis. Our numerical study makes use of a new computational Message Passing Interface (MPI) framework Mechanic which we developed to run massive numerical experiments on CPU clusters. The code is illustrated on a model Hamiltonian introduced by Froeschl\’e et al. We confirm that the \nu Octantis planet could reside in a retrograde orbit, according with a hypothesis of Eberle & Cuntz. It may be present in a zone of stable motions which has a structure of the Arnold web formed due to overlapping of low-order mean motion resonances and their sub-resonances. We also re-analyzed the available radial velocity data in terms of self-consistent Newtonian N-body model. We found stable best-fit solutions that obey the observational constraints. They correspond to retrograde, strictly anti-aligned orbits of the binary and the planet. However, these solutions are confined in very small stable regions of the phase space. The presence of a real planet in the system is still questionable, because its formation would be hindered by strong dynamical perturbations.

We investigate the orbital stability of a putative Jovian planet reported by Ramm et. al in a compact binary \nu Octantis. Our numerical study makes use of a new computational Message Passing Interface (MPI) framework Mechanic which we developed to run massive numerical experiments on CPU clusters. The code is illustrated on a model Hamiltonian introduced by Froeschl\’e et al. We confirm that the \nu Octantis planet could reside in a retrograde orbit, according with a hypothesis of Eberle & Cuntz. It may be present in a zone of stable motions which has a structure of the Arnold web formed due to overlapping of low-order mean motion resonances and their sub-resonances. We also re-analyzed the available radial velocity data in terms of self-consistent Newtonian N-body model. We found stable best-fit solutions that obey the observational constraints. They correspond to retrograde, strictly anti-aligned orbits of the binary and the planet. However, these solutions are confined in very small stable regions of the phase space. The presence of a real planet in the system is still questionable, because its formation would be hindered by strong dynamical perturbations.

Long baseline laser interferometers used for gravitational wave detection have proven to be very complicated to control. In order to have sufficient sensitivity to astrophysical gravitational waves, a set of multiple coupled optical cavities comprising the interferometer must be brought into resonance with the laser field. A set of multi-input, multi-output servos then lock these cavities into place via feedback control. This procedure, known as lock acquisition, has proven to be a vexing problem and has reduced greatly the reliability and duty factor of the past generation of laser interferometers. In this article, we describe a technique for bringing the interferometer from an uncontrolled state into resonance by using harmonically related external fields to provide a deterministic hierarchical control. This technique reduces the effect of the external seismic disturbances by four orders of magnitude and promises to greatly enhance the stability and reliability of the current generation of gravitational wave detector. The possibility for using multi-color techniques to overcome current quantum and thermal noise limits is also discussed.

We describe the use of Graphics Processing Units (GPUs) for speeding up the code NBODY6 which is widely used for direct $N$-body simulations. Over the years, the $N^2$ nature of the direct force calculation has proved a barrier for extending the particle number. Following an early introduction of force polynomials and individual time-steps, the calculation cost was first reduced by the introduction of a neighbour scheme. After a decade of GRAPE computers which speeded up the force calculation further, we are now in the era of GPUs where relatively small hardware systems are highly cost-effective. A significant gain in efficiency is achieved by employing the GPU to obtain the so-called regular force which typically involves some 99 percent of the particles, while the remaining local forces are evaluated on the host. However, the latter operation is performed up to 20 times more frequently and may still account for a significant cost. This effort is reduced by parallel SSE/AVX procedures where each interaction term is calculated using mainly single precision. We also discuss further strategies connected with coordinate and velocity prediction required by the integration scheme. This leaves hard binaries and multiple close encounters which are treated by several regularization methods. The present nbody6-GPU code is well balanced for simulations in the particle range $10^4-2 \times 10^5$ for a dual GPU system attached to a standard PC.

In this paper we explore the impact of including redshift information on cosmological applications with the forthcoming generation of large-scale, deep radio continuum surveys. By cross-matching these radio surveys with shallow optical to near-infrared surveys we can essentially separate the source distribution into a low redshift sample and the high-z tail of the radio sources that are unidentified, thus providing a constraint on the evolution of cosmological parameters such as those of dark energy. We examine two radio surveys, the Evolutionary Map of the Universe (EMU) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). A crucial advantage is their combined potential to provide a deep, full-sky survey. The surveys used for the cross-identifications are SkyMapper and SDSS, for the southern and northern skies, respectively. We concentrate on the galaxy clustering angular power spectrum as our benchmark observable and find that the possibility of including this low redshift information yields major improvements in the results. With this approach, we are able to put strict constraints on the dark energy parameters, i.e. w_0=-0.9+/-0.065(0.087) and w_a=-0.24+/-0.19(0.26) with(without) priors from Planck; this corresponds to a Figure of Merit (FoM) of circa 400(>200), which is two to three orders of magnitude times better than the case without any redshift information and more than three times better than what obtained by using only the cross-identified sources.

A new sampling theorem on the sphere has been developed recently, reducing the number of samples required to represent a band-limited signal by a factor of two for equiangular sampling schemes. For signals sparse in a spatially localised measure, such as in a wavelet basis, overcomplete dictionary, or in the magnitude of their gradient, for example, a reduction in the number of samples required to represent a band-limited signal has important implications for sparse signal reconstruction on the sphere. A more efficient sampling of the sphere improves the fidelity of sparse signal reconstruction through both the dimensionality and spatial sparsity of signals. To demonstrate this result we consider a simple inpainting problem on the sphere and consider signals sparse in the magnitude of their gradient. We develop a framework for total variation (TV) inpainting on the sphere by making a connection to the underlying continuous signal via a sampling theorem. Numerical simulations are performed, verifying the enhanced fidelity of sparse signal reconstruction due to the more efficient sampling of the sphere provided by the new sampling theorem.

Directional detection is a promising Dark Matter search strategy. Even though it could accommodate to a sizeable background contamination, electron/recoil discrimination remains a key and challenging issue as for direction-insensitive detectors. The measurement of the 3D track may be used to discriminate electrons from nuclear recoils. While a high rejection power is expected above 20 keV ionization, a dedicated data analysis is needed at low energy. After identifying discriminant observables, a multivariate analysis, namely a Boosted Decision Tree, is proposed, enabling an efficient event tagging for Dark Matter search. We show that it allows us to optimize rejection while keeping a rather high efficiency which is compulsory for rare event search.With respect to a sequential analysis, the rejection is about 20 times higher with a multivariate analysis, for the same Dark Matter exclusion limit.

The past decade has seen a dramatic improvement in the quality of data available at both high (HE: 100 MeV to 100 GeV) and very high (VHE: 100 GeV to 100 TeV) gamma-ray energies. With three years of data from the Fermi Large Area Telescope (LAT) and deep pointed observations with arrays of Cherenkov telescope, continuous spectral coverage from 100 MeV to $\sim10$ TeV exists for the first time for the brightest gamma-ray sources. The Fermi-LAT is likely to continue for several years, resulting in significant improvements in high energy sensitivity. On the same timescale, the Cherenkov Telescope Array (CTA) will be constructed providing unprecedented VHE capabilities. The optimisation of CTA must take into account competition and complementarity with Fermi, in particularly in the overlapping energy range 10$-$100 GeV. Here we compare the performance of Fermi-LAT and the current baseline CTA design for steady and transient, point-like and extended sources.

We develop an exact wavelet transform on the three-dimensional ball (i.e. on the solid sphere), which we name the flaglet transform. For this purpose we first construct an exact harmonic transform on the radial line using damped Laguerre polynomials and develop a corresponding quadrature rule. Combined with the spherical harmonic transform, this approach leads to a sampling theorem on the ball and a novel three-dimensional decomposition which we call the Fourier-Laguerre transform. We relate this new transform to the well-known Fourier-Bessel decomposition and show that band-limitness in the Fourier-Laguerre basis is a sufficient condition to compute the Fourier-Bessel decomposition exactly. We then construct the flaglet transform on the ball through a harmonic tiling, which is exact thanks to the exactness of the Fourier-Laguerre transform (from which the name flaglets is coined). The corresponding wavelet kernels have compact localisation properties in real and harmonic space and their angular aperture is invariant under radial translation. We introduce a multiresolution algorithm to perform the flaglet transform rapidly, while capturing all information at each wavelet scale in the minimal number of samples on the ball. Our implementation of these new tools achieves floating point precision and is made publicly available. We perform numerical experiments demonstrating the speed and accuracy of these libraries and illustrate their capabilities on a simple denoising example.

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is an airborne observatory consisting of a specially modified Boeing 747SP with a 2.7-m telescope, flying at altitudes as high as 13.7 km (45,000 ft). Designed to observe at wavelengths from 0.3 micron to 1.6 mm, SOFIA operates above 99.8 % of the water vapor that obscures much of the infrared and submillimeter. SOFIA has seven science instruments under development, including an occultation photometer, near-, mid-, and far-infrared cameras, infrared spectrometers, and heterodyne receivers. SOFIA, a joint project between NASA and the German Aerospace Center DLR, began initial science flights in 2010 December, and has conducted 30 science flights in the subsequent year. During this early science period three instruments have flown: the mid-infrared camera FORCAST, the heterodyne spectrometer GREAT, and the occultation photometer HIPO. This article provides an overview of the observatory and its early performance.

The closest examples of high-mass star birth occurs in deeply embedded environments at kiloparsec distances. Although much progress has been made, an observationally validated picture of the dominant processes which allows the central hydrostatic object to grow in mass has yet to be established. The observational technique of optical interferometry has demonstrated its potential in the field of high-mass star formation by delivering a milli-arcsecond infrared view on the complex accretion environment. We provide an overview of the scientific results obtained with multi-aperture telescope arrays and briefly discuss future instruments and their anticipated impact on our understanding of massive young stellar objects.

A direct WIMP (Weakly Interacting Massive Particle) detector with a neutron veto system is designed to better reject neutrons. An experimental configuration is studied in the present paper: 984 Ge modules are placed inside a reactor neutrino detector. The neutrino detector is used as a neutron veto device. The neutron background for the experimental design has been estimated using the Geant4 simulation. The result show that the neutron background can decrease to O(0.01) events per year per tonne of high purity Germanium. We calculate the sensitivity to spin-independent WIMP-nucleon elastic scattering. An exposure of one tonne $\times$ year could reach a cross-section of about 4$\times$$10^{-11}$ pb.

Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using forward magnetohydrodynamic simulation and synthetic white-light imaging. Identifiable as an inclined streak of enhanced brightness in a time-elongation map, the travelling shock can be readily imaged by an observer located within a wide range of longitudes in the ecliptic. Different parts of the shock front contribute to the imaged brightness pattern viewed by observers at different longitudes. Moreover, even for an observer located at a fixed longitude, a different part of the shock front will contribute to the imaged brightness at any given time. The observed brightness within each imaging pixel results from a weighted integral along its corresponding ray-path. It is possible to infer the longitudinal location of the shock from the brightness pattern in an optical sky map, based on the east-west asymmetry in its brightness and degree of polarization. Therefore, measurement of the interplanetary polarized brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging.

We present a framework to interactively volume-render three-dimensional data cubes using distributed ray-casting and volume bricking over a cluster of workstations powered by one or more graphics processing units (GPUs) and a multi-core CPU. The main design target for this framework is to provide an in-core visualization solution able to provide three-dimensional interactive views of terabyte-sized data cubes. We tested the presented framework using a computing cluster comprising 64 nodes with a total of 128 GPUs. The framework proved to be scalable to render a 204 GB data cube with an average of 30 frames per second. Our performance analyses also compare between using NVIDIA Tesla 1060 and 2050 GPU architectures and the effect of increasing the visualization output resolution on the rendering performance. Although our initial focus, and the examples presented in this work, is volume rendering of spectral data cubes from radio astronomy, we contend that our approach has applicability to other disciplines where close to real-time volume rendering of terabyte-order 3D data sets is a requirement.

We present EzGal, a flexible python program designed to easily generate observable parameters (magnitudes, colors, mass-to-light ratios) for any stellar population synthesis (SPS) model. As has been demonstrated by various authors, the choice of input SPS models can be a significant source of systematic uncertainty. A key strength of EzGal is that it enables simple, direct comparison of different models sets. EzGal is also capable of generating composite stellar population models (CSPs) and can interpolate between metallicities for a given model set. We have created a web interface to run EzGal and generate observables for a variety of star formation histories and model sets. We make many commonly used SPS models available from this interface; the BC03 models, an updated version of these models, the Maraston models, the BaSTI models, and finally the FSPS models. We use EzGal to compare magnitude predictions for the model sets as a function of wavelength, age, metallicity, and star formation history. We recover the well-known result that the models agree best in the optical for old, solar metallicity models, with differences at the ~0.1 magnitude level. The most problematic regime for SPS modeling is for young ages (7500 Angstroms) where scatter between models can vary from 0.3 mags (Sloan i) to 0.7 mags (Ks). We find that these differences are best understood as general uncertainties in SPS modeling. Finally we explore a more physically motivated example by generating CSPs with a star formation history matching the global star formation history of the universe. We demonstrate that the wavelength and age dependence of SPS model uncertainty translates into a redshift dependent model uncertainty, highlighting the importance of a quantitative understanding of model differences when comparing observations to models as a function of redshift.

We present EzGal, a flexible python program designed to easily generate observable parameters (magnitudes, colors, mass-to-light ratios) for any stellar population synthesis (SPS) model. As has been demonstrated by various authors, the choice of input SPS models can be a significant source of systematic uncertainty. A key strength of EzGal is that it enables simple, direct comparison of different models sets. EzGal is also capable of generating composite stellar population models (CSPs) and can interpolate between metallicities for a given model set. We have created a web interface to run EzGal and generate observables for a variety of star formation histories and model sets. We make many commonly used SPS models available from this interface; the BC03 models, an updated version of these models, the Maraston models, the BaSTI models, and finally the FSPS models. We use EzGal to compare magnitude predictions for the model sets as a function of wavelength, age, metallicity, and star formation history. We recover the well-known result that the models agree best in the optical for old, solar metallicity models, with differences at the ~0.1 magnitude level. The most problematic regime for SPS modeling is for young ages (7500 Angstroms) where scatter between models can vary from 0.3 mags (Sloan i’) to 0.7 mags (Ks). We find that these differences are best understood as general uncertainties in SPS modeling. Finally we explore a more physically motivated example by generating CSPs with a star formation history matching the global star formation history of the universe. We demonstrate that the wavelength and age dependence of SPS model uncertainty translates into a redshift dependent model uncertainty, highlighting the importance of a quantitative understanding of model differences when comparing observations to models as a function of redshift.

Over the past two decades, every scholarly publisher has migrated at least the mechanical aspects of their journal publishing so that they utilize digital means. The academy was comfortable with that for a while, but publishers are under increasing pressure to adapt further. At the American Astronomical Society (AAS), we think that means bringing our publishing program to the point of being fully digital, by establishing procedures and policies that regard the digital objects of publication primarily. We have always thought about our electronic journals as databases of digital articles, from which we can publish and syndicate articles one at a time, and we must now put flesh on those bones by developing practices that are consistent with the realities of article at a time publication online. As a learned society that holds the long-term rights to the literature, we have actively taken responsibility for the preservation of the digital assets that constitute our journals, and in so doing we have not forsaken the legacy pre-digital assets. All of us who serve as the long-term stewards of scholarship must begin to evolve into fully digital publishers.

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

We present a new computational scheme aimed at reducing the complexity of the chemical networks in astrophysical models, one which is shown to markedly improve their computational efficiency. It contains a flux-reduction scheme that permits to deal with both large and small systems. This procedure is shown to yield a large speed-up of the corresponding numerical codes and provides good accord with the full network results. We analyse and discuss two examples involving chemistry networks of the interstellar medium and show that the results from the present reduction technique reproduce very well the results from fuller calculations.

In this work, we have explored the advantages and drawbacks of using GPUs instead of CPUs in the calculation of a standard 2-point correlation function algorithm, which is useful for the analysis of Large Scale Structure of galaxies. Taking into account the huge volume of data foreseen in upcoming surveys, our main goal has been to accelerate significantly the analysis codes. We find that GPUs offer a 100-fold increase in speed with respect to a single CPU without a significant deviation in the results. For comparison’s sake, an MPI version was developed as well. Some issues, like code implementation, which arise from using this option are discussed.

We develop and demonstrate a classification system constituted by several Support Vector Machines (SVM) classifiers, which can be applied to select quasar candidates from large sky survey projects, such as SDSS, UKIDSS, GALEX. How to construct this SVM classification system is presented in detail. When the SVM classification system works on the test set to predict quasar candidates, it acquires the efficiency of 93.21% and the completeness of 97.49%. In order to further prove the reliability and feasibility of this system, two chunks are randomly chosen to compare its performance with that of the XDQSO method used for SDSS-III’s BOSS. The experimental results show that the high faction of overlap exists between the quasar candidates selected by this system and those extracted by the XDQSO technique in the dereddened i-band magnitude range between 17.75 and 22.45, especially in the interval of dereddened i-band magnitude < 20.0. In the two test areas, 57.38% and 87.15% of the quasar candidates predicted by the system are also targeted by the XDQSO method. Similarly, the prediction of subcategories of quasars according to redshift achieves a high level of overlap with these two approaches. Depending on the effectiveness of this system, the SVM classification system can be used to create the input catalog of quasars for the GuoShouJing Telescope (LAMOST) or other spectroscopic sky survey projects. In order to get higher confidence of quasar candidates, cross-result from the candidates selected by this SVM system with that by XDQSO method is applicable.

We report on observations of a free-Shercliff-layer instability in a Taylor-Couette experiment using a liquid metal over a wide range of Reynolds numbers, $Re\sim 10^3-10^6$. The free Shercliff layer is formed by imposing a sufficiently strong axial magnetic field across a pair of differentially rotating axial endcap rings. This layer is destabilized by a hydrodynamic Kelvin-Helmholtz-type instability, characterized by velocity fluctuations in the $r-\theta$ plane. The instability appears with an Elsasser number above unity, and saturates with an azimuthal mode number $m$ which increases with the Elsasser number. Measurements of the structure agree well with 2D global linear mode analyses and 3D global nonlinear simulations. These observations have implications for a range of rotating MHD systems in which similar shear layers may be produced.

We examine systematic differences in the derived X-ray properties of galaxy clusters as reported by three different groups: Vikhlinin et al. (2009a), Mantz et al. (2010b), and Planck Collaboration (2011b). The sample overlap between any two pairs of works ranges between 16 to 28 galaxy clusters in common. We find systematic differences in most reported properties, including the total cluster mass, M500. The most extreme case is an average 45% \pm 5% difference in cluster mass between the Planck Collaboration (2011b) and Mantz et al. (2010b), for clusters at z > 0.13 (averaged over 16 clusters). These mass differences induce differences in cluster observables defined within an R500 aperture. After accounting for aperture differences, we find very good agreement in gas mass estimates between the different groups. However, the soft-band X-ray luminosity, LX, core-excised spectroscopic temperature, TX, and gas thermal energy, YX = MgasTX display mean differences at the 5%-15% level. We also find that the low (z \leq 0.13) and high (z \geq 0.13) galaxy cluster samples in Planck Collaboration (2011b) appear to be systematically different: the YSZ/YX ratio for these two sub- samples is ln(YSZ/YX) = -0.06 \pm 0.04 and ln(YSZ/YX) = 0.08 \pm 0.04 for the low and high redshift sub-samples respectively.

As part of the DarkSide program of direct dark matter searches using liquid argon TPCs, a prototype detector with an active volume containing 10 kg of liquid argon, DarkSide-10, was built and operated underground in the Gran Sasso National Laboratory in Italy. A critically important parameter for such devices is the scintillation light yield, as photon statistics limits the rejection of electron-recoil backgrounds by pulse shape discrimination. We have measured the light yield of DarkSide-10 using the readily-identifiable full-absorption peaks from gamma ray sources combined with single-photoelectron calibrations using low-occupancy laser pulses. For gamma lines of energies in the range 122-1275 keV, we get consistent light yields averaging 8.887\pm0.003(stat)\pm0.444(sys) p.e./keV_ee. With additional purification, the light yield measured at 511 keV increased to 9.142\pm0.006(stat) p.e./keV_ee.

We have measured electron-ion recombination for Fe XII forming Fe XI using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. The measured merged beams recombination rate coefficient (MBRRC) for collision energies from 0 to 1500 eV is presented. This work uses a new method for determining the absolute MBRRC based on a comparison of the ion beam decay rate with and without the electron beam on. For energies below 75 eV, the spectrum is dominated by dielectronic recombination (DR) resonances associated with 3s-3p and 3p-3d core excitations. At higher energies we observe contributions from 3-N’ and 2-N’ core excitations DR. We compare our experimental results to state-of-the-art multi-configuration Breit-Pauli (MCBP) calculations and find significant differences, both in resonance energies and strengths. We have extracted the DR contributions from the measured MBRRC data and transformed them into a plasma recombination rate coefficient (PRRC) for temperatures in the range of 10^3 to 10^7 K. We show that the previously recommended DR data for Fe XII significantly underestimate the PRRC at temperatures relevant for both photoionized plasmas (PPs) and collisionaly ionized plasmas (CPs). This is to be contrasted with our MCBP PRRC results which agree with the experiment to within 30% at PP temperatures and even better at CP temperatures. We find this agreement despite the disagreement shown by the detailed comparison between our MCBP and experimental MBRRC results. Lastly, we present a simple parameterized form of the experimentally derived PRRC for easy use in astrophysical modelling codes.

The impact of the solar activity on the heliosphere has a strong influence on the modulation of the flux of low energy galactic cosmic rays arriving at Earth. Different instruments, such as neutron monitors or muon detectors, have been recording the variability of the cosmic ray flux at ground level for several decades. Although the Pierre Auger Observatory was designed to observe cosmic rays at the highest energies, it also records the count rates of low energy secondary particles (the scaler mode) for the self-calibration of its surface detector array. From observations using the scaler mode at the Pierre Auger Observatory, modulation of galactic cosmic rays due to solar transient activity has been observed (e.g., Forbush decreases). Due to the high total count rate coming from the combined area of its detectors, the Pierre Auger Observatory (its detectors have a total area greater than $16\,000$\,m$^2$) detects a flux of secondary particles of the order of $\sim 10^8$\,counts per minute. Time variations of the cosmic ray flux related to the activity of the heliosphere can be determined with high accuracy. In this paper we briefly describe the scaler mode and analyze a Forbush decrease together with the interplanetary coronal mass ejection that originated it. The Auger scaler data are now publicly available.

Interactive visualization and simulation of astrophysical phenomena help astronomers and enable digital planetariums and television documentaries to take their spectators on a journey into deep space to explore the astronomical wonders of our universe in 3D.

Interactive visualization and simulation of astrophysical phenomena help astronomers and enable digital planetariums and television documentaries to take their spectators on a journey into deep space to explore the astronomical wonders of our universe in 3D.

Interactive visualization and simulation of astrophysical phenomena help astronomers and enable digital planetariums and television documentaries to take their spectators on a journey into deep space to explore the astronomical wonders of our universe in 3D.

For pulsar projects it is often necessary to predict the pulse phase in advance, for example, when preparing for new observations. Interpolation of the pulse phase between existing measurements is also often required, for example, when folding X-ray or gamma-ray observations according to the radio pulse phase. Until now these procedures have been done using various ad hoc methods. The purpose of this paper is to show how to interpolate or predict the pulse phase optimally using statistical models of the various noise processes and the phase measurement uncertainty.

As part of NASA’s mission to explore habitable planets orbiting nearby stars, this paper explores the detection and characterization capabilities of a 4-m space telescope plus 50-m starshade located at the Earth-Sun L2 point, a.k.a. the New Worlds Observer (NWO). Our calculations include the true spectral types and distribution of stars on the sky, an iterative target selection protocol designed to maximize efficiency based on prior detections, and realistic mission constraints. We carry out both analytical calculations and simulated observing runs for a wide range in exozodiacal background levels ({\epsilon} = 1 – 100 times the local zodi brightness) and overall prevalence of Earth-like terrestrial planets ({\eta}\oplus = 0.1 – 1). We find that even without any return visits, the NWO baseline architecture (IWA = 65 mas, limiting FPB = 4\times10-11) can achieve a 95% probability of detecting and spectrally characterizing at least one habitable Earth-like planet, and an expectation value of ~3 planets found, within the mission lifetime and {\Delta}V budgets, even in the worst-case scenario ({\eta}\oplus = 0.1 and {\epsilon} = 100 zodis for every target). This achievement requires about one year of integration time spread over the 5 year mission, leaving the remainder of the telescope time for UV-NIR General Astrophysics. Cost and technical feasibility considerations point to a “sweet spot” in starshade design near a 50-m starshade effective diameter, with 12 or 16 petals, at a distance of 70,000-100,000 km from the telescope.

We report on the performance and commissioning of a cryogenic distillation column for low radioactivity underground argon at Fermi National Accelerator Laboratory. The distillation column is designed to accept a mixture of argon, helium, and nitrogen and return pure argon with a nitrogen contamination less than 10 ppm. In the first commissioning, we were able to run the distillation column in a continuous mode and produce argon that is 99.9% pure. After running in a batch mode, the argon purity was increased to 99.95%, with 500 ppm of nitrogen remaining. The efficiency of collecting the argon from the gas mixture was between 70% and 81%, at an argon production rate of 0.84-0.98 kg/day.

We report on the performance and commissioning of a cryogenic distillation column for low radioactivity underground argon at Fermi National Accelerator Laboratory. The distillation column is designed to accept a mixture of argon, helium, and nitrogen and return pure argon with a nitrogen contamination less than 10 ppm. In the first commissioning, we were able to run the distillation column in a continuous mode and produce argon that is 99.9% pure. After running in a batch mode, the argon purity was increased to 99.95%, with 500 ppm of nitrogen remaining. The efficiency of collecting the argon from the gas mixture was between 70% and 81%, at an argon production rate of 0.84-0.98 kg/day.

We present a new, approximate method of calculating the column density of ammonia in mapping observations of the 23 GHz inversion lines. The temperature regime typically found in star forming regions allows for the assumption of a slowly varying partition function for ammonia. It is therefore possible to determine the column density using only the (J=1,K=1) inversion transition rather than the typical combination of the (1,1) and (2,2) transitions, with additional uncertainties comparable to or less than typical observational error. The proposed method allows column density and mass estimates to be extended into areas of lower signal to noise ratio. We show examples of column density maps around a number of cores in the W3 and Perseus star-forming regions made using this approximation, along with a comparison to the corresponding results obtained using the full two-transition approach. We suggest that this method is a useful tool in studying the distribution of mass around YSOs, particularly in the outskirts of the protostellar envelope where the (2,2) ammonia line is often undetectable on the short timescales necessary for large area mapping.

We report on the first large-scale production of low radioactivity argon from underground gas wells. Low radioactivity argon is of general interest, in particular for the construction of large scale WIMP dark matter searches and detectors of reactor neutrinos for non-proliferation efforts. Atmospheric argon has an activity of about 1 Bq/kg from the decays of 39Ar; the concentration of 39Ar in the underground argon we are collecting is at least a factor of 100 lower than this value. The argon is collected from a stream of gas from a CO2 well in southwestern Colorado with a Vacuum Pressure Swing Adsorption (VPSA) plant. The gas from the well contains argon at a concentration of 400-600 ppm, and the VPSA plant produces an output stream with an argon concentration at the level of 30,000-50,000 ppm (3-5%) in a single pass. This gas is sent for further processing to Fermilab where it is purified by cryogenic distillation. The argon production rate is presently 0.5 kg/day.

We report on the first large-scale production of low radioactivity argon from underground gas wells. Low radioactivity argon is of general interest, in particular for the construction of large scale WIMP dark matter searches and detectors of reactor neutrinos for non-proliferation efforts. Atmospheric argon has an activity of about 1 Bq/kg from the decays of 39Ar; the concentration of 39Ar in the underground argon we are collecting is at least a factor of 100 lower than this value. The argon is collected from a stream of gas from a CO2 well in southwestern Colorado with a Vacuum Pressure Swing Adsorption (VPSA) plant. The gas from the well contains argon at a concentration of 400-600 ppm, and the VPSA plant produces an output stream with an argon concentration at the level of 30,000-50,000 ppm (3-5%) in a single pass. This gas is sent for further processing to Fermilab where it is purified by cryogenic distillation. The argon production rate is presently 0.5 kg/day.

A new code, named MAP, is written in Fortran language for magnetohydrodynamics (MHD) calculation with the adaptive mesh refinement (AMR) and Message Passing Interface (MPI) parallelization. There are several optional numerical schemes for computing the MHD part, namely, modified Mac Cormack Scheme (MMC), Lax-Friedrichs scheme (LF) and weighted essentially non-oscillatory (WENO) scheme. All of them are second order, two-step, component-wise schemes for hyperbolic conservative equations. The total variation diminishing (TVD) limiters and approximate Riemann solvers are also equipped. A high resolution can be achieved by the hierarchical block-structured AMR mesh. We use the extended generalized Lagrange multiplier (EGLM) MHD equations to reduce the non-divergence free error produced by the scheme in the magnetic induction equation. The numerical algorithms for the non-ideal terms, e.g., the resistivity and the thermal conduction, are also equipped in the MAP code. The details of the AMR and MPI algorithms are described in the paper.

We investigate the possibilities of wavelet-based RM Synthesis for the recognition of structures of regular and turbulent magnetic fields in extended magnetized objects, like galaxies and galaxy clusters. Wavelets allow to reformulate the RM Synthesis method in a scale-dependent way and to visualize the data as a function of Faraday depth and scale. We present observational tests to recognize regular magnetic fields without and with one or two reversals along the line of sight and imprints of turbulent magnetic fields. A region with a regular magnetic field generates a broad “disk” in Faraday space (“Faraday spectrum”), with two “horns” if the distribution of cosmic-ray electrons is broader than that of the thermal electrons. Each field reversal generates one asymmetric “horn” on top of the “disk”. A region with a turbulent field can be recognized by a “Faraday forest” of many components. We argue that the ratio of maximum to minimum wavelengths is an important parameter because it determines the range of scales which can be identified in Faraday space. Full recognition of magnetic field structures in spiral galaxies or galaxy clusters requires analysis of data cubes in position-position-Faraday depth (“PPF cubes”), observed over a wide and continuous wavelength range, from about 100 MHz to several GHz, providing good resolution as well as recognition of a wide range of scales in Faraday space. The planned SKA fulfils this condition and will be close to a perfect “Faraday telescope”. The combination of data from the present telescopes LOFAR (low frequencies) and EVLA (high frequencies) looks promising for the recognition of magnetic structures at all scales. The addition of WSRT or GMRT data at intermediate frequencies to those from LOFAR and EVLA data fills the gap between the LOFAR and EVLA wavelength ranges.

We investigate the possibilities of wavelet-based RM Synthesis for the recognition of structures of regular and turbulent magnetic fields in extended magnetized objects, like galaxies and galaxy clusters. Wavelets allow to reformulate the RM Synthesis method in a scale-dependent way and to visualize the data as a function of Faraday depth and scale. We present observational tests to recognize regular magnetic fields without and with one or two reversals along the line of sight and imprints of turbulent magnetic fields. A region with a regular magnetic field generates a broad “disk” in Faraday space (“Faraday spectrum”), with two “horns” if the distribution of cosmic-ray electrons is broader than that of the thermal electrons. Each field reversal generates one asymmetric “horn” on top of the “disk”. A region with a turbulent field can be recognized by a “Faraday forest” of many components. We argue that the ratio of maximum to minimum wavelengths is an important parameter because it determines the range of scales which can be identified in Faraday space. Full recognition of magnetic field structures in spiral galaxies or galaxy clusters requires analysis of data cubes in position-position-Faraday depth (“PPF cubes”), observed over a wide and continuous wavelength range, from about 100 MHz to several GHz, providing good resolution as well as recognition of a wide range of scales in Faraday space. The planned SKA fulfils this condition and will be close to a perfect “Faraday telescope”. The combination of data from the present telescopes LOFAR (low frequencies) and EVLA (high frequencies) looks promising for the recognition of magnetic structures at all scales. The addition of WSRT or GMRT data at intermediate frequencies to those from LOFAR and EVLA data fills the gap between the LOFAR and EVLA wavelength ranges.

We survey the job situation of women in astronomy in Germany and of German women abroad and review indicators for their career development. Our sample includes women astronomers from all academic levels from doctoral students to professors, as well as female astronomers who have left the field. We find that networking and human support are among the most important factors for success. Experience shows that students should carefully choose their supervisor and collect practical knowledge abroad. We reflect the private situation of female German astronomers and find that prejudices are abundant, and are perceived as discriminating.We identify reasons why women are more likely than men to quit astronomy after they obtain their PhD degree. We give recommendations to young students on what to pay attention to in order to be on the successful path in astronomy.

Modeling the radiation generated by accreting matter is an important step towards realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al. (2008), have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.

Modeling the radiation generated by accreting matter is an important step towards realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al. (2008), have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.

We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified interface. Our framework is intrinsically parallel, and conveniently separates all component numerical modules in memory. This strict separation allows automatic unit conversion, distributed execution of modules on different cores within a cluster or grid, and orderly recovery from errors. The framework can be efficiently implemented and incurs an acceptable overhead. In practice, we measure the time spent in the framework to be less than 1% of the wall-clock time. Due to the unified structure of the interface, incorporating multiple modules addressing the same physics in different ways is relatively straightforward. Different modules may be advanced serially or in parallel. Despite initial concerns, we have encountered relatively few problems with this strict separation between modules, and the results of our simulations are consistent with earlier results using more traditional monolithic approaches. This framework provides a platform to combine existing simulation codes or develop new physical solver codes within a rich “ecosystem” of interchangeable modules.

Brown has recently introduced a covariant formulation of the BSSN equations which is well suited for curvilinear coordinate systems. This is particularly desirable as many astrophysical phenomena are symmetric with respect to the rotation axis or are such that curvilinear coordinates adapt better to their geometry. However, the singularities associated with such coordinate systems are known to lead to numerical instabilities unless special care is taken (e.g., regularization at the origin). Cordero-Carrion will present a rigorous derivation of partially implicit Runge-Kutta methods in forthcoming papers, with the aim of treating numerically the stiff source terms in wave-like equations that may appear as a result of the choice of the coordinate system. We have developed a numerical code solving the BSSN equations in spherical symmetry and the general relativistic hydrodynamic equations written in flux-conservative form. A key feature of the code is that it uses a second-order partially implicit Runge-Kutta method to integrate the evolution equations. We perform and discuss a number of tests to assess the accuracy and expected convergence of the code, namely a pure gauge wave, the evolution of a single black hole, the evolution of a spherical relativistic star in equilibrium, and the gravitational collapse of a spherical relativistic star leading to the formation of a black hole. We obtain stable evolutions of regular spacetimes without the need for any regularization algorithm at the origin.