DI Herculis as a test of internal stellar structure and general relativity
New apsidal motion rate and evolutionary models
A. Claret^{1}  G. Torres^{2}  M. Wolf^{3}
1  Instituto de Astrofísica de Andalucía, CSIC, Apartado 3004, 18080 Granada, Spain
2 
HarvardSmithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA
3  Astronomical Institute, Faculty of Mathematics and Physics,
Charles University Prague, 180 00 Praha 8, V Holesovickách 2,
Czech Republic
Received 22 December 2009 / Accepted 9 February 2010
Abstract
Context. For the past three decades there has been a serious
discrepancy between the observed and theoretical values of the apsidal
motion rate
of
the eccentric eclipsing binary DI Her, which has even been
interpreted occasionally as a possible failure of general relativity
(GR). A number of plausible explanations have been put forward.
Recent observations of the RossiterMcLaughlin effect have shown
convincingly that the reason for the anomaly is that the rotational
axes of the stars and the orbital axis are misaligned, which changes
the predicted rate of precession significantly.
Aims. Although the disagreement is now drastically smaller as a
result of those measurements, it formally remains at the level
of 50%, owing possibly to errors in the measured apsidal motion
rate, outdated stellar models, or inaccuracies in the stellar
parameters. The aim of this paper is to address each of these issues to
improve the agreement even more.
Methods. New times of minimum have been collected and used for
redetermining of the apsidal motion rate. Based on the latest
determinations of the absolute dimensions of the binary, we computed
new stellar evolution models with updated physical inputs, and derived
improved apsidal motion constants for the components. We performed
Monte Carlo simulations to infer the theoretical distribution
of
,
including the contributions from GR, as well as tidal and
rotational distortions. All observational errors have been
accounted for.
Results. Our simulations yield a retrograde apsidal motion rate
due to the rotationallyinduced oblateness
of 0.00056 deg cycle^{1} (mode of the distribution), a GR contribution of +0.00068 deg cycle^{1}, and a tidal contribution of +0.00034 deg cycle^{1}, leading to a total predicted rate of +0.00046 deg cycle^{1}. This is in excellent agreement with the newly measured value of +0.00042 deg cycle^{1}. The formal difference is now reduced to 10%, a small fraction of the observational uncertainties.
Key words: binaries: eclipsing  stars: interiors  stars: rotation
1 Introduction
For some three decades, the eccentric Btype eclipsing binary DI Her (P = 10.55 days) has been known to have an apsidal motion rate that is too slow compared to theoretical predictions (for a brief review, see Claret 1998). Until very recently, this has presented a serious challenge to our understanding of stellar physics, and even the validity of general relativity (GR) has been called into question given that the GR contribution is dominant for DI Her. The discrepancy in the total rate of precession amounts to 400% or more. This anomaly has persisted even when considering hypothetical stellar models with infinite mass concentration (apsidal motion constant k_{2} = 0), which tends to reduce the predictions, and this was often interpreted as a failure of GR. The problem was all the more puzzling when considering that studies of other binary systems suitable for testing GR and with welldetermined absolute dimensions indicated good agreement with theory (see, e.g., Claret 1997; Wolf et al. 2010). This suggested to many that the problem may lie with DI Her itself. A number of hypotheses were put forward over the years to explain the disagreement, including
 1.
 a rapid ongoing circularization of the orbit, affecting the measured times of minimum;
 2.
 the presence of a circumstellar cloud between the components, which would affect the gravitational field of the stars;
 3.
 the presence of a distant third star in the system, which would contribute to the motion of the line of apsides;
 4.
 an alternative theory of gravitation.
