A&A 480, 193-198 (2008)
DOI: 10.1051/0004-6361:20078229
K. Geißler1,2 - G. Chauvin1,3 - M. F. Sterzik1
1 - European Southern Observatory, Alonso de Cordova 3107, Vitacura, Santiago, Chile
2 -
Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany
3 -
Laboratoire d'Astrophysique, Observatoire de Grenoble, UJF, CNRS, BP 53,
38041 Grenoble Cedex 9, France
Received 5 July 2007 / Accepted 8 December 2007
Abstract
Context. Brown dwarfs exhibit complex atmospheric signatures, and their properties are highly sensitive to effective temperature, surface gravity, and metallicity. Several physical properties of brown dwarfs in binary systems can be well inferred from the primary, and therefore allow one to better constrain their atmospheres.
Aims. We want to constrain atmospheric models of brown dwarfs in binary systems using narrow-band mid-infrared photometry.
Methods. High spatial resolution and sensitivity is required to resolve the components. Therefore we have obtained deep mid-infrared images of four close binary systems with brown dwarf companions using VISIR at the Very Large Telescope (VLT) in three narrow-band filters at 8.6, 10.5 and 11.25
m.
Results. Three brown dwarf companions (GJ 229 B, HD 130948 BC and HR 7329 B) were detected at 8.6
m. HD 130948 BC was also observed at 10.5
m. Finally, we establish upper flux limits for the other narrow band filters with null detections.
Conclusions. Our results, are in general, compatible with previous observations and model expectations for these objects. For HD 130948 BC, we conclude that there is photometric variability on a significance level of 2.8
based on repeated observations. The bandpass around 10.5
m appears specifically well suited for variability studies, and we speculate that either inhomogeneities in the atmospheric NH3 distribution, or silicate absorption, might cause its time-variability.
Key words: stars: low-mass, brown dwarfs - stars: fundamental parameters
Brown dwarfs (BDs) bridge the gap in mass between low-mass stars and giant planets. Hundreds of them have been discovered in the past decade, mainly in wide-field optical (SDSS, Stoughton et al. 2002) and near-infrared (e.g., 2MASS, DENIS; Cutri et al. 2003; Epchtein et al. 1997) surveys. Two main classes of BDs emerged based on their optical and infrared spectral properties, the L-dwarfs and T-dwarfs. They typically cover effective temperature ranges of
K and
K, respectively. The modeling of atmospheres cooler than
K is a challenge, because it must include an appropriate treatment of a plethora of molecular opacity's and dust processes (formation, condensation, size distribution and mixing). The most recent atmosphere models include additional properties such as age (gravity) and metallicity, and seem to reproduce the spectral signatures and the infrared colours of L and T dwarfs reasonably well. Only the L-T transition, occurring around a relatively narrow temperature range
K, remains problematic (for a discussion of state-of-the-art models see Burrows et al. 2006).
As an additional complication, the determination of distances, ages and metallicity for field BDs remains difficult and may be the cause of various uncertainties in their physical parameters. In that sense, BD companions of nearby stars are extremely interesting, as the system characteristics are well known thanks to the bright primary star. For these systems the influence of age and metallicity on the spectral signatures of substellar objects can be explored with high accuracy.
Table 1: Target properties.
Initially, most BD studies were conducted in the optical and near-infrared. However, with the recent development of mid-infrared observing techniques and instruments from the ground or from space, this spectral region became essential because numerous molecular lines, including dominant absorption bands of
and
,
are present. The Spitzer Space Telescope has provided a wealth of high-quality photometry and low- and medium resolution spectroscopy of L and T dwarfs in the 5-20
m regime (see, e.g. Leggett et al. 2007; Cushing et al. 2006; Mainzer et al. 2007). They enable one to probe vertical mixing, clouds and non-equilibrium chemistry around the L-T boundary in unprecedented detail.
The power of sensitive, ground-based, high spatial resolution, mid-infrared imaging has recently been demonstrated for the close BD binary companion
Indi Ba and Bb (see Sterzik et al. 2005). The relative photometry between both components allows one to derive effective temperatures independently from the determinations in near-infrared. As the distance is well known for the primary, the absolute photometry constrains radii and bolometric corrections, in contrast to Spitzer observations (Roellig et al. 2004) that suffered from insufficient angular resolution to resolve both components (Spitzer diffraction limit at 10
m is
3
).
In order to continue our effort to constrain atmospheric and evolutionary models of BDs, we have therefore conceived a mini-survey of close BD companions using the mid-infrared imager VISIR at the VLT. Section 2 reviews the target's properties and describes our observations. Section 3 describes our data reduction and analysis process. Section 4 presents the results of our survey, i.e. the photometry and the astrometry of the detected companions and the sensitivity limits obtained. Finally, Sect. 5 compares our results with the predictions of cool atmosphere models, to discuss their applicability and limitations.
For our target selection, we have only considered the confirmed members of binary (or multiple) systems with known distances, as their primaries are well characterized in terms of metallicity and age. Only BD companions with expected mid-IR fluxes stronger than 1 mJy and separations larger then 0.5
were selected, in order to fully adapt and exploit the sensitivity and spatial resolution of the mid-infrared imager VISIR at the VLT. We finished with a short list of four systems: GJ 229, HD 130948, HR 7329 and HR 7672. Their main characteristics are summarized in Table 1.
Table 2: Observing Log.
To ensure that the primary and the companion are within the FoV during chopping and nodding, the system was aligned horizontally on the detector. At the same time, this simplifies the reduction process, since the shifting and adding of the frames can be done using the brighter primary star. A summary of the observing log is given in Table 2, including the mean airmass during the observing run and the total integration time.
The reduction was done using self-written IDL routines for bad-pixel replacement and for the shifting and adding of the frames. Bad-pixels were replaced by the mean value of the surrounding pixels within a box of 9
9 pixel, before subtracting (A-B) nodding positions. In the following, the relative shifts between the frames of one data-set were calculated via cross-correlation of the bright primary, before the frames were averaged. Since the VISIR detector was affected by randomly triggered stripes during part of the observations, a destriping technique developed by Pantin et al. (2007) was applied to the final co-added images.
Standard aperture photometry was used to determine the relative photometry of all detected BDs. Using IDL ATV routines, a curve-of-growth method was applied to the brown dwarf companions to obtain the apertures where the signal-to-noise ratio is maximised. In the following, those apertures were used for the primary as well as for the standard stars. Thus the count-rate to flux conversion factor was determined and relative photometry obtained. The variation of the count-rate to flux conversion factors with aperture radius was screened for at least 3 consecutive aperture radii between 4 and 7 pixels (corresponding to radii of 0.3
to 0.525
). At 10
m the VISIR diffraction limit is
0.3
,
the chosen aperture radii are of the order of twice the diffraction limit. To calibrate the flux values different standard stars
, observed before and after the targets, were used. The error bar estimation of the flux calibration is derived from the flux variations of the source measured in different apertures and from different standard stars. In cases where two independent measurements were taken, an average of the measured fluxes is quoted in Table 4.
As already mentioned, a destriping technique was applied to the final images to clean it of random stripes and thereby improve the image quality. To estimate the impact of the destriping on the photometry of the BDs, we performed the aperture photometry before and after the destriping process. In those cases in which the source is not located close to a stripe, no influence is noticeable. In contrast in the cases where the BD is close to a stripe, a decrease in the measured count-rate, and consequently in flux, of the BD in the destriped images is perceivable. Nevertheless, this effect is expected, since the stripes in the ``raw'' images fall within the aperture radii and lead to an overestimation of the count-rate and therefore of the flux.
![]() |
Figure 1:
VISIR detection images of GJ 299B at |
| Open with DEXTER | |
![]() |
Figure 2:
VISIR images of HD 130948 in the SIV filter taken on the 5th of August ( left) and on the 3rd of August ( middle). The binary brown dwarf companion was detected in the data from the 5th of August and is marked by a box in the leftmost image. The flux was measured to be 5.7 |
| Open with DEXTER | |
Three of the four brown dwarfs were detected in PAH1, namely GJ 229 B, HR 7329 B and HD 130948 BC (see Fig. 1), while only HD 130948 BC could be detected in SIV. Note that HD 130948 BC, a binary brown dwarf, was not resolved in our observations. HR 7672 B could not be detected in any of the filters. While the resolution of VISIR is sufficient to separate the brown dwarf and the primary (
0.79
Liu et al. 2002; assuming negligible orbital motion), the data quality in PAH1 and SIV is low. The PSF of the primary is elongated, affecting the area in which the brown dwarf is expected, and thus adding noise.
Table 3: Separations and position angles of the detected brown dwarfs.
In Table 3 the measured separations and position angles of the detected brown dwarfs are given. To obtain the separation as well as the position angle of the brown dwarfs relative to their primaries, the pixelscale and N-orientation provided in the image header were used. Golimowski et al. (1998), used the HST's Wide Field Planetary Camera 2 (WFPC2) to observe GJ 229 B at three epochs, which were spread over approximately one year. Orbital motion of GJ 229 B was detected and a relative change of separation of (0.088Table 4: VISIR photometry of the primaries and brown dwarfs. In cases of non-detections, upper limits are provided. The fluxes are quoted in mJy.
As a final step we compare our obtained photometry to the models developed by Allard et al. (2001) and Burrows et al. (2006). Using their theoretical spectra, provided online
, absolute model fluxes were calculated by integrating the theoretical spectrum over the VISIR filter bandpasses. The object radii R, which determine the absolute spectral flux calibration, were obtained from evolutionary calculations by Burrows et al. (1997). In Table 5 the calculated model fluxes are listed. Furthermore, the age and effective temperature combinations for which the object radius was determined are given.
From Allard et al. (2001) we employed the AMES-cond and AMES-dusty models, representing the two extreme cases, in which either all dust has disappeared from the atmosphere (AMES-cond) or dust settling throughout the atmosphere is negligible (AMES-dusty). Following Allard et al. (2001) the AMES-dusty models should successfully describe dwarfs with effective temperatures greater than 1800 K, while the AMES-cond models are better suited to describe the atmospheres of dwarfs with
K.
Table 5:
Predicted mid-IR fluxes from different theoretical models. Values consistent within 3
or with the given upper limit were marked in boldface.
![]() |
Figure 3: Comparison of the limiting background obtained from the two data sets of HD 130948 taken in the SIV filter. The over-plotted point corresponds to the detection on the 5th of August, with a flux of 5.7 mJy. |
| Open with DEXTER | |
GJ 229 B:
Saumon et al. (2000), have used high-resolution infrared spectra to determine the metallicity, effective temperature and gravity of the T dwarf (see also Table 1). While using an age of 0.2 Gyr they derived an effective temperature of 950 K
80 K and a gravity of log g = 5
0.5. Later on, Leggett et al. (2002) compared the observed low- and high- resolution spectra of GJ 229 A and GJ 229 B to theoretical spectra (AMES-models). Their best fit yields an
= 1000
100 and a gravity of log g
3.5 for GJ 229 B, as well as a metallicity of
and an age of
30 Myr (range 16-45 Myr) for the system. While the metallicity is determined within the spectra fitting procedure, the age is derived by a comparison with evolutionary models, and mainly constrained by the observed luminosity and the derived effective temperature of the A component. Using VISIR mid-IR photometry, absolute model predictions of both Allard et al. (2001) and Burrows et al. (2006) can then be tested for different combinations of
,
[M/H] and gravity (when available, e.g. see Fig. 4). From Table 5, we see that the PAH1, SIV and PAH2 photometry are consistent with model predictions for a
= 900 K, log g = 5.0 and [M/H] = 0 companion. At solar metallicity, an effective temperature of
= 1000 K can be excluded at more than 2 sigma. At subsolar
[M/H] = -0.5 metallicity, low gravity (log g < 4.5) values remain excluded for
= 1000 K. Therefore, excluding young ages predictions of Leggett et al. (2002), the VISIR photometry clearly favours the initial physical parameters proposed by Saumon et al. (2000) for solar and subsolar metallicities.
![]() |
Figure 4: GJ 229 B: comparison of theoretical spectra from Allard et al. (2001) and Burrows et al. (2006) with the VISIR photometry. |
| Open with DEXTER | |
HD 130948 BC: As already mentioned, HD 130948 BC is a binary brown dwarf consisting of two L 4 dwarfs. To compare our observations to the model predictions, we assumed that both brown dwarfs contribute equally to the measured fluxes, as the simplest assumption. An unequal distribution of the measured flux on one of the two brown dwarfs would only increase the afterwards described effect. The suggested effective temperature of
1900 K places the binary in the regime of the AMES-dusty models. While the predicted PAH1 flux is in good agreement with our measurements, most models fail to reproduce the SIV flux detected on August 5, 2006, but would be consistent with the non-detection on August 3, 2007. Using the models provided by Burrows et al. (2006) we tested the influence of different metallicities ([M/H]=0.0 and
[M/H]=-0.5) on the predicted fluxes. The change is marginal, only about 0.05 mJy.
While analysing the SIV data obtained for HD 130948, we found that the binary companion was detectable in only one of the two data sets. A possible explanation of this result is an intrinsic variability around 10.5
m of one or both of the L dwarfs in the binary. Variability at 10.5
m could either be caused by ammonia (
)
or silicates. Since the
absorption features, which were first identified in the mid-infrared spectra taken with Spitzer/IRAS by Roellig et al. (2004), and Cushing et al. (2006), appear at roughly the L/T transition this is unlikely to be the cause of the observed variability.
A more favourable explanation may be an inhomogeneous distribution of silicate clouds, which characterise the atmospheres of L dwarfs with effective temperatures of roughly 1400-2000 K (Burrows et al. 2001). Future VISIR observations of HD 130948 BC at 10.5
m over different timescales should secure this photometric variability.
Using VISIR at the VLT, we performed a mini-survey of brown dwarfs in binary systems. The four selected brown dwarfs were imaged in three narrow band filters at 8.6, 10.5 and 11.25
m. At 8.6
m three of the brown dwarfs were detected and photometry was obtained. None of the brown dwarfs was detected at 11.25
m, and only HD 130948 BC was detected at 10.5
m. The observations of HD 130948 BC at 10.5
m indicate a possible variation of one or both brown dwarfs of the binary.
To constrain the atmospheric properties of the brown dwarfs we compared the mid-infrared photometry to theoretical model spectra by Allard et al. (2001) and Burrows et al. (2006). The measured mid-infrared fluxes and upper limit of GJ 229 B are consistent with the characteristic parameter obtained, by Saumon et al. (2000) (
950 K, log g
5, respectively), while values of the effective temperature and gravity as suggested by Leggett et al. (2002) (
K, log g
3.5, respectively) result in too high model fluxes. As for HD 130948 BC, the model fluxes for
1900 K, log g
5 fit the measurement at 8.6
m and the upper limits obtained at 10.5
m and 11.25
m. Nevertheless, the models are not in agreement with the flux measured for the detection of HD 130948 BC at 10.5
m during one observing epoch.