Abstract

We present new optical and near-infrared spectra of WD0137−349; a close white dwarf–brown dwarf non-interacting binary system with a period of ≈114 min. We have confirmed the presence of H α emission and discovered He, Na, Mg, Si, K, Ca, Ti and Fe emission lines originating from the brown-dwarf atmosphere. This is the first brown-dwarf atmosphere to have been observed to exhibit metal emission lines as a direct result of intense irradiation. The equivalent widths of many of these lines show a significant difference between the day-side and night-side of the brown dwarf. This is likely an indication that efficient heat redistribution may not be happening on this object, in agreement with models of hot Jupiter atmospheres. The H α line strength variation shows a strong phase dependency as does the width. We have simulated the Ca ii emission lines using a model that includes the brown-dwarf Roche geometry and limb darkening, and we estimate the mass ratio of the system to be 0.135 ± 0.004. We also apply a gas-phase equilibrium code using a prescribed drift-phoenix model to examine how the chemical composition of the brown-dwarf upper atmosphere would change given an outward temperature increase, and discuss the possibility that this would induce a chromosphere above the brown-dwarf atmosphere.

1 INTRODUCTION

Brown dwarfs occupy the mass range between low-mass stars and high-mass planets; they form like stars through molecular cloud fragmentation but never reach the hydrogen burning mass limit of ≈0.072 M (Burrows & Liebert 1993; Chabrier & Baraffe 1997). This lack of hydrogen fusion throughout their evolution results in different atmospheric properties to that of stars. In some respects there is a strong similarity between brown-dwarf and planetary atmospheres such as Jupiter (e.g. Helling & Casewell 2014). The dividing line on this continuum of objects is the deuterium burning mass limit, which occurs above 13 MJup. This rule was adopted by the Working Group on Extrasolar Planets of the International Astronomical Union in 2002 (Boss et al. 2007; Spiegel, Burrows & Milsom 2011). In order to draw clear parallels between planets and brown dwarfs, they need to be observed in similar environments i.e. close in to their host star.

Brown-dwarf companions that orbit their main-sequence star host within 3 au are rare compared to planetary or stellar companions to main-sequence stars (Grether & Lineweaver 2007). This is known as the brown-dwarf desert and does not appear to extend out to larger orbital radii (Metchev 2006). This desert is thought to exist due to the difficulties in forming extreme mass ratio binaries (q ∼ 0.02–0.1), although observational biases may play a part too (see fig. 2 in Burgasser et al. 2007). This is not the case for exoplanets that are thought to form through core accretion in a proto-planetary disc (Armitage 2010).

There have been a number of brown-dwarf companions detected around white-dwarf stars, but such systems are even rarer than brown dwarfs around main-sequence stars with only 0.5 per cent of white dwarfs having a brown-dwarf companion (Steele et al. 2011). These systems are thought to form through post-common envelope evolution. This is the process by which the white-dwarf progenitor evolves off the main sequence, expands and engulfs the brown-dwarf companion. The brown dwarf loses orbital angular momentum to the envelope thus causing the brown dwarf to spiral in towards the core and the envelope to be ejected (Politano 2004). This results in a detached system with a white dwarf and a close, tidally locked brown-dwarf companion known as a post-common envelope binary (PCEB).

Recent infrared (IR) studies have provided candidate systems (e.g. Debes et al. 2011; Girven et al. 2011; Steele et al. 2011) that have been identified by looking for IR excesses in SDSS spectra. Only a handful of the candidate systems have been confirmed as PCEBs: GD1400 (WD+L6, P = 9.98 h; Farihi & Christopher 2004; Dobbie et al. 2005; Burleigh et al. 2011), WD0137−349 (WD+L6-L8, P = 116 min; Burleigh et al. 2006; Maxted et al. 2006; Casewell et al. 2015), WD0837+185 (WD+T8, P = 4.2 h; Casewell et al. 2012), NLTT5306 (WD+L4-L7, P = 101.88 min; Steele et al. 2013), SDSS J141126.20+200911.1 (WD+T5, P = 121.73 min; Beuermann et al. 2013; Littlefair et al. 2014), SDSS J155720.77+091624.6 (WD+L3-L5, P = 2.27 h; Farihi, Parsons & Gänsicke 2017), SDSS J120515.80−024222.6 (WD+L0, P = 71.2 min) and SDSS J123127.14+004132.9 (WD+BD, P = 72.5 min; Parsons et al. 2017). These systems are thought to be progenitors of cataclysmic variables with a brown-dwarf donor (e.g. Littlefair, Dhillon & Martín 2003; Hernández Santisteban et al. 2016).

There are advantages to studying irradiated brown dwarfs over hot Jupiters. First, the atmospheres of brown dwarfs have been very well characterized (e.g. Cushing, Rayner & Vacca 2005) and the models associated with brown-dwarf atmospheres are thus more advanced than exoplanet models (e.g. Helling & Fomins 2013; Helling & Casewell 2014; Marley & Robinson 2015). Secondly, brown dwarfs are typically brighter than exoplanets and can be directly observed at wavelengths longwards of 1.2 μm. Thirdly, observing brown dwarf or planetary atmospheres around main-sequence stars is difficult due to the brightness of their host. A solution to this is to observe these systems in their more evolved form. A brown dwarf emits primarily in the near-infrared (NIR), whereas the white-dwarf contribution dominates at short wavelengths. Due to this, there is a high contrast between the brown dwarf and white dwarf, making it easier to separate their spectral components. Thus, the observation driven study of short period irradiated brown-dwarf atmospheres can give us insights into the chemistry of the atmospheres of hot Jupiters.

WD0137−349 has been the subject of several studies over the last decade e.g. (Burleigh et al. 2006; Maxted et al. 2006; Casewell et al. 2015). The white dwarf was first discovered using low resolution spectroscopy and added to the McCook & Sion catalogue (McCook & Sion 1999). Follow up high resolution spectroscopy by Maxted et al. (2006) revealed the presence of a brown-dwarf companion and concluded that it was not affected by the common envelope phase. 2MASS photometry and observations using the Gemini telescope and the Gemini Near-Infrared Spectrograph revealed an IR excess. This allowed Burleigh et al. (2006) to better characterize the brown dwarf and estimate its spectral type to be L6–L8, although these observations were taken of the unirradiated hemisphere of the brown dwarf. The brown dwarf intercepts 1 per cent of the primarily ultraviolet (UV) white-dwarf radiation and is likely to be tidally locked to its 16 500 K host. This results in one side being continually irradiated giving the brown dwarf a temperature difference of 500 K between its ‘night-’ and ‘day-’ side (Casewell et al. 2015).

Our observations and data analysis are summarized in the next section. In Section 3, we present our measured radial velocities and equivalent widths of all the emission lines detected, we also present the refined ephemeris of this system. In Section 4, we put WD0137−349B into the context of similar binaries. In Section 5, we discuss how we used a drift-phoenix model to do a first test of how the chemical composition of the brown dwarf upper atmosphere would change given an outward temperature increase. We discuss the possibility that this would induce a chromosphere above the brown-dwarf atmosphere.

2 OBSERVATIONS AND DATA REDUCTION

We obtained 78 spectra using the XSHOOTER instrument (Vernet et al. 2011) mounted at the ESO VLT-UT3 (’Melipal’) telescope in Paranal, Chile. We observed in the wavelength range 300.0–2480.0 nm using the three independent arms: VIS, UVB, NIR. The observations took place on the nights of 2014 August 28 between 04:45 and 09:45, and 2014 August 29 between 04:35 and 07:12 universal time. The 78 exposures, each 300 s long, were taken covering 2.4 orbits on the first night and 1.3 orbits on the second night. This gives us 21 spectra per full orbit. The weather conditions on the first night were good with an airmass ranging between 1.02 and 1.33 and seeing of 0.98 arcsec. The second night, the weather was not as good due to a high wind of 13 m s−1 at the start of the night and slowing to 8 m s−1. The airmass ranged between 1.02 and 1.58 and the seeing between 0.63 and 1.1 arcsec. The data and standards were taken in nodding mode and were reduced using gasgano (v2.4.3; Izzo et al. 2004) following ESO's xshooter pipeline (v2.5.2).