An attractive alternative is a configuration for the binary in which the spin axes are tilted relative to the orbital axis (Shakura 1985; Guinan & Maloney 1985; Company et al. 1988; Claret 1998). This would introduce a correction to the rotational term of the apsidal motion that can be negative (regression of the line of apsides), resulting in a smaller total . The young age of the system (a few Myr) certainly makes this plausible, as tidal forces would not have had time to operate and bring the spin axes of the stars into alignment with the axis of the orbit (see Sect. 3). A remarkable confirmation of this prediction has been obtained recently by Albrecht et al. (2009), who provided compelling observational evidence for a peculiar orientation of the stars by taking advantage of the RossiterMcLaughlin effect. This is a distortion of the spectral lines that occurs during the eclipses, as a result of the partial obscuration of the rotating stellar disks. Using this effect, they were able to measure the angles between the sky projections of the spin axes and the orbital axis, and show that there is indeed a strong misalignment for both components of the binary. It is the first documented case in which this has been shown convincingly for stars in a doublelined eclipsing binary^{}.
With this new information Albrecht et al. (2009) managed to reduce the discrepancy with the predicted apsidal motion rate of DI Her very significantly, so that theory and observation are now in much better agreement. Some of the angles involved in computing the effect of the tilted axes remain unknown, however, so the theoretical prediction necessarily assumes a statistical character in this case, which Albrecht et al. (2009) tackled by means of Monte Carlo simulations. The nominal difference between the measured apsidal motion rate and the mode of their theoretical probability distribution for is still about 50% at present, although the distribution is rather wide so theory and observation are formally consistent with each other. Some of this difference could come from the measurement of itself, from inaccuracies in the adopted absolute dimensions of the stars, or from the stellar evolution models used to infer the apsidal motion constants k_{2}, which are needed as input to the theoretical calculations.
In this paper we reexamine all the ingredients in order to provide the most accurate predictions and observational constraints. We incorporate new measurements of the eclipse times to derive an improved value of for DI Her that is significantly higher than the value adopted by Albrecht et al. (2009). We use the most recent determinations of the mass, radius, and other properties of the stars, and we revisit the temperature determinations. With these we provide new fits to stellar evolution models and derive new apsidal motion constants. We then repeat the comparison of the measured apsidal motion rate with theory, also examining the sensitivity of the predictions to several of the key observables.
Table 1: Observational parameters for DI Her.
2 Observational data
2.1 The absolute dimensions
The absolute masses and radii of the components of DI Her used here are adopted from the recent compilation by Torres et al. (2009a), in which these and other properties were reexamined or redetermined with a uniform methodology. The uncertainties are adopted from the same source. We summarize all observational results in Table 1, and describe some of them later. The effective temperatures reported by Torres et al. (2009a) are based on the individual B  V colors as given by Popper (1982) and on his color/temperature relation (Popper 1980). As a check, we used the Strömgren photometry obtained by Hilditch & Hill (1975) and the calibration of Crawford (1978), and derived a mean dereddened color for the binary that is essentially identical to Popper's. The color/temperature calibration by Paunzen et al. (2005), which combines several more recent temperature scales, then yields individual temperatures of 17 300 K and 15 400 K for the primary and secondary, with uncertainties of 800 K. These are the temperatures we adopt in this work, which are 300 K hotter than those listed by Torres et al. (2009a).
2.2 New apsidal motion rate
The apsidal motion rate of DI Her was determined by means of an diagram analysis. The method we followed is that described by Giménez & GarcíaPelayo (1983),
which is an iterative weighted leastsquares procedure including terms
in the eccentricity up to the fifth order. There are five adjustable
quantities to be determined: the rate of periastron advance
,
the sidereal period ,
the eccentricity e, and the longitude of periastron
at the reference epoch T_{0}. The periastron position
at any epoch E is then given by the linear equation
(1) 
and the relation between the sidereal and the anomalistic periods, and , is given by
(2) 
The value for the orbital inclination was adopted as from the lightcurve analysis of Guinan & Maloney (1985).
Table 2: New times of minimum light of DI Her.
We used all precise photoelectric and CCD timings recently recalculated by Kozyreva & Bagaev (2009, their Tables 1 and 2), as well as new eclipse times obtained by us at two observatories with the following instrumentation:
 Observatory and Planetarium of Johann Palisa, Ostrava, Czech Republic: 0.20m or 0.30m telescopes with an SBIG ST8XME CCD camera and Johnson V, R, and I filters;
 Valasské Mezirící Observatory, Czech Republic: 0.30m Celestron Ultima telescope with an SBIG ST7 CCD camera and Johnson R filter.
Figure 1: Differential light curve of DI Her during the primary minimum of 2005 May 20, as measured by Ms. Hana Kucáková at the Johann Palisa Observatory, Technical University Ostrava, Czech Republic. The observations were made through a Johnson R filter with a 20cm reflecting telescope and SBIG ST8XME CCD camera, with exposure times of 25 s. The comparison star was GSC 0210900247 (V = 10.04). The light curve contains more than 500 points obtained during nearly 6 h of monitoring. 

