We present models of alpha Centauri A and B implementing an entropy calibration of the mixing-length parameter alpha_MLT, recently developed and successfully applied to the Sun (Spada et al. 2018, ApJ, 869, 135). In this technique the value of alpha_MLT in the 1D stellar evolution code is calibrated to match the adiabatic specific entropy derived from 3D radiation-hydrodynamics simulations of stellar convective envelopes, whose effective temperature, surface gravity, and metallicity are selected consistently along the evolutionary track. The customary treatment of convection in stellar evolution models relies on a constant, solar-calibrated alpha_MLT. There is, however, mounting evidence that this procedure does not reproduce the observed radii of cool stars satisfactorily. For instance, modelling alpha Cen A and B requires an ad-hoc tuning of alpha_MLT to distinct, non-solar values. The entropy-calibrated models of alpha Cen A and B reproduce their observed radii within 1% (or better) without externally adjusted parameters. The fit is of comparable quality to that of models with freely adjusted alpha_MLT for alpha Cen B (within 1 sigma), while it is less satisfactory for alpha Cen A (within ~ 2.5 sigma). This level of accuracy is consistent with the intrinsic uncertainties of the method. Our results demonstrate the capability of the entropy calibration method to produce stellar models with radii accurate within 1%. This is especially relevant in characterising exoplanet-host stars and their planetary systems accurately.

Main sequence, solar-like stars (M < 1.5 Msun) have outer convective envelopes that are sufficiently thick to affect significantly their overall structure. The radii of these stars, in particular, are sensitive to the details of inefficient, super-adiabatic convection occurring in their outermost layers. The standard treatment of convection in stellar evolution models, based on the Mixing-Length Theory (MLT), provides only a very approximate description of convection in the super-adiabatic regime. Moreover, it contains a free parameter, alpha_MLT, whose standard calibration is based on the Sun, and is routinely applied to other stars ignoring the differences in their global parameters (e.g., effective temperature, gravity, chemical composition) and previous evolutionary history. In this paper, we present a calibration of alpha_MLT based on three-dimensional radiation-hydrodynamics (3D RHD) simulations of convection. The value of alpha_MLT is adjusted to match the specific entropy in the deep, adiabatic layers of the convective envelope to the corresponding value obtained from the 3D RHD simulations, as a function of the position of the star in the (log g, log T_eff) plane and its chemical composition. We have constructed a model of the present-day Sun using such entropy-based calibration. We find that its past luminosity evolution is not affected by the entropy calibration. The predicted solar radius, however, exceeds that of the standard model during the past several billion years, resulting in a lower surface temperature. This illustrative calculation also demonstrates the viability of the entropy approach for calibrating the radii of other late-type stars.

Sirius, the seventh-nearest stellar system, is a visual binary containing the metallic-line A1 V star Sirius A, brightest star in the sky, orbited in a 50.13-year period by Sirius B, the brightest and nearest white dwarf (WD). Using images obtained over nearly two decades with the Hubble Space Telescope (HST), along with photographic observations covering almost 20 years, and nearly 2300 historical measurements dating back to the 19th century, we determine precise orbital elements for the visual binary. Combined with the parallax and the motion of the A component, these elements yield dynamical masses of 2.063+/-0.023 Msun and 1.018+/-0.011 Msun for Sirius A and B, respectively. Our precise HST astrometry rules out third bodies orbiting either star in the system, down to masses of ~15-25 Mjup. The location of Sirius B in the H-R diagram is in excellent agreement with theoretical cooling tracks for WDs of its dynamical mass, and implies a cooling age of ~126 Myr. The position of Sirius B in the mass-radius plane is also consistent with WD theory, assuming a carbon-oxygen core. Including the pre-WD evolutionary timescale of the assumed progenitor, the total age of Sirius B is about 228+/-10 Myr. We calculated evolutionary tracks for stars with the dynamical mass of Sirius A, using two independent codes. We find it necessary to assume a slightly sub-solar metallicity, of about 0.85 Zsun, to fit its location in the luminosity-radius plane. The age of Sirius A based on these models is about 237-247 Myr, with uncertainties of +/-15 Myr, consistent with that of the WD companion. We discuss astrophysical puzzles presented by the Sirius system, including the probability that the two stars must have interacted in the past, even though there is no direct evidence for this, and the orbital eccentricity remains high.

We introduce the Yale-Potsdam Stellar Isochrones (YaPSI), a new grid of stellar evolution tracks and isochrones of solar-scaled composition. In an effort to improve the Yonsei-Yale database, special emphasis is placed on the construction of accurate low-mass models (Mstar < 0.6 Msun), and in particular of their mass-luminosity and mass-radius relations, both crucial in characterizing exoplanet-host stars and, in turn, their planetary systems. The YaPSI models cover the mass range 0.15 to 5.0 Msun, densely enough to permit detailed interpolation in mass, and the metallicity and helium abundance ranges [Fe/H] = -1.5 to +0.3, and Y = 0.25 to 0.37, specified independently of each other (i.e., no fixed Delta Y/Delta Z relation is assumed). The evolutionary tracks are calculated from the pre-main sequence up to the tip of the red giant branch. The isochrones, with ages between 1 Myr and 20 Gyr, provide UBVRI colors in the Johnson-Cousins system, and JHK colors in the homogeneized Bessell & Brett system, derived from two different semi-empirical Teff-color calibrations from the literature. We also provide utility codes, such as an isochrone interpolator in age, metallicity, and helium content, and an interface of the tracks with an open-source Monte Carlo Markov-Chain tool for the analysis of individual stars. Finally, we present comparisons of the YaPSI models with the best empirical mass- luminosity and mass-radius relations available to date, as well as isochrone fitting of well-studied ste

One of the largest sources of uncertainty in stellar models is caused by the treatment of convection in stellar envelopes. One dimensional stellar models often make use of the mixing length or equivalent approximations to describe convection, all of which depend on various free parameters. There have been attempts to rectify this by using 3D radiative-hydrodynamic simulations of stellar convection, and in trying to extract an equivalent mixing length from the simulations. In this paper we show that the entropy of the deeper, adiabatic layers in these simulations can be expressed as a simple function of og g and log T_{eff} which holds potential for calibrating stellar models in a simple and more general manner.

The nearby star Procyon is a visual binary containing the F5 IV-V subgiant Procyon A, orbited in a 40.84 yr period by the faint DQZ white dwarf Procyon B. Using images obtained over two decades with the Hubble Space Telescope, and historical measurements back to the 19th century, we have determined precise orbital elements. Combined with measurements of the parallax and the motion of the A component, these elements yield dynamical masses of 1.478 +/- 0.012 Msun and 0.592 +/- 0.006 Msun for A and B, respectively. The mass of Procyon A agrees well with theoretical predictions based on asteroseismology and its temperature and luminosity. Use of a standard core-overshoot model agrees best for a surprisingly high amount of core overshoot. Under these modeling assumptions, Procyon A's age is ~2.7 Gyr. Procyon B's location in the H-R diagram is in excellent agreement with theoretical cooling tracks for white dwarfs of its dynamical mass. Its position in the mass-radius plane is also consistent with theory, assuming a carbon-oxygen core and a helium-dominated atmosphere. Its progenitor's mass was 1.9-2.2 Msun, depending on its amount of core overshoot. Several astrophysical puzzles remain. In the progenitor system, the stars at periastron were separated by only ~5 AU, which might have led to tidal interactions and even mass transfer; yet there is no direct evidence that these have occurred. Moreover the orbital eccentricity has remained high (~0.40). The mass of Procyon B is somewhat lower than anticipated from the initial-to-final-mass relation seen in open clusters. The presence of heavy elements in its atmosphere requires ongoing accretion, but the place of origin is uncertain.

In an effort to better understand the details of the stellar structure and evolution of metal poor stars, the Gemini North telescope was used on two occasions to take speckle imaging data of a sample of known spectroscopic binary stars and other nearby stars in order to search for and resolve close companions. The observations were obtained using the Differential Speckle Survey Instrument, which takes data in two filters simultaneously. The results presented here are of 90 observations of 23 systems in which one or more companions was detected, and 6 stars where no companion was detected to the limit of the camera capabilities at Gemini. In the case of the binary and multiple stars, these results are then further analyzed to make first orbit determinations in five cases, and orbit refinements in four other cases. Mass information is derived, and since the systems span a range in metallicity, a study is presented that compares our results with the expected trend in total mass as derived from the most recent Yale isochrones as a function of metal abundance. These data suggest that metal-poor main-sequence stars are less massive at a given color than their solar-metallicity analogues in a manner consistent with that predicted from the theory.