We see H α absorption and emission features that move in anti-phase over the orbital period of the system. The data were normalized and phase-binned using Tom Marsh's molly software. The H α absorption feature has a narrow core and broad wings. To fit this line, we followed the method from Casewell et al. (2015) and used three Gaussian profiles: two for the wings and one for the core. The best-fitting widths were 0.149, 2.360 and 8.30 nm, and the depth of the narrowest profile was fixed. A fourth Gaussian was used to model the H α emission, height and width were free to roam. This allowed us to measure the radial velocities of the brown dwarf and white dwarf independently. For the other emission lines, there are no features coming from the white dwarf at those wavelengths therefore only one Gaussian was necessary to model these emission lines.

3 RESULTS

The data from the two nights were analysed together and separately and no significant difference was found between the two. Therefore, we will use the analysis of both nights together.

3.1 Refining the ephemeris

The radial velocities of the white and brown dwarf were extracted from our multiGaussian fits. These were then fitted with a sine function of the form,
\begin{equation} V_{\rm r} = \gamma + K\sin {2\pi \phi }, \end{equation}
(1)
where Vr is the measured radial velocity in km s−1, γ is the systemic velocity of the object in km s−1 and K is the semi-amplitude in km s−1. Phase, ϕ, is calculated by subtracting the zero-point of the ephemeris, T0, from the HJD, t, and dividing by the period, P,
\begin{equation} \phi = \frac{t - T_0}{P}. \end{equation}
(2)
To determine γ, K, P and T0, we used a least squares method to fit equations 1 and 2 to the H α absorption and emission lines in our 78 spectra. Any points more than 3σ away from the model were excluded to improve the fit and 1 km s−1 was added in quadrature to the standard statistical errors until a reduced χ2 of 1 was achieved. The error bars shown on Fig. 1 reflect this and the solid lines are the best fit from the least-squares fit. The period and zero ephemeris were then held constant when fitting the other emission lines and their semi-amplitudes and systemic velocities are listed in Table 1.
Figure 1.

The radial velocity values of the H α emission (red ‘+’) and absorption (black ‘x’). These have been fit with using equation (1) the values from table 1.

Table 1.

The best-fitting radial velocity parameters (cf. equation 1). The period and zero ephemeris were determined by fitting the white-dwarf radial velocity data and can be found in Table 2. P and T0 were then held constant for all subsequent emission lines. The equivalent widths of the emission lines were determined, where possible, from the line profiles.