Open with DEXTER 
The solution was found by minimizing the standard sum, and the errors of the observations were adopted as published. The resulting apsidal motion parameters and their uncertainties from the least squares fitting are presented at the bottom of Table 1, and the corresponding diagram is shown in Fig. 2. The value for the best fit was . This is considerably greater than the number of degrees of freedom (), which is most likely due to the formal errors of the eclipse timings being underestimated, as is often found to be the case. The uncertainty of reported in Table 1 already accounts for this effect. The extra scatter compared to the uncertainties is seen more clearly in Fig. 3, in which the influence of the apsidal motion terms has been subtracted from the measurements. There is no evidence in this diagram of any other significant variations, such as might be induced by the presence of a third body in the system.
Figure 2: The diagram of DI Her (eclipse timing residuals), together with our bestfit apsidal motion model. The solid line and circles correspond to the primary, and the dashed line and triangles to the secondary. 

Open with DEXTER 
Figure 3: The residuals from the photoelectric and CCD timings of DI Her, after subtraction of the apsidal motion terms. Symbols are as in Fig. 2. 

Open with DEXTER 
3 Stellar evolution models
As a prerequisite for investigating the theoretical apsidal motion of DI Her, the observed stellar properties such as the mass, radius, and temperature of each component must be well matched by stellar evolution models at a single age for the binary, since these models are then used to infer the apsidal motion constants k_{2} that factor into the rotational and tidal contributions to . The basic properties of the stellar evolution code used here have been described by Claret (2004). The models include overshooting beyond the formal boundary of the convective core, as determined by the Schwarzschild criterion. The distance over which this is considered to occur is specified as = , where is the pressure scale height taken at the edge of the convective core, as given by Schwarzschild's criterion, and is a free parameter. We adopt here , following the results of Claret (2007). The evolution code incorporates radiative opacities from the tables provided by Iglesias & Rogers (1996), and the calculations of Alexander & Ferguson (1994) for lower temperatures.
The chemical composition of doublelined eclipsing binaries is often difficult to determine observationally, and that is the case also for DI Her. This introduces some uncertainty in the comparison with stellar evolution models, since the metal abundance Z then becomes a free parameter. For young systems such as DI Her, Z is quite sensitive to the effective temperatures of the stars. We computed evolutionary tracks for the measured masses and explored a range of compositions, finding the best fit for (X, Z) = (0.71, 0.02), which suggests an abundance very near that of the Sun. The observations along with our best fit are displayed in Fig. 4. The agreement with the temperatures is quite satisfactory, given the error bars and other uncertainties described below. Precise age determinations for young systems are not easy to make and require very accurate measurements of the absolute radii. This is illustrated clearly in Fig. 5. As small as the radius errors already are in DI Her (<2% relative errors, which are among the best available for eclipsing binaries; see Torres et al. 2009a), the uncertainties allow for a wide range of ages between 2 Myr and 7 Myr. Nominally the best fit to the radii is achieved at about 4.5 Myr with these models. We note, in passing, that other models such as those from the YonseiYale series (Yi et al. 2001; Demarque et al. 2004) do not yield as good a fit, at least with the grid of standard compositions (X, Z) that are publicly available.
Figure 4: Stellar evolution models compared with the measurements for DI Her. The evolutionary tracks are computed for the measured masses (see Table 1). The best fit has (X, Z) = (0.71, 0.02), implying a chemical composition near solar. 

Open with DEXTER 
Figure 5: Radius evolution of DI Her, based on the models shown in Fig. 4. The shaded areas represent the measured radii and their uncertainties. The bestfit age of about 4.5 Myr is indicated by the vertical dotted line. 

