Showing 1–50 of 74 results for author: Demarque, P

Evidence for strong progenitor age dependence of type Ia supernova luminosity standardization process

Authors: Young-Wook Lee, Chul Chung, Pierre Demarque, Seunghyun Park, Junhyuk Son, Yijung Kang

Abstract: Supernova (SN) cosmology is based on the assumption that the width-luminosity relation (WLR) and the color-luminosity relation (CLR) in the type Ia SN luminosity standardization would not show {absolute magnitude differences} with progenitor age. Unlike this expectation, recent age datings of stellar populations in host galaxies have shown significant correlations between progenitor age and Hubble… ▽ More Supernova (SN) cosmology is based on the assumption that the width-luminosity relation (WLR) and the color-luminosity relation (CLR) in the type Ia SN luminosity standardization would not show {absolute magnitude differences} with progenitor age. Unlike this expectation, recent age datings of stellar populations in host galaxies have shown significant correlations between progenitor age and Hubble residual (HR). Here we show that this correlation originates from a strong progenitor age dependence of the zero-points of the WLR and the CLR, in the sense that SNe from younger progenitors are fainter each at given light-curve parameters $x_1$ and $c$. This $4.6σ$ result is reminiscent of Baade's discovery of the zero-point variation of the Cepheid period-luminosity relation with age, and, as such, causes a serious systematic bias with redshift in SN cosmology. Other host properties show substantially smaller and insignificant offsets in the WLR and CLR for the same dataset. We illustrate that the differences between the high-$z$ and low-$z$ SNe in the WLR and CLR, and in HR after the standardization, are fully comparable to those between the correspondingly young and old SNe at intermediate redshift, indicating that the observed dimming of SNe with redshift may well be an artifact of over-correction in the luminosity standardization. When this systematic bias with redshift is properly taken into account, there is little evidence left for an accelerating universe, in discordance with other probes, urging the follow-up investigations with larger samples at different redshift bins. △ Less

Submitted 29 September, 2022; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: 13 pages, 11 figures, 1 table, Accepted for publication in MNRAS

Stellar evolution models with entropy-calibrated mixing-length parameter: application to red giants

Authors: Federico Spada, Pierre Demarque, Friedrich Kupka

Abstract: We present evolutionary models for solar-like stars with an improved treatment of convection that results in a more accurate estimate of the radius and effective temperature. This is achieved by improving the calibration of the mixing-length parameter, which sets the length scale in the 1D convection model implemented in the stellar evolution code. Our calibration relies on the results of 2D and 3… ▽ More We present evolutionary models for solar-like stars with an improved treatment of convection that results in a more accurate estimate of the radius and effective temperature. This is achieved by improving the calibration of the mixing-length parameter, which sets the length scale in the 1D convection model implemented in the stellar evolution code. Our calibration relies on the results of 2D and 3D radiation hydrodynamics simulations of convection to specify the value of the adiabatic specific entropy at the bottom of the convective envelope in stars as a function of their effective temperature, surface gravity and metallicity. For the first time, this calibration is fully integrated within the flow of a stellar evolution code, with the mixing-length parameter being continuously updated at run-time. This approach replaces the more common, but questionable, procedure of calibrating the length scale parameter on the Sun, and then applying the solar-calibrated value in modeling other stars, regardless of their mass, composition and evolutionary status. The internal consistency of our current implementation makes it suitable for application to evolved stars, in particular to red giants. We show that the entropy calibrated models yield a revised position of the red giant branch that is in better agreement with observational constraints than that of standard models. △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: 12 pages, 10 figures; accepted for publication in MNRAS

Testing the entropy calibration of the radii of cool stars: models of alpha Centauri A and B

Authors: F. Spada, P. Demarque

Abstract: 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 stel… ▽ More 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. △ Less

Submitted 2 September, 2019; originally announced September 2019.

Comments: 9 pages, 7 figures; accepted for publication in Monthly Notices of the Royal Astronomical Society

Improved calibration of the radii of cool stars based on 3D simulations of convection: implications for the solar model

Authors: F. Spada, P. Demarque, S. Basu, J. D. Tanner

Abstract: 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 T… ▽ More 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. △ Less

Submitted 5 November, 2018; originally announced November 2018.

Comments: 16 pages, 14 figures, accepted for publication in the Astrophysical Journal

The Sirius System and its Astrophysical Puzzles: Hubble Space Telescope and Ground-Based Astrometry