Wavelength Semi-amplitude Systemic velocity Day-side EW Night-side EW
λ (Å) K (km s−1) γ (km s−1) (Å) (Å)
H α 6562.76 (ab) 28.3 ± 0.4 18.4 ± 0.3
6562.76 (em) −192.4 ± 1.2 4.5 ± 0.9 −0.79 ± 0.021 −0.03 ± 0.021
H β 4861.30 (ab) 28.9 ± 0.9 6.8 ± 0.3
4861.30 (em) −189.5 ± 5.2 2.1 ± 3.7 −0.29 ± 0.030 −0.09 ± 0.030
H i 10938.17 −184.8 ± 8.7 −14.8 ± 6.2
He i 21 655.34 −119.3 ± 50.3 1.7 ± 24.2 −0.66 ± 0.1 0.24 ± 0.1
Na i 5889.95 −188.2 ± 2.4 2.9 ± 1.7 −0.20 ± 0.021 −0.02 ± 0.021
5895.92 −185.7 ± 2.5 6.6 ± 1.8 −0.16 ± 0.021 0.04 ± 0.021
8183.27 −191.2 ± 7.0 1.0 ± 4.8 −0.13 ± 0.021 0.03 ± 0.021
8194.81 −184.8 ± 6.1 −4.4 ± 4.2 −0.19 ± 0.021 0.06 ± 0.021
Mg i 5167.32 −188.1 ± 2.4 −3.5 ± 1.7 −0.15 ± 0.030 −0.05 ± 0.030
5172.68 −186.7 ± 2.3 −3.8 ± 1.6 −0.17 ± 0.030 −0.09 ± 0.030
5183.60 −187.9 ± 1.7 −6.4 ± 1.2 −0.14 ± 0.030 −0.07 ± 0.030
8806.76 −191.4 ± 5.0 −23.0 ± 3.5 −0.40 ± 0.021 −0.03 ± 0.021
Si i 10 585.14 −171.8 ± 12.7 10.9 ± 9.0 −0.32 ± 0.1 0.04 ± 0.1
10 868.79 −182.5 ± 5.0 0.5 ± 3.7 −0.09 ± 0.1 0.23 ± 0.1
12 031.51 −187.1 ± 4.8 −13.2 ± 3.7 −0.29 ± 0.1 0.03 ± 0.1
K i 7698.96 −187.5 ± 2.0 7.1 ± 1.4 −0.14 ± 0.021 0.00 ± 0.021
Ca i 12 816.04 −184.8 ± 6.3 19.8 ± 4.2
Ca ii 3933.66 −185.0 ± 4.2 −1.2 ± 2.9 −0.08 ± 0.030 0.05 ± 0.030
8498.02 −186.5 ± 0.8 6.3 ± 0.6 −0.46 ± 0.021 −0.01 ± 0.021
8542.09 −185.6 ± 1.1 7.1 ± 0.8 −0.75 ± 0.021 0.00 ± 0.021
8662.14 −186.9 ± 1.0 5.1 ± 0.7 −0.73 ± 0.021 −0.05 ± 0.021
Ti i 11 752.31 −177.6 ± 4.4 13.1 ± 3.0 −0.05 ± 0.1 −0.06 ± 0.1
11 982.96 −191.0 ± 3.5 8.2 ± 2.7 −0.14 ± 0.1 −0.05 ± 0.1
19 860.98 −158.7 ± 8.5 −6.2 ± 4.6 −0.44 ± 0.1 0.05 ± 0.1
Fe i 8327.05 −186.3 ± 2.2 3.6 ± 1.6 −0.11 ± 0.021 0.003 ± 0.021
8334.53 −185.7 ± 4.6 19.7 ± 3.2 −0.10 ± 0.021 0.04 ± 0.021
8387.77 −185.2 ± 2.2 7.7 ± 1.6 −0.16 ± 0.021 −0.01 ± 0.021
8688.62 −188.9 ± 1.4 7.2 ± 1.0 −0.15 ± 0.021 0.01 ± 0.021
8824.22 −186.4 ± 2.5 6.6 ± 1.7 −0.05 ± 0.021 −0.01 ± 0.021
10 682.39 −193.9 ± 16.8 40.3 ± 12.1 −0.10 ± 0.1 −0.05 ± 0.1
11 970.50 −186.5 ± 11.6 2.1 ± 7.4 −0.04 ± 0.1 −0.10 ± 0.1
13 120.47 −160.6 ± 27.5 −31.3 ± 17.1 −0.07 ± 0.1 −0.18 ± 0.1
13 147.92 −210.1 ± 10.6 5.0 ± 6.5 −0.07 ± 0.1 0.07 ± 0.1
14 876.54 −219.6 ± 25.4 14.0 ± 15.5 −0.77 ± 0.1 −0.1 ± 0.1
15 022.67 −185.0 ± 5.3 6.7 ± 3.7 −0.13 ± 0.1 0.02 ± 0.1
15 036.47 −168.7 ± 18.0 3.9 ± 12.4 −0.14 ± 0.1 −0.24 ± 0.1
19 771.29 −180.5 ± 15.5 9.0 ± 7.9
Fe ii 8434.49 −198.2 ± 3.9 17.8 ± 2.8 −0.15 ± 0.021 −0.07 ± 0.021
10 690.51 −180.8 ± 4.7 −0.8 ± 3.5 −0.15 ± 0.1 0.02 ± 0.1
11 825.65 −182.2 ± 10.8 18.0 ± 7.2 −0.03 ± 0.1 0.05 ± 0.1
15 885.67 −171.0 ± 15.4 2.1 ± 10.1
17 106.54 −175.0 ± 12.7 10.6 ± 7.8 −0.17 ± 0.1 −0.20 ± 0.1
Wavelength Semi-amplitude Systemic velocity Day-side EW Night-side EW
λ (Å) K (km s−1) γ (km s−1) (Å) (Å)
H α 6562.76 (ab) 28.3 ± 0.4 18.4 ± 0.3
6562.76 (em) −192.4 ± 1.2 4.5 ± 0.9 −0.79 ± 0.021 −0.03 ± 0.021
H β 4861.30 (ab) 28.9 ± 0.9 6.8 ± 0.3
4861.30 (em) −189.5 ± 5.2 2.1 ± 3.7 −0.29 ± 0.030 −0.09 ± 0.030
H i 10938.17 −184.8 ± 8.7 −14.8 ± 6.2
He i 21 655.34 −119.3 ± 50.3 1.7 ± 24.2 −0.66 ± 0.1 0.24 ± 0.1
Na i 5889.95 −188.2 ± 2.4 2.9 ± 1.7 −0.20 ± 0.021 −0.02 ± 0.021
5895.92 −185.7 ± 2.5 6.6 ± 1.8 −0.16 ± 0.021 0.04 ± 0.021
8183.27 −191.2 ± 7.0 1.0 ± 4.8 −0.13 ± 0.021 0.03 ± 0.021
8194.81 −184.8 ± 6.1 −4.4 ± 4.2 −0.19 ± 0.021 0.06 ± 0.021
Mg i 5167.32 −188.1 ± 2.4 −3.5 ± 1.7 −0.15 ± 0.030 −0.05 ± 0.030
5172.68 −186.7 ± 2.3 −3.8 ± 1.6 −0.17 ± 0.030 −0.09 ± 0.030
5183.60 −187.9 ± 1.7 −6.4 ± 1.2 −0.14 ± 0.030 −0.07 ± 0.030
8806.76 −191.4 ± 5.0 −23.0 ± 3.5 −0.40 ± 0.021 −0.03 ± 0.021
Si i 10 585.14 −171.8 ± 12.7 10.9 ± 9.0 −0.32 ± 0.1 0.04 ± 0.1
10 868.79 −182.5 ± 5.0 0.5 ± 3.7 −0.09 ± 0.1 0.23 ± 0.1
12 031.51 −187.1 ± 4.8 −13.2 ± 3.7 −0.29 ± 0.1 0.03 ± 0.1
K i 7698.96 −187.5 ± 2.0 7.1 ± 1.4 −0.14 ± 0.021 0.00 ± 0.021
Ca i 12 816.04 −184.8 ± 6.3 19.8 ± 4.2
Ca ii 3933.66 −185.0 ± 4.2 −1.2 ± 2.9 −0.08 ± 0.030 0.05 ± 0.030
8498.02 −186.5 ± 0.8 6.3 ± 0.6 −0.46 ± 0.021 −0.01 ± 0.021
8542.09 −185.6 ± 1.1 7.1 ± 0.8 −0.75 ± 0.021 0.00 ± 0.021
8662.14 −186.9 ± 1.0 5.1 ± 0.7 −0.73 ± 0.021 −0.05 ± 0.021
Ti i 11 752.31 −177.6 ± 4.4 13.1 ± 3.0 −0.05 ± 0.1 −0.06 ± 0.1
11 982.96 −191.0 ± 3.5 8.2 ± 2.7 −0.14 ± 0.1 −0.05 ± 0.1
19 860.98 −158.7 ± 8.5 −6.2 ± 4.6 −0.44 ± 0.1 0.05 ± 0.1
Fe i 8327.05 −186.3 ± 2.2 3.6 ± 1.6 −0.11 ± 0.021 0.003 ± 0.021
8334.53 −185.7 ± 4.6 19.7 ± 3.2 −0.10 ± 0.021 0.04 ± 0.021
8387.77 −185.2 ± 2.2 7.7 ± 1.6 −0.16 ± 0.021 −0.01 ± 0.021
8688.62 −188.9 ± 1.4 7.2 ± 1.0 −0.15 ± 0.021 0.01 ± 0.021
8824.22 −186.4 ± 2.5 6.6 ± 1.7 −0.05 ± 0.021 −0.01 ± 0.021
10 682.39 −193.9 ± 16.8 40.3 ± 12.1 −0.10 ± 0.1 −0.05 ± 0.1
11 970.50 −186.5 ± 11.6 2.1 ± 7.4 −0.04 ± 0.1 −0.10 ± 0.1
13 120.47 −160.6 ± 27.5 −31.3 ± 17.1 −0.07 ± 0.1 −0.18 ± 0.1
13 147.92 −210.1 ± 10.6 5.0 ± 6.5 −0.07 ± 0.1 0.07 ± 0.1
14 876.54 −219.6 ± 25.4 14.0 ± 15.5 −0.77 ± 0.1 −0.1 ± 0.1
15 022.67 −185.0 ± 5.3 6.7 ± 3.7 −0.13 ± 0.1 0.02 ± 0.1
15 036.47 −168.7 ± 18.0 3.9 ± 12.4 −0.14 ± 0.1 −0.24 ± 0.1
19 771.29 −180.5 ± 15.5 9.0 ± 7.9
Fe ii 8434.49 −198.2 ± 3.9 17.8 ± 2.8 −0.15 ± 0.021 −0.07 ± 0.021
10 690.51 −180.8 ± 4.7 −0.8 ± 3.5 −0.15 ± 0.1 0.02 ± 0.1
11 825.65 −182.2 ± 10.8 18.0 ± 7.2 −0.03 ± 0.1 0.05 ± 0.1
15 885.67 −171.0 ± 15.4 2.1 ± 10.1
17 106.54 −175.0 ± 12.7 10.6 ± 7.8 −0.17 ± 0.1 −0.20 ± 0.1
Table 1.

The best-fitting radial velocity parameters (cf. equation 1). The period and zero ephemeris were determined by fitting the white-dwarf radial velocity data and can be found in Table 2. P and T0 were then held constant for all subsequent emission lines. The equivalent widths of the emission lines were determined, where possible, from the line profiles.

