ALMA’s view of the Sun’s nearest neighbours

ALMA’s view of the Sun’s nearest neighbours

The submm/mm SEDs of the Centauri binary and a new source
R. Liseau Department of Earth and Space Sciences, Chalmers University of Technology, Onsala Space Observatory, SE-439 92 Onsala, Sweden,    V. De la Luz SCiESMEX, Instituto de Geofisica, Unidad Michoacan, Universidad Nacional Autonoma de Mexico, Morelia, Michoacan, Mexico    E. O’Gorman Dublin Institute for Advanced Studies, Astronomy and Astrophysics Secton, 31 Fitzwilliam Place Dublin 2, Ireland    E. Bertone Instituto Nacional de Astrofísica, Óptica y Electrónica (INAOE), Luis Enrique Erro 1, Sta. María Tonantzintla, Puebla, Mexico    M. Chavez Instituto Nacional de Astrofísica, Óptica y Electrónica (INAOE), Luis Enrique Erro 1, Sta. María Tonantzintla, Puebla, Mexico    F. Tapia SCiESMEX, Instituto de Geofisica, Unidad Michoacan, Universidad Nacional Autonoma de Mexico, Morelia, Michoacan, Mexico
Received ; accepted
Key Words.:
stars: chromospheres – stars: solar-type – (stars:) binaries: general – stars: individual: tauri AB – submillimeter: stars – radio continuum: stars

Context:The precise mechanisms that provide the non-radiative energy for heating the chromosphere and the corona of the Sun and other stars are at the focus of intense contemporary research.

Aims:Observations at submm/mm wavelengths are particularly useful to obtain information about the run of the temperature in the upper atmosphere of Sun-like stars. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to study the chromospheric emission of the  Centauri binary system in all six available frequency bands during Cycle 2 in 2014-2015.

Methods:Since ALMA is an interferometer, the multi-telescope array is particularly suited for the observation of point sources. With its large collecting area, the sensitivity is high enough to allow the observation of nearby main-sequence stars at submm/mm wavelengths for the first time. The comparison of the observed spectral energy distributions with theoretical model computations provides the chromospheric structure in terms of temperature and density above the stellar photosphere and the quantitative understanding of the primary emission processes.

Results:Both stars were detected and resolved at all ALMA frequencies. For both  and B, the existence and location of the temperature minima, firstly detected from space with Herschel, are well reproduced by the theoretical models of this paper. For , the temperature minimum is deeper than for A and occurs at a lower height in the atmosphere, but for both stars, is consistently lower than what is derived from optical and UV data. In addition, and as a completely different matter, a third point source was detected in Band 8 (405 GHz, 740 m) in 2015. With only one epoch and only one detection, we are left with little information regarding that object’s nature, but conjecture that it might be a distant solar system object.

Conclusions:The submm/mm emission of the  stars is indeed very well reproduced by modified chromospheric models of the Quiet Sun. This most likely means that the non-radiative heating mechanisms of the upper atmosphere that are at work in the Sun are operating also in other solar-type stars.

1 Introduction

Outside the solar system, Alpha Centauri () is our nearest neighbour, only a little more than a parsec away (0742). It is a double star, and its primary  has the same spectral type and luminosity class as the Sun, viz. G2 V. The secondary, , is a somewhat cooler star, of spectral type K1 V. Using asteroseismology, the age of the main-sequence stars  and B has been determined to  Gyr by Thévenin et al. (2002), whereas statistical methods resulted in estimates of 8 to 10 Gyr, depending on the method used, the Ca II  index or the X-ray luminosity, respectively (see, e.g., Eiroa et al., 2013, and references therein).

The proximity of , the similarity of A, and the differences of B, compared to the Sun provide an excellent opportunity to study the stellar-solar relationship, as the understanding of the physics of the Sun and the stars is an iterative process that provides feed-back in both directions. For instance, an outstanding problem of modern solar physics is the heating of the outer atmospheric layers, i.e., of the chromosphere and the corona (Wedemeyer-Böhm et al., 2007). A few hundred kilometers above the solar photosphere, the temperature gradient changes sign at the location of the temperature minimum. From early theoretical models of the chromosphere, this phenomenon was already found also for  and B (and in addition, for  Boo and  CMi: Ayres et al., 1976). The primary observables were the wings of optical and UV resonance lines, e.g. Ca II H&K and Mg II h&k, the cores of which are formed higher up in the chromosphere. In addition, high temperature tracers also include high ionization lines and continua in the UV from the transition region and radio emission and X-rays from the corona.