Authors: Howard E. Bond, Gail H. Schaefer, Ronald L. Gilliland, Jay B. Holberg, Brian D. Mason, Irving W. Lindenblad, Miranda Seitz-McLeese, W. David Arnett, Pierre Demarque, Federico Spada, Patrick A. Young, Martin A. Barstow, Matthew R. Burleigh, Donald Gudehus

Abstract: 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 historica… ▽ More 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. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Comments: Accepted by Astrophysical Journal

The Yale-Potsdam Stellar Isochrones (YaPSI)

Authors: F. Spada, P. Demarque, Y. -C. Kim, T. S. Boyajian, J. M. Brewer

Abstract: 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-hos… ▽ More 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 △ Less

Submitted 11 March, 2017; originally announced March 2017.

Comments: 17 pages, 14 figures; accepted for publication in the Astrophysical Journal

Entropy in adiabatic regions of convection simulations

Authors: Joel D. Tanner, Sarbani Basu, Pierre Demarque

Abstract: 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… ▽ More 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. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: Accepted for publication in ApJ Letters

Hubble Space Telescope Astrometry of the Procyon System

Authors: Howard E. Bond, Ronald L. Gilliland, Gail H. Schaefer, Pierre Demarque, Terrence M. Girard, Jay B. Holberg, Donald Gudehus, Brian D. Mason, Vera Kozhurina-Platais, Matthew R. Burleigh, Martin A. Barstow, Edmund P. Nelan

Abstract: 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 compone… ▽ More 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. △ Less

Submitted 2 October, 2015; originally announced October 2015.

Comments: Accepted by Astrophysical Journal

Observations of Binary Stars with the Differential Speckle Survey Instrument. V. Toward an Empirical Metal-Poor Mass-Luminosity Relation

Authors: Elliott P. Horch, William F. van Altena, Pierre Demarque, Steve B. Howell, Mark E. Everett, David R. Ciardi, Johanna K. Teske, Todd J. Henry, Jennifer G. Winters

Abstract: 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 t… ▽ More 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. △ Less

Submitted 1 March, 2015; originally announced March 2015.

Stellar Diameters and Temperatures VI. High angular resolution measurements of the transiting exoplanet host stars HD 189733 and HD 209458 and implications for models of cool dwarfs

Authors: Tabetha Boyajian, Kaspar von Braun, Gregory A. Feiden, Daniel Huber, Sarbani Basu, Pierre Demarque, Debra A. Fischer, Gail Schaefer, Andrew W. Mann, Timothy R. White, Vicente Maestro, John Brewer, C. Brooke Lamell, Federico Spada, Mercedes López-Morales, Michael Ireland, Chris Farrington, Gerard T. van Belle, Stephen R. Kane, Jeremy Jones, Theo A. ten Brummelaar, David R. Ciardi, Harold A. McAlister, Stephen Ridgway, P. J. Goldfinger , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 20 November, 2014; originally announced November 2014.

Modeling Convective Core Overshoot and Diffusion in Procyon Constrained by Asteroseismic Data

Authors: D. B. Guenther, P. Demarque, M Gruberbauer

Abstract: 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… ▽ More 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. △ Less

Submitted 24 April, 2014; originally announced April 2014.

The Effect of Metallicity-Dependent T-tau Relations on Calibrated Stellar Models

Authors: Joel D. Tanner, Sarbani Basu, Pierre Demarque

Abstract: 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 r… ▽ More 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. △ Less

Submitted 31 March, 2014; v1 submitted 5 February, 2014; originally announced February 2014.

Journal ref: 2014, ApJ, 785, L13

Helium-Abundance and Other Composition Effects on the Properties of Stellar Surface Convection in Solar-like Main-sequence Stars

Authors: Joel D. Tanner, Sarbani Basu, Pierre Demarque

Abstract: We investigate the effect of helium abundance and $α$-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 gri… ▽ More We investigate the effect of helium abundance and $α$-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 $α$-element enhancement and find that it has a considerably smaller effect on the convective dynamics in the SAL compared to that of helium abundance. △ Less

Submitted 1 October, 2013; originally announced October 2013.

Comments: Accepted for publication in ApJ

The radius discrepancy in low mass stars: single vs. binaries

Authors: F. Spada, P. Demarque, Y. -C. Kim, A. Sills

Abstract: 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 c… ▽ More 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 α. 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. △ Less