Wavelength Semi-amplitude Systemic velocity Day-side EW Night-side EW
λ (Å) K (km s−1) γ (km s−1) (Å) (Å)
H α 6562.76 (ab) 28.3 ± 0.4 18.4 ± 0.3
6562.76 (em) −192.4 ± 1.2 4.5 ± 0.9 −0.79 ± 0.021 −0.03 ± 0.021
H β 4861.30 (ab) 28.9 ± 0.9 6.8 ± 0.3
4861.30 (em) −189.5 ± 5.2 2.1 ± 3.7 −0.29 ± 0.030 −0.09 ± 0.030
H i 10938.17 −184.8 ± 8.7 −14.8 ± 6.2
He i 21 655.34 −119.3 ± 50.3 1.7 ± 24.2 −0.66 ± 0.1 0.24 ± 0.1
Na i 5889.95 −188.2 ± 2.4 2.9 ± 1.7 −0.20 ± 0.021 −0.02 ± 0.021
5895.92 −185.7 ± 2.5 6.6 ± 1.8 −0.16 ± 0.021 0.04 ± 0.021
8183.27 −191.2 ± 7.0 1.0 ± 4.8 −0.13 ± 0.021 0.03 ± 0.021
8194.81 −184.8 ± 6.1 −4.4 ± 4.2 −0.19 ± 0.021 0.06 ± 0.021
Mg i 5167.32 −188.1 ± 2.4 −3.5 ± 1.7 −0.15 ± 0.030 −0.05 ± 0.030
5172.68 −186.7 ± 2.3 −3.8 ± 1.6 −0.17 ± 0.030 −0.09 ± 0.030
5183.60 −187.9 ± 1.7 −6.4 ± 1.2 −0.14 ± 0.030 −0.07 ± 0.030
8806.76 −191.4 ± 5.0 −23.0 ± 3.5 −0.40 ± 0.021 −0.03 ± 0.021
Si i 10 585.14 −171.8 ± 12.7 10.9 ± 9.0 −0.32 ± 0.1 0.04 ± 0.1
10 868.79 −182.5 ± 5.0 0.5 ± 3.7 −0.09 ± 0.1 0.23 ± 0.1
12 031.51 −187.1 ± 4.8 −13.2 ± 3.7 −0.29 ± 0.1 0.03 ± 0.1
K i 7698.96 −187.5 ± 2.0 7.1 ± 1.4 −0.14 ± 0.021 0.00 ± 0.021
Ca i 12 816.04 −184.8 ± 6.3 19.8 ± 4.2
Ca ii 3933.66 −185.0 ± 4.2 −1.2 ± 2.9 −0.08 ± 0.030 0.05 ± 0.030
8498.02 −186.5 ± 0.8 6.3 ± 0.6 −0.46 ± 0.021 −0.01 ± 0.021
8542.09 −185.6 ± 1.1 7.1 ± 0.8 −0.75 ± 0.021 0.00 ± 0.021
8662.14 −186.9 ± 1.0 5.1 ± 0.7 −0.73 ± 0.021 −0.05 ± 0.021
Ti i 11 752.31 −177.6 ± 4.4 13.1 ± 3.0 −0.05 ± 0.1 −0.06 ± 0.1
11 982.96 −191.0 ± 3.5 8.2 ± 2.7 −0.14 ± 0.1 −0.05 ± 0.1
19 860.98 −158.7 ± 8.5 −6.2 ± 4.6 −0.44 ± 0.1 0.05 ± 0.1
Fe i 8327.05 −186.3 ± 2.2 3.6 ± 1.6 −0.11 ± 0.021 0.003 ± 0.021
8334.53 −185.7 ± 4.6 19.7 ± 3.2 −0.10 ± 0.021 0.04 ± 0.021
8387.77 −185.2 ± 2.2 7.7 ± 1.6 −0.16 ± 0.021 −0.01 ± 0.021
8688.62 −188.9 ± 1.4 7.2 ± 1.0 −0.15 ± 0.021 0.01 ± 0.021
8824.22 −186.4 ± 2.5 6.6 ± 1.7 −0.05 ± 0.021 −0.01 ± 0.021
10 682.39 −193.9 ± 16.8 40.3 ± 12.1 −0.10 ± 0.1 −0.05 ± 0.1
11 970.50 −186.5 ± 11.6 2.1 ± 7.4 −0.04 ± 0.1 −0.10 ± 0.1
13 120.47 −160.6 ± 27.5 −31.3 ± 17.1 −0.07 ± 0.1 −0.18 ± 0.1
13 147.92 −210.1 ± 10.6 5.0 ± 6.5 −0.07 ± 0.1 0.07 ± 0.1
14 876.54 −219.6 ± 25.4 14.0 ± 15.5 −0.77 ± 0.1 −0.1 ± 0.1
15 022.67 −185.0 ± 5.3 6.7 ± 3.7 −0.13 ± 0.1 0.02 ± 0.1
15 036.47 −168.7 ± 18.0 3.9 ± 12.4 −0.14 ± 0.1 −0.24 ± 0.1
19 771.29 −180.5 ± 15.5 9.0 ± 7.9
Fe ii 8434.49 −198.2 ± 3.9 17.8 ± 2.8 −0.15 ± 0.021 −0.07 ± 0.021
10 690.51 −180.8 ± 4.7 −0.8 ± 3.5 −0.15 ± 0.1 0.02 ± 0.1
11 825.65 −182.2 ± 10.8 18.0 ± 7.2 −0.03 ± 0.1 0.05 ± 0.1
15 885.67 −171.0 ± 15.4 2.1 ± 10.1
17 106.54 −175.0 ± 12.7 10.6 ± 7.8 −0.17 ± 0.1 −0.20 ± 0.1
Wavelength Semi-amplitude Systemic velocity Day-side EW Night-side EW
λ (Å) K (km s−1) γ (km s−1) (Å) (Å)
H α 6562.76 (ab) 28.3 ± 0.4 18.4 ± 0.3
6562.76 (em) −192.4 ± 1.2 4.5 ± 0.9 −0.79 ± 0.021 −0.03 ± 0.021
H β 4861.30 (ab) 28.9 ± 0.9 6.8 ± 0.3
4861.30 (em) −189.5 ± 5.2 2.1 ± 3.7 −0.29 ± 0.030 −0.09 ± 0.030
H i 10938.17 −184.8 ± 8.7 −14.8 ± 6.2
He i 21 655.34 −119.3 ± 50.3 1.7 ± 24.2 −0.66 ± 0.1 0.24 ± 0.1
Na i 5889.95 −188.2 ± 2.4 2.9 ± 1.7 −0.20 ± 0.021 −0.02 ± 0.021
5895.92 −185.7 ± 2.5 6.6 ± 1.8 −0.16 ± 0.021 0.04 ± 0.021
8183.27 −191.2 ± 7.0 1.0 ± 4.8 −0.13 ± 0.021 0.03 ± 0.021
8194.81 −184.8 ± 6.1 −4.4 ± 4.2 −0.19 ± 0.021 0.06 ± 0.021
Mg i 5167.32 −188.1 ± 2.4 −3.5 ± 1.7 −0.15 ± 0.030 −0.05 ± 0.030
5172.68 −186.7 ± 2.3 −3.8 ± 1.6 −0.17 ± 0.030 −0.09 ± 0.030
5183.60 −187.9 ± 1.7 −6.4 ± 1.2 −0.14 ± 0.030 −0.07 ± 0.030
8806.76 −191.4 ± 5.0 −23.0 ± 3.5 −0.40 ± 0.021 −0.03 ± 0.021
Si i 10 585.14 −171.8 ± 12.7 10.9 ± 9.0 −0.32 ± 0.1 0.04 ± 0.1
10 868.79 −182.5 ± 5.0 0.5 ± 3.7 −0.09 ± 0.1 0.23 ± 0.1
12 031.51 −187.1 ± 4.8 −13.2 ± 3.7 −0.29 ± 0.1 0.03 ± 0.1
K i 7698.96 −187.5 ± 2.0 7.1 ± 1.4 −0.14 ± 0.021 0.00 ± 0.021
Ca i 12 816.04 −184.8 ± 6.3 19.8 ± 4.2
Ca ii 3933.66 −185.0 ± 4.2 −1.2 ± 2.9 −0.08 ± 0.030 0.05 ± 0.030
8498.02 −186.5 ± 0.8 6.3 ± 0.6 −0.46 ± 0.021 −0.01 ± 0.021
8542.09 −185.6 ± 1.1 7.1 ± 0.8 −0.75 ± 0.021 0.00 ± 0.021
8662.14 −186.9 ± 1.0 5.1 ± 0.7 −0.73 ± 0.021 −0.05 ± 0.021
Ti i 11 752.31 −177.6 ± 4.4 13.1 ± 3.0 −0.05 ± 0.1 −0.06 ± 0.1
11 982.96 −191.0 ± 3.5 8.2 ± 2.7 −0.14 ± 0.1 −0.05 ± 0.1
19 860.98 −158.7 ± 8.5 −6.2 ± 4.6 −0.44 ± 0.1 0.05 ± 0.1
Fe i 8327.05 −186.3 ± 2.2 3.6 ± 1.6 −0.11 ± 0.021 0.003 ± 0.021
8334.53 −185.7 ± 4.6 19.7 ± 3.2 −0.10 ± 0.021 0.04 ± 0.021
8387.77 −185.2 ± 2.2 7.7 ± 1.6 −0.16 ± 0.021 −0.01 ± 0.021
8688.62 −188.9 ± 1.4 7.2 ± 1.0 −0.15 ± 0.021 0.01 ± 0.021
8824.22 −186.4 ± 2.5 6.6 ± 1.7 −0.05 ± 0.021 −0.01 ± 0.021
10 682.39 −193.9 ± 16.8 40.3 ± 12.1 −0.10 ± 0.1 −0.05 ± 0.1
11 970.50 −186.5 ± 11.6 2.1 ± 7.4 −0.04 ± 0.1 −0.10 ± 0.1
13 120.47 −160.6 ± 27.5 −31.3 ± 17.1 −0.07 ± 0.1 −0.18 ± 0.1
13 147.92 −210.1 ± 10.6 5.0 ± 6.5 −0.07 ± 0.1 0.07 ± 0.1
14 876.54 −219.6 ± 25.4 14.0 ± 15.5 −0.77 ± 0.1 −0.1 ± 0.1
15 022.67 −185.0 ± 5.3 6.7 ± 3.7 −0.13 ± 0.1 0.02 ± 0.1
15 036.47 −168.7 ± 18.0 3.9 ± 12.4 −0.14 ± 0.1 −0.24 ± 0.1
19 771.29 −180.5 ± 15.5 9.0 ± 7.9
Fe ii 8434.49 −198.2 ± 3.9 17.8 ± 2.8 −0.15 ± 0.021 −0.07 ± 0.021
10 690.51 −180.8 ± 4.7 −0.8 ± 3.5 −0.15 ± 0.1 0.02 ± 0.1
11 825.65 −182.2 ± 10.8 18.0 ± 7.2 −0.03 ± 0.1 0.05 ± 0.1
15 885.67 −171.0 ± 15.4 2.1 ± 10.1
17 106.54 −175.0 ± 12.7 10.6 ± 7.8 −0.17 ± 0.1 −0.20 ± 0.1