We present direct radii measurements of the well-known transiting exoplanet host stars HD 189733 and HD 209458 using the CHARA Array interferometer. We find the limb-darkened angular diameters to be theta_LD = 0.3848 +/- 0.0055 and 0.2254 +/- 0.0072 milliarcsec for HD 189733 and HD 209458, respectively. HD 189733 and HD 209458 are currently the only two transiting exoplanet systems where detection of the respective planetary companion's orbital motion from high resolution spectroscopy has revealed absolute masses for both star and planet. We use our new measurements together with the orbital information from radial velocity and photometric time series data, Hipparcos distances, and newly measured bolometric fluxes to determine the stellar effective temperatures (T_eff = 4875 +/- 43, 6093 +/- 103 K), stellar linear radii (R_* = 0.805 +/- 0.016, 1.203 +/- 0.061 R_sun), mean stellar densities (rho_* = 1.62 +/- 0.11, 0.58 +/- 0.14 rho_sun), planetary radii (R_p = 1.216 +/- 0.024, 1.451 +/- 0.074 R_Jup), and mean planetary densities (rho_p = 0.605 +/- 0.029, 0.196 +/- 0.033 rho_Jup) for HD 189733 b and HD 209458 b, respectively. The stellar parameters for HD 209458, a F9 dwarf, are consistent with indirect estimates derived from spectroscopic and evolutionary modeling. However, we find that models are unable to reproduce the observational results for the K2 dwarf, HD 189733. We show that, for stellar evolutionary models to match the observed stellar properties of HD 189733, adjustments lowering the solar-calibrated mixing length parameter from 1.83 to 1.34 need to be employed.

We compare evolved stellar models, which match Procyons mass and position in the HR diagram, to current ground-based asteroseismic observations. Diffusion of helium and metals along with two conventional core overshoot descriptions and the Kuhfuss nonlocal theory of convection are considered. We establish that one of the two published asteroseismic data reductions for Procyon, which mainly differ in their identification of even versus odd l-values, is a significantly more probable and self-consistent match to our models than the other. The most probable models according to our Bayesian analysis have evolved to just short of turnoff, still retaining a hydrogen convective core. Our most probable models include Y and Z diffusion and have conventional core overshoot between 0.9 and 1.5 pressure scale heights, which increases the outer radius of the convective core by between 22% to 28%, respectively. We discuss the significance of this comparatively higher than expected core overshoot amount in terms of internal mixing during evolution. The parameters of our most probable models are similar regardless of whether adiabatic or nonadiabatic model p-mode frequencies are compared to the observations, although, the Bayesian probabilities are greater when the nonadiabatic model frequencies are used. All the most probable models (with or without core overshoot, adiabatic or nonadiabatic model frequencies, diffusion or no diffusion, including priors for the observed HRD location and mass or not) have masses that are within one sigma of the observed mass 1.497+/-0.037 Msun.

Mixing length theory is the predominant treatment of convection in stellar models today. Usually described by a single free parameter, alpha, the common practice is to calibrate it using the properties of the Sun, and apply it to all other stellar models as well. Asteroseismic data from Kepler and CoRoT provide precise properties of other stars which can be used to determine alpha as well, and a recent study of stars in the Kepler field of view found alpha to vary with metallicity. Interpreting alpha obtained from calibrated stellar models, however, is complicated by the fact that the value for alpha depends on the surface boundary condition of the stellar model, or T-tau relation. Calibrated models that use typical T-tau relations, which are static and insensitive to chemical composition, do not include the complete effect of metallicity on alpha. We use 3D radiation-hydrodynamic simulations to extract metallicity-dependent T-tau relations and use them in calibrated stellar models. We find the previously reported alpha-metallicity trend to be robust, and not significantly affected by the surface boundary condition of the stellar models.

We investigate the effect of helium abundance and $\alpha$-element enhancement on the properties of convection in envelopes of solar-like main-sequence stars stars using a grid of 3D radiation hydrodynamic simulations. Helium abundance increases the mean molecular weight of the gas, and alters opacity by displacing hydrogen. Since the scale of the effect of helium may depend on the metallicity, the grid consists of simulations with three helium abundances ($Y=0.1, 0.2, 0.3$), each with two metallicities ($Z=0.001, 0.020)$. We find that changing the helium mass fraction generally affects structure and convective dynamics in a way opposite to that of metallicity. Furthermore, the effect is considerably smaller than that of metallicity. The signature of helium differs from that of metallicity in the manner in which the photospheric velocity distribution is affected. \rev{We also find that helium abundance and surface gravity behave largely in similar ways, but differ in the way they affect the mean molecular weight}. A simple model for spectral line formation suggests that the bisectors and absolute Doppler shifts of spectral lines depends on the helium abundance. We look at the effect of $\alpha$-element enhancement and find that it has a considerably smaller effect on the convective dynamics in the SAL compared to that of helium abundance.

A long-standing issue in the theory of low mass stars is the discrepancy between predicted and observed radii and effective temperatures. In spite of the increasing availability of very precise radius determinations from eclipsing binaries and interferometric measurements of radii of single stars, there is no unanimous consensus on the extent (or even the existence) of the discrepancy and on its connection with other stellar properties (e.g. metallicity, magnetic activity). We investigate the radius discrepancy phenomenon using the best data currently available (accuracy about 5%). We have constructed a grid of stellar models covering the entire range of low mass stars (0.1-1.25 M_sun) and various choices of the metallicity and of the mixing length parameter \alpha. We used an improved version of the Yale Rotational stellar Evolution Code (YREC), implementing surface boundary conditions based on the most up-to-date PHOENIX atmosphere models. Our models are in good agreement with others in the literature and improve and extend the low mass end of the Yale-Yonsei isochrones. Our calculations include rotation-related quantities, such as moments of inertia and convective turnover time scales, useful in studies of magnetic activity and rotational evolution of solar-like stars. Consistently with previous works, we find that both binaries and single stars have radii inflated by about 3% with respect to the theoretical models; among binaries, the components of short orbital period systems are found to be the most deviant. We conclude that both binaries and single stars are comparably affected by the radius discrepancy phenomenon.

We examine how metallicity affects convection and overshoot in the superadiabatic layer of main sequence stars. We present results from a grid of 3D radiation hydrodynamic simulations with four metallicities ($Z=0.040$, 0.020, 0.010, 0.001), and spanning a range in effective temperature ($4950 < \rm T_{\rm eff} < 6230$). We show that changing the metallicity alters properties of the convective gas dynamics, and the structure of the superadiabatic layer and atmosphere. Our grid of simulations show that the amount of superadiabaticity, which tracks the transition from efficient to inefficient convection, \rev{is sensitive to changes in metallicity. We find that increasing the metallicity forces the location of the transition region to lower densities and pressures, and results in larger mean and turbulent velocities throughout the superadiabatic region.} We also quantify the degree of convective overshoot in the atmosphere, and show that it increases with metallicity as well.

We examine the effect of different radiative transfer schemes on the properties of 3D simulations of near-surface stellar convection in the superadiabatic layer, where energy transport transitions from fully convective to fully radiative. We employ two radiative transfer schemes that fundamentally differ in the way they cover the 3D domain. The first solver approximates domain coverage with moments, while the second solver samples the 3D domain with ray integrations. By comparing simulations that differ only in their respective radiative transfer methods, we are able to isolate the effect that radiative efficiency has on the structure of the superadiabatic layer. We find the simulations to be in good general agreement, but they show distinct differences in the thermal structure in the superadiabatic layer and atmosphere.

We present models of the components of the systems KOI-126 and CM Draconis, the two eclipsing binary systems known to date to contain stars with masses low enough to have fully convective interiors. We are able to model satisfactorily the system KOI-126, finding consistent solutions for the radii and surface temperatures of all three components, using a solar-like value of the mixing-length parameter \alpha in the convection zone, and PHOENIX NextGen 1D model atmospheres for the surface boundary conditions. Depending on the chemical composition, we estimate the age of the system to be in the range 3-5 Gyr. For CM Draconis, on the other hand, we cannot reconcile our models with the observed radii and T_eff using the current metal-poor composition estimate based on kinematics. Higher metallicities lessen but do not remove the discrepancy. We then explore the effect of varying the mixing length parameter \alpha. As previously noted in the literature, a reduced \alpha can be used as a simple measure of the lower convective efficiency due to rotation and induced magnetic fields. Our models show a sensitivity to \alpha (for \alpha < 1.0) sufficient to partially account for the radius discrepancies. It is, however, impossible to reconcile the models with the observations on the basis of the effect of the reduced \alpha alone. We therefore suggest that the combined effects of high metallicity and \alpha reduction could explain the observations of CM Draconis. For example, increasing the metallicity of the system towards super-solar values (i.e. Z = 2 Z_sun) yields an agreement within 2 \sigma with \alpha = 1.0.

