Astrometric search for a planet around VB 10^{⋆}
^{1}
Main Astronomical Observatory, National Academy of Sciences of the
Ukraine, Zabolotnogo
27, 03680
Kyiv, Ukraine
email: laz@mao.kiev.ua
^{2}
Observatoire de Genève, Université de Genève,
51 Chemin des Maillettes,
1290
Sauverny,
Switzerland
^{3}
INTACSIC Centro de Astrobiología, 28850 Torrejón de Ardoz, Madrid, Spain
^{4} Instituto de Astrofísica de Canarias, C/ Vía Láctea s/n,
38200 La Laguna (Tenerife), Spain
^{5}
Centro de Astrofísica, Universidade do Porto,
Rua das Estrelas,
4150762
Porto,
Portugal
^{6} Departamento de Física e Astronomia, Faculdade de Ciências,
Universidade do Porto, Portugal
Received:
13
July
2010
Accepted:
5
November
2010
We observed VB 10 in August and September 2009 using the FORS2 camera of the VLT with the aim of measuring its astrometric motion and of probing for the presence of the announced planet VB 10b. We used the published STEPS astrometric positions of VB 10 over a timespan of 9 years, which allowed us to compare the expected motion of VB 10 due to parallax and proper motion with the observed motion and to compute precise deviations. The singleepoch precisions of our observations are about 0.1 mas, and the data showed no significant residual trend, while the presence of the planet should have induced an apparent proper motion greater than 10 mas yr^{1}. Subtraction of the predicted orbital motion from the observed data produces a strong trend in position residuals of VB 10. We estimated the probability that this trend is caused by random noise. After taking all the uncertainties into account and using MonteCarlo resampling of the data, we are able to reject the existence of VB 10b with the announced mass of 6.4 M_{J} with a false alarm probability of only 5 × 10^{4}. A 3.2 M_{J} planet is also rejected with a false alarm probability of 0.023.
Key words: astrometry / techniques: high angular resolution / planetary systems / stars: individual: VB 10
© ESO, 2011
1. Introduction
The detection of a planetarymass companion to the nearby M8 ultracool dwarf VB 10 (GJ 752B) was announced by Pravdo & Shaklan (2009). Using groundbased astrometric observations over nine years with the STEPS camera of the 5 m Palomar telescope, the authors derived a full astrometric solution for the system that suggests a planetary companion with a 271 day orbital period and 6.4 M_{J} mass. Even though the VB 10b astrometric signature is large with a peaktopeak value of 8 mas, the detection is debatable because of imperfect orbital sampling and the comparably low astrometric accuracy of the individual measurements.
Zapatero Osorio et al. (2009) presented highresolution nearinfrared observations of VB 10 obtained between 2001 and 2008 with the NIRSPEC instrument on the Keck II telescope. They emphasize the necessity for a better sampling of the orbital phase to constrain the orbital parameters and the individual masses of the system precisely. Bean et al. (2010) used the nearinfrared CRIRES spectrograph and did not detect the planet’s signature in radialvelocity (RV). AngladaEscudé et al. (2010) performed a joint analysis of precise RV observations over 175 days with MIKE/Magellan and ESPaDOnS/CFHT spectrographs, RV measurements by Zapatero Osorio et al. (2009), and astrometric data by Pravdo & Shaklan (2009). They show that the observed astrometric motion does not come from an unseen companion.
We present results of the astrometric search for the planet using the FORS2 camera of the VLT, the only instrument that can achieve astrometric precision of individual measurements of 0.05–0.1 mas (Lazorenko et al. 2009). Despite the short observation period of two months, we obtained good enough data to verify the hypothesis that VB 10b exists.
2. Observations
Observations were made with the FORS2/VLT camera (Appenzeller et al. 1998) in imaging mode during five nights from 2 August to 25 September 2010. The images cover a 4.2 × 4.2′ field of view with a pixel (px) scale of ~0.1″/px. For each night, we obtained 21 to 72 frames of 7 s exposures with the R_{special} filter. Seeing varied from 0.55″ to 0.9″. To reduce differential chromatic refraction (DCR) of the atmosphere, we observed near meridian.
Because the observed motion of VB 10 in the sky is a long arc (over 2 px), we had to know the precise pixel scale to convert pixel to arcsec units without loss of accuracy. From the CMC14 and NOMAD catalogue positions of stars in our reference frame, we derived acceptably precise scales of 0.12538 ± 0.00007″/px and 0.12567 ± 0.00005″/px in RA and Decl, respectively.
3. Data reduction and analysis
Our null hypothesis was that the planet VB 10b exists. To test it, we computed the residuals of our measurements compared to the position of VB 10 expected from the combination of proper, parallactic, and orbital motions. By determining the probability that the measured residuals are compatible with the expected mean value of zero, we could confirm or reject the null hypothesis. We also performed the corresponding computations by assuming that the planet does not exist.
In all cases, we considered that the astrometric signal Ψ in VB 10 position induced by the planet is defined by the orbital elements of VB 10b given in the discovery paper by Pravdo & Shaklan (2009). The predicted orbital motion over the measurement timespan of 17 days is 0.48 mas in RA and 0.92 mas in Dec, and for 54 days it has 1.95 mas in RA and 3.28 mas in Dec. Compared to the astrometric precision of FORS2 of about 0.1 mas per epoch, these displacements are large and should be detected in this study.
3.1. Photocentre determination
Raw images were flatfielded and biassubtracted to exclude pixeltopixel variations in the CCD sensitivity. To increase the number of reference stars, we measured all star images with R = 14–21. At this faint end, the star field is very crowded. This forced us to improve our procedure of photocentre computation, which was initially developed for isolated stellar images (Lazorenko 2006). We scanned images obtained under various seeing conditions and made a detailed census of star positions at subpixel precision, of fluxes, and of the pointspread function variation across the CCD. This information was used to accurately model and subsequently subtract the contamination of background counts caused by distant stars. Wings of star profiles at distances up to 50 px were approximated by exponential function with free parameters smoothly varying over the CCD. The photocentres and ŷ were computed by fitting star profiles in a 11 × 11 px window with a model with 12 free parameters and an auxiliary oscillating function in the central 5 × 5 px window (Lazorenko 2006). Computations based on the LevenbergMarquardt numerical algorithm of the least squares fit were found to give good enough results for stars of approximately equal brightnesses at separations larger than about 10 px.
3.2. Astrometric model
Reductions were based on the method previously applied to FORS1/2 observations (Lazorenko et al. 2007, 2009). The method was shown to efficiently mitigate atmospheric image motion, geometric field distortion, DCR, and other effects, thus ensuring stable precision over time scales of a few days to a few years. For every star in the field, i.e. VB 10 and the reference stars, the measured photocentre positions and ŷ_{m} in frame m at time t were represented by the model (1)The lefthand side contains the free model parameters and the model function . The expression for ydata is similar but contains cosγ instead of sinγ. Here, x_{0} is a zero point, k the mode (an arbitrary even integer, usually from 4 to 16), a polynomial in x and y of order k/2 − 1, which models the sum of atmospheric image motion and geometric distortion for each frame m. The parameter μ_{x} is the proper motion, π the parallax, and p_{x} the parallax factor in x. The displacement of the star image due to DCR is modelled by a term with leading parameter ρ, which depends on the star’s colour, and contains the zenith distance z and the angle γ between a direction to zenith and yaxis. The next term describes an image displacement opposite to that of DCR and introduced by the longitudinal atmospheric dispersion compensator (LADC) of the VLT (Avila et al. 1997). This displacement also depends on the star colour via the parameter d ≈ −ρ. Both ρ and d are free model parameters, and their computation does not require external colour data. The LADC is automatically adjusted to the average zenith distance z_{L} over a given series of frames by setting a distance b between its two prisms to b ~ tanz_{L}. Finally, Ψ(t) represents the induced orbital motionof VB 10 if the planet exists. For reference stars, Ψ(t) = 0.
Equation (1) defines a system of equations in the combined {x, y, t}domain for which the unique solution is derived under the condition that the model parameters of reference stars are orthogonal to each basic function of Φ_{k}. It is solved iteratively for all reference stars available in a circular region of radius R_{k}, which increases with k and is centred on the target VB 10. The optimal field size R_{k} is the size at which the noise from a reference field is equal to the noise from atmospheric image motion. The number of reference stars depending on k varies from a few dozen to 500. Thanks to the large number of reference stars, the function Φ accurately reproduces the coordinate grid distortion introduced by image motion. Solutions for Φ_{km} at each mode k are then used to form equations for the target only in the time domain. The solution of this new set of equations yields the target’s model parameters and position deviations from the model. The final output is the average obtained from the solutions at all modes k.
3.3. Comparison of STEPS and FORS2 reference frames
Because of different zeropoints, STEPS proper motion cannot be directly applied to FORS2 astrometry. The STEPS reference frame is given by 15 bright stars. For the FORS2 reference frame, we used both these (except STEPS star Nr. 10 which is a binary) and a number of much fainter stars. The system of proper motions is defined under the condition ∑ μ_{i}w_{i} = 0 where summation is taken over all reference stars i with weightsw_{i} approximately proportional to their brightness (Lazorenko et al. 2009). In particular, this condition is valid for the proper motions of the STEPS stars, which due to their high brightness (they contribute to over 40% of the total light flux), are the basis of the FORS2 reference frame. One may therefore expect that ∑ μ_{i}w_{i}/ ∑ w_{i} (the weighted mean proper motion in the system of FORS2 proper motions) taken over these stars, is nearly zero. The computed value of this sum (the difference of proper motion zeropoints) was found to be small: − 0.1 ± 1.0 and + 0.3 ± 1.0 mas yr^{1} in RA and Decl. The low precision meant that we did not apply it but included it in the error budget assuming that the uncertainty in the proper motion zero point is ~1 mas yr^{1}. Thus, both proper motion systems are consistent, but only on time scales of ~0.1 yr.
3.4. Short and full dataset
The data reduction of the first epoch (2 Aug) was problematic because of incorrect telescope pointing, owing to which all stars in a 300 px wide area just below VB 10 (Fig. 1a) were imaged to another CCD chip and could not serve as a reference on this particular night. We dealt with this problem in two ways. First, we put aside all stars within this area from frames of all nights, using only stars outside of this area as reference objects. The sizes and shapes of the reference fields in this case are shown in Fig. 1a by five circular segments, each corresponding to k increasing from six (small radius) to 14 (large radius). Alternatively, we used circular reference areas, which for the 2 Aug. night were the same R_{k} size but with the lower half vignetted (Fig. 1b). We also examined a short dataset without the 2 Aug. epoch with circular reference area (Fig. 1c). The asymmetry of the first two cases degraded the precision obtained compared to the symmetric configurations of the short dataset.
We did obtain a solution for all three cases represented in Fig. 1. For the final result, we present two solutions. One is obtained from the short dataset (case “c”, MJD between 55082 and 55099 days), and the second obtained as the average of the results in cases “a” and “b”. The latter solution is referred to as having been obtained from the full dataset (MJD between 55045 and 55099 days).
Fig. 1 Circular reference fields increasing in size with increasing modal number k in the field of view of FORS2 in the case of: a) the full dataset when an area unavailable on 2 Aug. frames was cut away at all nights; b) the same but the problem area cut only at 2 Aug. frames; c) the short dataset not using 2 Aug. frames. 