3.2 Spectral features

We have identified several other emission lines: He i, Na i, Mg i, Si i, K i, Ca i and ii, Ti i, and Fe i and ii (see Table 1). The results for each line are very similar, so for simplicity, we will analyse H α and use Ca ii (866.214 nm) for comparison.

Fig. 2 shows the trailed spectra of H α, demonstrating the binary nature of the system. The adjacent trailed spectra of the Ca ii emission line highlight the similarity between the H α and Ca ii emission features. These emission features disappear briefly when the ‘night’-side of the brown dwarf is facing the observer (ϕ = 0.5).1 Therefore, it is likely that the emission is only coming from the heated side of the brown dwarf. This is happening in both trails; however, it is slightly less obvious in the H α trail due to the absorption feature.

Figure 2.

Trailed spectra demonstrating the binary nature of the system. The H α emission line is in anti-phase with the absorption feature. The Ca ii emission closely resembles the H α emission indicating they originate from the same place: the brown dwarf.

In addition to our search for emission features, we also looked for absorption lines that may be present in the system either from the brown dwarf – which has been seen in sdB systems (Wood & Saffer 1999) – or intrinsic to the white dwarf, indicating it is a DAZ. The white dwarf in this system is within the temperature range of 16 000–20 000 K suggested by Kilic et al. (2006) to be cool enough to potentially show metal absorption lines in its spectrum; however, it does not. Only the Balmer lines are present coming from the white dwarf. There is also no emission in phase with the white dwarf therefore no mass transfer or gas circulation can be inferred.

4 COMPARISON TO SIMILAR SYSTEMS

WD0137−349B is unique because it is the first irradiated brown dwarf to exhibit He, Na, Mg, Si, K, Ca, Ti and Fe emission lines in its atmosphere. The radial velocities of these lines confirm that they originate from the brown dwarf and are most likely a direct result of irradiation from the 16 500 K white-dwarf host. Below, we contrast the observable properties of WD0137−349B to other, similarly irradiated brown dwarf and low-mass stellar companions to white dwarfs in short period binaries.

A large number of PCEBs have low-mass stellar companions, and Balmer line emission has been detected in a number of these. For example, LTT 560 displays H α emission with dual components, one of which is due to stellar activity on the M dwarf (Tappert et al. 2007). Ca i and ii, Na i, K i, Fe i and Ti were observed in the secondary's atmosphere as absorption (Tappert et al. 2011). Irradiation induced emission lines outside of the Balmer series are less common but have been seen. For example, Ca ii and Na i emission in GD 448 (Maxted et al. 1998), and He i and Mg ii emission in NN Serpentis (Parsons et al. 2010).

A handful of PCEBs have substellar companions but emission from these is seen less frequently. H α emission has been detected in WD0137−349 before (Maxted et al. 2006); it has also been seen in two other systems: SDSS J120515.80−024222.6 and SDSS J155720.77+091624.6. The former is an eclipsing system with a short period (71.2 min) and a ∼23 500 K white dwarf (Parsons et al. 2017), and the latter has a longer period (2.27 h) and a 21 800 K white dwarf (Farihi et al. 2017).

Whilst WD0137−349B is the first brown-dwarf companion to exhibit irradiation induced metal emission lines, irradiation has been observed to induce temperature differences across the surfaces of brown-dwarf companions. For example, SDSS J141126.20+200911.1 displays a flux excess in the Ks band, indicative of a temperature increase on the irradiated side of the brown dwarf (Littlefair et al. 2014), and the cataclysmic variable SDSS J143317.78+101123.3 shows an ∼ 200 K temperature difference across the surface of the brown dwarf (Hernández Santisteban et al. 2016). WD0137−349 has a temperature difference of 500 K between the day-side and night-side of the brown dwarf (Casewell et al. 2015).

These irradiated white-dwarf–brown-dwarf systems do not all manifest the same observable properties. In order to understand this, we need to look at the atmospheres of isolated brown dwarfs to establish a base from which we can interpret irradiated brown dwarfs.

H α emission as been seen in the atmospheres of ‘hyperactive, brown dwarfs. Pineda et al. (2016) characterized the H α emission in three such objects objects: 2MASS J00361617+1821104 (2MASS 0036), 2MASS J17502484−0016151 (2MASS 1750) and SDSS J042348.57−041403.5 (SDSS 0423). 2MASS 0036 is an L dwarf with variable H α, which is thought to be linked through magnetic interaction with a companion. The H α emission is not as strong in 2MASS 1750, which has been classified as an L5.5 brown dwarf, and SDSS 0423 is an L6+T2 brown dwarf binary with the H α emission coming from the L dwarf. The detection of H α emission is linked to chromospheric activity in brown-dwarf atmospheres. The M9.5 dwarf PC 0025+0447 (M < 0.06 M) discovered by Schneider et al. (1991) displays consistent H α emission. This emission and the emission of He i and O i is thought to be due to coronal activity on the brown dwarf (Martín, Basri & Zapatero Osorio 1999). 2MASS J04183483+2131275, an L5±0.5 brown dwarf, is the first, in the Hyades open star cluster, to display chromospheric H α emission (Pérez-Garrido, Lodieu & Rebolo 2017). This activity is unusual and is not expected from a brown dwarf of this age. A double degenerate binary system, 2MASS J13153094-2649513 (L5 + T7; Burgasser et al. 2011), shows H α, Na i and K i emission coming from the L dwarf (Fuhrmeister, Schmitt & Hauschildt 2005). The lack of irradiation suggests the emission is non-thermal, the cause is thought to be the presence of a chromosphere driven by an unusually strong magnetic field (Burgasser et al. 2011).

WD0137−349B displays metal emission lines akin to irradiated low-mass stellar companions but is the first sub-stellar companion to do so. So-called hyperactive brown dwarfs have shown similar emission that has been taken as evidence for the presence of a chromosphere. The difference being hyperactive brown-dwarf chromospheres are thought to be generated by unusually strong magnetic fields; whereas, WD0137−349B is due to the intense irradiation from the white dwarf. We will discuss the presence of brown-dwarf chromospheres with respect to WD0137−349B in the next section.

5 DISCUSSION

There are some distinct atmospheric differences between WD0137−349 and field brown dwarfs. For example, molecules such as CrH and FeH are expected to be present in an L type dwarf, but we do not detect either of these molecules in our data, nor do we see any sign of Cr, Rb or Cs, which are also expected in a substellar atmosphere. These elements are raised in the atmosphere through the process of convection (Oppenheimer et al. 1998), so their non-detection suggests that irradiation is hindering the convective processes that would normally occur in a brown-dwarf atmosphere.