We have analysed oscillations of the red giant star HD 186355 observed by the NASA Kepler satellite. The data consist of the first five quarters of science operations of Kepler, which cover about 13 months. The high-precision time-series data allow us to accurately extract the oscillation frequencies from the power spectrum. We find the frequency of the maximum oscillation power, {\nu}_max, and the mean large frequency separation, {\Delta}{\nu}, are around 106 and 9.4 {\mu}Hz respectively. A regular pattern of radial and non-radial oscillation modes is identified by stacking the power spectra in an echelle diagram. We use the scaling relations of {\Delta}{\nu} and {\nu}_max to estimate the preliminary asteroseismic mass, which is confirmed with the modelling result (M = 1.45 \pm 0.05 M_sun) using the Yale Rotating stellar Evolution Code (YREC7). In addition, we constrain the effective temperature, luminosity and radius from comparisons between observational constraints and models. A number of mixed l = 1 modes are also detected and taken into account in our model comparisons. We find a mean observational period spacing for these mixed modes of about 58 s, suggesting that this red giant branch star is in the shell hydrogen-burning phase.

We present a brief overview of the history of attempts to obtain a clear detection of solar-like oscillations in cluster stars, and discuss the results on the first clear detection, which was made by the Kepler Asteroseismic Science Consortium (KASC) Working Group 2.

Asteroseismology of stars in clusters has been a long-sought goal because the assumption of a common age, distance and initial chemical composition allows strong tests of the theory of stellar evolution. We report results from the first 34 days of science data from the Kepler Mission for the open cluster NGC 6819 -- one of four clusters in the field of view. We obtain the first clear detections of solar-like oscillations in the cluster red giants and are able to measure the large frequency separation and the frequency of maximum oscillation power. We find that the asteroseismic parameters allow us to test cluster-membership of the stars, and even with the limited seismic data in hand, we can already identify four possible non-members despite their having a better than 80% membership probability from radial velocity measurements. We are also able to determine the oscillation amplitudes for stars that span about two orders of magnitude in luminosity and find good agreement with the prediction that oscillation amplitudes scale as the luminosity to the power of 0.7. These early results demonstrate the unique potential of asteroseismology of the stellar clusters observed by Kepler.

Understanding the reasons of the cyclic variation of the total solar irradiance is one of the most challenging targets of modern astrophysics. These studies prove to be essential also for a more climatologic issue, associated to the global warming. Any attempt to determine the solar components of this phenomenon must include the effects of the magnetic field, whose strength and shape in the solar interior are far from being completely known. Modelling the presence and the effects of a magnetic field requires a 2D approach, since the assumption of radial symmetry is too limiting for this topic. We present the structure of a 2D evolution code that was purposely designed for this scope; rotation, magnetic field and turbulence can be taken into account. Some preliminary results are presented and commented.

We extend the analysis of Penev et al. (2007) to calculate effective viscosities for the surface convective zones of three main sequence stars of 0.775Msun, 0.85Msun and the present day Sun. In addition we also pay careful attention to all normalization factors and assumptions in order to derive actual numerical prescriptions for the effective viscosity as a function of the period and direction of the external shear. Our results are applicable for periods that are too long to correspond to eddies that fall within the inertial subrange of Kolmogorov scaling, but no larger than the convective turnover time, when the assumptions of the calculation break down. We find linear scaling of effective viscosity with period and magnitudes at least three times larger than the Zahn (1966, 1989) prescription.

The GIII red giant star epsilon Oph has been found to exhibit several modes of oscillation by the MOST mission. We interpret the observed frequencies of oscillation in terms of theoretical radial p-mode frequencies of stellar models. Evolutionary models of this star, in both shell H-burning and core He-burning phases of evolution, are constructed using as constraints a combination of measurements from classical ground-based observations (for luminosity, temperature, and chemical composition) and seismic observations from MOST. Radial frequencies of models in either evolutionary phase can reproduce the observed frequency spectrum of epsilon Oph almost equally well. The best-fit models indicate a mass in the range of 1.85 +/- 0.05 Msun with radius of 10.55 +/- 0.15 Rsun. We also obtain an independent estimate of the radius of epsilon Oph using high accuracy interferometric observations in the infrared K' band, using the CHARA/FLUOR instrument. The measured limb darkened disk angular diameter of epsilon Oph is 2.961 +/- 0.007 mas. Together with the Hipparcos parallax, this translates into a photospheric radius of 10.39 +/- 0.07 Rsun. The radius obtained from the asteroseismic analysis matches the interferometric value quite closely even though the radius was not constrained during the modelling.

In the second paper of this series we pursue two objectives. First, in order to make the code more sensitive to small effects, we remove many approximations made in Paper I. Second, we include turbulence and rotation in the two-dimensional framework. The stellar equilibrium is described by means of a set of five differential equations, with the introduction of a new dependent variable, namely the perturbation to the radial gravity, that is found when the non-radial effects are considered in the solution of the Poisson equation; following the scheme of the first paper, we write the equations in such a way that the two-dimensional effects can be easily disentangled. The key concept introduced in this series is the equipotential surface. We use the underlying cause-effect relation to develop a recurrence relation to calculate the equipotential surface functions for uniform rotation, differential rotation, rotation-like toroidal magnetic fields and turbulence. We also develop a more precise code to numerically solve the two-dimensional stellar structure and evolution equations based on the equipotential surface calculations. We have shown that with this formulation we can achieve the precision required by observations by appropriately selecting the convergence criterion. Several examples are presented to show that the method works well. Since we are interested in modeling the effects of a dynamo-type field on the detailed envelope structure and global properties of the Sun, the code has been optimized for short timescales phenomena (down to 1 yr). The time dependence of the code has so far been tested exclusively to address such problems.

We compute the p-mode oscillation frequencies and frequency splittings that arise in a two-dimensional model of the Sun that contains toroidal magnetic fields in its interior.

We report on a new, wide field ($20^{\prime} \times 20^{\prime}$), multicolor ($UBVI$), photometric campaign in the area of the nearby old open cluster NGC 2112. At the same time, we provide medium-resolution spectroscopy of 35 (and high-resolution of additional 5) Red Giant and Turn Off stars. This material is analyzed with the aim to update the fundamental parameters of this traditionally difficult cluster, which is very sparse and suffers from heavy field star contamination. Among the 40 stars with spectra, we identified 21 {\it bona fide} radial velocity members which allow us to put more solid constraints on the cluster's metal abundance, long suggested to be as low as the metallicity of globulars. As indicated earlier by us on a purely photometric basis (Carraro et al. 2002), the cluster [Fe/H] abundance is slightly super-solar ([Fe/H] =0.16$\pm$0.03) and close to the Hyades value, as inferred from a detailed abundance analysis of 3 of the 5 stars with higher resolution spectra. Abundance ratios are also marginally super solar. Based on this result, we revise the properties of NGC 2112 using stellar models from the Padova and Yale-Yonsei groups. For this metal abundance, we find the cluster's age, reddening, and distance values are 1.8 Gyr, 0.60 mag, and 940 pc, respectively. Both the Yale-Yonsei and Padova models predict the same values for the fundamental parameters within the errors. Overall, NGC 2112 is a typical solar neighborhood, thin disk star cluster, sharing the same chemical properties of F-G stars and open clusters close to the Sun. This investigation outlines the importance of a detailed membership analysis in the study of disk star clusters.

We describe a web-based cgi calculator for constructing synthetic color-magnitude diagrams for a simple stellar population (SSP) using the Yonsei-Yale (YY) isochrone data base. This calculator is designed to be used interactively. It creates quick look CMD displays in (B-V) and (V-I) colors. Stochastic effects on the CMDs are included. Output in tabular form is also provided for special purpose displays, or for combining the CMDs of different stellar populations. This research tool has applications in studies of the stellar content of our Galaxy and external systems. It provides an easy way to interpret the CMDs in resolved stellar populations. It offers the means to explore the dependence of the integrated properties of unresolved stellar systems on stellar parameters (ages, chemical composition, binarity) and on the characteristics of their parent population (IMF slope and mass range).

