abstract: Time-Spatial data plays a crucial role for different fields such as traffic management. These data can be collected via devices such as surveillance sensors or tracking systems. However, how to efficiently analyze and visualize these data to capture essential embedded pattern information is becoming a big challenge today. Classic visualization approaches focus on revealing 2D and 3D spatial information and modeling statistical test. Those methods would easily fail when data become massive. Recent attempts concern on how to simply cluster data and perform prediction with time-oriented information. However, those approaches could still be further enhanced as they also have limitations for handling massive clusters and labels. In this paper, we propose a visualization methodology for mobility data using artificial neural net techniques. This method aggregates three main parts that are Back-end Data Model, Neural Net Algorithm including clustering method Self-Organizing Map (SOM) and prediction approach Recurrent Neural Net (RNN) for extracting the features and lastly a solid front-end that displays the results to users with an interactive system. SOM is able to cluster the visiting patterns and detect the abnormal pattern. RNN can perform the prediction for time series analysis using its dynamic architecture. Furthermore, an interactive system will enable user to interpret the result with graphics, animation and 3D model for a close-loop feedback. This method can be particularly applied in two tasks that Commercial-based Promotion and abnormal traffic patterns detection.
author: Zhenghao Chen; Jialong Zhou; Xiuying Wang
title: Visual Analytics of Movement Pattern Based on Time-Spatial Data: A Neural Net Approach

# Introduction

Time-spatial data consists of spatial and temporal contextual information, which is a complicated type of data set. Nowadays, with the development of technique, it is becoming more and more convenient to collect such data with detailed day and time information than it was ever before . Taking mobility data as an example which represents the objects movement in real time with recorded the geographical information , this type of data can be collected using Global Positioning System (GPS), sensors or tracking systems. The information embed from those data such as behavior pattern can be also implied in urban development, traffic transportation, and big event management. Therefore the understanding of this particular part of data is essential. To obtain meaningful results from those data, this paper proposes a methodology that utilizes neural network together with other visualization approaches. The framework of our method is shown as Figure. There are three majority parts. Firstly, this approach utilizes the relational database combined with dynamic non-relational database as back-end. Temporal database will be generated once user requests from the relational database. The second part is features extraction algorithm. We propose two approaches prediction and clustering. For clustering, SOM will be applied to cluster visiting pattern in order to group the visitors and find the abnormal groups that can be further suspected as crime or others. Moreover, the inspections will be put on such abnormal clusters. Regarding prediction, LSTM network is utilized to predict upcoming locations that travelers would like to go in next time-stamp. Eventually, the extracted information such as clustering and prediction result will be displayed using graphics and multi-media techniques including 3D time-spatial cube, 2D heat map, relational map and animation to users. Such displays will help users understand behavior pattern of moving objects better.

Using our approach, we participated in the public challenges that Visual Analytics Science and Technology Challenge 2017 (VAST 2017) as well as assist a mobile company in China for their promotion activity. In terms of commercial based visualization such as mobile users promotion , we used this method to obtain user behavior pattern. Based on those information, we can provide our suggestions to the company manger. On the other hand, we take traffic visualization and prediction in VAST 2017 to research how movement of vehicles affect wildlife via understanding their pattern by conducting the pipeline of visual analytics. Consequently, we provide suggestions to a Chinese mobile company for user promotion and successfully predict the affect of vehicles to wildlife in a nature preserve. Overall contributions can be summarized as:

-   We set up a time-spatial data visual analytics pipeline to understand movement behavior.

-   We modify and utilize neural net approach to propose a clustering and prediction methodology .

-   We use this approach to resolve two data analytics problems including commercial promotion as well as pattern detection and prediction.

# Related Work

The previous work are main concerned on two parts that clustering and aggregation as well as feature learning. .

## Clustering and Aggregation

Increasing amount of data collection, analysis of trajectories is becoming even more massive. Therefore, single mobility data analysis and visualization become un-meaningful. The state of the art research focus on aggregation, clustering and grouping in order to simplify time-spatial context information via reducing the size of data representation . With those advance techniques, the visualization of time-series data becomes much more elegant than before. The biggest challenge for those research is to generate an efficient algorithm for aggregation, grouping and clustering.

## Feature Learning

The big breakthrough was conducted with machine learning. Classic clustering methods such as Gaussian Mixture Model (GMM) with Expectation Maximization (EM) and K-means serves as effective tool. However, clustering raw trajectory with multiple-dimension is not sufficient for larger time-spatial data set . Recently, the attempts of using Neural Network Self-Organizing Map (SOM) to cluster the trajectory data is done in the 2012 . That research combines temporal-context information and other associated data such as social-economics background of users. Using such data, researchers tried to visualize the pattern of different groups of visitors and their potential to commit crime. In 2015, the T.Landersberger used spatial-simplification algorithm to dynamically reduce the flow information after aggregation and clustering which is the latest approach to achieve the aggregation task so far. State-of-art method deep learning enables the prediction of time-spatial for next timestamp. Researchers start looking at how to predict the upcoming location of movement and generate an inference or generative model. Recurrent Neural Net (RNN) is one of the popular deep learning models which makes it easier to deal with vectors where each dimension has strong dependency from previous one..

# A Dynamic to Static Data Model

There are two main types of data. For one, we have relational database as static database. That is all the data will be stored in the relational database model using relationship $Object - [Activity] -> StayPoint$. Also, both stay points and objects are stored using tree data structure. Specifically, stay point will have sub-stay points as generalization. Moreover, non-relational dynamic database will be generated from static database once user requests. Those data will consist of check-in frequency, time-spend, squence and Time-Oriented matrix.

**Visiting-Frequency Matrix** It is to indicate how many times a visitor visit in every location, we use visiting frequency matrix.

**Time-Spent Matrix** Simply using frequency count might lost a lot of information of visitor behavior information. Hence, we also need to consider the period that visitor spend in certain attraction.

**Sequence Vector** Sequence vector records the location that a traveler visits in certain period without considering time information. Such vector will simply track places a traveler check in using a one-dimensional vector.

**Time-Oriented Matrix** Time-oriented matrix have two dimensions for representing user-location-traveling-stay information. Specifically, X-dimension represents the users list and Y-Dimension represents the time that is divided into certain time duration. The value of matrix will be the discrete identification value of each location that a traveler stays in at certain time stamp.

# Neural Net Based Feature Extraction Algorithms

This section will introduce the Neural Net Based Feature Extraction Algorithms including SOM-based clustering and RNN-based Prediction, the algorithm is inspired from our previous work. We use these two different algorithms for generating predicted flow map and clustering Unified Matrix(U-matrix). Our main contribution is concerned on the post-training visualization rather than modifying the neural net.

## Self-organizing Map based Clustering

An unsupervised approach Self-organizing Map would be utilized here. Furthermore, most of previous work tried to cluster the movement trajectory directly. However, such clustering become very difficult for clustering large and irregular trajectories. Different from those, in this approach, I cluster moving people based on their preferences of locations (time spent and frequency matrix). This clustering method helps to find people who have the high potential to be in the same groups once they have similar or even the same interests. Based on those, we can further infer what category of visitors they are. For instance, a group of family with children moving in a park would probably have higher check-in frequency on kiddie rides than a group of teenagers.

Essentially, there are two parts for this work. For one, we use SOM to train and render a feature map after clustering. For another, after clustering result is rendered, we use associated visual analytics approach to recognize different patterns from the rendered Unified Distance Matrix from SOM. There are couple of essential steps including **Initialization**, **Competition**, **Cooperation** **Adaption** **Group Aggregation** and **Visiting Pattern Analysis**. The detail of those parts will be explained in the following sections.

In addition, associated visual analytics is crucial in this approach. Using those approaches to visualize the patterns such as visitor preferences on different locations will be useful to understand movement behaviors. Besides, other features including communication, time spent on a site, preference categories are also considered as essential features for clustering.

## Deep Learning based Prediction

As we process the time-oriented matrix as dynamic data which indicates where visitor stays at each timestamp. We will use this matrix as training data to forecast the travel flow of visitors. By using this, we are not going to lose any time information which results in better performance than using pure sequence matrix . That is because duration that individual stay in a location is an essential factor of behavior choice . Given an example, if we want to predict which location traveler would like to go at 18:00 pm and we have training time-oriented vector $(2, 2, 2, 3, 2, 41, 41, 2, 3, 5)$ for 10 hours and from 8:00 am to 18:00 pm. We can take vector $(2, 2, 2, 3, 2, 41, 41, 2, 3)$ as training set while 5 as targeted label as visitor is at location 5 at 18:00 pm. Therefore we can de-composite the time-oriented Matrix to be training matrix and labour vector