Date Start UTC End UTC Synthesized Beam
yyyy-mm-dd hh min sec hh min sec hh mm ss.s              hh mm ss.s                    
B3 2014-07-03 00 47 20.4 01 38 19.4 14  39  28.893   14  39  28.333      19
B7 2014-07-07 02 26 26.4 02 44 53.8 14  39  28.883   14  39  28.325      47
B9 2014-07-18 00 56 05.7 01 26 49.4 14  39  28.870   14  39  28.309      36
B6 2014-12-16 11 04 36.6 11 18 34.2 14  39  28.650   14  39  28.120      71
B4 2015-01-18 13 35 24.5 13 59 40.8 14  39  28.624   14  39  28.110      82
B8 2015-05-02 03 04 14.2 03 25 01.7 14  39  28.439   14  39  27.934     
Table 1: Positions of  and B with ALMA in Right Ascension and Declination (ICRS J 2000.0)
Primary beam corrected flux density, (mJy), and signal-to-noise [S/N]
Band 9 Band 8 Band 7 Band 6 Band 4 Band 3
679 GHz 405 GHz 343.5 GHz 233 GHz 145 GHz 97.5 GHz
442 m 740 m 873 m 1287 m 2068 m 3075 m
A [71] 1 [168] [137] [170] [83] [281]
B 1    [13]      [87] [34] 1 [124] [34]      [80]
Table 2: ALMA flux density data for the tauri binary

The temperature minimum of  was directly observed in the far-infrared spectral energy distribution (SED) by Liseau et al. (2013). However, the far-infrared data did not resolve the binary in its individual components and the interpretation had to rely on photometry at shorter wavelengths. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) at three frequencies finally resolved the pair and the individual SEDs were spectrally mapped throughout the sub-millimeter (submm), up to 3 mm (Liseau et al., 2015).  was observed with three more ALMA bands during Cycle 2. The stars themselves were unresolved and appeared as point sources to ALMA. With regard to the stellar-solar connection, these observations would refer to analogues of the Quiet Sun, for which the intensity is integrated over the solar disk.

The metallicity of  is slightly higher than that of the Sun, i.e. (Torres et al., 2010), a fact that could favor the existence of planets around the stars (e.g., Wang & Fischer, 2015). Examining a wealth of radial velocity data, Dumusque et al. (2012) announced the discovery of an Earth-mass planet around . That was however challenged by Hatzes (2013), Demory et al. (2015) and Rajpaul et al. (2016) who were unable to confirm the existence of this object.

Attempts to detect planets around  with direct imaging in the optical and the near infrared have hitherto been unsuccessful, see Kervella et al. (2006); Kervella & Thévenin (2007) and Kervella et al. (2016, in preparation). At these wavelengths, any feeble planetary signal within several arcseconds from the stars would be totally swamped by their overwhelming glare (-magnitude = ), alternatively be hidden behind the coronagraphic mask inside the inner working angle. This contrast problem would be naturally overcome for closeby faint objects with ALMA, an interferometer that for point sources in the reconstructed images generates a much cleaner point spread function (PSF), and our imaging results of  with ALMA are discussed toward the end of this paper.

The organization of this paper is as follows: Sect. 2 reports the observations and the data reduction. Sect. 3 briefly presents the results, which are discussed in Sect. 4. We round off with our conclusions in Sect. 5.

Figure 1: Measurements of the flux density of  (blue circles) and of  (red circles) with ALMA, with statistical error bars inside the symbols. Left: Assuming that , least-square fits to the Band 3 to 9 flux densities are shown by dashed lines, with the power law exponent shown next to them. Right: Shown by the solid lines are fits, performed as in the left panel, to the data above, and by dotted lines below 200 GHz ( 1.5 mm). The ALMA bands, with their central wavelengths, are identified at the bottom of the figure.
Figure 2: Brightness temperature in Kelvin at ALMA wavelengths in m, for Bands 3 to 9 of the G-star  (left, blue) and the K-star  (right, red). In addition to the observational rms-errors (solid bars), the estimated absolute errors, including calibration uncertaities, are shown as dashed error bars. The stellar photospheres are represented by extrapolations to PHOENIX model atmospheres of Brott & Hauschildt (2005) for the stars’ respective (, , [Fe/H]) and are shown as black dashed lines. The ALMA bands are indicated below. A solar model chromosphere (VAL IIIC, Vernazza et al., 1981) is shown as long dashes, with data for the Sun from Loukitcheva et al. (2004) as black open circles.

2 Observations and data reduction

The binary B was observed in all six ALMA continuum bands during the period July 2014 to May 2015 (Table 2). The field of view (primary beam) varied from about 10 for the shortest wavelength to about 1 for the longest. Similarly, the angular resolution (synthesized half power beam width) ranged from 02 to 15. With angular diameters of 0008 and 0006 for A and B at 2 (Kervella et al., 2003), the stars were-point like to the ALMA interferometer in all wave-bands (cf. Table 1). The ALMA program code is 2013.1.00170.S and the observations in Band 3, 7 and 9 have already been described in detail by Liseau et al. (2015) and will not be repeated here.

