c ESO 2014
Astronomy & Astrophysics manuscript no. harding_AA_2012_20865
November 7, 2014
Spin-orbit alignment in the very low mass binary regime
The L dwarf tight binary 2MASSW J0746425+200032AB
L. K. Harding1, 2 , G. Hallinan1 , Q. M. Konopacky3,⋆ , K. M. Kratter4,⋆⋆ , R. P. Boyle5 , R. F. Butler2 , and A. Golden6, 2
arXiv:1304.5290v1 [astro-ph.SR] 19 Apr 2013
1
2
3
4
5
6
Cahill Center for Astrophysics, California Institute of Technology, 1200 E. California Blvd., MC 249-17, Pasadena, CA 91125,
USA
e-mail: lkh@astro.caltech.edu; gh@astro.caltech.edu
Centre for Astronomy, National University of Ireland, Galway, University Road, Galway, Ireland
e-mail: ray.butler@nuigalway.ie
Dunlap Institute for Astronomy and Astrophysics, University of Toronto, 50 St. George Street, Toronto M5S 3H4, Ontario, Canada
e-mail: konopacky@di.utoronto.ca
JILA, University of Colorado, Boulder, CO 80309-0440, USA
e-mail: kaitlin.kratter@jila.colorado.edu
Vatican Observatory Research Group, Steward Observatory, University of Arizona, Tucson, AZ 85721, USA
e-mail: boyle@ricci.as.arizona.edu
Department of Genetics (Computational Genetics), Albert Einstein College of Medicine, Bronx NY 10461, USA
e-mail: aaron.golden@einstein.yu.edu
Accepted April 17, 2013
ABSTRACT
Studies of solar-type binaries have found coplanarity between the equatorial and orbital planes of systems with <40 AU separation.
By comparison, the alignment of the equatorial and orbital axes in the substellar regime, and the associated implications for formation
theory, are relatively poorly constrained. Here we present the discovery of the rotation period of 3.32 ± 0.15 hours from 2MASS
J0746+20A - the primary component of a tight (2.7 AU) ultracool dwarf binary system (L0+L1.5). The newly discovered period,
together with the established period via radio observations of the other component, and the well constrained orbital parameters and
rotational velocity measurements, allow us to infer alignment of the equatorial planes of both components with the orbital plane of the
system to within 10 degrees. This result suggests that solar-type binary formation mechanisms may extend down into the brown dwarf
mass range, and we consider a number of formation theories that may be applicable in this case. This is the first such observational
result in the very low mass binary regime. In addition, the detected period of 3.32 ± 0.15 hours implies that the reported radio period
of 2.07 ± 0.002 hours is associated with the secondary star, not the primary, as was previously claimed. This in turn refutes the
claimed radius of 0.78 ± 0.1 R J for 2MASS J0746+20A, which we demonstrate to be 0.99 ± 0.03 R J .
Key words. binaries: close – brown dwarfs – stars: formation – stars: low-mass – stars: rotation
1. Introduction
Investigation of the spin and orbital properties of solar-type binaries has been underway for several decades (Weis 1974; Abt
& Levy 1976; Bodenheimer 1978; Fekel 1981; Hale 1994, and
references therein). A large fraction of close (<40 AU) solartype binaries have been found to have coplanar spins and orbits (Hale 1994). Studies by Weis (1974) yielded similar results over a wider range of spectral classes [B - F]. More recent
studies have turned up a diversity of configurations. Jensen et
al. (2004) and Monin et al. (2006) have reported planar alignment and misalignment for wider-separation binaries. Similar
to these studies, in the case of very close binaries of semi-major
axis ∼0.3 AU, there have been some examples where systems
exhibit both aligned and misaligned axes (Albrecht et al. 2009,
2011). More recently, Wheelwright et al. (2011) have also reported coplanarity between HAe/Be binary systems and circumstellar disks. Thus, it appears that while the spin-orbit alignment
is common, there are exceptions at all orbital separations.
⋆
⋆⋆
Dunlap Fellow
Hubble Fellow
In this paper, we turn our attention to the spin-orbit alignment of stars at the bottom end of the initial mass function, where
binarity is less common (Lada 2006; Raghavan et al. 2010). In
the years following the detection of the first brown dwarf Gl
229B by Nakajima et al. (1995), a number of surveys yielded
the discovery of a large number of low mass star (< 0.1 M⊙ ) and
brown dwarf binary systems, e.g. Burgasser et al. (2007, and references therein). Following these discoveries, the introduction of
laser guide star (LGS) adaptive optics (AO) systems on groundbased telescopes provided the means of assessing the dynamical mass of such systems (Bouy et al. 2004; Dupuy et al. 2010;
Konopacky et al. 2010). More recently, Konopacky et al. (2012)
obtained resolved LGS AO spectroscopic measurements of individual component rotation velocities for a sample of eleven very
low mass (VLM) dwarf binaries (Mtot ≤ 0.185 M⊙ , as defined by
Close et al. (2003)). These data provided additional parameters
for sources separated by ∼1 - 10 AU, but could only be used to
tentatively investigate spin-orbit alignment, since other parameters such as individual component rotation periods, and system
properties inferred from evolutionary models, such as radii, were
still either unknown or poorly constrained.
Article number, page 1 of 9
Table 1. Properties of 2MASSW J0746425+200032AB
Parameter
Rotation period (hrs)
v sin i (km s−1 )
Equatorial Velocity (km s−1 )
Period ratio
v sin i ratio
Orbital period (yrs)
Semi-major axis (mas)
Inc.ORB (deg)‡
Inc.EQ (deg)‡
Age (log yrs)
Masstotal (M⊙ )
Mass (M⊙ )
Lithium?
Radius (R J )
Gravity (log g)
Lbol (log L/L⊙ )
Abs. mag (J)
Abs. mag (H)
Abs. mag (K)
References
2MASS J0746+20A
3.32 ± 0.15†
19 ± 2
36 ± 4†
0.62 ±0.02
0.03
0.57±0.13
0.10
12.71 ± 0.07
237.3+1.5
−0.4
41.8 ± 0.5
32 ± 4†
9.1 ± 0.1†
0.151 ± 0.003
0.078 ± 0.004†
No
0.99 ± 0.03†
5.34 ± 0.02†
-3.64 ± 0.02
11.85 ± 0.04
11.13 ± 0.02
10.62 ± 0.02
1, 3-5
2MASS J0746+20B
2.07 ± 0.002
33 ± 3
56 ± 2†
...
...
12.71 ± 0.07
237.3+1.5
−0.4
41.8 ± 0.5
36 ± 4†
9.1 ± 0.1†
0.151 ± 0.003
0.073 ± 0.004†
No
0.96 ± 0.02†
5.34 ± 0.02†
-3.77 ± 0.02
12.36 ± 0.10
11.54 ± 0.03
10.98 ± 0.02
1-5
Notes. † Rotation period discovered in this work; all other parameters derived here based on evolutionary models of Chabrier et al. (2000).
‡Inc.ORB is the system orbital inclination as measured by Konopacky et al. (2012), whereas Inc.EQ is the equatorial inclination of each component,
calculated in this work.
References. (1) This work; (2) Berger et al. (2009); (3) Konopacky et al. (2010); (4) Bouy et al. (2004); (5) Chabrier et al. (2000).
Previous studies have shown that the presence of magnetic fields can affect binary formation (Mestel 1977; Bodenheimer 1978; Fekel 1981; Li et al. 2004, and references therein),
whereby such fields could potentially contribute to a loss of angular momentum on large scales, or the tilting of stellar spins on
small scales. Are these concerns of special relevance in the ultracool dwarf regime? M dwarfs later than M3 are now associated
with intense magnetic activity, often possessing surface magnetic field strengths of a few kG and greater (Reiners & Basri
2007; Hallinan et al. 2008; Berger et al. 2009; West et al. 2011).
In this paper, we now have sufficient data to investigate the orbital properties for the magnetically-active VLM dwarf binary,
2MASS J0746+20AB (Mtot = 0.151 ± 0.003 M⊙ ). This detection can shed light on the formation of VLM binary stars, and
could signal that a scaled-down version of the formation mechanism for solar-type binary systems holds in this regime, despite
the presence of a ∼1.7 kG magnetic field in this case (Antonova
et al. 2008). Characterizing the fundamental properties of VLM
binary star formation is important in establishing a correlation,
if any, in the formation and evolution of all types of binary stars.
2. 2MASSW J0746425+200032AB
2.1. Properties
2MASS J0746+20AB is an L dwarf binary (L0+L1.5) that is
located at a distance of 12.20 ± 0.05 pc (Dahn et al. 2002). It
is a tight binary system, with a separation of ∼2.7 AU (Reid
et al. 2001) and an effective temperature of between 1900 2225 K (Vrba et al. 2004). The latest high-precision dynamical
mass measurements yielded a total mass of 0.151 ± 0.003 M⊙
(Konopacky et al. 2010), initially measured to be 0.146 ±0.016
0.006
M⊙ by Bouy et al. (2004), placing the dwarf in the VLM binary regime. Bouy et al. (2004) based their individual compoArticle number, page 2 of 9
nent masses on a model-derived age estimate. However, Gizis
& Reid (2006) questioned these individual mass estimates, and
argued against the secondary component being a brown dwarf as
predicted by Bouy et al. (2004), and instead favored a substellar
or low mass star classification.
The L dwarf was reported as an active radio source by
Antonova et al. (2008), who estimated magnetic field strengths
of ∼1.7 kG based on the detection of a single highly polarized
burst of emission. The dominant emission was quiescent, the nature of which is still debated. Following this observation, Berger
et al. (2009) reported radio emission with a rotation period of
2.07 ± 0.002 hours. Simultaneously, they detected periodic Hα
emission. The period was the same in both instances, and was
consistent with stellar rotation. They also estimated a stellar radius of 0.078 ± 0.010 R J for 2MASS J0746+20A; in order to
do this, they assigned an average v sin i of 27 ± 3 km s−1 (taken
from a number of previous studies (Reid et al. 2002; Bailer-Jones
2004; Reiners & Basri 2008)) to the primary. This assumption
was based on the observations of Bouy et al. (2004), who established that 70% of the light contribution was coming from the
primary star. Thus, they argued that both the Hα and radio emission were also emanating from the primary component of the
system. In addition, marginal evidence of periodicity (of a few
hours) has previously been reported by Clarke et al. (2002) and
Bailer-Jones (2004). Target properties are shown in Table 1.
2.2. Instrumentation and observations
We obtained observations over a ∼2 year baseline for 2MASS
J0746+20 to search for periodic photometric variability from either, or both components of the binary system. The observations encompassed 4 separate epochs in January 2009, February
2010, November 2010 and December 2010. These were carried
L. K. Harding et al.: Spin-orbit alignment in the very low mass binary regime
Table 2. Observation details
Source
Date
of obs.
(UT)
Length
of obs.
(∼hrs)
Exp.
time
(s × coadd)
Photometric
error
(%)
Band
Telescope
/ instrument
(1)
2MASS J0746+20AB
(2)
2009 Jan 25
2009 Jan 26
2009 Jan 28
2010 Feb 19
2010 Feb 20
2010 Nov 13
2010 Nov 14
2010 Dec 2
2010 Dec 12
2010 Dec 13
2010 Dec 14
(3)
6.0
6.8
7.4
4.5
4.0
4.6
5.5
6.0
3.0
6.8
7.0
(4)
25 × 1
25 × 1
25 × 1
5 × 12
5 × 12
5 × 12
5 × 12
5 × 12
5 × 12
5 × 12
5 × 12
(5)
0.21
0.28
0.24
0.27
0.30
0.31
0.33
0.25
0.29
0.32
0.34
(6)
I
I
I
I
I
I
I
I
I
I
I
(7)
VATT/4K
VATT/4K
VATT/4K
VATT/GUFI
VATT/GUFI
VATT/GUFI
VATT/GUFI
VATT/GUFI
VATT/GUFI
VATT/GUFI
VATT/GUFI
Notes. Column (1) Target of campaign. (2) Date of each observation in UT. (3) The total time observed per night for each observation, in hours.
(4) The exposure time for each observation in seconds, followed by the binning factor used to increase data point S/N. (5) The mean photometric
error for light curve data points per night. (6) The VATT Johnson photometric waveband used (I: ∼7200 - 9100 Å). (7) Telescope and detector
used for a given observation. These observations were taken over 4 separate epochs, spanning ∼62 hours of data over a 2 year baseline.
out with the GUFI (Galway Ultra-Fast Imager) and VATT 4K
photometers on the VATT telescope7 using the VATT I-Arizona
filter (∼7200 - 9100 Å). GUFI was first commissioned by astronomers at NUI Galway to operate on the 1.5 m Loiano telescope in Bologna, Italy (Sheehan & Butler 2008). It was modified thereafter for use on the VATT telescope, where it is currently stationed as a visitor instrument. GUFI was developed as
a dedicated high-speed photometer for a long-term observational
campaign that focuses on the periodic variability of a sample of
radio detected ultracool dwarfs, and the associated mechanisms
responsible for such periodic signals (Harding et al. 2013). It
provides a FOV of ∼ 3′ × 3′ , a plate scale of 0.35′′ pixel−1 , and
performs high-time resolution imaging (e.g. 34 frames per second full frame, with ∼2 ms readout rates). Data was taken with
exposure times of 5 seconds where frames were later summed
in image space to 1 minute to increase the signal to noise (S/N).
GUFI was used for the February 2010, November 2010 and December 2010 epochs.
We used the VATT 4K photometer for the first epoch of the
2MASS J0746+20 observations in January, 2009, before GUFI’s
commissioning in May 2009. The VATT 4K is the facility photometer of the VATT observatory, and has a native plate scale of
0.188′′ pixel−1 and a FOV of ∼12.5′ × 12.5′ . Details of observations, and the typical photometric error per given observation
are shown in Table 2, where the full sample includes ∼62 hours
of observations.
pending). We also carried out photometry for all reference stars
in order to measure their level of variability - this comparison
ensured that the periodicity was indeed intrinsic to the binary.
Reference stars were chosen after assessing their stability (in
terms of flux) throughout the night, their isolation on the CCD
chip, the properties of their seeing profiles, and their color index with respect to the target. Photometric apertures (in pixels)
which provided the highest S/N for the target star were selected
for aperture photometry. Since GUFI did not resolve each component of the binary system, the chosen aperture included the
combined flux from both stars. Differential photometry was obtained by dividing the target flux by the mean flux of selected
reference stars. Finally, after differential photometry was carried
out, we assessed the photometric light curves for evidence of
any residual systematic trends due to e.g. the effects of increasing/decreasing airmass. These effects were negligible throughout the campaign, where typical seeing at the VATT observatory
was ∼0.7 - 1.6′′ . However, where small trends were observed
due to such second order extinction effects, we fitted a linear (or
quadratic, trend depending) fit to the data and then applied an
inverse form of the chosen fit. Removing trends is important in
a variable light curve, since they can introduce aliasing effects
in the subsequent spectral analysis (discussed in Section 3). The
calculated photometric errors are shown in Table 2 (column 5),
and the procedure for calculating these errors is outlined in Section 3.
2.3. Data reduction and differential photometry
3. Periodic Variability
Standard data reduction techniques were employed via a GUFI
pipeline (Sheehan & Butler 2008), which performs bias subtraction (using zero-integration frames) and flat-fielding (using twilight flat-fields). The FOVs of GUFI and the VATT 4K provided
up to 20 reference stars for the target field (detector FOV de-
3.1. Variability analysis
7
The Vatican Advanced Technology Telescope (VATT) telescope (1.83 m) facility is operated by the Vatican Observatory, and is part of the Mount Graham International Observatory.
Details of VATT detectors can be found here:
http://vaticanobservatory.org/VATT/index.php/telescope-instruments/.
The Lomb-Scargle (LS) periodogram (Lomb 1976; Scargle
1982) was used to search for periodic variability in the binary
light curves. This technique is very effective for unevenly spaced
data, and uses the discrete Fourier Transform (DFT). The algorithms output power spectra with a range of peaks of different
amplitudes. The significance of these peaks are then analyzed to
a given sigma level as computed by the LS algorithm, which corresponds to possible periodic variability. We selected a range of
peaks of >5σ significance, and then inspected the various soluArticle number, page 3 of 9
HJD: 2454857.631957755
1.01
1.01
1.01
1
0.99
6
1
0.99
0.98
8
10
January 25 2009 UT (hrs)
Relative flux
1.02
0.98
4
6
8
January 26 2009 UT (hrs)
0.99
0.98
10
1.01
1.01
0.99
6
7
8
February 19 2010 UT (hrs)
1
0.99
0.98
9
Relative flux
1.01
1
5
HJD: 2455514.814986587
6
7
February 20 2010 UT (hrs)
1
0.99
0.98
8
1.01
1.01
0.99
8
9
10
11
12
November 14 2010 UT (hrs)
0.99
0.98
13
Relative flux
1.01
Relative flux
1.02
1
HJD: 2455543.815873595
8
1
0.99
0.98
10
12
December 2 2010 UT (hrs)
1.01
1.01
0.99
9
10
11
12
December 13 2010 UT (hrs)
13
Relative flux
1.01
Relative flux
1.02
8
11
12
December 12 2010 UT (hrs)
1
0.99
0.98
8
10
12
December 14 2010 UT (hrs)
13
Reference star lightcurve
1.02
0.98
10
HJD: 2455544.785973326
1.02
1
9
10
11
12
November 13 2010 UT (hrs)
HJD: 2455542.919808880
1.02
0.98
8
HJD: 2455532.794271032
1.02
1
6
8
10
January 28 2009 UT (hrs)
HJD: 2455513.820684682
1.02
5
4
HJD: 2455247.685167070
1.02
0.98
Relative flux
1
1.02
Relative flux
Relative flux
HJD: 2455246.705281595
Relative flux
HJD: 2454859.653165393
1.02
Relative flux
Relative flux
HJD: 2454856.710011647
1.02
1
0.99
0.98
8
9
10
11
UT (hrs)
12
Fig. 1. I-band light curves of 2MASS J0746+20AB. Each figure shows relative flux vs. UT times (and dates). We also mark the HJD time above
each figure corresponding to the start-point of each observation. All of the data presented here was taken in the ‘VATT Arizona’ I-band broadband
filter (∼7200 - 9100 Å), where the baseline extends over ∼2 years - as shown in Table 1. In this work, we report the discovery of consistent periodic
variability from 2MASS J0746+20A, with a period of 3.32 ± 0.15 hours, and an amplitude variability range of 0.40 - 1.52%. We note that the
full in-depth analysis investigating the stability of the target’s phase and amplitude is discussed in other work (Harding et al. 2013). Some poor
weather conditions during constituent epochs were present - e.g. January 25 & 26 2009 were taken under thin cloud and high wind conditions
and thus were binned by a factor of 2 compared to other data sets (Table 1). The vertical arrow marked on the November 14 2010 light curve
denotes a period of complete cloud cover, and so data was removed accordingly. The red light curve (bottom right) is an example of a reference
star plotted against all others (taken from November 13 2010), illustrating the stability of the chosen reference stars with respect to the target. The
mean amplitude variability of the reference stars is shown in Table 3.
tions by phase folding the light curves to these periods. Comparing these phase folded light curves from both individual nights,
and from different epochs, allowed us to assess the level of agreement of phase, for a given period.
The standard deviation of the phase fold was also considered.
We calculated this via Phase Dispersion Minimization (PDM)
techniques (Stellingwerf 1978). PDM is a least squares fit calculation where the correct period used in the routines produces
the least data point scatter in the resulting phase fold. This calArticle number, page 4 of 9
culation is called the ‘PDM Theta statistic’ (Θ), where the minimum value is the most likely period solution. We selected a
broad range of periods (e.g. 1 - 5 hours), and ran the PDM algorithm giving a corresponding range of Θ minima. The significance of each Θ was assessed via 105 Monte Carlo trials, which
tests whether any detected Θ minimum could be a result of noise
alone.
Finally, we assessed the peak to peak (PtP) amplitude variability of the target light curves (PtPtarget ) by fitting a model sinu-
L. K. Harding et al.: Spin-orbit alignment in the very low mass binary regime
Table 3. Periodic photometric variability from 2MASS J0746+20A
Phase fold of all epochs
1.02
2M J0746+20A
3.32 ± 0.15
3.318
>>5
3.32
0.40 - 1.52
∼0.36
Photometric
1
2M J0746+20B
2.07 ± 0.002
...
...
...
...
...
Radio/Hα
2
Raw
Bin factor 5
1.015
1.01
1.005
Relative flux
Parameter
Rotation period (hrs)
LS† period (hrs)
LS sign.† (σ)
PDM‡ period (hrs)
PtPtarget (%)
Mean σre f erence (%)
Variability
References
1
0.995
0.99
†
Notes. Period obtained from the Lomb-Scargle (LS) periodogram analysis, as shown in Figure 3 (left). We also include the significance/false
alarm probability of the detection as calculated by the LS algorithm.
‡
Period obtained from the Phase Dispersion Minimization (PDM) technique, shown in Figure 3 (right).
References. (1) This work; (2) Berger et al. (2009).
soidal signal to the raw data using the detected period, and varied
the phase and amplitude of the model. A Chi squared (χ2 ) minimization was performed giving the best fit for a given night’s
peak to peak amplitude variability. The scatter in the reference
star light curves (σre f erence ) were established by calculating their
standard deviation.
3.2. Photometric error and period uncertainty estimation
We estimated the error in the relative magnitude of the target star
using the iraf.phot8 task. The error in magnitude was then converted to an error in flux (since we plot relative flux vs. UT in
our photometric light curves). By using iraf.phot, we take both
formal (e.g. flat-fielding) and informal (e.g. fringing9) errors in
to account - these are usually quite difficult to assess independently.
The period uncertainty was calculated for individual nights
using a χ2 test. We could not run this test on the entire baseline
at once due to large gaps in the time series. Furthermore, we
did not achieve an accurate enough period in consecutive epochs
to phase connect the 2 year baseline together. The quoted period and error in this work was therefore derived within epochs,
which was found to be consistent for all.
3.3. Results
Although we do not resolve each component of the binary as a
point source, most intriguingly, we show optical periodic modulation of 3.32 ± 0.15 hours in VATT I-band (Figure 1). We show
phase folded light curves for this period in Figure 2. Therefore,
this optical periodic variability originates from the other component to that producing the radio emission - reported by Berger
et al. (2009) where the binary exhibited periodic bursts of radio emission of 2.07 ± 0.002 hours. By adopting the estimated
radii in this work (discussed in the next section), in addition to
the v sin i measurements shown in Table 1 (Konopacky et al.
2012), we derive maximum period values of ∼4.22 hours and
8
Image Reduction and Analysis. Found here: http://iraf.noao.edu/.
Fringing is an additive optical effect at red/NIR wavelengths that occurs in the thinned substrate of back-illuminated CCDs. This is due
to atmospheric spectral emission such as OH, and varies in amplitude
across the frame. Thus it needs to be removed if varying at a level
greater than the amplitude variability of the target star. In the case of
2MASS J0746+20, we assessed this effect and found it to be negligible.
9
0.985
0.98
0
0.2
0.4
0.6
0.8
1
Phase
1.2
1.4
1.6
1.8
2
Fig. 2.
Phase folded light curves of all epochs from the 2MASS
J0746+20A observation campaign. The black light curve is folded using the raw light curves as shown in Figure 1, and the red is the same
binned to a factor of 5. Each light curve is phase folded to the detected
period of 3.32 ± 0.15 hours.
∼2.38 hours for 2MASS J0746+20A and 2MASS J0746+20B,
respectively. This infers that we detect the optical rotation period of 2MASS J0746+20A, the primary component of the system. This result is contrary to what was claimed by Berger et
al. (2009), who attributed the radio period to the primary based
upon an incorrect assumption that the v sin i broadening was due
to the primary star, since it is more luminous (Bouy et al. 2004).
A LS periodogram showing the significance of the detection
in our work is shown in Figure 3 (left), and the PDM analysis in Figure 3 (right). Both analyses support a period of 3.32
hours for one binary component. The PtPtarget variability range
for all light curves was found to be ∼0.40 - 1.52%, where we
also observe changes in light curve morphology throughout the
campaign. The source of this behavior is discussed in other work
(Harding et al. 2013). However, in this paper we do investigate
the presence of the other binary component - since variable chromospheric emission and an active magnetosphere was reported
(Berger et al. 2009). Thus the presence of magnetic surface features such as hot or cool spots must be considered for both components. We highlight a cluster of reasonably low false alarm
probability power spectra around ∼12 days−1 (or ∼2 hours), approximating the radio period. To further assess the possible presence of an underlying period from the secondary component, we
fitted a sinusoidal signal with a period of 3.32 hours to each of
the raw light curves. We then subtracted this out, and re-ran
LS periodograms on each of the remaining data points. For each
night, those probabilities were no longer present and thus did not
support the presence of the 2.07 ± 0.002 hour period. We also
point out multiple power spectra and Θ minima around the reported rotation periods in Figure 3 (left and right). These peaks
and minima are present as a result of spectral leakage, which is
due to large gaps in the data between consecutive observations.
Binary system rotation periods are an important parameter in
the diagnosis of the system’s orbital coplanarity. A radius estimate, in addition to a v sin i measurement of a system can only
loosely infer spin-orbital alignments - the rotation period is an
essential variable in the calculation of ‘v’, which makes its discovery all the more pertinent. Furthermore, the period reported
here has allowed us to infer which period matches each binary
member, and thus allowed for an effective estimate of masses
and radii, which we discuss in the following section.
Article number, page 5 of 9
Phase Dispersion Minimization
Power spectrum for 2MASS J0746+2000A (all epochs)
1
450
p = 7.232−1 x 24 = 3.32 hours
400
0.95
350
0.9
Theta statistic ( Θ)
Power
300
250
200
0.85
0.8
150
0.75
100
α = 5.7ε−07
α = 0.0027
α = 0.3173
50
0
0
5
10
Frequency (days −1)
15
20
0.7
p = 3.32 hours
0.65
3
3.2
3.4
3.6
Period (hours)
3.8
4
Fig. 3.
[LEFT] Lomb-Scargle periodogram of all 2MASS J0746+20A epochs, calculated from the combined dataset from Figure 1. The
three red dashed-dotted horizontal lines (top-to-bottom) represent a 5σ, 3σ and 1σ false-alarm probability of the peaks, as determined by the LS
periodogram algorithm. The x-axis is plotted in days−1 , since each epoch was time-stamped in units of Heliocentric Julian Days (HJD) - for period
accuracy over such a lengthy baseline. [RIGHT] Phase Dispersion Minimization (PDM) plot of the entire ∼2 year observation baseline. Here we
plot of the PDM ‘Theta (Θ) statistic’ vs. Period (in hours), where the routines minimize the value at 3.32 hours.
4. An estimate of age, mass and radius
By adopting the established total system mass of 0.151 ± 0.003
M⊙ , the photometric J H K measurements and bolometric luminosity measurements of Konopacky et al. (2010), in addition
to the lack of detected lithium in the binary dwarf’s spectrum
(Bouy et al. 2004), we were able to place constraints on the evolutionary models of Chabrier et al. (2000), in determining the
mass range and radii of each component. We make no initial assumptions for the age of the system, however young ages were
ruled out based on this absence of lithium. Thus we had three
measured quantities to estimate the mass track (i.e. J H K colors, Lbol , and Li=0). Once we identified a range of ages that did
not contain lithium for either component, we interpolated over a
range of masses based on the correlation between the J H K colors of Konopacky et al. (2010), and those of the Chabrier et al.
(2000) models, and then over the bolometric luminosities - thus
establishing the best agreement between each quantity. Most
importantly, assuming the stars are coeval, the masses were ultimately constrained since the sum of each component’s mass
could not be more than the estimated total mass of 0.151 ± 0.003
M⊙ for the system.
We find an age of ∼1 - 1.5 Gyr for the binary based on this assessment, as well as individual mass estimates of 0.078 ± 0.004
M⊙ and 0.073 ± 0.004 M⊙ for 2MASS J0746+20A and 2MASS
J0746+20B, respectively. These mass estimates are largely consistent with Bouy et al. (2004), and in good agreement with Gizis
& Reid (2006) and Konopacky et al. (2010), and furthermore
infer that each component lies at, or just below, the substellar
boundary, supporting the prediction of a low mass star classification for the secondary member (Gizis & Reid 2006). The difference in rotational velocity between these stars is most intriguing,
Article number, page 6 of 9
considering the similarity of component mass estimates. Perhaps there is a larger difference in mass between each member
than what has been inferred from the evolutionary models in the
above studies.
By contrast, an age of ∼1 - 1.5 Gyr identifies the system as
a much older binary dwarf than originally predicted by Bouy
et al. (2004), who found the system to be ∼150 - 500 Myr old.
This is a large discrepancy. However, we point out that Bouy et
al. (2004) put forward these ages despite the absence of lithium
in the binary’s spectrum - which is expected to be present for
stars of this age. We note however, that some studies suggest
that lithium can indeed be depleted at younger low mass star
ages, due to the effect of episodic accretion (Baraffe & Chabrier
2010). Although the absence/presence of lithium was originally
used as a test for sub-stellarity (e.g. Rebolo et al. (1992)), it is
consequently not necessarily a robust indicator of stellar age. A
much younger age is inconsistent with the surface gravity estimates of Schweitzer et al. (2001), who compared spectra to the
models of Allard et al. (2001). Their temperature estimates do
agree with those of Bouy et al. (2004), but their estimated gravity is too high for a ∼150 - 500 Myr old object. These gravity
estimates (from high-resolution spectra: log g ∼5.0 [K I λ7685]
- 5.5 [Ca I λ5673]; and from low-resolution spectra: log g ∼ 6.0)
are more consistent with the ages we infer in this work.
Therefore, adopting the estimates of age and mass above
places each star just below 1 R J , as shown in Figure 4, where
we estimate radii of 0.99 ± 0.03 R J for 2MASS J0746+20A and
0.96 ± 0.02 R J for 2MASS J0746+20B. These predictions are
consistent with those of Konopacky et al. (2010). However, the
radius estimate for 2MASS J0746+20A is inconsistent with the
estimate of Berger et al. (2009), due to their assumption of a v
sin i of 27 km s−1 for 2MASS J0746+20A, and a period of 2.07
L. K. Harding et al.: Spin-orbit alignment in the very low mass binary regime
1.35
log 8.2 − 8.7 years
log 9.0 years
log 9.1 years
log 9.2 years
1.3
sin iA =
1.25
Radius (RJ)
1.2
1.15
Range given by Bouy et al. (2004)
for ages of ~150 − 500 Myr
1.1
1.05
1
2MASS J0746+20A
0.95
2MASS J0746+20B
0.9
0.05
0.055
0.06
0.065
0.07
0.075
0.08
0.085
0.09
0.095
Mass (MO. )
Fig. 4.
Isochrones of radii (R J ) vs. mass (M⊙ ) for 2MASS
J0746+20AB, for ages of ∼1 (log 9.0; red), ∼1.25 (log 9.1; blue) and
∼1.5 Gyr (log 9.2; green), derived in this work by using the evolutionary models of Chabrier et al. (2000). The errors associated with the
mass and radii estimates are based on the constraints of the Chabrier
et al. (2000) models for a given age, as outlined in Section 4. We also
include evolutionary tracks for ages of ∼150 - 500 Myr (log 8.2 - 8.7
years) as predicted by Bouy et al. (2004). These are shown by the grey
dash-dotted lines. Based on the work this paper, the binary appears to
be a much older age than the range predicted by Bouy et al. (2004),
which require much larger radii for both components.
± 0.002 hours. This rotational velocity is approximately 30%
larger than the established v sin i of 19 ± 2 km s−1 (Konopacky
et al. 2012), thus incorrectly placing their radius estimate much
lower than those predicted by evolutionary models, at 0.78 ± 0.1
R J . Our newly established radii are much more consistent with
model predictions, based on our identifying the primary component as the non-radio pulsing star, and strong evidence supporting a spin-orbit alignment to within 10 degrees for both stars.
See Table 1 for a summary of system properties, including our
estimates of surface gravity.
5. Inferred spin-orbit alignment
The discovery of the rotation period of 2MASS J0746+20A
in this paper, in addition to the spectroscopic observations of
Konopacky et al. (2012) and the radio observations of Berger
et al. (2009), yield a v sin i ratio of 0.57+0.13
−0.10 , and a period ra.
By
adopting
these
v
sin
i and rotation period
tio of 0.62+0.02
−0.03
measurements, the estimated radii in this work indicate that both
components in the binary have their equatorial spin axis aligned
with the orbital plane to within 10◦ - consistent with the trends
for solar-type binaries with separations ≤ 40 AU (Hale 1994).
Is this evidence for a common formation pathway across a wide
swath of stellar masses? Our best fit indicates that the axes are
orientated with respect to the observer at 32 ± 4◦ and 36 ± 4◦ ,
respectively (Figure 5 and Figure 6). If the measurements are
taken at face value, both objects could also be perfectly coplanar,
implying that evolutionary models over-predict their radii for an
assumed mass and age. The inclination angle of the equatorial
spin axis of each component, sin iA and sin iB , were derived as
follows:
v sin iA
v sin iB
and sin iB =
,
2 · π · rA /PA
2 · π · rB /PB
(1)
where v is the rotational velocity, r is the radius and P is the
period. We estimated the uncertainly in equatorial inclinations
above by using Monte Carlo simulations. Each of the variables
in Equation 1 had known errors as calculated by Konopacky et al.
(2012) (v sin i), this work (rA and rB , PA ) and Berger et al. (2009)
(PB), respectively. For each of these variables, we generated 105
copies, where each iteration had random noise (normally distributed) added based on the reported errors. For each of the 105
simulations, we measure the standard deviation of the variables
and then calculated a mean value for each. This mean standard
deviation over all of the simulations is thus taken as the uncertainty (randomly calculated) in each quantity and applied to the
final calculation of the equatorial inclinations in Equation 1 (as
shown in Table 2 and Figure 5 & Figure 6).
Finally, we highlight a geometric effect with respect to system spin axes inclinations. The inclinations of the rotation axes
of each component to our line of sight may be equal, but this
does not necessarily imply that the orbital planes are coplanar.
For example, two edge-on spins can still be orthogonal on the
sky, and could be coincidentally equal.
6. Discussion: implications for formation
This is the first study to assess the orbital alignment properties
of a VLM binary system. The observed orbit-spin alignment
is consistent with several different formation pathways. There
are numerous binary (and multiple) star system formation theories, the most prominent of which are turbulent core fragmentation, disk fragmentation, and formation via dynamical interactions both during the main accretion phase and after, e.g. (Hoyle
1953; Adams et al. 1989; Laughlin & Bodenheimer 1994; Bonnell 1994a; Tsuribe & Inutsuka 1999; Padoan & Nordlund 2002;
Bate et al. 2003; Bonnell et al. 2004; Clark et al. 2008; Kratter
et al. 2010; Kratter 2011). A nearly coplanar system might result from any of these models, although formation via a dynamical interaction at late times is unlikely; tidal interactions are too
weak at these distances, so alignment would be coincidental.
Disk fragmentation naturally produces aligned systems, and
can also drive components towards equal mass (Adams et al.
1989; Laughlin & Bodenheimer 1994; Bonnell 1994a; Bonnell
& Bate 1994b; Kratter et al. 2010). However, disk fragmentation is more likely for higher mass systems (>1 M⊙ ) (Kratter et
al. 2008; Offner et al. 2010). Notably, these models have not
been extended down to the brown dwarf regime. Stamatellos
& Whitworth (2009) proposed an alternative disk fragmentation
scenario, where binary brown dwarfs are born within the disk
of a more massive star. This scenario is somewhat inconsistent
with 2MASS J0746+20AB, as it produces only very high eccentricity binaries (they predict e > 0.6 compared to this system’s
e = 0.487 (Konopacky et al. 2010)). It is also unclear whether
ejection would preserve spin alignment.
Both core fragmentation and competitive accretion can also
produce aligned systems. In the former scenario, fragments may
share the core’s net angular momentum vector (Matsumoto &
Hanawa 2003). We note that while equal inclinations is a necessary condition for aligned angular momenta, it is not a sufficient
one, unless the inclinations are zero. Recent work by Jumper
& Fisher (2012) has shown that a straightforward extrapolation
of the turbulent core fragmentation model to lower masses naturally reproduces the separation distribution of brown dwarfs.
Article number, page 7 of 9
Spin−orbit inclination angle, Θ (degrees)
150
50
60
70
80
90
100
110
120
Spin−orbit inclination angle, Θ (Degrees)
130
50
60
70
80
90
100
110
120
130
4.0
v Sin i
v Sin i
2
−1
3
J
Jupiter radii (R )
100
1.5
36
0.99
10
20
30
56
0.96
0.5
36 +/− 4
0.5
41.8
50
60
70
80
90
Orbital inclination angle, i (degrees)
Fig. 5. Equatorial rotational velocity (y-axis, left) vs. inclination angle
of the orbital plane (x-axis, bottom). The measured v sin i of 2MASS
J0746+20A is shown by the red solid curve. The dashed lines in all
cases represent the associated errors. The green vertical solid line (xaxis, top) highlight the alignment of the spin-orbit axes. Y-axis, right
corresponds to the radius of the dwarf (RJ = ∼69550 km) as calculated
in this work, where we have marked the estimated radius of 0.99 ±
0.03 RJ . The black horizontal lines show the corresponding equatorial
velocity of 36 ± 4 km s−1 . Crucially, this equatorial velocity and the associated alignment only apply for the period of 3.32 ± 0.15 hours and a
radius of 0.99 ± 0.03 RJ . The measured equatorial inclination of 32 ± 4
degrees is marked with the downward arrow, which is within 10 degrees
of the orbital inclination angle of 41.8 ± 0.5 degrees (Konopacky et al.
2012).
However, an analysis by Dupuy & Liu (2011), finds that the eccentricity distribution of ultracool dwarf binaries is statistically
distinct from that of solar type systems, and more consistent
with the clustered, competitive accretion model of Bate (2009,
2012). Even if such systems are born misaligned, tidal torquing
between disks can re-align close systems (see Lubow & Ogilvie
(2000)). Interaction with a circumbinary disk, as seen in Bate
(2012), might also align stars with initially random orientations.
As noted above, magnetic interactions, which might be even
stronger for fully convective stars, can also alter spin-orbit alignment. The αω-dynamo, thought to be responsible for the generation of magnetic fields of main sequence stars, no longer functions for fully convective objects. Magnetic fields many orders
of magnitude greater in strength than solar-type systems have
been shown to exist in the low mass star regime (Hallinan et
al. 2007; Reiners & Basri 2007; Antonova et al. 2008; Berger
et al. 2009). These fields could have a physical effect on alignment. 2MASS J0746+20AB, at a separation of only ∼2.7 AU
possesses a large-scale ∼1.7 kG magnetic field (Antonova et al.
2008). Indeed, a strongly misaligned system would be more indicative of dynamical processing.
Although alignment in 2MASS J0746+20AB cannot be used
to distinguish between various formation models, more measurements of similar systems might elucidate trends with mass. For
Article number, page 8 of 9
1.5
o
32 +/− 4o
0
100
Jupiter radii (RJ)
Equatorial rotational velocity (km s )
2MASS J0746+20B: 33 +/− 3 km s−1
−1
Equatorial rotational velocity (km s )
2MASS J0746+20A: 19 +/− 2 km s−1
0
10
20
30
41.8
50
60
70
80
90
Orbital inclination angle, i (degrees)
Fig. 6.
Equatorial rotational velocity (y-axis, left) vs. inclination
angle of the orbital plane (x-axis, bottom) for 2MASS J0746+20B. As
before, the measured v sin i of 2MASS J0746+20B is represented by
the blue solid curve, where dash-dotted lines represent the errors in each
measurement. We highlight our estimated radius of 0.96 ± 0.02 RJ (yaxis, right), which, together with the radio period of 2.07 ± 0.002 hours
(Berger et al. 2009), infers an equatorial inclination angle of 36 ± 4
degrees.
example, one might expect the separation at which systems transition from aligned to misaligned to be smaller for lower mass
systems, particularly if magnetic fields play a role. Comparing
the orbital properties of stars across the mass spectrum may indicate where different formation pathways dominate.
7. Conclusion
We report on the orbital coplanarity of the L dwarf tight binary
2MASS J0746+20AB. Recently, high-precision dynamical mass
and individual rotation velocity measurements were obtained
(Bouy et al. 2004; Konopacky et al. 2010, 2012), as well a rotation period for one component (Berger et al. 2009). We present
the discovery of the rotation period from the primary component
in this work. From these data, we infer that the binary orbital
plane is oriented perpendicular to the stellar spin axes to within
10◦ . Such alignment has previously been observed in studies
of solar-type binaries (Hale 1994), and also more massive stars
(Wheelwright et al. 2011). This work is the first direct evidence
of spin-orbit alignment in the VLM binary regime. We outline
the numerous binary formation models that are consistent with
the observed alignment. Further theoretical work, and a larger
sample of VLM systems (a campaign that we have already commenced), will place tighter constraints on the most likely formation pathways.
The discovery of the 3.32 ± 0.15 hour period implies that the
radio period of 2.07 ± 0.002 hours is associated with the secondary star and not the primary star, as was suggested in Berger
et al. (2009). This in turns refutes the claimed radius of 0.78
L. K. Harding et al.: Spin-orbit alignment in the very low mass binary regime
± 0.1 R J for 2MASS J0746+20A. We find that the primary and
secondary have radii of 0.99 ± 0.03 R J and 0.96 ± 0.02 R J respectively.
Acknowledgements. This work was largely carried out under the National University of Ireland Traveling Studentship in the Sciences (Physics). LKH gratefully acknowledges the support of the Science Foundation Ireland (Grant Number 07/RFP/PHYF553) as well as Dr. Mark Lang of NUI Galway. We thank
the VATT team for their help and guidance. Finally, we would like to thank the
referee for their careful reading of our work, and for their valuable input and
suggestions on how to improve this manuscript.
References
Abt, H. A. & Levy, S. G. 1976, ApJS, 30, 273
Adams, F. C. & Ruden, S. P. & Shu, F. H. 1976, ApJ, 347, 959A
Albrecht, S. & Reffert, S. & Snellen, I. A. G. & Will, J. N. 2009, Nature, 461,
373
Albrecht, S. & Will, J. N. & Carter, J. A. & Snellen, I. A. G. & de Mooij, E. J.
W. 2011, ApJ, 726, 68
Allard, F. & Hauschildt, P. H. & Alexander, D. R. & Tamanai, A. & Schweitzer,
A. 2001, ApJ, 556, 357
Antonova, A. & Doyle, J. G. & Hallinan, G. & Bourke, S. & Golden, A. 2008,
ApJ, 487, 317
Bailer-Jones, C. A. L. 2004, A&A, 419, 703
Baraffe, I. & Chabrier, G. 2010, A&A, 521, 44
Bate, M. R. & Bonnell, I. A. & Bromm, V. 2003, MNRAS, 339, 577
Bate, R. B. 2009, MNRAS, 392, 590
Bate, R. B. 2012, MNRAS, 419, 3115
Berger, E. & Rutledge, R. E. & Phan-Bao, N. et al. 2009, ApJ, 695, 310
Bodenheimer, P. 1978, ApJ, 224, 488
Bonnell, I. A. 1994, MNRAS, 269, 837B
Bonnell, I. A. & Bate, M. R. 1994, MNRAS, 271, 999
Bonnell, I. A. & Vine, S. G. & Bate, M. R. 2004, MNRAS, 349, 735
Bouy, H. & Duchêne, G. & Köhler, R. et al. 2004, A&A, 423, 341
Burgasser, A. J. & Reid, I. N. & Siegler, N. et al. 2007, in Protostars and Planets
V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson, AZ: Univ. Arizona Press),
427
Clarke, F. J. & Oppenheimer, B. R. & Tinney, C. G. 2002, MNRAS, 335, 1158
Clark, P. C. & Bonnell, I. A. & Klessen, R. S. 2008, MNRAS, 386, 3
Chabrier, G. & Baraffe, I. & Allard, F. & Hauschildt, P. 2000, ApJ, 542, 464
Close, L. M. & Siegler, N. & Freed, M. & Biller, B. 2003, ApJ, 587, 407
Dahn, C. C. & Harris, H. C. & Vrba, F. J. et al. 2002, AJ, 124, 1170
Donati, J. F. & Morin, J. & Petit, P. et al. 2008, MNRAS, 390, 545
Dupuy, T. J. & Liu, M. C. & Bowler, B. P. et al. 2010, ApJ, 721, 1725
Dupuy, T. J. & Liu, M. C. 2011, ApJ, 733, 122
Fekel Jr., F. C. 1981, ApJ, 246, 879
Gizis, J. E. & Reid, N. I. 2006, AJ, 131, 638
Hale, A. 1994, AJ, 107, 1
Hallinan, G. & Bourke, S. & Lane, C. et al. 2007, ApJ, 663, L25
Hallinan, G. & Antonova, A. & Doyle, J. G. et al. 2008, ApJ, 684, 644
Harding, L. K. & Hallinan, G. & Boyle, R. P. et al. 2013, ApJ, submitted
Hoyle, F. 1953, ApJ, 118, 513
Jensen, E. L. N. & Mathieu, R. D. & Donar, A. X. & Dullighan, A. 2004, ApJ,
600, 789
Jumper, P. H. & Fisher, R. T. 2012, ApJ, submitted (arXiv1206.1045)
Konopacky, Q. M. & Ghez, A. M. & Barman, T. S. et al. 2010, ApJ, 711, 1087
Konopacky, Q. M. & Ghez, A. M. & Fabrycky, D. C. et al. 2012, ApJ, in press
Kratter, K. M. & Matzner, C. D. & Krumholz, M. R. 2008, ApJ, 681, 375
Kratter, K. M. & Matzner, C. D. & Krumholz, M. R. & Klein, R. I. 2010, ApJ,
708, 1585
Kratter, K. M. 2011, in ASP Conf. Ser. 447, Evolution of Compact Binaries, ed.
L. Schmidtobreick & M. R. Schreiber & C. Tappert, 47
Lada, C. J. 2006, ApJ, 640, L63
Laughlin, G. & Bodenheimer, P. 1994, ApJ, 436, 335L
Li, P. S. & Norman, M. L. & Mac Low, M. M. & Heitsch, F. 2004, ApJ, 605, 800
Lomb, N. R. 1976, ApJS, 39, 447
Lubow, S. H. & Ogilvie, G. I. 2000, ApJ, 538, 326
Matsumoto, T. & Hanawa, T. 2003, ApJ, 595, 913
Mestel, L. 1977, In IAU Symposium, 75, 213
Monin, J. & Ménard, F. & Peretto, N. 2006, A&A, 446, 201
Nakajima, T. & Oppenheimer, B. R. & Kulkarni, S. R. et al. 1995, Nature, 378,
463
Offner, S. S. R. & Kratter, K. M. & Matzner, C. D. & Krumholz, M. R. & Klein,
R. I. 2010, ApJ, 725, 1485
Padoan, P. & Nordlund, A. 2002, ApJ, 576, 870
Raghavan, D. & McAlister, H. A. & Henry, T. J. et al. 2010, ApJS, 190, 1
Rebolo, R. & Martin, E. L. & Magazzu, A. 1992, ApJ, 389, L83
Reid, I. N. & Burgasser, A. J. & Cruz, K. L. et al. 2001, AJ, 121, 1710
Reid, I. N. & Kirkpatrick, J. D. & Liebert, J. et al. 2002, AJ, 124, 519
Reiners, A. & Basri, G. 2007, ApJ, 656, 1121
Reiners, A. & Basri, G. 2008, ApJ, 684, 1390
Scargle, J. D. 1982, ApJ, 263, 835
Schweitzer, A. & Gizis, J. E. & Hauschildt, P. H. & Allard, F. & Reid, I. N. 2001,
ApJ, 555, 368
Sheehan, B. & Butler, R. F., in AIP Conference Proceedings, Volume 984, pp.
162 - 167 (2008)
Stamatellos, D. & Whitworth, A. P. 2009, MNRAS, 392 413
Stellingwerf, R. F. 1982, ApJ, 224, 953
Tsuribe, T. & Inutsuka, S. I. 1999, ApJ, 526, 307
Vrba, F. J. & Henden, A. A. & Luginbuhl, C. B. et al. 2004, AJ, 127, 2948
Weis, E. W. 1974, ApJ, 190, 331
West, A. A. & Morgan, D. P. & Bochanski, J. J. et al. 2011, AJ, 141, 97
Wheelwright, H. E. & Vink, J. S. & Oudmaijer, R. D. & Drew, J. E. 2011, A&A,
532, A28
Article number, page 9 of 9
View publication stats