Open with DEXTER 
The theoretical apsidal motion constants inferred from our models are 0.05 and 0.05 for the primary and secondary, respectively. These values indicate that the modern models favor stars that are slightly more centrally concentrated in mass than those adopted by Claret (1998) and Albrecht et al. (2009). This goes in the direction of decreasing the predicted apsidal motion rate. Tests were made to evaluate the influence on the inferred k_{2} of uncertainties in the age determination and in the chemical composition, and these uncertainties are included in the errors reported above.
Regarding the measured effective temperatures, it is worth mentioning that because the spin axes of the stars are tilted relative to the orbital axis, the temperature distribution over their disks could be highly nonuniform and asymmetric as viewed from the Earth because of gravity darkening, and it is even possible one or both components are viewed nearly poleon. There may therefore be a bias in the measured effective temperatures that could amount to several hundred degrees for DI Her, but because the exact orientation is not known (see next section), the effect is difficult to estimate more precisely. Rotation itself changes the internal structure of stars and affects their evolution, generally leading to cooler temperatures. To explore this effect more quantitatively for DI Her we have computed models accounting for rotation (rigidbody approximation), and we find differences of 200 K, depending on the inclination relative to the line of sight. Once again, since the precise orientation of the spin axes is unknown, we are unable to be more specific regarding this effect.
With the models computed in this section we have calculated the age at which the stellar spin axes in DI Her are expected to become aligned with the axis of the orbit due to the action of tidal forces. The procedure follows closely the one described by Torres et al. (2009b, and references therein) for the case of the eclipsing system LV Her. The result using the theory by Zahn (1977, 1989) is 10^{8} yr. With the more general equations by Hut (1981), valid for systems of high eccentricity such as this, we obtain 3 10^{7} yr, which is an order of magnitude longer than the evolutionary age. These estimates suggest, as is observed, that tidal forces have not had enough time in this system to make the rotation axes parallel to the axis of the orbit.
4 The apsidal motion of DI Her
The total rate of apsidal motion from theory to be compared with the
measured value is given by the sum of the individual contributions of
each component due to tidal and rotational distortions and the general
relativistic term. When accounting for the possibility of misalignment
(Shakura 1985), we have
where
The subindices ``tidal, j'' and ``rot, j'' in Eq. (3) denote the tidal and rotational contributions of component j, and ``GR'' represents the relativistic contribution to the periastron advance. The angle i is the measured inclination of the orbital plane relative to the line of sight (see Table 1), are the angles between the rotation axes and the normal to the orbital plane, and are the angles between the rotation axes and the line of sight to the observer. The symbols and represent the rotational angular velocity for component j and the Keplerian orbital angular velocity, A is the semimajor axis of the orbit, and R_{j} and m_{j} are the radius and mass of component j. The symbols k_{2j} represent the apsidal motion constant of component j, which we determined in the previous section. The auxiliary functions f(e) and g(e) depend only on the eccentricity e and can be written as f(e) = (1  e^{2})^{2} and g(e) = (8 + 12e^{2} + e^{4}) f(e)^{2.5}/8. The factors capture the dependence on the orientation of the stars. When the spin axes are perfectly aligned with the orbit, , , and , so that Eq. (3) reduces to the classical expression presented by Sterne (1939a,b).
The relativistic contribution to the apsidal motion, which is independent of stellar structure, was given by LeviCivita (1937) as
where m_{i} and A are in solar units.
By introducing numerical values in Eq. (4), we obtain for the tidal contribution = 0.000339^{+40}_{34} deg cycle^{1}, where the error bars are given in units of the last significant digit and account for all observational uncertainties. This is about 10% less than obtained by Albrecht et al. (2009), who adopted slightly different absolute dimensions for the stars. The relativistic contribution is = 0.000677^{+6}_{7} deg cycle^{1}.
Figure 6: Monte Carlo calculations for the apsidal motion of DI Her, based on 2 10^{5} simulations. From top to bottom, the first three panels show the separate contributions from GR, tidal distortions (2nd and 3rd terms of Eq. (3)), and rotational distortions (4th and 5th terms of Eq. (3)), based on the measurements by Albrecht et al. (2009). The distributions account for all observational errors. The bottom panel shows the predicted total apsidal motion rate on the left (sum of all the above contributions), and on the right the result obtained if the spin axes are assumed to be aligned with the orbital axes, rather than strongly misaligned, as observed. The measured rate of apsidal motion of = 0.00042 0.00012 deg cycle^{1} is indicated with the vertical dashed line and shaded area, and agrees remarkably well with the calculation for misaligned axes. The mode of the distribution for the case of aligned axes, on the other hand, is more than 5 times larger than the measured value. 