Submitted 26 August, 2013; originally announced August 2013.

Comments: Accepted for publication in the Astrophysical Journal

Variation of Stellar Envelope Convection and Overshoot with Metallicity

Authors: Joel D. Tanner, Sarbani Basu, Pierre Demarque

Abstract: 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… ▽ More 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. △ Less

Submitted 22 February, 2013; originally announced February 2013.

Comments: Accepted for publication in ApJ

Comparing the Effect of Radiative Transfer Schemes on Convection Simulations

Authors: Joel D. Tanner, Sarbani Basu, Pierre Demarque

Abstract: 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 moment… ▽ More 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. △ Less

Submitted 17 September, 2012; originally announced September 2012.

Modelling fully convective stars in eclipsing binaries: KOI-126 and CM Draconis

Authors: Federico Spada, Pierre Demarque

Abstract: 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 parame… ▽ More 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 α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 α. As previously noted in the literature, a reduced α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 α(for α< 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 αalone. We therefore suggest that the combined effects of high metallicity and α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 σwith α= 1.0. △ Less

Submitted 20 February, 2012; originally announced February 2012.

Comments: 7 pages, 4 figures, accepted for publication in MNRAS

Modeling Kepler Observations of Solar-like Oscillations in the Red-giant Star HD 186355

Authors: Chen Jiang, Biwei Jiang, Jorgen Christensen-Dalsgaard, Timothy R. Bedding, Dennis Stello, Daniel Huber, Søren Fandsen, Hans Kjeldsen, Christoffer Karoff, Benoit Mosser, Pierre Demarque, Michael N. Fanelli, Karen Kinemuchi, Fergal Mullally

Abstract: 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, ν_max, and the me… ▽ More 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, ν_max, and the mean large frequency separation, Δν, are around 106 and 9.4 μ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 Δν and ν_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. △ Less

Submitted 28 November, 2011; v1 submitted 5 September, 2011; originally announced September 2011.

Comments: 26 pages, 5 figures and 2 tables

Journal ref: 2011 ApJ 742 120

Solar-like oscillations in cluster stars

Authors: D. Stello, S. Basu, T. R. Bedding, K. Brogaard, H. Bruntt, W. J. Chaplin, J. Christensen-Dalsgaard, P. Demarque, Y. P. Elsworth, R. A. García, R. L. Gilliland, S. Hekker, D. Huber, C. Karoff, H. Kjeldsen, Y. Lebreton, S. Mathur, S. Meibom, J. Molenda-Żakowicz, A. Noels, I. W. Roxburgh, V. S. Aguirre, C. Sterken, R. Szabó

Abstract: 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. 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. △ Less

Submitted 30 June, 2010; originally announced July 2010.

Comments: 4 pages, 7 figures, accepted by Astronomische Nachrichten

Detection of solar-like oscillations from Kepler photometry of the open cluster NGC 6819

Authors: Dennis Stello, Sarbani Basu, Hans Bruntt, Benoit Mosser, Ian R. Stevens, Timothy M. Brown, Jorgen Christensen-Dalsgaard, Ronald L. Gilliland, Hans Kjeldsen, Torben Arentoft, Jerome Ballot, Caroline Barban, Timothy R. Bedding, William J. Chaplin, Yvonne P. Elsworth, Rafael A. Garcia, Marie-Jo Goupil, Saskia Hekker, Daniel Huber, Savita Mathur, Soren Meibom, Reza Samadi, Vinothini Sangaralingam, Charles S. Baldner, Kevin Belkacem , et al. (27 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 30 December, 2009; originally announced January 2010.

Comments: 5 pages, 4 figures, accepted by ApJ (Lett.)

2D solar modeling

Authors: P. Ventura, V. Penza, L. Li, S. Sofia, S. Basu, P. Demarque

Abstract: 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 sola… ▽ More 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. △ Less

Submitted 26 December, 2009; originally announced December 2009.

Comments: Accepted for publication on Astrophysics and Space Science

Journal ref: Astrophys.Space Sci.328:295-300,2010

Asteroseismology and interferometry of the red giant star epsilon Oph

Authors: A. Mazumdar, A. Merand, P. Demarque, P. Kervella, C. Barban, F. Baudin, V. Coude du Foresto, C. Farrington, P. J. Goldfinger, M. -J. Goupil, E. Josselin, R. Kuschnig, H. A. McAlister, J. Matthews, S. T. Ridgway, J. Sturmann, L. Sturmann, T. A. ten Brummelaar, N. Turner

Abstract: 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 measurement… ▽ More 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. △ Less

Submitted 18 June, 2009; originally announced June 2009.

Comments: 11 pages, accepted for publication in Astronomy & Astrophysics

Dissipation Efficiency in Turbulent Convective Zones in Low Mass Stars

Authors: Kaloyan Penev, Dimitar Sasselov, Frank Robinson, Pierre Demarque

Abstract: 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… ▽ More 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. △ Less

Submitted 25 September, 2009; v1 submitted 28 October, 2008; originally announced October 2008.

Comments: 13 pages, 3 figures Effective viscosity scaling changed by a factor of ~100. More details provided for the numerical models

Journal ref: Astrophys.J.704:930-936,2009

Calculation of Solar P-mode Oscillation Frequency Splittings Based on a Two-dimensional Solar Model

Authors: Linghuai Li, Sarbani Basu, Sabatino Sofia, Pierre Demarque

Abstract: 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 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. △ Less

Submitted 27 October, 2008; originally announced October 2008.

Comments: 4 pages, 3 figures, in proceedings of GONG2008/SOHO XXI: Solar-stellar dynamos as revealed by helio- and asteroseismology, August 11-15, 2008 - Boulder, Colorado

Two-Dimensional Stellar Evolution Code Including Arbitrary Magnetic Fields. II. Precision Improvement and Inclusion of Turbulence and Rotation

Authors: Linghuai Li, Sabatino Sofia, Paolo Ventura, Valentina Penza, Shaolan Bi, Sarbani Basu, Pierre Demarque

Abstract: 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 th… ▽ More 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. △ Less

Submitted 27 October, 2008; originally announced October 2008.

Comments: 56 pages, 8 figures, submitted to ApJS

The old open cluster NGC 2112: updated estimates of fundamental parameters based on a membership analysis

Authors: Giovanni Carraro, Sandro Villanova, Pierre Demarque, Christian Moni Bidin, Ginny McSwain

Abstract: 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 diff… ▽ More 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. △ Less

Submitted 21 February, 2008; originally announced February 2008.

Comments: 12 pages, 10 eps figure, in press in MNRAS

A cgi synthetic CMD calculator for the YY Isochrones

Authors: P. Demarque, S. N. Virani, E. J. Murphy, J. -H. Woo, Y. -C. Kim, S. K. Yi

Abstract: 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… ▽ More 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). △ Less

Submitted 2 January, 2008; originally announced January 2008.

Comments: 10 pages, 2 figures; cgi CMD calculator available at http://www.astro.yale.edu/demarque/yyiso.html Electronic preprint only

YREC: The Yale Rotating Stellar Evolution Code

Authors: P. Demarque, D. B. Guenther, L. H. Li, A. Mazumdar, C. W. Straka

Abstract: 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… ▽ More 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. △ Less

Submitted 22 October, 2007; originally announced October 2007.

Comments: 10 pages, 2 figures, ApSS accepted

Journal ref: Astrophys.SpaceSci.316:31-41,2008

Short-Period Variables in the Local Group Dwarf Galaxies Tucana and LGS3

Authors: E. J. Bernard, M. Monelli, C. Gallart, A. Aparicio, G. Bertelli, S. Cassisi, A. A. Cole, P. Demarque, A. E. Dolphin, I. Drozdovsky, H. C. Ferguson, S. Hidalgo, M. Mateo, L. Mayer, J. Navarro, F. Pont, E. D. Skillman, P. B. Stetson, E. Tolstoy

Abstract: 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… ▽ More 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. △ Less

Submitted 25 January, 2007; originally announced January 2007.

Comments: 4 pages, 3 figures. To appear in the proceedings for the VII Scientific Meeting of the Spanish Astronomical Society (SEA) held in Barcelona, September 12-15, 2006

Asteroseismic Diagnostics of Stellar Convective Cores

Authors: Anwesh Mazumdar, Sarbani Basu, Braxton L. Collier, Pierre Demarque

Abstract: 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 de… ▽ More 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. △ Less

Submitted 15 August, 2006; originally announced August 2006.

Comments: Accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.372:949-958,2006

On Dissipation inside Turbulent Convection Zones from 3D Simulations of Solar Convection

Authors: Kaloyan Penev, Dimitar Sasselov, Frank Robinson, Pierre Demarque

Abstract: 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… ▽ More 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. △ Less

Submitted 3 July, 2006; originally announced July 2006.

Comments: 13 pages, 4 figures

Journal ref: Astrophys.J.655:1166-1171,2007

The hot stars in orbit around the M31 central supermassive black hole: are they young or old?

Authors: Pierre Demarque, Shanil Virani

Abstract: 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 collision… ▽ More 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. △ Less

Submitted 5 October, 2006; v1 submitted 13 March, 2006; originally announced March 2006.

Comments: 8 pages, 2 figures, accepted by A&A for publication, in press

NGC 6791: an exotic open cluster or the nucleus of a tidally disrupted galaxy?

Authors: Giovanni Carraro, Sandro Villanova, Pierre Demarque, M. Virginia McSwain, Giampaolo Piotto, Luigi R. Bedin

Abstract: 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 abun… ▽ More 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. △ Less

Submitted 30 December, 2005; originally announced December 2005.

Comments: 29 pages, 10 eps figure, accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.643:1151-1159,2006

2-D Stellar Evolution Code Including Arbitrary Magnetic Fields. I. Mathematical Techniques and Test Cases

Authors: L. H. Li, P. Ventura, S. Basu, S. Sofia, P. Demarque

Abstract: 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 t… ▽ More 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. △ Less

Submitted 27 January, 2006; v1 submitted 8 November, 2005; originally announced November 2005.

Comments: 88 pages, 8 figures, accepted for publication in ApJS

The Non-Homologous Nature of Solar Diameter Variations

Authors: Sabatino Sofia, Sarbani Basu, Pierre Demarque, Linghuai Li, Gerard Thuillier

Abstract: 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 support… ▽ More 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. △ Less

Submitted 27 September, 2005; originally announced September 2005.

Comments: 11 pages, 3 figures, aastex style, accepted for publication by ApJL

Journal ref: Astrophys.J.623:L147-L150,2005

Space and Ground Based Pulsation Data of Eta Bootis Explained with Stellar Models Including Turbulence

Authors: Christian W. Straka, Pierre Demarque, D. B. Guenther, Linghuai Li, Frank J. Robinson

Abstract: 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 ou… ▽ More 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. △ Less

Submitted 14 September, 2005; originally announced September 2005.

Comments: 13 pages, 7 figures, ApJ in press

Journal ref: Astrophys.J.636:1078-1086,2006

Simulating the outer layers of Procyon A: a comparison with the Sun

Authors: F. J. Robinson, P. Demarque, D. B. Guenther, Y. -C. Kim, K. L. Chan

Abstract: 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… ▽ More 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. △ Less

Submitted 30 June, 2005; originally announced July 2005.

Comments: 18 pages 10 figures

Journal ref: Mon.Not.Roy.Astron.Soc.362:1031-1037,2005

Core Overshoot: An Improved Treatment and Constraints from Seismic Data

Authors: Christian W. Straka, Pierre Demarque, D. B. Guenther

Abstract: 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 o… ▽ More 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. △ Less

Submitted 1 May, 2005; originally announced May 2005.

Comments: 45 pages, 26 figures, accepted in ApJ

Journal ref: Astrophys.J. 629 (2005) 1075-1090

Super Helium-Rich Population and the Origin of Extreme Horizontal-Branch Stars in Globular Clusters

Authors: Young-Wook Lee, Seok-Joo Joo, Sang-Il Han, Chul Chung, Chang H. Ree, Young-Jong Sohn, Yong-Cheol Kim, Suk-Jin Yoon, Sukyoung K. Yi, Pierre Demarque

Abstract: 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 larg… ▽ More 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. △ Less

Submitted 23 January, 2005; originally announced January 2005.

Comments: 4 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters

Journal ref: Astrophys.J. 621 (2005) L57-L60

Y2 Isochrones with an Improved Core Overshoot

Authors: Pierre Demarque, Jong-Hak Woo, Yong-Cheol Kim, Sukyoung Yi

Abstract: 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 c… ▽ More 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. △ Less

Submitted 1 September, 2004; originally announced September 2004.

Comments: Accepted for publication in ApJS (December issue)

Journal ref: Astrophys.Space Sci. 291 (2004) 261-262

Asteroseismic determination of helium abundance in stellar envelopes

Authors: Sarbani Basu, Anwesh Mazumdar, H. M Antia, Pierre Demarque

Abstract: 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 env… ▽ More 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 Γ_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 σ_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}. △ Less