Open with DEXTER 
3.5. Recovering LADC positions
VB 10 is very red and differs considerably in colour from the reference stars. Therefore, the DCR displacement of this star is very large and about 30 mas in Decl and 10 mas in RA. This is taken into account by the free model parameters ρ and d, but their adjustment requires knowledge of the LADC separation b.
Fig. 2 Distribution of epoch position residuals ⟨ y ⟩ of field stars (dots of different type) as a function of d for four September nights and linear dependences of ⟨ y ⟩ on d (lines of different type) for each night. Upper panel: with initial LADC separations b. Lower panel: the same with recovered b values. Size of dots is proportional to the star brightness. 

Open with DEXTER 
Information on b was not part of the obtained fits headers and was not accessible to us. Therefore we had to solve the inverse problem of recovering LADC separations from the observations. Because b is fixed for a series of frames on a single night, we had to determine a small number of values for b (one value per night). Our approach is based on the observation that incorrect values of b bias the average night position residuals ⟨ x ⟩ and ⟨ y ⟩ of field stars. According to Eq. (1), this bias linearly depends on the star colour parameter d during one given night. This is illustrated by Fig. 2, which shows the distribution and linear dependencies of ⟨ y ⟩ on d for field stars in the case of the short dataset (the effect is greatest in y).
Because all observations were made at small hour angles within ± 0.7 h, we initially assumed that LADC was always set to the separation corresponding to the meridian. By applying small corrections to these initial values of b, we iteratively reached a solution without dependence of ⟨ y ⟩ on colours (lower panel of Fig. 2). For these computations, we used all stars within the entire field of view of FORS2, processing them as target objects relative to their own subsets of reference stars. In spite of being based on field stars only, this procedure produced a much smaller dispersion of ⟨ y ⟩ values for VB 10 and for reference stars with extreme colours. Corrections to b were small and within ±4% of the initial values and were found with a precision corresponding to ±0.025 mas error in ⟨ y ⟩ for VB 10. For the full dataset, precision degraded to ±0.068 mas becausewe could not use stars in the problematic 300 px area below VB 10 (Fig. 1a).
3.6. Treatment of colourdependent terms
An essential drawback of this procedure is that the recovery of b introduces small colourdependent terms μ′(d) and π′(d), which are similar to proper motion and parallax. This is because the above iterations do not garantee convergence to the actual values of b ~ tanz_{L}. For example, the restored value of tanz_{L} may differ from its actual value by a term that progressively changes in time as , where is an arbitrary constant and ⟨ sinγ ⟩ is the average value of sinγ during the given night. As a result, the term dtanz_{L}sinγ in Eq. (1) generates the term , which is linearly dependent on d and which therefore can be treated as an extra image motion . The term μ′(d) compensates for the linear change of tanz_{L} in time. Similarly, we may assume that the restored values of tanz_{L} contain terms proportional to p_{x} and p_{y}. In this case, the solution of Eq. (1) for parallax should contain the compensating colourdependent term π′(d).
When recovering b values, we cannot control the amplitude of and of the equivalent parameter related to parallax. However, they can be detected as a proper motion and parallax dependence on colours (i.e. on d), which we model as a linear trend in proper motions of field stars and statistically correct for it.
The strong correlation between t, p_{x}, and p_{y} does not allow us to determine μ′(d) and π′(d) separately. However, this is not required because, for short times p_{x} ~ t, we can approximate the sum of the proper motion and parallax displacement by , where μ′ is the new effective quantity that substitutes for both μ′ and π′. Thus by processing the short dataset, we find and use only a single colour term μ′(d). The treatment is similar for the y components. Subtracting μ′t from the measured positions thus simultaneously eliminates the colour term π′ (see next section). The value of this correction for VB 10 in Decl is μ′ = 29 ± 4.5 mas yr^{1} for the full and μ′ = 27 ± 2.9 mas yr^{1} for the short dataset, respectively. In RA, the corrections are an order of magnitude smaller.
For the full dataset, the above approximation of by is still valid for reference stars, most of which have moderate d and therefore very small colourinduced parallaxes π′. But it is not precise enough for VB 10. Because of its extreme colour term d, the expected π′ value is also large and cannot be accurately eliminated by applying the correction μ′(d) as in the case of the short dataset. Therefore, we had to treat the term π′ of VB 10 as a free model parameter.
Fig. 3 Orbital motion Ψ(t) of VB 10 in RA (circles) and Decl (triangles) for the full and short FORS2 datasets and the linear approximations (solid and dashed lines) of these motions. 