The metal emission lines that we have detected allow us to draw some parallels, for example, Na i and K i are expected to be in the atmosphere of a late L-type brown dwarf (Kirkpatrick et al. 1999; Helling & Casewell 2014). A brown dwarf of spectral type L7 or later would not be expected to display TiO in its spectrum (Kirkpatrick et al. 1999) and indeed we do not detect any molecules at all. We have however detected Ti i suggesting irradiation may have split this molecule leaving Ti in the upper-atmosphere.

Atmospheres of single or binary brown dwarfs are understood to have very low effective temperatures. The gas-phase is therefore dominated mostly by H2, as well as H2O, CO or CH4, depending on the effective temperature. If the brown dwarf is irradiated by an external source, the atmospheric temperature will rise similar to the outer atmospheres of irradiated planets (see fig. 1 in Barman, Hauschildt & Allard 2004). Schmidt et al. (2015) present observations from the BOSS survey that suggest that brown dwarfs have chromospheres and this is corroborated by Sorahana, Suzuki & Yamamura (2014), who reach a similar conclusion. Rodríguez-Barrera et al. (2015) demonstrate that the upper parts of a brown-dwarf atmosphere can be magnetically coupled (see their Sect.6.1) despite the low atmospheric temperatures that allow the formation of clouds (e.g. Helling & Casewell 2014). Brown dwarfs are known to be convectively active such that magnetic waves could be excited thus causing a heating of the upper, low-density atmosphere; this is similar to what occurs within the Sun (e.g. Mullan & MacDonald (2016)). Schmidt et al. (2015) and Sorahana et al. (2014) both invoke a simplified representation of the chromosphere temperature structure. A chromosphere in WD0137−349B could for example originate from MHD processes like Alfvén wave heating but also be a result of the high-energy irradiation from the white dwarf. We have therefore followed a similar approach to these papers to mimic the presence of a chromosphere. This allows us a simple first test of the expected chemical composition in the hot and diluted upper part of the brown-dwarf atmosphere.

Fig. 3 shows a first test of how ions would dominate the chemical composition in a hot upper atmosphere of WD0137−349B. We apply a gas-phase equilibrium code (see Bilger, Rimmer & Helling 2013) and use a prescribed drift-phoenix model (Tgas, gas)-structure for Teff = 2000 K, log(g) = 5.0, [M/H] = solar as input for the inner atmospheric part. We manipulate the temperature in the upper atmosphere to be indicative of a temperature increase caused by a chromosphere or external energy input (black line, top panel, Fig. 3). The gas phase temperature is now high enough to increase the local degree of thermal ionization (blue line, top panel) substantially as is suggested by the observations presented here. Panels 2–4 in Fig. 3 demonstrate which ions would emerge if WD0137−349B exhibited such a rise in local temperature in the upper atmosphere. Various ions dominate over their atomic counterpart at different gas temperatures i.e. at different atmospheric heights. For example, H+ (H ii) is the dominating H-binding species rather high in our example atmosphere, where Tgas ≈ 6000 K, while Fe ii (Fe+) dominates over Fe i (Fe) at Tgas ≈ 2100 K for the given density structure; the uppermost atmospheric layer is now fully ionized with n(H+) = n(e). While the calculation presented in Fig. 3 does consider only collisional ionization in LTE, it gives an indication of how the upper atmosphere of a brown dwarf would change if irradiation, or another mechanism for forming a chromosphere, caused an outward temperature increase, which is necessary to explain the emission features observed for WD0137−349B.

Figure 3.

Chemical equilibrium results for a drift-phoenix model atmosphere (Teff = 2000 K, log(g) = 5.0, solar metallicity) with an artificially added chromosphere-like temperature increase for atmospheric heights z <104 cm (z decreases outwards). The top panel shows the temperature of the brown dwarf (black, steeper line) due to a chromosphere or external energy input and the degree of thermal ionization (blue, flatter line) that these higher temperatures induce. The second, third and fourth panels show the chemical equilibrium results plotted for a selected number of atoms (solid lines) and ions (dashed lines) each with a corresponding label. As expected, the inner part of the atmosphere is dominated by the atomic form of an element but an outward increasing temperature (e.g. resulting from a chromosphere or corona) will cause different ions to appear height dependently.

Our simple approach suggests that atomic emission lines should be expected to occur in brown-dwarf atmospheres. Irradiated brown dwarfs will, however, have strong temperature gradients in the observable atmosphere similar to irradiated giant gas planets, e.g. HD 189733b (Lee et al. 2016). Such ‘surface’ temperature gradients are caused by hydrodynamic transport of atmospheric gas and will cause the local chemistry to change substantially (Helling et al. 2016). Therefore, a far hotter day-side of an irradiated brown dwarf would show higher ionized elements in emission than the terminator regions, for example. Our data supports this as the strength of the emission lines significantly weaken when the night-side of the brown dwarf is in view. This will be discussed further in the next section.

5.1 Characterizing the line profiles

The equivalent widths (EWs) listed in Table 1 demonstrate that there is a clear and dramatic difference between the day-side and night-side of the brown dwarf that fits with the 500 K temperature difference found by Casewell et al. (2015). Fig. 4 shows the variability in EW of the H α emission line with respect to the phase of the system. EW is effectively a measure of line strength and there is a clear weakening of the line at ϕ = 0.5: where the dark side of the brown dwarf is facing the observer. The same has been performed for the 8662 Å Ca ii line and is also plotted in Fig. 4; the shape of this variation for these two lines matches closely and the Ca ii line appears only marginally weaker. It is likely that the phase dependent variability of the brown dwarf indicates significantly different atmospheric properties between the day-side and night-side. This may give rise to extreme variations in weather (Lee et al. 2015). From Burleigh et al. (2006), we see that the brown dwarf dominates the continuum at wavelengths longwards of 1.95 μm therefore lines detected in the NIR are more largely affected and caution should be exercised when interpreting these EWs.

Figure 4.

Equivalent width of the H α emission (red ‘+’) and the Ca ii 8662 Å line (black ‘x’) as a function of phase. The data have been duplicated over two phases.

We have compared our measured EWs to published EWs of two ‘hyperactive’ brown-dwarf systems (2MASS J0418 and SDSS J0423) that display H α emission. There are other brown dwarfs with published EWs for this line; however, they are most often early L dwarfs on or around the hydrogen burning mass limit, or the H α emission is sporadic, probably due to flaring. The two included in Fig. 5 have sustained H α emission. We also compare to V471 Tau as it is the only other irradiated companion, to our knowledge, whose published EWs display a clear phase dependency. Comparisons to the hyperactive brown dwarfs is valid as the underlying continuum is unaffected by non-brown-dwarf components such as reflection from a companion. Setting the continuum level in V471 Tau was more complex due to line blending and rotational broadening but Bois, Mochnacki & Lanning (1991) accounted for this. Thus, we can compare our H α values to these literature values. Data are sparse thus any conclusions drawn are preliminary and more observations are required to build a better picture of irradiated atmospheres. That said, Fig. 5 suggests that the emission from WD0137−349B is more akin to irradiation-induced emission like that seen on V471 Tau than emission seen in ‘hyper-active’ brown dwarfs, which is driven by unusually strong magnetic fields.

Figure 5.

A comparison of our measured equivalent widths to other brown dwarfs that display this emission. V471 Tau has also been included as it has day- (red, lower) and night- (blue, higher) side equivalent widths with a clear phase dependence similar to WD0137−349B.

We also observe a phase-dependent change in the width of the emission lines as is demonstrated in Fig. 6. From this we can see that the Ca ii emission line is narrower at nearly every point in phase suggesting it may be coming from a region lower in the brown-dwarf atmosphere. The signal-to-noise ratio of the Ca ii line between ϕ = 0.35 and 0.7 was too low to fit.

Figure 6.