The stellar evolution code YREC is outlined with emphasis on its applications to helio- and asteroseismology. The procedure for calculating calibrated solar and stellar models is described. Other features of the code such as a non-local treatment of convective core overshoot, and the implementation of a parametrized description of turbulence in stellar models, are considered in some detail. The code has been extensively used for other astrophysical applications, some of which are briefly mentioned at the end of the paper.

We present preliminary results concerning the search for short-period variable stars in Tucana and LGS3 based on very deep HST/ACS imaging. In the fraction of the observed field we studied in each galaxy, a total of 133 and 30 variables were found, respectively. For Tucana, we identified 76 of them as RR Lyrae (RRL) stars pulsating in the fundamental mode (RRab) and 32 in the first-overtone mode (RRc), as well as 2 anomalous Cepheids (AC). The mean period of the RRab and RRc is 0.59 and 0.35 days, respectively. In the case of LGS3, we found 24 RRab and 4 RRc, with mean periods 0.61and 0.39 days, respectively, plus two candidate ACs. These values place both galaxies in the Oosterhoff gap.

AIMS: The cluster of hot stars observed in orbit around the central black hole of M31 has been interpreted as a 200 Myr starburst. The formation of a population of young stars in close proximity to a massive black hole presents a difficult challenge to star formation theory. We point out that in a high stellar density environment, the course of stellar evolution is modified by frequent collisions and mergers. METHODS: Blue stragglers, which are the results of mergers in globular clusters, occupy the same position in the color-magnitude diagram as the observed hot stars in M31. For confirmation, the integrated spectrum of P3 is shown to be compatible with the spectral energy distribution of a blue horizontal branch field star. RESULTS: We suggest an old stellar population of evolved blue horizontal-branch stars and of merger products cannot be ruled out on the basis of the available data. Observations are suggested that would help distinguish between a ``young'' and ``old'' stellar population interpretation of the observations.

We present a detailed study of the small frequency separations as diagnostics of the mass of the convective core and evolutionary stage of solar-type stars. We demonstrate how the small separations can be combined to provide sensitive tests for the presence of convective overshoot at the edge of the core. These studies are focused on low degree oscillation modes, the only modes expected to be detected in distant stars. Using simulated data with realistic errors, we find that the mass of the convective core can be estimated to within 5% if the total stellar mass is known. Systematic errors arising due to uncertainty in the mass could be up to 20%. The evolutionary stage of the star, determined in terms of the central hydrogen abundance using our proposed technique, however, is much less sensitive to the mass estimate.

The development of 2D and 3D simulations of solar convection has lead to a picture of convection quite unlike the usually assumed Kolmogorov spectrum turbulent flow. We investigate the impact of this changed structure on the dissipation properties of the convection zone, parametrized by an effective viscosity coefficient. We use an expansion treatment developed by Goodman & Oh 1997, applied to a numerical model of solar convection (Robinson et al. 2003) to calculate an effective viscosity as a function of frequency and compare this to currently existing prescriptions based on the assumption of Kolmogorov turbulence (Zahn 1966, Goldreich & Keeley 1977). The results match quite closely a linear scaling with period, even though this same formalism applied to a Kolmogorov spectrum of eddies gives a scaling with power-law index of 5/3.

A high-precision two-dimensional stellar evolution code has been developed for studying solar variability due to structural changes produced by varying internal magnetic fields of arbitrary configurations. Specifically, we are interested in modeling the effects of a dynamo-type field on the detailed internal structure and on the global parameters of the Sun. The high precision is required both to model very small solar changes (of order of $10^{-4}$) and short time scales (or order of one year). It is accomplished by using the mass coordinate to replace the radial coordinate, by using fixed and adjustable time steps, a realistic stellar atmosphere, elements diffusion, and by adjusting the grid points. We have also built into the code the potential to subsequently include rotation and turbulence. The current code has been tested for several cases, including its ability to reproduce the 1-D results.

We report on high resolution Echelle spectroscopy of 20 giant stars in the Galactic old open clusters NGC 6791 obtained with Hydra at the WIYN telescope. High precision radial velocity allow us to isolate 15 {\it bona fide} cluster members. From 10 of them we derive a global [M/H]=+0.39$\pm$0.05. We therefore confirm that NGC 6791 is extremely metal rich, exhibits a few marginally sub-solar abundance ratios, and within the resolution of our spectra does not show evidences of spread in metal abundance. With these new data we re-derive the cluster fundamental parameters suggesting that it is about 8 Gyr old and 4.3 kpc far from the Sun. The combination of its chemical properties, age, position, and Galactic orbit hardly makes NGC 6791 a genuine Population I open cluster. We discuss possible interpretations of the cluster peculiarities suggesting that the cluster might be what remains of a much larger system, whose initial potential well could have been sufficient to produce high metallicity stars, and which has been depopulated by the tidal field of the Galaxy. Alternatively, its current properties may be explained by the perturbation of the Galactic bar on an object originated well inside the solar ring, where the metal enrichment had been very fast.

We show in this paper that the changes of the solar diameter in response to variations of large scale magnetic fields and turbulence are not homologous. For the best current model, the variation at the photospheric level is over 1000 times larger than the variation at a depth of 5 Mm, which is about the level at which f-mode solar oscillations determine diameter variations. This model is supported by observations that indicate larger diameter changes for high degree f-modes than for low degree f-modes, since energy of the former are concentrated at shallower layers than the latter.

The space telescope MOST is now providing us with extremely accurate low frequency p-mode oscillation data for the star Eta Boo. We demonstrate in this paper that these data, when combined with ground based measurements of the high frequency p-mode spectrum, can be reproduced with stellar models that include the effects of turbulence in their outer layers. Without turbulence, the l=0 modes of our models deviate from either the ground based or the space data by about 1.5-4.0 micro Hz. This discrepancy can be completely removed by including turbulence in the models and we can exactly match 12 out of 13 MOST frequencies that we identified as l=0 modes in addition to 13 out of 21 ground based frequencies within their observational 2 sigma tolerances. The better agreement between model frequencies and observed ones depends for the most part on the turbulent kinetic energy which was taken from a 3D convection simulation for the Sun.

Compared to the Sun, the atmospheric structure and convective flow in Procyon A exhibit the following characteristics: (1) the highly superadiabatic transition layer (SAL) is located at much shallower optical depth; it is in a dynamically active region, and its outer region is located part of the time in the optically thin atmosphere; (2) the outer region of the SAL moves from an optically thin region to thick region and back again over a time of 20-30 minutes. This motion, which is driven by the granulation, takes place in a time approximately half the turnover time of the largest granules; The main reason for the radically different radiative-convective behaviour in Procyon A compared to the Sun is the role played by turbulent eddies in determining the overall flow/thermal structure. The turbulent pressure and turbulent kinetic energy can exceed 50 % of the local gas pressure (compared to about 10-20 % in the Sun). The Procyon A simulation thus reveals two distinct timescales - the autocorrelation time of the vertical velocity and the characteristic timescale of the SAL which is tied to granulation. Just below the surface the autocorrelation decay time is about 5 minutes in Procyon A, and the SAL motion timescale is 20-30 mins. When the SAL penetrates the optically thin region there are efficient radiative losses and the peak of the SAL is low. We speculate that these losses damp out the relative amplitudes in luminosity (temperature fluctuations) compared to velocity (Doppler). Although this will not affect the frequencies of the peaks in the power spectrum, it will probably lower the average amplitude of the peaks relative to the noise background.

We present a comprehensive set of stellar evolution models for Procyon A in an effort to guide future measurements of both traditional stellar parameters and seismic frequencies towards constraining the amount of core overshoot in Procyon A and possibly other stars. Current observational measurements of Procyon A when combined with traditional stellar modeling only place a large upper limit on overshoot of alphaOV < 1.1. By carrying out a detailed pulsation analysis, we further demonstrate, how p- and g-mode averaged spacings can be used to gain better estimates of the core size. For both p- and g-modes, the frequency spacings for models without overshoot are clearly separated from the models with overshoot. In addition, measurements of the l=0 averaged small p-mode spacings could be used to establish Procyon A's evolutionary stage. For a fixed implementation of overshoot and under favorable circumstances, the g-mode spacings can be used to determine the overshoot extent to an accuracy of +-0.05 Hp. However, we stress that considerable confusion is added due to the unknown treatment of the overshoot region. This ambiguity might be removed by analyzing many different stars. A simple non-local convection theory developed by Kuhfuss is implemented in our stellar evolution code and contrasted with the traditional approaches. We show that this theory supports a moderate increase of the amount of convective overshoot with stellar mass of Delta(alphaOV) = +0.10 between 1.5 Msun and 15 Msun. This theory places an upper limit on Procyon A's core overshoot extent of ~0.4 Hp which matches the limit imposed by Roxburgh's integral criterion.