The second step is to use a prediction method to generate a model and then input the test data into model to generate the test label which is the predicted result. Afterwards we can visualize the result using several visualization approaches.The prediction method is the most crucial part. In this paper we use Recurrent Neural Net to serve as our predictor as its advantage to deal with time series data. Also the classical classification approaches K-Nearest Neighborhood (KNN), Naive Bayes (NB), Decision Tree (DT), Random Forest (RF), Support Vector Machine (SVM) and Adaptive Boosting (AdaBoost) will be used for evaluation.

# Experiment

In this part, we will introduce two of our experiments. For one, it is a mobile user commercial promotion visualization. For another, it is visualization of abnormal pattern of a nature preserve.

## Visual Analytics for Mobile User Commercial Promotion

### Data Set

This data set is from a Chinese Mobile company records the users use the mobile phone stations to visit the IP address in the real time. The data set is formed in $<userId, TimeStamp, Station, positon, IP>$.

### Outstanding User Detection

Using SOM to cluster the visiting pattern, we can get the result. There are around 150 stations activated per day. We query one station and cluster the visiting pattern for a week (5 working days) to view the result. As shown in Figure , the down left, down right and up right corners are shown as the outstanding users while the others can be interpreted as regular users.By further visualizing the visiting pattern of both outstanding and regular users, we see that, outstanding users use the station with a large number of frequency and keep activity on the station in quite long time. While the patterns of regular users illustrate that they only use the station with a very small number and activity time is very short too. Therefore, we know that if we want to have a commercial promotion, we should put more efforts on those outstanding users as they are more likely see the advertisements with more time spent on using the mobile station.

### Station Usage Prediction

Figure shows the prediction accuracy among the different methods. From the result, except the uniform method, other methods achieve very outstanding performance. The highest one is still obtained by RNN with 98% following by KNN, most frequency, RF and SVM with 97.5%, 97.2%, 96.4% and 96.2%. It is also interesting to see that, in this data set, most methods can obtain excellent performance. From that, we can also see the accuracy changes through time. RNN (indicated by green line) achieves mostly the highest result among all the time. Especially, when the training set is 300 minutes (5 hours), RNN becomes the best predictor. Also besides RNN, we can see that Decision Tree is also a excellent predictor in this task. Figure shows the actual usage of mobile stations at 17:00 pm, 14 May, 2015 as well as the predicted usage pattern of stations at the same time. We can see that RNN can mostly forecast the usage even with some slight mistakes. In this task, by successfully predicting the usage pattern of stations, we can provide suggestions to the managers. For instance, which station may have a large number of users in next hour, therefore more technical support should be concerned on. For the commercial promotion sides, once we know which station will be busy in next hour, commercial company can put advertisements onto that station.

## Pattern Analysis for a Nature Preserve

The data is obtained from VAST Challenge 2017[^1] which is trajectory data in a Nature Preserve Boonsong Lekagul. Locations of this preserve can be divided into five categories: Entrances, General-gates, Gates, Ranger-stops, and Camping. The vehicle that traveling around can be divided into 2-axle car (or motorcycle), 2-axle truck, 3-axle truck, 4-axle (and above) truck, 2-axle bus and 3-axle bus. Data can be formatted as $<Timestamp,car-id,car-type,gate-name>$.

### Abnormal Pattern Detection

Using the SOM clustering. It is quite obvious to see there are four outliers on the U-matrix map shown in Figure . These four points are all in the distinct clusters especially right up one which have the largest difference to others.By further visiting the heat map, we can visualize the frequency that those vehicles check-in such relationship can also be indicated by relationship map. 3D time-cube can clearly show trajectories of these four outstanding patterns. We can see that all of these four vehicles travel around in this nature preserve in a long duration. The most outstanding vehicle keeps visiting certain locations for 92 days. We can assume this vehicle will have strongest effect on the wildlife as the longest time it spent in this nature preserve. In other words, as this vehicle occurred over multiple days in the nature preserve, the animals around the path it traveled through would most likely to have been affected for a long time. Such pattern would have much more impact on wildlife than other regular traveling vehicle.

### Traffic Pattern Prediction

From Figure , we see that RNN get the most outstanding performance with 92.3%, other three methods SVM, KNN and Random Forest get 87.6%, 84.2% and 79% respectively for the second, third and fourth outstanding performance. Overall, for this data set, the average accuracy for most methods are quite good. Figure shows the actual traffic movement from 19:00 pm to 20:00 pm in period from 01 Jan 2016 to 29 May 2016 as well as the RNN predicted movement in the same time period. In this task, I use the data set from 01 May 2015 to 31 Dec 2015 from 0:00 am to 20:00 pm movement as training data to predict the attendance at 20:00 pm from 01 Jan 2016 to 29 May 2016 . In this task, RNN can successfully outline the four majority trajectories that all start from General Gate 6 but to Camping 5, Ranger Stop 2, General Gate 1 and Camping 1 respectively. As the time is 19:00-20:00 pm, we know that the night time wildlife in those four paths are going get most effects. Correspondingly, we will suggest the manager to take corresponding activities to help the wildlife among these four paths so that mitigate the effects.

# Conclusion and Future Work

The advantage of this framework over conventional methods is that we focus the feature extraction that are clustering and prediction. Such results will indicate essential information to user such as behavior pattern. Moreover, using neural net to render the clustering and predciton rather than simply constructing a visual analytics system, users do not have to conduct expensive experiments and suffer from massive information because of large time-spatial data set. Future work will be concerned on the how to enhance current method and what new methods should be introduced to resolve the temporal-spatial data visual analytics.

[^1]: <http://www.vacommunity.org/VAST+Challenge+2017>


abstract: The Smart Focal Planes (SmartFP) activity is a European Joint Research Activity funded to develop novel optical technologies for future large telescope instrumentation . In this paper, we will discuss the image slicer developments being carried out as part of this initiative. Image slicing techniques have many applications in the plans for instrumentation on Extremely Large Telescopes and will be central to the delivery of the science case. A study of a virtual *multi-object multi-ifu spectrograph and imager* (MOMSI) for a hypothetical OWL-class telescope reveals the need for focal plane splitting, deployable imagers and very small beam steering elements like deployable IFUs. The image slicer workpackage, lead from Durham University in collaboration with LFM Bremen, TNO Delft, UKATC Edinburgh, CRAL Lyon, LAM Marseille, Padua University and REFLEX Prague, is evaluating technologies for manufacturing micro optics in large numbers to enable multi-object integral field spectroscopy.
address: Center for Advanced Instrumentation (CfAI), Netpark Research Institute,Joseph Swan Road, Netpark, Sedgefield TS21 3FB, United Kingdom; INAF, Vicolo dell'Osservatorio 5, I - 35122 Padova, Italy; Reflex s.r.o., Novodvorska 994, 142 21 Praha 4, Czech Republic; Laboratoire d'Astrophysique de Marseille, Traverse du Siphon-Les trois Lucs, 13012 Marseille, France; UKATC, Blackford Hill, Edinburgh EH9 3HJ, Scotland; LFM, Badgasteiner Strasse 2, D-28359 Bremen, Germany
author: J.Schmoll, D.J.Robertson, C.M.Dubbeldam, J.Yao; F.Bortoletto; L.Pina, R.Hudec; E.Prieto; C.Norrie, S.Ramsay-Howat; W.Preuss
title: Optical replication techniques for image slicers

# MOMSI - a virtual instrument study

The **MOMSI** (**m**ulti-**o**bject **m**ulti-ifu **s**pectrograph and **i**mager) project is a virtual instrument, being designed for a 100m class ELT featuring an f/6 focus (, ). Due to the physical size of the focal plane (2.92 mm/arcsec, or 350 mm for 2 arcmin FOV), together with the desire to have diffraction-limited sampling, the amount of information available forces a pre-selection of regions of interest. Hence, even for imaging purposes pickoff systems would have to be considered. To get a high multiplex advantage, about 100 deployable IFUs are planned with about 100 $\times$ 100 spatial elements each. The challenges from the micro optics point of view are:

**Number of sub-systems:**

:   Each IFU will require a reimaging pickoff system and at least three mirror arrays, being the slicer itself, the pupil and slit mirror. Assuming an array size of 100 elements, there would be 30,000 optical surfaces just for the basic function of a slicer-type IFU. Manufacture and alignment will be critical, as well as the production time scale and costs.

**Size of individual components:**

:   To reduce the complexity and mass of the whole system, there are no dedicated fore-optics planned for these slicers. Hence the slice width will be in the order of 100 to 300 $\mu$m, being about a factor of 5 less than the typical existing IFU slit widths of e.g. 800 $\mu$m for the GNIRS IFU .

# The replication experiments