The observations in Band 4, 6, and 8 were taken in the standard wideband continuum mode with 8 GHz effective bandwidth spread over four spectral windows in each of the bands. The Band 4 observations, taken on 2015 Jan 18 with 34 antennas, were centered on 145 GHz (2068 m), with  24 min of observing time with 5.5 min on-source. The Band 6 observations, taken on 2014 Dec 16 with 35 antennas, were centered on 233 GHz (1287 m), with  14 min of observing time with  2 min on-source. Finally, the Band 8 observations, taken on 2015 May 2 with 37 antennas, were centered on 405 GHz (740 m), with  21 min of observing time with  7 min on-source.

The visibilities were flagged and calibrated following standard procedures using the CASA package111CASA is an acronym for Common Astronomy Software Application. v4.2.2 for Band 4 and 6, and v4.3.1 for Band 8. The quasar J1617-5848 was used as complex gain calibrator in Band 4 and 8, while J1408-5712 was used in Band 6. The quasar J1427-4206 was used as bandpass calibrator in Band 6 and 8, while J1617-5848 was used in Band 4. Flux calibration was done using Ceres in Band 4 when at 74 elevation, while  was at 44. The quasar 1427-421 was used for flux calibration in Band 6 when it was at 57 elevation and  was at 46, while Titan was used in Band 8 when it was at 50 elevation and  was at 52.

Imaging was performed using natural weighting in Band 4, 6, and 8 with one round of phase-only self-calibration was carried out on all three images to improve the rms noise. The synthesized beam sizes are listed in Table 1 and the primary beam-corrected flux densities and the rms noise per synthesized beam in the pointing center are listed in Table 2. We also imaged the ALMA spectral windows separately in each band to assess the spectral index within each band and the resultant flux densities for  and B are listed in Table A.1. and A.2., respectively.

(m) (GHz) in-band in-band
9 1442 679
8 1740 405
7 1873 1343.5
6 1287 233
4 2068 145
3 3075 1197.5
Table 3: Stellar flux ratios and in-band (spw 1 - spw 4) spectral indices. Too small bandwidth or too large errors.

3 Results

The binary system is well resolved at all frequencies. The J2000-coordinates for  and B on the observational dates are presented in Table 1, together with the sizes of the synthesized beams (ellipses with semi-major axes and semi-minor axes in arcseconds) and their orientations (position angle in degrees). The frequencies of the bands are given in Table 2, where the primary beam corrected flux densities, , are reported together with their statistical errors. As can be seen, the signal-to-noise ratio, S/N, spans the range 10–100 for , and excels to nearly 300 for . The absolute flux calibration is quoted in terms of goals222, viz. better than 5% for bands B 3 and B 4, better than 10% for B 6 and B 7, and at best about 20% for B 8 and B 9. These goals are shown for  and B in Fig. 2.

3.1 Relative fluxes from 0.4 to 3.1 mm

The average flux ratio for the binary over the ALMA bands 3 through 9 is (Table 3). This would be close to the ratio of their respective solid angles , where the radii are those of their interferometrically measured photospheric disks of uniform brightness (Kervella et al., 2003). Comparison with the value for the range 0.09 m to 70 m, i.e., (Liseau et al., 2013), indicates an apparently remarkable constancy of the flux ratio over four orders of magnitude in wavelength, from the photospheric emission in the visible to that in the micro-wave regime.

3.2 Spectral slopes of the SEDs

A first order characterization of the emission mechanism(s) can be obtained from the spectral slope of the logarithmic SED. Assuming that , linear regression (Press et al., 1986)333, and . to the Band 3 to 9 data results in a spectral index with a for . For , the corresponding and , see Fig. 1. The goodness-of-fit is for both.

This apparent constancy of the slope close to a value of two over the entire ALMA range, from 0.4 to 3.1 mm, is perhaps surprising. A more careful inspection of the data reveals that the slopes at the shorter wavelengths appear marginally steeper, but that the long-wavelength data, not totally unexpected, seem to flatten out. Dividing the data into two sub-sets for both stars, i.e. below and above 1.5 mm (200 GHz), yields for the spectral indices of the -SED and . Similarly, for , and (Fig. 1). In these cases, the formal fit errors are considerably larger for both  and B. However with regard to the fits in the left panel, the observed Band 3 flux densities are in excess by more than for A and by more than for B. Therefore, the flattening of the SEDs towards lower frequencies is real.