Open with DEXTER 
3.7. Subtraction of parallax and proper motions
Since the FORS2 observations cover a small fraction (<20%) of the orbital period, Ψ(t) can be approximated as a linear function of t: (2)where Q_{α} = ∂Ψ_{α}(t)/∂t and Q_{δ} = ∂Ψ_{δ}(t)/∂t. With the orbital elements given by Pravdo & Shaklan (2009), we obtained Q_{α}(PS) = 13.4 and Q_{δ}(PS) = 22.4 mas yr^{1} (Fig. 3). Astrometric acceleration terms (deviations from the linearity) are smaller than 0.1 mas for the full dataset, so are negligible. For the short dataset, they are even smaller than 0.02 mas.
In addition, Ψ(t) is approximately a linear function of the parallax factors p_{x}(t) and p_{y}(t), which themselves have an approximately linear time dependence. This causes a strong correlation between parallax, proper motion, and orbital motion, making Eq. (1) degenerated. Therefore, we subtracted parallax and proper motion from the measured positions of VB 10. Precise values of these parameters were found based on the published STEPS astrometric measurements which cover a 9 year period. The best fit of STEPS data yielded μ_{α}cos(δ) = −586.8 ± 0.2 mas yr^{1}, μ_{δ} = −1361.0 ± 0.2 mas yr^{1}, and π = 168.0 ± 1.2 mas. These values are very close to the estimates given by Pravdo & Shaklan (2009) and AngladaEscudé et al. (2010).
3.8. Position residuals of models with and without a planet
Subtracting the contributions of parallax and proper motion from the measured positions of VB 10 resulted in simplifying Eq. (1), which for VB 10 took the form (3)with only four free parameters x_{0}, y_{0}, ρ, and d in the case of the short dataset. For the full dataset, it incorporates the extra parameter π′ (the term in curly braces, see Sect. 3.6). The righthand side of Eq. (3) contains the measured coordinates and ŷ_{m}, the inplate solutions and derived from the reference stars, and all other corrections.
There is a correlation between the orbital signal and the night’s average angles ⟨ sinγ ⟩ , which are not zero and tend to increase in time. Because of this correlation, the model Eq. (1) filters out any component of the signal (e.g. Ψ), which is linearly dependent on t and enters the righthand side of Eq. (1). Therefore, the output position residuals only contain a part of the initial signal amplitude. However, this does not hamper our statistical analysis, which we performed after the subtraction of Ψ(t) from the observed positions, thus assuming zero input signal and consequently zero output signal. In this way, the impact of correlations is minimized. If the planet does not exist but the subtraction of Ψ(t) was applied, we should detect the inverse signal −Ψ(t) reduced in amplitude because of the correlation between ⟨ sinγ ⟩ and Ψ(t). This is the case corresponding to the last row of Table 2 (Sect. 4), where the measured signal (expressed by the parameter Q) has about half its expected value −Q(PS).
Fig. 4 Frametoframe deviations in Decl of VB 10 from its model motion at zero (Ψ = 0, filled circles) and nonzero planetary signal (Ψ ≠ 0, open circles), the epoch average deviations (large diamonds), and the singleframe precisions σ_{1} (asterisks) for the cases a), b), and c) of reference field configurations shown in Fig.1. 