Recent observations for the color-magnitude diagrams (CMDs) of the massive globular cluster Omega Centauri have shown that it has a striking double main sequence (MS), with a minority population of bluer and fainter MS well separated from a majority population of MS stars. Here we confirm, with the most up-to-date Y2 isochrones, that this special feature can only be reproduced by assuming a large variation (Delta Y = 0.15) of primordial helium abundance among several distinct populations in this cluster. We further show that the same helium enhancement required for this special feature on the MS can by itself reproduce the extreme horizontal-branch (HB) stars observed in Omega Cen, which are hotter than normal HB stars. Similarly, the complex features on the HBs of other globular clusters, such as NGC 2808, are explained by large internal variations of helium abundance. Supporting evidence for the helium-rich population is also provided by the far-UV (FUV) observations of extreme HB stars in these clusters, where the enhancement of helium can naturally explain the observed fainter FUV luminosity for these stars. The presence of super helium-rich populations in some globular clusters suggests that the third parameter, other than metallicity and age, also influences CMD morphology of these clusters.

Convective core overshoot affects stellar evolution rates and the dating of stellar populations. In this paper, we provide a patch to the $Y^2$ isochrones with an improved treatment of convective core overshoot. The new tracks cover the transition mass range from no convective core to a fully developed convective core. We compare the improved isochrones to CMDs of a few well observed open star clusters in the Galaxy and the Large Magellanic Cloud. Finally we discuss future prospects for improving the treatment of core overshoot with the help of asteroseismology.

Intermediate degree modes of the solar oscillations have previously been used to determine the solar helium abundance to a high degree of precision. However, we cannot expect to observe such modes in other stars. In this work we investigate whether low degree modes that should be available from space-based asteroseismology missions can be used to determine the helium abundance, Y, in stellar envelopes with sufficient precision. We find that the oscillatory signal in the frequencies caused by the depression in \Gamma_1 in the second helium ionisation zone can be used to determine the envelope helium abundance of low mass main sequence stars. For frequency errors of 1 part in 10^4, we expect errors \sigma_Y in the estimated helium abundance to range from 0.03 for 0.8M_sun stars to 0.01 for 1.2M_sun stars. The task is more complicated in evolved stars, such as subgiants, but is still feasible if the relative errors in the frequencies are less than 10^{-4}.

We present a preliminary theoretical seismic study of the astronomically famous star 51 Peg. This is done by first performing a detailed analysis within the Hertzsprung-Russell diagram (HRD). Using the Yale stellar evolution code (YREC), a grid of stellar evolutionary tracks has been constructed for the masses 1.00 M_sun, 1.05 M_sun and 1.10 M_sun, in the metallicity range Z=0.024-0.044, and for values of the Galactic helium enrichment ratio DY/DZ in the range 0-2.5. Along these evolutionary tracks, we select 75 stellar model candidates that fall within the 51 Peg observational error box in the HRD (all turn out to have masses of 1.05 M_sun and 1.10 M_sun. The corresponding allowable age range for these models, which depends sensitively on the parameters of the model, is relatively large and is ~2.5 - 5.5 Gyr. For each of the 75 models, a non-radial pulsation analysis is carried out, and the large and small frequency spacings are calculated. The results show that just measuring the large and small frequency spacings will greatly reduce the present uncertainties in the derived physical parameters and in the age of 51 Peg. Finally we discuss briefly refinements in the physics of the models and in the method of analysis which will have to be included in future models to make the best of the precise frequency determinations expected from space observations.

This paper describes a series of 3D simulations of shallow inefficient convection in the outer layers of the Sun. The computational domain is a closed box containing the convection-radiation transition layer, located at the top of the solar convection zone. The most salient features of the simulations are that: i)The position of the lower boundary can have a major effect on the characteristics of solar surface convection (thermal structure, kinetic energy and turbulent pressure). ii)The width of the box has only a minor effect on the thermal structure, but a more significant effect on the dynamics (rms velocities). iii)Between the surface and a depth of 1 Mm, even though the density and pressure increase by an order of magnitude, the vertical correlation length of vertical velocity is always close to 600 km. iv) In this region the vertical velocity cannot be scaled by the pressure or the density scale height. This casts doubt on the applicability of the mixing length theory, not only in the superadiabatic layer, but also in the adjacent underlying layers. v) The final statistically steady state is not strictly dependent on the initial atmospheric stratification.

We present a database of the latest stellar models of the $Y^2$ (Yonsei-Yale) collaboration. This database contains the stellar evolutionary tracks from the pre-main-sequence birthline to the helium core flash that were used to construct the $Y^2$ isochrones. We also provide a simple interpolation routine that generates stellar tracks for given sets of parameters (metallicity, mass, and $\alpha$-enhancement).

This is the first of a series of three papers devoted to the calibration of a few parameters of crucial importance in the modeling of the evolution of intermediate-mass stars, with special attention to the amount of convective core overshoot. To this end we acquired deep V and R photometry for three globular clusters of the Large Magellanic Cloud (LMC), namely NGC 2173, SL 556 and NGC 2155, in the age interval 1-3 Gyr. In this first paper, we describe the aim of the project, the VLT observations and data reduction, and we make preliminary comparisons of the color-magnitude diagrams with both Padova and Yonsei-Yale isochrones. Two following papers in this series present the results of a detailed analysis of these data, independently carried out by members of the Yale and Padova stellar evolution groups. This allows us to compare both sets of models and discuss their main differences, as well as the systematic effects that they would have to the determination of the ages and metallicities of intermediate-age single stellar populations.

We present a new set of isochrones in which the effect of the alpha-element enhancement is fully incorporated. These isochrones are an extension of the already published set of YY Isochrones (Yi et al. 2001: Paper 1), constructed for the scaled-solar mixture. As in Paper 1, helium diffusion and convective core overshoot have been taken into account.The range of chemical compositions covered is 0.00001 < Z < 0.08. The models were evolved from the pre-main-sequence stellar birthline to the onset of helium burning in the core. The age range of the full isochrone set is 0.1 -- 20 Gyr, while younger isochrones of age 1 -- 80 Myr are also presented up to the main-sequence turn-off. Combining this set with that of Paper 1 for scaled-solar mixture isochrones, we provide a consistent set of isochrones which can be used to investigate populations of any value of alpha-enhancement. We confirm the earlier results of Paper 1 that inclusion of alpha-enhancement effects further reduces the age estimates of globular clusters by approximately 8 percent if [alpha/Fe]=+0.3. It is important to note the metallicity dependence of the change in age estimates (larger age reductions in lower metallicities). This reduces the age gap between the oldest metal-rich and metal-poor Galactic stellar populations and between the halo and the disk populations. The isochrone tables, together with interpolation routines have been made available via internet; http://www.astro.yale.edu/demarque/yyiso.html http://www-astro.physics.ox.ac.uk/~yi/yyiso.html http://csaweb.yonsei.ac.kr/~kim/yyiso.html

We present an analysis of CMDs of three intermediate-age LMC clusters, namely NGC 2173, SL 556 and NGC 2155. The main goal of our project is to investigate the amount of convective core overshoot necessary to reproduce the CMDs of relatively metal-poor, intermediate age stellar populations. We conclude that a moderate amount of overshoot and some fraction of binary stars are essential for reproducing the observed shapes around the turnoff in the CMD's of all three clusters: unresolved binary stars fill in the expected core contraction gap, and make a unique sequence near the gap, which cannot be reproduced by single stars alone, even with a larger amount of overshoot. From our overall analysis such as, shape of isochrones, star counts, color distribution, and synthetic CMD comparisons, we conclude that overshoot ~ 20% of the local pressure scale height best reproduces the CMD properties of all three clusters. The best age estimates are 1.5, 2.1 and 2.9 Gyr for NGC 2173, SL 556 and NGC 2155, respectively.

The Sun is believed to have been the recipient of a substantial amount of metal-rich material over the course of its evolution, particularly in the early stages of the Solar System. With a long diffusion timescale, the majority of this accreted matter should still exist in the solar convection zone, enhancing its observed surface abundance, and implying a lower-abundance core. While helioseismology rules out solar models with near-zero metallicity cores, some solar models with enhanced metallicity in the convection zone might be viable, as small perturbations to the standard model. Because of the reduced interior opacity and core temperature, the neutrino flux predicted for such models is lower than that predicted by the standard solar model. This paper examines how compatible inhomogeneous solar models of this kind are with the observed low and intermediate degree p-mode oscillation data, and with the solar neutrino data from the SNO Collaboration. We set an upper limit on how much metal-rich accretion took place during the early evolution of the Sun at about 2 Earth masses of iron (or about 40 Earth masses of meteoric material).