Observations at longer wavelengths would help to better constrain the run of the SED. Unfortunately, at declination south of  the number of sensitive observing facilities is limited. Trigilio et al. (2013) and (2014) proposed Australia Telescope Compact Array (ATCA) observations at 15 mm (17 GHz) and 16 cm (2 GHz). C. Trigilio privately communicated to us that both stars were recently detected at 17 GHz. However, having no further information, we provide here our own flux estimates for ATCA observations of the binary . These are based on extrapolations beyond ALMA-Band 3 and the sensitivity specifications of the 6 km compact array for the K-band (15 mm) and C/X-band (4 cm)444, resulting in estimates of the (0.27 mJy) and 13 (0.13 mJy) for  and (0.04 mJy) and 7 (0.02 mJy) for , respectively. These values refer to 12 hour on-source integrations (rms = 0.003 mJy). The corresponding brightness temperatures are shown below, in Fig. 4.

Spectral indices for flux integrations over the individual bands are shown in Table 3, except for Band 9, where the fractional bandwidth is too small for meaningful measurement, and for Band 7, where the relative errors are too large (negative slope within the band). Inside the individual bands, the data were collected through four spectral windows (spw; see Fig. 7), with the flux data for these provided in Appendix A.

Band (m) (GHz) (mJy) (mJy) (mJy) (km) (K)
3 3075 1197.5 11 1 1 2143:
4 2068 145 11 1 1 2140:
6 1287 233 1 1 1180
7 1873 1343.5 1 1 1965
8 1740 405 1 1 1950
9 1442 679 1050
Table 4: Brightness temperatures and chromospheric heights of

4 Discussion

4.1 The stellar brightness temperatures

The direct observation of the temperature minima of B at far infrared wavelengths indicated a clear kinship with the Sun’s chromosphere (Liseau et al., 2013, 2015). At these wavelengths, the continuum opacity is dominated by inverse bremsstrahlung, with some contribution due to free-free H processes (e.g., Dulk, 1985; Wedemeyer et al., 2015).

Fig. 2 displays the observed spectral energy distributions of both stars in terms of their brightness temperatures555The brightness temperature, or radiation temperature, is the temperature of a blackbody that emits the same amount of radiation as the observed flux at a given frequency.


where is the stellar radius, the height at which the observed radiation originates, is the radiation frequency, is the distance to the source, is the observed flux density, and the other symbols have their usual meaning.

For the Sun, , where refers to the height above the solar photosphere, where the optical depth in the visual and . We assume similar -values for the  stars and use their photospheric radii, i.e. , where refers to the values determined by Kervella et al. (2003). When (Rayleigh-Jeans regime), Eq. 1 simplifies to


Consequently, in the Rayleigh-Jeans regime (RJ), optically thick free-free emission (or Bremsstrahlung) will behave as , so that the spectral index, (Fig. 1). In that case, observed brightness temperatures correspond to actual physical temperatures. The data for the  stars reveal a positive temperature gradient, reminiscent of the solar chromosphere, and different frequencies probe the temperature stratification of the atmosphere. To determine the chromospheric height values , requires a structure model of the atmosphere, that details the run of density and fractional ionization of the gas (De la Luz et al., 2014; Loukitcheva et al., 2015, and references therein).

In Fig. 3, for the disk integrated  is compared with observed values for the Quiet Sun (Loukitcheva et al., 2004).

4.2 Theoretical model chromospheres for

The region close to the temperature minimum is optically thick in the FIR/submm (Liseau et al., 2015) which, as a consequence of the negative temperature gradient, limits our view to higher, cooler layers above the optical photosphere. Therefore, the received flux at a given frequency measures directly the temperature of the plasma at a particular atmospheric height. That can be used to construct analytically the temperature profile to first order and over a limited region, e.g. Liseau et al. (2015, and refereces therein).

A more sophisticated method is to build a theoretical model chromosphere that at its base is anchored in the photosphere. The result of this is displayed in Fig. 3, showing both the temperature minimum and the temperature increase that are retrieved by the semi-empirical non-LTE model chromosphere of , based on a modified hydrostatic equilibrium model (C7) of the solar chromosphere (Avrett & Loeser, 2008; De la Luz et al., 2014). C7 can be viewed as an average of the five most widely used solar chromosphere models (Vernazza et al., 1981; Loukitcheva et al., 2004; Fontenla et al., 2007; Avrett & Loeser, 2008; De la Luz et al., 2014).

The temperature profile is computed iteratively from the modified density/pressure structure, ionization balance and opacity (lines and continua). As the conditions in the chromosphere strongly deviate from thermodynamical equilibrium, both the ionization-excitation and the radiative transfer are treated in non-LTE (De la Luz & Tapia, in preparation). Figure 3 also shows the sharp drop in proton density and the increase of the turbulent speed , steepening into shocks. Although  is not a solar analogue like A, a modified solar model also provides an acceptable fit to the data. The modeled of the K-star  is also shown in Fig. 3.

For , the temperature profile is shallower than for the Sun and  K at  km, where the proton density  cm. The corresponding model parameters for  are 3407 K, 560 km and  cm, respectively.