Open with DEXTER 
Figure 4 shows the results in terms of model deviations in Decl, where the expected signal is strongest. The singleframe precision σ_{1} includes errors in the photocentre measurements, the reference frame noise, and the atmospheric noise. It varies from 0.4 to 0.7 mas depending on seeing. The effect of the vignetted reference field of the 2 Aug. epoch (configurations “a” and “b”) is seen as a degradation of σ_{1} to over 1 mas. At other epochs, σ_{1} is larger compared to configuration “c” because of a larger R_{k}.
Clear conclusions can be drawn from the short dataset when the model (Eq. (3)) is simplestand does not require incorporation of parallax (Sect. 3.6). We considered the cases with the predicted orbital motion subtracted (Ψ ≠ 0) and with Ψ = 0. The epoch average deviations ⟨ x ⟩ and ⟨ y ⟩ are very small and randomly scattered when assuming Ψ = 0, but display a negative trend in time if Ψ ≠ 0. Small position deviations do not correspond to a “zero” measurement. Instead, they demonstrate very precise position measurements, which track the proper motion and parallax displacement at the daily rate of 2.7 mas and 5.0 mas in RA and Dec, respectively. The motion of VB 10 over the CCD surface (Fig. 5) for the measurement timespans of 17 days is 46 mas in RA and 86 mas in Decl, and is dominated by parallax and proper motion. DCR effects induce a smallscale scatter in the measured positions of one night with an amplitude of about 2 mas. Their structure for a typical night is shown in 5fold magnification in Fig. 5.
Similar computations were performed for the full dataset, where we had to account for the parallax correction π′(d) (see Sect. 3.6), so used the five parameters x_{0}, y_{0}, ρ, d, and π′ for fitting of the model Eq. (3). This introduced a strong correlation between p_{x}, p_{y}, and Ψ, which reduced the amplitude of the detected signal.
Table 1 and Fig. 6 summarize the results for the epoch residuals ⟨ x ⟩ and ⟨ y ⟩ . The astrometric precision is described by a nominal precision σ_{N} based on errors in photocentre determination, the reference frame noise and atmospheric noise. The term σ_{sum} also includes error components that dominate at long time spans and come from the uncertainties in b, μ′, pixel scale, and proper motion of VB 10. For the full dataset, it also includes the uncertainty in the parallax of VB 10. The mathematical expectation of the rootmeansquare of ⟨ x ⟩ and ⟨ y ⟩ , derived from the least squares fit (Eq. (3)) is given by σ_{fit}.
4. Discussion
Because the signal Ψ was subtracted in Eq. (3), the epoch deviations ⟨ x ⟩ and ⟨ y ⟩ should have an expectation value of zero and a random scatter irrespective of the model (with or without planet). A small dispersion of ⟨ x ⟩ and ⟨ y ⟩ values indicates that the tested hypothesis is correct. The data shown in Fig. 6 and Table 1 is more consistent with Ψ = 0 (the hypothesis of no planetary companion), while incorporating a nonzero planetary signal Ψ ≠ 0 increases the scatter of data points, which does not support the existence of the planet VB 10b. We estimated the probability P that the planet exists in two ways.
The first estimate P(χ^{2}) was based on the values of , where the sum is taken over epochs. The number of degrees of freedom (d.o.f.) is equal to the number of epochs minus the number of model parameters related to the time domain. There are two (x_{0}, y_{0}) and three (x_{0}, y_{0}, π′) effective parameters for the short and the full datasets, respectively. The probabilities P(χ^{2}) corresponding to the χ^{2} values are given in Table 2. For both the short and the full datasets, they are equally small at the 2% level if the model assumes Ψ ≠ 0. Because the probability of this event is low, the large scatter of ⟨ x ⟩ and ⟨ y ⟩ is unlikely to be caused by random noise in the observations.
Fig. 5 Measured inframe positions , ŷ − Φ^{ { y } } (squares) of VB 10, the model motion (solid line) over the CCD between 8 and 25 Sept. (short dataset) defined primarily by proper motion and parallax, and the magnified DCR signature for a single night (upper right corner). 