Submitted 15 February, 2004; originally announced February 2004.

Comments: To appear in MNRAS. 11 pages, 11 figures

Journal ref: Mon.Not.Roy.Astron.Soc. 350 (2004) 277

A Preliminary Seismic Analysis of 51 Peg: Large and Small Spacings from Standard Models

Authors: Eric J. Murphy, Pierre Demarque, D. B. Guenther

Abstract: 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… ▽ More 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. △ Less

Submitted 19 December, 2003; originally announced December 2003.

Comments: 22 pages, 5 figures, 3 tables. Accepted for publicaton by ApJ

Journal ref: Astrophys.J. 605 (2004) 472-479

3D convection simulations of the outer layers of the Sun using realistic physics

Authors: F. J. Robinson, P. Demarque, L. H. Li, S. Sofia, Y. -C. Kim, K. L. Chan, D. B. Guenther

Abstract: 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 o… ▽ More 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. △ Less

Submitted 12 December, 2002; originally announced December 2002.

Comments: 16 pages, 22 figures, mnras style file (mn.sty), accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 340 (2003) 923

The Y^2 Stellar Evolutionary Tracks

Authors: Sukyoung K. Yi, Yong -Cheol Kim, Pierre Demarque

Abstract: 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 $α$-enha… ▽ More 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 $α$-enhancement). △ Less