Given the prospect of large amounts of similar shaped surfaces, optical replication appears to be a very efficient way of manufacturing tens to hundreds of similar IFUs. Mirror replication techniques are well known for x-ray optics, because for these mirror types classical polishing of the positive forms is very difficult. Several tests were developed to check if this technique is useful for the replication of larger amounts of micro-optics. The replications are performed by two different companies, REFLEX (Prague) and MEDIA LARIO (Lecco). Both groups use roughly the same replication methods, depositing Nickel onto a gold coated mandrel surface. The separation occurs by using the different thermal expansion of the mandrel and replica materials. After that the replica is mounted into a cell of a similar thermal coefficient. The gold layer, that also acts as release agent, sticks to the replica and can be used as an optical surface. Alternatively it can be overcoated with other reflective layers.

# Replication trials

One major issue to address for assessing the suitability of this replication technique was to explore the interface conditions for the master form (mandrel). Usuallly for protruding (male) parts the vertical sides require an angle of at least 1$^{\circ}$ to enable releasing, while the edges require 30 $\mu$m chamfers. While both of these conditions are easy to meet for pupil and slit mirrors, for small slicing mirror components these requirements present problems, especially the chamfer which would lead to unacceptably low fill factors. To determine how these conditions could be optimized for small aperture parts, "replication challengers" were produced. These diamond machined aluminium test parts are 10 mm $\times$ 10 mm square, flat surfaces protruding out of a 25mm diameter cylindric blank. Their four edges have chamfers of different size (0, 10, 20 and 30 $\mu$m were specified). Each set consists of four challenger parts, featuring release angles of 1$^{\circ}$, 0.5$^{\circ}$, 0.2$^{\circ}$ and 0.1$^{\circ}$ while having all four chamfer widths each. The split of the lot into four parts was made to assure that at least some parts can be replicated if other fail. Two sets of the four pieces were manufactured by LFM, Bremen, by diamond fly-cutting. Each set was sent to one of the replication companies. Figure shows a challenger part in the different stages of replication.

# Mandrel wear

The mandrels wear down during replication, making it impossible to reach the desired surface quality after a while. To find out, how long a tool stands, a 100 $\mu$m thick NiP coated challenger part is currently in the process to obtain subsequent replicas from the same mandrel. Additionally REFLEX works with bare aluminium mandrels, using a special type of release agent to detach the parts after replication.

# Master material

The material commonly used is aluminium, which has to be Nickel-Phosphor plated to harden the surface to allow repeated replications. While the aluminium part is easy to produce using diamond machining techniques, the NiP coating tends to introduce surface errors if it is thicker than a few $\mu$m. Alternatively an about 100 $\mu$m thick layer can be applied before machining. Both methods introduce a further complication to make the surface hard enough for subsequent replication. Solutions are

**Zerodur mandrels:**

:   The hardness of Zerodur allows it to be used as a mandrel without the need for NiP coating. Furthermore, the use of Zerodur means that the mandrel-replica sandwich has to be heated up for separation, producing a higher temperature difference than for the standard separation procedure. To test these advantages, two optically polished Zerodur mandrels of the "challenger" type have been made to compare their capabilities with the aluminium mandrels.

**Release agent:**

:   REFLEX tested a proprietary agent that is deposited onto the mandrel before the gold coating is applied. Due to the improved ease of separating the parts, the mandrel wear is expected to be much lower than for using bare aluminium.

**Double replication:**

:   For optical replication the mandrel has to have a negative shape, which can be very difficult to achieve for image slicer components. To solve this problem, the feasibility of a double replication (positive original mandrel, negative replica as intermediate mandrel, and second replica positive again) has been tested by Reflex. Profilometer measurements of the mandrel and both replicas indicate that the deviations introduced by second replication are of the same order as for the single replication.

All three tests are under investigation: While the release agent has been tested at REFLEX already, the Zerodur mandrel replication is currently in progress. The double replication will be checked while replicating pupil mirror prototype parts, which can be machined and tested much easier as a positive shape.

# Realistically shaped arrays

Moving on from the abstract "challenger" part to a more realistic facetted surface, a 7-segment array has been designed similar to an all-spherical array used for the GNIRS-IFU pupil mirrors . As a positive shape, this array can be compared directly to the replica. This is produced by double replication, testing this process as well as the feasibility of elongated parts with stepped surfaces. First results of REFLEX indicate that double replication of slicing mirror arrays seems to work, albeit quantitative results are not available yet.

# First results

**Chamfers and angles:**

:   All chamfer and angle sizes could be reproduced easily, and all parts could be separated.

**Surface roughness:**

:   The surface roughness does not deteroriate very much, in the case of the release agent used by REFLEX this agent even seems to improve surface roughness, giving rise to smoother replicas.

**Surface shape:**

:   The shape is subject to a large scale bend, maximum deviations are 2 to 6 $\mu$m. Hence the optical form specifications are not met.

**Process duration:**

:   1 mm of Ni layer deposit takes about one week, causing long production cycles. Double replication may be a promising alternative.

While the first two results look promising, shape deformation and process duration can prohibit a use of this technique for the slicer application. Current experiments try to reduce the bending by slightly changing the chemistry of the process. Further experiments are planned to clarify if double replication can be used to reduce the process duration.

# Conclusions & future work

The first results above indicate that the replication of micro-optical surfaces is feasible, even in situations where no large chamfers can be applied. Albeit these results are encouraging, there remain several problems to be solved: Ways to improve the surface roughness and to suppress the surface deformation have to be found. Apart from that, the feasibility to replicate large numbers of pieces within an acceptable time scale has to be proven to exploit the advantages of this technique for the mass production of slicer micro-optics. Currently several ways of double replication are tested to reduce production time.

# Acknowledgements

This work is funded by the EU FP6 Opticon Smart Focal Planes consortium.


author: M. J. Griffin; A. Abergel ; A. Abreu ; P. A. R. Ade ; P. André ; J.-L. Augueres ; T. Babbedge ; Y. Bae ; T. Baillie ; J.-P. Baluteau ; M. J. Barlow ; G. Bendo ; D. Benielli ; J. J. Bock ; P. Bonhomme ; D. Brisbin ; C. Brockley-Blatt ; M. Caldwell ; C. Cara ; N. Castro-Rodriguez ; R. Cerulli ; P. Chanial ; S. Chen ; E. Clark ; D. L. Clements ; L. Clerc ; J. Coker ; D. Communal ; L. Conversi ; P. Cox ; D. Crumb ; C. Cunningham ; F. Daly ; G. R. Davis ; P. De Antoni ; J. Delderfield ; N. Devin ; A. Di Giorgio ; I. Didschuns ; K. Dohlen ; M. Donati ; A. Dowell ; C. D. Dowell ; L. Duband ; L. Dumaye ; R. J. Emery ; M. Ferlet ; D. Ferrand ; J. Fontignie ; M. Fox ; A. Franceschini ; M. Frerking ; T. Fulton ; J. Garcia ; R. Gastaud ; W. K. Gear ; J. Glenn ; A. Goizel ; D. K. Griffin ; T. Grundy ; S. Guest ; L. Guillemet ; P. C. Hargrave ; M. Harwit ; P. Hastings ; E. Hatziminaoglou ; M. Herman ; B. Hinde ; V. Hristov ; M. Huang ; P. Imhof ; K. J. Isaak ; U. Israelsson ; R. J. Ivison ; D. Jennings ; B. Kiernan ; K. J. King ; A. E. Lange$\dagger$ ; W. Latter ; G. Laurent ; P. Laurent ; S. J. Leeks ; E. Lellouch ; L. Levenson ; B. Li ; J. Li ; J. Lilienthal ; T. Lim ; J. Liu ; N. Lu ; S. Madden ; G. Mainetti ; P. Marliani ; D. McKay ; K. Mercier ; S. Molinari ; H. Morris ; H. Moseley ; J. Mulder ; M. Mur ; D. A. Naylor ; H. Nguyen ; B. O'Halloran ; S. Oliver ; G. Olofsson ; H.-G. Olofsson ; R. Orfei ; M. J. Page ; I. Pain ; P. Panuzzo ; A. Papageorgiou ; G. Parks ; P. Parr-Burman ; A. Pearce ; C. Pearson ; I. Pérez-Fournon; F. Pinsard ; G. Pisano ; J. Podosek ; M. Pohlen ; E. T. Polehampton ; D. Pouliquen ; D. Rigopoulou ; D. Rizzo ; I. G. Roseboom ; H. Roussel ; M. Rowan-Robinson ; B. Rownd ; P. Saraceno ; M. Sauvage ; R. Savage ; G. Savini ; E. Sawyer ; C. Scharmberg ; D. Schmitt ; N. Schneider ; B. Schulz ; A. Schwartz ; R. Shafer ; D. L. Shupe ; B. Sibthorpe ; S. Sidher ; A. Smith; A. J. Smith; D. Smith; L. Spencer ; B. Stobie ; R. Sudiwala ; K. Sukhatme ; C. Surace ; J. A. Stevens ; B. M. Swinyard ; M. Trichas ; T. Tourette ; H. Triou ; S. Tseng ; C. Tucker ; A. Turner ; M. Vaccari ; I. Valtchanov ; L. Vigroux ; E. Virique ; G. Voellmer ; H. Walker ; R. Ward ; T. Waskett ; M. Weilert ; R. Wesson ; G. J. White ; N. Whitehouse ; C. D. Wilson ; B. Winter ; A. L. Woodcraft ; G. S. Wright ; C. K. Xu ; A. Zavagno ; M. Zemcov ; L. Zhang ; E. Zonca
date: Received March 31 2010; accepted April 21 2010
title: The *Herschel*-SPIRE instrument and its in-flight performance [^1]