Open with DEXTER 
Epoch deviations ⟨ x ⟩ and ⟨ y ⟩ in RA and Decl of VB 10 for the models with Ψ = 0 and Ψ ≠ 0 (Eq. (3)).
Fig. 6 Astrometric residuals in the scenario with planet (open blue symbols) and without planet (filled black symbols) for the short dataset (top panel) and the full dataset (bottom panel) as a function of time. Circles and triangles mark residuals in RA and Decl, respectively. For clarity, only one errorbar corresponding to σ_{sum} is shown at each epoch, but the uncertainties for the four respective measurements are identical. 

Open with DEXTER 
Probabilities P(χ^{2}), P(Q) and linear trends Q (mas yr^{1}).
Secondly, we noticed that the observed residuals ⟨ x ⟩ and ⟨ y ⟩ show a linear trend instead of a random distribution. This can be caused by random errors in the observations, but also may indicate a wrong value of the signal Ψ(t) subtracted in Eq. (3). In the latter case, we would expect that the residuals ⟨ x ⟩ and ⟨ y ⟩ show a linear dependence −Qt opposite to Ψ(t) (see Eq. (2)).
A strong negative trend is seen for the short dataset residuals computed with Ψ ≠ 0 (Fig. 6), whereas they should be near zero to support the planetary hypothesis. In contrast, a much weaker trend is observed when assuming Ψ = 0 (no planet). For instance, the measured values of Q given in Table 2 for the short dataset and Ψ ≠ 0 are –7.82 and –11.32 mas yr^{1} in RA and Decl, respectively. If the null hypothesis Ψ ≠ 0 is correct, their expected values are 0 and 0 mas yr^{1}, and –13.4 and –22.4 mas yr^{1} if the hypothesis is wrong (Ψ = 0). In the last case, the measured Q values are always lower than the values Q(PS) induced by the orbital motion, because their magnitude is damped by correlations between Ψ(t) and model parameters (cf. Sect. 3.8).
Because the χ^{2} criterion is not the most efficient one for characterizing linear signals, we developed another, more powerful approach to obtaining an alternative estimate of P. We considered Q as quantities describing the dispersion of ⟨ x ⟩ and ⟨ y ⟩ better than χ^{2} and performed Monte Carlo simulations to see if the observed features can be explained by random noise in the observations. We simulated each data frame by a random Gaussian noise with rootmeansquare of σ_{1} and added components modelling the errors in μ′, b, and in the proper motion and parallax of VB 10. In addition, we included a 1 mas yr^{1} error to account for the uncertainties in the zero point of proper motion (cf. Sect. 3.3).
After fitting the data with the model Eq. (3), we searched the residuals for linear trends in time and estimated their coefficients Q^{∗}. The falsealarm probabilities P(Q) that a random noise produces higher Q^{∗} values than the observed ones are given in Table 2. In the case of the short dataset and Ψ ≠ 0, we find that the observed linear trends can be explained by Gaussian
noise with a probability of P(Q) = 5 × 10^{4}, which does not support the existence of the planet. For the full dataset, we find P(Q) = 0.02, which is not sufficient to draw definite conclusions. In spite of the longer timespan, hence seemingly better conditions for characterizing a planetary signal, the full dataset does not provide a better constraint because of the uncertainty in μ′, whose contribution increases with time. Besides, the planetary signal is substantially damped because of its strong correlation with the colour correction π′ used in Eq. (3) as a free model parameter.
5. Conclusions
We conclude that the presence of the announced planet around VB 10 is not supported by astrometry. Even assuming half the planetary mass (i.e. 3.2 M_{J}), simulations give a low falsealarm probability P(Q) = 0.023, which still does not prove the existence of VB 10b. Our result obtained from astrometry alone agrees with the conclusion of Bean et al. (2010) and AngladaEscudé et al. (2010) based on RV data.
This study is the first application of the FORS2 camera in the search for exoplanets by means of optical astrometry. Because of the high astrometric precision of FORS2, the availability of external STEPSbased proper motion and parallax of VB 10, and the strong expected orbital signal, it was possible to verify the planetary companion hypothesis within the extremely short observation period of 17 days, which is unusual for astrometric works of this type. Successful use of STEPS data for reduction of FORS2 observations is only justified because we verified that the STEPS and FORS2 proper motion reference frames are consistent within 1 mas yr^{1} uncertainty.
We have demonstrated a mean nominal precision of 0.09 mas per epoch of FORS2/VLT observation for data of reasonable quality, despite the problems caused by the uncertainty in the LADC position. This precision is sufficient for astrometric detection of planets around ultracool dwarfs.
Acknowledgments
We thank Dr. G. AngladaEscudé whose comments have helped to improve the paper. P.F. and N.C.S. would like to acknowledge the support by the European Research Council/European Community under the FP7 through a Starting Grant, as well as the support from the Fundação para a Ciência e a Tecnologia (FCT), Portugal, in the form of a grant with reference PTDC/CTEAST/098528/2008. NCS would further like to acknowledge the support from Fundação para a Ciência e a Tecnologia (FCT), Portugal, through a Ciência 2007 contract funded by FCT/MCTES (Portugal) and POPH/FSE (EC).
References
 AngladaEscudé, G., Shkolnik, E. L., Weinberger, A. J., et al. 2010, ApJ, 711, L24 [NASA ADS] [CrossRef] (In the text)
 Appenzeller, I., Fricke, K., Fürtig, W., et al. 1998, The Messenger, 94, 1 [NASA ADS] (In the text)
 Avila, G., Rupprecht, G., & Beckers, J. M. 1997, in SPIE, 2871 (In the text)
 Bean, J. L., Seifahrt, A., Hartman, H., et al. 2010, ApJ, 711, L19 [NASA ADS] [CrossRef] (In the text)
 Lazorenko, P. F. 2006, A&A, 449, 1271 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Lazorenko, P. F., Mayor, M., Dominik, M., et al. 2007, A&A, 471, 1057 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Lazorenko, P. F., Mayor, M., Dominik, M., et al. 2009, A&A, 505, 903 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Pravdo, S. H., & Shaklan, S. B. 2009, ApJ, 700, 623 [NASA ADS] [CrossRef] (In the text)
 Zapatero Osorio, M. R., Martín, E. L., del Burgo, C., et al. 2009, A&A, 505, L5 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
