Electrode level Monte Carlo model of radiation damage effects on astronomical CCDs

Current optical space telescopes rely upon silicon Charge Coupled Devices (CCDs) to detect and image the incoming photons. The performance of a CCD detector depends on its ability to transfer electrons through the silicon efficiently, so that the signal from every pixel may be read out through a single amplifier. This process of electron transfer is highly susceptible to the effects of solar proton damage (or non-ionizing radiation damage). This is because charged particles passing through the CCD displace silicon atoms, introducing energy levels into the semi-conductor bandgap which act as localized electron traps. The reduction in Charge Transfer Efficiency (CTE) leads to signal loss and image smearing. The European Space Agency's astrometric Gaia mission will make extensive use of CCDs to create the most complete and accurate stereoscopic map to date of the Milky Way. In the context of the Gaia mission CTE is referred to with the complementary quantity Charge Transfer Inefficiency (CTI = 1-CTE). CTI is an extremely important issue that threatens Gaia's performances. We present here a detailed Monte Carlo model which has been developed to simulate the operation of a damaged CCD at the pixel electrode level. This model implements a new approach to both the charge density distribution within a pixel and the charge capture and release probabilities, which allows the reproduction of CTI effects on a variety of measurements for a large signal level range in particular for signals of the order of a few electrons. A running version of the model as well as a brief documentation and a few examples are readily available at http://www.strw.leidenuniv.nl/~prodhomme/cemga.php as part of the CEMGA java package (CTI Effects Models for Gaia).Comment: Accepted by MNRAS on 13 February 2011. 15 pages, 7 figures and 5 table

Comparing simulations and test data of a radiation damaged charge-coupled device for the Euclid mission

The visible imager instrument on board the Euclid mission is a weak-lensing experiment that depends on very precise shape measurements of distant galaxies obtained by a large charge-coupled device (CCD) array. Due to the harsh radiative environment outside the Earth’s atmosphere, it is anticipated that the CCDs over the mission lifetime will be degraded to an extent that these measurements will be possible only through the correction of radiation damage effects. We have therefore created a Monte Carlo model that simulates the physical processes taking place when transferring signals through a radiation-damaged CCD. The software is based on Shockley–Read–Hall theory and is made to mimic the physical properties in the CCD as closely as possible. The code runs on a single electrode level and takes the three-dimensional trap position, potential structure of the pixel, and multilevel clocking into account. A key element of the model is that it also takes device specific simulations of electron density as a direct input, thereby avoiding making any analytical assumptions about the size and density of the charge cloud. This paper illustrates how test data and simulated data can be compared in order to further our understanding of the positions and properties of the individual radiation-induced traps

The astrometric core solution for the Gaia mission. Overview of models, algorithms and software implementation

The Gaia satellite will observe about one billion stars and other point-like sources. The astrometric core solution will determine the astrometric parameters (position, parallax, and proper motion) for a subset of these sources, using a global solution approach which must also include a large number of parameters for the satellite attitude and optical instrument. The accurate and efficient implementation of this solution is an extremely demanding task, but crucial for the outcome of the mission. We provide a comprehensive overview of the mathematical and physical models applicable to this solution, as well as its numerical and algorithmic framework. The astrometric core solution is a simultaneous least-squares estimation of about half a billion parameters, including the astrometric parameters for some 100 million well-behaved so-called primary sources. The global nature of the solution requires an iterative approach, which can be broken down into a small number of distinct processing blocks (source, attitude, calibration and global updating) and auxiliary processes (including the frame rotator and selection of primary sources). We describe each of these processes in some detail, formulate the underlying models, from which the observation equations are derived, and outline the adopted numerical solution methods with due consideration of robustness and the structure of the resulting system of equations. Appendices provide brief introductions to some important mathematical tools (quaternions and B-splines for the attitude representation, and a modified Cholesky algorithm for positive semidefinite problems) and discuss some complications expected in the real mission data.Comment: 48 pages, 19 figures. Accepted for publication in Astronomy & Astrophysic

<i>Gaia</i> Data Release 1. Summary of the astrometric, photometric, and survey properties

Context. At about 1000 days after the launch of Gaia we present the first Gaia data release, Gaia DR1, consisting of astrometry and photometry for over 1 billion sources brighter than magnitude 20.7. Aims. A summary of Gaia DR1 is presented along with illustrations of the scientific quality of the data, followed by a discussion of the limitations due to the preliminary nature of this release. Methods. The raw data collected by Gaia during the first 14 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium (DPAC) and turned into an astrometric and photometric catalogue. Results. Gaia DR1 consists of three components: a primary astrometric data set which contains the positions, parallaxes, and mean proper motions for about 2 million of the brightest stars in common with the HIPPARCOS and Tycho-2 catalogues – a realisation of the Tycho-Gaia Astrometric Solution (TGAS) – and a secondary astrometric data set containing the positions for an additional 1.1 billion sources. The second component is the photometric data set, consisting of mean G-band magnitudes for all sources. The G-band light curves and the characteristics of ∼3000 Cepheid and RR-Lyrae stars, observed at high cadence around the south ecliptic pole, form the third component. For the primary astrometric data set the typical uncertainty is about 0.3 mas for the positions and parallaxes, and about 1 mas yr−1 for the proper motions. A systematic component of ∼0.3 mas should be added to the parallax uncertainties. For the subset of ∼94 000 HIPPARCOS stars in the primary data set, the proper motions are much more precise at about 0.06 mas yr−1. For the secondary astrometric data set, the typical uncertainty of the positions is ∼10 mas. The median uncertainties on the mean G-band magnitudes range from the mmag level to ∼0.03 mag over the magnitude range 5 to 20.7. Conclusions. Gaia DR1 is an important milestone ahead of the next Gaia data release, which will feature five-parameter astrometry for all sources. Extensive validation shows that Gaia DR1 represents a major advance in the mapping of the heavens and the availability of basic stellar data that underpin observational astrophysics. Nevertheless, the very preliminary nature of this first Gaia data release does lead to a number of important limitations to the data quality which should be carefully considered before drawing conclusions from the data