The general consensus is that in order to reproduce the observed solar p-mode oscillation frequencies, turbulence should be included in solar models. However, until now there has not been any well-tested efficient method to incorporate turbulence into solar modeling. We present here two methods to include turbulence in solar modeling within the framework of the mixing length theory, using the turbulent velocity obtained from numerical simulations of the highly superadiabatic layer of the sun at three stages of its evolution. The first approach is to include the turbulent pressure alone, and the second is to include both the turbulent pressure and the turbulent kinetic energy. The latter is achieved by introducing two variables: the turbulent kinetic energy per unit mass, and the effective ratio of specific heats due to the turbulent perturbation. These are treated as additions to the standard thermodynamic coordinates (e.g. pressure and temperature). We investigate the effects of both treatments of turbulence on the structure variables, the adiabatic sound speed, the structure of the highly superadiabatic layer, and the p-mode frequencies. We find that the second method reproduces the SAL structure obtained in 3D simulations, and produces a p-mode frequency correction an order of magnitude better than the first method.

We assess the potential of asteroseismology for determining the fundamental properties of individual $\delta$ Scuti stars. We computed a grid of evolution and adiabatic pulsation models using the Yale Rotating Evolution Code to study the systematic changes in low-order ($\ell = 0, 1, 2,$ and 3) modes as functions of fundamental stellar properties. Changes to the stellar mass, chemical composition, and convective core overshooting length change the observable pulsation spectrum significantly. In general, mass has the strongest effect on evolution and on pulsation, followed by the metal abundance. Changes to the helium content have very little effect on the frequencies until near the end of the main sequence. Changes to each of the four parameters change the $p$-mode frequencies more, both in absolute and relative terms, than they do the $g$- and mixed-mode frequencies, suggesting that these parameters have a greater effect on the outer layers of the star. We also present evolution and pulsation models of the well-studied star FG Virginis, outlining a possible method of locating favorable models in the stellar parameter space based upon a definitive identification of only two modes. Specifically, we plot evolution models on the (period-period ratio) and (temperature-period ratio) planes to select candidate models, and modify the core overshooting parameter to fit the observed star. For these tests, we adjusted only the mass, helium and metal abundances, and core overshooting parameter, but this method can be extended to include the effects of first-order rotational splitting and second-order rotational distortion of pulsation spectra.

In stellar evolution calculations, the local pressure scale height is often used to empirically constrain the amount of convective core overshoot. However, this method brings unsatisfactory results for low-mass stars (< 1.1 -1.2 Mo for Z= \Zs) which have very small cores or no convective core at all. Following Roxburgh's integral constraint,we implemented an upper limit of overshoot within the conventional method of alpha parameterization in order to remove an overly large overshoot effect on low-mass stars. The erroneously large effect of core overshoot due to the failure of alpha parameterization can be effectively corrected by limiting the amount of overshoot to < 15 % of the core radius. 15 % of the core radius would be a proper limit of overshoot which can be implemented in a stellar evolution code for intermediate to low mass stars. The temperature structure of the overshoot region does not play a crucial role in stellar evolution since this transition region is very thin.

We have constructed a new set of isochrones, called the Y^2 Isochrones, that represent an update of the Revised Yale Isochrones (RYI), using improved opacities and equations of state. Helium diffusion and convective core overshoot have also been taken into consideration. This first set of isochrones is for the scaled solar mixture. A subsequent paper will consider the effects of alpha-element enhancement, believed to be relevant in many stellar systems. Two additionally significant features of these isochrones are that (1) the stellar models start their evolution from the pre-main sequence birthline instead of from the zero-age main sequence, and (2) the color transformation has been performed using both the latest table of Lejeune et al., and the older, but now modified, Green et al. table. The isochrones have performed well under the tests conducted thus far. The reduction in the age of the Galactic globular clusters caused by this update in stellar models alone is approximately 15% relative to RYI-based studies. When the suggested modification for the alpha-element enhancement is made as well, the total age reduction becomes approximately 20%. When post-RGB evolutionary stages are included, we find that the ages of globular clusters derived from integrated colors are consistent with the isochrone fitting ages. Isochrones are available from http://achee.srl.caltech.edu.

We present results of three dimensional simulations of the uppermost part of the sun, at 3 stages of its evolution. Each model includes physically realistic radiative-hydrodynamics (the Eddington approximation is used in the optically thin region), varying opacities and a realistic equation of state (full treatment of the ionization of H and He). In each evolution model, we investigate a domain, which starts at the top of the photosphere and ends just inside the convection zone (about 2400 km in the sun model). This includes all of the super-adiabatic layer (SAL). Due to the different positions of the three models in the $log (g) $ vs $log T_{eff}$ plane, the more evolved models have lower density atmospheres. The reduction in density causes the amount of overshoot into the radiation layer, to be greater in the more evolved models.

We have constructed models for the sun at three stages of its evolution: a zero-age main sequence model, the present sun, and a subgiant model. For each model, the turbulent pressure and turbulent kinetic energy were calculated from 3-d radiative hydrodynamical simulations (described in the poster by Robinson et al.), and inserted into the 1-d stellar models. We note that in these simulations, the turbulent pressure is not a free parameter, but can be computed from the resulting velocity field. We show the calculated p-mode frequencies for the model of the present sun, with and without turbulent pressure, and compare them to the observed solar frequencies. When the turbulent pressure is included in the models, the calculated frequencies are brought closer to the observed frequencies in the sun by up to two $\mu Hz$, strictly from structural effects. The effect of including turbulent pressure on p-mode frequencies is also shown for the zero-age main sequence model. Our models also suggest that the importance of turbulent pressure increases as the star evolves into the subgiant region. We discuss the importance of also including realistic turbulence as well as radiation in the non-adiabatic calculation of oscillation frequencies.

The presence of a substantial number of hot stars in the extremely metal-rich} open cluster NGC6791 has been a mystery. If these hot stars are in their core helium burning phase, they are significantly bluer (hotter) than predicted by canonical stellar evolution theory. No obvious explanation is available yet. We consider the effects of mass loss during the evolution of horizontal branch (HB) stars as their possible origin. We find that the addition of mass loss causes HB stars to evolve to be hotter and fainter. Mass loss has a more pronounced effect for less massive stars and thus naturally widens the temperature (and color) distribution on the HB. If mass loss rates are higher for more metal-rich stars, this phenomenon would be even more pronounced in the metal-rich populations, such as NGC6791. We find that mass loss on the HB may be a viable method of forming subdwarf B (sdB) stars both in the field and in clusters, especially when the metallicity is high.

A grid of synthetic horizontal-branch (SHB) models based on HB evolutionary tracks with improved physics has been constructed to reconsider the theoretical calibration of the dependence of M_v(RR) on metallicity in globular clusters, and the slope of the mean <M_v(RR)>-[Fe/H] relation. The SHB models confirm Lee's earlier finding (Lee 1991) that the slope of the <M_v(RR)>-[Fe/H] relation is itself a function of the metallicity range considered, and that in addition, for a given [Fe/H], RR Lyrae luminosities depend on HB morphology. This is due to the fact that HB stars pass through the RR Lyrae instability strip at different evolutionary stages, depending on their original position on the HB. At [Fe/H]=-1.9, and for HB type 0, the models yield M_v(RR)=0.47 \pm 0.10. The mean slope for the zero-age HB models is 0.204. Since there is no simple universal relation between M_v(RR) and metallicity that is applicable to all globular clusters, the HB morphology of each individual cluster must be taken into account, in addition to [Fe/H], in deriving the appropriate M_v(RR). Taking HB morphology into account, we find that the slope of the mean <M_v(RR)>-[Fe/H] relation varies between 0.36 for the clusters with galactocentric distances R_gc less than 6 kpc and 0.22 for clusters with 6<R_gc<20 kpc. Implications for interpreting observations of field RR Lyrae variables and for absolute globular cluster ages and galactic chronology are briefly discussed.

A grid of stellar evolution models for Procyon A has been calculated. These models include the best physics available to us (including the latest opacities and equation of state) and are based on the revised astrometric mass of Girard et al (1996). Models were calculated with helium diffusion and with the combined effects of helium and heavy element diffusion. Oscillation frequencies for l=0,1,2 and 3 p-modes and the characteristic period spacing for the g-modes were calculated for these models. We find that g-modes are sensitive to model parameters which effect the structure of the core, such as convective core overshoot, the heavy element abundance and the evolutionary state (main sequence or shell hydrogen burning) of Procyon A. The p-modes are relatively insensitive to the details of the physics used to model Procyon A, and only depend on the evolutionary state of Procyon A. Hence, observations of p-mode frequencies on Procyon A will serve as a robust test of stellar evolution models.

