PIONIER: a 4-telescope visitor instrument at VLTI††thanks: Based on observations collected at the European Southern Observatory, Paranal, Chile (commissioning data and 087.C-0709).
Key Words.:Techniques: high angular resolution - Techniques: interferometric - Instrumentation: high angular resolution - Instrumentation: interferometers
Context:PIONIER stands for Precision Integrated-Optics Near-infrared Imaging ExpeRiment. It combines four m Auxilliary Telescopes or four m Unit Telescopes of the Very Large Telescope Interferometer (ESO, Chile) using an integrated optics combiner. The instrument has been integrated at IPAG starting in December 2009 and commissioned at the Paranal Observatory in October 2010. It provides scientific observations since November 2010.
Aims:In this paper, we detail the instrumental concept, we describe the standard operational modes and the data reduction strategy. We present the typical performance and discuss how to improve them.
Methods:This paper is based on laboratory data obtained during the integrations at IPAG, as well as on-sky data gathered during the commissioning at VLTI. We illustrate the imaging capability of PIONIER on the binaries Sco and HIP11231.
Results:PIONIER provides 6 visibilities and 3 independent closure phases in the H band, either in a broadband mode or with a low spectral dispersion (R=40), using natural light (i.e. unpolarized). The limiting magnitude is in dispersed mode under median atmospheric conditions (seeing, ms) with the m Auxiliary Telescopes. We demonstrate a precision of deg on the closure phases. The precision on the calibrated visibilities ranges from 3 to 15% depending on the atmospheric conditions.
Conclusions:PIONIER has been installed and successfully tested as a visitor instrument for the VLTI. It permits high angular resolution imaging studies at an unprecedented level of sensitivity. The successful combination of the four m Unit Telescopes in March 2011 demonstrates that VLTI is ready for 4-telescope operation.
The Very Large Telescope Interferometer (VLTI, Haguenauer et al. 2010) is equipped with four Unit Telescopes of m (UTs) and four relocatable Auxiliary Telescopes of m (ATs). It offers a unique combination of interferometric imaging capability and high sensitivity. A nice review of the scientific opportunities is being published by Berger et al. (2011, sub. to A&A Rev). It emphasizes the need to recombine a large number of telescopes to provide reliable snapshot imaging capabilities. However, the current instrumentation suite can only handle two or three telescopes simultaneously. The next generation of facility instruments that will combine four telescopes is not expected to be operational prior to 2014.
Therefore, IPAG111Institut de Planétologie et d’Astrophysique de Grenoble and its partners have proposed to use a visitor instrument222http://www.eso.org/sci/facilities/paranal/instruments/vlti-visitor that provides VLTI, since 2010, with a new observational capability combining imaging, sensitivity and precision. The principle of this Precision Integrated-Optics Near-infrared Imaging ExpeRiment (PIONIER) was approved by the ESO Science and Technical Commitee in spring 2009. The instrument has been integrated at IPAG starting in December 2009 and commissioned at the Paranal Observatory in October 2010. Only a few nights of commissioning were necessary before the instrument started routinely delivering scientific data with the ATs. On 17 March 2011 light collected by all four of the m UTs was successfully combined for the first time. A picture of PIONIER installed at the VLTI laboratory is shown in Fig. 1.
PIONIER relies on the integrated optics (IO) technology to combine four beams (Malbet et al. 1999; Berger et al. 2001; Benisty et al. 2009). In one observation, it provides the measurements of six visibilities and three independent closure phases with low spectral resolution () in the H band. In many respect, it is complementary to the current AMBER instrument that provides three visibilities and one closure phase with various spectral resolutions () in the H and K bands (Petrov et al. 2007). PIONIER paves the road for the forthcoming instruments GRAVITY and MATISSE (Gillessen et al. 2010; Lopez et al. 2008), by enabling the 4-telescope operation of VLTI.
The paper is organized as follow: section 2 presents the main astrophysical programs of PIONIER and the corresponding high-level specifications. Section 3 describes the instrument and its different subsystems. Section 4 presents the PIONIER operations, including regular day-time maintenance, flux injection, scientific fringe tracking and regular calibrations. Section 5 presents the data reduction strategy, focusing on the aspects that are specific to the instrument, or that have not been published elsewhere, and the current performance. Finally, Section 6 illustrates the imaging capability of PIONIER on the Be star Sco and the spectroscopic binary HIP11231. The paper ends with brief conclusions and remarks about possible evolutions of the instrument.
2 Science drivers and specifications
PIONIER is designed for aperture synthesis imaging with a specific emphasis on fast fringe recording to allow precision closure phases and visibilities to be measured. Among all the possible science cases, we defined three key astrophysical programs.
2.1 Young Stellar Objects
Young stars are surrounded by an accretion disk made of gas and dust, thought to be the birthplace of planets. Studying the innermost astronomical unit is of major importance to understand the physical conditions for star and planet formation.
Herbig Ae/Be stars are young stars of intermediate-mass, whose disks have been intensively studied with near-infrared interferometers in the last 13 years. Measurements enabled the determination of near-infrared sizes, found not to agree with classical accretion disk model predictions (see Millan-Gabet et al. (2007) for a review). Consequently the view on the inner structure of accretion disks has been substantially revised. The concept of “puffed-up” inner rim has been introduced (Dullemond et al. 2001; Isella & Natta 2005). It marks the boundary in the disk where dust can no longer exist due to sublimation and forms a barrier to the photons potentially causing an alteration of the vertical structure of the disk. Recently, strong evidence was found that the matter flowing through this limit to the star is a major contributor to the excess emission (Tannirkulam et al. 2008; Benisty et al. 2010). Revealing the nature of the matter inside the dust sublimation limit and the exact structure of the inner disk at this transition is now investigated with care. We have identified a small sample of HAe/Be stars whose environments have been clearly resolved by the VLTI and should therefore be ideal candidates for imaging with PIONIER in order to reveal the morphology of the inner astronomical unit near-infrared emission. The simultaneous combination of four telescopes is required to obtain reliable snapshot image reconstructions of these complex and time-variable objects (Benisty et al. 2011).
In addition, PIONIER will probe the structure of disks around the T Tauri stars, which are the low-mass analogs of Herbig stars. In these fainter objects, the inner edge of the dusty disk is closer to the star, and is not expected to be fully resolved, preventing us from direct imaging. However, the visibilities will give access to the inner disk location (Akeson et al. 2005; Eisner et al. 2009) and the closure phases will probe the morphology and surface brightness profile. In the best cases, it will be possible to constrain, through proper radiative transfer modeling, the inner rim position, ellipticity and thickness. All of these are of major importance to understand the conditions for star and planet formation. As T Tauri stars are the faintest objects of our key programs, they set the need for limiting magnitude of .
2.2 Debris disks around main sequence stars
The inner solar system contains a cloud of small dust grains created when small bodies –asteroids and comets– collide and outgas. This dust cloud, called the zodiacal disk, has long been suspected to have extrasolar analogs around main sequence stars. The very inner part of these exozodiacal disks have remained elusive until recently, when near-infrared interferometric observations revealed a resolved excess flux of about 1% around several nearby stars (Absil et al. 2006, 2008, 2009; di Folco et al. 2007; Akeson et al. 2009). These excesses have been interpreted as being due to hot, possibly transient exozodiacal dust within AU from these stars.
Our main goal is to survey main sequence stars to study the hot dust content during the late phases of terrestrial planet formation, and its connection with outer resevoirs of cold dust. We expect that about 100 main sequence stars of various ages between Myr and a few Gyr will be observed during the PIONIER lifetime, providing us with a solid statistical sample to derive the time dependence of the bright exozodi phenomenon. In addition, for the brightest cases, PIONIER could be used to characterize the grain properties thanks to multi-colour information, and to study the inner disk morphology. The challenge of this program is to obtain an accuracy of 1% or less on the calibrated visibility.
2.3 Binaries, faint companions and Hot Jupiters
Long baseline interferometric observations naturally complement the radial velocity and adaptive optics surveys to discover and characterize binary systems. For instance, radial velocity measurements on massive stars are challenged by the lack of spectral lines and their intrinsic broadening. On active stars, the intrinsic radial velocity jitter may easily hide the signal of a faint companion. Moreover, compared with adaptive optics, the access to smaller separations and therefore shorter periodes increases the possibility to determine the dynamical masses of the components.
Our first attempts with PIONIER (Absil et al. 2011, submitted) demonstrate that the angular separation range from to mas can be covered within one hour of observation, thanks to the simultaneous measurement of 4 closure phases at a low spectral resolution. Accordingly, we have initiated several surveys from massive to low mass stars, including young stars and stars in nearby moving groups. In order to reach a sufficient number of targets, these programs typically require a limiting magnitude of and a precision of on the closure phases (1:200 dynamic range).
Deep integration on few selected targets will permit to assess what is the best possible dynamic range achievable, until possibly reaching the planetary regime (1:2000). If successful PIONIER will have the means to separate the flux contribution of the planet from that of its parent star and potentially obtain low spectral resolution information. In order to achieve such a result, systematic biases on the closure phase estimate have to be hunted down and understood at the unprecedented level of deg (Zhao et al. 2010, 2011).
2.4 High level specifications
|Topic||Sp. Band||Mag.||error||CP error|
|Herbig AeBe disks||H,K||5%||5deg|
|T Tauri disks||H,K||5%||2deg|
|Demonstrated||H||7.5 (AT)||15 – 3%||0.5deg|
The specifications of PIONIER have been defined as a trade-off between the requirements for the key astrophysical programs and the technical contraints due to the project small size and fast timeline. A summary is presented in Table LABEL:tab:obs_log. We used the experience gained with IONIC-3 (Berger et al. 2003), VINCI (Le Bouquin et al. 2004, 2006), and AMBER to make the following strategic choices:
the combination of 4 telescopes (6 baselines) simultaneously,
the use of an integrated optics beam combiner,
operation in the H- and K-band (not simultaneously),
a small spectral dispersion across few channels,
a fast readout of the camera.
As built, PIONIER matches all these requirements except for the K band extension which is contemplated in the future. The instrument has been integrated and commissioned with the H-band integrated beam combiner that was already available at the IPAG laboratory (Benisty et al. 2009).
3 Instrument description
Figure 2 summarizes the key elements of the PIONIER instrument. Functionally, it consists of the following subsystems.
3.1 Injection and Optical Delay Unit (IOPDU)
The Injection and Optical Delay Unit (IOPDU) injects the free-space beams from the VLTI into optical fibers. It includes a tip-tilt correction, an OPD modulation and a polarization control. The IOPDU is made of four strictly identical arms, one per VLTI beam.
A dichroic mirror extracts the required band from the VLTI beams, while other wavelengths are used to feed the infrared VLTI guiding camera IRIS. The translation stage that supports these optics has two observing positions (the H-band and K-band dichroic) and a free position to let the beam pass unaffected when PIONIER is not in use.
A modulation of the optical path length is introduced by a mirror mount in a Physik Instrument (PI) piezo translation stage. It provides an OPD range of m.
The tip-tilt mirror, mounted on PI piezo devices, allows the beam angle to be corrected up to a frequency of Hz with a range of (laboratory angle, corresponding to times the size of the magnified PSF).
A Lithium-Niobate plate of 2 mm thickness is used to compensate for the polarization phase-shift arising in the single-mode fiber of the IOBC (see Sec. 3.2). The amount of phase-shift between the vertical and horizontal axes is accurately adjusted by tilting differentially the plates in the individual beams.
The off-axis parabola focuses the light into the single mode, polarization maintaining fiber.
PIONIER’s fast tip-tilt correction comes in addition to the one provided by the telescopes. It mainly compensates for the additional contributions coming from the tunnel turbulence. The IRIS guiding camera provides the beam angle of arrival measurements at a rate from Hz to Hz depending on the target brightness (Gitton et al. 2004).
3.2 The Integrated Optics Beam Combiner (IOBC)
The Integrated Optics Beam Combiner (IOBC) takes as input the signals from the IOPDU, and delivers 24 interferometric outputs.
The IOBC is fed by polarization maintaining single-mode H-band fibers, whose lengths have been equalized with an accuracy of m. Consequently, these fibers introduce a negligible differential chromatic dispersion. The polarization phase-shift between the vertical and horizontal axes (i.e. the neutral axes of the IO chip, aligned with those of the fibers) is of the order of 1 fringe. This polarization phase-shift can be compensated by adjusting the inclination of the Lithium-Niobate plate in the IOPDU (see Sec. 3.1).
The design of the IOBC is described in detail in Benisty et al. (2009). The four incoming beams are split in three and distributed in the circuit in a pairwise combination. A so-called “static-ABCD” combining cell is implemented for each baseline. It generates simultaneously four phase states (almost in quadrature). Consequently, 24 outputs have to be read. The low-chromaticity phase shift is obtained through the use of specific waveguides with carefully controlled refraction index. This combiner can be used both in fringe-scanning mode (VINCI-like) or ABCD-like mode. There is currently one combiner available in the H band and one that is being developed for the K band (Jocou et al. 2010).
3.3 Imaging Optics and Dispersion Unit (IMODU)
The Imaging Optics and Dispersion Unit (IMODU) images the 24 outputs of the IOBC onto the camera’s focal plane, with or without spectral dispersion in the perpendicular direction (1, 3 or 7 spectral channels across the H band). This unit is the one that was used with the IONIC-3 instrument at IOTA (Berger et al. 2003). Its image quality permits to focus about 80% of the flux into a single pixel of the PICNIC camera (pixel size m). A Wollaston prism can be inserted to acquire separately the fringes for the two linear polarizations (vertical and horizontal axes). This device is only used to internally adjust the compensation of the polarization phase-shift (see Sec. 3.2). Once this adjustment is done, the Wollaston prism is removed and the fringes are acquired in natural (unpolarized) light without loss of contrast. Consequently the instrument provides no polarization information in targets.
In practice, these prisms are inserted and removed manually from the optical beam. The positioning of the Wollaston prism is repeatable at pixel. The positioning of the dispersive prism is repeatable at pixel. Altogether, these measurements are compatible with a stability of for the effective wavelength and for the flux splitting ratio between the outputs.
3.4 The detector (DET/DFE)
The detector (DET) comprises the focal plane array, the cryostat, the filter wheels, and the internal electronics. More electronic functions, including the digitizing of the video signal, are provided by the Detector Frontend Electronics (DFE). The detector is the PICNIC camera previously used at PTI and IOTA (Traub et al. 2004; Pedretti et al. 2004). Its electronics has been partially upgraded to improve the frame rate. The array is clocked at MHz and performs analog-to-digital conversions every s. The main limitation for speed is the time for the video signal to stabilize once a pixel has been addressed (s). The resulting time for the elementary operations reset and read are presented in Table 2. Consequently, the detector achieves a frame-rate of kHz when reading 168 pixels (24 outputs times 7 spectral channels) in a non-destructive way.
We characterized the readout noise following the same procedure as described by Pedretti et al. (2004). We found a noise of per read while Pedretti et al. reported per read on this camera. This noise is reduced to by averaging 8 successive analog-to-digital conversions of the same pixel. It is still possible to reduce the noise by averaging several frames in a non-destructive way, although it has a strong impact on the minimum integration time possible.
|Detector elementary||Required time|
|Reset full quadrant||205|
|Read 24 pixels 1 line||190|
|Read 24 pixels 7 lines||1150|
3.5 The Calibration and Alignment Unit (CAU)
The Calibration and Alignment Unit (CAU) serves dual functions: (1) to inject, as a substitute for the astronomical beams, mutually coherent beams derived from a common Tungsten source, for laboratory verification and health check purposes; (2) to reflect, via a corner cube mirror, a reverse-propagating beam from the output of the IOBC towards the IRIS guiding camera, for the purpose of verifying the alignment of the PIONIER instrument with the IRIS reference positions.
3.6 Control System and electronics (CS)
Control System (CS) includes hardware control of the instrument units, detector readout, quicklook and interaction with VLTI. The hardware has been designed to follow the standard VLTI architecture including an instrument workstation (WS), an Instrument Control System (ICS), a Local Control Unit (LCU) and a Detector Control System (DCS) running on the detector WS. The instrument WS is a standard Linux PC platform. The detector WS is an industrial Linux PC platform equipped with a PCI-7300A board used to generate clock signals and to acquire detector data. Both stations are running under the VLT Software.
Instead of the traditional LCUs, the instrument is controlled by an Embedded PC and several EtherCAT modules like remote IO and stepper motion controllers from the company Beckhoff. Through a fruitful collaboration with ESO, PIONIER has been used as a pilot project for testing these commercial components and developing the associated VLT Software extensions. This success not only tries to simplify the coding of ICS but it also attempts to introduce technologies envisaged for E-ELT instruments (Kiekebusch et al. 2010).
PIONIER electronics is housed in a cooled cabinet close to the instrument optical table. It includes the detector front-end electronics, the DC detector power supplies together with the detector WS, tip-tilt and scanning piezo controllers, EtherCAT, calibration lamps and control electronics for electro-mechanical devices.
Dialog between PIONIER and VLTI (e.g. sending new coordinates) is achieved through the VLTI Interferometer Supervisor Software. The detector WS is also listening to the VLTI Reflective Memory Network (RMN) to get real-time centroid information from the IRIS guiding camera.
The raw output data follow the standard defined by ESO for its interferometric instruments (Ballester & Sabet 2002), including the general FITS header. The only difference is the addition of a new dimension (depth) in the column DATAj of the binary table IMAGING_DATA. The third dimension of DATAj represents the consecutive analog-to-digital conversions of each pixel. Consequently, a typical night provides several gigabytes of data.
The Observation Software running on the instrument workstation coordinates the execution of an exposure for a given observing mode. From the high level point of view, PIONIER is operated via the Broker of Observing Blocks that executes Observing Blocks (OBs) fetched from the standard p2pp ESO software. OBs can be conveniently generated by the aspro2333http://www.jmmc.fr/aspro_page.htm preparation software from the Jean-Marie Mariotti Center (JMMC).
4 Operating PIONIER
This section describes the PIONIER operations, including regular day-time optical alignment, flux injection procedure, scientific fringe recording and regular calibrations.
4.1 Optical alignment
As a visitor instrument, PIONIER is maintained and operated by the PIONIER team. It is only operated in visitor-mode and the observer is responsible for checking the health of the instrument at the beginning and during the observing run. The day-to-day stability of the beam angle inside PIONIER is , that is a quarter of the magnified PSF size. This is well within the range of the tip-tilt actuators (). The day-to-day position stability of the fringes acquired with the internal calibration is better than m in OPD. The same typical variations are observed with the VLTI calibration source MARCEL. Again, this is well within the range of the modulation mirror (m) and does not impact the operations.
However, the two following adjustments should be carefully checked at the beginning of each observing run:
the vertical/horizontal positioning of the IOBC, in order to optimize the coupling of the output spots into the detector pixels ;
the tilt angle of the Lithium-Niobate plates, as the polarisation phase-shifts of the fibers drift by typically over a month timescale.
These adjustments remain valid for the few following days of a typical run.
4.2 Flux injection
The stellar images are stabilized in tip-tilt by the IRIS guiding camera located in the VLTI laboratory (Gitton et al. 2004), only a few meters away from PIONIER. The image shift caused by the difference of atmospheric refraction index between the H-band (PIONIER fibers) and the K-band (IRIS guiding camera) is handled by the Interferometer Supervisor Software. The software automatically offsets the IRIS guiding point so that the stellar image in the H-band is kept fixed whatever the airmass and pointing direction.
The amount of flux injected in the fibers is optimized by a grid search with the internal tip-tilt of PIONIER. We use a square grid of 10x10 steps of each. The resulting image is fitted by a 2D gaussian. The procedure is applied sequentially on the four beams, during the evening twilight, using the beacon sources located on the VLTI telescopes. The resulting accuracy is better than of the PSF size. The precision is dramatically reduced by speckle noise when the procedure has to be performed on stellar light.
4.3 The interferometric signal
The fringes are temporally scanned across the coherence length by means of the long-range OPD modulation devices. Figure 3 shows examples of raw interferograms obtained in the high SNR and low SNR regimes. The four beams are modulated at non redundant velocities in order to generate the six fringe signals. The speed is defined by the ratio between the length of the smallest scan (set to m) and the duration of a scan. For faint stars, its optimal value is a compromise between sensitivity and the necessity to freeze the atmospheric effects for the fringes at the lowest frequency . For bright stars, is imposed by the detector frame rate and the need to correctly sample the fringe at the highest frequency .
The most used scan setups are summarized in Table 3. In mode “fowler”, the detector is reset once and is left integrating non destructively during the entire scan. This mode achieves the best frame rate and the lowest readout noise (when oversampling the fringes). In mode “double”, the detector performs a reset-read-read sequence for each step in scan. This mode is used on bright stars that would saturate the detector in “fowler” mode. Its main drawback is to slow down dramatically the frame rate. Consequently, it is not possible to use the mode “double” with the spectral dispersion over 7 channels as the resulting frame rate would be too slow compared with the atmospheric coherence time.
|Detector||# of spectral||steps in||steps per fringe|
4.4 Quick-look and fringe tracking
During the observation, each scan is processed by a quick-look algorithm that also implements a slow group-tracking technique:
The spectral channels are added together in a pseudo broad-band signal. The four outputs of each baseline are summed together taking into account their relative phase shift. The result is a single interferometric signal per baseline represented in red in Fig. 3.
At the end of the scan, the signal is processed by the algorithm presented by Pedretti et al. (2005). It provides a SNR and a group-delay position for each baseline.
Then these six group-delay measurements are converted into offsets for the four input beams. The inversion algorithm only takes into account the baselines with an SNR larger than 3. The redundancy is used to recover the offsets of a maximal number of beams. The results are an optimal estimate of the offsets to be applied to the beams and a detection flag per beam (0/1).
For all beams with a positive detection flag, the zero positions of the scanning piezo-electric devices is updated with the new offset value, before the next scan starts.
The typical repetition rate of the loop (Hz) is enough to keep the group-delay residuals within m under median atmospheric conditions. The redundancy of the 6 baselines to recover 3 optical path differences makes the tracking robust against temporary flux losses or strongly resolved baselines.
4.5 Observations and calibration
Each pointing typically produces 5 files of 100 fringe scans. Then a dark exposure is recorded in a file of 100 scans with the same detector setup but with all shutters closed. The flux splitting ratios of each input are finally calibrated by files of 50 scans with one beam illuminated at a time.
The sequence takes between min and min depending on the modulation speed. An additional 3 min are needed to preset the VLTI telescopes and delay lines and acquire the star. A single calibrated point (science + calibrator) takes between to min. It is possible to skip the optimisation of the telescope guiding system if the V magnitudes of the science and calibration stars are within mag. It significantly speeds up the target re-acquisition (min instead of min).
The spectral calibration is performed at the end of the night. It consists of long scans of m with high SNR recorded in the internal CAU. The effective wavelength is estimated by Fourier Transform Spectrometry. We checked the linearity of the OPD modulation by observing a m infrared laser. The day-to-day repeatability of the effective wavelength is , which is consistent with the optical stability presented in the previous section. Taking into account the systematic uncertainties, the final accuracy on the effective wavelength is . It matches the requirement of the key programs presented in Sec. 2. In particular, it meets the required 1% accuracy on the calibrated visibility associated to the most demanding program, for which the visibility models are not very sensitive to the wavelength calibration (Sec 2.2: faint debris disks around marginally resolved photospheres).
5 Data reduction and performances
In this section, we focus on the steps of the data reduction that are specific to the PIONIER instrument, or that have not already been published.
The data reduction software of PIONIER converts the raw FITS file produced by the instrument into calibrated visibilities and closure phase measurements written in the standard OIFITS format (Pauls et al. 2005). These output files are science-ready and can be directly handled by software such as LITpro (model fitting, Tallon-Bosc et al. 2008) or MIRA (image reconstruction, Thiébaut 2008). The PIONIER data reduction software is written as a yorick444http://yorick.sourceforge.net package called pndrs and publicly available555http://apps.jmmc.fr/~swmgr/pndrs.
Our strategy is directly inspired from the FLUOR, VINCI and IONIC-3 experiments. In these instruments, the scanning method associated with an estimate of the visibility in the Fourier space provided sensitive and accurate measurements, without the need for tricky internal calibrations. Basically, the intensity measured at the output (respectively B, C and D) of the baseline of the combiner is written as:
where are the baselines (12, 13, 14, 23, 24, 34), are the photometric flux from input beams, and are the instrumental visibilities and phases of the beam combiner, are the internal flux splitting ratio, and is the OPD modulation across the scan plus the atmospheric piston. The complex visibility of the observed target is .
The basic operations that are related to the PICNIC detector are performed first. The consecutive non-destructive reads are subtracted accordingly to the detector mode (“double” or “fowler”) to provide the intensity. The dark level is removed. In the “fowler” mode of the detector, a dark level is estimated independently for each step in scan.
The flux splitting ratios are estimated from the dedicated calibration files. They are computed independently for each wavelength. We use these flux splitting ratios to compute what we call the flat-fielded kappa matrix:
and the flat-fielded signal:
We use the instrumental visibilities coming from laboratory measurements of the beam combiner, actually close to one.
5.2 Disentangling the photometric and interferometric signals
The pairwise design allows the instantaneous photometric signal to be extracted by a global fit to the coherent data. The main idea is to add together the outputs and from a given baseline, after flat-fielding:
The opposite phase relation between and (respectively and ) due to energy conservation removes the interferometric part of the signal:
with . Equation 5 defines a linear system between the six signals (one per baseline) and the four photometric signals (one per beam). This system can be inverted and is well constrained thanks to the pairwise combination. We obtain an estimate of the for each step in scan. This method to recover the photometric signals was used in the PIONIER precursor instrument IONIC-3 (Monnier et al. 2004). Although it is not widely known, this method is a fundamental advantage of multi-telescope pairwise combination (see $ 4.1.2 in Blind et al. 2011).
We build a single interferometric signal per baseline combining the information of all outputs. We first inject the estimates into the flat-fielded kappa matrix in order to remove the continuum part of the interferometric signals. Then we co-add the ,, and outputs taking into account their relative phase shift (estimated from laboratory measurements):
5.3 Closure phase estimate
The interferometric signal is filtered in Fourier space in order to keep only the frequencies near the modulation frequency of each baseline. Then the baselines are combined together into bispectra, that are integrated over the length of the scan:
Four bispectra are formed , , and . Figure 4 shows examples of bispectrum statistics obtained for 100 scans in the high SNR and low SNR cases. The phase terms related to the fringe modulation and the atmospheric piston cancel out and all bispectrum estimates point toward the same phase: the closure phase . We take the argument of the bispectrum average to estimate the final closure phase. We use a bootstraping method to compute its standard deviation.
5.4 Visibility estimate
The critical step is the photometric normalization of the interferometric signal by the photometric signal . We implement two versions.
For bright stars, the photometric normalization is performed independently for each step in scan. This simultaneity is the key to obtain precisely calibrated interferograms on bright stars. As discussed by Coudé Du Foresto et al. (1997) and Kervella et al. (2004b), we use a Wiener filtering of . Still this is not enough to avoid numerical instabilities when get close to zero. We filter out the scans with such a behavior.
For fainter stars, the photometric signal is averaged over the length of a scan. When dominated by the detector noise, it dramatically enhances the precision on the photometric signal. In a sense this is an extreme case of the previous method, when the Wiener filtering is so strong that it is equivalent to a simple time average.
After the photometric normalization, we use the method detailled by Coudé Du Foresto et al. (1997), applied independently to each spectral channel and baseline, to obtain an estimate of the square visibility . This method is based on the Power Spectral Density (PSD) of the signal to disentangle the fringe power (peak) from the detector and photon noise biases (continuum). Typical PSD are presented in Fig. 6. The shape of the underlying bias power can be correctly fitted with the form where is the wavenumber. The first term corresponds to the photon noise and the second term corresponds to the detector noise. This underlying structure of the noise is due to the “fowler” readout scheme: the detector noise is high-pass filtered by the difference between two consecutive reads.
5.5 Transfer function calibration
The instrumental and atmospheric contributions to the visibilities and closure phases (the so-called transfer function) are monitored by interleaving the observation of science targets with calibration stars. We implemented several strategies to interpolate the transfer function at the observing time of the science targets. We found that the most robust is to simply average all the transfer function estimates of the night, with a weight inversely related to the time separation. The uncertainty on the transfer function at the time of the scientific observations takes into account both the error bars and the dispersion of the transfer function estimates: a sequence of precise (small error bars) but dispersed points leads to large uncertainties on the calibrated product. An example of a transfer function with its interpolation is displayed in Fig. 6.
In this section, the term median atmospheric conditions refers to seeing, ms, and wind speedm/s.
The limiting magnitude of PIONIER is defined by its capability to detect the fringes during the observations (SNR per scan). The typical measured sensitivity is for a dispersion over 7 spectral channels and using a scan of 512 steps with the detector in mode “fowler”. When the atmospheric conditions are better than the median, it is possible to reach a limiting magnitude of using a scan of 1024 steps. In broad band, we have been able to detect and track the fringes on several unresolved stars with under good conditions. These numbers are for the Auxiliary Telescopes.
The statistical uncertainties on the closure phases (error bars in Fig. 6, right) are compatible with the data dispersion. The final accuracy of the calibrated closure phases generally ranges from deg to deg. These uncertainties not only depend on the target brightness but also on the turbulence strength. Sequences with stable closure phases, down to deg, have been recorded on bright stars under atmospheric conditions better than median. At such a level, new possible biases should be studied such as the dependence of the instrumental closure phase on the spectral type. This has not been done yet.
The measured accuracy on calibrated visibilities ranges from to depending on the atmospheric conditions. On bright targets, the dispersion of the uncalibrated data is sometimes larger than the statistical uncertainties computed by pndrs (Fig. 6, left). It means that we are sometimes facing some non-stationary biases. Our prime suspects are the mechanical OPD perturbations that could be related to the wind strength, perhaps by shaking the telescope structure. Indeed, when the wind speed is higher than m/s, the PSDs recorded by PIONIER are clearly affected by vibrations, as seen in Fig. 6 left. These perturbations are typically in the range Hz, with a typical life-time of a few seconds. This should be further investigated before drawing definite conclusions.
The VINCI and IONIC-3 instruments achieved a visibility accuracy of , using a similar instrumental concept to PIONIER (see for instance Kervella et al. 2004a). We plan to dedicate a few observing nights to define the best observing setup to reach the same level of accuracy.
6 Image reconstruction of Sco and HIP11231
In this section we illustrate the imaging capability of PIONIER on the stars Sco and HIP11231. Sco is a well known bright Be star () accompanied by a faint companion on a highly eccentric orbit. HIP11231 is a G1V spectroscopic binary with apparent magnitude . The PIONIER observations are presented in Appendix A and B together with their plane coverage. Sco has been observed once with the intermediate configuration (D0-H0-G1-I1). HIP11231 has been observed at 7 epochs, and with the three AT quadruplets (D0-H0-G1-I1, E0-G0-H0-I1, A0-K0-G1-I1). Data have been processed with the pndrs package.
6.1 Image reconstruction and model fitting of Sco
We perform model-independent image reconstruction with the MIRA software from Thiébaut (2008). The starting point is a delta function at . The pixel scale is mas and the field-of-view is pixels. In practice, the MIRA algorithm minimizes a joint criterion which is the sum of (1) a likelihood term which measures the compatibility with the data, and (2) a regularization term which imposes priors on the image. The relative weight between these two terms is controlled by a multiplicative factor called “hyperparameter”. We use the “total variation” regularization associated with positivity constraint as recommended by Renard et al. (2011). We set the hyperparameter with a small value of 100, so that the weight of the regularization term is kept small with respect to the fit to the data. It brings some superresolution, at the cost of an increased level for the noise in the image. We use the information of all the spectral channels to improve the plane coverage (“gray image” hypothesis).
The reconstructed image is presented in Fig. 7. The companion is easily detected. The main component appears extended as the disk surrounding it is marginally resolved by the VLTI baselines. This result illustrates the snapshot imaging capability of PIONIER thanks to the simultaneous beam combination of 4 telescopes.
In parallel, we perform a fit with a binary model with one resolved component. Results are presented in Table LABEL:tab:binaireDeltaSco. The resulting visibility and closure phase curves are overlaid to the data in Appendix A. We checked that the position of the companion is consistent with the prediction by Tycner et al. (2011). This image reconstruction of Sco demonstrates that the pndrs package computes the sign of the closure phases and of the plane in a consistent manner. In other words, it means that the produced OIFITS files are compliant with the sky orientation defined by Pauls et al. (2005).
6.2 Image reconstruction, model fitting and orbital solution of HIP11231
We use the same MIRA parameters as for Sco except that the pixel scale is set to mas. The reconstruction is performed independently for each epoch. We use the information of all the spectral channels to improve the plane coverage (“gray image”). Reconstructed images are presented in Fig. 9. The spectroscopic binary is easily resolved by PIONIER+VLTI with the three AT configurations. The noise in the images is kept below 2.5% of the image maximum. The flux ratio between the primary and the secondary, measured by aperture photometry in each reconstructed image, ranges from to .
In parallel, we perform a global fit to the full data set with the following parameters: (1) the flux ratio, forced to be identical for all epochs and (2) the relative positions, that we let free between each epoch. We combined together all the spectral channels to improve the plane coverage. Such a global fit is conveniently performed with the LITpro software (Tallon-Bosc et al. 2008). The best-fit flux ratio is . At each epoch, the relative separation could be retrieved unambiguously. The fitted positions and the flux ratio are consistent with the results from model-independent image reconstruction. The resulting visibility and closure phase curves are overlaid to the data in Appendix B.
Figure 9 shows the relative separations derived from the fit overlaid with the best orbital solution. We adopt the spectroscopic period derived by Barker et al. (1967) from radial velocity, but we allowed all the other Campbell elements to vary, namely: the time passage through periastron , the eccentricity , the semi-major axis , the position angle of the ascending node , the argument of periastron and the orbital inclination . Parameters of the best orbital solution are presented in Table LABEL:tab:orbit together with the solution determined from Hipparcos astrometry by Jancart et al. (2005). The difference in semi-major axis is due to the fact that Jancart et al. (2005) reconstructed the orbit of the barycenter of light from unresolved astrometry, while we reconstruct the orbit of the secondary from resolved astrometry. The difference in time passage through periastron can be explained by the uncertainties on the period and by the fact that our observations have been taken several orbital periods later than the Hipparcos ones. The disagreement in is more striking as we cannot find a simple explanation.
7 Conclusions and perspectives
The 4-telescope PIONIER instrument has been successfully developed and commissioned at the VLTI interferometer. The use of integrated optics technologies combined with the expertise gained from projects such as FLUOR, AMBER, VINCI and IONIC-3 made it possible to design, build and commission the instrument in about 18 months. PIONIER is routinely used for science operations in the H-band and permits high angular resolution imaging studies at an unprecedented level of sensitivity and precision. The overall performances are in agreement with the requirements from the key programs, although some work is still needed to achieve the best possible accuracy on the calibrated visibilities and closure phases.
PIONIER is expected to stay a few years at VLTI as a visitor instrument, until the arrival of the second generation instruments MATISSE and GRAVITY. Within this lifetime, several improvements are contemplated:
The imaging and dispersion unit will be upgraded with a motorized translation stage for the dispersion prisms. It will allow to change the operational mode (broad band or dispersed) during the night. This should be installed in PIONIER around December 2011.
Changing the operational wavelength from the H-band to the K-band would benefit all the key science programs. A component is already available at IPAG. The design of the associated imaging optics is under study. Ideally, the changes will be installed in PIONIER early 2012.
The arrival of a new generation of infrared detectors should provide detector noise as low as at a frame-rate larger than kHz in destructive mode (Gach et al. 2009; Finger et al. 2010). Installed in PIONIER, such a detector would increase the sensitivity by mag and would solve the readout speed and saturation issues.
PIONIER has the status of a visitor-instrument at Paranal Observatory. The possibility to offer PIONIER to the community is currently investigated.
Acknowledgements.PIONIER is funded by the Université Joseph Fourier (UJF, Grenoble) through its Poles TUNES and SMING and the vice-president of research, the Institut de Planétologie et d’Astrophysique de Grenoble, the “Agence Nationale pour la Recherche” with the program ANR EXOZODI, and the Institut National des Science de l’Univers (INSU) with the programs “Programme National de Physique Stellaire” and “Programme National de Planétologie”. The integrated optics beam combiner is the result of a collaboration between IPAG and CEA-LETI based on CNES R&T funding. The authors want to warmly thank all the people involved in the VLTI project. This work is based on observations made with the ESO telescopes. It made use of the Smithsonian/NASA Astrophysics Data System (ADS) and of the Centre de Donnees astronomiques de Strasbourg (CDS). All calculations and graphics were performed with the freeware Yorick.
- Absil et al. (2006) Absil, O., di Folco, E., Mérand, A., et al. 2006, A&A, 452, 237
- Absil et al. (2008) Absil, O., di Folco, E., Mérand, A., et al. 2008, A&A, 487, 1041
- Absil et al. (2009) Absil, O., Mennesson, B., Le Bouquin, J.-B., et al. 2009, submitted in ApJ
- Akeson et al. (2005) Akeson, R. L., Boden, A. F., Monnier, J. D., et al. 2005, ApJ, 635, 1173
- Akeson et al. (2009) Akeson, R. L., Ciardi, D. R., Millan-Gabet, R., et al. 2009, ApJ, 691, 1896
- Ballester & Sabet (2002) Ballester, P. & Sabet, C. 2002, VLTI Data Interface Control Document, Tech. rep., ESO
- Barker et al. (1967) Barker, E. S., Evans, D. S., & Laing, J. D. 1967, Royal Greenwich Observatory Bulletin, 130, 355
- Benisty et al. (2009) Benisty, M., Berger, J.-P., Jocou, L., et al. 2009, Astronomy and Astrophysics, 498, 601
- Benisty et al. (2010) Benisty, M., Natta, A., Isella, A., et al. 2010, A&A, 511, A74+
- Benisty et al. (2011) Benisty, M., Renard, S., Natta, A., et al. 2011, A&A, 531, A84+
- Berger et al. (2003) Berger, J., Haguenauer, P., Kern, P. Y., et al. 2003, in Interferometry for Optical Astronomy II. Edited by Wesley A. Traub. Proc. of SPIE, Vol. 4838, pp. 1099-1106., 1099–1106
- Berger et al. (2001) Berger, J. P., Haguenauer, P., Kern, P., et al. 2001, A&A, 376, L31
- Blind et al. (2011) Blind, N., Absil, O., Le Bouquin, J.-B., Berger, J.-P., & Chelli, A. 2011, A&A, 530, A121+
- Bordé et al. (2002) Bordé, P., Coudé du Foresto, V., Chagnon, G., & Perrin, G. 2002, A&A, 393, 183
- Coudé Du Foresto et al. (1997) Coudé Du Foresto, V., Ridgway, S., & Mariotti, J.-M. 1997, A&AS, 121, 379
- di Folco et al. (2007) di Folco, E., Absil, O., Augereau, J.-C., et al. 2007, A&A, 475, 243
- Dullemond et al. (2001) Dullemond, C. P., Dominik, C., & Natta, A. 2001, ApJ, 560, 957
- Eisner et al. (2009) Eisner, J. A., Graham, J. R., Akeson, R. L., & Najita, J. 2009, ApJ, 692, 309
- Finger et al. (2010) Finger, G., Baker, I., Dorn, R., et al. 2010, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 7742, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Gach et al. (2009) Gach, J.-L., Balard, P., Daigle, O., et al. 2009, in EAS Publications Series, Vol. 37, EAS Publications Series, ed. P. Kern, 255–270
- Gillessen et al. (2010) Gillessen, S., Eisenhauer, F., Perrin, G., et al. 2010, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 7734, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Gitton et al. (2004) Gitton, P. B., Leveque, S. A., Avila, G., & Phan Duc, T. 2004, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 5491, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. W. A. Traub, 944–+
- Haguenauer et al. (2010) Haguenauer, P., Alonso, J., Bourget, P., et al. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7734, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Isella & Natta (2005) Isella, A. & Natta, A. 2005, A&A, 438, 899
- Jancart et al. (2005) Jancart, S., Jorissen, A., Babusiaux, C., & Pourbaix, D. 2005, A&A, 442, 365
- Jocou et al. (2010) Jocou, L., Perraut, K., Nolot, A., et al. 2010, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 7734, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Kervella et al. (2004a) Kervella, P., Coude Du Foresto, V., Segransan, D., & di Folco, E. 2004a, in SPIE Conf. Series, Vol. 5491, New Frontiers in Stellar Interferometry, Proc. of SPIE, Vol.me 5491. Edited by Wesley A. Traub. Bellingham, WA: The International Society for Optical Engineering, 2004., p.741, ed. W. A. Traub, 741–+
- Kervella et al. (2004b) Kervella, P., Ségransan, D., & Coudé du Foresto, V. 2004b, A&A, 425, 1161
- Kiekebusch et al. (2010) Kiekebusch, M. J., Chiozzi, G., Knudstrup, J., Popovic, D., & Zins, G. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7740, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Lafrasse et al. (2010) Lafrasse, S., Mella, G., Bonneau, D., et al. 2010, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 7734, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Le Bouquin et al. (2006) Le Bouquin, J.-B., Labeye, P., Malbet, F., et al. 2006, A&A, 450, 1259
- Le Bouquin et al. (2004) Le Bouquin, J.-B., Rousselet-Perraut, K., Kern, P., et al. 2004, A&A, 424, 719
- Lopez et al. (2008) Lopez, B., Antonelli, P., Wolf, S., et al. 2008, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 7013, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Malbet et al. (1999) Malbet, F., Kern, P., Schanen-Duport, I., et al. 1999, A&AS, 138, 135
- Mérand et al. (2005) Mérand, A., Bordé, P., & Coudé Du Foresto, V. 2005, A&A, 433, 1155
- Millan-Gabet et al. (2007) Millan-Gabet, R., Malbet, F., Akeson, R., et al. 2007, Protostars and Planets V, 539
- Monnier et al. (2004) Monnier, J. D., Traub, W. A., Schloerb, F. P., et al. 2004, ApJ, 602, L57
- Pauls et al. (2005) Pauls, T. A., Young, J. S., Cotton, W. D., & Monnier, J. D. 2005, PASP, 117, 1255
- Pedretti et al. (2004) Pedretti, E., Millan-Gabet, R., Monnier, J. D., et al. 2004, PASP, 116, 377
- Pedretti et al. (2005) Pedretti, E., Traub, W. A., Monnier, J. D., et al. 2005, Appl. Opt., 44, 5173
- Petrov et al. (2007) Petrov, R. G., Malbet, F., Weigelt, G., et al. 2007, A&A, 464, 1
- Renard et al. (2011) Renard, S., Thiébaut, E., & Malbet, F. 2011, ArXiv e-prints
- Tallon-Bosc et al. (2008) Tallon-Bosc, I., Tallon, M., Thiébaut, E., et al. 2008, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 7013, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Tannirkulam et al. (2008) Tannirkulam, A., Monnier, J. D., Harries, T. J., et al. 2008, ApJ, 689, 513
- Thiébaut (2008) Thiébaut, E. 2008, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 7013, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
- Traub et al. (2004) Traub, W. A., Berger, J.-P., Brewer, M. K., et al. 2004, in New Frontiers in Stellar Interferometry, Proc. of SPIE, Vol. 5491. Edited by Wesley A. Traub. Bellingham, WA: The International Society for Optical Engineering, 2004., p.482, 482–+
- Tycner et al. (2011) Tycner, C., Ames, A., Zavala, R. T., et al. 2011, ApJ, 729, L5+
- Zhao et al. (2011) Zhao, M., Monnier, J. D., Che, X., et al. 2011, PASP, 123, 964
- Zhao et al. (2010) Zhao, M., Monnier, J. D., Che, X., et al. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7734, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
Appendix A Observations of Sco
|Date||Average MJD||Baseline||# of||# closure||Time spend|
Appendix B Observations of HIP11231
|Date||Average MJD||Baseline||# of||# closure||Time spend|