Thermal and Optical Characterization of Undoped and Neodymium-Doped Y3ScAl4O12 Ceramics

Y3–3xNd3xSc1Al4O12 (x = 0, 0.01, and 0.02) ceramics were fabricated by sintering at high temperature under vacuum. Unit cell parameter refinement and chemical analysis have been performed. The morphological characterization shows micrograins with no visible defects. The thermal analysis of these ceramics is presented, by measuring the specific heat in the temperature range from 300 to 500 K. Their values at room temperature are in the range 0.81–0.90 J g1–K–1. The thermal conductivity has been determined by two methods: by the experimental measurement of the thermal diffusivity by the photopyroelectric method, and by spectroscopy, evaluating the thermal load. The thermal conductivities are in the range 9.7–6.5 W K–1 m–1 in the temperature interval from 300 to 500 K. The thermooptic coefficients were measured at 632 nm by the dark mode method using a prism coupler, and the obtained values are in the range 12.8–13.3 × 10–6 K–1. The nonlinear refractive index values at 795 nm have been evaluated to calibrate the nonlinear optical response of these materials.This work is supported by the Spanish Government under projects MAT2011-29255-C02-01-02, MAT2013-47395-C4-4-R, and the Catalan Government under project 2014SGR1358. It was also funded by the European Commission under the Seventh Framework Programme, project Cleanspace, FP7-SPACE-2010-1-GA No. 263044

Gaia Data Release 1: Testing parallaxes with local Cepheids and RR Lyrae stars

Context. Parallaxes for 331 classical Cepheids, 31 Type II Cepheids, and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). Aims. In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, which involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity (PL), period-Wesenheit (PW) relations for classical and Type II Cepheids and infrared PL, PL-metallicity (PLZ), and optical luminosity-metallicity (M V -[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. Methods. Classical Cepheids were carefully selected in order to discard known or suspected binary systems. The final sample comprises 102 fundamental mode pulsators with periods ranging from 1.68 to 51.66 days (of which 33 with σ Ω /Ω < 0.5). The Type II Cepheids include a total of 26 W Virginis and BL Herculis stars spanning the period range from 1.16 to 30.00 days (of which only 7 with σ Ω /Ω < 0.5). The RR Lyrae stars include 200 sources with pulsation period ranging from 0.27 to 0.80 days (of which 112 with σ Ω /Ω < 0.5). The new relations were computed using multi-band (V,I,J,K s ) photometry and spectroscopic metal abundances available in the literature, and by applying three alternative approaches: (i) linear least-squares fitting of the absolute magnitudes inferred from direct transformation of the TGAS parallaxes; (ii) adopting astrometry-based luminosities; and (iii) using a Bayesian fitting approach. The last two methods work in parallax space where parallaxes are used directly, thus maintaining symmetrical errors and allowing negative parallaxes to be used. The TGAS-based PL,PW,PLZ, and M V - [Fe/H] relations are discussed by comparing the distance to the Large Magellanic Cloud provided by different types of pulsating stars and alternative fitting methods. Results. Good agreement is found from direct comparison of the parallaxes of RR Lyrae stars for which both TGAS and HST measurements are available. Similarly, very good agreement is found between the TGAS values and the parallaxes inferred from the absolute magnitudes of Cepheids and RR Lyrae stars analysed with the Baade-Wesselink method. TGAS values also compare favourably with the parallaxes inferred by theoretical model fitting of the multi-band light curves for two of the three classical Cepheids and one RR Lyrae star, which were analysed with this technique in our samples. The K-band PL relations show the significant improvement of the TGAS parallaxes for Cepheids and RR Lyrae stars with respect to the Hipparcos measurements. This is particularly true for the RR Lyrae stars for which improvement in quality and statistics is impressive. Conclusions. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent the first Gaia-calibrated relations and form a work-in-progress milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia Data Release 2 (DR2) in 2018. © ESO, 2017

Gaia Data Release 1: Open cluster astrometry: performance, limitations, and future prospects

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information.Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters.Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed.Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier HIPPARCOS-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters.Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the HIPPARCOS data, with clearly increased luminosities for older A and F dwarfs