{'fraction_non_alphanumeric': 0.039028030531153036, 'fraction_numerical': 0.030622798863591733, 'mean_word_length': 4.507390583315307, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 3, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Absolute dimensions of eclipsing binaries. XXVI. Setting a new standard: Masses, radii, and abundances for the F-type systems AD Bootis, VZ Hydrae, and WZ Ophiuchi

We determine absolute dimensions and abundances for the three F-type main-sequence detached eclipsing binaries AD Boo, VZ Hya, and WZ Oph and perform a detailed comparison with results from recent stellar evolutionary models. uvby light curves and standard photometry were obtained at ESO,radial velocity observations at CfA facilities, and high-resolution spectra with ESO's FEROS spectrograph. State-of-the-art methods were applied for the analyses. Masses and radii that are precise to 0.5-0.7% and 0.4-0.9%, respectively, have been established for the components, which span the ranges of 1.1 to 1.4 M_sun and 1.1 to 1.6 R_sun. The [Fe/H] abundances are from -0.27 to +0.10, with uncertainties between 0.07 and 0.15 dex. We find indications of a slight alpha-element overabundance of [alpha/Fe] ~ +0.1$ for WZ Oph. The secondary component of AD Boo and both components of WZ Oph appear to be slightly active. Yale-Yonsai and Victoria-Regina evolutionary models fit the components of AD Boo and VZ Hya almost equally well, assuming coeval formation, at ages of about 1.75/1.50 Gyr (AD Boo) and 1.25/1.00 Gyr (VZ Hya). BaSTI models, however, predict somewhat different ages for the primary and secondary components. For WZ Oph, the models from all three grids are significantly hotter than observed. A low He content, decreased envelope convection coupled with surface activity, and/or higher interstellar absorption would remove the discrepancy, but its cause has not been definitively identified. We have demonstrated the power of testing and comparing recent stellar evolutionary models using eclipsing binaries, provided their abundances are known. The strongest limitations and challenges are set by T_eff and interstellar absorption determinations, and by their effects on and correlation with abundance results.


Introduction
Accurate stellar mass and radius data from eclipsing binary systems are valuable empirical test data for stellar evolution models, basically because they are free of any scale dependent calibrations. Additional tests become possible when the effective temperatures for the components can also be determined, most often from well-calibrated photometry. The binaries then also serve as primary distance indicators, e.g. for nearby stars, stellar clusters, and Local Group galaxies, and they define the empirical massluminosity relations over a broad mass range. However, the most stringent model tests require that spectroscopic element abundances are available as well, minimising the number of free parameters in the comparison with theory (e.g. Andersen 1991). Unfortunately, the number of binary components with complete and accurate data is still small.
In our efforts to provide such complete data for a larger sample of eclipsing binary systems, F-G type main-sequence stars have been given high priority for two reasons. First, their spectra offer favourable conditions for accurate radial velocity determinations and abundance analyses. Second, such stars cover the mass range within which convective cores begin to develop and affect the observational determination of stellar ages, in particular due to difficulties and shortcomings related to the physics of core overshoot. Moreover, these ages cover the range of interest in studies of the evolution of the Galactic disk (e.g. Holmberg et al. 2007).
In this paper we present new analyses of the three detached F-type double-lined eclipsing binaries AD Boo, VZ Hya, and WZ Oph. From earlier studies they are known to have masses in the range 1.15-1.41 M ⊙ and radii corresponding to the lower half of the main-sequence band, but better accuracy and -especially -abundance data are needed for meaningful tests of state-of-the art stellar models. We have therefore obtained new radial velocity observations, light curves, standard photometry, and highresolution spectra covering a wide wavelength range, enabling us to determine all relevant parameters of the systems with high precision.
As pointed out by Popper (1998a), and also seen in Fig. 1 1 , AD Boo is of particular interest for having components with the greatest differences in mass and radius among well-studied binaries in its mass range. Hence, it may provide stringent tests of stellar models, but the masses obtained in two recent analyses of the system disagree by slightly more than their stated errors (Lacy 1997a;Popper 1998a). Moreover, in both studies the radii were based on reanalyses of the relatively old B, V light curves (normal points) by Zhai et al. (1982), which have an accuracy of only about 0.02 mag and show out-of-eclipse variations by about 0.075 mag.
Previous determinations of the dimensions of VZ Hya and WZ Oph were also based on old material (Popper 1965) and therefore had an accuracy of only 2-4%; see Popper (1980). As in the case of AD Boo, the components of VZ Hya are fairly different, whereas those of WZ Oph are nearly identical.
The paper is structured as follows: First, Sect. 2 describes our observations and analysis techniques adopted. The spectroscopic and photometric analyses of the individual systems are then discussed in detail in Sect. 3. The resulting new, accurate absolute dimensions and abundances are presented in Sect. 4 and compared with three recent grids of stellar evolutionary models in Sect. 5.
Throughout the paper, the component eclipsed at the deeper eclipse at phase 0.0 is referred to as the primary (p) component, the other stars as the secondary (s).

Spectroscopic observations
The radial velocity observations for the three binaries were conducted at the Harvard-Smithsonian Center for Astrophysics (CfA) with nearly identical echelle spectrographs on the 1.5-m Wyeth reflector at the Oak Ridge Observatory (Harvard, Massachusetts), the 1.5-m Tillinghast reflector at the F. L. Whipple Observatory (Mt. Hopkins, Arizona), and the Multiple Mirror Telescope (also 1 http://www.astro.keele.ac.uk/∼jkt/ Fig. 1. Main-sequence eclipsing binaries with both components in the 1.1-1.5 M ⊙ range having masses and radii accurate to 2% or better. The nine systems shown are, in order of increasing mass of the primary component: EW Ori, UX Men, HS Hya, V505 Per, IT Cas, V1143 Cyg, CD Tau, DM Vir, and BW Aqr. WZ Oph, VZ Hya, and AD Boo, are shown with star symbols. The full drawn line is the Y 2 0.3 Gyr isochrone for Z = 0.0181 (Demarque et al. 2004).
on Mt. Hopkins) prior to its conversion to a monolithic mirror. A single echelle order 45Å wide was recorded with intensified photon-counting Reticon detectors, at a central wavelength of approximately 5187Å, which includes the Mg I b triplet. The resolving power of these instruments is λ/∆λ ≈ 35 000, and the signal-to-noise (S/N) ratios achieved per resolution element of 8.5 km s −1 range from 10 to 43 for AD Boo, 15 to 39 for VZ Hya, and 10 to 25 for WZ Oph. As seen in Figures 4,9,and 14, the observations cover the orbits well.
For detailed abundance studies, we have obtained highresolution spectra at higher S/N ratios with the FEROS fiber echelle spectrograph at the ESO 1.52-m and 2.2-m telescopes at La Silla, Chile (Kaufer et al. 1999(Kaufer et al. , 2000. The spectrograph, which resides in a temperature-controlled room, covers without interruption the spectral region from the Balmer jump to 8700Å, at a constant velocity resolution of 2.7 km s −1 per pixel (λ/∆λ = 48 000). An observing log is given in Table 1.
A modified version of the MIDAS FEROS package, prepared by H. Hensberge, was used for the basic reduction of the FEROS spectra 2 3 . Compared to the standard package, background removal, definition and extraction of the individual orders, and wavelength calibration are significantly improved. The observations were reduced night by night using calibration exposures (ThAr and flat field) obtained during the afternoon/morning. Standard (rather than optimal) extraction was applied, and no order merging was attempted. Typically, a standard error of the wavelength calibration of 0.002-0.003Å was obtained. Subsequently, dedicated IDL 4 programs were applied to remove cosmic ray events and other defects, and for normalization of the individual orders. For each order, only the central part with acceptable signal-to-noise ratios was kept for further analysis. A short region of sample FEROS spectra is shown in Fig. 2.

Radial velocity determinations
Radial velocities were obtained from the CfA spectra using the two-dimensional cross-correlation algorithm TODCOR (Zucker & Mazeh 1994). Templates were selected from a library of synthetic spectra based on model atmospheres by R. L. Kurucz (see also Nordström et al. 1994;Latham et al. 2002), calculated over a wide range in effective temperature (T eff ), rotational velocity (v sin i), surface gravity (log g), and metallicity ([M/H]). The templates producing the best match to each component of the binary were identified by running extensive grids of cross-correlations and selecting the combination giving the highest correlation value averaged over all exposures, weighted according to the strength of each spectrum (e.g. Torres et al. 2002). These grids were in general run in temperature and rotational velocity only, since those are the parameters that affect the velocities the most. Surface gravities were fixed according to preliminary estimates. If needed, grids for different assumed heavy element abundances were compared, leading indirectly to a rough estimate of [M/H] for the system, which was in turn compared with the results from the detailed abundance analysis. The stability of the zero-point of our velocity system was monitored by means of exposures of the dusk and dawn sky, and small run-to-run corrections were applied as described by Latham (1992). In addition to the radial velocities we derived the light ratio at the mean wavelength of our observations (5187Å), following Zucker & Mazeh (1994).
Because of the narrow wavelength coverage of the CfA spectra, the radial velocities were checked for systematic errors that might occur as a result of lines of the stars moving in and out of the spectral window with orbital phase (Latham et al. 1996). Numerical simulations were carried out as described in more detail by Torres et al. (1997), and corrections based on them were applied to the raw velocities. These corrections are included in the measurements reported in Tables A.1, A.2, and A.3.

Photometric observations
uvby light curves for AD Boo, VZ Hya, and WZ Oph were obtained with the Strömgren Automatic Telescope (SAT) at ESO, La Silla during several campaigns between March 1988 and March 1997. They contain 652, 1180, and 697 points in each band, respectively, and cover all phases well. The accuracy per point is 4-5 mmag in vby and 6-7 mmag in u, far better than previous photometry of these systems. The light curves are kept in the stable uvby instrumental system of the SAT in order not to introduce additional scatter from transformation to the standard uvby system. Improved linear ephemerides have been calculated for all three systems from the new and numerous published times of minima. No signs of period changes were detected, and the orbits were assumed to be circular since there is no evidence for eccentricity, either in previous data or in our own.
In addition to the light curves, homogeneous standard uvbyβ indices have been established from SAT observations on dedicated nights in which the binaries (and comparison stars) were observed together with a large sample of uvby and β standard stars.
We refer to Clausen et al. (2008, hereafter CVG08) for further details on the photometry and the new ephemerides.

Photometric analyses
The choice of binary model for the light curve analysis depends strongly on the relative radii of the components and the irradiation between them. In our case, the systems are well-detached with relative radii between 0.10 and 0.17, the components are little deformed (oblateness between 0.001 and 0.007), and reflection effects are small. Therefore, the simple Nelson-Davis-Etzel model (Nelson & Davis 1972, Etzel 1981, Martynov 1973, which represents the deformed stars as biaxial ellipsoids and applies a simple bolometric reflection model, is expected to be perfectly adequate. We have adopted the corresponding EBOP code (Popper & Etzel 1981), supplemented by an extended version JKTEBOP 5 including Monte Carlo simulations (see e.g. Southworth et al. 2004aSouthworth et al. , 2004b, for the analyses of the three systems and the assignment of realistic errors for the photometric elements. Since the accuracy of the light curves is high, we have also checked for model-dependent systematic errors in the photometric elements by performing parallel analyses with two more advanced binary models, WINK (Wood 1972) and the widely used Wilson-Devinney (WD) model (Wilson & Devinney 1971;Wilson 1993Wilson , 1994. WINK uses tri-axial ellipsoids, whereas WD is based on Roche geometry. Our WINK version is modified and extended as described by Vaz (1984Vaz ( , 1986, Vaz & Nordlund (1985), and Nordlund & Vaz (1990). Our WD code is a revised and extended version of the 1994 code by Wilson; see Vaz et al. (1995) and Casey et al. (1998). In both codes we use the predictions of ATLAS9 model atmospheres (Kurucz 1993) to connect the flux ratios of the two components in different passbands.
In their publicly available versions, EBOP, WINK, and WD carry out a least-squares optimization of the parameters by means of the classical differential corrections method. In JKTEBOP, this has been replaced by the Levenberg-Marquardt minimization algorithm (MRQMIN: Press et al. 1992).
Throughout this paper we have analysed the uvby light curves independently, and for each band identical weights have been assigned to all observations. A photometric scale factor (the magnitude at quadrature) was always included as an adjustable parameter, and the phase of primary eclipse was allowed to shift from 0.0. The mass ratio between the components was kept at the spectroscopic value, and synchronous rotation was assumed. Gravity darkening coefficients/exponents corresponding to convective atmo-spheres were applied. The bolometric reflection albedo coefficients were fixed at 0.5 in the WINK and WD analyses, again due to convection; for EBOP the simple bolometric reflection model was used.
Linear limb-darkening coefficients were either assigned from theoretical calculations (van Hamme 1993;Díaz-Cordovés et al. 1995;Claret 2000) according to the effective temperatures, surface gravities, and abundances, or included as adjustable parameters. For all components of the three systems, the linear coefficients by Claret (2000) are about 0.10 higher than those by van Hamme (1993) and about 0.03 higher than those by Díaz-Cordovés et al. (1995). In general, a change in the coefficients by 0.10 causes systematic differences between the optimum theoretical light curves of just a few mmag through the eclipses, and the overall quality of the fits to the observations is nearly identical. If adjusted, the coefficients for the components of VZ Hya and WZ Oph turn out closest to the values by van Hamme; for AD Boo the coefficients either become unrealistic or convergence is not achieved at all. In addition, a few WD analyses with the non-linear law by Klinglesmith & Sobieski (1970) and coefficients by van Hamme (1993) were included. No improvements in the light curve fits were seen. The influence on the other photometric elements of adopting different limb-darkening coefficients/laws are discussed below for the individual systems. In order to include uncertainties of the limb-darkening coefficients in the uncertainties of the final photometric elements, random variations within +/-0.05 have been introduced in the JKTEBOP Monte Carlo simulations mentioned above.
In tables and text on photometric solutions we use the following symbols: i orbital inclination; r relative radius; k = r s /r p ; u linear limb darkening coefficient; y gravity darkening coefficient; J central surface brightness; L luminosity; Ω surface potential; T eff effective temperature.

Effective temperatures and interstellar reddening
Although significantly improved (semi)empirical temperature scales and calibrations have become available during the last decade, the determination of effective temperatures of binary components remains a difficult task; see e.g. Popper (1980), Andersen (1991), andClausen (2004) for general remarks. For comparison with theoretical models, luminosity calculations, and distance determinations, it is on the other hand important to diminish both systematic and random errors as much as possible. In the analyses of AD Boo, VZ Hya, and WZ Oph we have chosen to follow three independent approaches: First, as part of the light curve analyses we calculate the individual uvby indices for the components from the combined standard uvby indices (CVG08) and the luminosity ratios between the components in the four bands (based on mean geometric elements). Interstellar reddening is determined from the intrinsic colour calibration by Olsen (1988), which has a precision of 0.009 mag, and the results are compared to those from other sources. Effective temperatures can then be obtained from a series of (semi)empirical calibrations such as those by Alonso et al. (1996), Ramírez & Meléndez (2005), and Holmberg et al. (2007). Also, as an important check on these calibrations, effective temperature differences between the components can be derived from the empirical flux calibration by Popper (1980) and the y surface flux ratios. They are obtained to high precision in the light curve analyses, since they depend mainly on the observed depths of primary and secondary eclipses. We have applied both the BC scale used by Popper (1980) and that by Flower (1996), which in the temperature range of AD Boo, VZ Hya, and WZ Oph, lead to nearly identical temperature differences with errors of only about 25 K (excluding uncertainties in Popper's flux scale).
Second, we attempt to obtain spectroscopic temperatures from the high-resolution and high-S/N spectra, either combined or disentangled, as part of the abundance analyses; see e.g. Santos et al. (2004) and Bruntt et al. (2004) for details. And thirdly, as mentioned in Sect. 2.2, temperatures are estimated by identifying the templates producing the best match to each component in the TODCOR analyses, and interpolating in log g or metallicity as necessary. As discussed below for the individual systems, the two latter methods however tend to be problematic.
Finally, the photometric indices of the combined light of a system may provide a first estimate of the effective temperature of its "average" component, and if the components are nearly identical, as is the case for WZ Oph, a direct determination is possible. In Table 2 we present such temperatures derived from the combined uvby indices and 2MASS 6 JHK s observations. As seen, the temperatures from the different uvby and V JHK s calibrations agree well for WZ Oph. The spread for AD Boo and VZ Hya presumably reflects the fact that their combined flux distribu-6 http://www.ipac.caltec.edu/2mass tions are distorted compared to a single star, because the components differ by 400-600 K in effective temperature.

Abundance analyses
For meaningful comparisons between theoretical stellar models and empirical absolute dimensions -as determined from detached double-lined eclipsing binaries -the chemical composition, and in particular the metal abundance, must be known for the binary. However, detailed abundance studies have only been done for a few systems. As part of the analyses of AD Boo, VZ Hya, and WZ Oph, we have used the new spectroscopic and photometric data to establish their metal abundance. Several methods were applied/attempted:  Simon & Sturm (1994) and a revised version of the corresponding original code developed by E. Sturm. It assumes a constant light level, so only spectra taken outside of eclipse can be included. We note that all three systems are in fact constant to within a few percent outside of eclipse. The orbital elements were fixed at the results based on the TODCOR analyses of the CfA spectra. For AD Boo, with only two FEROS spectra available, disentangling is not possible. Disentangling of the five VZ Hya spectra was attempted, but the quality of the resulting component spectra was not acceptable; more than five spectra and/or a better phase distribution seem to be needed. For both systems the observed double-lined spectra were therefore used directly for the abundance analyses.
The versatile semi-automatic IDL tool VWA 7 , developed by Bruntt, was applied for the abundance analyses. We refer to Bruntt et al. (2002Bruntt et al. ( , 2004Bruntt et al. ( , 2008 for further details. For our particular study, VWA has been extended to allow the analysis of double-lined spectra in addition to normal spectra from single stars -or disentangled binary component spectra. We have adopted the effective temperatures, surface gravities, and rotational velocities determined for the binary components in this study (Table 21), but for comparison analyses for higher or lower temperatures were also done. Microturbulence velocities were derived from the calibration by Edvardsson et al. (1993), which has a scatter of about 0.3 km s −1 , but have been tuned, if needed, as part of the analyses.
In each spectrum, VWA selects the least blended lines from atomic line lists extracted from the Vienna Atomic Line Database (VALD; Kupka et al. 1999). For each selected line, the synthetic spectrum is then calculated, and the input abundance is changed iteratively until the equivalent width of the observed and synthetic lines match. The synthetic spectra are generated with the SYNTH software Table 2. Effective temperatures (K) of the 'average' components of AD Boo, VZ Hya, and WZ Oph. The V 0 magnitudes and the (b−y) 0 and c 0 indices are based on the out-of-eclipse uvby standard indices from CVG08. The 2MASS observations (J, H, K s ) were obtained at phases 0.216 (AD Boo), 0.181 (VZ Hya), and 0.495 (WZ Oph), respectively. For WZ Oph, the V magnitude was calculated from the value outside eclipses and the y eclipse depth at phase 0.495. A V is the adopted visual interstellar absorption. References are: A96, Alonso et al. (1996); H07, Holmberg et al. (2007); RM05, Ramírez & Meléndez (2005); M06, Masana et al. (2006). The results from A96 are based on their uvby calibration, those from RM05 on their uvby, (V − J), (V − H), and (V − K s ) calibrations (in that order); the calibration by M06 is for (V − K s ).  (Valenti & Piskunov 1996). Atmosphere models were either interpolated from the grid of modified ATLAS9 models by Heiter et al. (2002) or calculated from the MARCS code.
We have applied an improved version of the original model by Gustafsson et al. (1975); see Jørgensen et al. (1992) for further details.
For the binary components, we calculate all abundances relative to the Sun, line by line. We have used either an atlas spectrum of the solar flux (Wallace et al. 1998) or a FEROS sky spectrum. In general, this approach gives better agreement between the abundances determined from the individual lines of a given element/ion.
If we adopt the VALD line data and the solar photosphere abundances by Grevesse & Sauval (1998), we obtain for the Sun the Fe abundances listed in Table 3 for the solar atlas and FEROS spectra. Only lines with equivalent widths between 10 and 80 mÅ measured in both spectra were included. It is gratifying to see that the results from the FEROS spectrum agree very well with those from the atlas spectrum. For both spectra, Fe I and Fe II lines give slightly different [Fe/H] results. They agree a bit better for the MARCS models than for ATLAS9 models, and the MARCS models also reproduce the solar result slightly better.
Although we will focus on Fe, we have also analysed other elements, and results from the solar atlas spectrum are given in Table 4. For most of the ions with many measured lines we obtain good agreement with Grevesse & Sauval (1998), but there are exceptions such as Si I, Cr II, and Mn I. This may indicate errors in the VALD line data (a problem with Mn was found by Bruntt et al. 2004) which should, however, be eliminated to a large extent in the binary abundance analyses performed relative to the Sun. Comparison with the most recent solar photosphere abundances by Grevesse et al. (2007), which for many elements are derived using updated line data, 3D models and departures from LTE, is less obvious, since we apply 1D, LTE atmosphere codes.

Disentangled CfA spectra
As seen in Tables 5, 11, and 16, a large number of CfA spectra are available for all three systems. We have selected those collected outside eclipses with sufficient S/N ratios for disentangling. The actual numbers of spectra used are 113 (AD Boo), 29 (VZ Hya), and 40 (WZ Oph). After careful normalization of the individual spectra, the disentangled component spectra were constructed by adopting the orbital elements listed in the tables mentioned above. The luminosity ratios between the components, which cannot be determined from the disentangling procedure, were calculated by interpolation between the b and y luminosity ratios of the photometric solutions. The disentangled component spectra were then compared to the CfA library of synthetic spectra mentioned in Sect. 2.2. In general we find that the short wavelength range covered (45Å), and the correspondingly small number of lines available, limits the success of this approach. Most importantly, it was difficult to break the correlation between metal abundance and effective temperature. Results for the individual systems will, however, be presented below.

TODCOR template matching
As mentioned in Sect. 2.2, rough estimates of the metal abundance for a system can be obtained directly from the observed composite spectra by running extensive grids of two-dimensional cross-correlations (TODCOR) using templates of different metallicity. Results from this exercise will also be discussed below for the individual systems.

uvby calibration
Finally, we have derived metal abundances from the uvby indices for the individual binary components (Table 21) and the recent calibrations by Holmberg et al. (2007), which are based on a carefully selected sample of spectroscopic [Fe/H] results for several hundred stars. The calibration adopted for 0.30 < b − y < 0.46 has a dispersion of only 0.07 dex; that for bluer stars has a dispersion of 0.10 dex. Taking into account the uncertainties of the uvby indices and the Table 4. Abundances for the Sun determined from Solar atlas spectrum (Wallace et al. 1998) and ATLAS9 models (Heiter et al. 2002). N is the number of lines used; only ions with three or more measured lines are listed. GS1998 is the logarithmic element abundances (H = 12.00) from Grevesse & Sauval (1998

Spectroscopic and photometric analyses
In this section we present in detail the spectroscopic and photometric analyses of AD Boo, VZ Hya, and WZ Oph, including abundances, and compare our results with those obtained in previous studies. Additional references to earlier work are given by CVG08.

Radial velocities and spectroscopic orbit for AD Boo
Radial velocities for the components of AD Boo were determined by TODCOR analyses of 124 usable CfA spectra outside of eclipse, observed from January 1989 to June 1997. Surface gravities (log g) were fixed at 4.0 (primary) and 4.5 (secondary), and the heavy element abundance was assumed to be solar, the closest values in our grid of synthetic spectra to the actual values determined in Sect. 4. We inferred temperatures of 6250±150 K for the primary of AD Boo (the more massive star), and 5680 ± 300 K for the secondary, both significantly lower than the photometric values, because the templates have lower metallicity than measured for AD Boo; see Sects. 3.1.3 and 4. The measured v sin i values are 38 ± 2 km s −1 and 37 ± 5 km s −1 , respectively, with the secondary being considerably more uncertain due to its much weaker lines. The synchronous values are 39 km s −1 for the primary and 30 km s −1 for the secondary. Lacy (1997a) reported 38 ± 1 and 31 ± 1 km s −1 , whereas Popper (1998a) did not give any values. Templates from our library with parameters closest to the measured values were used to derive the radial velocities.
In the case of AD Boo we found that the radial velocity corrections for systematics described in Sect. 2.2, which are shown in Fig. 3, were quite significant, reaching up to 10 km s −1 for the secondary star and up to 2 km s −1 for the primary. The effect on the absolute masses is large, ∼6% for the primary and ∼4% for the secondary. The final velocities including corrections are listed in Table A.1. The light ratio derived from the spectra, including also corrections for systematics analogous to those derived for the velocities, is ℓ s /ℓ p = 0.37 ± 0.02 at the mean wavelength of our observations (5187Å). The spectroscopic orbital solution for AD Boo is given in Table 5, and the observations and computed orbit are shown graphically in Figure 4 along with the residuals for each star. For comparison, we have also disentangled 113 Table 5. Spectroscopic orbital solution for AD Boo. T is the time of central primary eclipse.

Photometric elements for AD Boo
The uvby light curves of AD Boo contain 652 observations in each band and were observed on 42 nights during five periods between March 1988 and March 1992 (CVG08). The accuracy per point is about 0.005 mag (vby) and 0.007 mag (u), but comparison of the data from the five periods reveals that at some phases, systematic differences in the light level of 0.01-0.03 mag exist, increasing from y to u. In particular, observations during the 1989 eclipses are systematically fainter. The effect is shown in Fig. 5 and indicates surface activity (spots) at a rather low level on one or both components, the cooler secondary being the most likely candidate. Zhai et al. (1982) also noticed scatter in their B and V light curves, which were obtained from April to June 1981, but at a much higher level of about 0.075 mag compared to their observing error of about 0.02 mag.
There is very little information on activity in the literature, and AD Boo does not appear to be an X-ray source. The presence of chromospheric activity in terms of the Rossby number 8 was studied by Hall (1994). The secondary in AD Boo is perhaps cool enough to display significant activity (see above) despite having a mass as high as 1.21 M ⊙ , which would imply a relatively thin convective envelope. An estimate of the color index of B − V = 0.58 based on the temperature allows the convective turnover time to be inferred, and this together with the rotational period (assuming synchronous rotation) leads to a Rossby number of approximately 0.20, which is in the range where other stars tend to show spot activity (see Hall 1994). Furthermore, our two high-resolution FEROS spectra taken in January 1999 (see Table 1) show in fact a weak emission feature in the Ca II H and K lines at the position of the secondary component, supporting also some activity. Nothing is seen in the H α line.
EBOP and WD solutions, based on linear limbdarkening coefficients by van Hamme (1993) are presented in Tables 6 and 7, respectively, and O −C residuals of the v observations from the theoretical light curve computed for the photometric elements given in Table 6 are shown in Fig. 5. As seen, differences between the photometric elements obtained from the two models are small, and the mean values of r p + r s and k agree within their uncertainties. WINK solutions, which are not included, agree well with those from EBOP and WD. Slightly different orbital inclinations are obtained, however, for the four bands from all three codes. This could be due to either the systematic light level differences mentioned above (spots), which increase from y to u, to a small amount of (red) third light, or to a slightly incorrect theoretical wavelength dependence of the limb-darkening coefficients.
Excluding the 1989 observations (see Fig 5), we obtain the EBOP solutions listed in Table 8. The relative Table 6. Photometric solutions for AD Boo (all data) from the EBOP code adopting linear limb darkening coefficients by van Hamme (1993). The errors quoted for the free parameters are the f ormal errors determined from the iterative least squares solution procedure.  Table 7. Photometric solutions for AD Boo (all data) from the WD code adopting linear limb darkening coefficients by van Hamme (1993); see Table 6. Gravity darkening exponents of 0.33 and bolometric albedo coefficients of 0.5 were applied, as appropriate for convective envelopes. T eff,p was assumed to be 6575 K. radii change by about 0.7% in opposite directions and their sum by about −0.6%, whereas the orbital inclination is unchanged in all four bands. Including 2% of third light changes i by about 0. • 3, which is the maximum difference between the four bands, but we have not found any direct support for even such a small amount of third light in any of the analyses. If we adopt linear limb-darkening coefficients from Claret (2000), which are 0.04-0.07 higher, the orbital inclination decreases by 0. • 16, and the relative radius increases by about 0.0013 for the primary component but is unchanged for the secondary. If adjusted, the linear limb-darkening coefficients either become unrealistic and the orbital inclinations from the four bands disagree much more, or convergence fails. Finally, WD solutions including non-linear limb darkening (see Sect. 2.4) give photometric elements very close to those presented in Table 7. We are   Table 6. 1989 observations are shown with square symbols. therefore unable to point at a clear explanation for the small i discrepancies, which are, however, of no consequence for the determination of the absolute properties.
In order to assign realistic errors to the photometric elements, we have supplemented these analyses with 10 000 JKTEBOP Monte Carlo simulations in each band; see Fig. 6. Furthermore, as shown in Fig. 7, we have com-  pared EBOP solutions for a range of fixed k = r s /r p values near 0.75. The adopted photometric elements listed in Table 9 are the weighted mean values of the EBOP solutions adopting the linear limb-darkening coefficients by van Hamme. Errors are based on the Monte Carlo simulations and comparison between the uvby solutions, including those where the 1989 observations were excluded. We find that secondary eclipse is in fact total, whereas at phase 0.0 about 61% of the y light from the primary component is eclipsed. The spectroscopic light ratio of ℓ s /ℓ p = 0.37±0.02 at 5187Å agrees quite well with that from the light curve solutions (0.40). For comparison, Lacy (1997a), who applied the EBOP code to analyse the B and V normal points by Zhai et al. (1982)

Abundances for AD Boo
Two high S/N FEROS spectra are available for abundance analyses of AD Boo; see Table 1. We have concentrated on the one with the highest S/N, which gives significantly more accurate abundances. Due to the relatively high rotational velocities and low luminosity ratio, the number of suitable unblended lines is far lower than for VZ Hya and WZ Oph, especially for the secondary component. ATLAS9 models were used, and the effective temperatures, surface gravities and rotational velocities listed in Table 21 were adopted, together with microturbulence velocities of 1.5 km s −1 for both components. The temperatures, which are derived from the uvby indices of the components, could not be constrained further from the spectroscopic analysis.
As seen in Metallicities have also been derived by comparing the disentangled CfA spectra to the extensive library of synthetic spectra especially prepared for the 45Å region centered at 5187Å. Surface gravities were fixed as above, and synchronous rotation of the components was adopted. The best fits are obtained for [M/H] = +0.13 and T eff = 6450 K (primary), and [M/H] = +0.30 and T eff = 6200 K (secondary), but a strong correlation between metallicities and Within the fairly large errors, the results from the different methods are in reasonable agreement. Giving higher weight to the analysis of the FEROS spectra, we adopt [Fe/H] = +0.10 ± 0.15 for AD Boo. This result should, however, be checked and preferably improved via analyses of disentangled component spectra based on a sufficiently large number of high-dispersion spectra.

VZ Hya
VZ Hya = HD72257 = HIP41834 (P = 2. d 90) was discovered by O'Connell (1932) to be an eclipsing binary, and spectroscopic elements were obtained from 75Å/mm spectra by Struve (1945). The currently most reliable absolute dimensions of the components of VZ Hya are those given by Popper (1980) in his critical review on stellar masses.

Radial velocities and spectroscopic orbit for VZ Hya
Our TODCOR analysis of VZ Hya was based on 42 usable CfA spectra outside of eclipse, observed from January 1989 to April 1991. Templates of solar metallicity, log g = 4.5 (both components), and temperatures of 6750 K (primary) and 6250 K (secondary) were used to derive the radial velocities. As seen in Fig. 8, TODCOR corrections for systematics in the radial velocities are fairly small, typically under 1 km s −1 for both components. They make very little difference in the masses, changing them by only +0.4% (primary) and −0.1% (secondary), which is well within the errors. The scatter around the orbit is about the same with and without corrections applied. Interpolating between log g = 4.0 and 4.5 grids to the exact gravities given in Table 21, we obtain temperatures of 6590 ± 100 K for the primary and 6150±150 K for the secondary, both in excellent agreement with the photometric values; see Sect. 4. The corresponding light ratio is ℓ s /ℓ p = 0.53 ± 0.02 at the mean wavelength of our observations (5187Å). Finally, the grids of crosscorrelations give best v sin i values of 21 and 20 km s −1 , with errors of 2 km s −1 and 3 km s −1 , respectively; they are consistent with synchronous rotation.
The spectroscopic orbital solution for VZ Hya is given in Table 11, and the observations and computed orbit are shown graphically in Figure 9 along with the residuals for each star. Our velocity semiamplitudes agree well with those by Popper (1965), (K p , K s ) = (94.3 ± 2.3, 103.9 ± 1.8) but are significantly more accurate.
For comparison, we have also disentangled the 29 best CfA spectra outside of eclipse and obtain velocity semiamplitudes of K p = 94.62 and K s = 105.87, in good agreement with the results based on TODCOR velocities.

Photometric elements for VZ Hya
The uvby light curves of VZ Hya contain 1180 observations in each band and were observed on 44 nights during four periods between February 1989 and April 1992 (CVG08). The average observational accuracy per point is about 4 mmag (vby) and 6 mmag (u), and throughout all phases the points scatter at these levels, indicating that the components of VZ Hya are constant within the precision of our    Table 12. Gravity darkening exponents of 0.33 and bolometric albedo coefficients of 0.5 were adopted, as appropriate for convective envelopes. T eff,p was assumed to be 6650 K. Tables 12 and 13, respectively. As seen, the results from the different bands agree very well, and the two models lead to practically identical photometric elements. The same holds for WINK solutions not included here. If we adopt instead linear limb-darkening coefficients by Claret (2000), which are 0.05-0.08 higher, the orbital inclination decreases by 0. • 2, and the relative radius is unchanged for the primary component but increases by 0.0009 for the secondary. If adjusted, the linear limb-darkening coefficients by van Hamme are reproduced to within ±0.05 and the solutions are close to the results in Tables 12 and 13. WD solutions including non-linear limb darkening (see Sect. 2.4) give photometric elements very close to those presented in Table 13. O − C residuals of the b observations from the theoretical light curve computed for the photometric elements given in Table 12 are shown in Fig. 10.  Table 12. As done for AD Boo, in order to assign realistic errors to the photometric elements, we have supplemented these analyses with 10 000 JKTEBOP Monte Carlo simulations in each band, and EBOP solutions for a range of fixed k values near 0.85; see Fig. 11 and Fig. 12. The adopted photometric elements listed in Table 14 are the weighted mean values of the EBOP solutions adopting the linear limb-darkening coefficients by van Hamme. Errors are based on the Monte Carlo simulations and a comparison between the uvby solutions. Our results are significantly more accurate than obtained in previous studies of VZ Hya. At phase 0.0, about 76% of the y light from the primary component is eclipsed, and at phase 0.5 about 98% of the y light of the secondary is blocked out, meaning that secondary eclipse is almost total. The spectroscopic light ratio of ℓ s /ℓ p = 0.53 ± 0.02 at 5187Å agrees quite well with that from the light curve solutions (0.56).
Previous analyses based on the photoelectric U BV light curves by Walker (1970) gave 3-4% larger relative radii for the primary component, whereas those for the secondary component agree well with our result (Walker 1970, Wood 1971, Cester et al. 1978. Another set of photoelectric U BV light curves was published by Padalia & Srivastava (1975), who found primary eclipse to be an occultation (larger star in front) rather than a transit (smaller star in front). This picture of VZ Hya was soon shown by Popper (1976) to be  wrong but was nevertheless defended once more by Padalia (1986).

Abundances for VZ Hya
Five FEROS spectra are available for the detailed abundance analysis of VZ Hya; see Table 1. As mentioned in Sect. 2.6, disentangling of so few spectra was unsuccess-  ful, and we have instead used the three most suitable observed spectra and analysed them as double-lined. The results from the VWA analysis are presented in Table 15. ATLAS9 models were used, and the final effective temperatures, surface gravities and rotational velocities listed in Table 21 were adopted, together with microturbulence velocities of 2.0 and 1.6 km s −1 , respectively, calculated from the calibration by Edvardsson et al. (1993). Excluding ions with few measurements, we find that within the uncertainties, identical abundances are derived from Cr I and Ni I lines. The α-elements Si and Ca seem slightly enhanced, but the difference is barely significant.
For the metallicities derived by comparing the disentangled CfA spectra with the library of synthetic spectra, surface gravities were fixed as above, and synchronous rotation of the components was adopted. The best fits are obtained for [M/H] = +0.09 and T eff = 6630 K (primary), and [M/H] = −0.23 and T eff = 6130 K (secondary), but as for AD Boo, a strong correlation between metallicities and temperatures exist, and the determinations are therefore quite uncertain. Forcing the metallicities of the components to be identical, we derive [M/H] = −0.15, and temperatures of 6595 and 6185 K for the components, in good agreement with the FEROS result and the photometric temperatures. Estimated uncertainties of these CfA metallicities and temperatures are (at least) 0.1 dex and 100 K, respectively.
Comparing TODCOR analyses of the observed CfA spectra for a range of adopted template metallicities yields a metallicity estimate of [M/H] = +0.02, assuming temperatures of 6650 and 6300 K and forcing the metal abundance of the components to be identical.
Finally, the calibrations by Holmberg et al. (2007) and the de-reddened uvby indices in Table 21

WZ Oph
WZ Oph = HD154676 = HIP83719 (P = 4. d 18) was discovered by J. H. Metcalf to be an eclipsing binary (Pickering 1917), and Sanford (1937) obtained the first radial velocities. The currently best modern analysis of WZ Oph, leading to masses and radii accurate to 3-4%, was done by Popper (1965) from photoelectric light curves and 20 A/mm spectra.
For the spectroscopic and photometric analyses of WZ Oph we have adopted the well established ephemeris from CVG08, where the zero point has been shifted by half a period compared to that by Popper (1965). We note that the primary component eclipsed at phase 0.0 according to our adopted ephemeris, and corresponding to the slightly deeper minimum for the uvby light curves, is given by Popper (1965) as component 2. Popper found that his component 2, i.e. our primary component, has slightly stronger spectral lines. This is consistent with our measurements of equivalent widths in the disentangled FEROS spectra for lines stronger than about 100 mÅ.

Radial velocities and spectroscopic orbit for WZ Oph
Radial velocities for the components of WZ Oph were determined with TODCOR on the basis of 40 usable CfA spectra outside of eclipse, observed from February 1989 to September 1991. Identical templates of solar composition with log g = 4.0 and a temperature of 6250 K were adopted for the two stars. The TODCOR corrections for systematics shown in Fig. 13 are less than 0.5 km s −1 at all phases for both components. This makes a difference of only +0.3% in the minimum masses, about half their uncertainties, and the scatter of the orbit does not change much with corrections. v sin i values of 18 and 17 km s −1 are obtained for the primary and secondary components, respectively, with estimated uncertainties of 1 km s −1 . They are consistent with synchronous rotation.
Our best spectroscopic estimates of the effective temperatures of the components by interpolating to a log g value close to that listed in Table 21 are 6260 K for the primary and 6340 K for the secondary, assuming solar metallicity. Interpolating to the metallicity of [M/H] = −0.27 (see below) gives 6050 K and 6120 K. The lower values for the primary compared to the secondary disagree formally with the trend from photometry, which gives a slightly hotter primary based on the slightly deeper primary eclipse. However, the difference is well within the spectroscopic error estimate of 150 K. The light ratio is ℓ s /ℓ p = 0.98 ± 0.02 at the mean wavelength of our observations (5187Å).
The spectroscopic orbital solution for WZ Oph is given in Table 16 and the observations and computed orbit are shown graphically in Figure 14 along with the residuals for each star. For comparison, we have also disentangled the 40 CfA spectra outside of eclipse and obtain velocity semiamplitudes of K p = 88.62 and K s = 88.92, in good agreement with the results based on the TODCOR velocities.
Our velocity semiamplitudes are larger than those by Popper (1965), (K p , K s ) = (85.9 ± 1.8, 87.0 ± 2.0), and as in the case of VZ Hya they are of much higher precision.

Photometric elements for WZ Oph
The uvby light curves of WZ Oph contain 697 observations in each band and were observed on 41 nights during six   (1993); see Table 17. Gravity darkening exponents of 0.33 and bolometric albedo coefficients of 0.5 were adopted, as appropriate for convective envelopes. T eff,p was assumed to be 6165 K. periods between March 1991 and March 1997 (CVG08). The average observational accuracy per point is about 4 mmag (vby) and 6 mmag (u), but throughout all phases the points scatter by 1-4 mmag more than this, highest in u. At a given phase, the observations from different seasons do not differ systematically, and it is thus unclear if the extra scatter in the light curves is due to slight variability of one or both stars (activity), or perhaps to an unknown instrumental/observational effect not seen in other systems. The Rossby number inferred for the mean component of WZ Oph does not clearly place it in the group of stars that show variability due to spots (Hall 1994), and is thus inconclusive. On the other hand, the high S/N FEROS spectrum taken at phase 0.33 (see Table 1) indicates very weak emission from both components in the Ca II H and K lines, but not in H α . This suggests WZ Oph is in fact mildly active.
EBOP and WD solutions, based on linear limbdarkening coefficients by van Hamme (1993) are presented  Table 17. in Tables 17 and 18, respectively. As seen, the results from the different bands agree well, and the two models lead to nearly identical photometric elements. WINK solutions not included here lead to changes in the relative radii of −0.0008 (primary) and +0.0008 (secondary). Formally, the less massive secondary component is slightly larger and cooler, but more luminous, than the primary for all three models. In reality the two stars are, however, very nearly identical. If we adopt linear limb-darkening coefficients by Claret (2000), which are 0.10 higher than those by van Hamme, the orbital inclination decreases by 0. • 1, and the relative radius increases by about 0.0004 for the primary component but is unchanged for the secondary. If adjusted, the vby linear limb-darkening coefficients by van Hamme are reproduced to within ±0.05 and the solutions are close to the results in Table 17. The u coefficients, however, become unrealistically low (0.60), and i changes by 0. • 2. WD solutions including non-linear limb darkening (see Sect. 2.4) give photometric elements very close to those presented in Table 18 for the uvb bands, but the u solution now disagrees in having k ≈ 1.05. O − C residuals of the y observations from the theoretical light curve computed for the photometric elements given in Table 17 are shown in Fig. 15.
Realistic errors have been assigned to the photometric elements by performing 10 000 JKTEBOP Monte Carlo simulations in each band, and EBOP solutions for a range of fixed k values near 1.0; see Fig. 16 and Fig. 17. The adopted photometric elements listed in Table 19 are the weighted   find that about 92% of the light of the components is shadowed out at both central eclipses.

Abundances for WZ Oph
14 FEROS spectra are available of WZ Oph (see Table 1), which allowed us to obtain disentangled spectra of high quality for the abundance studies. Results from the VWA analysis are presented in Table 20. ATLAS9 models were adopted, and the effective temperatures, surface gravities and rotational velocities listed in Table 21 were used, together with microturbulence velocities of 1.7 km s −1 for both components, calculated from the calibration by Edvardsson et al. (1993). This choice of parameters is supported by the spectroscopic analyses, since we see no dependence on equivalent width and excitation potential for abundances determined from individual lines. The precision of the temperatures (±100 K) could not be improved from the abundance analysis.
As Excluding ions with only a few lines measured, we find that within the uncertainties identical abundances are derived for Cr and Ni, noting that Cr II deviates, as seen for the Sun (Table 4), although our abundances are relative so that the effect should cancel out. The α-elements Si and Ca seem to be enhanced by about 0.1 dex, whereas Na and Ti (with only few lines measured) appear normal. The Si result may perhaps be questioned, as in the case of Cr II above, since Si I gave a low abundance for the Sun. For Mn, which together with Al is expected to anticorrelate with α-elements, we actually derive a rather low abundance of −0.54±0.13, taking into account temperature and microturbulence uncertainties as above. Unfortunately, problems with the Mn I line were also seen for the Sun. We tend, however, to conclude, that a slight α-enhancement of about 0.1 dex may be present. Finally, we derive an abundance for V close to solar. As for AD Boo and VZ Hya, metallicities have also been derived from the disentangled CfA spectra. Surface gravities were fixed as above, and synchronous rotation of the components was adopted. The best fits are obtained for [M/H] = −0.14 and T eff = 6280 K (primary), and [M/H] = −0.29 and T eff = 6250 K (secondary), but as in the other two binaries there is a strong correlation between metallicity and temperature, and the determinations are therefore quite uncertain. Forcing the metallicities of the components to be identical, we derive [M/H] = −0.23 and temperatures of 6230 and 6310 K for the components, in good agreement with the FEROS abundance result and in reasonable agreement with the photometric temperatures. Estimated uncertainties of metallicities and temperatures from this method are (at least) 0.1 dex and 100 K, respectively.
TODCOR analyses of the observed (composite) CfA spectra for a range of fixed template metallicities yield an estimate of [M/H] = −0.27, assuming temperatures of 6165 and 6115 K and constraining the results for the components to be identical.
Finally, the calibration by Holmberg et al. (2007) and the de-reddened uvby indices yield [Fe/H] = −0.17 ± 0.11 and −0.24 ± 0.12 for the primary and secondary components, respectively, also in good agreement with the FEROS result, which we adopt for WZ Oph.

Absolute dimensions
Absolute dimensions for AD Boo, VZ Hya, and WZ Oph are presented in Table 21, as calculated from the new spectroscopic and photometric elements given in Tables 5, 9 , 11, 14, 16, and 19.
As seen, masses and radii precise to 0.5-0.7% and 0.4-0.9%, respectively, have been established for the binary components. For AD Boo, the masses agree well with those determined by Popper (1998a) but are significantly more accurate, whereas those by Lacy (1997a) are about 2% higher. Both studies report radii in agreement with ours. Popper (1965Popper ( , 1980 reached masses and radii accurate to about 3% for VZ Hya, and within this uncertainty his results agree with ours. The less massive secondary component of WZ Oph appears slightly larger and more luminous but cooler than the primary. However, within the uncertainties the two components can be considered identical. On average, we obtain masses and radii for WZ Oph which are about 8% and 5% larger, respectively, than those given by Popper (1965), and they are far more accurate.
The uvby indices for the components of AD Boo, VZ Hya, and WZ Oph, included in Table 21, were calculated from the the combined indices of the systems at phase 0.25 (CVG08) and the luminosity ratios between their components (Tables 9, 14, 19). The E(b − y) interstellar reddenings, also given in Table 21, were determined from uvbyβ standard photometry and other sources, as described in more detail in Sect. 4.1.
The adopted effective temperatures listed in Table 21 were calculated from the calibration by Holmberg et al. (2007) Alonso et al. (1996) and the b − y calibration by Ramírez & Meléndez (2005) agree, but for AD Boo we note that both calibrations lead to a temperature difference between the components that is about 120 K higher. A similar trend was noticed by Holmberg et al. (2007). The empirical flux scale by Popper (1980) and the y flux ratios between the components (Tables 9, 14, 19) yield well established temperature differences of 450 K (AD Boo), 335 K (VZ Hya), and 20 K (WZ Oph), supporting the temperature calibration we have chosen. Details on temperatures derived as part of the TODCOR analyses of the observed CfA spectra (assuming the final metallicities) and the abundance analyses of the disentangled CfA spectra are given in Sect. 3. In general, the results agree well with those included in Table 21, except for WZ Oph where slightly higher temperatures were obtained from the analysis of the disentangled spectra.
The distances to the binaries were calculated from the "classical" relation, (1) + 10 log(T eff /T eff ⊙ ) + BC adopting the solar values and bolometric corrections given in Table 21, and A V /E(b − y) = 4.27. As seen, the distances have been established to 5-6 %, accounting for all error sources and including the use of other BC scales (e.g. Code et al. 1976, Bessell et al. 1998, Girardi et al. 2002. For AD Boo and VZ Hya, the empirical K surface brightness -T eff relation by Kervella et al. (2004) leads to nearly identical and perhaps even more precise distances; see Southworth et al. (2005) for details. K photometry outside eclipses is unfortunately not available for WZ Oph. VZ Hya and WZ Oph are included in the Hipparcos catalog (ESA 1997), which gives parallaxes of 5.03 ± 1.34 mas and 7.99 ± 1.37 mas, respectively, corresponding to distances of 199 ± 53 and 125 ± 21 pc. Matching the shorter Hipparcos distance for WZ Oph would require lower tem- Table 21. Astrophysical data for AD Boo, VZ Hya, and WZ Oph. T ef f ⊙ = 5780 K, M bol⊙ = 4.74, and bolometric corrections (BC) by Flower (1996) have been applied. v sync is the equatorial velocity for synchronous rotation. peratures, whereas higher interstellar absorption would imply higher temperatures and therefore a nearly unchanged distance. WZ Oph therefore does not formally belong to the group of eclipsing binaries within 125 pc discussed by Popper (1998b), which could be of use to improve the radiative flux scale, provided more precise astrometric parallaxes were available.

Interstellar reddening and absorption
The interstellar reddenings of AD Boo, VZ Hya, and WZ Oph have been determined as follows: For AD Boo, we adopt an interstellar reddening of E(b − y) = 0.025 ± 0.015, which is derived from the calibration by Olsen (1988) and the standard uvbyβ indices at phase 0.0, where a significant fraction of the primary component is eclipsed. Ideally, indices at phase 0.5, where the secondary component is totally eclipsed, should be used, but such β observations are unfortunately not available. Combining instead theoretical β values (Moon & Dworetsky 1985) and the individual uvby indices, reddenings of E(b − y) ≈ 0.01 ± 0.02 are obtained for both components. Alternatively, 2MASS JHK s photometry for stars within 1 • of AD Boo on the sky and closer than 200 pc show a mean absorption of A V = 0.18±0.04, corresponding to E(b − y) = 0.04 ± 0.01 (Knude, private communication). Finally, the model by Hakkila et al. (1997) yields A V = 0.10 or E(b − y) = 0.02, in the direction of and at the distance of AD Boo, and the Schlegel et al. (1998) dust maps give the same amount of total extinction.
For VZ Hya, the calibration by Olsen and the standard uvbyβ indices at phase 0.5, where the secondary is nearly totally eclipsed, lead to negligible interstellar reddening, as do the indices at phase 0.0, where a significant fraction of the primary component is eclipsed. On the other hand, 2MASS JHK s photometry for stars within 30 ′ of VZ Hya on the sky and closer than 150 pc show a mean absorption of A V = 0.17 ± 0.06, corresponding to E(b − y) = 0.040 ± 0.014. Furthermore, the model by Hakkila et al. also yields A V = 0.17 in the direction of and at the distance of VZ Hya, whereas the Schlegel et al. dust maps give a total extinction of A V = 0.12, corresponding to E(b − y) = 0.028. To encompass the different results we have adopted E(b − y) = 0.02 ± 0.02.
For WZ Oph we adopt an interstellar reddening of E(b − y) = 0.033 ± 0.012, which is derived from the calibration by Olsen and the standard uvbyβ indices outside of eclipse, which represent the nearly identical components well. 2MASS JHK s photometry for stars within 30 ′ of WZ Oph on the sky and closer than 180 pc show a mean absorption of A V = 0.10 ± 0.04, corresponding to E(b − y) = 0.024 ± 0.009. In the direction and at the distance of WZ Oph, the model of Hakkila et al. yields a very high absorption of A V = 0.52, whereas the Schlegel et al. dust maps give of total extinction of A V = 0.28 or E(b − y) = 0.066.

Discussion
In the following, we compare the new dimensions obtained for AD Boo, VZ Hya, and WZ Oph with properties predicted from some of the latest stellar evolutionary models. Below, we first briefly present the key ingredients of the selected grids, referring to the original papers for full descriptions. The confrontation between observations and theory for these three systems can be constrained better than in most previous studies, first of all because our [Fe/H] determinations and abundance results for other heavy elements are probably among the most detailed done to date for binaries. This means that we can select models not only with respect to mass, but also with respect to Z. On the other hand, for a given mass and Z, the observable properties of models at a given age such as radius, effective temperature, and luminosity still depend on the adopted input physics, including the treatment of core and envelope convection, diffusion, etc., and on the assumed element mixture and Y ,Z relation as well. In order to constrain such free model parameters, accurate results from many binaries are required, supplemented by other results such as those from asteroseismology and cluster analysis. We will focus on the three binaries studied here; a larger sample will be dealt with in forthcoming papers, cf. Fig. 1.
From the binary perspective, masses and radii are the most direct parameters available, free of any scale dependent calibrations. The M − R plane is therefore well suited for isochrone tests, especially when the binary components are significantly different. In addition, the T eff − R or T eff − log(g) planes allow tests of model temperatures, which for a given mass and Z depend, e.g., on the abundance mixture, Y , core overshoot treatment, and surface convection efficiency. In our case, the binary temperatures are based on several calibrations and methods, and furthermore the temperature difference between the components is well constrained from the light curves; cf. Sect. 2.5. Adopting the usual assumption of coeval formation from the same raw material, identical ages must be reached in these planes, as well as in the M − log(L) plane.

Summary of models
As a reference, we have adopted the Y 2 (Yonsei-Yale) evolutionary tracks and isochrones by Demarque et al. (2004) 9 , which are based on up-to-date input physics and include an improved core overshoot treatment where Λ OS = λ ov /H p depends on mass, and which also takes into account the composition dependence of M conv crit 10 The mixing length parameter in convective envelopes is calibrated using the Sun, and is held fixed at l/H p = 1.7432. Helium diffusion has been included. All tracks are started at the pre-mainsequence birth line. Tracks and isochrones are available for a broad range of masses, ages, and heavy element abundances. The enrichment law Y = 0.23 + 2Z is adopted, together with the solar mixture by Grevesse et al. (1996), leading to (X,Y ,Z) ⊙ = (0.71564,0.26624,0.01812). Models with enhanced α-element abundances are also included; see Kim et al. (2002) for further details. We have used the abundance, mass, and age interpolation routines provided by the Y 2 group.
Another recent example of stellar evolution calculations is given by the extensive Victoria-Regina grids 9 http://www.astro.yale.edu/demarque/yystar.html 10 Defined as "the mass above which stars continue to have a substantial convective core even after the end of the pre-MS phase." (VandenBerg et al., 2006) 11 , which adopt a somewhat different physically based description of the core overshoot with F over depending on mass and metallicity and calibrated observationally via cluster CMD's. Diffusive processes are not treated. All tracks have been computed from the pre-main-sequence stage. The VRSS grids (scaled-solar) used here adopt Y = 0.23544 + 2.2Z and the solar mixture by Grevesse & Noels (1993), leading to (X,Y ,Z) ⊙ = (0.7044,0.2768,0.0188). The mixing length parameter in convective envelopes is calibrated using the Sun and held at α MLT = 1.90. For further details on the chemical abundances adopted for the VR0A, VR2A, and VR4A grids (different [α/Fe]) we refer the reader to VandenBerg et al. (2000). We have used the isochrone interpolation routines provided by the Victoria-Regina group; track interpolations have not been done.
A third source of modern evolutionary tracks and isochrones for specific masses and ages is the BaSTI database (Pietrinferni et al., 2004) 12 , available for several compositions. It is based on a recent version of the FRANEC evolutionary code (Degl'Innocenti et al., 2007). Models without and with core overshoot are available; in the latter case the choice of λ ov depends on mass but not on metallicity. Y = 0.245 + 1.4Z and the solar mixture by Grevesse & Noels (1993) are used, and the mixing length parameter in convective envelopes is calibrated using the Sun and held fixed at l/H p = 1.913. Mass loss is included, but atomic diffusion is taken into account only in the standard solar model. All tracks have been computed from the pre-main-sequence phase.
A comparison of the main sequence models between 1.1 and 1.5 M ⊙ for identical [Fe/H] values close to those of our stars reveals some differences between the three grids described above. This is exemplified in Figs. 18 and 19. In the lower main-sequence band, the Y 2 and the Victoria-Regina tracks agree quite well. The isochrones are slightly shifted, and for part of the mass interval the effect increases in the upper part of the main-sequence band due to differences in the treatment of core overshoot and the ramping procedure adopted. The BaSTI tracks differ somewhat in shape and location, and consequently significant differences are seen in the isochrones as well.
Before focusing on the confrontation between binary observations and models, we point out that we have not considered here models based on scaling of the most recent solar mixture (Grevesse et al. 2007), which implies a significantly lower metallicity of Z ⊙ ≈ 0.013 and a different conversion between [Fe/H] and Z. At present, solar models computed with these abundances disagree with precise results inferred from helioseismology, and the cause(s) have yet to be identified (e.g. Basu & Antia 2007).

Comparison with the binary results
For AD Boo, comparisons between the absolute dimensions and Y 2 models calculated for the measured [Fe/H] are shown in Figs. 20-22. As seen in Fig. 21, an isochrone corresponding to ∼1.75 Gyr would fit both components very well, and for the primary component this age is also deduced from Figs. 20    ±0.15 dex translates to an age uncertainty of about ±0.15 Gyr, comparable to that from mass and radius uncertainties. The model tracks fit the components well but predict a ∼ 100 K lower temperature difference between the components than established from both photometry and the empirical flux scale; see Sect. 4. Consequently, the predicted model luminosity of the secondary is also marginally higher than observed (Fig. 22). If real, the slight discrepancy may be related to the activity of the secondary component. The same picture is seen for Victoria-Regina models, which give a slightly lower age of ∼1.5 Gyr. On the other hand, the BaSTI models fit AD Boo less well. For the observed masses and radii, they predict quite different ages of about 2.25 and 1.75 Gyr for the primary and secondary, respectively. This is probably related to the different shape of the BaSTI tracks ( Fig. 19) in the lower main-sequence region.
Figs. 23-25 compare the absolute dimensions for VZ Hya with the Y 2 models calculated for the measured [Fe/H]. As seen in Fig. 24, an isochrone for an age of 1.25 Gyr fits both components within errors, although a slightly steeper isochrone slope would formally improve the fit still more. The abundance uncertainty of ±0.12 dex corresponds to about ±0.25 Gyr in age, similar to the contribution from mass and radius uncertainties. The model tracks are slightly hotter than observed, and consequently the ages derived from Figs. 23 and 25 become too high and low, respectively. Both tracks would fit well for a slightly higher [Fe/H] of −0.15, safely within the uncertainty of the measured value. Adopting [Fe/H] = −0.15, identical ages of 1.25 Gyr are obtained from all three diagrams. Victoria-Regina tracks for the observed [Fe/H] fit the components well, being slightly cooler than the Y 2 tracks. For the observed masses and radii, ages become lower, about 1.0 and 0.7 Gyr for primary and secondary, respectively; they agree within errors. The BaSTI models yield about 1.25 and 0.75 Gyr, and in this case the difference is significant.
Whereas the Y 2 models in general match the observed properties of AD Boo and VZ Hya well, the situation for  Alternatively, the models would need to have a significantly lower helium content and/or an envelope mixing length parameter well below l/H p = 1.7432. The latter can not a priori be ruled out, since the photometry and spectroscopy indicate the presence of some degree of surface activity  Baraffe et al. (1998) 13 but only for solar metallicity. They suggest that l/H p ∼ 1 or even lower would be needed, although models computed specifically for this case would be required in order to evaluate the combined effects on the stellar radii and temperatures in detail. The possible connection between activity and convection, as indicated from studies of solar-type binaries, has been discussed by Clausen et al. (1999), Torres et al. (2006), and others.
As stated in Sect. 4.1, we have adopted an interstellar reddening of E(b − y) = 0.033 ± 0.012 for WZ Oph, equivalent to A V = 0.14 ± 0.05 for a standard extinction law. There were, however, some indications of higher absorption from other sources. If we therefore somewhat arbitrarily double it and assume A V = 0.28 and E(b − y) = 0.066, the effective temperatures derived from the uvby indices increase by about 250 K, i.e. to about 6400 K. This has an effect on the abundances. From the Fe I lines we would derive [Fe/H] ∼ −0.13 for both components, whereas the results from the Fe II lines would be practically unchanged at −0.27. The photometric determinations from the calibration by Holmberg et al. (2007) would increase by about 0.12 dex for both components. The discrepancy between neutral and ionized Fe lines argues against the higher temperatures, but assuming them, together with the higher abundance from the Fe I lines, does on the other hand re-  move the model discrepancies, as seen in Fig. 27. The age of WZ Oph then becomes about 3 Gyr.
An independent check would be to derive intrinsic colours and/or effective temperatures from line-depth ratios measured in the disentangled FEROS spectra, see e.g. Gray (1994), but WZ Oph is unfortunately slightly too hot for the calibrations available. Adopting a subset of those by Kovtyukh et al. (2003), valid to 6100 or 6150 K, we derive mean effective temperatures close to 6200 K for both components. The individual results from the 9 (primary) and 16 (secondary) measured line-depth ratios scatter by  Although we cannot pinpoint the reason for the discrepancy between the observations and the models for WZ Oph, we suspect they are connected to composition (He) or con-vection anomalies in the stars and/or the determination of interstellar reddening, rather than to physical or numerical problems with the model calculations. This is supported by the general good agreement seen for AD Boo and especially VZ Hya, but detailed studies of additional systems similar to these are needed for full clarification. Work on several new F-G type binaries is in progress.

Summary and conclusions
We have presented precise absolute dimensions and abundances for the three F-type double-lined eclipsing binaries AD Boo, VZ Hya, and WZ Oph. The results, which are based on state-of-the-art analyses of new photometric and spectroscopic observations, were used to test the recent Y 2 , Victoria-Regina, and BaSTI stellar evolutionary models.
Masses and radii precise to 0.5-0.7% and 0.4-0.9% have been established for the binary components. This level can only be reached if excellent observations are available and are analysed carefully using adequate tools. Special care was taken to identify and avoid possible systematic errors. Several methods were applied in order to determine reliable and accurate (100-150 K) effective temperatures and interstellar extinctions, still a challenging task whose importance can not be overestimated.
System [Fe/H] abundances were based on numerous lines and were derived relative to the Sun. We have reached precisions between 0.07 and 0.15 dex. Even having many Fe I and Fe II lines of different strength and excitation potential available, the strong correlation between abundance and effective temperature is difficult to break. T eff and microturbulence uncertainties contribute significantly to the final abundance precision. Abundances for other heavy elements, based on fewer lines, were also derived; in one case (WZ Oph) indications of a slight α-element overabundance is seen. Comparing the results for the three binaries clearly demonstrates that analyses of disentangled spectra are far superior to those of composite spectra, essentially because many more lines become available and blending becomes less of an issue.
AD Boo has components with the greatest difference in mass and radius among well-studied F-type binaries (1.41 M ⊙ , 1.61 R ⊙ ; 1.21 M ⊙ , 1.22 R ⊙ ) and an abundance of [Fe/H] = +0.10. Slight surface activity is present on the secondary component. The system is fitted well by Y 2 and Victoria-Regina evolutionary models for identical component ages (as determined from masses and radii) of 1.75 and 1.5 Gyr, respectively. The fact that the secondary appears marginally cooler than the corresponding evolutionary track may be related to mild surface activity. BaSTI models predict different ages of 2.25 and 1.75 Gyr for the components, and therefore do not fit AD Boo as well.
VZ Hya has somewhat different components (1.27 M ⊙ , 1.31 R ⊙ ; 1.14 M ⊙ , 1.11 R ⊙ ) and an abundance of [Fe/H] = −0.20. No signs of surface activity are seen. Y 2 and Victoria-Regina evolutionary tracks fit VZ Hya well. Nearly identical component ages (as determined from masses and radii) of 1.25 Gyr are obtained from the Y 2 isochrones, whereas the Victoria-Regina calculations give 1.0 and 0.7 Gyr for the primary and secondary. An even larger difference, 1.25 versus 0.75 Gyr, is found when using the BaSTI models.
WZ Oph consists of two nearly identical components (1.22 M ⊙ , 1.41 R ⊙ ). The [Fe/H] abundance is −0.27, and as mentioned above a slight α-element overabundance is possible. Photometry and spectroscopy suggest low-level surface activity on both components. Contrary to the case for AD Boo and VZ Hya, neither of the three model grids are able to fit the observations for WZ Oph; models are significantly hotter than observed. We tentatively conclude that this is caused either by WZ Oph anomalies such as a low He content or decreased envelope convection, or by underestimated interstellar reddening. Further observations will be required to resolve the issue.
We believe our study is among the most detailed carried out to date for main-sequence eclipsing binaries, and that it sets new standards for critical tests of stellar evolutionary models based on accurate binary data. With only three systems available, it is premature to draw firm conclusions and suggest specific model shortcomings and/or improvements. We find, however, better agreement with the Y 2 and Victoria-Regina models than with the BaSTI models. In contradiction to this, Tomasella et al. (2007) recently found that dedicated BaSTI models, including core overshoot and both helium and heavy element diffusion, represent the components of the young ∼ 1.25M ⊙ system V505 Per well. Their study is based on new absolute dimensions of extraordinary high formal precision.
As shown in Fig. 1, several additional F-type binaries are available, and we are presently working on re-analyses for some of them, including abundance determinations. We are also conducting a larger program on new F-G type systems exhibiting various levels of surface activity (Clausen et al. 2001), and we expect some very useful insights to come from these studies, and from a parallel program on eclipsing binaries in clusters spanning a wide range in age and metallicity (e.g. Grundahl et al. 2006).