Source author record

Pierre Demarque

Pierre Demarque appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2016arXiv

Entropy in adiabatic regions of convection simulations

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.

preprint2015arXiv

Hubble Space Telescope Astrometry of the Procyon System

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.

preprint2015arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

Variation of Stellar Envelope Convection and Overshoot with Metallicity

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.

preprint2012arXiv

Comparing the Effect of Radiative Transfer Schemes on Convection Simulations

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.

preprint2012arXiv

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

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.

preprint2011arXiv

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

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.

preprint2009arXiv

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

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.

preprint2009arXiv

Dissipation Efficiency in Turbulent Convective Zones in Low Mass Stars

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.

preprint2008arXiv

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

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.

preprint2006arXiv

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

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.

preprint2005arXiv

The Non-Homologous Nature of Solar Diameter Variations

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.

preprint2002arXiv

The YY Isochrones for alpha-element Enhanced Mixtures

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

preprint1997arXiv

On the Origin of the UV Upturn in Elliptical Galaxies. I. Sensitivity of UV Population Synthesis to Various Input Parameters

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

preprint1997arXiv

On the Origin of the UV Upturn in Elliptical Galaxies. II. Test of the HB Hypothesis

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

preprint1997arXiv

The Relative Ages of Galactic Globular Clusters

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