Submitted 9 October, 2002; originally announced October 2002.

Comments: 7 pages, TeX, 1 eps figure. ApJS, 2003, vol.144 (Feb), in press

Testing intermediate-age stellar evolution models with VLT photometry of LMC clusters. I. The data

Authors: C. Gallart, M. Zoccali, G. Bertelli, C. Chiosi, P. Demarque, L. Girardi, E. Nasi, J. -H. Woo, S. Yi

Abstract: 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 t… ▽ More 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. △ Less

Submitted 8 August, 2002; originally announced August 2002.

Comments: 27 pages with 10 figures. Accepted by the Astronomical Journal

The YY Isochrones for alpha-element Enhanced Mixtures

Authors: Yong -Cheol Kim, Pierre Demarque, Sukyoung K. Yi, David R. Alexander

Abstract: 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… ▽ More 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 △ Less

Submitted 7 August, 2002; originally announced August 2002.

Comments: LaTeX, 11 Figures, 3 Tables. The Astrophysical Journal Suppl, 2002, accepted

Journal ref: Astrophys.J.Suppl.143:499-512,2002

Testing intermediate-age stellar evolution models with VLT photometry of LMC clusters. II. Analysis with the Yale models

Authors: J. -H. Woo, C. Gallart, P. Demarque, S. Yi, M. Zoccali