The Spectral and Photometric Imaging Receiver (SPIRE), is the *Herschel* Space Observatory's submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 $\mu$m, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 $\mu$m (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4 x 8, observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6 . The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the *Herschel* cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5 –- 2.

# Introduction

The SPIRE instrument is designed to exploit the particular advantages of the *Herschel* Space Observatory (Pilbratt et al. 2010) for observations in the submillimetre region: its large (3.5-m), cold ($\sim$ 85-K), low-emissivity ($\sim$ 1%) telescope; unrestricted access to the poorly explored 200-700 $\mu$m range; and the large amount (over 20,000 hrs) of observing time. In this paper we summarise the key design features of the instrument, outline its main observing modes, and present a summary of its measured in-flight performance and scientific capabilities. Further details of the instrument calibration are given in Swinyard et al. (2010).

# SPIRE instrument design

SPIRE consists of a three-band imaging photometer and an imaging Fourier Transform Spectrometer (FTS). The photometer carries out broad-band photometry ($\lambda$/$\Delta\lambda$ $\sim$ 3) in three spectral bands centred on approximately 250, 350 and 500 $\mu$m, and the FTS uses two overlapping bands to cover 194-671 $\mu$m (447-1550 GHz).

Figure 1 shows a block diagram of the instrument. The SPIRE focal plane unit (FPU) is approximately 700 x 400 x 400 mm in size and is supported from the 10-K *Herschel* optical bench by thermally insulating mounts. It contains the optics, detector arrays (three for the photometer, and two for the spectrometer), an internal $^{3}$He cooler to provide the required detector operating temperature of $\sim$ 0.3 K, filters, mechanisms, internal calibrators, and housekeeping thermometers. It has three temperature stages: the *Herschel* cryostat provides temperatures of 4.5 K and 1.7 K via high thermal conductance straps to the instrument, and the $^{3}$He cooler serves all five detector arrays.

The photometer and FTS both have cold pupil stops conjugate with the *Herschel* secondary mirror, which is the telescope system pupil, defining a 3.29-m diameter used portion of the primary. Conical feedhorns (Chattopadhyay et al. 2003) provide a roughly Gaussian illumination of the pupil, with an edge taper of around 8 dB in the case of the photometer. The same $^{3}$He cooler design (Duband et al. 2008) is used in SPIRE and in the PACS instrument (Poglitsch et al. 2010). It has two heater-controlled gas gap heat switches; thus one of its main features is the absence of any moving parts. Liquid confinement in zero-$g$ is achieved by a porous material that holds the liquid by capillary attraction. A Kevlar wire suspension system supports the cooler during launch whilst minimising the parasitic heat load. The cooler contains 6 STP litres of $^{3}$He, fits in a 200 x 100 x 100-mm envelope and has a mass of $\sim$ 1.7 kg. Copper straps connect the 0.3-K stage to the five detector arrays, and are held rigidly at various points by Kevlar support modules (Hargrave et al. 2006). The supports at the entries to the 1.7-K boxes are also light-tight.

All five detector arrays use hexagonally close-packed feedhorn-coupled spider-web Neutron Transmutation Doped (NTD) bolometers (Turner et al. 2001). The bolometers are AC-biased with frequency adjustable between 50 and 200 Hz, avoiding 1/*f* noise from the cold JFET readout readouts. There are three SPIRE warm electronics units: the Detector Control Unit (DCU) provides the bias and signal conditioning for the arrays and cold electronics, and demodulates and digitises the detector signals; the FPU Control Unit (FCU) controls the cooler and the mechanisms, and reads out all the FPU thermometers; and the Digital Processing Unit (DPU) runs the on-board software and interfaces with the spacecraft for commanding and telemetry.

# Photometer design

Figure 2 shows the opto-mechanical layout of the photometer side of the FPU. It is an all-reflective design (Dohlen et al. 2000) except for the dichroics used to direct the three bands onto the bolometer arrays, and the filters used to define the passbands (Ade et al. 2006). The image is diffraction-limited over the 4 x 8 field of view, which is offset by 11 from the centre of the *Herschel* telescope's highly curved focal surface. The input mirror M3, lying below the telescope focus, receives the *f*/8.7 telescope beam and forms an image of the secondary at the flat beam steering mirror (BSM), M4. A calibration source (Pisano et al. 2005) placed behind a hole in the centre of the BSM, is used to provide a repeatable signal for the bolometers. It occupies an area contained within the region of the pupil obscured by the hole in the primary. Mirror M5 converts the focal ratio to *f*/5 and provides an intermediate focus at M6, which re-images the M4 pupil to a cold stop. The input optics are common to the photometer and spectrometer and the separate spectrometer field of view is directed to the other side of the optical bench panel by a pick-off mirror close to M6. The 4.5-K optics are mounted on the SPIRE internal optical bench. Mirrors M7, M8 and a subsequent mirror inside the 1.7-K box form a one-to-one optical relay to bring the M6 focal plane to the detectors. The 1.7-K enclosure also contains the three detector arrays, and two dichroic beam splitters to direct the same field of view onto the arrays so that it can be observed simultaneously in the three bands.

The bolometer array modules are bolted to the outside wall of the 1.7-K box. Inside each one, the $^{3}$He stage, accommodating the detectors, feedhorns and final filters, is thermally isolated from the 1.7-K mount by tensioned Kevlar threads, and cooled by a thermal strap to the $^{3}$He cooler. The three arrays contain 43 (500 $\mu$m), 88 (350 $\mu$m) and 139 (250 $\mu$m) detectors. The relative merits of feedhorn-coupled detectors, as used by SPIRE, and filled array detectors, which are used by *Herschel*-PACS (Poglitsch et al. 2010) and some ground-based instruments such as SCUBA-2 (Audley et al. 2007) and SHARC-II (Dowell et al. 2003), are discussed in detail in Griffin et al. (2002). In the case of SPIRE, the feedhorn-coupled architecture was chosen as the best option given the achievable sensitivity, the requirements for the largest possible field of view and high stray light rejection, and limitations on the number of detectors imposed by spacecraft resource budgets. The detector feedhorns are designed for maximum aperture efficiency, requiring an entrance aperture close to 2$F\lambda$, where $\lambda$ is the wavelength and $F$ is the final optics focal ratio. This corresponds to a beam spacing on the sky of 2$\lambda/D$, where $D$ is the telescope diameter. The array layouts are shown schematically in Fig. 3 which also shows a photograph of an array module.

The photometric passbands are defined by quasi-optical edge filters (Ade et al. 2006) located at the instrument input, at the 1.7-K cold stop, and directly in front of the detector arrays, the reflection-transmission edges of the dichroics, and the cut-off wavelengths of the feedhorn output waveguides. The filters also serve to minimise the thermal loads on the 1.7-K and 0.3-K stages.

The AC-biased bolometer signals are de-modulated by individual lock-in amplifiers (LIAs) in the DCU. Each LIA comprises a bandpass filter and a square wave demodulator, followed by a 5-Hz low-pass filter. The output of the LIA is a DC voltage proportional to the RMS value of the voltage at the bolometer output. The LIA outputs are multiplexed and sampled at nominally 18.6 Hz for telemetry to the ground. In order to achieve the necessary 20-bit sampling using a 16-bit analogue-to-digital converter, an offset subtraction scheme is implemented. After multiplexing, a suitable 4-bit DC offset (generated on board and transmitted to the ground in the science data stream) is subtracted from each signal prior to the final gain stage before digitisation. The 130 kbs available data rate allows the data to be transmitted to the ground with no on-board processing.

# Spectrometer design

The FTS (Swinyard et al. 2003; Dohlen et al. 2000) uses two broadband intensity beam splitters in a Mach-Zehnder configuration which has spatially separated input and output ports. One input port views a 2.6 diameter field of view on the sky and the other an on-board reference source. Two bolometer arrays at the output ports cover overlapping bands of 194-313 $\mu$m (SSW) and 303-671 $\mu$m (SLW). As with any FTS, each scan of the moving mirror produces an interferogram in which the spectrum of the entire band is encoded with the spectral resolution corresponding to the maximum mirror travel.

