Tidal effects on brown dwarfs: application to the eclipsing binary 2MASS J053521840546085
The anomalous temperature reversal in the context of tidal heating
R. Heller^{1}  B. Jackson^{2}  R. Barnes^{3,4}  R. Greenberg^{5}  D. Homeier^{6}
1  Hamburger Sternwarte (Universität Hamburg), Gojenbergsweg 112,
21029 Hamburg, Germany
2  Lunar and Planetary Laboratory, University of Arizona, Tucson, AZ
85721, USA
3  University of Washington, Dept. of Astronomy, Seattle, WA 98195,
USA
4  Virtual Planetary Laboratory, NASA, USA
5  Lunar and Planetary Laboratory, University of Arizona, Tucson, AZ
85721, USA
6  Institut für Astrophysik, GeorgAugustUniversität Göttingen,
FriedrichHundPlatz 1, 37077 Göttingen, Germany
Received 4 July 2009 / Accepted 8 February 2010
Abstract
Context. 2MASS J053521840546085
(2M053505) is the
only known eclipsing brown dwarf (BD) binary, and so may serve as a
benchmark for models of BD formation and evolution. However,
theoretical predictions of the system's properties seem inconsistent
with observations: i) the more massive (primary) component is
observed to be cooler than the less massive (secondary) one; ii) the
secondary is more luminous (by 10^{24} W)
than expected. Previous explanations for the temperature reversal have
invoked reduced convective efficiency in the structure of the primary,
connected to magnetic activity and to surface spots, but these
explanations cannot account for the enhanced luminosity of the
secondary. Previous studies also considered the possibility that the
secondary is younger than the primary.
Aims. We study the impact of tidal heating to the
energy budget
of both components to determine if it can account for the observed
temperature reversal and the high luminosity of the secondary. We also
compare various plausible tidal models to determine a range of
predicted properties.
Methods. We apply two versions of two different,
wellknown
models for tidal interaction, respectively: i) the
``constantphaselag'' model; and ii) the
``constanttimelag''
model and incorporate the predicted tidal heating into a model of BD
structure. The four models differ in their assumptions about the
rotational behavior of the bodies, the system's eccentricity and
putative misalignments
between the bodies' equatorial planes and the orbital plane of the
system.
Results. The contribution of heat from tides in
2M053505 alone
may only be large enough to account for the discrepancies between
observation and theory in an unlikely region of the parameter space.
The tidal quality factor
of BDs would have to be 10^{3.5} and the secondary
needs a spinorbit misalignment of .
However, tidal synchronization time scales for 2M053505 restrict the
tidal dissipation function to
and rule out intense tidal heating in 2M053505. We provide the first
constraint on Q for BDs.
Conclusions. Tidal heating alone is unlikely to be
responsible
for the surprising temperature reversal within 2M053505. But an
evolutionary embedment of tidal effects and a coupled treatment with
the structural evolution of the BDs is necessary to corroborate or
refute this result. The heating could have slowed down the BDs'
shrinking and cooling processes after the birth of the system 1 Myr
ago, leading to a feedback between tidal inflation and tidal heating.
Observations of old BD binaries and measurements of the
RossiterMcLaughlin effect for 2M053505 can provide further
constraints on .
Key words: celestial mechanics  binaries: eclipsing  stars: evolution  stars: individual: 2MASSJ053521840546085  brown dwarfs
1 Introduction
2MASS J053521840546085 (2M053505) is a benchmark object for brown dwarf (BD) science since it offers the rare opportunity of independent radius and mass measurements on substellar objects. The observed values constrain evolutionary and structural models (Chabrier et al. 2007; D'Antona & Mazzitelli 1997; Baraffe et al. 2002; Chabrier & Baraffe 2000; Baraffe et al. 1998). 2M053505 is located in the Orion Nebulae, a starforming region with an age of 1 (0.5) Myr. If both components formed together, as commonly believed, then this system allows for effective temperature ( ) and luminosity (L) measurements of two BDs at the same age.
However, this system is observed to have an unexpected temperature reversal (Stassun et al. 2006), contravening theoretical simulations: the more massive component (the primary) is the cooler one. From the transit light curve, the ratio of the effective temperatures can be accurately determined to (Mohanty et al. 2009; Gómez Maqueo Chew et al. 2009). From spectroscopic measurements then, the absolute values can be constrained. The primary, predicted to have K (Baraffe et al. 1998), has an observed value of 2700 K, whereas the surface temperature of the secondary, predicted to be K, is most compatible with K.
One explanation for the temperature discrepancies is suppression of convection due to spots on the surface of the primary. If a portion of a BD's surface is covered by spots, its apparent temperature will be reduced, resulting in an increase in the estimated radius in order for the measured and expected luminosities to agree (Chabrier et al. 2007). With a spot coverage of 3050% and a mixing length parameter most of the mismatches between predicted and observed radii for lowmass stars (LMS) can be explained (Ribas et al. 2008). Observations of spots on both of the 2M053505 components (Gómez Maqueo Chew et al. 2009), as inferred from periodic variations in the light curve, and measurements on the H line of the combined spectrum during the radial velocity maxima (Reiners et al. 2007) suggest that enhanced magnetic activity and the accompanying spots on the primary indeed play a key role for its temperature deviation. But even if the spot coverage on the primary serves as an explanation for the primary's reduced , the secondary's luminosity overshoot of W, as compared to the Baraffe et al. (1998) models, suggests some additional processes may be at work.
The temperature reversal between the primary and secondary may result from a difference between their ages. The secondary could be 0.5 Myr older than the primary, as proposed by Stassun et al. 2007 (see also D'Antona & Mazzitelli 1997). A difference of 0.5 Myr could allow the secondary to have converted the necessary amount of gravitational energy into heat^{}, which would explain its luminosity excess. But evolutionary models are very uncertain for ages 1 Myr (Wuchterl 2005; Marley et al. 2007; Baraffe et al. 2002; Mohanty et al. 2007) and, in any case, the age determination and physical natures of these very young objects is subject to debate (Stassun et al. 2009,2008). Furthermore, the mutual capture of BDs and LMS into binary systems after each component formed independently is probably too infrequent to account for the large number of eclipsing LMS binaries with either temperature reversals or inflated radii (Morales et al. 2009; Guenther et al. 2001; Ribas et al. 2008; Çakirli et al. 2009; Coughlin & Shaw 2007).
Table 1: Orbital and physical parameters of 2M053505.
Here, we consider the role that tidal heating may play in
determining the temperatures of the BDs. In Table 1
we show the parameters of 2M053505 necessary for our calculations. The
computed energy rates will add to the luminosity of the BDs in some way
(Sect. 2.3)
and will contribute to a temperature deviation compared to the case
without a perturbing body (Sect. 3). All these
energy rates must be seen in the context of the luminosities of the
BDs: W
(luminosity of the primary) and W
(luminosity of the secondary). At a distance a
to the primary component, its luminosity is distributed onto a sphere
with area .
The secondary has an effective  i.e. a 2Dprojected  area of .
With
as the flux of the primary at distance a, the
irradiation from the primary onto the secondary
is thus given by
Using that equation, we calculate the mutual irradiation of the BDs: W and W. These energy rates are two and three orders of magnitude lower, respectively, than the observed luminosity discrepancy. Hence, we assume that mutual irradiation can be ignored. This simplification is in contrast to the cases of the potentially inflated transiting extrasolar planets WASP4b, WASP6b, WASP12b, and TrES4, where stellar irradiation (Ibgui et al. 2010) dominates tidal heating by several magnitudes.
Various tidal models haven been used to calculate tidal heating in exoplanets (Barnes et al. 2009; Bodenheimer et al. 2001; Jackson et al. 2008a,b), which may in fact be responsible for previous discrepancies between interior models and radii of transiting exoplanets (Ibgui & Burrows 2009; Jackson et al. 2008a,b). This success in exoplanets motivates our investigation into BDs. While many different tidal models are available, there is no consensus as to which is the best. For this reason, we apply a potpourri of wellestablished models to the case of 2M053505 in order to compare the different results. As we show, tidal heating may account for the temperature reversal and it may have a profound effect on the longerterm thermal evolution of the system.
The coincidence of , with as the orbital and P_{1} primary's rotation period, has been noted before but we assume no resonance between the primary's rotation and the orbit for our calculations. These resonances typically occur in systems with rigid bodies where a fixed deformation of at least one body persists, such as in the SunMercury configuration with Mercury trapped in a 3/2 spinorbit resonance. We assume that, in the context of tides, BDs may rather be treated as fluids and the shape of the body is not fixed.
With this paper, we present the first investigation of tidal interaction between BDs. In Sect. 2 we introduce four models for tidal interaction and discuss how we convert the computed energy rates into an increase in effective temperature. Section 3 is devoted to the results of our calculations, while we deal with the observational implications in Sect. 4. We end with conclusions about tidal heating in 2M053505, and in BDs in general, in Sect. 5.
2 Tidal models
Two qualitatively different models of tidal dissipation and evolution have been developed over the last century: The ``constantphaselag'' (Goldreich & Soter 1966; FerrazMello et al. 2008; Wisdom 2008, Wis08 and FM08 in the following), and the ``constanttimelag'' model (Hut 1981, Hut81 in the following). In the former model, the forces acting on the deformed body are described by a superposition of a static equilibrium potential and a disturbing potential (FM08). The latter model assumes the time between the passage of the perturbing body overhead and the passage of the tidal bulge is constant. Although both models have been used extensively, it is not clear which model provides a more accurate description of the effects of tides, so we apply formulations of both models.
In the ``constantphaselag'' model of FM08, quantitative expressions have been developed to second order in eccentricity e while the others include also higher orders. Higher and higher order expansions require assumptions about the dependence of a body's tidal response to an increasing number of tidal frequencies, which involves considerable uncertainty. Therefore higher order expansions do not necessarily provide more accuracy (FM08; Greenberg 2009). In the ``constantphaselag'' model of Wis08, expressions in e are developed to order. The ``constanttimelag'' model of Hut81 does not include possible obliquities, while an enhanced version of that model by Levrard et al. (2007) (Lev07) does.
Tidal dissipation in BDs has not been observed or even considered previously, and hence, neither model should take precedence when calculating their tidal dissipation, especially since neither tidal model is definitive (Greenberg 2009). As our investigation is the first to consider tidal effects on BDs, we will employ several applicable, previously published models to 2M053505. By surveying a range of plausible models and internal properties, usually encapsulated in the ``tidal dissipation function'' Q (Goldreich & Soter 1966), we may actually be able to determine which model is more applicable to the case of BDs  assuming, of course, that tidal dissipation contributes crucially to the observed temperature inversion.
2.1 Constant phase lag
2.1.1 Tidal model #1
The potential of the perturbed body can be treated as the
superposition of periodic contributions of tidal frequencies at
different phase lags and the expression for the potential can be
expanded to first order in those lags (FM08). Those phase lags
of the ith body that we will need for our equations
are given by
where is the algebraic sign of x, thus , is the orbital frequency and are the rotational frequencies of the primary (i = 1) and secondary (i = 2), P_{i} being their rotational periods. The tidal frequencies are functions of the tidal quality factor Q of the deformed object, which parametrizes the object's tidal response to the perturber. It is defined as
(3) 
where E_{0} is the maximum energy stored in the tidal distortion and the integral over the energy dissipation rate dE/dt is the energy lost during one orbital cycle (Goldreich & Soter 1966). Although Ogilvie & Lin (2004) conclude that tidal dissipation rates of giant planets are not adequately represented by a constant Qvalue, many parameterized tidal models rely on this quantity. Measurements of the heat flux from Jupiter's moon Io during the flyby of the Voyager 1 spacecraft, combined with a specific model of the history of the orbital resonance, allowed for an estimate for the quality factor of Jupiter to be (Yoder 1979) while Aksnes & Franklin (2001) used historical changes in Io's orbit to infer that is around 10^{5.3}. However, Greenberg et al. (2008) pointed out that is not ruled out (see also Ioannou & Lindzen 1993; Peale & Greenberg 1980). Tides raised by Neptune on its moons help to constrain the planet's quality factor to (Zhang & Hamilton 2008). For M dwarfs, is assumed to be of order 10^{5}, whereas for rigid bodies like Earth (Mardling & Lin 2004; Ray et al. 2001, and references therein). For BDs, however, Q is even more uncertain, thus we will handle it as a free parameter in our procedures.
FM08 allows for the tidal amplitude to be different from what it would be if the tideraising body were fixed in space. This concern is met by the dynamical Love number under the assumption that the tidally disturbed body had infinite time to respond. Without better knowledge of a body's response to tides, we assume the dynamical Love number is the same as the potential Love number of degree 2, k_{2}. For the gas planets of the solar system, this number has been calculated by Gavrilov & Zharkov (1977). BDs may rather be treated as polytropes of order (Baraffe, private communication). We infer the Love number from the relation (Mardling & Lin 2002) and use the tables of apsidal motion constants given in Brooker & Olle (1955). These authors provide numerical calculations for for a polytrope of . We find and thus . This places k_{2} for BDs well in the regime spanned by the gas giants of the solar system: Jupiter ( k_{2} = 0.379), Saturn ( k_{2} = 0.341), Uranus ( k_{2} = 0.104) and Neptune ( k_{2} = 0.127) (Gavrilov & Zharkov 1977).
Before we proceed to the equations for the tidal heating
rates, we sum
up those for the orbital evolution of the semimajor axis a,
the eccentricity e and the putative
obliquity .
The latter parameter is the angle between the equatorial plane of one
of the two bodies in a binary system and the orbital plane (Winn et al. 2005),
frequently referred to as spinorbit misalignment. We use
Eqs. (56), (60) and (61) from FM08 but our equations
for a
binary system with comparable masses need slight modifications since
both constituents contribute significantly to the evolution of a
and e. We add both the terms for the secondary
being the perturber of the primary (i = 1, j
= 2) and vice versa, since only spinorbit coupling is relevant,
whereas spinspin interaction can be neglected. This
results in
where is the dynamical Love number, M_{i} the mass and R_{i} the radius of the deformed BD, , with as the obliquity of the perturbed body, and are the tidal phase lags, given in Eq. (2).
The total energy that is dissipated within the perturbed body,
its
tidal energy rate, can be determined by summing the work done by tidal
torques (Eqs. (48) and (49) in FM08). The change in orbital
energy
of the i
body due to the j
body is given by
and the change in rotational energy is deduced to be
where G is Newton's gravitational constant. The total energy released inside the body then is
The greaterthan sign in this equation is true, since either and orbital energy is converted into rotational energy, or and the body is decelerated by a transfer of rotational energy into orbital energy. In both cases, the dynamical energy of the system is released within the distorted body. For , e.g., Eqs. (8) and (9) yield and .
The approach for the calculation of tidal energy rates with tidal model
#1 depends on processes due to nonsynchronous rotation via
and includes a putative obliquity
and terms of e up to the second order.
After inserting the orbital and rotational periods for 2M053505, these
equations reduce to
Interestingly, for these particular values of , and n, the S_{1}terms for cancel each other, so that it is not a function of , whereas does depend on .
2.1.2 Tidal model #2
The model of Wis08 includes terms in eccentricity up to the 8th order,
predicting higher tidal energy rates than for the equations of
model #1. Equations for the evolution of the orbital
parameters
are not given in Wis08. Furthermore, in his theory the perturbed body
is assumed to be synchronously rotating with the orbital period. Since
this is not the case for either of the BDs in 2M053505, the following
equations will only yield lower limits for the tidal heating. The tidal
heating rates are given by
with
where we used and
following the nomenclature of Hut (1981) and Wisdom (2008) as indicated. Furthermore, k_{2,i} is the potential Love number of degree 2 for the ith component of the binary system and is a measure of the longitude of the node of the body's equator on the orbit plane with respect to the pericenter of its orbit. In order to estimate the impact of in the last term in Eq. (13), we assume this impact to be as large as possible, , and compare it to the preceding terms. We find that for the case of 2M053505 the first three terms are of order 1, whereas the term connected to varies between 10^{2} and 10^{5}, depending on . These irrelevant contributions give us a justification to neglect the unknown values of in 2M053505 for our computations, facilitating the comparisons to the other models.
2.2 Constant time lag
2.2.1 Tidal model #3
Instead of assuming phase lags and superposition of frequencydependent
potentials, the ``equilibrium tide'' model by Hut (1981) invokes a
constant time lag
between the line joining the centers of the two bodies and the
culmination of the tidal bulge on the distorted object. With that
assumption, the model of Hut81 is mutually exclusive with the
assumption of a fixed angle lag (Goldreich
& Soter 1966): in general, a fixed time lag and a
fixed angle lag result in very different behaviors of the tidal bulge^{}.
As for the case of the ``constanttimelag'' model, we first sum up the
equations governing the behavior of the orbital evolution. With the
purpose of easing a comparison between Hut81's equations
(Eqs. (9)(11) therein) and Eqs. (4)(7) from this
paper for the theory of the ``constantphaselag'' model #1,
we transform the former into
with as the apsidal motion constant of the perturbed body (see Sect. 2.1.1), as the radius of gyration of the i body, which is defined by the body's moment of inertia , and
Hut81 then calculates the energy dissipation rate within a binary system, caused by the influence of one of the two bodies on the other, as the change in the total energy . Here, and are the orbital and rotational energies of the body (Eqs. (A28)(A35) in Hut81). For the tidal heating rates of the i constituent within the binary, this yields
where
Unfortunately, with these equations for the tidal energy rates model #3 neglects a potential obliquity of the body, which prevents us from a direct comparison with the other tidal models.
2.2.2 Tidal model #4
Lev07 extended Hut81's formula for the tidal energy rate to the case of
an object in equilibrium rotation^{}
and they included possible obliquities (see also Neron de Surgy & Laskar
1997), though they do not give the equations for the orbital
evolution. Lev07's equations are equivalent to
where
(22) 
The ``annual tidal quality factor'' is given as . Even though Lev07's equations invoke and their equations resemble those of the models with constant phase lag, their approach still assumes a constanttimelag. Since Lev07 do not explicitly connect their to the Q of FM08 (model #1) and Wis08 (model #2), we keep Q and as two different constants for our further treatment.
With these expansions, Eq. (21) involves terms in eccentricity up to order e^{8}. But since model #4 assumes tidal locking, i.e. is not a function of , this model also yields just a lower limit for the heating rates (Wisdom 2008).
2.3 Converting tidal heating into temperature increase
Now that we have set up four distinct models for the calculations of
the additional tidal heating term for the BDs, there are two physical
processes that will be driven by these energy rates: tidal inflation
and temperature increase. Let's take
as the luminosity of either of the two 2M053505 BDs that it would have
if it were a single BD and
and
as its corresponding radius and effective temperature. Then, by the
StefanBoltzmann law (Boltzmann 1884; Stefan 1879)
where is the StefanBoltzmann constant. The radial expansion in the binary case is given by and the temperature increase by . In its present state, the BD has a luminosity
where is some additional internal energy rate. Solving Eq. (24) for the temperature increase yields:
In the next step, we quantify the amount of tidal energy that is converted into internal energy, leading to an increase in effective temperature. Since we will use the virial theorem for an ideal, monoatomic gas to estimate the partition between internal and gravitational energy, we first have to assess the adequacy of treating the 2M053505 BDs as ideal gases. We therefore show the degeneracy parameter as a function of radius in Fig. 1 (Chabrier & Baraffe 2000; Baraffe, private comm.). Here, is the Boltzmann constant, T is the local temperature within the gas and is the Fermi energy of a partially degenerate electron gas with an electron Fermi temperature . With respect to and , g being the body's gravitational acceleration at the surface, the BD structure model corresponds to that of the primary, but with an age of 4.9 Myr. We find that for most of the BD, i.e. that portion of the structure in which the majority of the luminosity is released, is of order 1. This means that we may indeed approximate the BDs as ideal gases.
Figure 1: Degeneracy parameter (solid line) with model parameters similar to those of the 2M053505 primary and radiusintegrated luminosity L (dashed line) as a function of radius. To fit into the plot, L is normalized to 10. 

Open with DEXTER 
With the time derivative of the virial theorem for an ideal monoatomic
gas (Kippenhahn &
Weigert 1990, Sect. 3.1 therein),
(26) 
where is the temporal change in gravitational energy, we find that half of the additional tidal energy is converted into internal energy and the other half causes an expansion of the BD. There are currently no models for tidal inflation in BDs and the treatment is beyond the scope of this paper. Instead of including the modeled BD radii into Eq. (25) we avoid further uncertainties and fix (see Sect. 5 for a discussion of tidal inflation in the evolutionary context). The increase in effective temperature due to tidal heating then becomes
For we took the values predicted by the Baraffe et al. (1998) models (see Table 1).
Our neglect of tidal inflation makes this temperature increase an upper limit. Given that this neglect is arbitrary, we estimate how our constraints for and would change, if tidal inflation played a role in 2M053505. Comparing the observed radii of both BDs with the model predictions (see Table 1), radial expansions of 10% for the primary and 20% for the secondary seem realistic. Theoretical investigations of tidal heating on the inflated transiting planet HD209458b by (Ibgui & Burrows 2009) support an estimate of tidal inflation by 20%. As a test, we assumed that the secondary BD in 2M053505 is tidally inflated, where its radius in an isolated scenario would be 80% of its current value, i.e. in Eq. (25). In the noninflated scenario with , the BD would reach a temperature increase of dT = 60 K at and with model #2 (see Sect. 3.3). With the inflation, however, is needed to achieve the same heating at , whereas no obliquity at would yield significant heating. Thus, if tidal inflation in the secondary BD increases its radius by 20%, then the value for the dissipation function required to yield the same would be about 0.8 smaller in than in the case of no inflation. Therefore, the temperature we report in Sect. 3 may, at worst, correspond to that is smaller by 0.8.
3 Results
3.1 Orbital evolution
In order to get a rough impression of how far the orbital configuration of the system has evolved, we used the equations given in FM08, to compute the change of its eccentricity e and of a possible obliquity of the secondary within the last 1.5 Myr. Since this time span is the upper bound for the system's age, confined by its localization within the Orion Nebulae and indicated by comparison with BD evolutionary tracks, we thus get the strongest changes in e and . If any initial obliquity would be washed out already, could be neglected in the calculations of tidal heating. Furthermore, the measured eccentricity e could give a constraint to the tidal dissipation function Q. Computations based on the theory of ``constanttimelag'' yield qualitatively similar results.
Figure 2: Orbital evolution of 2M053505 after model #1 going back in time for 1.5 Myr. Left: eccentricity evolution. Depending on and on the age of the system, its initial eccentricity has not been smaller than 0.3133, which is 97.4% of its current value. Right: obliquity evolution of the secondary BD for three different values of Q_{2}. Simulations started at ``time = 0'' for and were evolved backwards in time. For there is no significant change in . For all the treated values of Q_{2}, the obliquity of the 2 M053505 secondary is still close to its natal state. 

Open with DEXTER 
For the evolution of e, we relied on Eq. (6). We took the observed eccentricity e = 0.3216 as a starting value and evolved it backwards in time. To evolve the system into the past, we changed the sign of the right side of the equation. Furthermore, we assumed that the quality factors Q_{1} and Q_{2} of the primary and secondary are equal, leading to and , because we are merely interested in a tentative estimate so far. This assumption should be a good approximation due to the similarity of the both components in terms of composition, temperature, mass, and radius.
The observed eccentricity of the system might give a constraint to the possible values for since depends on via . Certain regimes could be incompatible with the observed eccentricity of the system at a maximum age of 1.5 Myr, if these values would have caused the eccentricity to decay rapidly to 0 within this time. However, our simulations (Fig. 2) show that the system has not yet evolved very far for the whole range of and that the eccentricity of 2M053505 is in fact increasing nowadays. In this system, circularization does not occur. The observed eccentricity of 0.3216 consequently does not constrain . In this first estimate, we fixed all other parameters in time, i.e. we neglected an evolution of the semimajor axis a, of possible obliquities and we used constant radii R_{i} and rotational frequencies . We did this because we cannot yet incorporate the evolutionary behavior of the components' radii R_{i} in the context of tides and furthermore, there is no knowledge about possible misalignments between the orbital plane and the equatorial planes of the primary and secondary, respectively. A consistent evolution of R_{i}, however, is necessary to evolve as a function of and , as given by Eq. (4). Such a calculation was beyond the scope of this study.
The relative spingeometry of the two BD rotational axes with
respect
to the orbital plane and with respect to each other is unknown in
2M053505. Anyhow, we can estimate if a possible obliquity that once
existed for one of the BDs would still exist at an age of
1.5 Myr
or if it would have been washed out up to the present. We used a
numerical integration of Eq. (7) to evolve
forward in time (Fig. 2).
For the secondary's initial obliquity ,
we plot the state of
as a function of the quality factor Q_{2}
after an evolution of 1.5 Myr. We see that even for a very
small quality factor of 10^{3}
and high initial obliquities the secondary is basically in its natal
configuration today. Thus, it is reasonable to include a putative
misalignment of the secondary with respect to the orbital plane in our
considerations. As shown below, this is crucial for the calculations of
the tidal heating and the temperature reversal.
3.2 Tidal heating in 2M053505 with model #1
In Fig. 3, we show the results for the tidal heating rates as computed after tidal model #1. As given by Eq. (11), the tidal heating of the primary does not depend on a putative obliquity, whereas that of the secondary does. Using this model, we find that the luminosity gain of the secondary is, over the whole Q range, smaller than that of the primary, which mainly results from the relation . Figure 3 also shows that a growing obliquity shifts the gain in thermal energy towards higher values for a fixed Q_{2}. The observed overshoot of 10^{24} W in the secondary's luminosity can be reproduced with very small quality factors of and high obliquities up to .
Figure 3: Tidal heating after model #1. Left: (Primary) Projection of onto the  plane. The stepsize between contour lines is chosen to be in . Right: (Secondary) Projection of onto the  plane. Although there is a dependence on , the energy rates at a fixed value for the quality factor are smaller than those for the primary. 

Open with DEXTER 
Figure 4: Temperature increase after model #1. Left: (Primary) Projection of dT_{1} onto the  plane. For a significant temperature increase, Q_{1} would have to be much smaller than 10^{3.5}, but such a temperature increase is not observed in the primary. Right: (Secondary) Projection of dT_{2} onto the  plane. Even for very low values of Q_{2} and high obliquities the observed temperature increase cannot be reconstructed. For any given point in the  plane, , which does not support the observed temperature reversal. 

Open with DEXTER 
Figure 5: Tidal heating after model #2. Left: (Primary) Projection of onto the  plane. The stepsize between contour lines is chosen to be in . The tidal energy rates strongly depend on a putative obliquity, different from model #1. Right: (Secondary) Projection of onto the  plane. For the three models (#1, #2, and #4) invoking Q and , these rates are the highest of all for any given point in the  plane  for the primary as well as for the secondary. 

Open with DEXTER 
Figure 6: Temperature increase after model #2. Left: (Primary) Projection of dT_{1} onto the  plane. Right: (Secondary) Projection of dT_{2} onto the  plane. For any given location in the  plane, model #2 yields the strongest temperature increase compared to the other models  both for the primary and the secondary, respectively. For a given spot in Q space there is an inversion in temperature increase: , i.e. the less massive BD is heated more. 

Open with DEXTER 
In Fig. 4, we show the results for the temperature increase as per Eq. (27) with the tidal energy rates coming from model #1. These rates yield only a slight temperature increase for both constituents. Even for low Q values of order 10^{4} and high obliquities of the secondary, the heating only reaches values 10 K. We also see that the heating for the primary is computed to be greater than that for the secondary and no temperature reversal would be expected. If both BDs have the same Q values, then model #1 is unable to explain the temperature reversal. We cannot rule out a system in which, e.g., Q_{1} = 10^{5} and Q_{2} = 10^{3}, for which model #1 could explain the reversal. However, there is no reason to expect that similar bodies have Q values that span orders of magnitude. Hence, we conclude that model #1 can neither reproduce the luminosity overshoot of the secondary nor the system's temperature reversal.
3.3 Tidal heating in 2M053505 with model #2
This model yields the highest heating rates and hence temperature increases. The contrast between the absolute energy rates within the primary and the secondary is very small. In fact, for any given point in Q space, the heating rates differ only by (Fig. 5). The tidal energy rates of the secondary become comparable to the observed luminosity overshoot at and , where W. A comparison of the heating rates from model #2 with those of model #1 for either of the BDs shows that model #2 provides higher rates, with growing contrast for increasing obliquities.
The temperature increase arising from the comparable heating rates is inverted for a given spot on the plane. If both BDs had the same obliquity and the same dissipation factor, the secondary would experience a higher temperature increase. As presented in Fig. 6, the temperature increase after model #2 is significant only in the regime of very low Q and high obliquities. Neglecting any orbital or thermal evolution of the system, the observed temperature reversal could be reproduced by assuming an obliquity for the secondary while the primary's rotation axis is nearly aligned with the normal of the orbital plane. We note that the real heating will probably be greater since model #2 assumes synchronous rotation, which is not the case for both BDs in 2M053505 (see Table 1). The values of Q_{2} and necessary to account for the observed increase in L_{2} and may thus be further shifted towards more reasonable numbers, i.e. Q_{2} might also be higher than 10^{3.5} and the obliquity might be smaller than . Thus, for a narrow region in the  plane, model #2 yields tidal energy rates for the secondary comparable to its observed luminosity overshoot and in this region the computed temperature increase can explain the observed temperature reversal.
3.4 Tidal heating in 2M053505 with model #3
Since the only free parameter in this model is the putative fixed time lag , we show the tidal heating rates for both the primary and the secondary only as a function of in Fig. 7 with . For this range, model #3 yields energy rates and temperature rises that are compatible with the observed luminosity and temperature overshoot of the secondary. For s the heating rate for the secondary becomes comparable to the observed one, namely W. However, assuming a similar time lag for the primary, the luminosity gain of the primary BD would be significantly higher than that of the secondary, which is not compatible with the observations. The assumption of should be valid since both BDs are very similar in their structural properties, such as mass, composition, temperature, and radius.
The corresponding temperature increase is plotted in Fig. 8. It shows that the more massive BD would experience a higher temperature increase than its companion, assuming similar time lags. Since tidal heating is underway in 2M053505 and was probably similar in the past (see Sect. 2.1.1), tidal heating after model #3 would have been more important on the primary, forcing it to be even hotter than it would be without the perturbations of the secondary. The temperature difference between the primary and the secondary, which is anticipated by BD evolutionary models, would be even larger. Thus, the temperature inversion cannot be explained by tidal model #3.
3.5 Tidal heating in 2M053505 with model #4
The calculations based on model #4 yield significant heating rates in both BDs. Like in the case of models #1 and #2, the luminosity gain of the secondary at a fixed obliquity is, over the whole range, smaller than that of the primary (Fig. 9). As for model #2, the difference between and is less pronounced than in model #1. Assuming spinorbit alignment for the primary and a pronounced obliquity of the secondary, tidal heating rates of W can be reached with and .
Like model #2, #4 produces a reversal in temperature increase by means of the modified StefanBoltzmann relation in Eq. (27), due to the comparable heating rates of both BDs and the significantly smaller radius of the secondary (Fig. 10). We find a reversal in tidal heating, i.e. for any given point in  space. A temperature increase of 40 K can be reached with and . Since the equations of model #4 provide merely a lower limit due to the assumption of asymptotic nonsynchronous rotation, might also be higher than 10^{3.5} and the obliquity might be smaller than . Similar to model #2, tidal model #4 can reproduce the observed temperature reversal in a narrow region of the  parameter space.
4 Discussion
We employed several tidal models to explore the tidal heating in
2M053505. We found that, assuming similar tidal quality factors Q
and obliquities
for both BDs, the constantphaselag model #2 and the constanttimelag
model #4 yield a stronger increase in effective temperature on the
secondary mass BD than on the primary. For certain regimes of Q_{2}
and ,
the tidal energy rates in the secondary are of the correct amount to
explain the larger temperature in the smaller BD. A comparison
between our computations based on the models #1 and #2 on the one hand
and #3 and #4 on the other hand is difficult. The reference to a fixed
tidal time lag might only be reconciled with the assumption of
as done by Lev07, which is at least questionable since the assumption
of a fixed time lag is not compatible with a fixed phase lag.
Furthermore, model #3 does not invoke obliquities, which also
complicates direct comparisons of the model output.
Figure 7: Tidal heating within the primary (solid line) and secondary (dashed line) after model #3. While the tidal heating rate of the secondary becomes comparable to its observed luminosity overshoot for s, if the same is applied to the primary, heating within the primary would lead to a larger luminosity than is observed. 

Open with DEXTER 
Figure 8: Temperature increase of the primary (solid line) and secondary (dashed line) after model #3. Contrary to what is observed, the primary would be hotter than the secondary. 

Open with DEXTER 
4.1 Constraints on the tidal dissipation function for BDs,
4.1.1 The RossiterMcLaughlin effect in 2M053505
The geometric implication of the most promising tidal models #2 and #4 is that the obliquity of the 2M053505 primary is negligible and that of the secondary is  provided tidal heating accounts for the reversal and the luminosity excess of the secondary. There does exist an observational method to measure the geometric configuration of eclipsing systems, called the RossiterMcLaughlin effect (RME) (McLaughlin 1924; Rossiter 1924).
Figure 9: Tidal heating after model #4. Left: (Primary) Projection of onto the  plane. The stepsize between contour lines is chosen to be in . Right: (Secondary) Projection of onto the  plane. 

Open with DEXTER 
Figure 10: Temperature increase after model #4. Left: (Primary) Projection of dT_{1} onto the plane. Right: (Secondary) Projection of dT_{2} onto the plane. For and an obliquity of the temperature increase of the secondary becomes similar to the observed one. For the whole range of Q and there is an inversion in temperature increase, similar to model #2: . 

Open with DEXTER 
The RME appears during transits in front of rotating stars. Hiding a fraction of the star's surface results in the absence of some corresponding rotational velocity contribution to the broadening of the stellar lines. Thus, the changes in the line profiles become asymmetric (except for the midpoint of the transit) and the center of a certain stellar line is shifted during a transit, which induces a change of the star's radial velocity. The shape of the resulting radial velocity curve depends on the effective area covered by the transiting object and its projected path over the stellar surface with respect to the spin axis of the covered object (for a detailed analysis of the RME see Ohta et al. 2005).
Using a code originally presented in Dreizler et al. (2009), we have undertaken simulations of the RME for various geometric configurations of 2M053505 during the primary eclipse^{} as it would be seen with the Ultraviolet and Visual Echelle Spectrograph (UVES) at the Very Large Telescope (VLT) (see Fig. 11). For the data quality we assumed the constraints given by the UVES at the VLT exposure time calculator^{} in version 3.2.2. The computations show that, using ThAr reference spectra and also the telluric A and B bands as benchmarks, a time sampling with one spectrum every 1245 s and a S/N of 7 around 8600 are necessary to get 21 measurements during the primary eclipse and an accuracy of 100 m/s.
Figure 11: Simulations for the RossiterMcLaughlin effect as it would be seen with UVES during the primary eclipse of 2M053505, which occurs when the secondary mass BD is occulted by the primary. The S/N is 7. Left: the orbital inclination i is fixed at (see Table 1) and , which means the transiting primary BD follows a path parallel to the secondary's equator. The alignment of the secondary's spin axis varies between (perpendicular to the line of sight) and . Right: with i fixed at and , varies between (primary path parallel to the secondary's equator) and (primary path strongly misaligned with the secondary's equator). 

Open with DEXTER 
In principle, there are four parameters for the background object of the transit to be fitted in our simulations of the RME: the rotational velocity , the inclination of the spin axes with respect to the line of sight , the angle between the projection of the spin and the projection of the orbital plane normal onto the celestial plane , and the orbital inclination with respect to the line of sight i. From light curve analyses, both rotational velocities in 2M053505 and the orbital inclination i are known. Thus, for the simulation of the primary eclipse and are the remaining free parameters.
The obliquities ,
i.e. the real 3dimensional angle between the orbital normal and the
spin axis of the occulted object, is related to the other
angles as
While the two obliquities are intrinsic angles of the system, they cannot be measured directly. They can only be inferred from i, and , which depend on the position of the observer with respect to the system. Since we are only interested in the possible options for the measurement of the obliquities in 2M053505, we refer the reader to the paper by Winn et al. (2005) for a discussion of Eq. (28) and the geometrical aspects of the RME. With the first term in Eq. (28) degrades to insignificance, which yields .
At low values for and the fitted solutions to the RME are degenerate and there are multiple solutions within a certain confidence interval. But our simulations for the transit show that the error due to the observational noise is on the same order as the error due to degeneracy and thus we find standard deviations in and of and , respectively. The uncertainty in depends not only on the uncertainties in and but also on the actual values of and . But in all cases, the standard deviation in the secondary's obliquity .
If present in 2M053505, a considerable misalignment of the secondary BD of could be detected with a 1 accuracy of or less. Thus, an observed value of would be a 2.5 detection of spinorbit misalignment. Unless RME measurements suggest , RME observations alone are unlikely to provide definitive evidence that any of the tidal models we consider is responsible for the temperature reversal.
4.1.2 Further observations of BD binaries
Besides the option of RME measurements for testing the geometric implications, there does exist a possibility to verify our estimate of for BDs in general. Comparison of observed orbital properties with values constrained by the equations that govern the orbital evolution might constrain the free parameters, here Q. Using Eq. (6), we find that, assuming only a slight initial eccentricity of 0.05, the eccentricity of a BD binary system similar to 2M053505, in terms of masses, radii, rotational frequencies, and semimajor axis would increase to 1 after 500 Myr if the quality factors of the two BDs are 10^{3.5} (see left panel in Fig. 12). A measurement of e in such an evolved state could not constrain Q in a 2M053505 analog since either the initial eccentricity could have been relatively large while the orbit evolved rather slowly due to high Q values or a small initial value of e could have developed to a large eccentricity due to small values of Q.
Figure 12: Orbital evolution of a 2M053505 analog after model #1. Left: eccentricity evolution for different values of for the next 500 Myr. The initial eccentricity was arbitrarily chosen: . For this binary will be disrupted within 500 Myr. Right: eccentricity evolution of a 2M053505 analog but with P_{1} = P_{2} = 14.05 d for different values of . Contrary to the scenario in the left figure, the changed rotational period of the primary BD now leads to circularization of the system. Measurements of e in LMS binaries with known ages can give lower limits to . 

Open with DEXTER 
We also simulate the evolution of a 2M053505 analog but with a different rotational frequency of the primary constituent in order to let the eccentricity decrease with time. We neglected the evolution of all the other physical and orbital parameters since we are merely interested in a tentative estimate. For a given candidate system the analysis would require a selfconsistent coupled evolution of all the differential equations. For the arbitrary case of P_{1} = P_{2} = 14.05 d we find that, even for the most extreme but unrealistic case of an initial eccentricity equal to 1, this fictitious binary would be circularized on a timescale of 100 Myr for (see right panel in Fig. 12). Findings of old, eccentric BD binaries with rotational and orbital frequencies that yield circularization in the respective system would set lower limits to Q.
4.1.3 Rotational periods in 2M053505
Another, and in fact a crucial, constraint on Q for BDs comes from the synchronization time scale of the two BDs in 2M053505. Following the equation given in Lev07 and taking the initial orbital mean motion and semimajor axis of the system as calculated with an uncoupled system of differential equations from model #1, we derive Myr for the primary and Myr for the secondary with . Since the rotation in both BDs is not yet synchronized with the orbit and the age of the system is about 1 Myr, is not consistent with the age of 2M053505. Both components should have synchronous rotation rates already. We find the consistent value for Q to be 10^{4.5}, yielding synchronization time scales Myr and Myr.
To make this estimate for Q more robust, we present the evolution of the BDs' rotational periods in Fig. 13 and compare it to the critical period for a structural breakup . The evolutionary tracks are calculated with model #1 and Eq. (30) in FM08. As a rough approach we do not couple this equation with those for the other orbital parameters. The left panel of Fig. 13 shows that for and the primary's initial rotation period 1 Myr ago is 0.3 d. The initial rotation period for the secondary, for and , is about 0.2 d, where the algebraic sign contributes for a retrograde revolution (right panel in Fig. 13). For most of its lifetime, the secondary would have had a retrograde rotation and just switched the rotation direction within the last few 10 000 yr, which is very unlikely in statistical terms. Since the orbital momentum is on the order of 10^{43} kg m^{2} / s and the individual angular momenta are about 10^{41} kg m^{2} / s, the shrinking process might not have had a serious impact on the rotational evolution. Tides have dominated the spin evolutions.
Figure 13: Rotational evolution of the two BDs in 2M053505 after model #1 for different values of Q_{1} and Q_{2}. Left: (Primary) Going backwards in time, the rotation period decreases. For , P_{1} drops below the critical period for structural breakup of 0.5 d already before the date of birth around 1 Myr ago. Right: (Secondary) For we show the tracks for and for comparison. For the rotation direction switches at about 0.18 Myr and for at roughly 10 000 yr. 

Open with DEXTER 
Following Scholz & Eislöffel (2005), the critical breakup period depends only on the body's radius and its mass. The radius evolution for BDs is very uncertain for the first Myr after formation but we estimate their initial radii to be as large as the solar radius. This yields d for both the primary and the secondary BD. As stated above, the moduli of the initial rotation periods of both BDs would have been smaller than 0.5 d for Q values of 10^{3.5}. This inconsistency gives a lower limit to Q_{1} and Q_{2} since values of 10^{3.5} would need an initial rotation periods of both BDs which are smaller than their critical breakup periods. Obliquities larger than would accelerate the (backwards) evolution and yield even larger lower limits for Q_{1} and Q_{2}. Thus, our simulations of the rotational period evolution of both BDs require , whereas the tidal synchronization timescale even claims .
4.2 Evolutionary embedment of tidal heating
Tidal heating must be seen in the evolutionary context of the system. On the one hand, the tidal energy rates generate a temperature increase on the KelvinHelmholtz timescale, which is 2 Myr for the BDs in 2M053505  and thus on the order of the system's age, as per Eq. (27). On the other hand, tidal heating will affect the shrinking and cooling process of young BDs in terms of an evolutionary retardation. As models show (D'Antona & Mazzitelli 1997; Chabrier et al. 2000; Baraffe et al. 1998; Chabrier & Baraffe 2000, single BDs cool and shrink significantly during their first Myrs after formation. Adding an energy source comparable to the luminosity of the object will slow down the aging processes such that the observed temperature and luminosity overshoot at some later point is not only due to the immediate tidal heating but also due to its past evolution. Consequently, the luminosity and temperature overshoot in the secondary might not (only) be due to presentday tidal heating, but it could be a result of an evolutionary retardation process triggered by the presence of the primary as a perturber. Coupled radiusorbit evolutionary models have already given plausible explanations for the inflated radii of some extrasolar planets (Miller et al. 2009; Ibgui & Burrows 2009; Ibgui et al. 2009; Gu et al. 2003; Ibgui et al. 2010).
For a consistent description of the orbital and physical history of 2M053505, one would have to include the evolution of obliquities , BD radii R_{i}, eccentricity e, semimajor axis a, and rotational frequencies . Note that there is a positive feedback between radial inflation and tidal heating: as tidal heating inflates the radius, the tidal heating rate can increase and  in turn  may cause the radius to inflate even more. In a selfconsistent orbital and structural simulation of 2M053505, tidal inflation, neglected in our computations of the increase in Eq. (27), will result naturally from the additional heating term introduced by tides.
In conjunction with 2M053505 that means the actual heating rates necessary to explain the and luminosity excess in the secondary are lower than they would have to be if there would be no historical context. Relating to Figs. 4, 6, 8, and 10, the implied obliquity and Q factor for the secondary are  again  shifted towards lower and higher values, respectively. Embedded in the historical context of tidal interaction in 2M053505, and may also explain the temperature reversal and the luminosity excess of the secondary.
These trends, however, are contrary to that induced by tidal inflation. If tidal heating is responsible for a radial expansion of 10 and 20% in the primary and secondary, the values of the dissipation factor necessary to explain the reversal would be 0.8 smaller in (see Sect. 2.3).
5 Conclusions
We surveyed four different published tidal models, but neglect any evolutionary background of the system's orbits and the components' radii to calculate the tidal heating in 2M053505. Our calculations based on models #2 and #4, which are most compatible with the observed properties of the system, require obliquities , and a quality factor in order to explain the luminosity excess of the secondary. Additionally, the observed temperature reversal follows naturally since we may reproduce a reversal in temperature increase due to tides: . In model #2, synchronous rotation of the perturbed body is assumed. Since this is not given in 2M053505, the actual heating rates will be even higher than those computed here. Our results for the heating rates as per model #2 are thus lower limits, which shifts the implied obliquity of the secondary and its Q factor to lower and higher values, respectively.
Considerations of the synchronization time scale for the BD duet and the individual rotational breakup periods yield constraints on for BDs. We derive a lower limit of . This is consistent with estimates of Qvalues for M dwarfs, , and the quality factors of Jupiter, , and Neptune, (see Sect. 2.1.1). With tidal heating alone can neither explain the temperature reversal in the system nor the luminosity excess of the secondary.
An obliquity of , however, would be reasonable in view of recent results from measurements of the RME in several transiting exoplanet systems^{}. Currently, out of 18 planets there are 7 with significant spinorbit misalignments and some of them are even in retrograde orbits around their host stars. A substantial obliquity might cause an enhanced heating in the 2M053505 secondary, while the primary's spin could be aligned with the orbital spin, leading to negligible heating in the primary.
Despite the advantages of distanceindependent radius and luminosity measurements of close, lowmass binaries, the comparison of fundamental properties of the constituents with theoretical models of isolated BDs must be taken with care. This applies also to the direct translation from the discrepancies between observed and modeled radii for a fixed metallicity into an apparent age difference as a calibration of LMS models (Stassun et al. 2009). Tidal heating might be a crucial contribution to discrepancies between predicted and observed radii in other eclipsing lowmass binary systems (Ribas et al. 2008). As recently shown by Ibgui & Burrows (2009), tidal heating in extrasolar giant planets in close orbits at AU with modest to high eccentricities of can explain the increased radii of some planets, when embedded in the orbital history with its host star.
Improvement of tidal theories is necessary to estimate the relation between tides and the observed radii of LMS being usually too large as compared to models. A tidal model is needed for higher orders of arbitrary obliquities and eccentricities that also accounts for arbitrary rotation rates. As stated by Greenberg (2009), a formal extension of the simple ``lagandadd'' procedure of tidal frequencies the theory of constant phase lag is questionable. Besides the extension, conciliation among the various models is needed. The results from the models applied here should be considered preliminary but are suggestive and indicate the possible importance of tides in binary BD systems.
Several issues remain to be addressed for a more detailed assessment of tidal heating in 2M053505: i) reconciliation and improvement of tidal theories; ii) selfconsistent simulations of the orbital and physical evolution of the system and the BDs; iii) measurements of the system's geometric configuration; iv.) constraints on the tidal quality factors of BDs.
Acknowledgements
Our sincere thanks go to J. L. Bean for initiating this collaboration. The advice of S. Dreizler on the computations of the RM effect and inspirations from A. Reiners on 2M053505 and BDs in general have been a valuable stimulation to this study. We acknowledge the help of Y. G. M. Chew and K. G. Stassun on the parametrization of the 2M053505 BD binary and we appreciate the contribution from Baraffe to the modeling of the BDs' structures. The referee JeanPaul Zahn deserves our honest gratitude for his crucial remark on the tidal synchronization time scale. R. Heller is supported by a PhD scholarship of the DFG Graduiertenkolleg 1351 ``Extrasolar Planets and their Host Stars''. R. Barnes acknowledges funding from NASA Astrobiology Institute's Virtual Planetary Laboratory lead team, supported by NASA under Cooperative Agreement No. NNH05ZDA001C. R. Greenberg, B. Jackson, and R. Barnes were also supported by a grant from NASA's Planetary Geology and Geophysics program. This research has made use of NASA's Astrophysics Data System Bibliographic Services.
References
 Aksnes, K., & Franklin, F. A. 2001, AJ, 122, 2734 [NASA ADS] [CrossRef] [Google Scholar]
 Baraffe, I., Chabrier, G., Allard, F., et al. 1998, A&A, 337, 403 [NASA ADS] [Google Scholar]
 Baraffe, I., Chabrier, G., Allard, F., et al. 2002, A&A, 382, 563 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barnes, R., Jackson, B., Raymond, S. N., West, A. A., & Greenberg, R. 2009, ApJ, 695, 1006 [NASA ADS] [CrossRef] [Google Scholar]
 Bodenheimer, P., Lin, D. N. C., & Mardling, R. A. 2001, ApJ, 548, 466 [NASA ADS] [CrossRef] [Google Scholar]
 Boltzmann, L. 1884, Annalen der Physik, 258, 291 [NASA ADS] [CrossRef] [Google Scholar]
 Brooker, R. A., & Olle, T. W. 1955, MNRAS, 115, 101 [NASA ADS] [Google Scholar]
 Çakirli, Ö., Ibanoglu, C., & Güngör, C. 2009, New Astron., 14, 496 [NASA ADS] [CrossRef] [Google Scholar]
 Chabrier, G., & Baraffe, I. 2000, ARA&A, 38, 337 [NASA ADS] [CrossRef] [Google Scholar]
 Chabrier, G., Baraffe, I., Allard, F., et al. 2000, ApJ, 542, 464 [NASA ADS] [CrossRef] [Google Scholar]
 Chabrier, G., Gallardo, J., & Baraffe, I. 2007, A&A, 472, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Coughlin, J. L., & Shaw, J. S. 2007, Journal of the Southeastern Association for Research in Astronomy, 1, 7 [NASA ADS] [Google Scholar]
 D'Antona, F., & Mazzitelli, I. 1997, Mem. Soc. Astron. Ital., 68, 807 [NASA ADS] [Google Scholar]
 Dreizler, S., Reiners, A., Homeier, D., et al. 2009, A&A, 499, 615 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 FerrazMello, S., Rodríguez, A., & Hussmann, H. 2008, Celest. Mechan. Dyn. Astron., 101, 171 (FM08) [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Gavrilov, S. V., & Zharkov, V. N. 1977, Icarus, 32, 443 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., & Soter, S. 1966, Icarus, 5, 375 [NASA ADS] [CrossRef] [Google Scholar]
 Gómez Maqueo Chew, Y., Stassun, K. G., Prsa, A., et al. 2009, ApJ, 699, 1196 [NASA ADS] [CrossRef] [Google Scholar]
 Greenberg, R. 2009, ApJ, 698, L42 [NASA ADS] [CrossRef] [Google Scholar]
 Greenberg, R., Barnes, R., & Jackson, B. 2008, in BAAS, 40, 391 [Google Scholar]
 Gu, P., Lin, D. N. C., & Bodenheimer, P. H. 2003, ApJ, 588, 509 [NASA ADS] [CrossRef] [Google Scholar]
 Guenther, E. W., Torres, G., Batalha, N., et al. 2001, A&A, 366, 965 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hut, P. 1981, A&A, 99, 126 (Hut81) [NASA ADS] [Google Scholar]
 Ibgui, L., & Burrows, A. 2009, ApJ, 700, 1921 [NASA ADS] [CrossRef] [Google Scholar]
 Ibgui, L., Spiegel, D. S., & Burrows, A. 2009, ApJ, submitted[arXiv:0910.5928] [Google Scholar]
 Ibgui, L., Burrows, A., & Spiegel, D. S. 2010, ApJ, 713, 751 [NASA ADS] [CrossRef] [Google Scholar]
 Ioannou, P. J., & Lindzen, R. S. 1993, ApJ, 406, 266 [NASA ADS] [CrossRef] [Google Scholar]
 Jackson, B., Barnes, R., & Greenberg, R. 2008a, MNRAS, 391, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Jackson, B., Greenberg, R., & Barnes, R. 2008b, ApJ, 681, 1631 [NASA ADS] [CrossRef] [Google Scholar]
 Kippenhahn, R., & Weigert, A. 1990, Stellar Structure and Evolution (Stellar Structure and Evolution, XVI (Berlin, Heidelberg, New York: SpringerVerlag), Also Astronomy and Astrophysics Library [Google Scholar]
 Levrard, B., Correia, A. C. M., Chabrier, G., et al. 2007, A&A, 462, L5 (Lev07) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mardling, R. A., & Lin, D. N. C. 2002, ApJ, 573, 829 [NASA ADS] [CrossRef] [Google Scholar]
 Mardling, R. A., & Lin, D. N. C. 2004, ApJ, 614, 955 [NASA ADS] [CrossRef] [Google Scholar]
 Marley, M. S., Fortney, J. J., Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2007, ApJ, 655, 541 [NASA ADS] [CrossRef] [Google Scholar]
 McLaughlin, D. B. 1924, ApJ, 60, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Miller, N., Fortney, J. J., & Jackson, B. 2009, ApJ, 702, 1413 [NASA ADS] [CrossRef] [Google Scholar]
 Mohanty, S., Baraffe, I., & Chabrier, G. 2007, in IAU Symp., 239, ed. F. Kupka, I. Roxburgh, & K. Chan, 197 [Google Scholar]
 Mohanty, S., Stassun, K. G., & Mathieu, R. D. 2009, ApJ, 697, 713 [NASA ADS] [CrossRef] [Google Scholar]
 Morales, J. C., Ribas, I., Jordi, C., et al. 2009, ApJ, 691, 1400 [NASA ADS] [CrossRef] [Google Scholar]
 Neron de Surgy, O., & Laskar, J. 1997, A&A, 318, 975 [NASA ADS] [Google Scholar]
 Ogilvie, G. I., & Lin, D. N. C. 2004, ApJ, 610, 477 [NASA ADS] [CrossRef] [Google Scholar]
 Ohta, Y., Taruya, A., & Suto, Y. 2005, ApJ, 622, 1118 [NASA ADS] [CrossRef] [Google Scholar]
 Peale, S. J., & Greenberg, R. J. 1980, in Lunar and Planetary Institute Conference Abstracts, 11, 871 [Google Scholar]
 Ray, R. D., Eanes, R. J., & Lemoine, F. G. 2001, Geophys. J. Int., 144, 471 [NASA ADS] [CrossRef] [Google Scholar]
 Reiners, A., Seifahrt, A., Stassun, K. G., Melo, C., & Mathieu, R. D. 2007, ApJ, 671, L149 [NASA ADS] [CrossRef] [Google Scholar]
 Ribas, I., Morales, J. C., Jordi, C., et al. 2008, Mem. Soc. Astron. It., 79, 562 [NASA ADS] [Google Scholar]
 Rossiter, R. A. 1924, ApJ, 60, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Scholz, A., & Eislöffel, J. 2005, A&A, 429, 1007 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stassun, K. G., Mathieu, R. D., & Valenti, J. A. 2006, Nature, 440, 311 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Stassun, K. G., Mathieu, R. D., & Valenti, J. A. 2007, ApJ, 664, 1154 [NASA ADS] [CrossRef] [Google Scholar]
 Stassun, K. G., Mathieu, R. D., Cargile, P. A., et al. 2008, Nature, 453, 1079 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Stassun, K. G., Hebb, L., LópezMorales, M., et al. 2009, in IAU Symp. 258, ed. E. E. Mamajek, D. R. Soderblom, & R. F. G. Wyse, 161 [Google Scholar]
 Stefan, J. 1879, Sitzungsberichte der mathematischnaturwissenschaftlichen Classe der kaiserlichen Akademie der Wissenschaften, 79, 391 [Google Scholar]
 Winn, J. N., Noyes, R. W., Holman, M. J., et al. 2005, ApJ, 631, 1215 [NASA ADS] [CrossRef] [Google Scholar]
 Wisdom, J. 2008, Icarus, 193, 637 (Wis08) [NASA ADS] [CrossRef] [Google Scholar]
 Wuchterl, G. 2005, Astron. Nachr., 326, 905 [NASA ADS] [CrossRef] [Google Scholar]
 Yoder, C. F. 1979, Nature, 279, 767 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, K., & Hamilton, D. P. 2008, Icarus, 193, 267 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
 ... heat^{}
 In contrast to the Baraffe et al. (1998) tracks, the models by D'Antona & Mazzitelli (1997) predict a temperature increase in BDs for the first 30 Myr of their existence.
 ... bulge^{}
 If e = 0 and , then there is a single tidal lag angle and the tidal dissipation funtion can be written as . For the course of an orbit, where the tidal evolution of n is negligible, both Q and can be fixed. However, in a general case where is constant in time, Q will decrease as the orbital semimajor axis decays and n increases. So Q would not be constant.
 ... rotation^{}
 Wis08 calls this ``asymptotic nonsynchronous rotation''.
 ... eclipse^{}
 The ``primary eclipse'' refers to the major flux decrease in the system's light curve. Due to the significantly higher effective temperature of the secondary mass BD the primary eclipse occurs when the primary mass component transits in front of the secondary companion, as seen from Earth.
 ... calculator^{}
 http://www.eso.org/observing/etcwww.eso.org/observing/etc
 ... systems^{}
 See www.hs.unihamburg.de/EN/Ins/Per/Heller for an overview.
All Tables
Table 1: Orbital and physical parameters of 2M053505.
All Figures
Figure 1: Degeneracy parameter (solid line) with model parameters similar to those of the 2M053505 primary and radiusintegrated luminosity L (dashed line) as a function of radius. To fit into the plot, L is normalized to 10. 

Open with DEXTER  
In the text 
Figure 2: Orbital evolution of 2M053505 after model #1 going back in time for 1.5 Myr. Left: eccentricity evolution. Depending on and on the age of the system, its initial eccentricity has not been smaller than 0.3133, which is 97.4% of its current value. Right: obliquity evolution of the secondary BD for three different values of Q_{2}. Simulations started at ``time = 0'' for and were evolved backwards in time. For there is no significant change in . For all the treated values of Q_{2}, the obliquity of the 2 M053505 secondary is still close to its natal state. 

Open with DEXTER  
In the text 
Figure 3: Tidal heating after model #1. Left: (Primary) Projection of onto the  plane. The stepsize between contour lines is chosen to be in . Right: (Secondary) Projection of onto the  plane. Although there is a dependence on , the energy rates at a fixed value for the quality factor are smaller than those for the primary. 

Open with DEXTER  
In the text 
Figure 4: Temperature increase after model #1. Left: (Primary) Projection of dT_{1} onto the  plane. For a significant temperature increase, Q_{1} would have to be much smaller than 10^{3.5}, but such a temperature increase is not observed in the primary. Right: (Secondary) Projection of dT_{2} onto the  plane. Even for very low values of Q_{2} and high obliquities the observed temperature increase cannot be reconstructed. For any given point in the  plane, , which does not support the observed temperature reversal. 

Open with DEXTER  
In the text 
Figure 5: Tidal heating after model #2. Left: (Primary) Projection of onto the  plane. The stepsize between contour lines is chosen to be in . The tidal energy rates strongly depend on a putative obliquity, different from model #1. Right: (Secondary) Projection of onto the  plane. For the three models (#1, #2, and #4) invoking Q and , these rates are the highest of all for any given point in the  plane  for the primary as well as for the secondary. 

Open with DEXTER  
In the text 
Figure 6: Temperature increase after model #2. Left: (Primary) Projection of dT_{1} onto the  plane. Right: (Secondary) Projection of dT_{2} onto the  plane. For any given location in the  plane, model #2 yields the strongest temperature increase compared to the other models  both for the primary and the secondary, respectively. For a given spot in Q space there is an inversion in temperature increase: , i.e. the less massive BD is heated more. 

Open with DEXTER  
In the text 
Figure 7: Tidal heating within the primary (solid line) and secondary (dashed line) after model #3. While the tidal heating rate of the secondary becomes comparable to its observed luminosity overshoot for s, if the same is applied to the primary, heating within the primary would lead to a larger luminosity than is observed. 

Open with DEXTER  
In the text 
Figure 8: Temperature increase of the primary (solid line) and secondary (dashed line) after model #3. Contrary to what is observed, the primary would be hotter than the secondary. 

Open with DEXTER  
In the text 
Figure 9: Tidal heating after model #4. Left: (Primary) Projection of onto the  plane. The stepsize between contour lines is chosen to be in . Right: (Secondary) Projection of onto the  plane. 

Open with DEXTER  
In the text 
Figure 10: Temperature increase after model #4. Left: (Primary) Projection of dT_{1} onto the plane. Right: (Secondary) Projection of dT_{2} onto the plane. For and an obliquity of the temperature increase of the secondary becomes similar to the observed one. For the whole range of Q and there is an inversion in temperature increase, similar to model #2: . 

Open with DEXTER  
In the text 
Figure 11: Simulations for the RossiterMcLaughlin effect as it would be seen with UVES during the primary eclipse of 2M053505, which occurs when the secondary mass BD is occulted by the primary. The S/N is 7. Left: the orbital inclination i is fixed at (see Table 1) and , which means the transiting primary BD follows a path parallel to the secondary's equator. The alignment of the secondary's spin axis varies between (perpendicular to the line of sight) and . Right: with i fixed at and , varies between (primary path parallel to the secondary's equator) and (primary path strongly misaligned with the secondary's equator). 

Open with DEXTER  
In the text 
Figure 12: Orbital evolution of a 2M053505 analog after model #1. Left: eccentricity evolution for different values of for the next 500 Myr. The initial eccentricity was arbitrarily chosen: . For this binary will be disrupted within 500 Myr. Right: eccentricity evolution of a 2M053505 analog but with P_{1} = P_{2} = 14.05 d for different values of . Contrary to the scenario in the left figure, the changed rotational period of the primary BD now leads to circularization of the system. Measurements of e in LMS binaries with known ages can give lower limits to . 

Open with DEXTER  
In the text 
Figure 13: Rotational evolution of the two BDs in 2M053505 after model #1 for different values of Q_{1} and Q_{2}. Left: (Primary) Going backwards in time, the rotation period decreases. For , P_{1} drops below the critical period for structural breakup of 0.5 d already before the date of birth around 1 Myr ago. Right: (Secondary) For we show the tracks for and for comparison. For the rotation direction switches at about 0.18 Myr and for at roughly 10 000 yr. 

Open with DEXTER  
In the text 
Copyright ESO 2010