The temperature minimum in the -scale of the  model, , is as low as what has been observed in CO lines from the Sun (, Avrett, 2003, and references therein). This is lower than what traditionally has been derived from the wings of resonance lines, viz. for both  and the Sun (Ayres et al., 1976; Avrett, 2003).

At the longest wavelengths the exponent of the observed SED changes, likely because the free-free emission is turning from optically thick to thin beyond 1.5 mm (frequency exponent tends from about 2 to 0). Especially at 3 mm, the Band 3 data are not well reproduced by the model, the density of which is too low to generate sufficient free-free and H opacity for the required flux. However, from Table 4, it can be seen that the radiation from  in Band 4 and 3 probably originates rather high up, at about 2000 km and near the base of the transition region (TR) into the hot corona, which is seen in the X-rays from the  binary (DeWarf et al., 2010; Ayres, 2014). The X-ray emission is particularly strong from the more active companion .

It is likely that it is in these thin layers of the TR base, where wave energy is dumped and dissipated (Soler et al., 2015; Shelyag et al., 2016). Therefore, this region is critical to the understanding of the heating processes of the outer atmospheres of the stars and the Sun. Given the available evidence, ALMA Band 5 observations will eventually be particularly crucial for the observation of these layers in the  stars. These stars deserve continued monitoring, including observations at longer wavelengths.

 and B are known to be variable on both short and long time scales (DeWarf et al., 2010; Ayres, 2015). In X-rays and the FUV, both stars show flickering but also solar-like magnetic cycles, with  being the more active one. Repeat observations would assess the level of activity in the submm/mm regime. Between 2014 and 2017,  is expected to go through its broad shallow maximum of its  year cycle, whereas B will presumably pass through a minimum of its 8 year cycle. Thus, perhaps in contrast to the solar case, changes of the chromospheric emission from the active K-dwarf could occur over a period of a few years, although such behavior, by analogy with the Sun, would not be expected for the less active .