The FTS focal plane layout is shown in Fig. 4. A single back-to-back scanning roof-top mirror serves both interferometer arms. It has a frictionless mechanism using double parallelogram linkage and flex pivots, and a Moiré fringe sensing system. A filtering scheme similar to the one used in the photometer restricts the passbands of the detector arrays at the two ports, defining the two overlapping FTS bands.

The FTS spectral resolution element is given by 1/(2$d$) where $d$ is the maximum optical path difference. The highest unapodised resolution available is 0.04 cm$^{-1}$ (1.2 GHz), corresponding to a FWHM of the instrument spectral response function of 0.048 cm$^{-1}$. For this resolution, $\lambda$/$\Delta$$\lambda$ varies between $\sim$ 1300 at the short-wavelength end and 370 at the long-wavelength end. The standard unapodised resolution adopted for spectrophotometry is 0.83 cm$^{-1}$ for which $\lambda$/$\Delta$$\lambda$ varies from approximately 20 to 60 between the long- and short-wavelength ends of the range.

The two hexagonally close-packed spectrometer arrays contain 37 detectors in the short-wavelength array and 19 in the long-wavelength array. The array modules are similar to those used for the photometer, with an identical interface to the 1.7-K enclosure. The feedhorn and detector cavity designs are optimised to provide good sensitivity across the whole wavelength range of the FTS. The SSW feedhorns are sized to give 2$F\lambda$ pixels at 225 $\mu$m and the SLW horns are 2$F\lambda$ at 389 $\mu$m. This arrangement has the advantage that there are many co-aligned pixels in the combined field of view. The SSW beams on the sky are 33 apart, and the SLW beams are separated by 51. A thermal source, SCal (Hargrave et al. 2003), at the second input port is available to allow the background power from the telescope to be nulled, thereby reducing the dynamic range needed in recording the interferograms (the amplitude of the interferogram central maximum is proportional to the difference in the radiant power from the two ports).

The FTS on-board signal chain is similar to that for the photometer, but has a higher low-pass-filter cut-off frequency of 25 Hz to accommodate the required signal frequency range, and the sampling rate is correspondingly higher (80 Hz). As with the photometer, the 130 kbs available data rate is sufficient to allow data transmission to Earth with no on-board processing.

# Observing modes

The photometer was designed to have three principal observing modes: point source photometry, field (jiggle) mapping, and scan mapping. The first two modes involve chopping and jiggling with the SPIRE beam steering mirror and nodding using the telescope. In practice, because of the excellent performance and simplicity of the scan-map mode, it is the optimum for most observations. Although its efficiency for point source and small map observations is low, it provides better data quality, and over a larger area, than the chopped modes, and can produce a sky background confusion limited map in a total observing time that is still dominated by telescope slewing overheads.

In scan-map mode, the telescope is scanned across the sky at 30 or 60/s. The scan angle is chosen to give the beam overlap needed for full spatial sampling over a strip defined by one scan line, and to provide a uniform distribution of integration time over the area covered by the scan. One map repeat is normally constituted by two scans in orthogonal directions to provide additional data redundancy and cross-linking. As well as SPIRE-only and PACS-only scan-map modes, *Herschel* can implement simultaneous five-band photometric scan mapping in SPIRE-PACS Parallel mode. Data are taken simultaneously in the three SPIRE bands and two of the PACS bands, providing a highly efficient means of making multi-band maps of large areas (at least 30 x 30). Scan speeds of 20 or 60/s are available in Parallel mode. The SPIRE detector sampling rate in this mode is reduced to 10 Hz to keep within the data rate budget, with minimal impact on the data quality.

The standard operating mode for the FTS is to scan the mirror at constant speed (nominally 0.5 mm/s) with the telescope pointing fixed, providing an optical path rate of change of 2 mm/s due to the factor of four folding in the optics. Radiation frequencies of interest are encoded by the scanning motion as detector output electrical frequencies in the range 3-10 Hz. The maximum scan length is 3.5 cm, corresponding to an optical path difference of 14 cm. Three values of spectral resolution are adopted as standard: high, medium, and low, giving unapodised resolutions of 0.04, 0.24, and 0.83 cm$^{-1}$ respectively (corresponding to 1.2, 7.2, and 25 GHz).

For point source observations with the FTS, the source is positioned at the overlapping SSW and SLW detectors at the array centres; but data are acquired for all of the detectors, providing at the same time a sparsely-sampled map of the emission from the region around the source. Likewise, a single pointing generates a sparse map of an extended object. Full and intermediate spatial sampling can be obtained using the SPIRE BSM to provide the necessary pointing changes between scans while the telescope pointing remains fixed. Regions larger than the field of view are mapped using telescope rastering combined with spatial sampling with the BSM. Full details of the SPIRE observing modes are given in the SPIRE Observers' Manual (2010).

# Instrument in-flight optimisation and achieved performance

## In-flight functional performance

All SPIRE subsystems are fully functional and meet their key requirements. The $^{3}$He cooler has a hold-time in excess of 46 hrs, allowing two full days of SPIRE operation after a cooler cycle. The achieved temperatures are $\sim$ 287 mK at the $^{3}$He cold tip, 310 mK for the photometer arrays and 315 mK for the FTS arrays. These values are in line with on-ground measurements and slightly better than used in the pre-launch instrument sensitivity model. Detector yield is excellent, with totals of only six unusable photometer bolometers and two unusable FTS bolometers. In-flight testing, carried out before the *Herschel* cryostat lid was removed, allowed the main results of system-level pre-flight testing (with *Herschel* in the Large Space Simulator facility at ESTEC) to be reproduced. The radiant background on the detectors, from the thermal emission of the *Herschel* telescope and any additional stray light, is consistent with on-ground characterisation of the *Herschel* telescope emissivity (Fischer et al. 2004), with a stray light level somewhat less than that adopted in pre-launch estimation of the performance. It is clear that the launch of *Herschel* proved to be a rather benign event for the system.

## Instrument optimisation and AOT characterisation

SPIRE has many parameters which are adjustable in flight to achieve best performance. These include: cooler recycle settings and detailed timing; JFET pre-amplifier supply voltage settings; detector array bias frequencies and voltages; detector lock-in amplifier phase settings; PCal and SCal applied power settings; telescope angular scan speeds and scan angles for scan-map observations; FTS mirror and BSM servo system parameters; and the FTS mirror scan speed. A good part of the Commissioning and Performance Verification (PV) phases of the *Herschel* mission involved systematic measurements and analysis to optimise the values of these parameters. In the case of the detector-related settings, separate optimisations were done for nominal source strength ($<$ 200 Jy for the photometer) and bright source settings. Most parameter optimisations resulted in values close to or identical to those predicted before launch. Important differences were: (i) the necessary frequency of PCal operation to track detector responsivity, planned for at least once per hour pre-launch, has been relaxed to once per observation regardless of the observation length; (ii) because of the low telescope and stray light background, it is possible to operate the FTS without switching on the SCal unit. This has the advantages that the operation and data reduction are simplified, and that the absence of SCal thermal emission reduces photon noise.

## Glitches due to ionising radiation

In space, ionising radiation hits deposit thermal energy in the bolometers causing spikes in the output voltage timelines. While these glitches are a nuisance, the effects for bolometric detectors are not as severe as in the case of low-background photoconductive detectors, which are prone to non-linear and memory effects as a result of charge build-up in the semiconductor material. Photoconductors can have operating currents of hundreds of electrons/s or less, but the much higher bias current of a typical cryogenic bolometer ($\sim$ 1 nA or $\sim$ 10$^{9}$ electrons/s) results in the deposited charge being swept out very rapidly, so that after the initial spike has decayed (on a timescale determined by the bolometer thermal time constant) the bolometer has no memory of the event, and its quiescent properties (responsivity and noise) are unaffected by its history of exposure to ionising radiation.

Two types of glitches are seen in the SPIRE detector timelines: large events associated with single direct hits on individual detectors, and smaller co-occurring glitches, seen simultaneously on many detectors in a given array, that are believed to be due to ionising hits on the silicon substrate that supports all the detectors in an array (frame hits). The glitch decay time constant is dictated by the bolometer thermal time constant (typically 6 ms), but the on-board analogue signal chains have low-pass electrical filters (designed to eliminate aliased high frequency noise from the sampled timelines), which have the effect of prolonging the duration of a glitch and reducing its amplitude. Because the photometer channels have a lower cut-off frequency (5 Hz) than the spectrometer channels (25 Hz), spectrometer glitches are detected more easily above the detector noise and so have a higher rate. With the currently adopted glitch recognition scheme, the detected rates per photometer detector are (0.6, 0.7, 1.5) glitches/detector/minute for (250, 350, 500) $\mu$m, leading to a fractional data loss of less than 1% in all cases. For the FTS, the equivalent numbers are $\sim$ 3 and 2/detector/minute/ for SSW and SLW respectively. At present, glitches are identified and flagged in the automatic processing pipelines. The corresponding data samples are rejected in the photometer map-making stage. In the FTS pipeline, the flagged sections of data are replaced by linear interpolation to ensure that artefacts are not introduced by subsequent Fourier processing in the pipeline. Work continues to improve the deglitching in both pipelines, but already glitches are not imposing a significant sensitivity penalty.