The FWHM of the H α line (red ‘+’) and the Ca ii line (black ‘x’) with respect to phase. The horizontal lines indicate the FWHM calculated from the vsin (i) for each line in the corresponding colour. This has been duplicated over two orbits for clarity.

To investigate this further, we have simulated the Ca ii emission lines using a model that takes into account the Roche geometry of the brown dwarf and limb darkening / brightening. We simulate an intrinsic spectrum at each point on the brown dwarf that is equal to Gaussians with the same full width at half-maximum (FWHM) as the resolution of the spectra. The emission flux is proportional to the incident flux from the white dwarf and we assume no intrinsic emission or absorption from the brown dwarf i.e. pure emission due to irradiation. The exposure time of 5 min is an appreciable fraction of the orbital period; to account for this, we have included phase smearing by numerical integration over three phase points.

We set up a grid of semi-amplitude, K, filling factor, f, and linear limb darkening, X, as defined in Maxted et al. (1998). We use Levenburg–Marquadt minimization to optimize γ and the ratio of the Ca ii intrinsic line strengths plus an overall scaling factor to match the line strengths of the observations. From Fig. 7, we can see that these simulations provide a good fit to the data, the best-fitting values are listed in Table 2 along with our calculated mass ratio. We also estimate the v sin (i) of this system to be 39.3 km s−1, which can be translated into FWHM using FWHM |$= \frac{v\,\sin\, (i)}{c} \lambda$|⁠, where λ is the wavelength of the line. This has been indicated with horizontal lines in Fig. 6 for H α and Ca ii.

Figure 7.

Trail plots of the Ca ii data and our simulated Ca ii 8662 Å emission line for comparison.

Table 2.

The values obtained from fitting equation (1) to the measured radial velocities of the H α absorption and emission lines. The error on the final digit of the period is shown in parentheses. We list the corrected semi-amplitude for the brown dwarf used to correct for the light from the emission being offset from the centre of mass of the brown dwarf; we used this value to calculate the masses of the white dwarf and brown dwarf.

Properties of the system
P (days) 0.079 429 939(1)
T0 2454 178.6761 ± 0.0003
Kc (km s−1) 210 ± 5
Mass ratio 0.135 ± 0.004
Properties of the system
P (days) 0.079 429 939(1)
T0 2454 178.6761 ± 0.0003
Kc (km s−1) 210 ± 5
Mass ratio 0.135 ± 0.004
Table 2.

The values obtained from fitting equation (1) to the measured radial velocities of the H α absorption and emission lines. The error on the final digit of the period is shown in parentheses. We list the corrected semi-amplitude for the brown dwarf used to correct for the light from the emission being offset from the centre of mass of the brown dwarf; we used this value to calculate the masses of the white dwarf and brown dwarf.

Properties of the system
P (days) 0.079 429 939(1)
T0 2454 178.6761 ± 0.0003
Kc (km s−1) 210 ± 5
Mass ratio 0.135 ± 0.004
Properties of the system
P (days) 0.079 429 939(1)
T0 2454 178.6761 ± 0.0003
Kc (km s−1) 210 ± 5
Mass ratio 0.135 ± 0.004

The work presented in this paper demonstrates the first brown-dwarf atmosphere to exhibit a range of emission lines as a direct result of intense UV irradiation from the 16 500 K white-dwarf primary. Ability to reproduce the line strength variation that is seen in WD0137−349 will be important in any future models of irradiated atmospheres.

5.2 Comparison to hot Jupiter models