Abstract: 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 reprod… ▽ More 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. △ Less

Submitted 6 August, 2002; originally announced August 2002.

Comments: 19 pages, 13 figures, accepted for the publication in AJ

Seismic Test of Solar Models, Solar Neutrinos and Implications for Metal-Rich Accretion

Authors: R. A. Winnick, Pierre Demarque, Sarbani Basu, D. B. Guenther

Abstract: 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 helioseismol… ▽ More 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). △ Less

Submitted 29 May, 2002; v1 submitted 5 November, 2001; originally announced November 2001.

Comments: Revision: 23 pages, 8 figures, 3 tables, to be published in ApJ, Sept 2002. Used newer neutrino cross sections, modified text accordingly, other minor revisions as well

Journal ref: Astrophys.J. 576 (2002) 1075-1084

Inclusion of turbulence in solar modeling

Authors: L. H. Li, F. J. Robinson, P. Demarque, S. Sofia, D. B. Guenther

Abstract: 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 tu… ▽ More 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. △ Less

Submitted 7 November, 2001; v1 submitted 5 September, 2001; originally announced September 2001.

Comments: 10 pages, 12 figures

Journal ref: Vol 567, no. 1, p. 1, ApJ March 1, 2002

Asteroseismology of delta Scuti stars - a parameter study and application to seismology of FG Virginis

Authors: M. R. Templeton, S. Basu, P. Demarque

Abstract: We assess the potential of asteroseismology for determining the fundamental properties of individual $δ$ 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, an… ▽ More We assess the potential of asteroseismology for determining the fundamental properties of individual $δ$ 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. △ Less

Submitted 28 August, 2001; originally announced August 2001.

Comments: To appear in Astrophysical Journal