Figure 3: Top: The SED of the model chromosphere of the G2 star , based on the modified solar C 7 model, is shown by the blue curve. Data are from Spitzer, Herschel and APEX (Liseau et al., 2013) (small blue sysmbols and dotted error bars) and from ALMA (big blue squares). The ALMA bands are indicated at the bottom of the figure and the stellar photosphere is shown as RJ(. For comparison, data for the Quiet Sun from Loukitcheva et al. (2004) are shown as black open circles. Middle: The run of with height with symbols as above. For comparison, also the corresponding model for the K 1 star  is shown in red, and, for reference, the solar C 7 model as black dots. Bottom: The run of density (H) and turbulent velocity with height is shown for the solar analogue .
Figure 4: Brightness temperatures for six solar-type stars at wavelengths from 0.5 mm to 6 cm (see the text). Detections were obtained below 1 cm and merely upper limits above that wavelength. The color coding and stellar identifications are given in the upper left corner of the figure. The open circles denote estimates of future ATCA detections of B in 12 hours at 20 and 6 GHz, respectively (see the text).

4.3 Comparison with other stars

4.3.1 Solar-type

In addition to , a handful of other solar-type stars (late F to early K) have been observed at long wavelengths. These stars are all within 6 pc. For  Eridani (K2 V) measurements have been made at 1.3 mm and 7 mm (MacGregor et al., 2015) and at 3.6 cm (Güdel, 1992) and 6 cm (Bower et al., 2009); for  Ceti (G8.5 V) at 1.3 mm (MacGregor et al., 2016) and at 8.7 mm and 2 cm (Villadsen et al., 2014). Further, 40 Eridani A (K0.5 V) and  Cassiopeiae A (F9 V) at 8.7 mm, and the latter also at 2 and 6 cm, have also been observed by Villadsen et al. (2014).

As seen in Fig. 4, there is only limited overlap with the wavelength domain of the  binary and upper limits, rather than detections, dominate at cm-wavelengths. However, for all detected cases (4 stars in addition to  and B), the fluxes were not consistent with photospheric values but significantly higher. Therefore, it was generally concluded that this excess emission originates in stellar chromospheres, similar to those in the Sun and B.

Figure 5: Left: Band 8 observation of tauri on 2 May 2015, with the color bar for the intensity shown below. Apart from the well known binary  and , a previously unknown source U was discovered less than 6 north of the primary A. The displayed image is primary beam corrected and the slightly oval synthesized beam (Table 1) is shown in grey in the lower left corner. In the FITS image, the FK5 (J2000.0) coordinates at mid-integration, i.e. JD 2457144.632219328, are R.A.=143928491, Dec.= 49 5183. Right: The logarithmic submm/mm-SED of the unidentified source U near  is consistent at the level with that of a blackbody, as indicated by the dashed line of slope 2.0 (cf. Table 5).

4.3.2 Giants

Harper et al. (2013, and references therein) discuss ongoing observational and theoretical work on giants (luminosity class III), addressing the possibility to observationally separate acoustic from MHD heating processes in the upper atmospheres due to the large scale heights in these stars. Their convective cells and envelopes are much larger than those of main-sequence stars, which may make possible to observationally distinguish between these effects. In addition and in contrast to the smaller and more compact main-sequence stars (class V), giants are relatively bright and hence offer themselves as possible candidates for calibration purposes for observations in the submm/mm/cm regime (see also Cohen et al., 2005).

4.4 A new, unidentified point-like source near

In May 2015, an unidentified object was detected in the Band 8 observations of  (Fig. 5). This point source, designated U and with integrated flux over the band of about 4 mJy (Table 5), was within a few arcseconds of the binary. As this object was not detected in any other data set (including UV, VIS and NIR with HST and VLT, see Kervella et al., 2006; Kervella & Thévenin, 2007), other epoch data are lacking and hence its nature is unknown.

Figure 5 also displays the SED of this object, consisting of one detection and five upper limits at the level. However, the data could be consistent with blackbody emission, viz. , and may be due to a submm galaxy, a stellar object, a brown dwarf or a planetary object. A companion star of the  system does not present a viable explanation, as any star would be brighter than 10th magnitude in the V-band, and hence must be discarded.

The submm galaxy option would imply that the proper motion of U would be minuscule, and that it would be quickly left behind the  stars as they pace, at the rate of 37 yr, through the sky. As  is in close projection to the plane of the Galaxy, a stellar nature of U may perhaps appear more natural. However, this putative star remained undetected in recent deep searches, implying that U is either a distant heavily extinguished background star or a nearby, very cold object, i.e. a brown dwarf or a planetary object. The parallax and proper motion would clearly distinguish among these possibilities.

Verly low-temperature brown dwarves like the T 8.5-type ULAS  with an estimated temperature of 575 K, or the even cooler Y 2 object WISE J with  K (Tinney et al., 2014; Leggett et al., 2015), can serve as known examples, i.e. an extremely cool brown dwarf at a distance of nearly 20 000 AU may be a viable candidate for the identification of source U. However, like the Y2 object, the Wide-field Infrared Survey Explorer (WISE) should have picked it up. Unless close to the very bright B, the moderate angular resolution of WISE (60) presented an obstacle to a clean detection.

In the solar system, the projected offset of  55 would correspond to a distance between Jupiter and Saturn666The accuracy of the absolute stellar positions will be addressed by Kervella et al. 2016 (in preparation). However, the identification of U as a planetary companion of  would be totally unrealistic, because the observed 740 m-flux would be too high by several orders of magnitude. If a body of planetary dimensions, U would possibly be bound to the solar system, but its distance would presently be undetermined. Fig. 6 shows the distances and flux densities at 740 m estimated for several known dwarf planets with the diameter as parameter. From the figure it is evident that U is likely more distant than Pluto, since an  1000 km body at roughly 40 AU would have been known for a long time, i.e. for at least ten years. For example, when examining a total of 766 925 known solar-system objects777 for being within 15 around  at the time of observation, we found no source down to the limiting -magnitude of 26.0. Therefore, a low-albedo, thermal Extreme Trans Neptunian Object (ETNO), would clearly be consistent with our data (see Fig. 6).

Band 9 Band 8 Band 7 Band 6 Band 4 Band 3
679 GHz 405 GHz 343.5 GHz 233 GHz 145 GHz 97.5 GHz
442 m 740 m 873 m 1287 m 2068 m 3075 m
Table 5: Primary beam corrected flux density and upper limits for the U-source in mJy
Figure 6: Band 8 flux density as function of the distance from the Sun with diameters as parameter, in 10 km and next to or atop the curves and arbitrarily limited to 6000 km, i.e. slightly smaller than the diameter of Mars. Both the surface temperature and the radius are a priori undetermined. A few known TNOs with their names are shown by the red dots ( In parentheses, the apparent diameter in milli-arcseconds and the estimated blackbody temperature are given. The size of the ALMA synthesized beam,  mas, is given in the upper right corner, confirming that these objects would be point-like to ALMA. The observed Band 8 flux density of the unidentified object is indicated by the horizontal blue-shaded dashed line (). The distance to the U-source remains to be determined.

5 Conclusions

Below, we briefly summarize our main conclusions.

  • ALMA observations of tauri at 0.44, 0.74, 0.87, 1.3, 2.1 and 3.1 mm clearly resolved the binary, but not the stellar disks, at all wavelengths. The spectral energy distributions of these continuum measurements are consistent with radiation that follows , except at the lowest frequencies where the SEDs appear to flatten. This is particularly pronounced for the more active secondary, a K 1 star, possibly indicative of time variability within half a year or, perhaps more likely, of optically thin free-free emission.

  • The ALMA data have been modeled with modified solar chromosphere models which result in the physical structure of the stellar chromospheres. This adapted solar model works very well for the solar analog  (G2 V), but also for the K1 V star . Comparison with the data indicates that the temperature minima of both  and B are deeper than on the Quiet Sun. These correspond to the low temperatures seen in lines of the CO molecule on the Sun and occur at atmospheric heights of 615 km and 560 km, respectively.

  • The ALMA data for B can be put into context with observations of other nearby solar-type stars that show that chromospheric mm-wave emission is a common feature among these stars and that an increase in the sample size can be expected in the near future.

  • The ALMA imaging at 0.74 mm led to the discovery of a previously unknown point source within a projected distance of 7.5 AU from B. The ALMA observations were performed at different occasions during one year (2014 - 2015), but this source was clearly detected only on one date. At the three sigma level, the SED of this object is consistent with that of a blackbody and we speculate about its nature. Unless it is a highly variable background source, we find it most likely that it is a distant member of our solar system.

We thank the referee for his/her valuable comments on the manuscript. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00170.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ.


  • Avrett (2003) Avrett, E. H. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 286, Current Theoretical Models and Future High Resolution Solar Observations: Preparing for ATST, ed. A. A. Pevtsov & H. Uitenbroek, 419
  • Avrett & Loeser (2008) Avrett, E. H. & Loeser, R. 2008, ApJS, 175, 229
  • Ayres (2014) Ayres, T. R. 2014, AJ, 147, 59
  • Ayres (2015) Ayres, T. R. 2015, AJ, 149, 58
  • Ayres et al. (1976) Ayres, T. R., Linsky, J. L., Rodgers, A. W., & Kurucz, R. L. 1976, ApJ, 210, 199
  • Bower et al. (2009) Bower, G. C., Bolatto, A., Ford, E. B., & Kalas, P. 2009, ApJ, 701, 1922
  • Brott & Hauschildt (2005) Brott, I. & Hauschildt, P. H. 2005, in ESA Special Publication, Vol. 576, The Three-Dimensional Universe with Gaia, ed. C. Turon, K. S. O’Flaherty, & M. A. C. Perryman, 565
  • Cohen et al. (2005) Cohen, M., Carbon, D. F., Welch, W. J., et al. 2005, AJ, 129, 2836
  • De la Luz et al. (2014) De la Luz, V., Chavez, M., Bertone, E., & Gimenez de Castro, G. 2014, Sol. Phys., 289, 2879
  • Demory et al. (2015) Demory, B.-O., Ehrenreich, D., Queloz, D., et al. 2015, MNRAS, 450, 2043
  • DeWarf et al. (2010) DeWarf, L. E., Datin, K. M., & Guinan, E. F. 2010, ApJ, 722, 343
  • Dulk (1985) Dulk, G. A. 1985, ARA&A, 23, 169
  • Dumusque et al. (2012) Dumusque, X., Pepe, F., Lovis, C., et al. 2012, Nature, 491, 207
  • Eiroa et al. (2013) Eiroa, C., Marshall, J. P., Mora, A., et al. 2013, A&A, 555, A11
  • Fontenla et al. (2007) Fontenla, J. M., Balasubramaniam, K. S., & Harder, J. 2007, ApJ, 667, 1243
  • Güdel (1992) Güdel, M. 1992, A&A, 264, L31
  • Harper et al. (2013) Harper, G. M., O’Riain, N., & Ayres, T. R. 2013, MNRAS, 428, 2064
  • Hatzes (2013) Hatzes, A. P. 2013, ApJ, 770, 133
  • Kervella & Thévenin (2007) Kervella, P. & Thévenin, F. 2007, A&A, 464, 373
  • Kervella et al. (2006) Kervella, P., Thévenin, F., Coudé du Foresto, V., & Mignard, F. 2006, A&A, 459, 669
  • Kervella et al. (2003) Kervella, P., Thévenin, F., Ségransan, D., et al. 2003, A&A, 404, 1087
  • Leggett et al. (2015) Leggett, S. K., Morley, C. V., Marley, M. S., & Saumon, D. 2015, ApJ, 799, 37
  • Liseau et al. (2013) Liseau, R., Montesinos, B., Olofsson, G., et al. 2013, A&A, 549, L7
  • Liseau et al. (2015) Liseau, R., Vlemmings, W., Bayo, A., et al. 2015, A&A, 573, L4
  • Loukitcheva et al. (2004) Loukitcheva, M., Solanki, S. K., Carlsson, M., & Stein, R. F. 2004, A&A, 419, 747
  • Loukitcheva et al. (2015) Loukitcheva, M., Solanki, S. K., Carlsson, M., & White, S. M. 2015, A&A, 575, A15
  • MacGregor et al. (2016) MacGregor, M. A., Lawler, S. M., Wilner, D. J., et al. 2016, ArXiv e-prints
  • MacGregor et al. (2015) MacGregor, M. A., Wilner, D. J., Andrews, S. M., Lestrade, J.-F., & Maddison, S. 2015, ApJ, 809, 47
  • Martí-Vidal et al. (2014) Martí-Vidal, I., Vlemmings, W. H. T., Muller, S., & Casey, S. 2014, A&A, 563, A136
  • Press et al. (1986) Press, W. H., Flannery, B. P., & Teukolsky, S. A. 1986, Numerical recipes. The art of scientific computing
  • Rajpaul et al. (2016) Rajpaul, V., Aigrain, S., & Roberts, S. 2016, MNRAS, 456, L6
  • Shelyag et al. (2016) Shelyag, S., Khomenko, E., de Vicente, A., & Przybylski, D. 2016, ApJ, 819, L11
  • Soler et al. (2015) Soler, R., Carbonell, M., & Ballester, J. L. 2015, ApJ, 810, 146
  • Thévenin et al. (2002) Thévenin, F., Provost, J., Morel, P., et al. 2002, A&A, 392, L9
  • Tinney et al. (2014) Tinney, C. G., Faherty, J. K., Kirkpatrick, J. D., et al. 2014, ApJ, 796, 39
  • Torres et al. (2010) Torres, G., Andersen, J., & Giménez, A. 2010, A&A Rev., 18, 67
  • Trigilio et al. (2013) Trigilio, C., Buemi, C., Umana, G., & Leto, P. 2013, Is alpha Centauri revealed by the new generation of radio telescopes?, ATNF Proposal
  • Vernazza et al. (1981) Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635
  • Villadsen et al. (2014) Villadsen, J., Hallinan, G., Bourke, S., Güdel, M., & Rupen, M. 2014, ApJ, 788, 112
  • Wang & Fischer (2015) Wang, J. & Fischer, D. A. 2015, AJ, 149, 14
  • Wedemeyer et al. (2015) Wedemeyer, S., Bastian, T., Brajsa, R., et al. 2015, ArXiv e-prints
  • Wedemeyer-Böhm et al. (2007) Wedemeyer-Böhm, S., Steiner, O., Bruls, J., & Rammacher, W. 2007, in Astronomical Society of the Pacific Conference Series, Vol. 368, The Physics of Chromospheric Plasmas, ed. P. Heinzel, I. Dorotovič, & R. J. Rutten, Heinzel

Appendix A Sub-band fluxes

Figure 7: Measurements of the flux density of  (blue circles) and of  (red circles) in the sub-band windows (spw), see Table 6. The error bars represent the rms-values. Band 9 is too narrow to allow meaningful measurement in sub-windows and only a single value is given.

The flux densities of the spectral windows per band are provided in Table 6 () and the data are plotted in Fig. 7. For Band 9, only a single value is given, as the windows are too narrow for meaningful individual measurement.

For , the relative drop in intensity in the second spw of Band 4 is conspicuous. This is not evident for , and the glitch can therefore not be caused by different calibrations. B are point sources and were observed simultaneously. Hence, simultaneous visibility fitting, with fixing the positions to reduce the noise (Martí-Vidal et al. 2014), should not result in such large differences, unless there is something in the data, e.g. a spectral feature, in  that is not present in the SED of . New observations of B, at higher S/N in Band 4, would be necessary to resolve this issue.

B (A) rms(A) (B) rms(B)
(GHz) (mJy) (mJy) (mJy) (mJy)
3 190.49459 3.03 0.04 1.42 0.06
3 192.43209 3.00 0.11 1.41 0.12
3 102.4946 3.75 0.06 1.64 0.07
3 104.4946 3.81 0.06 1.83 0.05
4 138.7133 5.92 0.15 2.50 0.15
4 140.6508 5.96 0.14 2.32 0.14
4 149.2758 6.74 0.15 2.62 0.15
4 151.2758 6.82 0.16 2.98 0.16
6 224.000 13.12 0.22 5.37 0.14
6 226.000 13.75 0.17 6.21 0.09
6 240.000 14.33 0.14 7.03 0.18
6 242.000 14.64 0.46 6.43 0.11
7 336.4946 26.75 0.53 10.50 0.59
7 338.4321 25.18 0.39 10.81 0.58
7 348.4946 25.69 0.38 12.22 0.27
7 350.4946 26.25 0.48 12.61 0.90
8 397.9946 37.17 0.57 15.66 0.28
8 399.9321 35.47 0.65 15.69 0.51
8 409.9946 38.65 0.65 17.39 0.57
8 411.9946 38.19 0.38 17.84 0.51
9 678.9600 107.20 1.50 57.60 4.50
Table 6: Sub-band (spw) flux densities for B
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description