## Photometer performance

Beam profiles: The photometer beam profiles have been measured using scan maps of Neptune, which provides high S/N and appears point-like in all bands (angular diameter $\approx$ 2). The (250, 350, 500)  $\mu$m beams are well described by 2-D Gaussians down to $\sim$ 15 dB, with mean FWHM values of (18.1, 25.2, 36.6) and mean ellipticities of (7%, 12%, 9%). These values, which are very similar to pre-launch predictions, represent an average over each array; there are small systematic variations at the level of $\sim$ 5% across the arrays. Beam maps are available (via the ESA *Herschel* Science Centre) allowing the detailed beam shape to be deconvolved from map data if desired. Each beam map constitutes an averaging in the map over all of the individual bolometers crossing the source, and represents the realistic point source response function of the system, including all scanning artefacts or astrometric uncertainties.

Flux calibration: The primary calibrator for the photometer is Neptune, for which we adopt the model of Moreno (1998, 2010), which has an estimated absolute uncertainty of 5%. The calibration scheme is defined in detail in Swinyard et al. (2010), Griffin (2010), and in the SPIRE Observers' Manual (2010). Neptune was not observable during most of the PV phase, and an initial calibration was established based on observations of the asteroid Ceres. At the time of writing, the Neptune-based calibration has yet to be incorporated into the automatic data pipeline, but interim cross-checks have been made to assess the accuracy of the Ceres-based and Neptune-based calibrations, and the overall total uncertainty of the current calibration is estimated as within 15%. Further improvement on this figure is expected to be made once the Neptune-based calibration is implemented and through incremental refinements during the course of the mission. Flux densities are quoted at wavelengths of 250, 350 and 500 $\mu$m, based on the *Herschel* convention of an assumed source SED with flat $\nu$S($\nu$) (i.e., spectral index -1).  Colour correction factors to convert to a different spectral index are small (a few percent) and within the current 15% overall uncertainty.

Sensitivity: The photometer sensitivity has been estimated from repeated scan maps of dark regions of extragalactic sky. A single map repeat is constituted by two orthogonal scans as implemented in the SPIRE-only scan-map AOT. Multiple repeats produce a map dominated by the fixed-pattern sky confusion noise, with the instrument noise having integrated down to a negligible value. This sky map can then be subtracted from individual repeats to estimate the instrument noise. The extragalactic confusion noise levels for SPIRE are assessed in detail by Nguyen et al. (2010), who define confusion noise as the standard deviation of the flux density in the map in the limit of zero instrument noise. Measured confusion noise levels in the (250, 350, 500) $\mu$m bands are (5.8, 6.3, 6.8) mJy. For the nominal scan speed of 30/s, the instrument noise is estimated at (9.0, 7.5, 10.8) mJy/sqrt($N_{reps}$) where $N_{reps}$ is the number of map repeats. These values are comparable to pre-launch estimates at 250 and 500 $\mu$m, and better by a factor of $\sim$ 1.7 at 350 $\mu$m. For the fast scan speed, the instrument noise levels correspond to (12.7, 10.6, 15.3) mJy/sqrt($N_{reps}$), precisely as expected given the factor of two reduction in integration time per repeat. For the nominal scan speed, the overall noise is within a factor of sqrt(2) of the (250, 350, 500) $\mu$m confusion levels for (3, 2, 2) repeats.