Open with DEXTER 
Figure 7: Sensitivity of the Monte Carlo calculations for the predicted apsidal motion rate of DI Her to changes in some of the input quantities measured with the RossiterMcLaughlin effect. The hatched histogram is the result using the values exactly as reported by Albrecht et al. (2009), and the vertical dashed line and shaded region marks the new measured apsidal motion rate reported in this work. The shaded histogram on the left is for = 160 km s^{1} instead of the nominal value of 108 km s^{1}, with all other quantities left at their measured values. The shaded histogram on the right is for instead of the observed value of 72. Similar effects are seen when varying those parameters for the secondary. 

Open with DEXTER 
The calculation of the rotational contribution is more involved because
of to the dependence on the orientation of the system. The angles
and
in Eq. (6) are not directly measured. However, they are related through the expression
in which are accessible to observation through the RossiterMcLaughlin effect and represent the angles between the projections of the spin axes and the orbital axis on the plane of the sky for each star. Albrecht et al. (2009) determined their values to be and for the primary and secondary of DI Her, respectively. The angular velocities depend on the projected linear rotational velocities, which were measured spectroscopically by Albrecht et al. (2009) as 4 km s^{1} and 4 km s^{1}, and on the unknown inclination angles needed to deproject the measured velocities in order to infer the equatorial values. We may therefore write
in which we represent in square brackets to indicate that it is a measured quantity.
We follow Albrecht et al. (2009) and perform Monte Carlo simulations in order to determine and , and with these the total rate of apsidal motion in Eq. (3). The result is a probability distribution rather than a single value. For this we assume that the angles are distributed completely at random^{} or, equivalently, that has a uniform distribution. Knowledge of is sufficient to determine through Eq. (8), and through Eq. (9). As in Albrecht et al. (2009), we exclude angles that lead to equatorial rotational velocities exceeding the breakup velocity, for which we adopt 600 km s^{1}.
In the first three panels of Fig. 6, we show separately the GR, tidal, and rotational contributions to the apsidal motion rate. The distributions for the first two were also generated with Monte Carlo simulations in which we have assumed that all observational quantities have normally distributed errors. The distribution for the rotational term in the third panel is much wider and asymmetrical, but also sharply peaked, and has a negative modal value of = 0.00056 deg cycle^{1}, as proposed by Shakura (1985) and others, indicating a retrograde motion of the line of apsides due to the misalignment of the spin and orbital axes. This result is numerically some 10% larger than that obtained by Albrecht et al. (2009), on account of the slightly different absolute dimensions and theoretical k_{2j} values adopted. The bottom panel of Fig. 6 displays the predicted distribution for the total apsidal motion rate (left), which has a mode of = 0.000462^{+7}_{473} deg cycle^{1} (asymmetric 1 errors enclosing 68.3% of the results). The agreement with the new measured value of = 0.00042 0.00012 deg cycle^{1} from Sect. 2.2 is excellent, although the theoretical distribution has a long tail toward the left that could certainly accomodate lower values of . Still, it is interesting to note that the measurement now agrees much more closely with the most probable value from theory, because the formal difference is only 10%. The improvement compared to the results by Albrecht et al. (2009) is due mostly to the refinement of the measured apsidal motion rate, and to some extent also to the updated absolute dimensions and stellar models leading to more realistic apsidal motion constants. In particular, our models imply stars that are somewhat more centrally concentrated in mass.
If the spin axes of DI Her are assumed to be perfectly aligned with the orbit, instead of strongly misaligned as observed, the resulting Monte Carlo distribution for the total apsidal motion rate is the one shown on the righthand side of the bottom panel of Fig. 6. The mode of this distribution is = 0.00220^{+18}_{13} deg cycle^{1}, which is 5.2 times greater than the measured value. The distribution is completely inconsistent with the observation. This result for the case of aligned axes is nearly a factor of two more than that presented by Albrecht et al. (2009). The reason for this significant difference is that theirs is based on calculations by Guinan & Maloney (1985) in which the projected rotational velocities of both stars were taken to be 45 km s^{1}, as measured spectroscopically in the 1970s, whereas we used the values reported by Albrecht et al. (2009) themselves, obtained in 2008, which are considerably higher (Table 1).
The sensitivity of our Monte Carlo calculations for the misaligned case to some of the quantities measured by Albrecht et al. (2009), such as and , is illustrated in Fig. 7. We show the changes in the resulting distributions when these values are altered one at a time for the primary star, with all other observational quantities held at their nominal values.
The DI Her observations of Albrecht et al. (2009) exploiting the RossiterMcLaughlin effect represent a breakthrough in our understanding of this system, and the accuracy of their measurements has played a central role in achieving the now excellent agreement with the measured apsidal motion rate of the binary, to within 10% (a small fraction of the measurement error). The main factor leading to this improvement has been the new measurement of the apsidal motion rate presented here, followed by other refinements including the consideration of stellar models with updated input physics, and more accurate absolute dimensions for the binary components. We now know that the spin axes of the DI Her components are strongly misaligned, and the mystery of its slow apsidal motion rate is solved. However, questions remain. During the process of binary star formation the stellar spins and the orbit presumably derive their angular momentum from the same source, which is the parent molecular cloud. One would then naturally expect that the axes would be closely aligned, yet in this system they are not. Thus, interest in this object is sure to continue as theorists try to answer these questions about DI Her, and observers look for other systems like it.
AcknowledgementsThe Spanish MEC (AYA200606375) is gratefully acknowledged for its support during the development of this work. G.T. acknowledges partial support through grant AST0708229 from the US National Science Foundation. The research of M.W. was supported by the Research Program MSM0021620860 of the Ministry of Education of the Czech Republic. We would like to thank Ms. Hana Kucáková, Ostrava Observatory, and Mr. Ladislav Smelcer, Valasské Mezirící Observatory, for their valuable assistance with the photometric observations, and the anonymous referee for helpful suggestions. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France, and of NASA's Astrophysics Data System Abstract Service.
References
 Albrecht, S., Reffert, S., Snellen, I. A. G., & Winn, J. N. 2009, Nature, 461, 373 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437, 879 [NASA ADS] [CrossRef] [Google Scholar]
 Claret, A. 1997, A&A, 327, 11 [NASA ADS] [Google Scholar]
 Claret, A. 1998, A&A, 330, 533 [NASA ADS] [Google Scholar]
 Claret, A. 2004, A&A, 424, 919 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Claret, A. 2007, A&A, 475, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Company, R., Portilla, M., & Giménez, A. 1988, ApJ, 335, 962 [NASA ADS] [CrossRef] [Google Scholar]
 Crawford, D. L. 1978, AJ, 83 [Google Scholar]
 Demarque, P., Woo, J.H., Kim, Y.C., & Yi, S. K. 2004, ApJS, 155, 667 [NASA ADS] [CrossRef] [Google Scholar]
 Giménez, A., & GarcíaPelayo, J. M. 1983, Ap&SS, 92, 203 [NASA ADS] [CrossRef] [Google Scholar]
 Guinan, E. F., & Maloney, F. P. 1985, AJ, 90, 1519 [NASA ADS] [CrossRef] [Google Scholar]
 Hilditch, R. W., & Hill, G. 1975, Mem. R. Astr. Soc., 79, 101 [NASA ADS] [Google Scholar]
 Hut, P. 1981, A&A, 99, 126 [NASA ADS] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Khodykin, S. A. 2007, IBVS, 5788 [Google Scholar]
 Kozyreva, V. S., & Bagaev, L. A. 2009, Astron. Lett., 35, 483 [NASA ADS] [CrossRef] [Google Scholar]
 Kwee, K. K., & van Woerden, H. 1956, Bull. Astron. Inst. Netherlands, 12, 327 [NASA ADS] [Google Scholar]
 LeviCivita, T. 1937, Amer. J. Math., 59, 225 [CrossRef] [Google Scholar]
 Motl, D. 2007, CM UNIPACK, http://cmunipack.sourceforge.net [Google Scholar]
 Paunzen, E., Schnell, A., & Maitzen, H. M. 2005, A&A, 444, 941 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Popper, D. M. 1980, ARA&A, 18, 115 [NASA ADS] [CrossRef] [Google Scholar]
 Popper, D. M. 1982, ApJ, 254, 203 [NASA ADS] [CrossRef] [Google Scholar]
 Runkle, R. C. 2003, Ph.D. Thesis, University of North Carolina, unpublished [Google Scholar]
 Shakura, N. I. 1985, Sov. Astr. Lett., 11(4), 225 [NASA ADS] [Google Scholar]
 Sterne, T. E. 1939a, MNRAS, 99, 451 [NASA ADS] [Google Scholar]
 Sterne, T. E. 1939b, MNRAS, 99, 662 [NASA ADS] [Google Scholar]
 Torres, G., Andersen, J., & Giménez, A. 2009a, A&A Rev., 18, 67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Torres, G., Lacy, C. H. S., & Claret, A. 2009b, AJ, 138, 1622 [NASA ADS] [CrossRef] [Google Scholar]
 Wolf, M., Claret, A., Kotkova, H., et al. 2010, A&A, 509, 18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Yi, S. K., Demarque, P., Kim, Y.C., et al. 2001, ApJS, 136, 417 [NASA ADS] [CrossRef] [Google Scholar]
 Zahn, J.P. 1977, A&A, 57, 383 [NASA ADS] [Google Scholar]
 Zahn, J.P. 1989, A&A, 220, 112 [NASA ADS] [Google Scholar]