We present a review of the present state of knowledge regarding the relative ages of Galactic globular clusters. First, we discuss the relevant galaxy formation models and describe the detailed predictions they make with respect to the formation timescale and chemical evolution of the globular clusters. Next, the techniques used to estimate globular cluster ages are described and evaluated with particular emphasis on the advantages and disadvantages of each method. With these techniques as a foundation, we present arguments in favor of the following assertions: 1) The age of a globular cluster is the likeliest candidate to be the global second parameter, which along with metal abundance, controls the morphology of the horizonal branch. 2) A total age range of as much as $\sim$5 Gyr exists among the bulk of the Galactic globulars. 3) There is a significant relation between age and metallicity among the Galactic globular clusters if the slope of the \mvrr-\feh relation is less than $\sim$0.23. These conclusions along with other supporting evidence favor a formation scenario in which the inner regions of the Galactic halo collapsed in a monotonic fashion over a short time period much less than 1 Gyr. In contrast, the outer regions of the halo fragmented and collapsed in a chaotic manner over several Gyrs.

A grid of stellar evolution models for Procyon A has been calculated. These models include the best physics available to us (including the latest opacities and equation of state) and are based on the revised astrometric mass of Girard et al. (1996). The long standing discrepancy between the evolutionary mass and the astrometric mass is now resolved, a result of the newly determined astrometric mass. Models were calculated with helium diffusion and with the combined effects of helium and heavy element diffusion. Oscillation frequencies for $\ell=0,1,2$ and 3 $p$-modes (and $g$-modes) were calculated for these models. The predicted $p$-mode frequencies are relatively unaffected by heavy element diffusion and convective core overshoot. The inclusion of a modest stellar wind which effectively suppresses the helium diffusion in the surface layers has a modest effect on the $p$-mode frequencies. The evolutionary state (main sequence or shell hydrogen burning) of Procyon A has the largest effect on the predicted $p$-mode frequencies. The $g$-modes show a greater sensitivity to the various model parameters.

Whether or not metal-rich HB stars are the dominant UV source in giant elliptical galaxies (gEs) is an important question in current astronomical research. We follow up our previous evolutionary population synthesis study with quantitative tests to answer this question affirmatively under the following three conditions: (1) Reimers' empirical mass loss formula is proper, (2) the mass loss efficiency parameter in metal-rich stars is somewhat larger than the value estimated from the metal-poor star studies, and (3) the true value of the helium enrichment parameter is positive. All three important empirical characteristics of the UV upturn (i.e., the fact that strong UV upturns are restricted to gEs, the positive UV upturn-metallicity correlation, and the narrow range of the temperature of the UV sources) are closely reproduced for reasonable ranges of input parameters. We discuss the major sources of uncertainties in the models, such as the production and role of hot horizontal-branch stars in gEs, and the importance of galactic nucleosynthesis.

We review five independent techniques which are used to set the distance scale to globular clusters, including subdwarf main sequence fitting utilizing the recent Hipparcos parallax catalogue. These data together all indicate that globular clusters are farther away than previously believed, implying a reduction in age estimates. This new distance scale estimate is combined with a detailed numerical Monte Carlo study designed to assess the uncertainty associated with the theoretical age-turnoff luminosity relationship in order to estimate both the absolute age and uncertainty in age of the oldest globular clusters. Our best estimate for the mean age of the oldest globular clusters is now $11.5\pm 1.3 $Gyr, with a one-sided, 95% confidence level lower limit of 9.5 Gyr. This represents a systematic shift of over 2 $\sigma$ compared to our earlier estimate, due completely to the new distance scale---which we emphasize is not just due to the Hipparcos data. This now provides a lower limit on the age of the universe which is consistent with either an open universe, or a flat, matter dominated universe (the latter requiring $H_0 \le 67 \kmsmpc$). Our new study also explicitly quantifies how remaining uncertainties in the distance scale and stellar evolution models translate into uncertainties in the derived globular cluster ages. Simple formulae are provided which can be used to update our age estimate as improved determinations for various quantities become available.

We present models of the late stages of stellar evolution intended to explain the UV upturn phenomenon in elliptical galaxies. Such models are sensitive to values of a number of poorly-constrained physical parameters, including metallicity, age, stellar mass loss, helium enrichment, and the distribution of stars on the zero age horizontal branch (HB). We explore the sensitivity of the results to values of these parameters, and reach the following conclusions. Old, metal rich galaxies, such as giant ellipticals, naturally develop a UV upturn within a reasonable time scale - less than a Hubble time - without the presence of young stars. The most likely stars to dominate the UV flux of such populations are low mass, core helium burning (HB and evolved HB) stars. Metal-poor populations produce a higher ratio of UV-to-V flux, due to opacity effects, but only metal-rich stars develop a UV upturn, in which the flux increases towards shorter UV wavelengths. Model color-magnitude diagrams and corresponding integrated spectra (for various values of age, metallicity, helium enrichment, mass loss efficiency, initial mass function, and the HB mass dispersion factor) are available on S.Y.'s world wide web site http://shemesh.gsfc.nasa.gov/model.html.

We consider the origin of the UV bright phase of metal-rich helium-burning stars, the slow blue phase (SBP), that was predicted by various earlier works. Based on improved physics including OPAL opacities, we confirm the existence of the SBP. In addition to our grid of evolutionary tracks, we provide an analytical understanding of the main characteristics of the SBP phenomenon. The SBP is slow because it is a slow evolving helium-shell-burning phase which is analogous to the early AGB phase. The SBP of a more metal-rich star is slower than a metal-poor counterpart if their Teff's are the same because a more metal-rich helium-burning star has a smaller mass than a metal-poor one and because lifetime increases as mass decreases. Metal-rich helium-burning stars easily become hot because the luminosity from the hydrogen-burning shell is extremely sensitive to the mean molecular weight whereas the luminosity from the helium-burning core is not. We suggest that the SBP phenomenon is a major cause of the UV upturn phenomenon in giant elliptical galaxies as will be shown in subsequent papers. The new HB tracks can be retrieved from S.Y.'s web site http://shemesh.gsfc.nasa.gov/astronomy.html.

We explore the evolution of collisionally merged stars in the blue straggler region of the HR diagram. The starting models for our stellar evolution calculations are the results of the smoothed particle hydrodynamics (SPH) simulations of parabolic collisions between main sequence stars performed by Lombardi, Rasio and Shapiro (1996). Since SPH and stellar evolution codes employ different and often contradictory approximations, it is necessary to treat the evolution of these products carefully. The mixture and disparity of the relevant timescales (hydrodynamic, thermal relaxation and nuclear burning) and of the important physical assumptions between the codes makes the combined analysis of the problem challenging, especially during the initial thermal relaxation of the star. In particular, the treatment of convection is important, and semiconvection must be modeled in some detail. The products of seven head-on collisions are evolved through their initial thermal relaxation, and then through the main sequence phase to the base of the giant branch. Their evolutionary tracks are presented. In contrast to the assumptions in previous work, these collision products do not develop substantial convective regions during their thermal relaxation, and therefore are not mixed significantly after the collision.