An important aspect of the photometer noise performance is the knee frequency that characterises the 1/*f* noise of the detector channels. Pre-launch, a requirement of 100 mHz with a goal of 30 mHz had been specified. In flight, the major contributor to low frequency noise is temperature drift of the $^{3}$He cooler. Active control of this temperature is available via a heater-thermometer PID control system, but has not yet been used in standard AOT operation. The scan-map pipeline (see Sect. 7) includes a temperature drift correction using thermometers, located on each of the arrays, which are not sensitive to the sky signal but track the thermal drifts. This correction works well and will be improved with a forthcoming update of the flux calibration parameters. Use of the thermometer signals to de-correlate thermal drifts in the detector timelines over a complete observation (E. Pascale, priv. comm.) can produce a 1/$f$ knee of as low as a 1-3 mHz. This corresponds to a spatial scale of several degrees at the nominal scan speed.

Observing overheads: Scan-map mode is very efficient for large area observations, but somewhat less so for smaller fields due to the time needed to turn the telescope around at the end of each scan leg. The HSpot observation preparation tool (version 4.4.4 at the time of writing) provides a detailed summary of the on-source integration time and the various overheads for any particular observation. For example, a 2$^\circ$ x 2$^\circ$ map with two repeats takes a total of 7.4 hrs (including all instrument and telescope overheads) with an the overall efficiency of 86%. For a 20 x 20 map with two repeats, the duration is 31 min. with an overall efficiency of 51% due to the larger fraction of time taken up by telescope turn-arounds. For the Small Map AOT ($\sim$ 5 x 5 map area), with two repeats, the duration is $\sim$ 8 min. with an efficiency of 15% (but the total observing time still has a large relative contribution from the 3-min. telescope slewing overhead).

The photometer mapping AOTs are designed to ensure that the area requested by the observer is observed with full sampling and highly uniform coverage at the nominal scan speeds. Maps also contain a region around the target area corresponding to the telescope turn-around periods, in which data are available to the user, but with a lower than nominal scan speed and with less complete spatial sampling.

## Spectrometer performance

Wavelength coverage and calibration: The wavelength coverage of the spectrometer is as described in Sect. 4, with the short wavelength SSW band covering 32.0-51.5 cm$^{-1}$ (194-313 $\mu$m) and the long wavelength SLW band covering 14.9-33.0 cm$^{-1}$ (303-671 $\mu$m). Wavelength calibration, verified using CO lines in galactic sources, is accurate within 1/10 of a high-resolution spectral resolution element across both bands.

Beam profiles: The beamwidths vary across the FTS bands, and this has been characterised using spectral mapping of Neptune as a suitable strong point source. The FWHM varies between 17 and 21 between the short- and long-wavelength ends of the SSW band. The SLW FWHM exhibits a dip within the band, with values of 37 at the short-wavelength end, 42 at the long-wavelength end, and a minimum of 29 at $\sim$ 425 $\mu$m. More details are given in Swinyard et al. (2010).

Flux calibration: The spectrometer calibration scheme and methods are described in detail by Swinyard et al. (2010). At the time of writing, the FTS flux calibration is still based on the asteroid Vesta, but will be updated in the near future based on observations of Neptune and Uranus using the planetary atmosphere models of Moreno (1998, 2010). Current overall calibration accuracy is estimated at 15-30%, depending on position in the bands, for frequencies above 20 cm$^{-1}$, and will be improved when the planet-based calibration is implemented. Correct subtraction of the thermal emission from the telescope (and from the $\sim$ 5 K instrument enclosure at the lower frequencies) currently requires expert interactive analysis, especially for faint sources, and work continues to establish an automatic pipeline calibration. The current FTS calibration is adequate for point or extended sources, but for sources partly extended with respect to the beam, calibration is difficult because it requires some a priori knowledge or modelling of the source intensity distribution.

Sensitivity: The line sensitivity (high-resolution mode; 5 $\sigma$; 1 hr) of the FTS, achieved to date, is typically 1.5 x 10$^{-17}$ W m$^{-2}$ for the SSW band and 2.0 x 10$^{-17}$ W m$^{-2}$ for the SLW band. The sensitivity varies slightly across the band in the manner expected from the shape of the Relative Spectral Responsivity Function (RSRF), as given in the SPIRE Observer's Manual (2010). At the time of writing, the noise level continues to integrate down as expected for at least 20 repeats (20 forward and 20 reverse scans; total on-source integration time of $\sim$ 45 min.). The FTS sensitivities are better than pre-launch (HSpot) predictions by a factor of 1.5 – 2. This is attributable to the low telescope background, the fact that the SCal source is not used, and to a conservatism factor that was applied to the modelled sensitivities to account for various uncertainties in the model.

Observing overheads: The telescope pointing is static for most of the time during typical FTS observations, so the overheads are very low. For example, a single pointing high-resolution observation with 20 repeats takes approximately one hour and has an overall efficiency of 86%. A fully sampled spatial map, requiring 16 jiggle positions of the BSM, with four repeats per position, takes about 2.4 hrs and is 91% efficient.

# Data-processing pipelines and data quality

The architecture and operation of the photometer pipeline is described by Griffin et al. (2008), and the spectrometer pipeline is described by Fulton et al. (2008).

Photometer data are processed in a fully automatic manner, producing measured flux density timelines (Level-1 products) and maps (Level-2 products). The pipeline includes conversion of telemetry packets into data timelines, calculation of the bolometer voltages from the raw telemetry, association of a sky position for each detector sample, glitch identification and flagging, corrections for various effects including time constants associated with the detectors and electronics, conversion from detector voltage to flux density, and correction for detector temperature drifts. The Level-1 products are calibrated timelines suitable for map-making. Maps can be made either with the maximum likelihood map-making algorithm MADmap (Cantalupo et al. 2010) or by "naïve" mapping, involving simple binning of the measured flux densities. The pipeline can also be run in an interactive mode, with selectable parameters associated with deglitching, baseline removal and map pixel size. Standard map pixel sizes of (6, 10, 14) are adopted for the (250, 350, 500) $\mu$m bands. The most important aspect of the Level-1 data quality that must be addressed prior to map-making is the amount of residual thermal baseline drift on the timelines. At the time of writing, the standard data products are naïve maps generated with pre-treatment of the Level-1 timelines to remove these residual drifts. The results are high-quality maps, with a low level of scan-related artefacts and with sky structures preserved over large spatial scales (degrees). Work continues to improve further the map-making process, in particular to automate the Level-1 timeline treatment step prior to creation of the maps. Once that has been completed, the potential advantages of MADmap in producing further enhancements will be investigated.

The FTS pipeline processes the telemetry data producing calibrated spectra. It shares some elements with the photometer pipeline, including the conversion of telemetry into data timelines and the calculation of bolometer voltages. Steps unique to the spectrometer are: temporal and spatial interpolation of the stage mechanism and detector data to create interferograms, apodisation, Fourier transformation, and creation of a hyperspectral data cube. Corrections are made for various instrumental effects including first- and second-level glitch identification and removal, interferogram baseline correction, temporal and spatial phase correction, non-linear response of the bolometers, variation of instrument performance across the focal plane arrays, and variation of spectral efficiency.

# Conclusions

The SPIRE instrument is fully functional with performance and scientific capabilities matching or exceeding pre-launch estimates. Flux calibration is currently accurate to 15% for the photometer and 15-30% for the FTS, and the current pipelines are already producing high-quality data. Major changes to the AOT implementation are not needed or planned, but further enhancements to the calibration and pipelines will continue as the mission progresses.

[^1]: *Herschel* is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.


abstract: White dwarfs represent the endpoint of stellar evolution for stars with initial masses between approximately 0.07 $M_{\odot}$ and 8-10 $M_{\odot}$, where $M_{\odot}$ is the mass of the Sun (more massive stars end their life as either black holes or neutron stars). The theory of stellar evolution predicts that the majority of white dwarfs have a core made of carbon and oxygen, which itself is surrounded by a helium layer and, for $\sim$`<!-- -->`{=html}80 per cent of known white dwarfs, by an additional hydrogen layer$^{1-3}$. All white dwarfs therefore have been traditionally found to belong to one of two categories: those with a hydrogen-rich atmosphere (the DA spectral type) and those with a helium-rich atmosphere (the non-DAs). Here we report the discovery of several white dwarfs with atmospheres primarily composed of carbon, with little or no trace of hydrogen or helium. Our analysis shows that the atmospheric parameters found for these stars do not fit satisfactorily in any of the currently known theories of post-asymptotic giant branch evolution, although these objects might be the cooler counterpart of the unique and extensively studied PG1159 star H1504+65 (refs 4-7). These stars, together with H1504+65, might accordingly form a new evolutionary sequence that follow the asymptotic giant branch.
author: P. Dufour, James Liebert, G. Fontaine, and N. Behara
title: Rare White dwarf stars with carbon atmospheres

Traces of carbon are typically observed as either neutral carbon lines or molecular C$_2$ Swan bands (defining the DQ spectral type) in cool helium-rich white dwarfs with effective temperatures ($T_{\rm eff}$) below $\sim$ 13,000 K. The presence of carbon in the atmospheres of these objects has been explained successfully by a model in which carbon is dredged-up from the underlying carbon/oxygen core by the deep helium convection zone$^8$. This model predicts a maximum contamination of carbon at an effective temperature of $\sim$`<!-- -->`{=html}10,000 K (corresponding approximately to the temperature at which the surface convection zone is maximal) before gradually decreasing with lower temperature, in agreement with atmospheric analysis determinations$^9$. We note that although some of these objects show a very high level of carbon pollution, helium always remains the dominant constituent of the atmosphere (the highest carbon abundances by number (N) found are around logN(C)/N(He) $\approx -3$).

Atmospheric abundance determination for DQ white dwarfs has been most successful for stars on the cool side of the maximum contamination and very few analyses have been done for hotter objects located on the ascending side of the curve$^{10-12}$. Recently$^{13}$, several new hot DQ white dwarfs have been discovered in the Sloan Digital Sky Survey (SDSS), providing a unique opportunity to test the dredge-up theory on the "hot side". However, appropriate atmospheric models were not available at that time. Because the atmospheric compositions of these white dwarfs were expected to be helium-dominated, effective temperature estimates (T$_{He}$) were obtained by comparing the SDSS colours with those of pure helium composition models. Although it was expected that such an approximation would overestimate the effective temperature by several thousand kelvin, some of these stars were so hot (T$_{He} \ge$ 25,000 K) that an overlap with most of the helium-rich star (the DB spectral type) temperature range seemed unavoidable. It was hypothesized that these stars had thinner outer helium envelopes so that dredge-up occurred earlier in the cooling sequence. These highly carbon-polluted white dwarfs are expected to be massive, so it was argued that they might represent the missing high-mass tail of the DB mass distribution$^{14}$.

Thus, it is with this scientific rationale in mind that we proceeded with the calculation of the appropriate atmospheric models for these objects. Because the continuum opacity of heavy elements might not be negligible for these objects, these new models have been updated with the latest C and O photoionization cross-sections from the Opacity Project$^{15}$. Although the analysis of the coolest objects ($T_{\rm eff}\le$`<!-- -->`{=html}15,000 K) is straightforward (results will be presented elsewhere; manuscript in preparation), we quickly realized that no combination of carbon and helium could successfully reproduce the observed features (mostly Cii lines) in the optical spectrum of the hottest ones by assuming a helium-dominated atmosphere. Indeed, such models predict the presence of a strong Hei $\lambda$=4,471 line that is not observed spectroscopically in our sample of hot DQ stars. We thus concluded that a good fit to both the spectra and the energy distribution was possible only by considering atmospheres made primarily of carbon, with little or no trace of hydrogen and helium.

Figure 1 shows our fits to the optical spectrum and photometric energy distribution of one such star using a grid of pure carbon atmosphere models, as well as models containing traces of hydrogen and helium. The fitting method is similar to that used for cooler DQ and needs not to be repeated here$^9$. This particular case shows a small trace of hydrogen which we believe is more likely to be the result of accretion from the interstellar medium than of primordial origin. We found eight more similar objects in the SDSS white dwarf catalogue$^{16}$, all of which were found to have a carbon-dominated atmosphere and a temperature between 18,000 and 23,000 K, although for some faint stars, detailed analysis is more uncertain because of the poor signal-to-noise ratio of the observations. Only crude upper limits for hydrogen and helium (of the order of logN(C)/N(H) $\ge$ 1.5 and logN(C)/N(He) $\ge$ 1.5) can be obtained from the other spectra. Our detailed model atmosphere fits provide fair estimates of the effective temperature and chemical composition, but higher-quality data (which we hope to obtain in December 2007, weather permitting) are essential for a precise measure of the surface gravity (and thus mass) of these objects.

These new stars are too hot to be explained by the standard convective dredge-up scenario and so we sought another explanation. Natural progenitor candidates to consider are the hot PG1159 stars. The latter are probably the result of a late He-shell flash$^{17}$ (the so-called "born-again" scenario) at the end of the post-asymptotic-giant-branch phase, which almost completely eliminates the remaining hydrogen and mixes material from the interior with the helium envelope. As a result, these objects re-enter the white dwarf cooling track, but this time they have a surface composition that is a mixture of helium, carbon, oxygen, and little or no hydrogen ($T_{\rm eff}$ between $\sim$`<!-- -->`{=html}75,000 and 200,000 K; typical abundances$^{18}$, in mass fraction, of He, C, O and Ne are 33$\%$, 50$\%$, 15$\%$ and 2$\%$). However, as these stars cool down, gravitational diffusion rapidly separates the heavier elements from the helium that tends to float to the surface and by the time they have reached 25,000 K, helium completely dominates the surface composition.

Models exploring the evolution of stars 9-11 $M_{\odot}$ have produced white dwarfs with O-Ne-Mg cores and CO envelopes and possibly little He and H; such models may explain our peculiar stars$^{19-21}$. However, it is unclear what quantity of He and H should remain for such stars. Furthermore, the resulting white dwarfs have high surface gravities ($\log g\approx 9$, in centimetre-gram-second units) and would thus show extremely broad carbon lines that are incompatible with the observed widths of the lines (Fig. 1 shows that $\log g
\approx 8$ is more likely, although analysis of a better spectrum is needed for more precision). The most likely scenario to explain the very existence of hot DQ white dwarfs with carbon-dominated atmospheres is that these stars are the progenies of objects such as H1504+65. The latter is a unique object among the known PG1159 stars so far. It is the hottest specimen of its class at $T_{\rm eff}\approx$ 200,000 K, and its atmospheric composition very unusual, with a mass fraction of $\sim 50\%$ C and $\sim 50\%$ O plus small traces of heavier elements, but no detectable helium or hydrogen. It is currently believed that H1504+65 is essentially a bare stellar nucleus produced by a particularly violent post-asymptotic-giant-branch very late thermal pulse that has destroyed in large part the remaining stellar envelope containing helium and hydrogen. The best available simulations of such an event, those of ref. 17, suggest indeed that, following the last He flash, the full envelope becomes convective and extends to deep enough layers for H and He to be completely consumed, except that very small traces of He may survive.

The uniformity of the chemical composition of the envelope of H1504+65 is believed to be maintained for a while by a residual wind but, with time, diffusion becomes dominant and He, C and O must separate under the influence of gravitational settling. The idea is that the wind slowly dies out with cooling (decreasing luminosity) and that the braking effect it provides on element separation becomes less and less efficient. Hence, the small residual trace of helium believed to exist in the envelope of H1504+65 eventually diffuses upward to form a thin layer above a C-enriched and O-depleted mantle. The total mass of an atmosphere is tiny ($\sim 10^{-14}$ to $10^{-15}$ M$_*$, where M$_*$ is the mass of the star), so there is ultimately enough accumulated He to form a full atmosphere and the descendant of H1504+65 would now be "disguised" as a He-atmosphere white dwarf after the PG1159 evolutionary phase. With further cooling, a convection zone develops in the C-enriched mantle owing to the recombination of that element while the overlaying He layer remains radiative. At some $T_{\rm eff}$, whose exact value is not yet known due to a current lack of proper models, it is conjectured that the subphotospheric C convection zone becomes active enough to be able to dilute from below the overlaying He layer. At that point the star would undergo a dramatic spectral change, transforming itself from a He-dominated atmosphere white dwarf to a star with a C-dominated atmosphere, because the mass in the C convection zone is orders of magnitude larger than the mass of the He layer.

Hence, the former PG1159 star H1504+65, showing initially a mixed C and O atmosphere, would now show a C-dominated atmosphere after an intermediate phase in which it would have been observed as a He-atmosphere white dwarf. We note, however, that helium must ultimately reappear at the surface in enough quantity to form a helium-rich atmosphere because no carbon-rich object has ever been discovered at lower temperatures ($T_{\rm eff}\le$ 15,000 K). The remaining helium probably floats again at the surface when the star cools down, perhaps turning the carbon-rich objects into DQ stars belonging to the second sequence of high carbon abundances$^{11,22}$. We are at present making evolutionary calculations to test the above scenario.

To conclude, we roughly estimate that the space density of carbon-rich white dwarfs, assuming a 50$\%$ completeness for the targeting in SDSS$^{16}$, is between $\sim 2.1\times10^{-6}$ pc$^{-3}$ and $\sim
7.0\times10^{-8}$ pc$^{-3}$ depending on the preferred value of $\log g$ (and thus radius) and $T_{\rm eff}$ for the stars. For comparison, the local space density of white dwarfs$^{23}$ is estimated to be $\sim
5.0\times10^{-3}$ pc$^{-3}$. White dwarfs with carbon dominated atmospheres are thus intrinsically rare and a relatively large volume of space (as was done by the Sloan Digital Sky Survey) had to be surveyed before many of them could be found.

. Iben, I., Jr. On the frequency of planetary nebula nuclei powered by helium burning and on the frequency of white dwarfs with hydrogen-deficient atmospheres. Astrophys. J. **277**, 333-354 (1984).
2. Koester, D. $\&$ Schoenberner, D. Evolution of white dwarfs. Astron. Astrophys. **154**, 125-134 (1986).
3. D'Antona, F. $\&$ Mazzitelli, I. Evolutionary times of white dwarfs-Long or short? IAU Colloq. **95**, 635-637 (1987).
4. Nousek, J. A. et al. H1504+65-an extraordinarily hot compact star devoid of hydrogen and helium. Astrophys. J. **309**, 230-240 (1986).
5. Werner, K. NLTE analysis of the unique pre-white dwarf H 1504 + 65. Astron. Astrophys. **251**, 147-160 (1991).
6. Werner, K. $\&$ Wolff, B. The EUV spectrum of the unique bare stellar core H1504+65. Astron. Astrophys. **347**, L9-L13 (1999).
7. Werner, K., Rauch, T., Barstow, M. A. $\&$ Kruk, J. W. Chandra and FUSE spectroscopy of the hot bare stellar core H 1504+65. Astron. Astrophys. **421**, 1169-1183 (2004).
8. Pelletier, C., Fontaine, G., Wesemael, F., Michaud, G. $\&$ Wegner, G. Carbon pollution in helium-rich white dwarf atmospheres. Time-dependent calculations of the dredge-up process. Astrophys. J. **307**, 242-252 (1986).
9. Dufour, P., Bergeron, P. $\&$ Fontaine, G. Detailed spectroscopic and photometric analysis of DQ white dwarfs. Astrophys. J. **627**, 404-417 (2005).
10. Wegner, G. $\&$ Koester, D. Atmospheric analysis of the carbon white dwarf G227-5. Astrophys. J. **288**, 746-750 (1985).
11. Thejll, P., Shipman, H. L., MacDonald, J. $\&$ Macfarland, W. M. An atmospheric analysis of the carbon-rich white dwarf G35 - 26. Astrophys. J. **361**, 197-206 (1990).
12. Desharnais, S., Wesemael, F., Chayer, P. $\&$ Kruk, J. W. FUSE observation of cool DB white dwarfs. ASP Conf. Ser. **372**, 265-268 (2007).
13. Liebert, J. et al. SDSS white dwarfs with spectra showing atomic oxygen and/or carbon lines. Astron. J. **126**, 2521-2528 (2003).
14. Beauchamp, A. Détermination des Paramètres Atmosphériques des étoiles Naines Blanches de Type DB. PhD thesis, Montréal (1995).
15. Behara, N. $\&$ Jeffery, C. S. LTE model atmosphere with new opacities. 1. Methods and general properties. Astron. Astrophys. **451**, 643-650 (2006).
16. Eisenstein, D. J. et al. A catalog of spectroscopically confirmed white dwarfs from the Sloan Digital Sky Survey data release 4. Astrophys. J. **167** (Suppl.), 40-58 (2006).
17. Herwig, F., Blöcker, T., Langer, N. $\&$ Driebe, T. On the formation of hydrogen-deficient post-AGB stars. Astron. Astrophys. **349**, L5-L8 (1999).
18. Werner, K. $\&$ Herwig, F. The elemental abundances in bare planetary nebula central stars and the shell burning in AGB stars. Publ. Astron. Soc. Pacif. **118**, 183-204 (2006).
19. Garcia-Berro, E. $\&$ Iben, I. On the formation and evolution of super-asymptotic giant branch stars with cores processed by carbon burning. 1: SPICA to Antares. Astrophys. J. **434**, 306-318 (1994).
20. Garcia-Berro, E., Ritossa, C. $\&$ Iben, I. On the evolution of stars that form electron-degenerate cores processed by carbon burning. III. The inward propagation of a carbon-burning flame and other properties of a 9 $M_{\odot}$ model star. Astrophys. J. **485**, 765-784 (1997).
21. Ritossa, C., Garcia-Berro, E. $\&$ Iben, I. On the evolution of stars that form electron-degenerate cores processed by carbon burning. V. Shell convection sustained by helium burning, transient neon burning, dredge-out, URCA cooling, and other properties of an 11 $M_{\odot}$ population I model star. Astrophys. J. **515**, 381-397 (1999).
22. Koester, D. $\&$ Knist, S. New DQ white dwarfs in the Sloan Digital Sky Survey DR4: confirmation of two sequences. Astron. Astrophys. **454**, 951-956 (2006).
23. Holberg, J. B., Oswalt, T. D. $\&$ Sion, E. M. A determination of the local density of white dwarf stars. Astrophys. J. **571**, 512-518 (2002).
This work has been partially supported by the NSF for work on SDSS white dwarfs. This work was also supported in part by the NSERC (Canada).

**Author Information** Reprints and permissions information is available at www.nature.com/reprints. Correspondence and requests for materials should be addressed to P.D. (email@example.com).