Footnotes
 ... binary^{}
 Spinorbit misalignment has been reported previously for several host stars of extrasolar transiting planets.
 ... random^{}
 In other words, Prob() .
All Tables
Table 1: Observational parameters for DI Her.
Table 2: New times of minimum light of DI Her.
All Figures
Figure 1: Differential light curve of DI Her during the primary minimum of 2005 May 20, as measured by Ms. Hana Kucáková at the Johann Palisa Observatory, Technical University Ostrava, Czech Republic. The observations were made through a Johnson R filter with a 20cm reflecting telescope and SBIG ST8XME CCD camera, with exposure times of 25 s. The comparison star was GSC 0210900247 (V = 10.04). The light curve contains more than 500 points obtained during nearly 6 h of monitoring. 

Open with DEXTER  
In the text 
Figure 2: The diagram of DI Her (eclipse timing residuals), together with our bestfit apsidal motion model. The solid line and circles correspond to the primary, and the dashed line and triangles to the secondary. 

Open with DEXTER  
In the text 
Figure 3: The residuals from the photoelectric and CCD timings of DI Her, after subtraction of the apsidal motion terms. Symbols are as in Fig. 2. 

Open with DEXTER  
In the text 
Figure 4: Stellar evolution models compared with the measurements for DI Her. The evolutionary tracks are computed for the measured masses (see Table 1). The best fit has (X, Z) = (0.71, 0.02), implying a chemical composition near solar. 