In an attempt to explore the role irradiation and opacity play in heat redistribution, Perna, Heng & Pont (2012) model the atmosphere of hot Jupiters using 3D atmospheric circulation models. The models are cloud-free and this results in an atmosphere opacity that is too low. They do find however that heat redistribution breaks down at irradiation temperatures greater than 2200–2400 K and attribute this to the advective time-scale becoming much longer than the radiative time-scale. Showman et al. (2013) use a circulation model to compare the atmospheres of hot Jupiters with and without intense irradiation. They found that strong insolation (stellar radiation that reaches the planet's surface) produces circulation dominated by high-altitude air flow from the day-side to the night-side and suggest a return flow at deeper levels. It is clear that irradiation of gas-giant planets introduces additional dynamical effects within their atmospheres. Younger, lower gravity (log g ≈ 3) brown dwarfs share many properties with these gas giant exoplanets (Faherty et al. 2013).

No 3D simulations have been created for higher gravity irradiated brown-dwarfs-like WD0137−349B, but we may expect comparable atmospheric effects. A recent study by Hernández Santisteban et al. (2016) looks at the interacting binary system of a brown dwarf being irradiated by a 13 200 K white dwarf. They model the atmosphere using a simple geometric reprocessing model and find poor heat redistribution from the day-side to the night-side. With a slightly hotter white dwarf such as WD0137−349 (Teff = 16 500 K), one would expect this effect to be similar and the strength variation over the orbit of all our observed emission lines does support this. Although the longer period of WD0137−349B is likely to make heat redistribution more effective. Full analysis of phase curves would need to be done to be sure of this effect.

6 CONCLUSION

We have used 78 XSHOOTER observations of WD0137−349AB to refine the ephemeris of the system using the emission and absorption lines of H α. We show that this system is detached as we see no indication of mass transfer in the trailed spectra and no lines other than Balmer lines coming from the white dwarf. We have detected new emission lines of He i, Na i, Mg i, Si i, K i, Ca i and ii, Ti i, and Fe i and ii originating from the brown-dwarf atmosphere as a direct result of irradiation.

WD0137−349 is unique as it is the first brown dwarf to exhibit irradiation induced metal emission lines. One would expect that shorter orbits with hotter white dwarfs, thus more intense irradiation, would lead to stronger Balmer line emission and metal emission like we see in the atmosphere of WD0137−349B. We have discussed comparable systems and find that despite their similarities, they do not all manifest the same observable properties.

To establish what to expect from an irradiated brown dwarf atmosphere, we use a drift-phoenix model and manipulate the temperature of the upper atmosphere to give an outward temperature increase. This provides a simple preliminary test of the expected chemical composition of a brown-dwarf atmosphere at different gas temperatures or atmospheric heights. The elements we have detected and their ionization states loosely match the predictions of this model, which supports the possibility of a brown-dwarf chromosphere being present.

We measured the equivalent widths of the emission lines in WD0137−349B and found the line strengths to weaken from the day-side to the night-side of the brown dwarf. This indicates that the irradiation is likely driving dramatic temperature and weather variations between the two sides of this tidally locked brown dwarf. We found the H α emission line to have a strong phase dependence and compared our H α day-side and night-side equivalent width to that of similar systems. We find that WD0137−349B has equivalent widths more akin to the irradiation induced H α emission of low-mass stellar companion atmospheres.

Any future models of irradiated atmospheres will need to reproduce the line strength variation that is seen here.

Acknowledgments

We would like to thank our helpful reviewer for their comments and suggestions. We also thank Tom Marsh for the use of molly. ESL acknowledges the support of STFC studentship. SLC acknowledges support from the University of Leicester College of Science and Engineering. CH highlights the financial support of the European community under the FP7 ERC starting grant 257431. This work was supported by the Science and Technology Facilities Council [ST/M001040/1]. This work is based on observations made with ESO telescopes at La Silla Paranal Observatory under programme ID 093.C-0211(A).

1

For consistency, we follow previous studies of WD0137−349 in adopting the unconventional phasing of this system and define phase 0 as when the brown dwarf is furthest from the observer.

REFERENCES

Armitage
P. J.
,
2010
,
Astrophysics of Planet Formation

Barman
T. S.
,
Hauschildt
P. H.
,
Allard
F.
,
2004
,
ApJ
,
614
,
338

Beuermann
K.
et al. ,
2013
,
A&A
,
558
,
A96

Bilger
C.
,
Rimmer
P.
,
Helling
C.
,
2013
,
MNRAS
,
435
,
1888

Bois
B.
,
Mochnacki
S. W.
,
Lanning
H. H.
,
1991
,
AJ
,
102
,
2079

Boss
A. P.
et al. ,
2007
,
Trans. IAU A
,
26
,
183

Burgasser
A. J.
Reid
I. N.
Siegler
N.
Close
L.
Allen
P.
Lowrance
P.
Gizis
J.
,
2007
,
Protostars and Planets V
, p.
427

Burgasser
A. J.
,
Sitarski
B. N.
,
Gelino
C. R.
,
Logsdon
S. E.
,
Perrin
M. D.
,
2011
,
ApJ
,
739
,
49

Burleigh
M. R.
et al. ,
2006
,
MNRAS
,
373
,
1416

Burleigh
M. R.
et al. ,
2011
, in
Schuh
S.
Drechsel
H.
Heber
U.
, eds,
AIP Conf. Proc. Vol. 1331., Planetary Systems Beyond The Main Sequence
.
Am. Inst. Phys.
,
New York
, p.
262

Burrows
A.
,
Liebert
J.
,
1993
,
Rev. Mod. Phys.
,
65
,
301

Casewell
S. L.
et al. ,
2012
,
ApJ
,
759
,
L34

Casewell
S. L.
et al. ,
2015
,
MNRAS
,
447
,
3218

Chabrier
G.
,
Baraffe
I.
,
1997
,
A&A
,
327
,
1039

Cushing
M. C.
,
Rayner
J. T.
,
Vacca
W. D.
,
2005
,
ApJ
,
623
,
1115

Debes
J. H.
,
Hoard
D. W.
,
Wachter
S.
,
Leisawitz
D. T.
,
Cohen
M.
,
2011
,
ApJS
,
197
,
38

Dobbie
P. D.
,
Burleigh
M. R.
,
Levan
A. J.
,
Barstow
M. A.
,
Napiwotzki
R.
,
Hubeny
I.
,
2005
,
A&A
,
439
,
1159

Faherty
J. K.
,
Rice
E. L.
,
Cruz
K. L.
,
Mamajek
E. E.
,
Núñez
A.
,
2013
,
AJ
,
145
,
2

Farihi
J.
,
Christopher
M.
,
2004
,
AJ
,
128
,
1868

Farihi
J.
,
Parsons
S. G.
,
Gänsicke
B. T.
,
2017
,
Nature Astron.
,
1
,
0032

Fuhrmeister
B.
,
Schmitt
J. H. M. M.
,
Hauschildt
P. H.
,
2005
,
A&A
,
439
,
1137

Girven
J.
,
Gänsicke
B. T.
,
Steeghs
D.
,
Koester
D.
,
2011
,
MNRAS
,
417
,
1210

Grether
D.
,
Lineweaver
C. H.
,
2007
,
ApJ
,
669
,
1220

Helling
C.
,
Casewell
S.
,
2014
,
A&AR
,
22
,
80

Helling
C.
,
Fomins
A.
,
2013
,
Phil. Trans. R. Soc. A
,
371
,
10581

Helling
C.
et al. ,
2016
,
MNRAS
,
460
,
855

Hernández Santisteban
J. V.
et al. ,
2016
,
Nature
,
533
,
366

Izzo
C.
,
Kornweibel
N.
,
McKay
D.
,
Palsa
R.
,
Peron
M.
,
Taylor
M.
,
2004
,
The Messenger
,
117
,
33

Kilic
M.
,
von Hippel
T.
,
Leggett
S. K.
,
Winget
D. E.
,
2006
,
ApJ
,
646
,
474

Kirkpatrick
J. D.
et al. ,
1999
,
ApJ
,
519
,
802

Lee
G.
,
Helling
C.
,
Dobbs-Dixon
I.
,
Juncher
D.
,
2015
,
A&A
,
580
,
A12

Lee
G.
,
Dobbs-Dixon
I.
,
Helling
C.
,
Bognar
K.
,
Woitke
P.
,
2016
,
A&A
,
594
,
A48

Littlefair
S. P.
,
Dhillon
V. S.
,
Martín
E. L.
,
2003
,
MNRAS
,
340
,
264

Littlefair
S. P.
et al. ,
2014
,
MNRAS
,
445
,
2106

Marley
M. S.
,
Robinson
T. D.
,
2015
,
ARA&A
,
53
,
279

Martín
E. L.
,
Basri
G.
,
Zapatero Osorio
M. R.
,
1999
,
AJ
,
118
,
1005

Maxted
P. F. L.
,
Marsh
T. R.
,
Moran
C.
,
Dhillon
V. S.
,
Hilditch
R. W.
,
1998
,
MNRAS
,
300
,
1225

Maxted
P. F. L.
,
Napiwotzki
R.
,
Dobbie
P. D.
,
Burleigh
M. R.
,
2006
,
Nature
,
442
,
543

McCook
G. P.
,
Sion
E. M.
,
1999
,
VizieR Online Data Catalog
,
3210
,
0

Metchev
S. A.
,
2006
,
PhD thesis
,
California Institute of Technology
,
California

Mullan
D. J.
,
MacDonald
J.
,
2016
,
ApJ
,
818
,
154

Oppenheimer
B. R.
,
Kulkarni
S. R.
,
Matthews
K.
,
van Kerkwijk
M. H.
,
1998
,
ApJ
,
502
,
932

Parsons
S. G.
,
Marsh
T. R.
,
Copperwheat
C. M.
,
Dhillon
V. S.
,
Littlefair
S. P.
,
Gänsicke
B. T.
,
Hickman
R.
,
2010
,
MNRAS
,
402
,
2591

Parsons
S. G.
et al. ,
2017
,
MNRAS
,
preprint (arXiv:1705.05856)

Pérez-Garrido
A.
,
Lodieu
N.
,
Rebolo
R.
,
2017
,
A&A
,
599
,
A78

Perna
R.
,
Heng
K.
,
Pont
F.
,
2012
,
ApJ
,
751
,
59

Pineda
J. S.
,
Hallinan
G.
,
Kirkpatrick
J. D.
,
Cotter
G.
,
Kao
M. M.
,
Mooley
K.
,
2016
,
ApJ
,
826
,
73

Politano
M.
,
2004
,
ApJ
,
604
,
817

Rodríguez-Barrera
M. I.
,
Helling
C.
,
Stark
C. R.
,
Rice
A. M.
,
2015
,
MNRAS
,
454
,
3977

Schmidt
S. J.
,
Hawley
S. L.
,
West
A. A.
,
Bochanski
J. J.
,
Davenport
J. R. A.
,
Ge
J.
,
Schneider
D. P.
,
2015
,
AJ
,
149
,
158

Schneider
D. P.
,
Greenstein
J. L.
,
Schmidt
M.
,
Gunn
J. E.
,
1991
,
AJ
,
102
,
1180

Showman
A. P.
,
Fortney
J. J.
,
Lewis
N. K.
,
Shabram
M.
,
2013
,
ApJ
,
762
,
24

Sorahana
S.
,
Suzuki
T. K.
,
Yamamura
I.
,
2014
,
MNRAS
,
440
,
3675

Spiegel
D. S.
,
Burrows
A.
,
Milsom
J. A.
,
2011
,
ApJ
,
727
,
57

Steele
P. R.
,
Burleigh
M. R.
,
Dobbie
P. D.
,
Jameson
R. F.
,
Barstow
M. A.
,
Satterthwaite
R. P.
,
2011
,
MNRAS
,
416
,
2768

Steele
P. R.
et al. ,
2013
,
MNRAS
,
429
,
3492

Tappert
C.
,
Gänsicke
B. T.
,
Schmidtobreick
L.
,
Aungwerojwit
A.
,
Mennickent
R. E.
,
Koester
D.
,
2007
,
A&A
,
474
,
205

Tappert
C.
,
Gänsicke
B. T.
,
Schmidtobreick
L.
,
Ribeiro
T.
,
2011
,
A&A
,
532
,
A129

Vernet
J.
et al. ,
2011
,
A&A
,
536
,
A105

Wood
J. H.
,
Saffer
R.
,
1999
,
MNRAS
,
305
,
820