Globular cluster age estimates based on the absolute magnitude of the main sequence turn-off (\mvto) are generally considered to be the most reliable from a theoretical viewpoint. However, the difficulty in determining \mvto\ in observed colour-magnitude diagrams leads to a large error in the derived age. In this paper, we advocate the use of the absolute magnitude of the point which is brighter than the turn-off and 0.05 mag redder (\mvb) as a precision age indicator. It is easy to measure this point on observed colour-magnitude diagrams, leading to small observational error bars. Furthermore, an extensive Monte Carlo calculation indicates that the theoretical uncertainty in \mvb\ is similar to \mvto. As a result, ages derived using \mvb\ are at least a factor of 2 more precise than those derived using \mvto. This technique is applied to the globular cluster M68 and an age of $12.8\pm 0.3\,$Gyr is derived (assuming $\mvrr = 0.20\,\feh + 0.98$), indicating that M68 is a `young' globular cluster. A homogeneous set of globular cluster age estimates with this precision would provide unprecedented insight into the formation of the Galactic halo.

We report the results of a detailed numerical study designed to estimate both the absolute age and the uncertainty in age (with confidence limits) of the oldest globular clusters. Such an estimate is essential if a comparison with the Hubble age of the universe is to be made to determine the consistency, or lack thereof, of various cosmological models. Utilizing estimates of the uncertainty range (and distribution) in the input parameters of stellar evolution codes we produced 1000 Monte Carlo realizations of stellar isochrones, with which we could fit the ages of the 18 oldest globular clusters. Incorporating the observational uncertainties in the measured color-magnitude diagrams for these systems and the predicted isochrones, we derived a probability distribution for the mean age of these systems. The one-sided 95$\%$ C.L.~ lower bound for this distribution occurs at an age of 12.07 Gyr. This puts interesting constraints on cosmology which we discuss. Further details, including a description of the distributions, covariance matrices, dependence upon individual input parameters, etc.\ will appear in a future article.

Main sequence turnoff magnitudes from the recent set of Yale isochrones (Chaboyer \ea 1995) have been combined with a variety of relations for the absolute magnitude of RR Lyr stars (\mvrr) to calibrate age as a function of the difference in magnitude between the main sequence turn-off and the horizontal branch (\dv). A best estimate for the calibration of \mvrr is derived from a survey of the current literature: $\mvrr = 0.20 \feh + 0.98$. This estimate, together with other calibrations (with slopes ranging from 0.15 to 0.30) has been used to derive \dv ages for 43 Galactic globular clusters. Independent of the choice of \mvrr, there is no strong evidence for an age-Galacto\-centric distance relationship among the 43 globular clusters. However, an age-metallicity relation exists, with the metal-poor clusters being the oldest. A study of the age distribution reveals that an age range of 5 Gyr exists among the bulk of the globular clusters. In addition, about 10\% of the sample are substantially younger, and including them in the analysis increases the age range to 9 Gyr. Once again, these statements are independent of the \mvrr ~relation. Evidence for age being the second parameter governing horizontal branch morphology is found by comparing the average \dv ~age of the second parameter clusters to the normal clusters. The second parameter clusters are found to be on average 2 - 3 Gyr younger than the other clusters, which is consistent with age being the second parameter. These results suggest that globular clusters were formed over an extended period of time, with progressively more metal-rich globular clusters ($\feh \ga -1.7$) being formed at later times.

This paper provides estimates of convective turnover time scales for Sun-like stars in the pre-main sequence and early post-main sequence phases of evolution, based on up-to-date physical input for the stellar models. In this first study, all models have solar abundances, which is typical of the stars in the Galactic disk where most of the available data have been collected. A new feature of these models is the inclusion of rotation in the evolutionary sequences, thus making it possible to derive theoretically the Rossby number for each star along its evolutionary track, based on its calculated rotation rate and its local convective turnover time near the base of the convection zone. Global turnover times are also calculated for the complete convection zone. This information should make possible a new class of observational tests of stellar theory which were previously impossible with semi-empirical models, particularly in the study of stellar activity and in research related to angular momentum transfer in stellar interiors during the course of stellar evolution.

We have studied the importance of the combined effects of rotation, diffusion, and convective overshoot on the p-mode oscillation spectrum and the neutrino flux of the standard solar model. To isolate the various physical affects included in the new rotation plus diffusion models we also constructed solar models to test the significance of diffusion and of overshoot by themselves. In previous studies, models that include helium diffusion during solar evolution were found to improve the predicted p-mode frequencies for some modes and worsen the agreement for others (Guenther \ea 1993). Here we verify this result for both the Bahcall and Loeb (1990) formulation of diffusion and the Proffitt and Michaud (1991) formulation of diffusion. We find that the effects of rotation on the Sun's structure in the outer layers perturbs the $p$-mode frequencies only slightly when compared to the more substantial effects due to diffusion. In the thin overshoot layer (taken here to be $0.1\, H_p$), we have compared the effect of overmixing in a radiative layer versus convective (adiabatic) penetration. Neither radiative overmixing nor adiabatic penetration has any significant effect on the $p$-modes, probably in part because the overshoot layer is constrained to be thin. The predicted neutrino flux in our diffusion plus rotation model is 7.12 SNU for Cl detectors, 127 SNU for Ga detectors and $5.00\times 10^6\,{\rm erg \, cm^{-2} }$ for the $^8$B neutrinos; this is approximately half-way between the standard solar model without diffusion, and the standard solar model with diffusion alone.

A new grid of theoretical isochrones based on the Yale stellar evolution code using the OPAL and Kurucz opacities has been constructed. The grid of isochrones spans a wide range of metallicities, helium abundances and masses. The construction of the isochrones is described and the isochrones are compared to galactic globular cluster observations. A solar calibrated mixing length ($\alpha = 1.7$) yields a good fit to globular cluster colour-magnitude diagrams. Ages for 40 globular clusters are determined using the \dv method and the formation of the halo is discussed.

The Yale stellar evolution code has been modified to include the combined effects of diffusion and rotational mixing on $^1$H, $^4$He and the trace elements $^3$He, $^6$Li, $^7$Li, and $^9$Be. The interaction between rotational mixing and diffusion is studied by calculating a number of calibrated solar models. The rotational mixing inhibits the diffusion in the outer parts of the models, leading to a decrease in the envelope diffusion by 25 -- 50\%. Conversely, diffusion leads to gradients in mean molecular weight which can inhibit the rotational mixing. The degree to which gradients in mean molecular weight inhibits the rotational mixing is somewhat uncertain. A comparison to the observed solar oblateness suggests that gradients in the mean molecular weight play a smaller role in inhibiting the rotational mixing previously believed. This is reinforced by the fact that the model with the standard value for the inhibiting effect of mean molecular weight on the rotational mixing depletes no Li on the main sequence. This is in clear contrast to the observations. A reduction in the inhibiting effect of mean molecular weight gradients by a factor of ten leads to noticeable main sequence Li depletion.

Stellar models with masses ranging from 0.5 to $1.3~M_\odot$ were constructed in order to compare to young cluster observations of Li and of rotation velocities. The amount of Li depletion in cool stars is sensitive to the amount of overshoot at the base of the surface convection zone, and the exact metallicity of the models. Even when this uncertainty is taken into account, the Li observations are a severe constraint for the models and rule out standard models, and pure diffusion models. Stellar models which include diffusion and rotational mixing in the radiative regions of stars are able to simultaneously match the Li abundances observed in the Pleiades, UMaG, Hyades, Praesepe, NGC 752 and M67. They also match the observed rotation periods in the Hyades. However, these models are unable to simultaneously explain the presence of the rapidly rotating late G and K stars in the Pleiades and the absence of rapidly rotating late F and early G stars.

In an attempt to understand the properties of convective energy transport in the solar convection zone, a numerical model has been constructed for turbulent flows in a compressible, radiation-coupled, non-magnetic, gravitationally stratified medium using a realistic equation of state and realistic opacities. The time-dependent, three-dimensional hydrodynamic equations are solved with minimal simplifications. The statistical information obtained from the present simulation provides an improved understanding of solar photospheric convection. The characteristics of solar convection in shallow regions is parameterized and compared with the results of Chan and Sofia's simulations of deep and efficient convection (Chan and Sofia 1989). We assess the importance of the zones of partial ionization in the simulation, and confirm that the radiative energy transfer is negligible throughout the region except in the uppermost scale heights of the convection zone, a region of very high super-adiabaticity. When the effects of partial ionization are included, the dynamics of flows are altered significantly. However, we confirm the Chan and Sofia result that kinetic energy flux is non-negligible and can have a negative value in the convection zone.

A large number of stellar evolution models with [Fe/H] = -2.3 and -3.3 have been calculated in order to determine the primordial $^7$Li abundance and to test current stellar evolution models by a comparison to the extensive database of Li abundances in extremely metal poor halo stars observed by Thorburn (1994). Standard models do a good job of fitting the observed Li abundances in stars hotter than 5600 K. They predict a primordial $^7$Li abundance of Log N(Li) = 2.24\pm 0.03$. Models which include microscopic diffusion predict a downward curvature in the $^7$Li destruction isochrones at hot temperatures which is not present in the observations. Thus, the observations clearly rule out models which include uninhibited microscopic diffusion of $^7$Li from the surface of the star. The [Fe/H] = -2.3 stellar models which include both diffusion and rotational mixing provide an excellent match to the observations. Both the plateau stars and the heavily depleted cool stars are well fit by these models. The rotational mixing leads to considerable $^7$Li depletion in these models and the primordial $^7$Li abundance inferred from these models is Log N(Li) = $3.08\pm 0.1$.