Open with DEXTER  
In the text 
Figure 5: Radius evolution of DI Her, based on the models shown in Fig. 4. The shaded areas represent the measured radii and their uncertainties. The bestfit age of about 4.5 Myr is indicated by the vertical dotted line. 

Open with DEXTER  
In the text 
Figure 6: Monte Carlo calculations for the apsidal motion of DI Her, based on 2 10^{5} simulations. From top to bottom, the first three panels show the separate contributions from GR, tidal distortions (2nd and 3rd terms of Eq. (3)), and rotational distortions (4th and 5th terms of Eq. (3)), based on the measurements by Albrecht et al. (2009). The distributions account for all observational errors. The bottom panel shows the predicted total apsidal motion rate on the left (sum of all the above contributions), and on the right the result obtained if the spin axes are assumed to be aligned with the orbital axes, rather than strongly misaligned, as observed. The measured rate of apsidal motion of = 0.00042 0.00012 deg cycle^{1} is indicated with the vertical dashed line and shaded area, and agrees remarkably well with the calculation for misaligned axes. The mode of the distribution for the case of aligned axes, on the other hand, is more than 5 times larger than the measured value. 

Open with DEXTER  
In the text 
Figure 7: Sensitivity of the Monte Carlo calculations for the predicted apsidal motion rate of DI Her to changes in some of the input quantities measured with the RossiterMcLaughlin effect. The hatched histogram is the result using the values exactly as reported by Albrecht et al. (2009), and the vertical dashed line and shaded region marks the new measured apsidal motion rate reported in this work. The shaded histogram on the left is for = 160 km s^{1} instead of the nominal value of 108 km s^{1}, with all other quantities left at their measured values. The shaded histogram on the right is for instead of the observed value of 72. Similar effects are seen when varying those parameters for the secondary. 

Open with DEXTER  
In the text 
Copyright ESO 2010