All Tables
Epoch deviations ⟨ x ⟩ and ⟨ y ⟩ in RA and Decl of VB 10 for the models with Ψ = 0 and Ψ ≠ 0 (Eq. (3)).
All Figures
Fig. 1 Circular reference fields increasing in size with increasing modal number k in the field of view of FORS2 in the case of: a) the full dataset when an area unavailable on 2 Aug. frames was cut away at all nights; b) the same but the problem area cut only at 2 Aug. frames; c) the short dataset not using 2 Aug. frames. 

Open with DEXTER  
In the text 
Fig. 2 Distribution of epoch position residuals ⟨ y ⟩ of field stars (dots of different type) as a function of d for four September nights and linear dependences of ⟨ y ⟩ on d (lines of different type) for each night. Upper panel: with initial LADC separations b. Lower panel: the same with recovered b values. Size of dots is proportional to the star brightness. 

Open with DEXTER  
In the text 
Fig. 3 Orbital motion Ψ(t) of VB 10 in RA (circles) and Decl (triangles) for the full and short FORS2 datasets and the linear approximations (solid and dashed lines) of these motions. 

Open with DEXTER  
In the text 
Fig. 4 Frametoframe deviations in Decl of VB 10 from its model motion at zero (Ψ = 0, filled circles) and nonzero planetary signal (Ψ ≠ 0, open circles), the epoch average deviations (large diamonds), and the singleframe precisions σ_{1} (asterisks) for the cases a), b), and c) of reference field configurations shown in Fig.1. 

Open with DEXTER  
In the text 
Fig. 5 Measured inframe positions , ŷ − Φ^{ { y } } (squares) of VB 10, the model motion (solid line) over the CCD between 8 and 25 Sept. (short dataset) defined primarily by proper motion and parallax, and the magnified DCR signature for a single night (upper right corner). 

Open with DEXTER  
In the text 
Fig. 6 Astrometric residuals in the scenario with planet (open blue symbols) and without planet (filled black symbols) for the short dataset (top panel) and the full dataset (bottom panel) as a function of time. Circles and triangles mark residuals in RA and Decl, respectively. For clarity, only one errorbar corresponding to σ_{sum} is shown at each epoch, but the uncertainties for the four respective measurements are identical. 

Open with DEXTER  
In the text 